El análisis de la interacción entre la radiación electromagnética y un sistema cuántico compuesto por un gas de átomos hidrogenoides se realiza inicialmente determinando el hamiltoniano de interacción desde una perspectiva cuántica entre un sistema de espines y bosones. Lo novedoso de la determinación realizada aquí es que consideramos la polarización circular para los bosones. El hamiltoniano clásico de la interacción espín-bosón se determina como el valor esperado del hamiltoniano cuántico, donde los vectores de estados están definidos como los estados coherentes del sistema de espínes.
El análisis de la interacción entre la radiación electromagnética y un sistema cuántico compuesto por un gas de átomos hidrogenoides se realiza inicialmente determinando el hamiltoniano de interacción desde una perspectiva cuántica entre un sistema de espines y bosones. Lo novedoso de la determinación realizada aquí es que consideramos la polarización circular para los bosones. El hamiltoniano clásico de la interacción espín-bosón se determina como el valor esperado del hamiltoniano cuántico, donde los vectores de estados están definidos como los estados coherentes del sistema de espínes.
La realización experimental del condensado de Bose-Einstein ha resaltado la importancia de controlar y manipular, con facilidad, sistemas cuánticos cada vez más complejos. En este contexto, el estudio de átomos ultra-fríos cargados en redes ópticas ha demostrado cumplir con estas exigencias, proporcionando escenarios ideales para investigar sistemas cuánticos muy complejos, que solo pueden ser analizados a través de simuladores cuánticos. En este trabajo, se implementa numéricamente el teorema de Floquet usando el método del espacio de Hilbert extendido. El objetivo es estudiar la dinámica de uno y dos bosones de espín 1 en una red óptica unidimensional, en el régimen de superfluido, utilizando como base el hamiltoniano de Bose-Hubbard de espín 1. Además, se aborda el estudio de dicho sistema bajo la influencia de campos de Zeeman lineales o cuadráticos externos, conducidos periódicamente de manera tanto lineal como cuadrática. En primer lugar, se presenta el modelo a estudiar. Comenzamos con el modelo de Bose-Hubbard sin espín para, posteriormente, incorporar el grado de libertad de espín 1 junto con sus interacciones e incluir campos de Zeeman, tanto lineales como cuadráticos, acoplados a un driving periódico, lo que sitúa a nuestro sistema en un marco dinámico. Posteriormente, se introduce la representación bosónica de SU(3) para espín 1 como grupo de simetría, junto con las bases, estados y observables relevantes necesarios para caracterizar completamente nuestro sistema. Para analizar eficazmente este sistema, se emplea una técnica numérica basada en la teoría de Floquet. Se describe el teorema de Floquet, sus implicaciones y sus analogías con el teorema de Bloch. Luego, se presenta el método numérico principal utilizado en nuestro análisis: el método del espacio de Hilbert extendido, que aprovecha la naturaleza periódica de nuestro sistema para transferir el problema al espacio armónico de Fourier, donde la evolución del estado se puede determinar fácilmente. Este método se valida resolviendo analíticamente, de manera exacta, nuestro sistema de una partícula en dos sitios y se compara con nuestra implementación numérica, demostrando una coincidencia sin errores discernibles. Entre los resultados obtenidos, se muestra que, independientemente de la presencia de interacciones cuadrupolares o incluso de interacciones del tipo densidad entre partículas, cualquier estado de espín 1 muestra inherentemente tanto orden cuadrupolar como espinorial. De hecho, se muestra que, cuando el estado inicial es un estado local, sin importar el campo aplicado, la dinámica del orden cuadrupolar, junto con la espinorial, puede condensarse en la desigualdad 10/9 ≤ ⟨Ŝ⟩²t + ⟨Q⟩²t ≤ 4/3, donde el límite superior e inferior resultan ser los valores propios de los operadores cuadrático y cúbico de Casimir del grupo SU(3), respectivamente. Adicionalmente, se muestra que la presencia de campos magnéticos genera nuevas fases y fenómenos en los sistemas, tales como la ruptura de la simetría de invarianza temporal en el caso de una sola partícula bajo un campo lineal, lo que da lugar a fenómenos como la precesión del espín o la resonancia del spin-mixing en el caso de dos partículas bajo un campo cuadrático, donde la intensidad del campo y la interacción entran en fase, mostrando oscilaciones máximas y rápidas que no decaen en el tiempo. Finalmente, se concluye que el método numérico del espacio de Hilbert extendido, basado en el teorema de Floquet, es una herramienta robusta y de alto nivel que proporciona una forma directa y sencilla de estudiar sistemas cuánticos de manera dinámica, incluso en presencia de interacciones complejas y espacios de Hilbert grandes. Además, demuestra una versatilidad notable para analizar los sistemas desde distintos regímenes impuestos por la periodicidad del driving acoplado, tales como el régimen de altas frecuencias, bajo la aproximación de onda rotante, el régimen fuerte, entre otros.
En esta tesis doctoral estudiamos el fenomeno de la polarizacion para campos cuanticos. En el caso de campos clasicos el estado de polarizacion es elegantemente visualizado por medio de la esfera de Poincare, y la polarizacion se cuantifica por medio de los parametros de Stokes. Estos son cantidades medibles que permiten la ordenacion de los estados de la luz segun su grado de polarizacion. Hasta cierto punto, este formalismo es extensible al caso cuantico, donde se definen los operadores de Stokes, cuyos valores esperados coinciden con los parametros de Stokes. Por analogia con el caso clasico, en el cual el grado de polarizacion se define como la longitud del vector de Stokes, podriamos definir en optica cuantica el grado de polarizacion de esta misma forma sustituyendo los parametros de Stokes por el valor promedio de los operadores de Stokes. Sin embargo, esta definicion tiene importantes inconvenientes. Se trata de una cantidad clasica que no distingue entre campos cuanticos con significativas diferencias en cuanto a su polarizacion. La apreciacion de los sutiles efectos que surgen en el mundo cuantico requiere del escrutinio de las fluctuaciones de polarizacion de mayor orden (mas alla de las correlaciones de segundo orden en los campos que son las unicas tenidas en cuenta en el tratamiento estandar de la polarizacion en optica clasica). La forma mas inmediata de conseguir esto es expandiendo el sector de polarizacion del estado (la parte del operador densidad que es accesible mediante medidas de polarizacion) en multipolos. Una vez hecho el desarrollo en multipolos, vemos como estos estan relacionados con las cantidades medibles en el laboratorio. La distribucion de probabilidad para estos multipolos proporciona asi una informacion completa sobre las propiedades de polarizacion de cualquier estado y, en terminos de estos, proponemos una medida adecuada para la evaluacion cuantitativa de las fluctuaciones de la polarizacion. En particular, resulta conveniente la introduccion de lo que denominamos distribucion acumulativa, que contiene la informacion de polarizacion hasta el orden M. Es remarcable que, dentro de un subespacio con un numero definido de fotones, los estados coherentes de SU(2) son los estados que alcanzan el valor maximo de la distribucion acumulativa para cualquier orden M. El caso contrario en el que esta distribucion se anula hasta dicho orden es el de los estados despolarizados a M orden. Existe un amplio consenso en considerar los estados coherentes como los estados con propiedades mas clasicas. Es irresistible preguntarse por los estados que minimizan esta distribucion acumulativa, ya que estos estados son lo opuesto a los estados coherentes de SU(2) y, por tanto, son los estados mas no-clasicos permitidos por la teoria cuantica. En esta tesis calculamos esos estados, a los que denominamos Reyes Cuanticos, por medio de metodos computacionales que hacen uso de las bases de Grobner. Ademas, empleamos la representacion de Majorana para obtener una valiosa imagen geometrica de estos estados. Mientras la representacion de Majorana para un estado coherente de SU(2) corresponde a un unico punto (degenerado) en la esfera unidad, los Reyes Cuanticos estan representados por constelaciones de Majorana con los puntos situados de una forma muy simetrica alrededor de la esfera. Estudiamos las relaciones entre estas constelaciones de Majorana y los t-disenos esfericos. La evidencia numerica nos hace conjeturar que el maximo orden de despolarizacion posible para un estado de N fotones coincide con el grado t de un t-diseno esferico optimo con N puntos. Finalmente, estudiamos las posibles aplicaciones de estos estados extremos con especial atencion a su utilidad en metrologia cuantica.
En este trabajo, desarrollado en el marco teórico detrás de los Matrix Product States (MPS) y Matrix Product Operators (MPO) con un enfoque en su implementación numérica en un lenguaje de programación de bajo nivel como lo es C++, se obtiene como resultado una librería que integra ambos métodos para computar eficientemente observables locales, correladores y operadores no triviales que involucran múltiples sitios, tales como los Hamiltonianos que incluyen términos de hopping o de interacción. La implementación es general y puede ser empleada para una variedad de sistemas. En este trabajo se hace el planteamiento teórico de un sistema conformado por un arreglo unidimensional de bosones de espín-1 en el régimen superfluido, descrito por el Hamiltoniano de Bose-Hubbard en presencia de un acoplamiento espín-órbita efectivo (SOC) modelado por un campo magnético externo. Expresar el sistema en la base del operador Ŝᵧ para alinear el sistema con la orientación del campo externo introduce una nueva base de estados helicoidales, haciendo emerger propiedades físicas novedosas cuyo estudio es de interés en el campo de los átomos ultrafríos. Adicionalmente, se implementan las técnicas de MPS y MPO en el estudio del modelo de Heisenberg cuadrático y de un sistema de espín-1/2 sujeto a interacción espín-órbita como prueba preliminar para la futura implementación en el sistema de espín-1 en el régimen superfluido. El análisis de estados base de múltiples sistemas dentro de distintos regímenes de parámetros provee información física relevante a través de correlaciones, operadores locales y energías de estado base, permitiendo indagar en los detalles característicos de cada una de las fases y modelos estudiados. Este trabajo demuestra que la inclusión del formalismo del MPO a la técnica convencional del MPS permite representar de manera eficiente operadores no triviales que exhiben diferentes rangos de complejidad, todo esto mientras se reduce el coste computacional de las simulaciones. Este formalismo abre direcciones prometedoras para futuros desarrollos, incluyendo algoritmos paralelizados de time evolution block decimation (TEBD) y protocolos de contracción optimizados para sistemas de mayor dimensión.
Los corrales cuánticos presentan propiedades interesantes debido a la interesante \ncombinación de confinamiento y focalización algunas de las cuales fueron observadas experimentalmente. Estudiamos el comportamiento estático y dinámico de estos sistemas. \nEn primer lugar analizamos el caso del corral cuántico no interactuante (sin impurezas), y con parámetros realistas obtuvimos los autoestados la elipse con el modelo de pared dura. Posteriormente estudiamos el sistema en presencia de dos impurezas de Kondo de espín S = 1/2 localizadas en los focos de la elipse. Éstas interactúan con los electrones itinerantes en la elipse vía un término de superintercambio J. Diagonalizamos este sistema numéricamente y estudiamos propiedades tales como correlación de espín y correlaciones dinámicas. Encontramos que, para valores de J chicos comparado con el ancho de banda, \nlos espines están en un estado singlete para número par de electrones en la elipse o triplete para número impar (régimen RKKY). En ese límite es posible describir el comportamiento a bajas energías con un hamiltoniano efectivo entre las impurezas. Para valores grandes de J más estados electrónicos están involucrados y los espines se descorrelacionan formando un estado de Kondo local con los electrones itinerantes. Estudiamos algunas magnitudes de la teoría de la información cuántica que es una herramienta alternativa para el análisis de sistemas cuánticos. Calculamos la entropía de von Neumann para diversas particiones del \nsistema y caracterizamos el cambio de régimen de RKKY al de Kondo ante el incremento de J, de acuerdo a los resultados obtenidos con el cálculo de funciones de correlaciones estáticas y dinámicas. \nAdicionalmente consideramos fluctuaciones de carga en los estados de las impurezas y analizamos propiedades estáticas y dinámicas en función de los parámetros relevantes del sistema. El modelo estudiado consiste de dos impurezas de Anderson localizadas en los focos e hibridizadas con magnitud V a los electrones en la elipse. Con el estudio de la correlación de espín entre las impurezas, y entre una impureza y los electrones en un mismo foco, caracterizamos los regímenes RKKY y Kondo. Esto está de acuerdo con los resultados obtenidos con la densidad de estados electrónica calculada usando un modelo más realista que incluye la hibridización con estados del continuo. \nEncontramos una relación entre el entrelazamiento entre las impurezas de espín y \nobservables físicos del sistema. Con ésto obtuvimos el diagrama de fases de entrelazamiento entre las impurezas para llenado par e impar, que presenta regiones en las que están entrelazadas y otras de estados separables.Mostramos resultados numéricos que sustentan las relaciones encontradas. \nComo aplicación adicional de los conceptos de información cuántica al estudio de las propiedades físicas de sistemas correlacionados, analizamos una cadena de espín S = 1 con condiciones de contorno abiertas y con interacción a primeros vecinos bilineal y bicuadrática (β). En particular estudiamos los modelos de Heisenberg (β = 0) y AKLT (β = 1/3). Realizamos diversos estudios en sistemas bipartitos puros y mixtos. Para caracterizar el entrelazamiento en un sistema bipartito puro usamos la entropía de von Neumann, en cambio para un sistema bipartito mixto, usamos la negatividad. Estos estudios nos permitieron estudiar las excitaciones de espín 1/2 localizadas en los bordes de la cadena [87, 88]. \nUn resultado interesante aparece cuando analizamos el entrelazamiento entre dos \npartes del sistema; que es máximo entre los extremos de la cadena y disminuye hacia el centro, reflejando el carácter fraccionario y la localización de las excitaciones de espín del modelo, que han sido observadas experimentalmente. \nEn conclusión, usando tanto técnicas tradicionales (correlaciones estáticas y dinámicas) como novedosas (información cuántica) para el estudio de sistemas interactuantes, investigamos el comportamiento de algunos sistemas nanoscópicos paradigmáticos. Para ésto utilizamos técnicas tanto analíticas como numéricas de actualidad. Esperamos que los resultados obtenidos en esta tesis sean útiles para una mejor comprensión de estos sistemas, así como para posibles aplicaciones nanotecnológicas. \n
El tema central de la tesis doctoral es el análisis de los sistemas cuánticos abiertos. Estos sistemas se caracterizan por estar sometidos a la interacción con el entorno, lo que provoca que su evolución deje de ser unitaria. Es por tanto necesario considerar modelos más allá de la ecuación de Schrödinger. Los sistemas cuánticos abiertos aparecen en numerosos campos, como la Física del Estado Sólido y la Dinámica Molecular. Por este motivo, un análisis detallado de sus propiedades y su dinámica es un tema digno de estudio con un gran abanico de aplicaciones. El enfoque elegido en esta tesis es el desarrollo de un formalismo geométrico que describa de forma adecuada las características de los sistemas cuánticos abiertos. La geometría diferencial ha demostrado ser una herramienta muy útil en el análisis de sistemas físicos. Desde mediados del siglo XX, se ha desarrollado con gran éxito una descripción geométrica de la Mecánica Clásica, principalmente en torno a las mecánicas lagrangiana y hamiltoniana. Por este motivo, resulta natural describir también los sistemas cuánticos en términos geométricos. Las ventajas de un formalismo geométrico resultan claras. Cuando tanto los sistemas clásicos como los cuánticos se describen en los mismos términos, es sencillo describir situaciones en las que existan interacciones clásico-cuánticas. Éste es el caso, por ejemplo, de muchos modelos de Dinámica Molecular, en los que los núcleos y los electrones son considerados respectivamente como partículas clásicas y cuánticas. Por otra parte, un formalismo geométrico de la Mecánica Cuántica posibilita un mejor entendimiento de las diferencias intrínsecas entre las teorías clásicas y las cuánticas. Dada su relevancia a lo largo de la tesis, el Capítulo 1 está enfocado al resumen y el análisis de la descripción geométrica de la imagen de Schrödinger de la Mecánica Cuántica. En formulación usual (algebraica), esta imagen se basa en la representación de los estados de los sistemas cuánticos mediante vectores en un espacio de Hilbert complejo. La transición a una formulación basada en geometría diferencial es inmediata para sistemas finito-dimensionales, dado que los espacios lineales de dimensión finita son casos triviales de variedades diferenciables. Las estructuras adicionales, en concreto el producto hermítico propio de los espacios de Hilbert y los escalares complejos, se describen mediante campos tensoriales en dichas variedades diferenciables, formando lo que se conoce como una estructura Kähler. Todos los ingredientes necesarios para el análisis de sistemas cuánticos pueden describirse en estas variedades de Kähler. Los observables se representan mediante funciones diferenciables, mientras que la dinámica se describe mediante curvas integrales de campos vectoriales hamiltonianos respecto a la forma simpléctica de la estructura Kähler. Esta caracterización puede llevarse a cabo en el espacio de Hilbert asociado a cualquier sistema cuántico finito-dimensional. Además, es posible analizar las propiedades geométricas del espacio proyectivo de Hilbert, el cual constituye el conjunto de estados puros del sistema. Su estructura Kähler puede ser deducida mediante un proceso de reducción de la estructura previamente obtenida en el espacio de Hilbert. Como aspecto novedoso, la tesis presenta en detalle este proceso de reducción, dotándolo de una descripción matemática adecuada. Todas las características de la imagen de Schrödinger pueden describirse de forma geométrica, lo que permite utilizar nuevas herramientas en el análisis de los sistemas cuánticos. Éste es precisamente el tema principal del Capítulo 2, en el cual se utilizan los sistemas de Lie-Kähler para resolver la ecuación de Schrödinger dependiente del tiempo. En geometría diferencial, un sistema de Lie es un sistema no homogéneo de ecuaciones diferenciales que admite una regla de superposición. En general, la obtención de esta regla de superposición es una ardua tarea, la cual puede aligerarse en presencia de estructuras adicionales que sean preservadas por la acción del sistema de Lie. De esta forma, según la estructura preservada, es posible hablar de sistemas de Lie-Hamilton, Lie-Dirac, etc. En el caso de la Mecánica Cuántica, una ecuación de Schrödinger dependiente del tiempo es un sistema de Lie que preserva la estructura de Kähler previamente descrita. Es por tanto un nuevo tipo de sistema de Lie, al que resulta natural denominar de Lie-Kähler. El Capítulo 2 presenta las propiedades de estos nuevos sistemas y describe un método riguroso para la obtención de sus reglas de superposición. El formalismo geométrico puede extenderse más allá de la imagen de Schrödinger, lo que resulta necesario en el contexto de los sistemas cuánticos abiertos, dado que tanto estados puros como estados mezcla son necesarios para su descripción. Por este motivo, el Capítulo 3 resume la imagen de Heisenberg de la Mecánica Cuántica y su representación de los estados puros y mezcla como funcionales lineales en el álgebra de Lie-Jordan de observables. Se explica también como las estructuras algebraicas de los observables pueden representarse geométricamente en el espacio dual de funcionales lineales en el álgebra. Este es el punto de partida de una de las principales contribuciones de la tesis. Un proceso de reducción, similar al realizado en el análisis de la imagen de Schrödinger, permite describir las propiedades geométricas de la variedad de estados puros y mezcla del sistema. De esta forma, se obtienen dos campos tensoriales en la variedad de estados, los cuales representan correctamente las estructuras algebraicas de los observables. Puede por tanto concluirse que el formalismo geométrico presentado en la tesis es completamente equivalente a la descripción algebraica tradicional, ya que se logra describir adecuadamente todas las propiedades de los sistemas cuánticos. Además, el formalismo geométrico ofrece una visión más clara de las propiedades intrínsecas de la Mecánica Cuántica, lo que facilita una mejor comprensión de la teoría. Por otra parte, un análisis geométrico de la variedad de estados permite estudiar su estratificación y sus propiedades. La tesis demuestra que se trata de una variedad con borde, cuyos puntos extremales son precisamente los estados puros del sistema. La estratificación de esta variedad resulta importante a la hora de considerar la dinámica inducida por campos gradiente y hamiltonianos. Con el objetivo de ilustrar todas estas propiedades, se analizan unos casos sencillos pero con relevancia física. A lo largo de la tesis, se muestran diversas aplicaciones del formalismo geométrico al análisis de sistemas cuánticos abiertos. El Capítulo 4 presenta la descripción de la evolución markoviana de sistemas cuánticos abiertos. Se dice que una evolución es markoviana si depende únicamente en el estado actual del sistema y no de los estados en instantes anteriores, es decir, si el sistema "no tiene memoria". En Mecánica Cuántica, la evolución markoviana se obtiene a partir de la ecuación de Kossakowski-Lindblad, una ecuación diferencial de primer orden en la variedad de estados puros y mezcla de un sistema cuántico abierto. El formalismo geométrico describe esta ecuación como un campo tensorial en esta variedad, lo que permite analizar las propiedades de sus curvas integrales. De esta manera, es posible considerar diversos aspectos de la evolución markoviana desde un punto de vista geométrico. Cualquier evolución no-unitaria determina un cambio en las propiedades algebraicas de los observables cuánticos, lo que puede resultar en una contracción del álgebra. En términos geométricos, esta contracción puede entenderse mediante el límite de una familia de campos tensoriales definida por el flujo del campo tensorial de Kossakowski-Lindblad. Otra característica importante de esta evolución es la existencia de variedades límite. Sus propiedades pueden determinarse gracias a la estructura afín existente, lo que a su vez permite investigar su relación con las contracciones de álgebras de observables. Por último, se ofrece una descripción geométrica de los problemas de control de sistemas cuánticos abiertos. Un análisis geométrico de la Mecánica Cuántica permite aplicar a estos problemas los resultados de la teoría de control de grupos de Lie. Como consecuencia, es posible realizar una clasificación de los sistemas cuánticos abiertos según sus propiedades de controlabilidad. Otro ejemplo de sistemas cuánticos abiertos aparece en el contexto de la Dinámica Molecular. En el estudio de sistemas moleculares, debido al gran número de partículas presentes, la ecuación de Schrödinger no puede ser resuelta ni siquiera por métodos numéricos. Por tanto, resulta útil considerar aproximaciones a la ecuación de Schrödinger. En particular, existen muchos modelos que consideran un comportamiento clásico de algunas de las partículas, normalmente los núcleos. El Capítulo 5 resume las propiedades estos modelos moleculares, y en particular del conocido como modelo de Ehrenfest. Es posible llevar a cabo una descripción geométrica de este modelo, basándose en las descripciones de los subsistemas clásico y cuántico. Como resultado, las ecuaciones del modelo pueden escribirse como ecuaciones hamiltonianas en una variedad de Poisson. A partir de estas propiedades, la tesis presenta una generalización del modelo de Ehrenfest a distribuciones estadísticas. Este es un paso importante, ya que se demuestra que esta descripción estadística predice la aparición de efectos relacionados con el fenómeno de decoherencia, algo que no ocurre en el modelo de Ehrenfest estándar. Se han realizado simulaciones numéricas, cuyos resultados respaldan la descripción de sistemas moleculares mediante el modelo estadístico de Ehrenfest. Por último, en este contexto resulta posible considerar distribuciones estadísticas con temperatura. La tesis presenta estas dist
El objetivo de este trabajo es realizar un estudio sobre las transiciones de fase cuánticas que se presentan en el modelo de Rabi, un modelo de interacción radiación-materia entre un átomo de dos niveles y un campo bosónico. Presentaremos y analizaremos su hamiltoniano para, a partir del mismo, caracterizar las transiciones de fase cuánticas que afectan tanto al estado fundamental como a los estados excitados del sistema. Primero analizaremos la transición en el estado fundamental, realizando el estudio cuántico mediante cálculos numéricos y luego haciendo un estudio analítico dentro de una aproximación semiclásica. Ello nos permitirá comparar los resultados y comprobar la validez de esta aproximación. Después, analizaremos la transición en los estados excitados, comparando de nuevo el estudio analítico con los cálculos numéricos, exactos. Veremos que la transición de fase en el estado fundamental se debe al cambio del potencial que observa el sistema, causando la rotura de la simetría de dicho estado. Por otro lado, comprobaremos que la transición en los estados excitados solo puede ocurrir cuando ya se ha dado la transición en el estado fundamental, y veremos que se caracteriza por la degeneración en niveles de energía inferiores a cierto valor crítico. En todos los casos, comprobaremos que los cálculos semiclásicos y los exactos coinciden cuando el cociente entre las energías que caracterizan al átomo y al campo bosónico tiene un valor grande comparado con 1. Para conseguirlo, hemos desarrollado desde cero diversos códigos MATLAB con el objetivo de poder realizar los cálculos necesarios para caracterizar las transiciones. En concreto, para estudiar la transición de fase en el estado fundamental, hemos desarrollado un código que permite obtener la energía del estado fundamental del sistema y otro que permite calcular, para dicho estado, el valor esperado del número de fotones y de σz, observables que se discutirán en el texto. También hemos creado un código que muestra la forma del potencial que afecta al sistema en su estado fundamental. Del mismo modo, al estudiar la transición en los estados excitados, hemos creado un código que obtiene todos los niveles de energía del sistema, con el objetivo de calcular numéricamente la densidad de estados. Otro código permite comparar este cálculo con la expresión analítica que obtendremos de la densidad de estados, y otro último calcula el valor de σz en todos los estados del sistema. Por último, dos códigos nos permitirán crear imágenes que sirvan para ilustrar lo que sucede en las transiciones; el primero de ellos representa la forma del potencial que afecta al sistema en diversas situaciones; el segundo estudia la variación de los diez primeros niveles energéticos del sistema. Todo esto se desarrolla en el trabajo. Los códigos pueden consultarse en el apéndice.