Abdiel-Ecuacion de Dirac

Materiales de Dirac: Un vistazo al futuro
C. Abdiel de Jesús Espinosa Champo
Facultad de Ciencias en Fı́sica y Matemáticas

Universidad Autónoma de Chiapas
Tuxtla Gutiérrez, Chiapas; a 14 de marzo de 2019
C. Abdiel Espinosa Ch. (FCFM-UNACH) Materiales de Dirac Tuxtla Gtz., Chis. 06/03/19 1 / 52
Índice
1 Introducción Propiedades
Ecuación de Dirac Aplicaciones
Materiales de Dirac 3 Straintrónica y optotrónica
2 Grafeno Straintrónica en grafeno
Historia 4 Conclusiones
Introducción
Outline
Introducción Ecuación de Dirac
Reseña historica
A principios del siglo pasado aparecieron dos de las teorı́as fı́sicas que
cambiaron nuestra forma de comprender la naturaleza: la RG a grandes
escalas y la mecánica cuántica a pequeñas escalas.
Las ecuaciones que rigen a cada una de ellas son:
8πG
Gµν = Tµν (Ecuación de Einstein) (1)
c4
∂ h̄2 2
i h̄ Ψ(r , t) = − ∇ Ψ(r , t) + V (r ) (Ec. de Schrödinger) (2)
∂t 2m
Reseña historica
A principios del siglo pasado aparecieron dos de las teorı́as fı́sicas que
cambiaron nuestra forma de comprender la naturaleza: la RG a grandes
escalas y la mecánica cuántica a pequeñas escalas.
Las ecuaciones que rigen a cada una de ellas son:
8πG
Gµν = Tµν (Ecuación de Einstein) (1)
c4
∂ h̄2 2
Ψ(r , t) = −
i h̄ ∇ Ψ(r , t) + V (r ) (Ec. de Schrödinger) (2)
∂t 2m
¿Qué pasa con las partı́culas cuánticas relativistas?
Actualmente
Versión resumida
Figura 1: “Lagrangiano del modelo estándar”

1er intento: Klein-Gordon
En 1926, Oskar Klein y Walter Gordon propusieron la ecuación

h mc 2 i
+( ) Ψ(r , t) = 0 (3)
h̄
donde
1 ∂2
≡ − ∇2
c2 ∂ t2
que es invariante
p bajo las transformaciones de Lorentz, cumpliendo la
relación E = p 2 c 2 + m2 c 4 ; no obstante Ψ no puede ser interpretada en
relación a una densidad de probabilidad como en la Ec. de Schrödinger.
2o intento: Dirac
En 1928, Paul Dirac formuló una ecuación que describe las partı́culas
elementales de espı́n 12 , como el electrón y es completamente consistente
con los principios de la mecánica cuántica y de la teorı́a de la relatividad
especial.
(i h̄γ µ ∂µ − mc)ψ = 0 (4)
donde las γ µ satisfacen las relaciones {γ µ , γ ν } = 2η µν 1.
Si consideramos la energı́a
E 2 = p 2 c 2 + m2 c 4 (5)
en el regimen ultra-relativista cp mc 2 , la energı́a llega a ser lı́neal, y con

la elección apropiada de las matrices γ µ , la ecuación de Dirac se desacopla
en dos ecuaciones que pueden ser descritos como
i h̄∂0 ψR/L = ±c(σ · p̂)ψR/L (6)
donde p̂ = −i h̄∇, σ es el vector de matrices de Pauli, y ψR/L son espinores

con dos componentes. El Hamiltoniano,
H = ±c(σ · p̂) (7)
es conocido como el Hamiltoniano de Weyl, y los espinores asociados

como los espinores de Weyl.
Figura 2: Espectro para una partı́cula libre.
Entonces, ¿los modelos relativistas de Dirac y Weyl tienen un poco que ver
con lo que sucede con los electrones en un sólido?
Entonces, ¿los modelos relativistas de Dirac y Weyl tienen un poco que ver
con lo que sucede con los electrones en un sólido?
La respuesta es no siempre; esto es porque los sólidos presentan un
entorno complejo donde los electrones interactúan con los núcleos y entre
sı́, y estas interacciones modifican efectivamente las propiedades
electrónicas.
Consideremos dos átomos A y B; y los estados ψA ± ψB
a) b)
Figura 3: a) 5 átomos b) gran número de átomos.
Figura 4: a)semiconductor b)semimetal c) metal.
Introducción Materiales de Dirac
Materiales de Dirac/Weyl
Son materiales sólidos cristalinos cuya estructura de bandas cerca de la

