Cálculo de la concurrencia para el modelo de Heisenberg

Calculation of concurrence to the Heisenberg model

Resúmenes

La concurrencia es una cantidad que nos permite medir el grado de entreveramiento que presenta un sistema cuántico y se puede calcular a partir de la matriz densidad reducida. En este artículo mostramos explicitamente como calcular la concurrencia para una cadena finita de espines s =1/2 descrita por el modelo de Heisenberg anistrópico. Nosotros mostramos que para cadenas finitas la concurrencia tiene un máximo en el punto crítico Δ = 1, la cual es una de las principales características en el límite termodinámico. Se mostró que es posible obtener información sobre las transiciones de fase del modelo calculando la concurrencia para cadenas de por lo menos 12 sitios

concurrencia; entreveramiento cuántico; modelo de Heisenberg


The concurrence is a quantity that allows us to measure the entanglement of a quantum system and can be calculated from the reduced density matrix. In this paper we show explicitly how to calculate the concurrence for a finite chain of spins s =1/2 described by the anisotropic Heisenberg model. We show that for finite chains the concurrence has a maximum at the critical point Δ = 1, which is the main feature in the thermodynamic limit. We observe that it is possible to obtain information about quantum phase transitions of the model by calculating the concurrence for chains of at least 12 sites

concurrence; quantum entanglement; Heisenberg model


ARTIGOS GERAIS

Cálculo de la concurrencia para el modelo de Heisenberg

Calculation of concurrence to the Heisenberg model

R. Castellanos; R. Franco; J. Silva-Valencia1 1 E-mail: jsilvav@unal.edu.co.

Departamento de Física, Universidad Nacional de Colombia, Bogotá, Colombia

RESUMEN

La concurrencia es una cantidad que nos permite medir el grado de entreveramiento que presenta un sistema cuántico y se puede calcular a partir de la matriz densidad reducida. En este artículo mostramos explicitamente como calcular la concurrencia para una cadena finita de espines s =1/2 descrita por el modelo de Heisenberg anistrópico. Nosotros mostramos que para cadenas finitas la concurrencia tiene un máximo en el punto crítico Δ = 1, la cual es una de las principales características en el límite termodinámico. Se mostró que es posible obtener información sobre las transiciones de fase del modelo calculando la concurrencia para cadenas de por lo menos 12 sitios.

Palavras-chave: concurrencia, entreveramiento cuántico, modelo de Heisenberg.

ABSTRACT

The concurrence is a quantity that allows us to measure the entanglement of a quantum system and can be calculated from the reduced density matrix. In this paper we show explicitly how to calculate the concurrence for a finite chain of spins s =1/2 described by the anisotropic Heisenberg model. We show that for finite chains the concurrence has a maximum at the critical point Δ = 1, which is the main feature in the thermodynamic limit. We observe that it is possible to obtain information about quantum phase transitions of the model by calculating the concurrence for chains of at least 12 sites.

Keywords: concurrence, quantum entanglement, Heisenberg model.

1. Introducción

El entreveramieto cuántico (EC) es la característica más intrigante de la teoría cuántica, el sorprendente comportamiento de los estados entreverados había sido acentuado por Einstein, Podolsky y Rosen en 1931 [1] quienes consideraron un sistema de dos partículas en un estado entreverado y mostraron que las suposiciones de localidad, realismo y completés de la mecánica cuántica fallan al dar una descripción física consistente de los sistemas cuánticos. Varios años después, Bell encontró sus famosas desigualdades [2] con las cuales mostró que el entreveramiento no puede ser consistente con ninguna teoría local que contenga variables ocultas. Desde entonces el entreveramiento se ha vuelto un recurso de mucho valor en la teoría de información y comunicación cuántica.

Debido al isomorfismo existente entre los espacios de Hilbert de un espín (|↓), |↑) y los de un qubit (|0), |1), se ha considerado que los sistemas de espines serian la implementación natural de los dispositivos de la teoría cuántica de información [3].

Los progresos experimentales en la síntesis de nuevos materiales han permitido observar que las propiedades magnéticas de muchos materiales reales pueden ser descritas mediante cadenas o escaleras de espines debilmente interactuantes, cuyas propiedades físicas pueden ser explicadas usando el modelo de Heisenberg o generalizaciones de este [4]. Estas cadenas o escaleras de espines presentan propiedades inusuales debidas a la presencia de fuertes fluctuaciones cuánticas y efectos topológicos, ejemplos de estas propiedades inusuales son los valles en las curvas de magnetización y las transiciones de fase cuánticas a temperatura cero que están caracterizadas por el cambio en las propiedades del estado base causada por modificaciones en los parámetros de interacción presentes en el Hamiltoniano.

