DESCRIPCIÓN TEÓRICA DEL DOPADO MOLECULAR DE MOLÉCULAS ORGÁNICAS CONJUGADAS USANDO LA TEORÍA DEL FUNCIONAL DENSIDAD (DFT) Lina Marcela Castellanos Daza Carlos Andrés Pinzón Osorio Universidad Distrital Francisco José de Caldas Facultad de ciencias y educación, Licenciatura en Física Bogotá, Colombia 2017
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
DESCRIPCIÓN TEÓRICA DEL DOPADO MOLECULAR DE MOLÉCULAS
ORGÁNICAS CONJUGADAS USANDO LA TEORÍA DEL FUNCIONAL DENSIDAD
(DFT)
Lina Marcela Castellanos Daza
Carlos Andrés Pinzón Osorio
Universidad Distrital Francisco José de Caldas
Facultad de ciencias y educación, Licenciatura en Física
Bogotá, Colombia
2017
II
Descripción teórica del dopado molecular de moléculas orgánicas
conjugadas usando la teoría del funcional densidad (DFT)
Lina Marcela Castellanos Daza
Carlos Andrés Pinzón Osorio
Informe de pasantía presentada como requisito parcial para optar al título de:
Licenciado en Física
Directores
Henry Alberto Méndez Pinzón PhD
Sol Milena Mejía Chica PhD
José Manuel Pérez Flórez Pérez PhD
Línea de investigación:
Semiconductores Orgánicos
Grupo de investigación:
Grupo de películas delgadas y nanofotónica
Grupo de Fotoquímica
Universidad Distrital Francisco José de Caldas
Facultad de ciencias y educación, Licenciatura en Física
Bogotá, Colombia
2017
III
Lista de Figuras
Figura 1. HOMO y LUMO de un semiconductor. .............................................................. 4
Figura 2. Ancho de banda de un metal (conductor), aislante y semiconductor................... 5
Figura 3. Estructura electrónica del tiofeno. ....................................................................... 5
Figura 4. Proceso para la obtención de átomos de carbono con hibridación 𝑠𝑝2. .............. 6
Figura 5. Orbitales de enlace molecular a partir de átomos de carbono con hibridación
Figura 17. Representaciones orbitales moleculares HOMO de los oligomeros y
LUMO de FxTCNQ. ......................................................................................................... 29
Figura 18. Representación del GAP calculado mediante el método OSC. ....................... 29
Figura 19. Comparación entre los niveles de energía de la molécula 4T+TCNQ. ........ 30
IV
Lista de Tablas
Tabla 1. Cronograma de actividades a realizar durante la pasantía. .................................. 3
Tabla 2. Cálculos de HOMO para 4T. ............................................................................. 25 Tabla 3. Calculo de GAP para el complejo conjugado 4T+F2TCNQ. ............................ 26
Tabla 4. Optimización y corrección de energía para 4T-TCNQ. ..................................... 27 Tabla 5. Energía del orbital más alto ocupado y del orbital más bajo desocupado. ........ 28 Tabla 6. Cálculos de GAP haciendo uso del método OSC. ............................................. 28
Tabla 7. Niveles de energia de los complejos conjugados. .............................................. 30
1.2 Resumen ............................................................................................................. 2 1.3 Cronograma de actividades ................................................................................ 2
CAPITULO II. Antecedentes
2.1 Aspectos generales de los polímeros conductores ............................................. 4 2.2 Estructura electrónica de los polímeros conjugados .......................................... 5 2.3 Dopaje y conducción electrónica en polímeros .................................................. 7 2.4 Tetracianoquinodimetano (TCNQ) como aceptor ............................................. 11
CAPITULO III. Metodología 3.1 Química computacional ..................................................................................... 13
3.1.1 Métodos teóricos ........................................................................................ 13 3.2 Teoría del funcional densidad (DFT) .................................................................. 14 3.3 Detalles computacionales ................................................................................. 17
3.3.1 Gausview ..................................................................................................... 17 3.3.2 Gaussian 09 ................................................................................................. 20 3.3.3 Optimización geométrica de las moléculas ............................................... 22 3.3.4 Niveles de energía HOMO y LUMO ............................................................ 23 3.3.5 Espectro RAMMAN UV/VIS e IR .................................................................. 23
CAPITULO IV. Resultados y Discusión
4.1 Optimización geométrica de moléculas dopantes/aceptoras y complejos de transferencia de carga .............................................................................................. 5 4.2 cálculo de niveles de energía de moléculas individuales y complejos de transferencia de carga ............................................................................................... 5 4.3 comparación de GAP entre dopantes/aceptores y complejos conjugados ........ 5
Conclusiones y recomendaciones
Bibliografía
A. Anexo
VI
Apéndice 1. Gap calculado con el método OSC
Apéndice 2. Espectros Raman, IR y UV-VIS de oligotiofenos y FxTCNQ
Apéndice 3. Participación en Tercera escuela de Química Computacional de las
Américas (PREFALC)
Apéndice 4. Participación en VI encuentro Nacional de Químicos teóricos y
computacionales
Apéndice 5. III Escuela Colombiana de teoría y computación en las ciencias
moleculares.
CAPITULO I . Introducción
1.1 Objetivos
1.1.1 objetivo general
Caracterizar moléculas semiconductoras orgánicas, con base en la combinación de cálculos de estructura electrónica DFT y la medición de propiedades optoelectrónicas efectuadas sobre oligotiofenos dopados con aceptores.
1.2.1 objetivos específicos
realizar la construcción de las estructuras moleculares de oligotiofenos y
tetracianoquinodimetanos fluorinados (FxTCNQ). desarrollar los cálculos de optimizaciones con el fin de obtener la geometría
basal de tanto los oligotiofenos desde n=1 a n=8 como de la familia de los
FxTCNQ. Obtener la energía basal de cada uno de los oligomeros y de los FxTCNQ a
través de una corrección con el fin de compararlos con los datos presentes en
la literatura.
2
1.2 Resumen
Los materiales orgánicos se han constituido en una alternativa valiosa para
aplicaciones en optoelectrónica. Estos materiales se clasifican entre moléculas
orgánicas conjugadas (COM) y polímeros conjugados (CP), los cuales han sido usados
como capas activas en dispositivos optoelectrónicas. Como el dopado de
semiconductores es un factor clave en el funcionamiento de los dispositivos
electrónicos modernos, se han hecho muchos esfuerzos desde la década pasada
para obtener capas orgánicas dopadas estables.
Dos modelos han surgido para explicar la baja eficiencia de la generación de carga
debida al dopado de semiconductores orgánicos: el modelo de transferencia entera
de carga (ICT) y el modelo de transferencia híbrida de carga (HCT). El primer modelo
concuerda con el dopado molecular adoptado de la física de semiconductores
inorgánicos: transferencia entera de carga desde el máximo orbital molecular
ocupado (HOMO) del semiconductor orgánico hacia el mínimo orbital molecular
desocupado (LUMO) de la molécula dopante tipo P, y viceversa para dopado tipo N.
Este proyecto, tratará desde una perspectiva teórica, apoyada por cálculos
computacionales, puntos relevantes como: formación de polarones, transferencia
interna de carga, generación de portadores libres, mediante un uso extenso del
cálculo de estructura electrónica DFT, del dopado de la familia de los tiofenos:
iniciando con monómeros, cuyo número se incrementa gradualmente (n=1, 2, 3,….,
8).Como moléculas dopantes aceptoras, se emplearán derivados fluorinados de la
molécula prototipo tetracianoquinodimetano (TCNQ). Este estudio suministrará
evidencia muy fuerte, soportada en cálculos de estructura electrónica basada en la
teoría del funcional densidad DFT, sobre la validez de los modelos que fundamentan
el mecanismo de dopado molecular: hibridización intermolecular ó transferencia
entera de carga, y proporcionará elementos que permitan conjugar ambas visiones
para lograr un modelo único que describa adecuadamente el proceso de generación
de carga inducido por el dopado molecular de semiconductores orgánicos.
3
1.3 Cronograma de actividades
ACTIVIDADES TIEMPO(H
ORAS)
Capacitación en fundamentos de mecánica cuántica 140
Capacitación en fundamentos de teoría de orbitales moleculares 120
Capacitación en fundamentos del sistema operativo Linux 10
Capacitación en manejo de software Gaussian 09 para cálculos de
estructura electrónica DFT
10
Capacitación en manejo de software Gaussview 09 para formulación y
visualización de resultados de cálculos de estructura
20
Capacitación en fundamentos de semiconductores orgánicos 100
Optimización de geometrías de moléculas conjugadas de la familia de
los oligotiofenos y de moléculas aceptoras fluorinadas derivadas del
(Tetracianoquinodimetano) TCNQ
80
Calculo de niveles de energía de orbitales de frontera y espectros
Raman, UV/VIS, e IR de oligotiofenos y aceptores
60
Optimización de geometrías de complejos de transferencia de carga
formados a partir de combinaciones donor/aceptor de
oligotiofeno/TCNQ
120
Calculo de niveles de energía de orbitales de frontera y espectros
Raman, UV/VIS e IR de complejos de transferencia de carga
60
Tiempo total planeado para el desarrollo de las actividades 720
Tabla 1. Cronograma de actividades a realizar durante la pasantía
4
CAPITULO II. Antecedentes
2.1 Aspectos generales de los polímeros conductores
Los polímeros conductores abarcan actualmente una amplia rama en la
investigación de la electrónica molecular, debido a que las características físicas y
químicas empleados para su construcción brindan grandes ventajas a nivel
tecnológico y electrónico, puesto que se pueden unificar algunas de las propiedades
de los polímeros y los metales, entre ellas “la resistencia mecánica, flexibilidad,
ligereza, alta conductividad eléctrica, resistencia a la corrosión, bajo coste de
fabricación y fácil preparación”, [1] optimizando de esta forma su uso y ampliando
el área de aplicaciones. Los primeros indicios que se tienen en este campo de
investigación se dieron a conocer con el premio nobel de química en el año 2000 por
AlanG. Heeger, Alan G. MacDiarmid y Hideki Shirakawa por el descubrimiento y
desarrollo de los polímeros conductores , quienes descubrieron que “al dopar una
película de poli acetileno con vapores de yodo, cloro o bromo la conductividad del
material aumentaba 1000 millones de veces.” [2]. Al igual que en los
semiconductores inorgánicos, el proceso de dopaje en las moléculas orgánicas
reduce la separación energética (Band Gap), necesaria para que se dé la conducción
eléctrica. “La estructura molecular de los polímeros (alternancia de enlaces simples
y dobles) tiene la propiedad fundamental de poseer orbitales electrónicos π
extendidos sobre toda la estructura, un electrón situado en uno de estos niveles
estaría muy deslocalizado y tendría una gran libertad de movimiento, de manera
que sería posible conseguir la conducción de electricidad.” [3] En estos sistemas
existe una banda de energía ocupada por electrones (HOMO) y una banda vacía que
contendrá a los electrones excitados (LUMO) véase Figura 1.
Figura 1. HOMO y LUMO del semiconductor
5
El valor de la separación energética entre HOMO y LUMO o ancho de banda (Band
Gap) define el tipo de polímero que se tiene: aislantes, conductores y
semiconductores. En la medida en que el ancho de banda tenga un valor de energía
menor, la conducción de los portadores de carga será mayor. Véase Figura 2.
Figura 2. Ancho de banda de un metal (conductor), aislante y semiconductor
2.2 Estructura electrónica en polímeros conjugados Los polímeros son cadenas carbonadas de una gran cantidad de monómeros a diferencia de los oligomeros que son cadenas finitas ( 2 a 13) de monómeros. El tiofeno hace parte de la familia de polímeros debido a que su estructura heterocíclica y aromática cumple con las características de los polímeros.
Figura 3. Estructura electrónica del tiofeno
6
La teoría de orbitales moleculares (OTM) y la hibridación están estrechamente relacionadas con la estructura electrónica de los polímeros conjugados, debido a que la conducción eléctrica en los polímeros se fundamenta en la OMT. “La estructura de doble enlace de carbono-carbono en los polímeros se fundamenta en una hibridación del tipo 𝑠𝑝2 entre los átomos de carbono. Según la OMT cuando “el orbital 2𝑠 se encuentra hibridado con dos de los orbitales 2𝑝, uno de los orbitales 2𝑝 queda sin hibridar, entonces se ubica un electrón en cada uno de los orbitales híbridos 𝑠𝑝2 y un electrón permanece en el orbital 2𝑝” Véase figura 4. [4]
Figura 4. Proceso para la obtención de átomos de carbono con hibridación 𝑠𝑝2 [4] “Desde el punto de vista de la OMT, el doble enlace carbono-carbono consiste de dos tipos de enlace, un enlace sigma σ y un enlace π. El enlace σ resulta del traslape de orbitales con hibridación sp2 y semeja a un eje enlazando los dos átomos de carbono. El enlace π resulta del traslape lateral de dos orbitales p; este tiene un plano nodal que semeja a un orbital p. En el estado base los electrones del enlace π están localizados entre los dos átomo de carbono, generalmente arriba y abajo del plano de la estructura de enlace σ”.[4] Véase figura 5.
Figura 5. Orbitales de enlace molecular a partir de átomos de carbono con
hibridación 𝑠𝑝2 [4]
7
Los electrones del enlace π tienen mayor energía que los electrones del enlace σ. Las energías relativas de los orbitales σ y π (con los electrones en estado basal) se muestran en la figura rr donde en enlace 𝜎∗ es el orbital de anti-enlace” [4] Figura rr
Figura 6. (a)Combinación de dos orbitales p para formar dos orbitales π. (b) energías relativas de los orbitales σ y π
En los polímeros conjugados, la hibridación 𝑠𝑝2𝑝𝑧 resulta, entonces en un electrón no apareado por cada átomo de carbono; en el eteno hay, en consecuencia, dos electrones sin aparear, o deslocalizados. Por tanto la conjugación en un polímero conduce, por tanto, a una estructura electrónica formada por múltiples estados discretos de energía, los cuales, al llegar a ser continuos, forman bandas” [5]. Cada banda tiene una capacidad finita de ser ocupada por electrones, así mismo la banda puede encontrase vacía. En los polímeros conjugados, los orbitales de enlace y anti enlace π generan bandas de energía, las que se encuentran ya sea totalmente ocupadas (banda π- π) o vacías (banda π-π ∗). La banda con mayor grado de ocupación se denomina HOMO ( Highest Occupied Molecular Orbital) y la banda con menor grado de ocupación se denomina LUMO (Lowest Unocuppied Molecular Orbital), estas bandas se encuentran separadas debido a la indeterminación de la función de onda en esta zona, esta separación energética se conoce como banda prohibida, siendo el valor de la diferencia energética existente entre las bandas HOMO y LUMO [6]
8
2.3 Dopaje y conducción electrónica en polímeros
La naturaleza de los polímeros es ser aislantes, pero debido a la adición de alguna sustancia se logra el dopado, la existencia de enlaces dobles y simples convirtiendo al polímero en un metal sintético. El metal sintético está compuesto por una cadena de carbonos que están enlazados mediante una sucesión intercalada de enlaces dobles y sencillos que de una manera muy eficaz permite la aparición de orbitales electrónicos con electrones que poseen una movilidad significativa. Véase figura 7
Figura 7. Estructura de tetra-tiofeno
La modificación de las propiedades del polímero mediante la introducción de “impurezas” se realiza de dos maneras distintas: el dopado tipo N que consiste en generar un exceso de electrones en la molécula anfitrión de tal manera que hay menos cantidad de huecos, este es menos usado debido a que no se puede generar en un ambiente regular abierto a la atmosfera debido a su característica oxidante por la alta concentración de oxígeno. Por otro lado está el dopado tipo P que genera un exceso de huecos sin que haya electrones asociados a este, esta característica facilita que sea el procedimiento más usado. Véase figura 8.
9
Figura 8. Representación del dopado tipo N (izquierda) y tipo P (derecha) [7]
Para explicar cómo es el proceso de conducción en los polímeros debemos acudir a la teoría de bandas en los materiales inorgánicos con el fin de realizar una analogía. La representación de la estructura electrónica de un átomo o molécula mediante bandas tomando como base el solapamiento de los orbitales de un átomo produciendo un numero discreto de orbitales moleculares es lo que se denomina teoría de bandas. La cual permite observar tres componentes que determinan la naturaleza del material: la banda de valencia en la cual van ubicados los electrones que se encargan de realizar los enlaces con otros átomos, la banda de conducción ocupada por los electrones libres que no están ligados a los átomos y que van a ser los encargados de realizar el proceso de conducción, y por último la banda prohibida o gap en la cual no pueden encontrarse electrones. Así pues tenemos que “cuando la banda de valencia esta solo parcialmente llena por los electrones permisibles, o si las dos bandas se superponen tal que no existe barrera energética entre ambas, la aplicación de un potencial permitirá el paso de algunos de los electrones hacia niveles vacantes donde estarán libres para desplazarse a través del sólido produciendo corriente eléctrica. Esto sería la descripción de un conductor. Cuando la banda de valencia está completa y separada de la banda de conducción vacante por una barrera energética, no se presenta el flujo de electrones bajo la influencia de un campo externo y para hacerlo se requiere una considerable cantidad energética. Tales materiales serán semiconductores o aisladores, dependiendo de la barrera energética. El modelo de banda asume que los electrones son deslocalizados y pueden extenderse sobre toda la estructura” [8] Ahora bien tomando los polímeros orgánicos y ya teniendo una aproximación a la
10
conducción eléctrica en los materiales inorgánicos, se debe hacer la aclaración de que no es posible para estos explicar el fenómeno anteriormente mencionado debido a que “los átomos son enlazados covalentemente entre sí, formando cadenas poliméricas que experimentan interacciones intermoleculares débiles. Luego la conducción macroscópica requerirá movimiento de los electrones, no solo a lo largo de la cadena sino de una cadena a otra. Esta teoría no explica el hecho de que los transportadores de cargas, generalmente electrones o huecos, en poliacetileno y polipirrol, no tienen espín” [8] pero la teoría de bandas ayudara a comprender a grandes rasgos lo que sucede en el polímero orgánico. Además para el fenómeno de la conducción eléctrica en polímeros se deben introducir ciertos conceptos, tales como soliton, polaron y bipolarones, que explican algunos fenómenos que ocurren dentro del material. “El movimiento de un electrón desde el máximo de la banda de valencia de un polímero conjugado, como el poliacetileno o polipirrol, origina una vacante (hueco o radical catión) que no se deslocaliza completamente como se esperaría de la teoría de banda clásica. Se produce una deslocalización parcial, extendiéndose sobre varias unidades monoméricas y causando una deformación estructural. El nivel de energía asociado con este catión radical representa un orbital no estable y tiene una mayor energía que las energías en la banda de valencia, es decir su nivel energético está en la barrera entre las bandas. Este radical catión está parcialmente deslocalizado sobre un segmento del polímero y es definido por la Física del Estado Sólido como polarón. Si otro electrón se desplaza desde el polímero con el polarón, pueden darse una posibilidad: -originarse otro polarón independiente en otro segmento del polímero o la salida desde el nivel del primer polarón (salida del electrón desapareado) dando lugar a un dicatión llamado por los físicos bipolarón”
Figura 9. Representación esquemática y diagramas energéticos de a) polaron y b) bipolaron [9]
De manera análoga como en los inorgánicos aparecen dos conceptos para definir la banda de valencia (HOMO) y la banda de conducción (LUMO) y la diferencia entre estos dos determinara el gap de la molécula, el dopado tipo N y tipo P respectivamente se caracteriza por que “los polímeros orgánicos en los que se realiza inyección de huecos, los niveles de Fermi de los electrodos poseen energías próximas al orbital HOMO, mientras que cuando se produce la inyección de electrones estos deben poseer energías próximas a la del orbital LUMO” [9] Para los polímeros es conveniente usar el dopado tipo P , debido a las facilidades que brinda para sus aplicación tecnológicas en ambientes abiertos a la atmosfera, de esta forma se tiene que estimar dos factores más que son importantes , la afinidad electrónica (EA) y el potencial de ionización (IP) . “El IP está relacionado con la energía necesaria para sustraer un electrón al sistema
12
y la EA está relacionada con la energía desprendida al adicionar un electrón al sistema. Por tanto, el IP está relacionado con la energía del orbital HOMO y la EA con el orbital LUMO. De este modo se suelen obtener valores teóricos que reproducen casi a la perfección los valores experimentales” [9].
2.4 TCNQ como aceptor
El dopado tipo P sobre polímeros orgánicos, en este caso el politiofeno se realiza usando el grupo de moléculas de la familia de TCNQ que hará de aceptor. Véase figura 10.
Figura 10. Familia de moléculas del FxTCNQ [10] Como aceptor dentro del complejo la molécula posee ciertas características que lo hacen ideal , tal como se describe a continuación “La molécula TCNQ es un sistema quinoide plano de closed-shell que tiene cuatro grupos de cianuro de alto electrón situados en las terminales de los átomos de carbono de metileno. Estas características explican su afinidad electrónica muy grande y el comportamiento asociado de tomar un electrón cuando se pone en contacto con virtualmente cualquier donante de electrones para formar el radical monociclo de TCNQ de open-shell. El único electrón no apareado en TCNQ ocupa la energía más baja vacía y se espera que resida principalmente localizada en los grupos terminales de diciano-metileno carbono con el fin de aprovechar al máximo la fuerte afinidad electrónica de los cianuros” [11]
13
Experimentalmente el valor de la densidad de carga local añadido por nitrógeno fue de 0.14 , lo que mostro que el exceso de electrones residía casi completamente sobre el componente C≡N mediante “14N nuclear quadrupolar resonance (NQR)” [12].
CAPITULO III. Metodología
3.1 Química Computacional
La química computacional simula estructuras y reacciones químicas numéricamente, usando las leyes fundamentales de la física. Algunos métodos no solo permiten explorar moléculas inestables, si no que se pueden estudiar estados de transición difíciles de aislar experimentalmente. Hay dos áreas principales en la química computacional orientadas al estudio de las moléculas y su reactividad, estas son la mecánica molecular (MM) y la teoría de estructura electrónica (EE). Ambas tienen en común los siguientes aspectos:
Cálculos de energía de una estructura particular
Optimización de geometría, encuentra el mínimo de energía para una estructura molecular cercana a la estructura de partida.
Cálculos de frecuencias vibracionales, resultado del movimiento interatómico en la molécula. [13]
3.1.1 Métodos teóricos
Mecánica molecular: los cálculos en MM no trata los electrones de un sistema molecular explícitamente, pero esta optimizada en el cálculo basado en la interacción electrón -núcleo. Los efectos electrónicos se encuentran implícitos en los campos de fuerza, que se componen de la suma de una serie de términos que modelizan las tensiones de los enlaces, de sus ángulos, sus torsiones, las interacciones de V. der W., las electrostáticas y las de los puentes de hidrogeno. [14]
Métodos de estructura electrónica: basado en las leyes de la mecánica cuántica, resolviendo la Ec. De Schrödinger se obtienen los estados cuánticos, energía y otras propiedades relacionadas. [15]
14
Los métodos de estructura electrónica se caracterizan por las aproximaciones matemáticas para resolver esta ecuación. Estos métodos pueden ser empíricos o ab-inito:
Métodos semi-empiricos: utilizan parámetros experimentales para simplificar el cálculo computacional.
Métodos ab-inito: a diferencia del MM no utilizan parámetros experimentales para el cálculo, están basados únicamente en las leyes de la mecánica cuántica.
Métodos de funcional de la densidad: se basa en el modelo molecular tomando en cuenta la correlación electrónica como funcional general de la densidad electrónica. [16]
3.2 Teoría del Funcional Densidad
El método DFT es una forma de calcular la correlación de electrones a través de funcionales generales de la densidad de electrones. La DFT funcionales dividen la energía electrónica en varios componentes que se calculan por separado: la energía cinética un término de correlación de intercambio que representa el resto de la interacción electrón-electrón (que se divide en componentes de intercambio y correlación en la mayoría de las formulaciones de la DFT reales)
Los métodos basados en la DFT derivan en última instancia de la investigación de la mecánica cuántica de 1920, especialmente del modelo Thomas-Fermi-Dirac y de la química cuántica del trabajo fundamental de Slater en 1950. El enfoque de la DFT se basa en una estrategia de modelar la correlación electrónica mediante de funcionales generales de la densidad electrónica.
Tales métodos deben sus orígenes modernos al teorema de Hohenberg-Kohn, publicado en 1964, que demostraba la existencia de una función única que determina exactamente la energía y la densidad del estado base. El teorema no proporciona la forma de este funcional, generalmente.
Siguiendo el trabajo de Kohn y Sham, los funcionales aproximados empleados por los métodos de la DFT actuales dividen la energía electrónica en varios términos:
𝐸 = 𝐸𝑇 + 𝐸𝑉 + 𝐸𝐽 + 𝐸𝑋𝐶 (1)
Donde 𝐸𝑇 es el término de energía cinética (derivado del movimiento de los electrones), 𝐸𝑉 incluye términos que describen la energía potencial de la atracción de electrones y de la repulsión entre pares de núcleos, 𝐸𝐽 es el término de repulsión electrón-electrón (es También descrita como la auto-interacción coulombiana de la
15
densidad electrónica) y 𝐸𝑋𝐶 es el término de correlación de intercambio e incluye la parte remanente de las interacciones electrón-electrón.
Todos los términos excepto la repulsión nuclear-nuclear son funciones de ρ, la densidad electrónica es dada por la siguiente expresión:
𝐸𝐽 =1
2∬ 𝜌(��1) (��12)−1𝜌(��2)𝑑��1𝑑��1 (2)
𝐸𝑇 + 𝐸𝑉 + 𝐸𝐽 Corresponde a la energía clásica de la distribución de la carga ρ. El término 𝐸𝑋𝐶 en la ecuación (2) explica los términos restantes en la energía:
La energía de intercambio que surge de la anti simetría de la función de onda de la mecánica cuántica.
Dinámica de correlación en los movimientos de los electrones individuales.
Hohenberg y Kohn demostraron que 𝐸𝑋𝐶 es determinado enteramente por (es un funcional de) la densidad de electrones. En la práctica, 𝐸𝑋𝐶 se suele aproximar como una integral que implica sólo las densidades de espín y posiblemente sus gradientes:
[8] Universidad del Pais Vasco. (n.d.). POLÍMEROS CONDUCTORES COMO ALTERNATIVA
TECNOLÓGICA EN LA PROTECCIÓN DE NUESTRO AMBIENTE. [Online] Available at:
http://POLÍMEROS CONDUCTORES COMO ALTERNATIVA TECNOLÓGICA EN LA PROTECCIÓN
DE NUESTRO AMBIENTE [Accessed 12 Jan. 2017].
[9] Marín Pérez, M. (2014). POLÍMEROS SEMICONDUCTORES -CONJUGADOS CON
UNIDADES ACEPTORAS Y DADORAS DE ELECTRONES DE INTERÉS EN CÉLULAS
FOTOVOLTAICAS. ESTUDIO TEÓRICO DE SU ESTRUCTURA MOLECULAR Y PROPIEDADES
OPTOELECTRÓNICAS. Pregrado. Universidad de Jaén, facultad de ciencias experimentales.
[10] Méndez, H. et al. Charge-transfer crystallites as molecular electrical dopants. Nat. Commun. 6:8560 doi: 10.1038/ncomms9560 (2015). [11] Lundqvist, B., Lundqvist, S. and Runnstrom-Reio, V. (1973). Collective properties of physical systems. 1st ed. Stockholm: Nobel Foundation, p.130. [12] J. Murgich and S. Pissanetzky, 14N nuclear quadrupolar resonance (NQR) Chem. Phys. Lett. 18, 420 (1972) . [13] Foresman, J. and Frisch, A. (1996). Exploring chemistry with electronic structure methods. 1st ed. Pittsburgh, PA: Gaussian, Inc., pp.3-11.
33
[14] Font, M. (n.d.). Química computacional aplicada-mecánica molecular Departamento de química orgánica sección de modelado molecular, Universidad de Navarra [15] Font, M. (n.d.). Química computacional aplicada-mecánica cuántica. Departamento de química orgánica sección de modelado molecular, Universidad de Navarra [16] Foresman, J. and Frisch, A. (1996). Exploring chemistry with electronic structure methods. 1st ed. Pittsburgh, PA: Gaussian, Inc., pp 111-139.
[17] Foresman, J. and Frisch, A. (1996). Exploring chemistry with electronic structure methods. 1st ed. Pittsburgh, PA: Gaussian, Inc., apéndices pp.272 -278.
[18] Adamo, C. & Barone, V. Toward reliable density functional methods without adjustable parameters: the PBE0 model. J. Chem. Phys. 110, 6158–6170 (1999).
[19] Méndez, H., Heimel, G.,(2013). Doping of Organic Semiconductors: Impact of Dopant Strength and Electronic Coupling. Angewandte Chemie, 125(30), pp.7905-7909.
[20] Zhu, L. Y., Kim, E. G., Yi, Y. P. & Bredas, J. L. Charge transfer in molecular complexes with
2,3,5,6-tetrafluoro-7,7,8,8-tetracyanoquinodimethane (F-4TCNQ): a density functional
theory study. Chem. Mater. 23, 5149–5159 (2011).
Anexos
34
Apéndice 1. Niveles de energía calculado para donor/aceptor
35
Apéndice 2. Espectros Raman, IR y UV-VIS de oligotiofenos y FxTCNQ
1T
2T
36
3T
37
4T
5T
38
6T
7T
39
8T
2a
40
2b
2c
41
2d
42
Apéndice 3. . Participación en Tercera escuela de Química
Computacional de las Américas (PREFALC)
43
Apéndice 4. Participación en VI encuentro Nacional de Químicos
teóricos y computacionales
44
Apéndice 5. III Escuela Colombiana de teoría y computación en las