energı́a de Fermi tienen relación de dispersión igual a las obtenidas a partir
de la ecuación de Dirac, además, el Hamiltoniano utilizado para describir a
los grados de libertad electrónicos es similar al Hamiltoniano Dirac-Weyl.
Materiales de Dirac/Weyl
Son materiales sólidos cristalinos cuya estructura de bandas cerca de la

energı́a de Fermi tienen relación de dispersión igual a las obtenidas a partir
de la ecuación de Dirac, además, el Hamiltoniano utilizado para describir a
los grados de libertad electrónicos es similar al Hamiltoniano Dirac-Weyl.
Figura 5: Materiales de Dirac hasta 2014. [T. O. Wehling et al. Dirac Materials,
Advances in Physics 63 (1), 2014.]
a)
b)
Figura 6: a) Vista lateral de siliceno b) Vista frontal de borofeno
Grafeno
Outline
Grafeno Historia
¿Qué es el grafeno?
El grafeno, fue el primer material verdaderamente bidimensional (2D) y cristalino

sintetizado por el hombre en el año 2004 [K. Novoselov et al. Science 306(5696), 2004].
a) b)
Figura 7: 1 a) Crecimiento de grafeno prı́stino sobre Cu(111). Imagen STM,
Vbias = +0.8 Volts, I = 0.8 nA, T = 80 K. b) Imagen STM estereográfica de una capa
de grafeno sobre Cu(111). I = 0.8 nA, Vbias = 0.8 Volts, T = 80 K. Rango de Z ∼ 0.06
nm. [George W. Flynn. The Journal of Chemical Physics 135(5), 2011].
Grafeno Historia
Grafeno Historia
Ignobel 2000
Figura 8: Por hacer flotar a una rana en un campo electromagnético de ∼ 16 T.
Grafeno Historia
¿De dónde y cómo se obtiene el grafeno?
La palabra grafeno hace referencia a una monocapa de carbono del grafito.

Estas monocapas pueden ser sintetizadas mediante la técnica conocida como
exfoliación micromecánica.
a) b)
Figura 9: 2 a) Grafito. b)Método de exfoliación micromecánica.
Grafeno Propiedades
Propiedades mecánicas
[Tim J. Booth et al. Nano Letters 8(8), 2008. ].
Se utilizó la técnica de
Microscopı́a de Fuerza Atómica
para realizar esfuerzos sobre
grafeno.
Se encontró que es 300 veces
más fuerte que el acero.
Es elástico hasta un 23 %
La tensión de ruptura excede 1
TPa.
Grafeno Propiedades
Propiedades mecánicas
[Tim J. Booth et al. Nano Letters 8(8), 2008. ].
Se utilizó la técnica de
Microscopı́a de Fuerza Atómica
para realizar esfuerzos sobre
grafeno.
Se encontró que es 300 veces
más fuerte que el acero.
Es elástico hasta un 23 %
La tensión de ruptura excede 1
TPa.
Grafeno Propiedades
Propiedades electrónicas y térmicas [Neil Savage. Nature, 483 (2012).]
Figura 10: Propiedades térmicas y eléctricas comparadas con otros materiales .
Grafeno Propiedades
La relación de dispersión del grafeno es lineal a bajas energı́as y forma un

cono (Cono de Dirac) alrededor del punto que coincide con el nivel de
Fermi del grafeno e implica que los electrones en grafeno se comportan
como fermiones de Dirac sin masa que se mueven a la velocidad de Fermi.
[K. S. Novoselov et al. Nature 438, 2005.].
Grafeno Aplicaciones
Circuitos integrados [Inanc Meric et al. Nature Nanotechnology 3(2008).]
El grafeno tiene alta movilidad