El entreveramiento es una cualidad de la mecánica cuántica, que corresponde a la presencia de correlaciones no locales entre diferentes partes de un sistema que no pueden ser explicadas de manera clásica. Un estado puro de un par de sistemas cuánticos se llama entreverado si no es factorizable (es decir, si el estado total no puede escribirse como producto de estados de una partícula), y un estado mezclado es entreverado si no se puede representar como la mezcla de estados puros factorizables [5].

Dada la importancia del concepto de entreveramiento en la teoría cuántica de la información (QIT) [6-8], la atención de este campo ha sido puesta, en los últimos años, en los sistemas estudiados en la materia condensada, donde es posible observar estados entreverados. Además se ha encontrado un especial vínculo entre el entreveramiento cuántico y las transiciones de fase cuánticas, las cuales describen cambios en el estado fundamental de un sistema a temperatura cero (T = 0), cuando se varía algún parámetro en el Hamiltoniano del sistema, como el campo magnético, el parámetro de anisotropía u otros. Wu [9] mostró de forma general que las transiciones de fase cuánticas son señaladas por medio de una discontinuidad en alguna medida del entreveramiento en un sistema cuántico. Algunas de las medidas del entreveramiento son por ejemplo la concurrencia y la entropía de von Neumann. Aunque el entreveramiento es una nueva herramienta para estudiar las transiciones de fase cuánticas en diversos modelos, aún no se sabe por que algunas medidas del entreveramiento son utiles y otras no.

Wootters and Hill [5, 10] definieron la concurrencia de un sistema bipartita de la siguiente forma

donde λi son los autovalores en orden decreciente de la matriz

donde ρr es la matriz densidad reducida del sistema considerado y σy es la matriz de Pauli usual en 2 dimensiones. Varios autores han encontrado que la concurrencia y/o otras medidas del entreveramiento son herramientas utiles para identificar transiciones de fase en diversos sistemas correlacionados ya que estas cantidades presentan valores extremos y/o singularidades en los puntos críticos [11-16]. Para los sistemas correlacionados de muchos cuerpos, cómo cuantificar apropiadamente el entreveramiento y cuál es su conexion intrínseca con las transiciones de fase cuánticas hansido preguntas tanto desafiantes como fundamentales en la ultima decada. Proximamente se espera que el entreveramiento cuántico juege un papel importante en la determinación de las transiciones de fase cuánticas, ya que ésto aportará a un conocimiento más profundo de éste tipo de fenómenos.

En este artículo calculamos de manera explicita la concurrencia para cadenas finitas de espines descritas por el modelo de Heisenberg unidimensional. Al mostrar como se calcula la concurrencia para este modelo se contribuye con la formación de los estudiantes de física de pregrado y postgrado, ya que el concepto de entreveramiento y sus medidas esta siendo ampliamente abordado en los cursos de mecánica cuántica y ya es incluido en algunos textos [17].

El artículo esta organizado de la siguiente manera, en la sección 2 mencionamos el modelo de Heisenberg, las diferentes fases que presenta y su diagrama de fases cuando un campo magnético externo es aplicado. En la sección 3 mencionamos las simetrias del modelo de Heisenberg, discutimos como determinar el Hamiltoniano para una cadena de espines, calculamos la matriz densidad reducida explicitamente para un sistema conformado por dos espines, calculamos la concurrencia y finalmente mostramos como se puede incrementar el tamaño del sistema. En las secciones 4 y 5 se presentan los resultados para la concurrencia en función del parámetro de anisotropía para redes de varios tamaños y las conclusiones del artículo respectivamente.

2. El modelo

Nosotros estudiaremos una cadena de espines s =1/2 interactuantes, descrita por el modelo de Heisenberg cuyo Hamiltoniano viene dado por

donde ,, denotan los operadores de espín y la suma es sobre todos los sitios de la cadena, la cual tiene una longitud L. J es conocida como la integral de intercambio y Δ es el parámetro de anisotropía, ambos son parámetros de interacción propios del modelo. El Hamiltoniano (3) también puede ser escrito en terminos de los operadores = + i y = - i así

