Modelo de bosones sin espín en interacción en una red
El modelo de Bose-Hubbard proporciona una descripción de la física de los bosones sin espín que interactúan en una red . Está estrechamente relacionado con el modelo de Hubbard que se originó en la física del estado sólido como una descripción aproximada de los sistemas superconductores y el movimiento de los electrones entre los átomos de un sólido cristalino . El modelo fue introducido por Gersch y Knollman [1] en 1963 en el contexto de los superconductores granulares. (El término " Bose " en su nombre se refiere al hecho de que las partículas en el sistema son bosónicas ). El modelo saltó a la fama en la década de 1980 después de que se descubrió que capturaba la esencia de la transición superfluido-aislante de una manera que era mucho más manejable matemáticamente que los modelos fermiónicos metal-aislante. [2] [3] [4]
El modelo de Bose-Hubbard se puede utilizar para describir sistemas físicos como átomos bosónicos en una red óptica , [5] así como ciertos aislantes magnéticos. [6] [7] Además, se puede generalizar y aplicar a mezclas de Bose-Fermi, en cuyo caso el hamiltoniano correspondiente se denomina hamiltoniano de Bose-Fermi-Hubbard.
Hamiltoniano
La física de este modelo viene dada por el hamiltoniano de Bose-Hubbard:
Aquí, denota la suma de todos los sitios reticulares vecinos y , mientras que y son operadores de creación y aniquilación bosónica tales que da el número de partículas en el sitio . El modelo está parametrizado por la amplitud de salto que describe la movilidad de los bosones en la red, la interacción en el sitio que puede ser atractiva ( ) o repulsiva ( ), y el potencial químico , que esencialmente establece el número de partículas. Si no se especifica, normalmente la frase "modelo de Bose-Hubbard" se refiere al caso en el que la interacción en el sitio es repulsiva.
Este hamiltoniano tiene simetría global , lo que significa que es invariante (sus propiedades físicas no cambian) ante la transformación . En una fase superfluida , esta simetría se rompe espontáneamente .
Espacio de Hilbert
La dimensión del espacio de Hilbert del modelo de Bose-Hubbard está dada por , donde es el número total de partículas, mientras que denota el número total de sitios reticulares. En fijos o , la dimensión del espacio de Hilbert crece polinomialmente, pero en una densidad fija de bosones por sitio, crece exponencialmente como . Se pueden formular hamiltonianos análogos para describir fermiones sin espín (el modelo de Fermi-Hubbard) o mezclas de diferentes especies de átomos (mezclas de Bose-Fermi, por ejemplo). En el caso de una mezcla, el espacio de Hilbert es simplemente el producto tensorial de los espacios de Hilbert de las especies individuales. Normalmente se incluyen términos adicionales para modelar la interacción entre especies.
Diagrama de fases
A temperatura cero, el modelo de Bose-Hubbard (en ausencia de desorden) está en un estado aislante de Mott en pequeño , o en un estado superfluido en grande . [8] Las fases aislantes de Mott se caracterizan por densidades de bosones enteros, por la existencia de una brecha de energía para excitaciones de partículas-huecos y por compresibilidad cero . El superfluido se caracteriza por una coherencia de fase de largo alcance, una ruptura espontánea de la simetría continua del hamiltoniano, una compresibilidad distinta de cero y una susceptibilidad superfluida. A temperatura distinta de cero, en ciertos regímenes de parámetros aparece una fase fluida regular que no rompe la simetría y no muestra coherencia de fase. Ambas fases se han observado experimentalmente en gases atómicos ultrafríos. [9]
En presencia de desorden, existe una tercera fase, la de "vidrio de Bose". [4] El vidrio de Bose es una fase de Griffiths y puede considerarse como un aislante de Mott que contiene raros "charcos" de superfluido. Estos depósitos de superfluido no están interconectados, por lo que el sistema sigue siendo aislante, pero su presencia cambia significativamente la termodinámica del modelo. La fase de vidrio de Bose se caracteriza por una compresibilidad finita, la ausencia de un espacio y una susceptibilidad infinita a los superfluidos . [4] Es aislante a pesar de la ausencia de un espacio, ya que la tunelización baja impide la generación de excitaciones que, aunque cercanas en energía, están separadas espacialmente. El vidrio de Bose tiene un parámetro de orden Edwards-Anderson distinto de cero [10] [11] y se ha sugerido (pero no demostrado) que muestra una ruptura de simetría de réplica . [12]
Teoría del campo medio
Las fases del modelo limpio de Bose-Hubbard se pueden describir utilizando un hamiltoniano de campo medio : [13] donde es el número de coordinación reticular . Esto se puede obtener a partir del hamiltoniano de Bose-Hubbard completo estableciendo donde , ignorando los términos cuadráticos en (asumiblemente infinitesimales) y reetiquetando . Debido a que este desacoplamiento rompe la simetría del hamiltoniano inicial para todos los valores distintos de cero de , este parámetro actúa como un parámetro de orden superfluido . Para simplificar, se supone que este desacoplamiento es el mismo en cada sitio, lo que excluye fases exóticas como supersólidos u otras fases no homogéneas. (Son posibles otros desacoplamientos). El diagrama de fases se puede determinar calculando la energía de este hamiltoniano de campo medio utilizando la teoría de perturbación de segundo orden y encontrando la condición para la cual . Para hacer esto, el hamiltoniano se escribe como una pieza local del sitio más una perturbación: donde los términos bilineales y su conjugado se tratan como la perturbación. Se supone que el parámetro de orden es pequeño cerca de la transición de fase . El término local es diagonal en la base de Fock , lo que da la contribución de energía de orden cero: donde es un entero que etiqueta el llenado del estado de Fock. La pieza perturbativa se puede tratar con la teoría de perturbación de segundo orden, que conduce a: La energía se puede expresar como una expansión en serie en potencias pares del parámetro de orden (también conocido como el formalismo de Landau ): Después de hacerlo, la condición para la transición de fase de segundo orden de campo medio entre el aislante de Mott y la fase superfluida viene dada por: donde el entero describe el llenado del lóbulo aislante de Mott. Trazar la línea para diferentes valores enteros de genera el límite de los diferentes lóbulos de Mott, como se muestra en el diagrama de fases. [4]
Implementación en redes ópticas
Los átomos ultrafríos en redes ópticas se consideran una realización estándar del modelo de Bose-Hubbard. La capacidad de ajustar los parámetros del modelo mediante técnicas experimentales simples y la falta de dinámica de red presente en los sistemas electrónicos de estado sólido significan que los átomos ultrafríos ofrecen una realización limpia y controlable del modelo de Bose-Hubbard. [14] [5] La mayor desventaja de la tecnología de red óptica es la vida útil de la trampa, ya que los átomos suelen quedar atrapados solo durante unas pocas decenas de segundos.
Para entender por qué los átomos ultrafríos ofrecen una realización tan conveniente de la física de Bose-Hubbard, el hamiltoniano de Bose-Hubbard se puede derivar a partir del segundo hamiltoniano cuantizado que describe un gas de átomos ultrafríos en el potencial reticular óptico. Este hamiltoniano viene dado por:
- ,
donde es el potencial de red óptica, es la amplitud de interacción (de contacto) y es el potencial químico. La aproximación de enlace fuerte da como resultado la sustitución , que conduce al hamiltoniano de Bose-Hubbard la física está restringida a la banda más baja ( ) y las interacciones son locales al nivel del modo discreto. Matemáticamente, esto se puede expresar como el requisito de que excepto para el caso . Aquí, es una función de Wannier para una partícula en un potencial de red óptica localizado alrededor del sitio de la red y para la banda de Bloch . [15]
Sutilezas y aproximaciones
La aproximación de enlace fuerte simplifica significativamente el segundo hamiltoniano cuantificado, aunque introduce varias limitaciones al mismo tiempo:
- En el caso de estados de un solo sitio con varias partículas en un solo estado, las interacciones pueden acoplarse a bandas de Bloch más altas, lo que contradice los supuestos básicos. Aun así, un modelo de banda única puede abordar la física de baja energía de un entorno de este tipo, pero con parámetros U y J que se vuelven dependientes de la densidad. En lugar de un parámetro U, la energía de interacción de n partículas puede describirse como cercana, pero no igual a U. [15]
- Al considerar la dinámica reticular (rápida), se agregan términos adicionales al hamiltoniano para que la ecuación de Schrödinger dependiente del tiempo se cumpla en la base de la función de Wannier (dependiente del tiempo). Los términos provienen de la dependencia temporal de las funciones de Wannier. [16] [17] De lo contrario, la dinámica reticular se puede incorporar haciendo que los parámetros clave del modelo dependan del tiempo, variando con el valor instantáneo del potencial óptico.
Resultados experimentales
Greiner et al. observaron experimentalmente transiciones de fase cuántica en el modelo de Bose-Hubbard [9] y el grupo de Immanuel Bloch observó parámetros de interacción dependientes de la densidad [18] . Desde 2009 es posible obtener imágenes con resolución atómica del modelo de Bose-Hubbard utilizando microscopios de gas cuántico [19] [20] [21]
Otras aplicaciones
El modelo de Bose-Hubbard es de interés en el campo de la computación cuántica y la información cuántica. Mediante este modelo se puede estudiar el entrelazamiento de átomos ultrafríos. [22]
Simulación numérica
En el cálculo de estados de baja energía, el término proporcional a significa que es improbable que haya una gran ocupación de un único sitio, lo que permite truncar el espacio de Hilbert local a estados que contengan partículas como máximo. Entonces, la dimensión del espacio de Hilbert local es La dimensión del espacio de Hilbert completo crece exponencialmente con el número de sitios de la red, lo que limita las simulaciones exactas por computadora de todo el espacio de Hilbert a sistemas de 15 a 20 partículas en 15 a 20 sitios de la red. [ cita requerida ] Los sistemas experimentales contienen varios millones de sitios, con un llenado promedio por encima de la unidad. [ cita requerida ]
Las redes unidimensionales pueden estudiarse utilizando el grupo de renormalización de la matriz de densidad (DMRG) y técnicas relacionadas, como la decimación de bloques con evolución temporal (TEBD). Esto incluye el cálculo del estado fundamental del hamiltoniano para sistemas de miles de partículas en miles de sitios de la red y la simulación de su dinámica gobernada por la ecuación de Schrödinger dependiente del tiempo . Recientemente, [ ¿cuándo? ] se han estudiado redes bidimensionales utilizando estados de pares entrelazados proyectados , una generalización de los estados de productos matriciales en dimensiones superiores, tanto para el estado fundamental [23] como para la temperatura finita. [24]
Las dimensiones más altas son significativamente más difíciles debido al rápido crecimiento del entrelazamiento . [25]
Todas las dimensiones pueden ser tratadas mediante algoritmos cuánticos de Monte Carlo , [ cita requerida ] que proporcionan una forma de estudiar las propiedades de los estados térmicos del hamiltoniano, y en particular el estado fundamental.
Generalizaciones
Se pueden derivar hamiltonianos de tipo Bose-Hubbard para diferentes sistemas físicos que contienen gas atómico ultrafrío en el potencial periódico. Entre ellos se incluyen:
- sistemas con interacciones densidad-densidad de mayor alcance de la forma , que pueden estabilizar una fase supersólida para ciertos valores de parámetros
- imanes dimerizados, donde los electrones de espín 1/2 están unidos en pares llamados dímeros que tienen estadísticas de excitación bosónica y se describen mediante un modelo de Bose-Hubbard
- interacción dipolar de largo alcance [26]
- sistemas con términos de tunelización inducidos por interacción [27]
- estructura de espín interna de los átomos, por ejemplo debido a la captura de una variedad degenerada completa de estados de espín hiperfinos (para F=1 conduce al modelo de espín 1 de Bose-Hubbard) [28] [ aclaración necesaria ]
- situaciones en las que el gas experimenta un potencial adicional, por ejemplo, en sistemas desordenados. [29] El desorden podría realizarse mediante un patrón de moteado o utilizando una segunda red óptica inconmensurable y más débil. En este último caso, la inclusión del desorden equivale a incluir un término adicional de la forma: .
Véase también
Referencias
- ^ Gersch, H.; Knollman, G. (1963). "Modelo de célula cuántica para bosones". Physical Review . 129 (2): 959. Bibcode :1963PhRv..129..959G. doi :10.1103/PhysRev.129.959.
- ^ Ma, M.; Halperin, BI; Lee, PA (1986-09-01). "Superfluidos fuertemente desordenados: fluctuaciones cuánticas y comportamiento crítico". Physical Review B . 34 (5): 3136–3143. Bibcode :1986PhRvB..34.3136M. doi :10.1103/PhysRevB.34.3136. PMID 9940047.
- ^ Giamarchi, T.; Schulz, HJ (1988-01-01). "Localización de Anderson e interacciones en metales unidimensionales". Physical Review B . 37 (1): 325–340. Bibcode :1988PhRvB..37..325G. doi :10.1103/PhysRevB.37.325. PMID 9943580.
- ^ abcd Fisher, Matthew PA; Grinstein, G.; Fisher, Daniel S. (1989). "Localización de bosones y la transición superfluido-aislante" (PDF) . Physical Review B . 40 (1): 546–70. Bibcode :1989PhRvB..40..546F. doi :10.1103/PhysRevB.40.546. PMID 9990946.,
- ^ ab Jaksch, D.; Zoller, P. (2005). "La caja de herramientas de Hubbard del átomo frío". Anales de Física . 315 (1): 52. arXiv : cond-mat/0410614 . Código Bibliográfico :2005AnPhy.315...52J. CiteSeerX 10.1.1.305.9031 . doi :10.1016/j.aop.2004.09.010. S2CID 12352119.
- ^ Giamarchi, Thierry; Rüegg, Christian; Tchernyshyov, Oleg (2008). "Condensación de Bose-Einstein en aislantes magnéticos". Nature Physics . 4 (3): 198–204. arXiv : 0712.2250 . Código Bibliográfico :2008NatPh...4..198G. doi :10.1038/nphys893. S2CID 118661914.
- ^ Zapf, Vivien; Jaime, Marcelo; Batista, CD (15 de mayo de 2014). "Condensación de Bose-Einstein en imanes cuánticos". Reseñas de Física Moderna . 86 (2): 563–614. Bibcode :2014RvMP...86..563Z. doi :10.1103/RevModPhys.86.563.
- ^ Kühner, T.; Monien, H. (1998). "Fases del modelo unidimensional de Bose-Hubbard". Physical Review B . 58 (22): R14741. arXiv : cond-mat/9712307 . Código Bibliográfico :1998PhRvB..5814741K. doi :10.1103/PhysRevB.58.R14741. S2CID 118911555.
- ^ ab Greiner, Markus; Mandel, Olaf; Esslinger, Tilman; Hänsch, Theodor W.; Bloch, Immanuel (2002). "Transición de fase cuántica de un superfluido a un aislante de Mott en un gas de átomos ultrafríos". Naturaleza . 415 (6867): 39–44. Código Bib :2002Natur.415...39G. doi :10.1038/415039a. PMID 11780110. S2CID 4411344.
- ^ Morrison, S.; Kantian, A.; Daley, AJ; Katzgraber, HG; Lewenstein, M.; Büchler, HP; Zoller, P. (2008). "Réplicas físicas y el vidrio de Bose en gases atómicos fríos". New Journal of Physics . 10 (7): 073032. arXiv : 0805.0488 . Bibcode :2008NJPh...10g3032M. doi :10.1088/1367-2630/10/7/073032. ISSN 1367-2630. S2CID 2822103.
- ^ Thomson, SJ; Walker, LS; Harte, TL; Bruce, GD (3 de noviembre de 2016). "Medición del parámetro de orden de Edwards-Anderson del vidrio de Bose: un enfoque de microscopio de gas cuántico". Physical Review A . 94 (5): 051601. arXiv : 1607.05254 . Código Bibliográfico :2016PhRvA..94e1601T. doi :10.1103/PhysRevA.94.051601. S2CID 56027819.
- ^ Thomson, SJ; Krüger, F. (2014). "Ruptura de simetría de réplica en el vidrio Bose". EPL . 108 (3): 30002. arXiv : 1312.0515 . Bibcode :2014EL....10830002T. doi :10.1209/0295-5075/108/30002. S2CID 56307253.
- ^ Sachdev, Subir (2011). Transiciones de fase cuánticas. Cambridge University Press. ISBN 9780521514682.OCLC 693207153 .
- ^ Jaksch, D.; Bruder, C.; Cirac, J.; Gardiner, C.; Zoller, P. (1998). "Átomos bosónicos fríos en redes ópticas". Physical Review Letters . 81 (15): 3108. arXiv : cond-mat/9805329 . Código Bibliográfico :1998PhRvL..81.3108J. doi :10.1103/PhysRevLett.81.3108. S2CID 55578669.
- ^ ab Lühmann, DSR; Jürgensen, O.; Sengstock, K. (2012). "Tunelización multiorbital e inducida por densidad de bosones en redes ópticas". New Journal of Physics . 14 (3): 033021. arXiv : 1108.3013 . Código Bibliográfico :2012NJPh...14c3021L. doi :10.1088/1367-2630/14/3/033021. S2CID 119216031.
- ^ Sakmann, K.; Streltsov, AI; Alon, OE; Cederbaum, LS (2011). "Modelos reticulares dependientes del tiempo óptimos para dinámicas de no equilibrio". New Journal of Physics . 13 (4): 043003. arXiv : 1006.3530 . Bibcode :2011NJPh...13d3003S. doi :10.1088/1367-2630/13/4/043003. S2CID 118591567.
- ^ Łącki, M.; Zakrzewski, J. (2013). "Dinámica rápida para átomos en redes ópticas". Physical Review Letters . 110 (6): 065301. arXiv : 1210.7957 . Código Bibliográfico :2013PhRvL.110f5301L. doi :10.1103/PhysRevLett.110.065301. PMID 23432268. S2CID 29095052.
- ^ Will, S.; Best, T.; Schneider, U.; Hackermüller, L.; Lühmann, DSR; Bloch, I. (2010). "Observación resuelta en el tiempo de interacciones coherentes de múltiples cuerpos en reactivaciones de fase cuántica". Nature . 465 (7295): 197–201. Bibcode :2010Natur.465..197W. doi :10.1038/nature09036. PMID 20463733. S2CID 4382706.
- ^ Bakr, Waseem S.; Gillen, Jonathon I.; Peng, Amy; Fölling, Simon; Greiner, Markus (2009). "Un microscopio cuántico de gas para detectar átomos individuales en una red óptica de régimen Hubbard". Nature . 462 (7269): 74–77. arXiv : 0908.0174 . Bibcode :2009Natur.462...74B. doi :10.1038/nature08482. PMID 19890326. S2CID 4419426.
- ^ Bakr, WS; Peng, A.; Tai, ME; Ma, R.; Simon, J.; Gillen, JI; Fölling, S.; Pollet, L.; Greiner, M. (30 de julio de 2010). "Investigación de la transición de superfluido a aislante de Mott a nivel de átomo único". Science . 329 (5991): 547–550. arXiv : 1006.0754 . Bibcode :2010Sci...329..547B. doi :10.1126/science.1192368. ISSN 0036-8075. PMID 20558666. S2CID 3778258.
- ^ Weitenberg, Christof; Endres, Manuel; Sherson, Jacob F.; Cheneau, Marc; Schauß, Peter; Fukuhara, Takeshi; Bloch, Emmanuel; Kuhr, Stefan (2011). "Direccionamiento de un solo giro en un aislante atómico Mott". Naturaleza . 471 (7338): 319–324. arXiv : 1101.2076 . Código Bib :2011Natur.471..319W. doi : 10.1038/naturaleza09827. PMID 21412333. S2CID 4352129.
- ^ Romero-Isart, O; Eckert, K; Rodó, C; Sanpera, A (2007). "Generación de transporte y entrelazamiento en el modelo Bose-Hubbard". Revista de Física A: Matemática y Teórica . 40 (28): 8019–31. arXiv : quant-ph/0703177 . Código Bib : 2007JPhA...40.8019R. doi :10.1088/1751-8113/40/28/S11. S2CID 11673450.
- ^ Jordan, J; Orus, R; Vidal, G (2009). "Estudio numérico del modelo de núcleo duro de Bose-Hubbard en una red cuadrada infinita". Phys. Rev. B . 79 (17): 174515. arXiv : 0901.0420 . Código Bibliográfico :2009PhRvB..79q4515J. doi :10.1103/PhysRevB.79.174515. S2CID 119073171.
- ^ Kshetrimayum, A.; Rizzi, M.; Eisert, J.; Orus, R. (2019). "Algoritmo de recocido de redes tensoras para estados térmicos bidimensionales". Phys. Rev. Lett . 122 (7): 070502. arXiv : 1809.08258 . Código Bibliográfico :2019PhRvL.122g0502K. doi :10.1103/PhysRevLett.122.070502. PMID 30848636. S2CID 53125536.
- ^ Eisert, J.; Cramer, M.; Plenio, MB (2010). "Coloquio: Leyes de área para la entropía del entrelazamiento". Reseñas de Física Moderna . 82 (1): 277. arXiv : 0808.3773 . Bibcode :2010RvMP...82..277E. doi :10.1103/RevModPhys.82.277.
- ^ Góral, K.; Santos, L.; Lewenstein, M. (2002). "Fases cuánticas de bosones dipolares en redes ópticas". Physical Review Letters . 88 (17): 170406. arXiv : cond-mat/0112363 . Código Bibliográfico :2002PhRvL..88q0406G. doi :10.1103/PhysRevLett.88.170406. PMID 12005738. S2CID 41827359.
- ^ Sowiński, T.; Dutta, O.; Hauke, P.; Tagliacozzo, L.; Lewenstein, M. (2012). "Moléculas dipolares en redes ópticas". Cartas de revisión física . 108 (11): 115301. arXiv : 1109.4782 . Código Bib : 2012PhRvL.108k5301S. doi : 10.1103/PhysRevLett.108.115301. PMID 22540482. S2CID 5438190.
- ^ Tsuchiya, S.; Kurihara, S.; Kimura, T. (2004). "Transición superfluida-aislante de Mott de bosones de espín 1 en una red óptica". Physical Review A . 70 (4): 043628. arXiv : cond-mat/0209676 . Código Bibliográfico :2004PhRvA..70d3628T. doi :10.1103/PhysRevA.70.043628. S2CID 118566913.
- ^ Gurarie, V.; Pollet, L.; Prokof'Ev, NV; Svistunov, BV; Troyer, M. (2009). "Diagrama de fases del modelo desordenado de Bose-Hubbard". Physical Review B . 80 (21): 214519. arXiv : 0909.4593 . Código Bibliográfico :2009PhRvB..80u4519G. doi :10.1103/PhysRevB.80.214519. S2CID 54075502.