de portadores de carga, con
poco ruido, lo cual lo hace
adecuado para ser usado en la
fabricación de trnasistores.
Circuitos integrados basados en
grafeno manejan frecuencias
arriba de los 10 GHz.
Transistores impresos en
plásticos flexibles que pueden
operar a 25 GHz.
Membranas impermeables
Figura 12: El grafeno es una membrana impermeable a gases estándares, incluyendo el

Helio. J. Scott Bunch et al. Impermeble Atomic Membranes from Graphene Sheets. Nano
Lett. 8(8):2458-2462 (2008).
Nanomedicina
Figura 13: El grafeno es una material que promete mucho en cuanto a las aplicaciones
a nanomedicina.Hong Ying Mao et al. Graphene: Promises, Facts, Opportunities, and
challenges in Nanomedicine. Chemical Review 113:3407-3424 (2013).
Nanomedicina
Figura 14: Hay mucha investigación en cuanto a aplicaciones del grafeno en

nanomedicina, tales como: material antibacteial, biosensores, tratamiento de cáncer.
En 2010, la Academia China de Ciencias encontró que hojuelas de oxido de

grafeno son altamente efectivas para matar bacterias tales como la
Escherichia Coli.
Aplicaciones
Straintrónica y optotrónica
Outline
Straintrónica y optotrónica Straintrónica en grafeno
Consideremos una nanobanda con terminación armchair (AGN) unida a dos

conductores metálicos como se indica en la figura. Y aplicamos un campo de
deformación, u(xi , yi ) ≈ (u(xi ), 0) a lo largo de la dirección armchair.
Supondremos que los conductores metálicos no se ven afectados por dicho campo.
Figura 16: Representación esquemática de un AGN conectado a electrodos en las

interfases L (izquierda) y R (derecha).
Modelo
La ecuación de Schrödinger que debe ser resuelta ĤΨ(E ) = E Ψ(E ), donde

Ĥ es el Hamiltoniano de amarre fuerte a primeros vecinos definido por:
†
Ĥ = − ∑ tr0i,j +δi,n
0 c 0 c 0
r r
0 + ε0 ∑ cr†0 cr0i,j (8)
i,j i,j +δ i,n i,j
ri,j r0i,j
y Ψ(E ) es la función de onda del sistema para una energı́a E dada. En el

formalismo de Landauer-Büttiker la función de onda del sistema puede ser
representado como Ψ(E ) = ∑r0i,j αi,j |ψi,j i donde los coeficientes complejos
αi,j pueden ser determinados usando el método de matrices de
transferencia.
Matriz de transferencia
Matriz de transferencia
Reescribiendo la ecuación de Schrödinger podemos obtener la ecuación:
− ∑ tr0 0 αi +τ,j +δ = (E − ε0 )αi,j , (9)

i,j +δi,n
τ,γ
donde τ and γ especifican los sitios de los primeros vecinos de r0i,j . El siguiente
paso es expresar los M coeficientes αi,j para un i dado como un vector ~αi .
Reescribiendo la ecuación de Schrödinger podemos obtener la ecuación:
− ∑ tr0 0 αi +τ,j +δ = (E − ε0 )αi,j , (9)

i,j +δi,n
τ,γ
donde τ and γ especifican los sitios de los primeros vecinos de r0i,j . El siguiente
paso es expresar los M coeficientes αi,j para un i dado como un vector ~αi .
Ahora definamos la matriz de transferencia 2M × 2M,P̂i , que conecta α ~ i con las
de las columnas vecinas,de tal forma que la matriz total cumple:
 L   R 
~α0 ~αN
  = P̂  , (10)
L
~α1 R
~αN +1
Ahora supongamos que tenemos M ondas viajeras de la derecha y de la

izquierda con amplitud unitaria para la energı́a dada obteniendo que la
forma de las funciones de onda en la interface es :
s0 s0 0
αijL = ∑s 0 (δs 0 s e ikx xi + rs 0 s e −ikx xi ) sin(kys yj )
(11)
s 0 0
αijR = ∑s 0 ts 0 s e ikx xi sin(kys yj ),
donde rs 0 s y ts 0 s son los coeficientess de reflexión y transmisión del s-ésimo
al s 0 -ésimo canal, y δs 0 s es la delta de Kronecker. Esto permite reducir el
problema a un sistema de dos ecuaciones lineales para cada canal. El
conjunto de ecuaciones a resolver es:
    