El modelo de Heisenberg para espines s = ½ en una dimension(1D) tiene solución exacta usando el anstaz de Bethe. En la Fig. 1 presentamos el diagrama de fases para el modelo de Heisenberg cuando un campo magnético h es aplicado, vemos que el modelo presenta tres fases: ferromaganética, paramagnética y antiferromagnética. Se observa que el campo crítico para el cual se da la transición ferro-para viene dado por hc/J = Δ +1 y que en ausencia de campo magnético el modelo tiene una fase paramagnética para -1 < Δ < 1, una fase ferromagnética para Δ < - 1 y una fase antiferromagnética para delta > 1 (ver Fig. 1). En la fase paramagnética los espines apuntan aleatoriamente en todas las direcciones (ver Fig. 1(a)) conduciendo a magnetización total igual a cero. En la fase ferromagnética todos los espines apuntan en la misma dirección y la magnetización total por sitio es 1/2 (ver Fig. 1(b)). En la fase antiferromagnética los espines vecinos están ordenados antiparalelamente y la magnetización total es igual a cero(ver Fig. 1(c)).



3. Diagonalización exacta, matriz densidad y concurrencia

La expresión "diagonalización exacta" es usada para describir un número de diferentes aproximaciones las cuales conducen a resultados numéricamente exactos para sistemas de redes finitas, diagonalizando directamente la respresentación matricial del Hamiltoniano del sistema en una base apropiada de muchas partículas. El proceso más simple, pero a su vez el que más consume tiempo y memoria es la diagonalización completa de la matriz Hamiltoniana que permite calcular todas las propiedades deseadas del sistema. Sin embargo, la dimensión de la base para un sistema cuántico fuertemente interactuante crece exponencialmente con el tamaño del sistema, así que es posible tratar a lo más con unos pocos sitios.

A continuación mostramos cómo se calculan la matriz densidad reducida y el entreveramiento para una cadena finita de espines a través de la diagonalización exacta del Hamiltoniano del modelo. Vamos a consid erar un bloque con dos espines y analizaremos su entreveramiento con el resto de la cadena que actúa como un baño térmico. Aquí se usan condiciones de frontera períodicas sobre el sistema y por lo tanto este es invariante bajo traslaciones, es decir que podemos escoger cualquier par de vecinos más cercanos como el bloque con dos espines.

Es posible mostrar que el Hamiltoniano (3) conmuta con la componente z del operador espín total, = ∑i , es decir

[H, ] = 0

Esto significa que el Hamiltoniano puede ser dividido en sectores o bloques con valores fijos del operador y que además los autoestados de H son a la vez autoestados del operador .

La base del espacio de Hilbert de cada sitio de la red está compuesta por los siguientes kets |↓ |↑, es decir el espacio de Hilbert de cada sitio tiene dimensión 2, por lo tanto 2L será la dimensión del espacio de Hilbert para una red con L sitios. La representación matricial de los operadores Sz y S+ para cada sitio de la red, en la base escogida es

Sz = S+ = ,

donde la primera y segunda columna de las matrices están asignadas a los estados |↓ y |↑ respectivamente. Considerando lo anterior, el Hamiltoniano para una cadena de dos sitios está dado por la expresión

H = .

Debido a la simetría del modelo explicada anteriormente es posible sectorizar por espín en la base escogida, asignando por ejemplo, al sector de espín Sz = 0 las columnas y filas de los estados | ↓↑ y | ↑↓, así el Hamiltoniano de dos sitios puede ser organizado en subbloques con espin definido -2, 0 y 2 (ver (3.)), así

H =

note que las entradas de la matriz que no están en los subbloques con espin definido son iguales a cero, lo cual hace que la matriz sea esparsa y permite diagonalizarla mas facilmente.

Ahora mostraremos como calcular el Hamiltoniano, la matriz densidad y la concurrencia para una cadena de cuatro espines s = ½ antiferromagnética tipo Heisenberg con condiciones de frontera periódicas y considerando J = 1 y Δ = 1, es decir isotropica. La cadena de espines será dividida en dos bloques cada uno con dos espines inicialmente, como se aprecia en la Fig. 3. El espacio de Hilbert correspondiente a cada bloque tendrá dimensión 4 y su base será la base de ocupación |↓↓, |↓↑, |↑↓, |↑↑ . De esta manera podemos escribir (3) para el caso períodico como


donde Hi, y con i = 1; 2 es el Hamiltoniano, el operador de proyección de espín en el eje z y el operador escalera respectivamente para cada espacio compuesto por dos sitios. En la ecuación (5) "h.c" significa el conjugado hermitico.

En el Hamiltoniano (5) tenemos el producto directo entre operadores de dos espacios de Hilbert cada uno de dimensión 4, el resultado será un espacio de Hilbert de dimensión 16 y los operadores de este nuevo espacio serán matrices de 16 × 16 como es de esperarse ya que la dimensión aumenta de forma exponencial con el número de sitios. Los elementos de matriz de los operadores del nuevo espacio se calculan de la siguiente forma: si A (dimensión m = 4) es la representación matricial de un operador del espacio 1 y B (dimensión n = 4) es la representación matricial de un operador del espacio 2, entonces la matriz D que es el producto directo de A y B tiene una dimensión mn y elementos [18]

donde

Se puede notar también que el producto de un elemento de en el subbloque con espín 2, con un elemento de en el subbloque con espín -2, aparecera en el subbloque con espin 0 del Hamiltoniano del espacio de Hilbert resultante del producto directo entre los espacios 1 y 2. Esto se evidencia, por ejemplo, en el siguiente caso

|↑↓

|↓↑ = |↑↓, ↓↑,

es decir el espacio resultante tendrá igualmente una estructura de bloques con los sectores de espin siendo obtenido a partir de los sectores de espin de los subbloques de y .

De la solución exacta del modelo de Heisenberg obtenida por medio del ansantz de Bethe [19], se sabe que el estado fundamental se caracteriza por tener Sz = 0, por tanto al considerar el Hamiltoniano (3) y sus simetrias encontramos que no es necesario diagonalizar toda la matriz Hamiltoniana (dimension 16), solo necesitamos el sector con espin = 0, el cual tiene una dimensión de 6 para cadena con 4 espines. El Hamiltoniano para el sector = 0 para J = 1 y Δ=1 es

Al diagonalizar (8) obtenemos que el estado fundamental tiene una energía E0 = -2.0 y su autoestado es

donde los

ij que aparecen al lado derecho representan las entradas numéricas reales de .

Para calcular el entreveramiento cuántico es necesario calcular la matriz densidad reducida (ρr) para un bloque con dos espines, esta matriz se puede hallar a partir del estado fundamental del sistema por medio de la expresión

donde

ij son las entradas numéricas del estado con i (j) corriendo a través de los estados del bloque (baño térmico), considerando que [H, ] = 0 tenemos que [ρr,] = 0, por lo tanto la matriz densidad puede ser escrita de forma matricial como

ρr =

cuyos elementos de matriz pueden ser expresados en términos de las funciones de correlación de espín [13] de la forma

A continuación mostramos cómo calcular algunos elementos de la matriz densidad a partir de la ec. (9) haciendo uso de la ec. (10),

que significa que el elemento ρ↓↓ estará ubicado en la fila correspondiente al estado ↓↓| y en la columna del estado |↓↓ . También para un elemento fuera de la diagonal tendremos

De esta manera obtenemos la matriz densidad reducida del bloque con dos sitios que es

ρr =

Ahora usando la matriz densidad anterior y la ma- triz de Pauli σy vamos a calcular la matriz

y obtenemos

S =

,

diagonalizando esta matriz nosotros obtenemos los autovalores λ1 = 0:56249, λ2, λ3, λ4 = 0:006944; con los cuales se calcula concurrencia que viene dada por [13]

Este es valor de la concurrencia para Δ = 1 para una red de cuatro sitios, esta cantidad mide el entreveramiento entre un bloque de dos sitios y el resto del sistema que actúa como un baño térmico, éste tipo de entreveramiento es de tipo local porque exhibe las correlaciones entre un estado local y la otra parte del sistema. Al variar el parámetro de anisotropía, los autovalores de la matriz densidad cambian y por tanto la medida del entreveramiento C(S) cambia, este hecho nos permite estudiar las transiciones de fase en el modelo, determinando como cambia C(S) en función del parámetro anisogtropía Δ.

Hasta aqui hemos calculado la concurrencia para una red de cuatro sitios, sin embargo la mecánica cuántica nos dice que la energía del estado fundamental en un pozo de potencial infinito depende del tamaño, por tanto al aumentar el tamaño de la red los valores de energía, los autovalores de la matriz densidad y la concurrencia van a cambiar, este efecto es comunmente denominado efecto de tamaño finito. Por tanto es necesario incrementar el tamaño de la red para obtener resultados más reales, para esto es necesario definir nuevos operadores que contengan toda la información de la red anterior y denominarlos como los operadores del nuevo espacio 1, así