−1 gs rss 1
  =  (12)
−ik s ik s
−e x hs tss e x
con    iNk s 
gs (N/4)+1 e x
  = ∏ ς̂r  , (13)
i(N+1)k s
hs r =0 e x
Obtenemos que la reflectancia, rss , de nuestro sistema es dado por

s
1 − (gs /hs )e ikx
|rss | = s
1 − (gs /hs )e −ikx
Del resultado previo, podemos obtener la conductividad (G) y el factor de

Fano (F), en el formalismo de Landauer-Büttiker [?], está dado por:
2e 2 M
G= |tss |2 , (14)
h ∑ s
Figura 17: Conductividad de un AGN deformado como función de E − ε0 para

grafeno prı́stino σ = 0 (color azul), y con deformación sinusoidal periódica con
σ = 5/7 (color rojo), y σ = (3 ∗ 89)/(4 ∗ 144) (color verde). Todas las curvas
fueron obtenidas para E = 0, λ = 0,04,
C. Abdiel Espinosa Ch. (FCFM-UNACH)
M = 1000, y N Tuxtla
Materiales de Dirac
= 100. Gtz., Chis. 06/03/19 41 / 52
Figura 18: Conductividad a bajas energı́as de un AGN deformado, como función

de E − ε0 para diferentes valores de longitud de onda de deformación (σ ). En
todos los casos hemos usado E = 0, λ = 0,08, M = 10000, y N = 1000.
Borofeno
Recientemente, existe un intenso interés en la sı́ntesis e investigación de
estructuras cristalinas 2D de Boro, conocidas como borofenos. El Boro es
un elemento fascinante debido a su complejidad quı́mica y estructural, y
los nanomateriales a base de boro de varias dimensiones han llamado
mucho la atención.
Figura 19: Lattice de borofeno 8-Pmmn
Iniciamos con el Hamiltoniano efectivo a bajas energı́as tipo Dirac
Ĥ = ρ(vx P̂x σx + vy P̂y σy + vt P̂y σ0 ) (15)
cuya relación de dispersión es

q
ρ
Eλ ,k = ρ h̄vt ky + λ h̄ vx2 kx2 + vy2 ky2 (16)
Figura 20: Gráfica de la energı́a E (k) como función de k, en la región alrededor

de un cono de Dirac.
Para borofeno irradiado por un campo electromagnetico perpendicular al

plano,
vt Π̂y vx Π̂x − ivy Π̂y
Ĥ = (17)
vx Π̂x + ivy Π̂y vt Π̂y
donde Π̂ = P̂ − (e/c)A y A el potencial vectorial dado

porA = EΩ0 cos(Gz − Ωt)(cos(θ ), sin(θ )).
Para campos intenso o longitud de onda larga (E0 /h̄Ω2 >> 300) la
ecuación para los espinores se transforman en una ecuación de Mathieu
d2
χ(φ ) + [a − 2q cos(2φ )] χ(φ ) = 0 (18)
dφ 2
donde 2
ζ
q=− (19)
2h̄Ω
y
ε 2
a= − 2q (20)
h̄Ω
Figura 21: Diagrama de estabilidad como función de los parámetros

adimensionales a y q para las soluciones de Mathieu. Las regiones de estabilidad
(región gris) y la de inestabilidad (región blanca) están divididas por las curvas
caracteristicas an (q) y bn (q) (lı́neas discontinuas).
Figura 22: Acercamiento del espectro de energı́a para el caso de longitud de onda
larga cerca de Υ = 0, i.e. ε = 0. Para un campo fijo (E0 = 2 V/m), la figura
muestra dos gaps de primer orden para θ = 0 y θ = π/2, respectivamente; y dos
bandas de energı́a. Observamos que el tamaño del gap decrece mientras θ
aumenta, con 0 ≤ θ ≤ π/2.
Conclusiones
Outline
Conclusiones
Conclusión
Iniciando con el descubrimiento del grafeno como el primer ejemplo

prominente, estos materiales han atraı́do un tremendo interés por sus
notables propiedades fundamentales, ası́ como su gran potencial y
viabilidad para aplicaciones, haciendo de los materiales de Dirac los
materiales del futuro.
Conclusiones
¡Gracias por su atención!