En la Fig. (4) se muestra la manera esquemática de realizar el aumento en la longitud del sistema. Por tanto el primer espacio queda con las propiedades de una cadena de cuatro espines con una dimensión de 16. Note que el operador se puede hallar simplemente como el adjunto del operador S+.


Usando rutinas de diagonalización estandard es posible repetir el procedimiento anterior hasta alcanzar aproximadamente una longitud de 12 sitios, sin embargo es posible alcanzar redes mayores usando algoritmos como el Lanczos o el Davidson para obtener el estado fundamental de matrices grandes [20, 21].

4. Resultados y Discusión

En la Fig. 5 se muestra el comportamiento de la concurrencia C en función del parámetro de anisotropía Δ para cadenas de espines de tamaños L = 4; 6 ; 12. Se observa que para L = 4 la concurrencia es cero hasta que Δ = -3.5, a partir de este valor es diferente de cero, note que un comportamineto similar se observa para L = 6 y L = 12, pero el valor de Δ para el cual deja de ser cero aumenta. En la Fig. 6(a) observa-mos la evolución del valor de Δ donde deja de ser cero la concurrencia en función de 1/L, en la figura vemos que efectivamente Δ tiende a -1 valor para el cual en ausencia de campo magnético tiene lugar la transición de la fase ferromagnética a la fase paramagnética (ver Fig. 2).




También observamos en la Fig. 5 que la concurrencia crece hasta alcanzar un máximo, el cual ocurre en Δ = 1 para todos los tamaños de redes finitas L = 4; 6 ; 12. A partir de este máximo la concurrencia decrece monotonamente en función del parámetro de anisotropía.

La curva continua en la Fig. 5 corresponde a la concurrencia exacta, la cual se puede obtener a partir de las funciones de correlación [6]

donde Gαα = Sα Sα , (α = x; y; z), estas funciones de correlación fueron calculadas usando el ansatz de Bethe [12] para una cadena de L = 1280 sitios. Vemos que la concurrencia deja de ser cero en Δ = -1 indicando la transición ferro-para en el modelo de Heisenberg, lo cual verifica la tendencia observada para redes finitas (ver Fig. 6(a)). Cuando Δ = 1 la energía del estado base es 1/4 - ln 2 [19] y se satisface que

con esto obtenemos que C = 0.3862, que es el valor máximo de la concurrencia.

Para una cadena de L = 4 sitios se mostro que la concurrencia para Δ = 1 es C = 0.499999, el cual es su valor máximo para esta cadena. En la Fig. 5 ve mos que a medida que el tamaño de la red aumenta la posición del máximo de la concurrencia no cambia, pero el valor de la concurrencia en este punto disminuye tendiendo hacia su valor en infinito C = 0.3862, la evolución del valor máximo de la concurrencia a medida que el tamaño de la cadena crece puede observarse en la Fig. 6(b).

Para una cadena de L = 12 se obtuvo que para anisotropías mayores a Δ = -1.01 la concurrencia es diferente de cero, que discrepa en 1.0% del valor en el límite termodinámico que es Δ = -1. El valor máximo de la concurrencia para esta cadena es C =0.398, para el cual obtuvimos un error de 3.06% al compararlo con el valor en el límite termodinámico. Observamos como para este modelo se puede obtener información sobre el entreveramiento y las transiciones de fase aún estudiando cadenas finitas de hasta L = 12 sitios.

Otra cantidad que nos permite estimar el entreveramiento cuántico en un sistema es la entropía de formación, la cual puede calcularse a partir de la concurrencia por medio de la expresión

donde x = 1/2+ √1 - C2/2.

En la Fig. 7 presentamos el comportamiento de la entropía de formación en función del parámetro de anisotropía Δ para redes finitas de tamaños L = 4, 6, 12. La curva continua corresponde a la entropía de formación exacta obtenida a partir de las funciones de correlación. En esta figura observamos que la entropía de formación tiene un máximo en Δ = 1 para todos los tamaños de red, lo cual indica la transición de la fase antiferromagnética a la fase paramagnética. Al igual que en el caso de la concurrencia se observa que a medida que crece el tamaño de la red, el valor del parámetro de anisotropía para el cual la entropía de formación deja de ser cero va tendiendo hacia Δ = -1 donde tiene lugar la transición de fase ferro-para.