Abdiel-Ecuacion de Dirac

Cargado por

Copyright:

Formatos disponibles

Abdiel-Ecuacion de Dirac

Cargado por

Información del documento

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Abdiel-Ecuacion de Dirac

Cargado por

Copyright:

Formatos disponibles

Materiales de Dirac: Un vistazo al futuro

C. Abdiel de Jesús Espinosa Champo

Facultad de Ciencias en Fı́sica y Matemáticas

Tuxtla Gutiérrez, Chiapas; a 14 de marzo de 2019

Figura 1: “Lagrangiano del modelo estándar”

1er intento: Klein-Gordon

En 1926, Oskar Klein y Walter Gordon propusieron la ecuación

en el regimen ultra-relativista cp  mc 2 , la energı́a llega a ser lı́neal, y con

i h̄∂0 ψR/L = ±c(σ · p̂)ψR/L (6)

donde p̂ = −i h̄∇, σ es el vector de matrices de Pauli, y ψR/L son espinores

H = ±c(σ · p̂) (7)

es conocido como el Hamiltoniano de Weyl, y los espinores asociados

Figura 2: Espectro para una partı́cula libre.

Consideremos dos átomos A y B; y los estados ψA ± ψB

Figura 4: a)semiconductor b)semimetal c) metal.

Son materiales sólidos cristalinos cuya estructura de bandas cerca de la

Son materiales sólidos cristalinos cuya estructura de bandas cerca de la

El grafeno, fue el primer material verdaderamente bidimensional (2D) y cristalino

Figura 8: Por hacer flotar a una rana en un campo electromagnético de ∼ 16 T.

¿De dónde y cómo se obtiene el grafeno?

La palabra grafeno hace referencia a una monocapa de carbono del grafito.

Propiedades electrónicas y térmicas [Neil Savage. Nature, 483 (2012).]

Figura 10: Propiedades térmicas y eléctricas comparadas con otros materiales .

La relación de dispersión del grafeno es lineal a bajas energı́as y forma un

Circuitos integrados [Inanc Meric et al. Nature Nanotechnology 3(2008).]

El grafeno tiene alta movilidad

Figura 12: El grafeno es una membrana impermeable a gases estándares, incluyendo el

Figura 14: Hay mucha investigación en cuanto a aplicaciones del grafeno en

En 2010, la Academia China de Ciencias encontró que hojuelas de oxido de

Consideremos una nanobanda con terminación armchair (AGN) unida a dos

Figura 16: Representación esquemática de un AGN conectado a electrodos en las

La ecuación de Schrödinger que debe ser resuelta ĤΨ(E ) = E Ψ(E ), donde

y Ψ(E ) es la función de onda del sistema para una energı́a E dada. En el

Reescribiendo la ecuación de Schrödinger podemos obtener la ecuación:

− ∑ tr0 0 αi +τ,j +δ = (E − ε0 )αi,j , (9)

Reescribiendo la ecuación de Schrödinger podemos obtener la ecuación:

− ∑ tr0 0 αi +τ,j +δ = (E − ε0 )αi,j , (9)

Ahora supongamos que tenemos M ondas viajeras de la derecha y de la

Obtenemos que la reflectancia, rss , de nuestro sistema es dado por

Del resultado previo, podemos obtener la conductividad (G) y el factor de

Figura 17: Conductividad de un AGN deformado como función de E − ε0 para

Figura 18: Conductividad a bajas energı́as de un AGN deformado, como función

Figura 19: Lattice de borofeno 8-Pmmn

Iniciamos con el Hamiltoniano efectivo a bajas energı́as tipo Dirac

Ĥ = ρ(vx P̂x σx + vy P̂y σy + vt P̂y σ0 ) (15)

cuya relación de dispersión es

Figura 20: Gráfica de la energı́a E (k) como función de k, en la región alrededor

Para borofeno irradiado por un campo electromagnetico perpendicular al

donde Π̂ = P̂ − (e/c)A y A el potencial vectorial dado

Figura 21: Diagrama de estabilidad como función de los parámetros

Iniciando con el descubrimiento del grafeno como el primer ejemplo

¡Gracias por su atención!

También podría gustarte

en el regimen ultra-relativista cp mc 2 , la energı́a llega a ser lı́neal, y con