5. Conclusiones

En este artículo hemos mostrado explicitamente como calcular la concurrencia, la cual es una de las principales medidas del entreveramiento cuántico. Este topico es de gran actualidad y ya se incluye en varios libros de texto de mecánica cuántica, por lo que mostrar como calcular medidas del entreveramiento como la concurrencia es fundamental y representa el principal aporte de este artículo.

El estudio del entreveramiento cuántico en cadenas de espines 1/2 descritas por el modelo de Heisenberg nos permite ganar información acerca de las transiciones de fase cuánticas presentadas por estos sistemas y compararlas con los resultados conocidos en el límite termodinámico. Nosotros mostramos que para cadenas finitas la concurrencia tiene un máximo en el punto crítico Δ = 1 indicando la transición antiferro-para. También se observo el punto crítico Δ = -1 donde la concurrencia deja de ser cero puede ser obtenido a partir de un analisis de efectos de tamaño finito.

Aunque los modelos estudiados ya tienen solución exacta, pequeñas variaciones en éstos conducen a modelos no solubles analíticamente, es allí donde un estudio numérico como el realizado sería fundamental para obtener información sobre los nuevos modelos.

Agradecimientos

Este trabajo fue financiado por la Division de Investigación Bogotá (DIB) de la Universidad Nacional de Colombia.

Recebido em 26/9/2009; Aceito em 26/1/2010; Publicado em 15/2/2011

  • 1
    E-mail:
    • [1] A. Einstein, B. Podolsky and N. Rosen, Physical Review 47, 777 (1935).
    • [2] J.S. Bell, Speakable and Unspeakable in Quantum Mechanics (Cambridge University Press, Cambridge, 1987).
    • [3] E. Plekhanov, A. Avella and F. Mancini, Physica B 403, 1282 (2008).
    • [4] D.C. Cabra, M.D. Grynberg, A. Honecker, P. Pujol, Condensed Matter Theories, vol. 16, Nova Science Publishers, New York, 2001.
    • [5] S. Hill and W. K. Wootters, Physical Review Letters 78, 5022 (1997)
    • [6] L. Amico, et al, Review of Modern Physics 80, 517 (2008).
    • [7] C T González-H, R Franco and J Silva-Valencia, Eur. J. Phys. 31, 681 (2010).
    • [8] J. J. Mendoza-Arenas, R. Franco, and J. Silva-Valencia, Phys. Rev. A 81, 062310 (2010).
    • [9] L. -A. Wu, et al, Physical Review Letters 93, 250404 (2004).
    • [10] W.K. Wootters, Physical Review Letters 80, 2245 (1998)
    • [11] Yan Chen et al., New Journal of Physics 8, 97 (2006).
    • [12] Shi-Jian Gu, Guang Shan-Tian and Hai-Qing Lin, Chin. Phys. Lett. 24, 2737 (2007).
    • [13] Shi-Jian Gu, H.Q. Lin and Y.Q. Li , New Journal of Physics 8, 61 (2006).
    • [14] Shi-Jian Gu, H.Q. Lin and Y.Q. Li Physical Review A 68, 042330 (2003).
    • [15] Shi-Jian Gu, G.S. Tian and H.Q. Lin Physical Review A 71, 052322 (2005).
    • [16] Shi-Jian Gu, Y.Q. Li and H.Q. Lin, Physical Review Letters 93, 086402 (2004).
    • [17] Stephen Gasiorowicz, Quantum Physics (John Wiley y Sons, Inc, New York, 2003).
    • [18] George B. Arfken and Hans J. Weber, Mathematical Methods for Physicists (Academic Press, New York, 1995).
    • [19] H.A. Bethe, Z. Physics 71, 205 (1931).
    • [20] E.R. Davidson. Computational Physics 7, 519 (1993).
    • [21] G.P. Zhang. Computational Physics Communications 109, 27 (1997).

    1 E-mail: jsilvav@unal.edu.co.

    Fechas de Publicación

    • Publicación en esta colección
      15 Mar 2011
    • Fecha del número
      Set 2010

    Histórico

    • Revisado
      26 Ene 2010
    • Recibido
      26 Set 2009
    • Acepto
      15 Feb 2011
    Sociedade Brasileira de Física Caixa Postal 66328, 05389-970 São Paulo SP - Brazil - São Paulo - SP - Brazil
    E-mail: marcio@sbfisica.org.br