Tese Doutorado Fernanda Bettanin
Tese Doutorado Fernanda Bettanin
Tese Doutorado Fernanda Bettanin
QUÍMICA
UNIVERSIDADE FEDERAL DO ABC
FERNANDA BETTANIN
TESE DE DOUTORADO
2015
UNIVERSIDADE FEDERAL DO ABC
CENTRO DE CIÊNCIAS NATURAIS E HUMANAS
Programa de Pós-Graduação em Ciência & Tecnologia - Quı́mica
FERNANDA BETTANIN
Santo André
05/2015
AGRADECIMENTO
Aos meus pais por terem me apoiado nessa fase tão desafiadora da minha vida, me
dando suporte psicológico e financeiro.
À minha orientadora, Paula Homem de Mello, por todo o conhecimento passado, sendo
sempre atenciosa, motivadora e compreensiva.
Aos meus colegas do Laboratório ABCSim por todo carinho e por me aturar até mesmo
nos piores dias.
Aos professores colaboradores Maurı́cio D. Coutinho Neto, Ronei Miotto, Rodrigo
Cordeiro, Anderson Orzari Ribeiro e Luı́s Gustavo Dias.
À Profa. Dra. Benedetta Mennucci por ter me possibilitado o estágio em Pisa, na
Itália, e ao grupo Molecolab por todo o apoio.
Às amigas-irmãs que fiz na Itália, “As lindas”, sem vocês seria muito mais triste e
difı́cil.
À UFABC e Cenapad pela infra-estrutura e a FAPESP pelo apoio financeiro.
“Se fracassar, ao menos que
fracasse ousando grandes feitos,
de modo que a sua postura não
seja nunca a dessas almas frias e
tı́midas que não conhecem nem
a vitória nem a derrota.”
Theodore Roosevelt
Claude Lévi-Strauss
RESUMO
Figura 1 – Processos fı́sicos que podem ocorrer após uma molécula absorver um
fóton. S0 é o estado eletrônico fundamental da molécula. S1 e T1 são
os primeiros estados excitados singleto e tripleto, respectivamente. As
setas vermelhas representam transições radiativas e as setas verdes
são as transições não-radiativas. A é absorção, F , fluorescência, P ,
fosforescência, CI, conversão interna, CIS, cruzamento intersistemas e
R, relaxação vibracional (Diagrama de Jablonski adaptado de Ochsner,
1997) [47]. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
Figura 2 – Esquema do processo da PDT, em que E é energia. . . . . . . . . . . . 20
Figura 3 – Mecanismos I e II de ação em PDT. S0 é estado fundamental singleto,
T0 é estado fundamental tripleto, S1 é estado singleto excitado e T1 é
estado tripleto excitado. . . . . . . . . . . . . . . . . . . . . . . . . . . 21
Figura 4 – Figura esquemática de uma membrana de bicamada lipı́dica. . . . . . . 23
Figura 5 – Figura esquemática do fosfolipı́deo palmitoil-oleil-fosfatidil-glicerol (POPG). 23
Figura 6 – Estrutura de ftalocianina livre. . . . . . . . . . . . . . . . . . . . . . . 24
Figura 7 – Estrutura de porfirina livre. . . . . . . . . . . . . . . . . . . . . . . . . 24
Figura 8 – Estrutura de clorina livre. Observar ausência de uma dupla ligação no
anel esquerdo inferior, fato que difere uma clorina de uma porfirina. . . 25
Figura 9 – Espectro de absorção de uma ftalocianina livre (fornecido pelo Prof.
Anderson Orzari Ribeiro). . . . . . . . . . . . . . . . . . . . . . . . . . 25
Sumário . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
1 Introdução . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
2 Revisão da Literatura . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
2.1 Terapia Fotodinâmica . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
2.1.1 Mecanismos de Ação em PDT . . . . . . . . . . . . . . . . . . . . . 18
2.2 Membrana Fosfolipı́dica . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
2.3 Estrutura e Propriedades de Ftalocianinas, Porfirinas e Clorinas . . . . . . 24
2.4 Estudos Teóricos sobre os Sistemas em Questão . . . . . . . . . . . . . . . 27
3 Métodos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
3.1 Mecânica Quântica . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
3.2 Funções de Onda para Sistemas Moleculares . . . . . . . . . . . . . . . . . 32
3.3 Teoria do Funcional da Densidade . . . . . . . . . . . . . . . . . . . . . . . 34
3.4 Forças de dispersão . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
3.5 Modelo Contı́nuo de Solvatação . . . . . . . . . . . . . . . . . . . . . . . . 38
3.5.1 Membrana simulada com modelo contı́nuo polarizável . . . . . . . . 45
3.5.1.1 Construção do Perfil de Constante Dielétrica . . . . . . . . 49
3.6 Dinâmica Molecular . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50
3.6.1 Condições Iniciais . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51
3.6.2 Cálculo dos potenciais . . . . . . . . . . . . . . . . . . . . . . . . . 52
3.6.3 Equações de movimento . . . . . . . . . . . . . . . . . . . . . . . . 53
3.6.3.1 Algoritmo de Verlet . . . . . . . . . . . . . . . . . . . . . 54
3.6.3.2 Algoritmo de leap-frog . . . . . . . . . . . . . . . . . . . . 55
3.6.4 Cálculo de propriedades . . . . . . . . . . . . . . . . . . . . . . . . 56
3.6.5 Limitações e aspectos importantes da dinâmica molecular . . . . . . 56
3.6.6 Dinâmica Molecular de sistemas lipı́dicos . . . . . . . . . . . . . . . 57
3.7 Método hı́brido: Mecânica Quântica e Mecânica Molecular . . . . . . . . . 58
4 Objetivos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62
5 Metodologia . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63
6 Resultados e Discussão . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 75
6.1 Estudo Metodológico . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 75
6.1.1 Otimização de geometria . . . . . . . . . . . . . . . . . . . . . . . . 75
6.1.2 Espectro Vibracional . . . . . . . . . . . . . . . . . . . . . . . . . . 76
6.1.3 Espectro de absorção . . . . . . . . . . . . . . . . . . . . . . . . . . 78
6.2 Análise de propriedades . . . . . . . . . . . . . . . . . . . . . . . . . . . . 81
6.2.1 Propriedades de Solvatação . . . . . . . . . . . . . . . . . . . . . . 81
6.2.2 Propriedades eletrônicas . . . . . . . . . . . . . . . . . . . . . . . . 91
6.2.2.1 Espectro de absorção . . . . . . . . . . . . . . . . . . . . . 103
6.2.3 Análise de componentes principais (PCA) e hierárquica (HCA) . . . 110
6.3 Dı́meros . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 122
6.4 Interação entre fotossensilizador e membrana de bicamada lipı́dica . . . . . 131
6.4.1 Dinâmica Molecular . . . . . . . . . . . . . . . . . . . . . . . . . . 131
6.4.1.1 Equilı́brio da membrana . . . . . . . . . . . . . . . . . . . 131
6.4.1.2 Equilı́brio da molécula . . . . . . . . . . . . . . . . . . . . 134
6.4.2 Método hı́brido QM/MMPol . . . . . . . . . . . . . . . . . . . . . . 136
6.4.3 Membrana descrita por IEFPCM (Modelo contı́nuo eletrostático) . 137
6.5 Fosforescência . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 140
7 Conclusões . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 143
Referências . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 146
A Apêndice . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 165
16
1 Introdução
Quando uma molécula no seu estado fundamental (S0 ) absorve um fóton, ocorre
transição eletrônica, em que o elétron é promovido para um estado excitado mais energético.
Juntamente à transição eletrônica da molécula, podem ocorrer também transições rotacio-
nais e vibracionais. Assim, bandas de absorção eletrônica geralmente são muito largas
porque vários nı́veis rotacionais e vibracionais diferentes correspondem a nı́veis de energia
muito próximos [48, 49].
Quando uma molécula, que no seu estado eletrônico fundamental é singleto (possui
dois elétrons com spins opostos em um orbital molécular de mais baixa energia) absorve
um fóton, ela é promovida deste estado fundamental S0 para um estado eletrônico singleto
excitado S1 , onde um dos elétrons do orbital de mais baixa energia é impulsionado para
um orbital de energia mais alta, mantendo o mesmo spin. Geralmente, o primeiro processo
que ocorre após a absorção é a relaxação vibracional ao nı́vel vibracional mais baixo de S1 .
Essa transição sem radiação é chamada de R1 , em que a energia vibracional é transferida
para o meio em forma de calor [48, 49].
Outro processo que pode ocorrer sem radiação, é chamado de conversão interna (CI),
quando a molécula entra em um nı́vel vibracional S0 altamente excitado que possui a mesma
energia de S1 . Neste caso, quando a molécula volta ao estado fundamental transferindo
sua energia para o meio através de colisões, sem radiação, o processo é chamado de R2
[48, 49].
A molécula também pode passar de S1 para o nı́vel excitado tripleto T1 , processo
chamado de cruzamento intersistemas (CIS), pois é uma transição onde o elétron que foi
19
Célula(Cancerosa(
LUZ(
(2º)(
Estado(Excitado( Estado(Excitado(
Corante( (4º)(
(3º)( (5º)( Oxigênio(
fotossensível( E(
(1º)(
Estado( Estado(
Fundamental( Fundamental(
elevado para um nı́vel de mais alta energia, inverte seu spin, tornando a molécula tripleto
(estado eletrônico com um elétron desemparelhado). Depois disso, a molécula pode sofrer
a relaxação vibracional R3 , indo para um nı́vel vibracional de T1 de menor energia. A
molécula ainda pode sofrer outro cruzamento interssistemas para voltar para S0 , e por fim
passar pela relaxação vibracional R4 . A molécula sofre transição radiativa ao passar de S1
para S0 (entre estados de mesmo spin, chamada de fluorescência) bem como de T1 para S0
(chamada de fosforescência, transição radiativa entre estados de diferentes spins) [48, 49].
As taxas relativas a CI, CIS, fluorescência e fosforescência dependem da molécula,
do solvente e de condições de temperatura e pressão. O tempo de vida da fluorescência é
considerado curto (10≠8 a 10≠4 s), e o tempo de vida da fosforescência é relativamente
longo (10≠4 a 102 s). Em virtude das relaxações, a fluorescência e a fosforescência têm
menor energia do que a radiação absorvida, emitindo em comprimento de onda maior do
que absorvido [26, 48, 49].
Neste trabalho, os FS estudados absorvem luz sofrendo uma excitação de S0 indo para
o estado singleto excitado S1 . Depois disso, sofrem relaxação não-radiativa realizando
cruzamento intersistemas populando o estado excitado T1 , cuja duração é maior que o
estado S1 . Quando sofrem decaimento para o estado fundamental S0 , podem transferir
energia para a molécula de oxigênio ou moléculas intermediárias ou até mesmo o solvente,
que se encontram nos tecidos e organelas celulares [47]. Essas moléculas, ao receber essa
energia vinda do FS, formam as chamadas espécies reativas de oxigênio: ânion radical
superóxido (O2•≠ ), peróxido de hidrogênio (H2 O2 ), dioxigênio singleto (1 O2 ), e radical
hidroxilo (HO• ). Alguns autores também consideram o peroxinitrito (ON OO≠ ) como
espécie reativa de oxigênio também.
A Figura 2 esquematiza o processo da PDT resumidamente incluindo os estados do
fotossensibilizador e uma das moléculas alvo, o oxigênio.
Existem dois mecanismos em PDT. O mecanismo l (fı́sico) consiste na transferência de
20
Figura 3: Mecanismos I e II de ação em PDT. S0 é estado fundamental singleto, T0 é estado
fundamental tripleto, S1 é estado singleto excitado e T1 é estado tripleto excitado.
21
da formação de agregados deve ser avaliada, pois a agregação pode alterar o comprimento
de onda de absorção do FS, ou seja, o mesmo deve ser solúvel no meio em que vai ser
administrado e no meio em que vai agir no organismo ou, mesmo agregado, apresentar
comprimento de onda de absorção compatı́vel com a terapia. O FS também deve ter
baixo custo de produção (sı́ntese fácil e com alto rendimento) para que haja interesse de
industrias farmacêuticas em produzi-lo.[27, 29, 30, 32, 49, 53].
A PDT é um tratamento aprovado em muitos paı́ses como Canadá, Alemanha, França,
Japão, Brasil, entre outros [54]. Porfirinóides, chamada de primeira geração de fotos-
R
sensibilizadores, tem alta produção de oxigênio singleto e o Photofrin• foi o primeiro
medicamento aprovado (1990) como fotossensibilizador para uso em PDT, possuindo o
maior rendimento quântico (0,89) encontrado até hoje. Outros medicamentos derivados
de porfirinas, segunda geração de fotossensibilizadores, foram aprovados pela FDA (Food
and Drug Administration) e orgãos regulatórios de outros paı́ses como por exemplo a
R
benzoporfirina DBP-MA, a protoporfirina PpIX, o Levulan• (ácido 1,5-aminolevulı́nico -
ALA), o aminolevulinato de metilo (MAL), o hexaminolevulinato (HAL) e o Visudyne
•R R
(Vertiporfin). Na Rússia, foi aprovado o medicamento Photosens• , uma ftalocianina.
R R
Existem também três medicamentos feitos de clorinas aprovados: Foscan• , LS11• e
R
Photochlor• . Além disso, alguns outros medicamentos não-porfirı́nicos foram aprova-
dos ou estão em fase de testes, dentre antraquinonas, fenotiazinas, xantenos, cianinas e
curcuminóides [53, 55].
Com intuito de comparar as moléculas a serem estudadas com os medicamentos
já aprovados, foram incluı́das as moléculas referentes aos medicamentos no presente
R
estudo: Levulan• , produzido pela DUSA Pharmaceuticals, Inc., tem grande eficiência
no tratamento de lesões causadas por acne, [56] e carcinoma de células basais e células
R
escamosas, cânceres urotelial e da bexiga urinária [57]; Visudyne • , produzido pela
Novartis Pharmaceuticals, é um medicamento aprovado para o tratamento de distúrbios
vasculares oculares, como miopia patológica e sı́ndrome de histoplasmose ocular, além de
neovascularização subforeal coróide devido degeneração macular relacionada à idade [58];
R
e Foscan• , produzido pela Biolitec Pharma Ltd., é considerado um medicamento de alta
eficiência para tratamento de cânceres de cabeça e pescoço [53, 59].
22
Figura 4: Figura esquemática de uma membrana de bicamada lipı́dica.
HO O O
P O
O
O O≠ Na+
23
Figura 6: Estrutura de ftalocianina livre.
24
Figura 8: Estrutura de clorina livre. Observar ausência de uma dupla ligação no anel esquerdo
inferior, fato que difere uma clorina de uma porfirina.
Fernanda Bettanin fez ótima apresentação e está aprovada.
Qx
Soret Qy
Figura 9: Espectro de absorção de uma ftalocianina livre (fornecido pelo Prof. Anderson Orzari
6
Ribeiro).
25
para o vermelho da principal banda de absorção causado pela saturação da dupla ligação
[84].
Para aumentar a eficiência das porfirinas livres, estudos propõem modificações em
seu centro metálico e em seus ligantes quanto ao número, posição e natureza de forma
a alterar suas propriedades espectrofotométricas, fotoquı́micas e de agregação [2]. Por
exemplo, é possı́vel aumentar a intensidade de absorção na região do visı́vel nas porfirinas
livres mudando sua estrutura básica com substituintes, mudanças na via de conjugação
ou adicionando um substituinte central [85]. O acréscimo de ligantes axiais para uti-
lização em sensores quı́micos tem grande influência em processos redox e propriedades
fotovoltaicas [7–9] e, portanto, em funções biológicas e catalı́ticas. Estes ligantes po-
dem trazer propriedades desejáveis como grande momento de dipolo, polarizabilidade
ou hiperpolarizabilidade [7]. Porfirinas com manganês têm aplicações como agentes de
aprimoramento de imagem de RMN, agentes de ligação e de radiodiagnóstico, clivagem de
DNA e oxidantes alimentı́cios [9]. A otimização de coeficientes de extinção das bandas Q
e B são um benefı́cio para as propriedades de captação de luz [21]. Os estados eletrônicos
mais elevados das porfirinas, acessados através do bombeamento da banda B e o efeito da
agregação nas suas propriedades ópticas têm despertado grande interesse [2, 3, 7–10, 86].
Ftalocianinas e clorinas possuem grande capacidade redox e possibilidade de se coorde-
nar a diferentes metais e, portanto, podem ser utilizadas em muitas aplicações citadas para
porfirinas [27, 40, 73, 74]. As metaloftalocianinas levam vantagem em relação às porfirinas
pela consideravelmente alta estabilidade térmica e por isso têm sido propostos como co-
rantes industriais, como catalisadores, em eletrodos de difusão de gás, em polı́meros eletro
e fotocondutores, em dispositivos semi-condutores, em eletrofotografia, em sensores de gás
semicondutores, em filmes de Langmuir-Blodgett, em cristal lı́quido, em catalisadores para
cátodos de célula combustı́vel e biossensores [3, 73]. Por sua atividade antimicrobiana,
clorinas também vêm sendo estudadas para terapia fotodinâmica sendo mais susceptı́veis
a mudanças na banda Q e red-shift quando com certos substituintes que as porfirinas
[27, 40]. Além disso, essas três classes de compostos são muito estudadas por participarem
de sistemas naturais, o que permite seu uso biológico na geração de oxigênio para aplicação
em PDT, por não apresentarem toxicidade na ausência de luz [87].
Apesar da grande quantidade de estudos experimentais, ainda pouco se sabe sobre os
efeitos dos substituı́ntes e do solvente nos estados excitados e sobre os mecanismos de
transferência eletrônica e de energia mediada por essas espécies [2, 3, 7, 8, 86].
26
crescido o número de estudos ab initio, principalmente com DFT (Teoria do Funcional
da Densidade) que apresenta precisão considerável e relativo baixo custo computacional
[2, 7, 8, 10, 21, 22, 24, 25, 27, 73, 86, 89–92]. Além das porfirinas, a DFT também tem sido
utilizada para o estudo das estruturas eletrônicas e obtenção de geometrias de ftalocianinas
e clorinas [73, 74].
Dentre as principais propriedades analisadas pelos trabalhos pesquisados estão: es-
trutura, espectro de absorção [2, 10, 21, 22, 24, 25, 27, 86, 89–92], propriedades eletrônicas
e emissão [2, 7–9, 21, 24, 25, 27, 73, 86, 90, 91], influência dos ligantes e do metal central
[2,7–9,21,22,25,27,33,73,85,89–91], potencial de ionização e afinidade eletrônica [8,9,24,25],
dipolo de transição e cargas atômicas [73, 92].
A maioria dos trabalhos sobre porfirinas e ftalocianinas com DFT utilizaram B3LYP
como funcional [2, 7, 10, 21, 22, 24, 25, 73, 74, 86, 89–92]. Outros funcionais utilizados foram
BLYP [27, 73], VWN-B-P [8], além de BHandHLYP [27]. T. E. Shubina e T. Clark [7]
realizaram um estudo metodológico para porfirinas de F e e Co, substituı́das com N O
e CO, e compararam os resultados obtidos com os funcionais B3LYP, OLYP, B3PW91,
M06L, com o método MP2 (as bases utilizadas foram 6-31G (d,p) e 6-311G (d,p)). Os
autores concluı́ram que B3LYP e OLYP/6-31G(d,p) apresentaram os melhores resultados
até mesmo do que o método MP2.
Zhang et al. [12] realizaram cálculos ab initio de espectros IV e Raman [93] e estimaram
o efeito da presença de isótopos no espectro IV para ftalocianinas e seus análogos.
As bases mais utilizadas são as de Pople (6-31G ou 6-311G com diferentes combinações
de funções de polarização e difusas) [2, 7, 21, 22, 24, 25, 73, 74, 86, 89–92]. No caso especı́fico
de clorinas, os funcionais BLYP e B3LYP, com a base 6-31G(d), forneceram resultados
satisfatórios em tempo computacional moderado [27, 74]. Esta famı́lia de bases tem sido
utilizada vastamente na literatura para sistemas similares aos tratados neste trabalho
[2, 7, 21, 23–25, 27, 73, 74, 86, 89–92, 94].
Para obtenção dos espectros e absorção, diversas metodologias já foram empregadas,
como TD-DFT (Teoria do Funcional da Densidade Dependente do Tempo), o método
mais comum para essa aplicação, [21, 24, 86, 89, 91, 92], TD-HF (Hartree-Fock Dependente
do Tempo) [8, 24, 25, 86, 89], ZINDO/S (método semi-empirico parametrizado para obter
espectros de absorção) [24,86], CI (Interação de Configurações) [24,25] e CASSCF/CASPT2
(Complete Active Space – Self-Consistent Field/ Second Order Pertubation Order) [7].
Porém, esses últimos métodos demandam um alto custo computacional.
A geometria planar de porfirinas, ftalocianinas e clorinas, principalmente quando
complexadas, lhes causam empilhamento devido às interações de dispersão geradas pelos
elétrons fi conjugados em seus sistemas aromáticos [95]. Esse empilhamento, ou seja,
a formação de dı́meros ou trı́meros, modificam o espectro de absorção em relação ao
monômero, e essa modificação pode levar à diminuição da eficiência da PDT [53, 87, 96].
Portanto, estudos vem sendo realizados com o intuito de entender essas estruturas e de
27
diminuir a agregação e seus efeitos [97,98]. No caso dos estudos teóricos, a DFT vem sendo
usada sempre com uma correção de longo alcance, já que o entendimento das forças de
dispersão (ligações fracas como van der Waals) dependem disso, e não são bem descritas
pelos funcionais desse método [79, 99, 100]. Uma das correções utilizadas são as DCACP’s
(do inglês, Dispersion Corrected Atom Centered Pseudopotentials) e existem trabalhos com
dı́meros utilizando-as [101–103] apresentando bons resultados.
Analisando os trabalhos presentes na literatura, é clara a importância de adicionar os
efeitos do meio na simulação [104–110]. Hoje os modelos de solvatação mais utilizados
em cálculos de mecânica quântica são os contı́nuos, como PCM (Polarizable Continuum
Model), que descrevem o solvente como um contı́nuo dielétrico. Apesar de ser uma
aproximação limitada, pois o solvente é modelado em termos de propagação de carga
superficial aparente numa cavidade de forma molecular, sua facilidade de utilização e
resposta razoável o tornam o modelo escolhido para este trabalho [104, 107].
Outra propriedade importante para desenvolvimento de novos fármacos para PDT,
dificilmente obtida experimentalmente e pouco relacionada com a eficiência fotodinâmica
nos trabalhos atuais, é a afinidade do corante com a membrana celular, onde a concentração
de oxigênio é maior que no meio aquoso [72, 111–113]. Alguns estudos indicam que,
preferencialmente, a molécula deve ter uma estrutura anfifı́lica: parte hidrofı́lica para
garantir solubilidade em água, viabilizando a administração intravenosa; e parte hidrofóbica
que permite a interação com a membrana fosfolipı́dica [111, 114–119]. Para simulação
computacional dessa propriedade, é necessário considerar também a posição e orientação
do corante em relação à membrana [72].
É desafiante estudar computacionalmente sistemas grandes e complexos como membra-
nas, biomoléculas e interfaces, por causa do alto custo e dificuldade de descrever efeitos de
equilı́brio relacionados a sistemas fı́sico-quı́micos como processos dinâmicos e flutuações
termodinâmicas. Simulações computacionais têm sido capazes de gerar informações sobre
estrutura e interações a nı́vel molecular e diversos pesquisadores tem utilizado diferentes
técnicas para a interação de corantes fotossensı́veis com membranas fosfolipı́dicas [120–125].
A aproximação atomı́stica tem sido utilizada em vários estudos [30, 126, 127] porém se é
necessária a adição de moléculas de solvente, como é o caso de simulação de membranas, o
custo computacional aumenta significativamente [128–130].
Mas nem sempre é necessária a aproximação atomı́stica. Por exemplo, Iozzi et al.
tem empregado o modelo contı́nuo de solvatação PCM (Polarizable Continuum Model)
[106,109] para simular sistemas biológicos [131]. Além de auxiliar na simulação de espectros
de absorção, o PCM é o modelo contı́nuo de solvatação mais utilizado para aplicações
envolvendo cálculos acurados de mecânica-quantica, como estudos de termodinâmica,
cinética, estrutura e espectroscopia em meios tanto isotrópicos quanto anisotrópicos e
interfaces entre dois fluidos [106, 131–136].
Nesse sistema soluto-solvente feito por PCM, é possı́vel obter a energia de interação de
28
solvatação, e assim a distribuição da probabilidade de posição, orientação e conformação do
soluto, no seu estado fundamental e excitado (não acessı́vel experimentalmente), levando
em consideração a não uniformidade e anisotropia do meio nos dois nı́veis. E, unindo
cálculos quânticos de cargas do soluto, também é possı́vel analisar a estabilidade do soluto
através da bicamada lipı́dica, além de deslocamento das bandas de emissão e absorção,
assim como perfis de transferência de energia livre dos mesmos, ao longo do meio [137,138].
O uso de PCM para simulação de membranas se tornou possı́vel a partir de Ster e
Feller, que explicaram parte do problema eletrostático de uma membrana de bicamada
[139]. Tanizaki e Feig [140], além de Im e colaboradores [141], resolveram o mesmo
problema para membranas biológicas. Este modelo não promove detalhes sobre a estrutura
atômica dos lipı́dios da membrana, mas descreve os efeitos da bicamada e das interações
eletrostáticas [142]. Usando a aproximação PCM para descrever membranas, o cálculo
se torna relativamente mais rápido que utilizando a aproximação atomı́stica, e promove
propriedades eletrônicas que não são alcançadas com cálculos de dinâmica molecular nem
mesmo com o custoso método QM/MM.
Além disso, vários estudos com PCM mostraram boa performance ao estudar peptı́deos
hidrofóbicos interagindo com a região de hidrocarbonetos de membrana lipı́dica [106, 143–
148]. Kessel e colaboradores comparam uma simulação de dinâmica molecular com um
modelo contı́nuo de solvatação e obtiveram bons resultados interagindo ácido valpróico
(molécula anfifı́lica) com membrana de bicamada lipı́dica [147].
29
3 Métodos
Ĥ =E (3.1)
C D
2 2
≠h 2 Ze
Ò ≠ (x, y, z) = E (x, y, z) (3.2)
8fi 2 m r
Na Eq. 3.2, o primeiro termo entre colchetes descreve a energia cinética e o segundo
termo descreve a energia potencial de um elétron de massa m, a uma distância r de uma
carga nuclear Z. E é a energia do elétron, e é sua carga e h é constante de Planck. O
quadrado de é a probabilidade de densidade eletrônica, que representa a probabilidade
de encontrar o elétron numa dada região [149].
Para um sistema molecular, o operador Hamiltoniano contém termos de energia cinética
dos núcleos e elétrons e termos de interação eletrostática em que os ı́ndices N e e dizem
respeito a núcleos e elétrons respectivamente 3.3.
≠h2 ÿ
M
1 2 h2 ÿ n
Ĥ = ÒA ≠ 2 Ò2a
8fi A=1 MA
2 8fi m a=1
¸ ˚˙ ˝ ¸ ˚˙ ˝
TN Te
M ÿ
ÿ n ÿM ÿ M ÿn ÿ n
ZA Z A ZB Za Z b
≠e2 + e2 + e2 (3.3)
A=1 a=1 rAa A=1 B>A RAB a=1 b>a Rab
¸ ˚˙ ˝ ¸ ˚˙ ˝ ¸ ˚˙ ˝
VN e VN N Vee
C D ⁄
~2 Ò2 2 “(r, rÕ ) HF Õ
≠ + ‚(r) + ‚H (r) „i (r) ≠ q
HF
d3 rÕ „ (r ) = ‘HF
i „i (r)
HF
(3.4)
2m |r ≠ rÕ | i
Com a aproximação, o primeiro termo da Equação 3.4, energia cinética dos núcleos,
pode ser negligenciado e o quarto termo, a repulsão entre núcleos, é considerado uma
constante. Como qualquer constante somada a um operador, são apenas alterados os
auto-valores e não há efeito sobre as auto-funções. Deriva-se, então, um novo Hamiltoniano
eletrônico (Equação 3.5) [149].
N
ÿ 1 2 ÿ N ÿ M
ZA ÿ N ÿ N
1
Ĥ = ≠ Òi ≠ + (3.5)
i=1 2 i≠1 A=1 riA i=1 j>i rij
A Equação 3.6 mostra a equação de Schrödinger envolvendo o Hamiltoniano eletrônico
e a sua solução é dada pela solução da função de onda eletrônica (Equação 3.7), que
descreve o movimento dos elétrons e depende parametricamente das coordenadas nucleares,
assim como a energia eletrônica (Equação 3.8), ou seja, para cada arranjo nuclear existirá
uma nova função de onda eletrônica. A energia total para núcleos fixos deve incluir a
constante de repulsão nuclear (Equação 3.9) [149].
31
Figura 10: Diagrama de superfı́cie potencial para molécula diatômica, dependente da energia
total da molécula e a distância entre seus núcleos RA [149].
Da mesma forma que é calculado o Hamiltoniano eletrônico para uma dada distância
entre os núcleos e substituem-se as coordenadas eletrônicas pelos seus auto-valores médios,
é necessário calcular um Hamiltoniano nuclear (Equação 3.10) para descrever o movimento
dos núcleos, dentro da média dos elétrons. O Hamiltoniano nuclear é descrito pela Equação
3.11, em que a energia total fornece um potencial para o movimento nuclear, que constitui
uma superfı́cie potencial [149].
M
ÿ 1 2
ÿM ÿ M
Z A ZB
Ĥnuclear = ÒA + (3.10)
A=1 2MA A=1 B>A RAB
M
ÿ 1 ÿM ÿ M
ZA Z B
Hnuclear = Ò2A + +
A=1 2MA A=1 B>A RAB
N
ÿ 1 2 ÿ N ÿ M
ZA ÿ N ÿ N
1
Òi ≠ + (3.11)
i=1 2 i≠1 A=1 riA i=1 j>i rij
32
Inicialmente, esses orbitais eram chamados de orbitais de Slater, que correspondiam a um
conjunto de funções que decaiam exponencialmente em relação à distância do núcleo. Mais
tarde, Frank Boys propôs que estes orbitais poderiam ser aproximados como combinações
lineares de orbitais Gaussianos, resolvidos computacionalmente com mais facilidade. Os
conjuntos de base descritos desta forma são chamados de GTOS (gaussian-type orbitals).
O menor destes conjuntos é chamado de conjunto de função de base mı́nima e contém um
número mı́nimo de funções de base necessário para representar todos os elétrons do átomo
[151, 152].
Primeiramente são feitas combinações de funções radiais para cada um dos átomos. A
combinação das funções radiais gera as funções atômicas que descrevem o orbital atômico.
Posteriormente, na Teoria dos Orbitais Moleculares (TOM), a aproximação de um orbital
molecular pode ser obtida a partir da combinação linear de seus orbitais atômicos (LCAO
– Linear Combination of Atomic Orbitals) para um sistema de N funções monoeletrônicas
„1 (x, y, z), . . . , „N (x, y, z), de acordo com a Equação 3.12 em que cµi (µ = 1, 2, . . . , N ) são
coeficientes a serem determinados e µu são funções de base monoeletrônicas conhecidas
como conjunto de base [151, 152].
N
ÿ
Âi = cµi „µ (3.12)
µ=1
É possı́vel estender a base para melhor adequá-la ao sistema a ser estudado. Entre as
formas mais comuns de estender as funções de base estão: (l) dobrar a função mı́nima
para cada orbital (ex. double-zeta - base mı́nima duplicada, triple-zeta – base mı́nima
triplicada), (ll) considerar apenas os orbitais de valência (split-valence basis set), (lll)
adicionar funções de polarização (p, d, . . .) e (lV) adicionar funções difusas (+, ++, . . .)
[152].
Nas bases de Pople, funções de polarização são denominadas por um asterisco “*” (dois
asteriscos significam que funções d e p foram adicionadas). Estas funções têm elevado
número quântico angular no conjunto de bases. Para melhor descrever a camada de
valência no caso de moléculas no estado excitado ou com carga negativa, são adicionadas
funções difusas, que possuem expoentes menores que as outras funções de base. Quando
em sua representação (+) é adicionado mais um sinal de soma (++) significa que também
foi incluida uma função gaussiana difusa do tipo s para átomos de hidrogênio [152].
Existe um grupo de bases chamadas de DGauss [153, 154] que foram otimizadas para
resolver as equações de Kohn-Sham e apresentam maior precisão e eficiência para estudos
de geometrias e energias moleculares [155]. A menor delas é chamada de DGauss Double
Zeta Valence with Polarization (DGDZVP).
33
3.3 Teoria do Funcional da Densidade
Uma vez que a função de onda envolve as coordenadas de todas as partı́culas, a
resolução da equação de Schrödinger se torna computacionalmente muito custosa, mesmo
com as aproximações de Born-Oppenheimer e Hartree-Fock. Uma alternativa a essa
abordagem é a Teoria da Funcional da Densidade (DFT, da sigla em inglês).
Conceitos como potencial quı́mico, dureza e maciez, ácidos e bases duros e macios,
equalização da eletronegatividade, entre outros, foram definidos a partir da DFT [151, 156,
157].
Os dois teoremas de Hohenberg e Kohn, formulados em 1964, mostram como determinar
o estado fundamental de um sistema, com um potencial externo, utilizando a densidade
eletrônica tridimensional como variável básica. Tais teoremas são uma reformulação da
mecânica quântica, baseada na densidade eletrônica em vez de funções de onda. Essa
densidade eletrônica fl(r) mede a probabilidade de encontrar um elétron no ponto de
coordenada r [158, 159].
O primeiro teorema diz que o potencial externo é um funcional único de fl(r) além
de uma constante aditiva. Ou seja, a densidade eletrônica determina o Hamiltoniano do
sistema, já que determina o número de elétrons e o potencial externo. A energia de um
sistema eletrônico é determinada pela densidade eletrônica (Equação 3.13, onde o v em
subscrito mostra a dependência do potencial externo) [158, 159].
E = E‚ [fl] (3.13)
Então, fl̃(r) define seu próprio ‚˜(r), o Hamiltoniano H̃ e ˜ (r1 , r2 , . . . , rn ). Esta última
função ˜ (r1 , r2 , . . . , rn ), pode ser utilizada como uma tentativa para o sistema com potencial
externo ‹(r). A energia E0 , de acordo com o princı́pio variacional, é dada pela Equação
3.15 [158, 159].
⁄ ⁄
E0 = E‚ [fl] = F [fl] + fl(r)‚(r)dr Æ E‚ [fl̃] = F [fl̃] + fl̃(r)‚(r)dr (3.15)
34
initio Hartree-Fock e pós-Hartree-Fock, além de apresentar um custo computacional mais
baixo do que métodos correlacionados tradicionais como interação de configurações, MP2
e coupled cluster [151, 156, 157].
Para a densidade eletrônica do estado fundamental, que satisfaz o princı́pio estacionário
requisitado pelo principio variacional (Equação 3.15), o potencial quı́mico dos elétrons, µ,
mede a tendência do elétron sair de um sistema em equilı́brio [151, 156, 157].
; 5⁄ 6<
” E‚ [fl] ≠ µ fl(r)dr ≠ N =0 (3.16)
A energia total pode ser escrita separando o potencial externo em um funcional trivial
da densidade eletrônica (Equação 3.17) [156].
⁄
E0 = È T̂ + Vˆe Í + fl(r)‚(r)dr (3.17)
1 ⁄ ⁄ fl(r1 )‚(r2 ) ⁄
E‚ [fl] = G[fl] + dr1 dr2 + fl(r)‚(r)dr (3.18)
2 |r1 ≠ r2 |
em que, G[fl] = Ts [fl] + Exc [fl]
O funcional da energia cinética de um sistema de elétrons que não interagem, Ts [fl],
tem a mesma densidade eletrônica do sistema de elétrons que interagem. Exc , além de
incluir o termo de troca e correlação, inclui a parte residual da energia cinética T [fl] ≠ Ts [fl],
sendo que T [fl] é a energia cinética exata para o sistema de elétrons que interagem [156].
A partir de um sistema de referência de elétrons que não interagem com um Hamiltoni-
ano que tenha um potencial local efetivo, ‚ef (r), é possı́vel escrever um novo Hamiltoniano
(Equação 3.19) [151, 156, 157].
1
H KS = ≠ Ò2 + ‚ef (r) (3.19)
2
Utilizando a mesma aproximação utilizada no método Hartree-Fock, é possı́vel escrever
a função KS , aproximada por um produto anti-simetrizado de N funções de onda de um
elétron, Âi (ri ), representado pelo determinante de Slater (para um sistema de camada
fechada) na Equação 3.20 [151, 156, 157].
- -
- Â1 (r1 ) Â2KS (r1 ) · · · ÂN (r1 ) --
- KS KS
- -
- Â1 (r2 ) Â2KS (r2 ) · · · ÂN (r2 ) --
- KS
1 KS
KS
=Ô -
- .. .. .. .. -
- (3.20)
N! -
- . . . . -
-
- KS -
-Â1 (rN ) Â2KS (rN ) · · · ÂN
KS
(rN )-
35
A partir da equação de Schrödinger de um elétron são obtidos os orbitais de Kohn-Sham
(KS), ÂiKS (Equação 3.21) [151, 156, 157].
3 4
1
≠ Ò2 + ‚ef ÂiKS = Ái ÂiKS (3.21)
2
Estabelecendo o potencial efetivo, é possı́vel estabelecer uma conexão entre o sistema
real e o sistema hipotético, tornando a densidade eletrônica resultante igual à densidade
eletrônica real (Equação 3.22), e assim se tem a energia cinética precisa (Equação 3.23)
[151, 156, 157].
N
ÿ
fls (r) = 2|ÂiKS (r)|2 = fl0 (r) (3.22)
i
N
ÿ 1
Ts [fl] = ÈÂiKS ≠ Ò2i ÍÂiKS (3.23)
i 2
A partir da expressão da energia, é obtido o potencial efetivo (Equação 3.24) [151, 156,
157].
⁄
fl(r1 )
‚ef (r) = ‚(r) + dr1 + ‚xc (r) (3.24)
|r ≠ r1 |
”Exc [fl]
em que, ‚xc (r) =
”fl(r)
36
Funcional Origem correlação Origem Troca Hartree-Fock (%)
B3LYP LYP B88 20
B3PW91 PW91 B88 25
OLYP LYP OPTX 0
M06-L M06-L M06-L 0
de van der Waals [159, 163]. Visando a necessidade de incluir a troca (para descrever a
indistinguibilidade eletrônica), foram desenvolvidos os funcionais hı́bridos, que incluem
diferentes porcentagens de troca Hartree-Fock [61,164,165]. Na Tabela 1 estão apresentadas
as origens da troca e da correlação de alguns funcionais e a porcentagem de Hartree-Fock
para cada um.
Em 1965, o método KS (Kohn-Sham) tornou possı́vel encontrar o funcional da energia
cinética exata (Ts ) para um sistema de uma partı́cula, dada pela Equação 3.25 [151].
N = - - >
ÿ - 1 2-
Ts = „i -≠ Ò -- „i
- (3.25)
i=1 2
A letra “s” subscrita mostra que a energia cinética foi calculada a partir de um
determinante de Slater. É possı́vel também escrever esta equação dependente da densidade
exata (Equação 3.26), a partir dos orbitais naturais (N O), vindos da matriz densidade
exata. Os números de ocupação no orbital, ni , correspondem ao número de elétrons em
cada orbital, portanto ficam entre 0 e 1 [151, 156, 157].
= - - >
{
ÿ - 1 2- N O
T[flexact ] = ni „N O -≠ Ò - „ (3.26)
i=1
i -
2 - i
{
ÿ - -2 {
ÿ
em que, flexact = ni --„N O -- eN= ni
i=1 i=1
Como a matriz exata não é conhecida, a densidade aproximada é escrita em funções
de onda auxiliares de um elétron, funções essas que descrevem os chamados orbitais de
Kohn-Sham (Equação 3.27) [151, 156, 157].
N
ÿ
fl(r) = |„i (r)|2 (3.27)
i=1
A inclusão deste termo significou uma grande melhoria na descrição do gap eletrônico
em diversos sistemas em que as aproximações LDA e GGA deixavam a desejar, assim
como na descrição de comprimentos de ligação [166]. A grande desvantagem deste método
é que os elementos da matriz de troca são não locais, o que gera um custo computacional
muito grande.
37
3.4 Forças de dispersão
Interações fracas, como forças de van der Waals, não são bem descritas pela DFT.
Os funcionais GGA não reproduzem bem interações em torno do mı́nimo de van der
Waals e o funcional LDA superestima essas interações. De qualquer forma, os dois tipos
de funcionais não são capazes de reproduzir a curva de longo alcance das interações de
dispersão [103, 167]. Tendo em vista a necessidade de incluir interações de van der Waals
em muitos estudos, como por exemplo a formação de agregados, tem-se desenvolvido
aproximações para descrever melhor a energia de dispersão.
Parâmetros empı́ricos ab initio chamados DCACPs (Dispersion Corrected Atom-
Centered Potentials) foram desenvolvidos recentemente desejando corrigir as limitações
dos funcionais GGA (generalized gradient approximation) [103, 167].
A correção DCACP Õ s consiste na introdução de potenciais efetivos centrados no núcleo
de cada átomo. Cada potencial é determinado por um processo de calibração dado por
cálculos de referência em sistemas representativos. Já foram parametrizados os átomos de
H, N, C, O, S, Si, P, halogênios e gases nobres [167]. O método é baseado no funcional
representado na Equação 3.28 de forma minimizada, onde rmin é a distância de equilı́brio
molecular e rmid é a distância no ponto médio dos cálculos de referência. Eref (r) é a
energia de interação de referência em r, FIˆ é termo de força para átomo Iˆ e wIˆ é um
fator de ponderação para que FIˆ e as contribuições de energia tenham a mesma ordem de
magnitude [103].
5 6
q
min
{‡i }
P (r)= {‡
min
i}
|E ref (r min )≠E(rminÕ ≠{‡i })|2 +|E ref (r mid )≠E(rmid {‡i })|2 + w1 |F1 (rminÕ {‡i })|2 (3.28)
i
O potencial DCACP, definido na Equação 3.29, é definido pelas variações dos parâmetros
‡i , que minimizam o funcional da Equação 3.28 e definem os potenciais não-locais do
átomo centrado em RI , onde r = |r ≠ RI | é a distância da posição do núcleo I, o vetor r é
o vetor unitário na direção de r ≠ RI e Yl m denota um harmônico esférico [103].
+1
ÿ
V1DCACP (r, rÕ ) = Ylm (r̂)p1 (r)‡1 p1 (rÕ )Ylm
ú
(r̂Õ ) (3.29)
m=≠1
38
concordância com resultados experimentas. Muitos modelos a utilizam [105, 106, 168] e
vem-se desenvolvendo melhorias no algoritmo computacional e parametrização de outros
solventes [108, 169–176].
Sabe-se que as interações eletrostáticas são as mais importantes para descrever um
solvente, além das forças de dispersão e repulsão de troca. As eletrostáticas, num modelo
contı́nuo dielétrico, são representadas pelo mesmo. As outras forças, quando incluı́das,
vem de modelos empı́ricos ou abordagens aproximadas [104].
Numa representação por contı́nuo dielétrico, cada molécula de solvente explicita
é substituı́da por dipolos pontuais, com movimento de dipolo µ ˛ , que interage com a
distribuição de carga do soluto, gerando a contribuição eletrostática. Assim, numa região
de volume V , contendo n moléculas de solvente, tem-se o vetor de polarização P˛ [14],
conforme Equação 3.30, onde µ ˛ i é o momento de dipolo da molécula i e <> significa uma
média térmica.
1 ÿ n
P˛ = <µ
˛> (3.30)
V i=1
Ou seja, P˛ é o vetor dipolo médio por unidade de volume. Os dipolos discretos então
são substituı́dos por um campo vetorial contı́nuo representado por P˛ . Mas esse, não
descreve da melhor forma a interação soluto solvente, apenas descreve macroscopicamente
o dielétrico, mas o modelo contı́nuo dielétrico ignora essa questão.
Esse modelo também assume que P˛ segue a aproximação de resposta linear, ou seja,
na presença de um campo elétrico E, ˛ o vetor P˛ segue a relação P˛ = ‰E, ˛ sendo chi a
suceptibilidade elétrica do meio. Este último introduz o comportamento macroscópico do
dielétrico, uma forma mais simples de representar a Equação 3.30. O vetor E ˛ é gerado
pelas cargas do soluto e pelos dipolos do solvente, que são representados por P˛ .
Além disso, nesse modelo é necessário definir uma cavidade para o volume que o soluto
ocupa, e entorno dessa cavidade fica o contı́nuo dielétrico. O tamanho da cavidade é
definida por esferas centradas em cada átomo do soluto, se sobrepondo, formando assim
uma cavidade realista [104].
O potencial eletrostático Ïdip (R),
˛ gerado pelo contı́nuo dielétrico em um ponto R, ˛
para um dipolo µ ˛ localizado na posição ˛r, é representado pela Equação 3.31.
1 µ ˛ (R
˛ ≠ ˛r)
Ïdip (R)
˛ = (3.31)
4fiÁ0 |R˛ ≠ ˛r|3
Se a Equação 3.31 for integrada considerando todo o espaço que o contı́nuo dielétrico
ocupa, se obtém o potencial eletrostático Ïd gerado pelo dielétrico (Equação 3.32).
1 P˛ (R
˛ ≠ ˛r)
Ïd (R)
˛ = dVr (3.32)
4fiÁ0 |R
˛ ≠ ˛r|3
39
Essa forma do potencial eletrostático está envolvendo dipolos, porém seria mais útil se
envolvesse cargas ao invés de dipolos. Para fazer essa transformação, usa-se a relação da
Equação 3.33.
S T
P˛ V˛r .P˛ (R
˛ ≠ ˛r)
V˛r U V = + P. (3.33)
˛ ≠ ˛r|
|R ˛ ≠ ˛r|
|R ˛ ≠ ˛r|3
|R
Substituindo a Equação 3.33 na Equação 3.32, tem-se a Equação 3.34:
S T
⁄
˛ = 1 P˛ V 1 ⁄ Vr .P
Ïd (R) V˛r U -- - dVr ≠ - - dV (3.34)
4fiÁ0 -R˛ ≠ ˛r-- 4fiÁ0 d --R˛ ≠ ˛r-- r
Sendo o contı́nuo dielétrico limitado a duas superfı́cies: S1, que contém o soluto, e S2,
que delimita o contı́nuo e o meio exterior, pode-se aplicar o teorema da divergência na
primeira integral da Equação 3.34, obtendo a Equação 3.35.
1 j P˛ .˛n 1 ⁄ Vr .P
Ïd (R)
˛ = - - dAr + - - dV (3.35)
4fiÁ0 S1 +S2 -˛
-R ≠ ˛ r-- 4fiÁ0 d --R˛ ≠ ˛r-- r
Podemos assim definir ‡S = P˛ .˛n e ‡V = Vr .P , assim a Equação 3.35 pode ser reescrita
dessa forma:
1 j ‡ (˛r) 1 ⁄ ‡V (˛r)
Ïd (R)
˛ = - S - dAr + - - dV (3.36)
4fiÁ0 S1 +S2 -˛
-R ≠ ˛
r-
- 4fiÁ0 d --R˛ ≠ ˛r-- r
j ⁄ j
(EÁ
˛ 0 ).˛ndAr = fl(˛r)dVr + P˛ .˛ndAr (3.39)
j ⁄
(EÁ
˛ 0 + P˛ ).˛ndAr = fl(˛r)dVr (3.40)
40
O vetor da primeira integral é o deslocamento elétrico (D),
˛ e sabendo que a permissivi-
dade do meio é Á = Á0 + ‰, então D˛ = ÁE.˛ A Equação 3.40 pode ser reescrita conforme
Equação 3.41.
j ⁄
˛ ndAr =
D.˛ fl(˛r)dVr (3.41)
˛ D
Ò. ˛ = fl(˛r) (3.43)
Sabendo que E
˛ é igual ao negativo do gradiente do potencial eletrostático, E
˛ = ÒÏ,
chegando na Equação 3.44.
fl(˛r)
Ò2 Ï(˛r) = ≠ (3.44)
Á
Onde, Ï é o potencial eletrostático total, obtido em vários pontos, ou seja, a soma
do potencial eletrostático do soluto (ÏS ) e do dielétrico (Ïd ), e fl é a densidade de carga
externa. ÏS é obtido da distribuição de carga do soluto,depois, sabendo Ï total, é possı́vel
calcular ÏSd . A Equação de Poisson serve tanto dentro do dielétrico quanto para o vácuo.
1 j ‡ (˛r) 1 j ‡ (˛r) 1 ⁄ Ò. ˛ P˛
Ïd (R)
˛ = - S - dAr + - S - dAr ≠ - - dV (3.45)
4fiÁ0 S1 -˛
r--
-R ≠ ˛ 4fiÁ0 S2 -˛
r--
-R ≠ ˛ 4fiÁ0 d --R˛ ≠ ˛r-- r
‡S (˛r) A
lim -- - Æ lim - - (3.46)
xæŒ ˛ ≠ ˛r-
-R
- xæŒ -
-R˛ ≠ ˛r--3
˛ = 1 j ‡ (˛r)
Ï(R) - S - dAr (3.47)
4fiÁ0 S1 -˛
-R ≠ ˛
r-
-
41
Outra expressão importante é a que dá as cargas de polarização (Equação 3.48).
dÏ
‡S (˛r) = (K ≠ 1)Á0 (3.48)
d˛n
onde, K = Á
Á0
=1+ ‰
Á0
[A]sol
= e≠ Gúsolv /RT
(3.49)
[A]g
Sabendo que o potencial quı́mico µ(A) do soluto, em fase lı́quida e gasosa, respectiva-
mente, é representada pela Equação 3.5.
No equilı́brio tem-se µsol (A) ≠ µg (A) = 0, e µúsol (A) ≠ µúsol (A) = Gúsolv (A). Assim
tem-se então a equação do potencial quı́mico de A em função do potencial quı́mico do
padrão de A como um gás ideal a 1 mol.L≠1 , do termo de concentração do soluto e da
energia livre de solvatação na Equação 3.51.
42
O processo de transferência pode ocorrer, reversı́vel, do soluto em uma posição fixa
em fase gasosa para uma posição fixa em solução. Ou seja, a energia livre de Gibbs é a
diferença entre os potenciais quı́micos em fase gasosa e lı́quida, que também é igual ao
trabalho (W) realizado sobre o sistema (Equação 3.52).
Carregando todos os pontos do espaço ocupado pelo soluto, faz-se a integral de todos
os elementos de carga dq:
⁄
dW = Ï(˛r, ⁄)d⁄fl(˛r)dVr (3.56)
1⁄
W = Ï(˛r)fl(˛r)dVr (3.58)
2
Quando usada no vácuo, a Equação 3.58, tem o potencial eletrostático gerado pela
molécula (Ï = ÏS ). Na fase lı́quida, esse potencial tem contribuição do dielétrico (Ï =
ÏS + Ïd ).Agora, sabendo o trabalho, pode-se escrever a expressão que dá a energia livre
de solvatação (Equação 3.59).
43
Utilizando a Equação 3.58, obtém-se:
1⁄ 1⁄
=
Gúsolv (ÏS (˛r) + Ïd (˛r))fl(˛r)dVr ≠ ÏS (˛r)fl(˛r)dVr (3.60)
2 2
que rearranjada dá a Equação 3.61:
1⁄
=
Gúsolv Ïd (˛r)fl(˛r)dVr (3.61)
2
Esta última equação mostra que a energia livre de solvatação é a metade da energia de
interação entre o dielétrico (representado pelas cargas de polarização de superfı́cie) e a
distribuição de carga da molécula [104].
O modelo de Born promove informações sobre como a energia de solvatação varia em
relação a carga, o raio da cavidade e a constante dielétrica do meio, conforme Equação
3.62, onde r é o raio da cavidade e q é a carga.
3 4
1 1 q2 1
Gúsolv = ≠ 1≠ (3.62)
2 4fiÁ0 r k
Observando a equação acima, tiram-se relações vistas experimentalmente como: dupli-
cando o valor da carga, a energia de solvatação quadruplica. Por outro lado, o aumento
do raio da cavidade, diminui a solvatação. Porém, uma falha desse cálculo se dá na
pequena variação da energia livre de solvatação diante a uma constante dielétrica alta. Isso
acontece pois, para descrever tal efeito, é necessário levar em conta interações especı́ficas
soluto-solvente que não são descritas neste modelo. Seria necessária uma cavidade que
dependesse do solvente ou até de acrescentar algumas moléculas explı́citas do solvente
[104].
Existe mais um termo necessário para o cálculo da energia livre de solvatação, o termo
que descreve a polarização do soluto, ou seja, que descreve como ocorre a deformação da
nuvem eletrônica do soluto diante das cargas de polarização do dielétrico. Sendo (0) a
função de onda do soluto em fase gasosa e (S) em solução, a Ep a energia de polarização
molecular é:
e - - f e - - f
(S) - 0- (S) (0) - 0- (0)
Ep = -Ĥ - ≠ -Ĥ - (3.63)
Ĥ S = Ĥ 0 + Vd (˛r, R)
˛ (3.64)
44
Então, incluindo a solvatação da molécula polarizada pelo dielétrico, dá a equação da
energia de solvatação (Equação 3.65),onde o primeiro termo é a interação do soluto com o
dielétrico e o segundo termo é a energia de polarização molecular.
1 e (S) -- - f
˛ -- (S) + Ep
Gúsolv = -Vd (˛
r, R) (3.65)
2
Dessa forma, substituindo a Equação 3.63 na Equação 3.65, tem-se a equação final da
energia de solvatação.
1 e (S) -- 0 - f e - -
˛ -- (S) ≠ (0) --Ĥ 0 -- (0)
f
Gúsolv = -Ĥ + (1/2)Vd (˛
r, R) (3.66)
2
O modelo contı́nuo de solvatação PCM [106, 109] é o mais conhecido e utilizado, além
de ser o modelo de solvatação mais acurado[106, 132, 133]. Neste modelo, a cavidade do
soluto é descrita por esferas sobrepostas, centradas em cada átomo do soluto, e o raio
depende da natureza de cada átomo. O potencial eletrostático gerado pelo dielétrico é
obtido pelo método de cargas aparentes de superfı́cie, e a densidade eletrônica, por meio
de cálculo ab initio, incluindo a polarização do soluto. A distribuição de cargas do soluto
é feita de forma contı́nua, sendo assim um método altamente realista [104].
Neste trabalho, é utilizado o modelo IEFPCM (Integral Equation Formulation Polari-
zable Continuum Model) [177, 178], que permite tratar em uma única abordagem, meios
complexos com dielétricos de natureza diferentes como lı́quidos isotrópicos, cristais lı́quidos
intrinsicamente anisotrópicos, matrizes sólidas ou soluções iônicas. Mais a fundo, este
método pode resolver o problema eletrostático de distribuição de carga da molécula dentro
da cavidade construı́da no meio complexo com uma equação linear envolvendo apenas
quantidades no limite da cavidade. Essas quantidades são definidas em termos da função
de Green [134].
Ë È
Òr. Á(˛r).ÒrG
˛ E (˛r, r˛Õ ) = 4fi”(˛r, r˛Õ ) (3.67)
45
A fim de estudar interfaces difusas planares, é necessário assumir Á(˛r) variando no eixo
z, perpendicular ao plano. Para explorar a simetria de translação ao longo dos eixos x e y,
se reescreve a função de Green pra o meio dielétrico:
⁄ Œ
G (˛r, r˛Õ ) =
E
GE (q1 , z, z Õ )J0 (q1 |(r˛|| ≠ r˛||Õ )q1 .dq1 (3.68)
0
Onde r˛|| e r˛||Õ são projeções de ˛r e r˛Õ no plano de interface e J0 é a função de Bessel de
ordem zero. Juntando as Equações 3.67 e 3.68, usando as propriedades das funções de
Bessel, obtém-se uma equação diferencial monodimensional para GE (q1 , z, z Õ ):
C D
d d
Á(z). GE (q1 , z, z Õ ) ≠ q12 Á(z)GE (q1 , z, z Õ ) = ≠2”(z ≠ z Õ ) (3.69)
dz dz
A solução desse tipo de equação pode ser escrita como a Equação 3.70,onde U1 e U2 são
duas soluções linearmente independentes da equação homogênea associada (Equação 3.71)
Y
2
_
] U (z)U2 (z Õ )
Á(z Õ )W (z Õ ) 1
æ z > zÕ
GE (q1 , z, z Õ ) = (3.70)
_ 2
[ U (z Õ )U2 (z)
Á(z Õ )W (z Õ ) 1
æ z < zÕ
C D
d d
Á(z). Un (z) ≠ q12 Á(z)Un (z) = 0 (3.71)
dz dz
Com z æ +Œ, U1 (z) æ 0 e z Õ æ ≠Œ, U2 (z) æ 0. W é a Wronskiana de U1 e U2 . Essa
função de Green obtida representa bem a singularidade de Coulomb. Deve-se considerar
esta singularidade da GE (z, z Õ ) já que ˛r e r˛Õ estão muito próximos, portanto os efeitos de
variação de Á estão escondidos pelo dominando 1/|˛r ≠ r˛Õ |. Portanto é necessário reescrever
GE (q1 , z, z Õ ) como na Equação 3.72, onde está incluı́da a contribuição da imagem da
transformada de Fourier da função Green, GE img (q1 , z, z ).
Õ
1
GE (q1 , z, z Õ ) = e≠q1 .|z≠z | + GE
img (q1 , z, z ) (3.72)
Õ Õ
C(z, z )q1
Õ
1
GE (˛r, r˛Õ ) = + GE
img (˛
r, r˛Õ ) (3.73)
C(z, z ).|˛r, r |
Õ ˛Õ
1
= lim q1 .eq1 |z≠z | .GE (q1 , z, z Õ ) (3.74)
Õ
C(z, z ) q1 æ+Œ
Õ
O grande ponto do IEF é que, o potencial eletrostático dentro de uma dada região I
(dentro da cavidade), pode ser descrito em termos de potencial e campo elétrico no limite
46
q
de uma região ( ) através de operadores integrais que envolvem a função de Green (G’)
adequada para aquela região. A expressão do potencial também contém uma integral de
volume na densidade de cargas localizada no interior dessa região (Equação 3.75)
⁄
4fiV (˛r) = (S .ˆ V )(˛r) ≠ (D .V )(˛r) + 4fi
I I I
GI (˛r, r˛Õ )fl(r˛Õ ).dr˛Õ (3.75)
q
Onde ˆ I u indica ÁI .n̂.Òu, e n̂ é o vetor normal a apontando para fora de I, e ÁI é a
constante dielétrica na região I. Os operadores S I e DI são definidos abaixo.
⁄
(S I .u)(˛r) = q GI (˛r, ˛s)u(˛s).d˛s (3.76)
⁄
(D .u)(˛r) = q ˆ˛sI GI (˛r, ˛s)u(˛s).d˛s
I
(3.77)
Esse potencial dentro da cavidade (região I), no vácuo, é uma simples função de Green.
Mas se há um meio (região E), é necessário levar em conta as contribuições de S I .‡ e
V, pois são a polarização molecular induzida do meio dielétrico, e será chamado de VR .
Para obter ‡ é necessário impor condições de contorno na superfı́cie da cavidade, ou seja,
a continuidade do potencial e do componente do vetor normal do vetor deslocamento
dielétrico (Equação 3.79).
Ë È
S E (2fi + DI ) + (2fi ≠ DE )S I .‡ = ≠S E .ˆV 0 ≠ (2fi ≠ DE )V 0 (3.79)
ú
Ë È
S E (2fiI + DI ) + (2fiI ≠ DE )S I A≠1 q = S E .E 0 ≠ (2fiI ≠ DE )V 0 (3.80)
ú
47
cada porção. As expressões de S I e DI estão em Mennucci, 1998 [179] e de S E e DE estão
em Frediani, 2004 [134].
Retomando a equação de energia eletrostática Eelec entre a densidade de carga da
molécula fl e o meio dielétrico, escrita em termos de fl e do potencial do campo de reação:
⁄
Eelec = d˛rVR (˛r)fl(˛r) (3.81)
Pode-se definir um Hamiltoniano Hef f para a molécula, adicionando uma energia de
interação eletrostática, apropriada para a densidade de carga molecular flM (˛r), para o
Hamiltoniano no vácuo. Em solução, a energia mais baixa do autovalor de Hef f , minimiza
um funcional energético com a forma da Equação 3.82.
1
G[flM ] = E0 [flM ] + Eelec [flM ] (3.82)
2
Onde E0 [flM ] é a energia da molécula em vácuo, calculado para um flM genérico. G[flM ]
é a energia livre. A contribuição eletrostática do solvente Gelec é dada pela Equação 3.83,
onde flsol
M é a densidade de carga molecular em solução e flM , em vácuo.
vac
Gelec = G[flsol
M ] ≠ E0 [flM ]
vac
(3.83)
Importante observar que não importa saber ao certo o que o meio dielétrico é de fato.
Uma vez que o modelo contı́nuo dielétrico é assumido, a única quantidade que resta saber
sobre o meio é o campo de reação.
Assumindo a aproximação de Born-Oppenheimer e a formulação de Kohn-Sham para
E0 [flM ] e, expandindo os orbitais à uma função de base finita, a minimização do funcional
da Equação 3.82 dá uma equação tipo Kohn- Sham que inclui termos relativos a interação
eletrostática molécula-dielétrico F Õ C = SCÁ , onde F Õ = F (P ) + 1/2(j + y + X(P )). Sendo
P a matriz densidade, F(P) é o operador Kohn-Sham no vácuo. A interação dos elétrons
da molécula com a polarização induzida do núcleo (y) é representada por j e a interação
eletrônica, induzida por elétrons, é representada por X(P). Essas matrizes são escritas em
termos de carga aparente de superfı́cie q. As expressões para j, y e X(P) são as mesmas
pra moléculas em meios homogêneos e podem ser achadas em Cammi, 1995 [180].
Como visto até agora, a diferença entre o cálculo para o estado fundamental entre o
meio homogêneo e um dielétrico não homogêneo, é nos valores da carga aparente q. Isso
também vale para o cálculo de energias de excitação, no mesmo esquema de TD-DFT,
resolvendo a equação do pseudo-autovalor:
SQ R Q RT Q R
A Bb 1 0 bV aXw b
Ua ≠wa =0 (3.84)
B A 0 ≠1 Yw
onde as matrizes A e B são dadas por:
0
Aai,bj = ”ab ”ij (Áa ≠ Ái ) + Kai,bj + Kai,bj
P CM
(3.85)
48
0
Bai,bj = Kai,jb + Kai,jb
P CM
(3.86)
ÿ
P CM
Kst,uv = qd (˛
si , [Âsú Ât ])V (˛
si , [Âuú Âv ] (3.87)
i
49
foram publicados [185, 186] e com esses perfis, foi possı́vel montar uma função que os
descrevessem (Equação 3.88).
5 3 46
dbulk z ≠ z0
d(z) = 1 ± tanh (3.88)
2 D
em que dbulk é a densidade aparente e o sinal ± depende da região do espaço que
o lı́quido considerado ocupa. A princı́pio, a densidade interfacial pode ser gerada pelo
cruzamento gradual entre os dois lı́quidos dependente do tempo, mas em muitos casos, a
interface permanece quase nı́tida. Já o modelo PCM leva em conta a média das várias
configurações de solventes na interação soluto-solvente, uma vez que o modelo descreve o
solvente através de uma média da constante dielétrica macroscópica [134].
Levando em consideração que para d=0, a constante dielétrica Á deve ser igual a 1,
a melhor forma de representar a dependência de Á em d é Á = 1 + kd, onde k é uma
constante que pode ser determinada assumindo Ábulk = 1 + kdbulk . Quando a interface está
entre dois lı́quidos diferentes, deve adicionar a suscetibilidade elétrica ‰ dos dois solventes
[‰ = (Á ≠ 1)/4fi].
5 3 46 5 3 46
Á1 ≠ 1 z ≠ z01 Á2 ≠ 1 z ≠ z02
Á(z) = 1 + 1 + tanh + 1 + tanh (3.89)
2 D1 2 D2
Assumindo z01 = z02 = z0 e D1 = D2 = D:
3 4
Á1 + Á2 Á1 ≠ Á2 z ≠ z0
Á(z) = + tanh (3.90)
2 2 D
Essa função é praticamente uma constante com o valor de Á1 para z >> z0 ou de Á2 para
z << z0 , dessa forma, ela descreve corretamente as constante dielétrica aparente dos dois
meios. Na região de interface, a permissividade muda suavemente de Á2 para Á1 . A largura
da zona interfacial pode ser ajustada através do parâmetro D. Pode ser verificado que, a
região em que a permissividade é sensivelmente diferente das suas constantes dielétricas
aparentes, é ¥ D/6. No caso, utiliza-se W=D/6 como uma estimativa da largura da
interface [134].
Para membranas, o perfil de permissividade é um pouco mais complexo (Equação 3.10).
Consideram-se três constantes dielétricas aparentes: Á1 para o lı́quido fora da membrana
(água), Á2 para a parte da monocamada e Á3 para o final desta (meio da bicamada). Duas
interfaces sucessivas definem a largura da monocamada, W12 , sendo que W12 = |z12 ≠ z23 |.
5 3 4 3 46
1 z ≠ z12 z ≠ z23
Á(z) = Á1 + Á3 (Á2 ≠ Á1 ) tanh ≠ (Á2 ≠ Á3 ) tanh (3.91)
2 D12 D23
Este último, foi o perfil utilizado para simular uma membrana de bicamada lipı́dica
durante o Programa Doutorado Sanduı́che realizado pela aluna junto ao grupo da Profa.
Mennucci na Università di Pisa, Itália.
50
3.6 Dinâmica Molecular
A Dinâmica Molecular (MD) foi introduzida por Alder e Wainright na década de 50,
mas apenas no inı́cio dos anos 70 cientistas computacionais começaram a obter resultados
similares aos experimentais. Rahman, em 1969, estudou propriedades do argônio lı́quido,
aplicando o potencial de Lennard-Jones e, a partir deste, muitos estudiosos procuraram
desenvolver um método para diferente tipos de sistemas a nı́vel atômico molecular [187].
Na década de 1980, uma série de campos de força para água, proteı́nas e DNA,
e programas de simulação por dinâmica molecular (AMBER [188], CHARMM [189],
GROMOS [190–192] e OPLS [193]) se tornaram disponı́veis, o que estimulou a aplicação
do método no campo da biofı́sica. Na sequência, estudos de sistemas lipı́dicos surgiram,
mesmo que no inı́cio com modelos simplificados de lipı́dios e ausência da influência de
solvente. Com o passar do tempo, os modelos foram estendidos para representar todos os
átomos do sistema e incluir efeito de solvente (água) e, a partir desses estudos, foi visto
que a dinâmica molecular, mesmo com algumas limitações, pode dar ideias detalhadas de
como se movimentam lipı́dios e proteı́nas [194].
Propriedades de transporte e equilı́brio, de um sistema clássico de muitos corpos,
são possı́veis de serem calculadas por simulação de dinâmica molecular. O sistema é
considerado clássico quando o movimento nuclear das partı́culas do sistema obedece a
mecânica clássica, o que é uma aproximação viável para muitos sistemas [187].
A MD é um método de simulação computacional que utiliza o movimento individual
de cada partı́cula de um sistema (átomos ou moléculas) para descreve-lo, podendo ele
ser sólido, lı́quido ou gasoso. Este movimento é descrito pela mudança de posições e
velocidades das partı́culas ao longo do tempo e esse conjunto de informações constroem
uma trajetória que é gerada através da integração das equações de movimento de Newton,
Lagrange ou Hamilton. Essa integração é feita por algoritmos que funcionam segundo o
esquema da Figura 11, em seguida, cada passo do algoritmo será comentada [187].
51
Figura 11: Fluxograma de um algoritmo genérico de Dinâmica Molecular [187].
ÿ ÿ ÿ ÿ ÿ
VT (r1 , r2 , ..., rn ) = VE + V◊ + VÏ + VLJ + VC (3.92)
Classicamente, pode-se representar uma ligação quı́mica como uma mola unindo dois
átomos. A lei de Hooke permite calcular o potencial harmônico linear, ou potencial
harmônico de estiramento (Equação 3.93), onde kE representa a constante de força de
estiramento, r é a distância entre os átomos e r0 é a distância no equilı́brio [187].
1
VE = kE (r ≠ r0 )2 (3.93)
2
O potencial harmônico angular representa a abertura ou fechamento do ângulo ◊ de
ligação formado por três átomos. Para obter o potencial pode-se adaptar a formula anterior
segundo Equação 3.94, onde K◊ é a constante de força e ◊0 é o ângulo entre os três átomos
no equilı́brio [187].
1
V◊ = k◊ (◊ ≠ ◊0 )2 (3.94)
2
O potencial de diedro descreve a torção de uma ligação, onde os ângulos de diedro são
„, ângulo entre os planos formados pelos átomos abc e bcd (Figura 12), sendo que o zero
corresponde a configuração CIS, em que a e d estão do mesmo lado [187].
Existem dois tipos de potencial de diedro: o impróprio e o próprio. O impróprio é
um potencial harmônico cuja função é manter a geometria planar de grupos como anéis
52
Figura 12: Representação do potencial de diedro [187].
aromáticos ou impedir que grupos girem. O potencial de diedro próprio permite a rotação
de grupos e pode ser descrito por uma série de Fourier truncada no terceiro termo, como
representado na Equação 3.95, onde n é o numero de mı́nimos para a torção de uma
determinada ligação quı́mica, „ é o ângulo diedral indicado na Figura 12 e é a defasagem
do ângulo diedral [187].
C 3
D
ÿ ÿ
V0 + Vn cod(nÏ + ) (3.95)
Diedro n+1
qa qb
VC = (3.97)
4fiÁ0 rab
53
Figura 13: Representação do potencial de Lennard-Jones [187].
Fi d2 ri (t)
Fi = mi ai , onde : ai = ∆ ai = (3.99)
mi dt2
Integrando a equação proposta, encontra-se as coordenadas da partı́cula passado o
tempo t. Este processo é repetido várias vezes até o tempo requerido, porém esta integração
não é possı́vel ser feita analiticamente, deve ser resolvida de forma numérica, por algoritmos
conhecidos. Aqui falaremos de dois deles [187].
54
Figura 14: Representação esquemática do algoritmo leap-frog, sendo r a posição e v as veloci-
dades no tempo [187].
Ë È
r(t + ”t) = 2r(t) + a(”t)2 ≠ r(t ≠ t”) + O (”t)4 (3.101)
O algoritmo leap-frog, uma variação do algoritmo de Verlet, foi utilizado neste trabalho.
Este calcula as posições das partı́culas no tempo t e suas velocidades na metade do intervalo
de integração, como mostrado na Figura 14 [187].
Fazendo a mesma expansão mostrada na seção anterior, chega-se na equação de
velocidade (Equação 3.102).
A B A B
”t ”t F (t)
v t+ =v t≠ + ”t (3.102)
2 2 m
A partir da velocidade pode ser calculada a posição da partı́cula (Equação 3.103).
A B
”t
v (t + ”t) = r(t) + v t ≠ ”t (3.103)
2
Assim como o algoritmo anterior, leap-frog também pode ser repetido até o tempo
requerido [195]. Este método foi escolhido por algumas vantagens que o mesmo retém:
rapidez na realização dos cálculos e aparição natural da velocidade no próprio algoritmo.
55
Entretanto, segundo a Figura 14, a velocidade encontrada sempre estará em defasagem com
a posição e consequentemente não será possı́vel ter a posição e a velocidade da partı́cula
ao mesmo tempo [187].
Por fim, a etapa de saı́da escreve (no output) as posições e velocidades no tempo de
integração requerido. A partir do output gerado, é possı́vel calcular a média de várias
propriedades do sistema, como energia, temperatura, pressão, entre outras [187].
N
ÿ mi vi2 (t)
T (t) = (3.105)
i=1 kB Nf
Ò
A flutuação da temperatura do sistema é da ordem de 1/ Nf , e Nf sendo da ordem
de 102 ≠ 103 , essa flutuação estatı́stica da temperatura é da ordem de 5-10%. Para uma
estimativa acurada da temperatura, são necessárias várias medições das flutuações da
mesma [194].
56
de campo de força, assim o efeito médio é retido porem não bem representado; é difı́cil
determinação da interação de van der Waals e potenciais de diedros; II) máximo de tempo
gasto em cada passo para tornar a dinâmica viável; e III) tratamento clássico do sistema
impossibilita considerar reações quı́micas sem tratar parte do sistema quanticamente,
porém isso não é uma limitação para simulação de sistemas lipı́dicos [196].
É importante que sejam levados em consideração o tempo e a escala em que os
fenômenos a serem estudados acontecem. Movimentos vibracionais estão dentre os mais
rápidos, poucos picosegundos. Isomerizações dos diedros da cauda lipı́dica ocorrem em
dezenas de picosegundos, e ficam mais devagar quando se aproximam da “cabeça” do
lipı́dio, demorando centenas de picosegundos. E os diedros dos grupos da “cabeça” dos
lipı́dios ainda demoram mais por conta da forte interação entre eles. Ou seja, a dinâmica
molecular de sistemas lipı́dicos são demorados, na escala de nanosegundos para rotações
de fosfolipı́dios, difusão lateral ou inserção de moléculas, por exemplo [196].
Outro ponto importante a ser levado em consideração é a adição de condições periódicas
de contorno. Em experimentos, é necessário utilizar um grande número de moléculas,
geralmente da ordem do número de Avogadro, o que é inviável para cálculos computacionais.
Porém, se reduzido esse número, efeitos de superfı́cie podem alterar os resultados. Portanto,
são utilizadas as condições periódicas de contorno, capazes de evitar interações com os
planos-limite. Isso acontece pois o método promove a interação das moléculas da superfı́cie
do sistema com as imagens da própria caixa, replicando o sistema infinitas vezes em todas
as direções do espaço. Dessa forma, se uma partı́cula “sair” da caixa, esta “entrará” do
lado oposto à sua saı́da, com a mesma velocidade que saiu [187].
Porém, o uso das condições periódicas de contorno, faz com que uma partı́cula possa
interagir com ela mesma através da sua imagem, ou então interagir mais de uma vez com
a mesma partı́cula. Além disso, para sistemas muito grandes, não é eficiente calcular
todas as interações. Portanto é necessário o uso de um raio de corte (cutoff ) R 6 L/2 que
impeça as interações duplas, onde L é o comprimento da aresta da caixa cúbica. Assim as
interações só ocorrem quando a distância entre as partı́culas é menor que R, interações
com outras partı́culas mais distantes são negligenciadas. Interações de van der Waals
rapidamente diminuem com a distância, já para interações coulombicas, principalmente
entre dipolos e cargas, essa distância é maior. O formato mais comum do cutoff é de dupla
esfera (uma pra todas as interações e uma maior para as eletrostáticas). Um método mais
acurado é o de Ewald que calcula as interações eletrostáticas resolvendo a equação de
Poisson [197–199].
57
Uma das primeiras considerações a serem feitas quando se deseja realizar uma dinâmica
molecular de sistemas lipı́dicos é sobre a escala de tempo: vibrações de ligações e ângulos,
movimentos libracionais e pequenas flutuações de diedro, assim como difusão e correlação
orientacional das moléculas de água e de pequenas moléculas, acontecem na mesma escala,
alguns picosegundos. Outros movimentos, como dos diedros da “cabeça” da membrana e
a interação delas, são mais lentos. A rotação da “cauda” no seu próprio eixo, é da escala
de nanosegundos e a difusão lateral de moléculas, dezenas de nanosegundos. Portanto, é
necessária uma avaliação sobre o dado a ser obtido com a simulação para saber o quanto
ela deverá durar. Quanto ao tamanho, hoje existe a limitação entre 100 e 200 moléculas
de lipı́deo, mas normalmente são utilizadas de 50 a 100 [196].
Normalmente a temperatura T e o número de partı́culas N são mantidos constantes
durante a simulação. Porém, o volume e a pressão devem ser ajustados portanto deve-se
definir as dimensões e pressão nas direções lateral e perpendicular. Mantendo a pressão
constante, é possı́vel obter a área por lipı́deo, um parâmetro (já mencionado anteriormente)
que pode também ser medido experimentalmente e auxilia a simulação para essas escolhas
[196].
A temperatura do sistema é dada pela energia cinética de todos os átomos, e ela tende
a se afastar do valo de partida devido o aquecimento gerado pelas forças externas ou ficção,
flutuações no processo de equilı́brio e eventuais erros de integração, portanto é necessário
um controle de temperatura. A pressão depende das forças e posições de todos os átomos
do sistema, determinando se o sistema expande ou contrai e assim o tamanho da caixa.
Os métodos de controle de temperatura e pressão são muitos, os dois mais utilizados são
weak coupling scheme [204], que é relativamente simples e não permite interpretações sobre
flutuações termodinâmicas no sistema, e o segundo é o extended system approach [205–207]
onde graus de liberdade são adicionados para simular um êmbolo para acoplamento de
pressão (barostato) e termostato para temperatura, sendo assim mais complexo, com
um conjunto mecânico-estatı́stico bem definido [152]. Neste trabalho foi utilizado como
barostato o algoritmo de Berensen [204], que altera as coordenadas e o tamanho da caixa
a cada passo. Quando a pressão muda em relação a uma pressão de referência, a pressão
do sistema relaxa de acordo com uma cinética de primeira ordem (Equação 3.106).
dP P0 ≠ P
≠ (3.106)
dT ·p
Para controle da temperatura, foi utilizado o “banho” termostático Langevin, que
corrige as alterações de temperatura do sistema lentamente, tendendo a temperatura
inicial.
58
3.7 Método hı́brido: Mecânica Quântica e Mecânica Molecular
Os métodos contı́nuos, como o PCM aqui descrito, são adequados para estudar propri-
edades eletrônicas de moléculas em meios homogêneos, e, até por esse motivo, também
foi descrito nesse trabalho uma forma de simular meios heterogêneos criando um gra-
diente dielétrico com PCM. Porém, para descrever propriedades eletrônicas em meios
heterogêneos, em que essas propriedades dependem da posição do soluto, o método hı́brido
chamado QM/MM (do inglês, quantum mechanical and molecular mechanical) é o método
mais conhecido para essa problemática e vem sendo bastante utilizado para descrever
efeitos de solvente em espectros de absorção [208–212]. Inclusive, encontra-se na literatura
trabalhos com membranas tratados por QM/MM [213–216].
O método QM/MM consiste em descrever explicitamente todo o meio que se deseja
estudar, tratando por mecânica quântica (QM) apenas a região do sistema em que os
fenômenos dependentes de elétrons acontecem e, todo o resto do sistema, por mecânica
molecular (MM). A energia e distribuição de carga da região tratada por MM é descrita
por funções clássicas de potencial e distribuição de carga. A energia das duas regiões do
sistema estão acopladas pelas mudanças nas cargas quânticas e no potencial eletrostático
das cargas e dipolos induzidos da região clássica e considerando as interações eletrostáticas
e de van der Waals entre os átomos das duas regiões. Assim, é incluı́do no cálculo o
efeito dielétrico e a polarização da região clássica e a orientação do momento de dipolo
do solvente. Os primeiros a descrever o tratamento de sistemas por QM/MM foram A.
Warshel e M. Levitt [217].
Para chegar na energia total do sistema não basta simplesmente somar as energias
das duas regiões. Termos acoplados devem ser considerados com atenção especial para
onde existem ligações quı́micas entre as regiões, principalmente se o limite entre elas inclui
ligações covalentes. Neste limite, existem átomos que devem ser considerados no cálculo de
energia das duas regiões. Existem vários formalismos para descrever a energia do sistema.
Neste trabalho, foi utilizada a forma aditiva, que é descrita pela Equação 3.107, em que
os subscritos indicam o tipo de cálculo e entre parênteses, a região do sistema que foi
considerada, sendo S o sistema inteiro, I, a região interna (tratada por QM), L são os
átomos que fazem a fronteira entre as duas regiões, também tratado por QM, e O a região
externa (tratada por MM) [218]. Para entender melhor L, ver as referências [218–220].
59
Equação 3.108 [218].
ÿ 1 ÿ ind soluto
ĤQM/M M = ĤQM/M
el
M + ĤQM/M M =
pol
qm V̂ (rm ) ≠ µ Ê (ra ) (3.111)
m 2 a a a
ÿ ÿ qm qn 1 ÿ ind ÿ qm (ra ≠ rm )
ĤM M = ĤM
el
M + ĤM M =
pol
≠ µ (3.112)
m n>m rmn 2 a a m |ra ≠ rm |3
60
no operador Fock por causa da dependência explicita dos dipolos induzidos na função de
onda QM. Além disso, em QM/MMPol não são consideradas contribuições de dispersão
e repulsão em HQM/M M e HM M como na maioria dos métodos QM/MM, descritos por
potenciais empı́ricos independentes dos graus de liberdade QM. Os dipolos induzidos
em cada sı́tio MM polarizável é dado pela Equação 3.113, em que a aproximação linear
foi assumida, negligenciando qualquer contribuição de caracter magnético em relação ao
campo de força, usando uma polarizabilidade isotrópica –a para cada ponto selecionado
na região externa. O termo Easolvente se refer ao campo elétrico total do solvente calculado
no sı́tio a e contem a soma das contribuições dos pontos de carga e momentos de dipolo
induzido no solvente. Ou seja, a Equação 3.113 deve ser resolvida interativamente, em
cada passo do SCF [221].
Ó Ô
a = –a (Ea
µind soluto
+ Esolvente
a q; µind ) (3.113)
µind = BE (3.114)
61
4 Objetivos
Este trabalho tem como objetivo geral estudar, utilizando métodos de quı́mica quântica,
derivados de porfirinas, ftalocianinas e clorinas, livres e complexadas com ı́ons de metais
de transição, e com diferentes ligantes, visando sua utilização em PDT.
Como objetivos especı́ficos realizar:
• uma avaliação de diversas propriedades como área, volume, dipolo, energia livre de
solvatação e dos estados excitados dos fotossensibilizadores;
64
F 4H F Al 4H
F Zn 4H F 1nh3
F 1nh3 4c F 2nh3
65
F 2nh3t F 3nh3
F 4nh3 F nh2ch3
F n3ch3 F 1oacc0
66
F n2acc0 F Zn ac2e
F Zn bc2e F Zn cd2e
F Zn c2ef F Zn c2eg
67
F Zn c2eh
Figura 15: Estruturas das ftalocianinas estudadas. Esfera cinza = carbono, branca = hidrogênio,
azul = nitrogênio, vermelha = oxigênio, rosa = alumı́nio, roxa = zinco e verde =
cloro.
membrana que mais contém moléculas de oxigênio dissolvidas, ajudando assim, a eficiência
da PDT [51].
A forma tradicional de simular essa interação é por cálculos de dinâmica molecular
(MD). Na literatura existem vários trabalhos que simulam membranas de bicamada lipı́dica
com diferentes pacotes computacionais de MD, inclusive a membrana POPG em questão
[61, 63, 64, 66–69]. Porém, é um método complexo pois a penetração da molécula induz
uma distorção na membrana, que descrita de forma discreta, assim como o solvente
que a circunda, torna o método custoso. Devido a essa complexidade e elevado custo
computacional deste método, foram feitos também cálculos com o modelo de solvatação
IEFPCM, para o qual a membrana foi descrita como um gradiente dielétrico, ao invés de
ser discreta, o que diminui consideravelmente o tempo de cálculo. Assim sendo, a quarta e
última parte deste trabalho, consistiu em simular a interação entre corante e membrana
de bicamada lipı́dica de duas formas, afim de obter o espectro UV-visı́vel do corante:
(I-a) por dinâmica molecular clássica, utilizando o pacote computacional AMBER12
[223], com o programa AmberTools13 [224] que inclui Lipid11 [200], campo de força pra
lipı́dio modular, descreve de forma mais robusta a dinâmica dos fosfolipı́dios (desenvolvido
para ser compatı́vel com outros campos de força do Amber). Nele está incluso também
parâmetros do General Amber Force Field (GAFF) [202] utilizado para o fotossensibilizador,
sendo que as cargas RESP do mesmo foram calculadas com HF/6-31g(d) no pacote
computacional Gaussian 09 (simulando solvente água por IEFPCM), pois, apesar de
ser um método custoso, esta é a abordagem de carga padrão aplicado nos campos de
68
P 4H P Al 4H
P Zn 4H P Si 2Cl
P 4fenil P 4b1cA
69
P 4b1cB P 4b2cAC
P 4b2cAT P 4b3cA
P 4b4cA P 4b4cB
70
P Si 2ch3 P Si 2fenil
Levulan Visudyne
Figura 16: Estruturas das porfirinas estudadas. Esfera cinza = carbono, branca = hidrogênio,
azul = nitrogênio, vermelha = oxigênio, rosa = alumı́nio, roxa = zinco, verde =
cloro e amarelo=silı́cio.
força GAFF e tem funcionado bem tanto para proteı́nas quanto para pequenas moléculas
[225, 226]. Alguns trabalhos na literatura utilizaram AMBER para estudos de membranas
de bicamada lipı́dica, inclusive a POPG [227–229]. Este método fornece, além de dados
sobre a penetração/estabilidade do corante na membrana, como fornece conformações para
realização de cálculos QM/MMPol (quantum mechanical/ molecular mechanical polarized)
[208, 230] que gerará o espectro UV-visı́vel. Os cálculos QM/MMPol permitem que exista
interação dos dipolos polarizáveis da mecânica molecular com a função de onda completa
da mecânica quântica, mantendo a conservação de energia da dinâmica molecular. Foram
selecionadas apenas dois corantes (F 2nh3 e P 4b2cAT) para realizar a interação com
membrana de bicamada lipı́dica e apenas um corante, F 2nh3 para realizar os cálculos
QM/MMPol, devido à complexidade e elevado custo computacional desses métodos.
O modelo da membrana de bicamada utilizada é aniônica, foi feita com 128 lipı́dios (64
em cada camada) do tipo palmitoil-oleil-fosfatidil-glicerol (POPG), sendo as dimensões da
membrana 6 x 6 x 4 nm e foram adicionados ao sistema moléculas de água e contra-ı́ons
utilizando gerador de input CHARMM-GUI Website Tool [231, 232]. Esses contra ı́ons
71
Foscan C 4H
C Zn 4H
Figura 17: Estruturas das clorinas estudadas. Esfera cinza = carbono, branca = hidrogênio,
azul = nitrogênio e roxa = zinco.
são 128 átomos de sódio como contra ı́ons da membrana e 2 átomos de cloreto para
neutralizar o fotossensibilizador, já que a molécula possui 2 cargas positivas. A hidratação
da membrana foi feita com 8078 moléculas de água, gerando uma espessura de 6 nm,
descrita pelo modelo TIP3P [193]. Para os fotossensibilizadores, previamente otimizados
com DFT/B3LYP/LanL2DZ, foram utilizadas as cargas RESP calculadas com HF/6-31g(d)
no pacote computacional Gaussian 09;
(I-b) Utilizando as conformações do sistema geradas pela metodologia “I-a”, é possı́vel
fazer um cálculo QM/MMPol, onde a molécula é tratada com mecânica quântica (QM, do
inglês: Quantum Mechanics) e, dentro de um raio de corte de 20 Å centrado na molécula,
os átomos são descritos por uma carga e uma polarizabilidade. Fora desse raio de corte, os
átomos são tratados por MM (Mecânica Molecular), ou seja, apenas por cargas parciais. A
partir disso, é extraı́do o espectro UV-visı́vel de várias conformações extraı́das da dinâmica
molecular;
(II) utilizando o modelo contı́nuo de solvatação (IEFPCM) para simular a membrana
por meio de um contı́nuo dielétrico que abrange toda a variação da constante dielétrica ao
longo da membrana fosfolipı́dica. Já com este método foi possı́vel realizar o estudo de
72
Etapa Parâmetros
Volume constante
Pressão não fixada
Minimização
Algorı́timo SHAKE = desligado
Cutoff = 10 Å
Número de passos = 2500; tempo de cada passo = 0,002 ps
Algorı́timo SHAKE = ligado
Banho termostático = Langevin
Temperatura inicial = 0 K; Temperatura final = 100 K
1 Aquecimento
Volume constante
Pressão não fixada
Cutoff = 10 Å
Lipı́deos fixados
Número de passos = 50000; tempo de cada passo = 0,002 ps
Algorı́timo SHAKE = ligado
Banho termostático = Langevin
Temperatura inicial = 100 K; Temperatura final = 300 K
2 Aquecimento
Pressão constante (condições periódicas de contorno)
Pressão anisotrópica simulada pelo barostato de Berendsen [164]
Cutoff = 10 Å
Lipı́deos fixados
Número de passos = 500000; Tempo de cada passo = 0,002 ps
Algorı́timo SHAKE = ligado
Banho termostático = Langevin
Produção Temperatura Constante = 300 K
Pressão Constante (condições periódicas de contorno)
Pressão anisotrópica simulada pelo barostato de Berendsen [164]
Cutoff = 10 Å
Tabela 2: Parâmetros utilizados em cada passo da dinâmica molecular realizada
interação com mais moléculas, pelo seu baixo custo computacional e complexidade.
Para as duas dinâmicas realizadas, em caixa de tamanho 90 Å3 , utilizando condições
periódicas de contorno, primeiramente foi feita uma minimização, com o intuito de
encontrar um mı́nimo de energia local, para definir as velocidades iniciais e remover
os “maus contatos”. Depois da minimização, foram realizadas 2 etapas de aquecimento,
mantendo os lipı́deos fixos, a primeira mantendo o volume constante e segunda mantendo
a pressão constante. O segundo aquecimento foi programado para aumentar a temperatura
lentamente até a temperatura necessária para realizar a etapa de produção. E, finalmente,
a etapa de produção foi realizada, somando 200 ns, sendo o tempo de cada passo 2 fs. Os
parâmetros de cada etapa estão na Tabela 2.
Além disso, a mesma metodologia (II) foi utilizada para simular a energia de emissão
dessas moléculas. Primeiramente, as moléculas foram otimizadas no estado tripleto. Depois
foram calculadas as energias dos estados tripleto excitado e singleto fundamental (estados
73
envolvidos no processo P da Figura 1) mantendo a geometria do estado tripleto. A
diferença de energia entre esses dois estados então, seria correspondente a energia de
emissão dessas moléculas. Esse cálculo foi realizado para cada ponto da membrana feita
com PCM para algumas das moléculas estudadas nesse trabalho.
74
6 Resultados e Discussão
Molécula
Metodologia
F 4H F Zn 4H
1 B3LYP/6-31g(d) 4,77 x 10≠02 4,61 x 10≠03
2 B3PW91/6-31g(d) 4,75 x 10≠02 8,24 x 10≠03
3 M06L/6-31g(d) 4,78 x 10≠02 4,14 x 10≠03
4 B3LYP/DGDZVP 4,83 x 10≠02 1,07 x 10≠02
5 B3LYP/LanL2DZ 5,90 x 10≠02 3,83 x 10≠02
6 PM6 9,10 x 10≠02 6,49 x 10≠02
Tabela 3: Valores de RMDS gerados pela sobreposição das geometrias calculadas pelas me-
todologias (funcional/base) de 1 à 6, com a estrutura cristalográfica para F 4H e
F Zn 4H
Qx Qy Soret Qy -Soret
Metodologia
⁄máx (nm) Erro (%) ⁄máx (nm) Erro (%) ⁄máx (nm) Erro (%) ⁄máx (nm) Erro (%)
1 B3LYP/6-31g(d) 622,0 9,7 621,5 5,0 335,5 0,6 286,0 9,6
5 B3LYP/LanL2DZ 638,0 7,4 637,0 2,6 342,5 -1,5 294,5 7,0
6 PM6 644,0 6,4 630,0 3,7 333,0 1,3 297,0 6,2
Experimental 689,0 - 654,0 - 337,5 - 283,0 -
Tabela 4: Comprimentos de onda máximo (⁄máx ) de absorção (nm) referentes às principais
bandas do espectro UV-visı́vel (Qx , Qy e Soret), diferença entre os ⁄máx das bandas
Qy e Soret e erros (%) destes em relação aos dados experimentais (em THF) para
F 4H calculados com B3LYP/6-31g(d), com geometrias obtidas pelos métodos 1, 5 e
6, em THF (IEFPCM)
Qx Soret Qx -Soret
Metodologia
⁄máx (nm) Erro (%) ⁄máx (nm) Erro (%) ⁄máx (nm) Erro (%)
1 B3LYP/6-31g(d) 614,5 7,9 333,0 3,1 281,5 13,0
5 B3LYP/LanL2DZ 630,0 5,5 342,0 0,4 288,0 11,0
6 PM6 657,5 1,4 337,0 1,9 320,5 0,9
Experimental 667,0 - 343,5 - 323,5 -
Tabela 5: Comprimentos de onda máximo (⁄máx ) de absorção (nm) referentes às principais
bandas do espectro UV-visı́vel (Qx e Soret), a diferença entre os ⁄máx das bandas
Qx e Soret e erros (%) destes em relação aos dados experimentais (em THF) para
F Zn 4H calculados com B3LYP/6-31g(d), com geometrias obtidas pelos métodos 1,
5 e 6, em THF (IEFPCM)
76
F 4H F Zn 4H
Metodologias Def. ang. N–H Est. C–N Est. C–N Est. C–C
‹ (cm )≠1
‹ (cm≠1 ) Erro ‹ (cm ) ‹ (cm≠1 ) Erro
≠1
Tabela 6: Número de onda (‹) em cm≠1 do principal pico do espectro vibracional experimentais
(fase sólida) e calculados (fase gasosa) por diversas metodologias (funcional/base),
de F 4H e F Zn 4H e erros (%) destes em relação aos dados experimentais
números de onda, e o erro (%) foi calculado a partir destes. São apresentados na Tabela
6 os números de onda (‹) experimentais (obtidas em pastilhas de KBr) e calculados
(fase gasosa) com as mesmas metodologias empregadas para obter as geometrias: B3LYP,
B3PW91, M06L utilizando a base 6-31g(d), B3LYP/ DGDZVP e B3LYP/LanL2DZ, em
fase gasosa.
Tanto para F 4H quanto para F Zn 4H, os menores erros de números de onda calculados
em relação ao experimental, foram obtidos pela metodologia B3LYP/LanL2DZ, sendo
assim, este foi o método escolhido para obter o espectro vibracional doravante. Entre-
tanto, é importante ressaltar que o pico com carácter muito forte obtido com os cálculos
B3LYP/LanL2DZ (1038 cm≠1 ), corresponde ao estiramento C–N no anel piridı́nico e não
deformação angular N–H, como Seoudi propõe [?]. Esta última vibração foi encontrada
nos cálculos em 889 cm≠1 . O mesmo acontece com F Zn 4H, o pico com carácter forte,
que Seoudi atribui a um estiramento C–N do anel pirrol, nos cálculos mostram uma
deformação angular C–C no anel benzil, em 1079 cm≠1 . De qualquer forma a vibração
escolhida para a caracterização de cada molécula, feita por Seoudi, 2005 entre outros, é
considerada válida considerando a diferença entre os espectros experimentais das duas
moléculas (Figura 18), e portanto o trabalho vai focar em reproduzir tal vibração. Na
Figura 19 e 20, são apresentados os espectros experimental e calculado com a metodologia
escolhida sobrepostos para F 4H e F Zn 4H, respectivamente.
Visando estabelecer apenas uma metodologia para otimizar todas as moléculas a
serem estudadas no presente trabalho, foi escolhida a metodologia 5 (B3LYP/LanL2DZ),
considerando que: com as metodologias de 1 a 4 o custo computacional é mais elevado;
dentre as metodologias 5 e 6, a última foi a que forneceu maiores valores de RMSD e
espectro vibracional com maiores erros; a base LanL2DZ inclui efeitos relativı́sticos, o que
pode ser importante na otimização de ftalocianinas complexadas com metais de transição.
77
Figura 18: Espectros vibracionais experimentais (em pastilha de KBr) de F 4H e F Zn 4H
sobrepostos.
78
Figura 19: Espectros vibracionais experimental (em pastilha de KBr) e calculado
(B3LYP/LanL2DZ, em fase gasosa), de F 4H, sobrepostos.
Qx Qy Soret Qy ≠ Soret
Metodologia
⁄máx Erro ⁄máx Erro ⁄máx Erro ⁄máx Erro
B3LYP/6-31g(d) 622,0 9,7 621,5 5,0 335,5 0,6 286,5 9,6
B3PW91/6-31g(d) 616,0 10,6 615,0 6,0 341,0 -1,0 274,0 13,4
M06L/6-31g(d) 629,0 8,7 628,5 3,9 314,0 7,4 314,5 0,6
OLYP/6-31g(d) 644,0 6,5 642,0 1,8 330,0 2,1 311,5 1,6
B3LYP/DGDZVP 634,0 8,0 634,0 3,1 348,0 -3,1 286,0 9,6
B3LYP/LanL2DZ 628,0 8,9 620,0 5,2 359,0 -6,4 261,0 17,5
ZINDO 853,0 -24,3 798,0 -28,1 324,0 4,7 474,0 -67,5
OLYP/6-31g(d) com
663,5 3,8 661,0 -1,1 374,0 -10,8 287,0 9,3
geometria B3LYP/LanL2DZ
Experimental 689,0 - 654,0 - 337,5 - 316,5 -
às bandas Qx (147æ 148 e 147æ 149) e Soret (134æ 148 e 134æ 149) para F Zn 4H,
calculados com o método escolhido, em THF.
Em relação às Figuras 21 e 22, observa-se que os orbitais HOMO (highest occupied
molecular orbital), LUMO (lowest unoccupied molecular orbital) e LUMO+1 apresentam
caráter mais anti-ligante que os orbitais mais internos (que contribuem para a banda
Soret). O HOMO é mais distribuı́do por toda a molécula, embora tenha participação
apenas de carbonos e LUMO e LUMO+1 apresentam participação de 2 nitrogênios cada
79
Figura 20: Espectros vibracionais experimental (em pastilha de KBr) e calculado
(B3LYP/LanL2DZ, em fase gasosa), de F Zn 4H, sobrepostos.
Q Soret Q ≠ Soret
Metodologia
⁄máx Erro ⁄máx Erro ⁄máx Erro
B3LYP/6-31g(d) 614,5 7,9 333,0 3,2 281,5 5,0
B3PW91/6-31g(d) 609,0 9,5 331,0 3,8 278,0 6,3
M06L/6-31g(d) 622,5 7,1 371,0 -7,4 251,5 17,5
OLYP/6-31g(d) 636,0 4,9 362,0 -5,1 274,0 7,8
B3LYP/DGDZVP 626,5 6,5 338,0 1,6 288,5 2,4
B3LYP/LanL2DZ 616,0 8,3 347,0 -1,0 269,0 9,8
ZINDO 822,0 -19,3 324,0 4,0 498,0 -57,3
OLYP/6-31g(d) com
654,0 2,0 336,0 2,2 318,0 7,1
geometria B3LYP/LanL2DZ
Experimental 667,0 - 343,5 - 295,5 -
além dos carbonos. Além disso, LUMO e o LUMO+1 apresentam participação de átomos
ao longo de um eixo (sendo que um deles tem seus orbitais no eixo x e outro no eixo
y). Como esperado, sempre as duas transições referentes à banda Q, são HOMO æ
LUMO e/ou HOMO æ LUMO+1, e como os orbitais virtuais estenderam-se cada um
ao longo de um eixo, tem-se as bandas Qx e Qy . A banda Soret envolve um orbital
preenchido bastante interno (120 para F 4H e 134 para F Zn 4H), com caráter ligante,
envolvendo apenas carbonos, que no processo de transição eletrônica, passa a concentrar a
80
120 133 (HOMO) 134 (LUMO) 135
Figura 21: Orbitais envolvidos nas transições relativos as bandas Qx , Qy e Soret, calculados
com OLYP/6-31g(d) em THF (IEFPCM) para F 4H.
Figura 22: Orbitais envolvidos nas transições relativos as bandas Qx e Soret, calculados com
OLYP/6-31g(d) em THF (IEFPCM) para F Zn 4H.
carga ao longo de um eixo que corresponde a LUMO ou LUMO+1. Com base no que foi
exposto até aqui, o procedimento adotado para as demais moléculas foi otimização com
DFT/B3LYP/LanL2DZ, assim como cálculo de espectros infravermelho, e simulação dos
espectros de absorção com TD-DFT/OLYP/6-31g(d).
81
o quanto o meio interfere nas propriedades eletrônicas do sistema. Foi estudada a variação
de solubilidade de acordo com a mudança dos substituintes e do metal complexado. No
Apêndice A, a Tabela 24 mostra a G em fase gasosa ( Gf g ), o G em solução aquosa
( Gag ) e a diferença entre essas energias ( Gsolv ) para todas as moléculas estudadas.
Para auxiliar na análise, foi feito o gráfico apresentado na Figura 23 com o Gsolv de cada
molécula. Como não foi possı́vel convergir alguns cálculos referentes a molécula F ac2e,
os dados de Gsolv da mesma não foram obtidos, portanto a mesma não é incluı́da em
nenhum dos gráficos que envolve Gsolv .
82
83
ΔGsolv'
0,000%
!0,700%
!0,650%
!0,600%
!0,550%
!0,500%
!0,450%
!0,400%
!0,350%
!0,300%
!0,250%
!0,200%
!0,150%
!0,100%
!0,050%
Figura 23:
F_4H%(0)%
P_4H%(0)%
Livres%
C_4H%(0)%
F_Al_4H%(0)%
P_Al_4H%(0)%
F_Zn_4H%(0)%
P_Zn_4H%(0)%
C_Zn_4H%(0)%
Complexadas%
P_Si_4H%(0)%
P_Si_2ch3%(0)%
P_Si_2fenil%(0)%
F_1oacc0%(0)%
F_n2acc0%(0)%
F_Zn_bc2e%(0)%
F_Zn_cd2e%(0)%
F_Zn_c2eg%(0)%
F_Zn_c2eh%(0)%
P_4fenil%(0)%
Foscan%(0)%
Levulan%(0)%
Visudyne%(0)%
F_1nh3%(1)%
F_1nh3_4c%(1)%
F_2nh3%(2)%
Com%subsOtuintes%
F_2nh3t%(2)%
F_3nh3%(3)%
F_4nh3%(4)%
F_nh2ch3%(1)%
F_n3ch3%(1)%
P_4b1cA%(1)%
P_4b1cB%(1)%
P_4b2cAC%(2)%
P_4b2cAT%(2)%
P_4b3cA%(3)%
P_4b4cA%(4)%
P_4b4cB%(4)%
Gsolv (eV) de todas as moléculas estudadas. Colunas em azul representam ftalocianinas, as em vermelho, porfirinas e as em verde,
O gráfico da Figura 23 mostra claramente que a adição de grupos carregados tem efeito
maior sobre a solubilidade em água do que a complexação por ı́ons metálicos ou a adição
de substituintes neutros. Entretanto, como o modelo utilizado não explicita as ligações de
hidrogênio, a Gsolv calculada para compostos com grupos carboxı́licos, como no caso da
F n2oacc0, deve estar subestimada.
Mais especificamente analisando as contribuições das diferentes mudanças nas estruturas
dos compostos:
• dentre as famı́lias estudadas não observa-se muita diferença no Gsolv , sendo que a
ftalocianina livre o possui mais elevado, seguido da porfirina livre e depois a clorina
livre. O mesmo acontece quando as mesmas complexadas com os ı́ons metais são
comparadas entre si;
• a adição de grupos fenila à porfirina livre (P 4fenil) não acusa diminuição drástica
de solvatação;
• dentre os compostos neutros, observa-se que os que possuem ligantes com éter são
os mais favorecidos (F Zn bc2e, F Zn cd2e, F Zn c2ef, F Zn c2eg e F Zn c2eh) ,
seguido do composto Foscan, que possui 4 anéis fenı́licos com uma substitução
álcool em cada um deles, e depois os compostos com ácidos carboxı́licos (F 1oacc0 e
F n2oacc0), Levulan que possui como ligante duas cadeias carbônicas terminadas
em álcool e o Visudyne, com ligantes esteres;
• as posições cis e trans não mostram diferença considerável no Gsolv , assim como as
A e B do ligante metil no anel piridı́nico nas moléculas P 4b1cA, P 4b1cB, P 4b4cA
e P 4b4cB.
84
Propriedades como dipolo, volume e área relacionam a estrutura molecular com a
atividade biológica da mesma [233], podendo indicar formas de interação desses fotos-
sensibilizadores com a membrana celular. Muito utilizadas em quı́mica medicinal, essas
propriedades ajudam na sı́ntese e planejamento de novos fármacos. Tanto o volume quanto
a área grid foram calculados pelo método em que é utilizada uma sonda do tamanho de
uma molécula de solvente (no caso, água) para montar uma espécie de malha (grid) ao
redor da molécula) e a partir dessa malha são calculados o volume e área grid. Já a área
aproximada é calculada pelo método de Hasel et al. [234] onde o raio de van der Waals é
tomado como base para o cálculo da área.
Ao que importa neste trabalho, sabe-se que quanto menor o dipolo da molécula (obtido
por cálculo DFT/B3LYP/LanL2DZ) e maior seu volume (assim como a área), maior a
capacidade da mesma interagir com meios hidrofóbicos como as membranas. Na Tabela 9
estão os valores de carga total da molécula, área grid, área aproximada, volume e dipolo
para todas as moléculas estudadas.
Com intuito de correlacionar propriedades, foram construı́dos gráficos correlacionando
o Gsolv com o volume (Figura 24) e com o dipolo (Figura 25 e 26).
85
Molécula Carga Área grid Área aproximada Volume Dipolo
F 4H 0 480,9 731,3 1328,5 0,0023
F Al 4H 0 501,3 758,3 1390,5 3,7214
F Zn 4H 0 482,7 736,2 1331,0 0,0009
F 1nh3 1 572,6 829,2 1482,8 33,3466
F 1nh3 4c 1 645,1 882,5 1592,4 44,6118
F 2nh3 2 665,2 909,6 1637,4 41,0710
F 2nh3t 2 665,5 909,3 1637,6 0,0081
F 3nh3 3 759,0 991,5 1792,8 26,0069
F 4nh3 4 1086,2 854,0 1947,7 0,0046
F nh2ch3 1 663,2 874,3 1579,6 32,6599
F n3ch3 1 708,3 886,8 1618,9 32,0383
F 1oacc0 0 594,0 835,3 1491,4 2,7238
F n2acc0 0 647,5 909,6 1664,3 1,8813
F Zn ac2e 0 883,2 1091,4 1967,8 3,3644
F Zn bc2e 0 920,9 1126,2 2027,7 3,0185
F Zn cd2e 0 936,9 1144,8 2055,7 2,0904
F Zn c2ef 0 964,4 1209,3 2179,2 6,3298
F Zn c2eg 0 987,1 1238,2 2210,7 7,5199
F Zn c2eh 0 1062,4 1297,0 2326,8 2,6030
P 4H 0 345,1 497,1 871,6 0,0000
P Al 4H 0 356,1 523,2 928,4 4,4625
P Zn 4H 0 348,5 505,8 870,0 0,0005
P Si 2Cl 0 360,8 521,8 955,7 0,0030
P 4fenil 0 638,3 905,2 1683,9 0,0001
P 4b1cA 1 690,7 928,1 1717,5 16,4964
P 4b1cB 1 684,2 929,2 1723,3 17,9573
P 4b2cAC 2 742,1 953,4 1761,3 25,0918
P 4b2cAT 2 742,4 949,4 1764,5 0,0352
P 4b3cA 3 796,2 973,4 1813,2 17,9188
P 4b4cA 4 851,5 1008,8 1871,3 0,0483
P 4b4cB 4 835,1 1023,9 1881,3 0,4345
P Si 2ch3 0 349,2 548,2 1007,6 0,9967
P Si 2fenil 0 425,5 683,9 1305,5 0,0160
Levulan 0 735,8 826,1 1559,6 4,0947
Visudyne 0 875,1 961,8 1827,7 3,3767
Foscan 0 707,4 957,1 1774,7 3,8808
C 4H 0 506,2 355,2 888,3 2,3260
C Zn 4H 0 504,7 358,2 886,4 2,3125
Tabela 9: Carga, área grid (Å2 ), área aproximada (Å2 ), volume (Å3 ) e dipolo (Debye) de todas
as moléculas estudadas.
86
2400,00%
F_Zn_c2eh%
F_Zn_c2eg%
2200,00%
F_Zn_c2ef%
F_Zn_cd2e%
F_Zn_bc2e%
2000,00%
F_4nh3%
P_4b4cB% P_4b4cA%
P_4b3cA% Visudyne%
1800,00% F_3nh3% Foscan%
Volume'(Å3)'
P_4b2cAC% P_4b2cAT%
P_4b1cA% P_4b1cB%
P_4fenil%
F_2nh3t% F_2nh3% F_n2acc0%
F_n3ch3%
1600,00%
F_1nh3_4c% Levulan%
F_nh2ch3%
F_1oacc0%
87
F_1nh3%
F_Al_4H%
1400,00%
F_Zn_4H% F_4H%
P_Si_2fenil%
1200,00%
P_Si_2ch3%
1000,00%
P_Si_4H% P_Al_4H%
C_Zn_4H% C_4H%
P_Zn_4H% P_4H%
800,00%
G0,750% G0,650% G0,550% G0,450% G0,350% G0,250% G0,150% G0,050% 0,050%
ΔGsolv'
88
50,00)
45,00) F_1nh3_4c)
F_2nh3)
40,00)
35,00)
F_1nh3) F_nh2ch3)
F_n3ch3)
30,00)
Dipolo&(Debye)&
F_3nh3) P_4b2cAC)
25,00)
89
20,00)
P_4b3cA) P_4b1cB)
P_4b1cA)
15,00)
10,00)
F_Zn_c2eg)
F_Zn_c2ef)
5,00) F_Zn_bc2e)
F_Zn_c2eh)
P_4b4cB)
F_4nh3) P_4b4cA) F_2nh3t) P_4b2cAT) F_Zn_cd2e)
0,00)
:0,700) :0,600) :0,500) :0,400) :0,300) :0,200) :0,100) 0,000)
ΔGsolv&
Figura 25: Gsolv (eV) e dipolo (Debye) de todas as moléculas estudadas. As moléculas no quadro vermelho estão representadas na Figura 26.
20,00%
18,00% P_4b1cB%
16,00% P_4b1cA%
14,00%
12,00%
Dipolo&(Debye)&
10,00%
90
8,00%
6,00%
P_Al_4H%
Foscan% Levulan%
4,00%
Visudyne% F_Al_4H%
F_1oacc0% C_Zn_4H%
F_n2acc0% C_4H%
P_Si_2fenil%
2,00%
P_Zn_4H%
P_Si_4H%
P_4fenil%
P_Si_2ch3%
ΔGsolv&
Figura 26: Gsolv (eV) e dipolo (Debye) das moléculas no quadro vermelho da Figura 25.
Analisando os gráficos das Figuras 25 e 26, observa-se que assim como para o volume,
a variação do Gsolv em relação ao dipolo é muito baixa frente à adição de cargas. As
moléculas neutras ficam agrupadas, sem apresentar muita variação no dipolo e no Gsolv .
Já entre as carregadas é importante ressaltar que, mesmo as moléculas com dipolo mais
alto como F 1nh3 4c e F 2nh3, não são as mesmas que apresentam maior Gsolv , e sim as
mais carregadas como F 4nh3 e P 4b4cB, que tem dipolo praticamente nulo.
91
#
e#
e#
h#
e#
g#
4c
f#
#
#
h3
c2
d2
2e
2e
2e
c2
c0
c0
H#
3#
#
3t#
3_
4H
3#
3#
3#
3#
_b
_4
2c
ch
_a
ac
ac
_c
_c
_c
_c
h
h
Al_
nh
1o
1n
1n
2n
2n
3n
4n
n3
n2
Zn
Zn
Zn
Zn
Zn
Zn
Zn
H#
F4
F_
F_
F_
F_
F_
F_
F_
F_
F_
F_
F_
F_
F_
F_
F_
F_
F_
F_
0,080#
0,070#
0,060#
0,050#
ΔH>L#(u.a.)##
0,040#
92
0,030#
0,020#
0,010#
0,000#
Fase#gasosa# Clorofórmio# Água#
Figura 27: H≠L (u.a.) de todas as ftalocianinas estudadas. As colunas em vermelho representam dados simulados em fase gasosa, as verdes em
clorofórmio e as azuis em água.
#
nil
#
#
h3
AC
AT
e#
H#
H#
A#
A#
A#
B#
B#
#
2fe
4H
#
4H
nil
yn
2c
lan
n#
_4
_4
1c
1c
2c
2c
3c
4c
4c
#
#
4fe
Al_
sca
ud
Si_
Si_
Si_
4H
4H
4b
4b
4b
4b
4b
4b
4b
Zn
Zn
vu
Vis
C_
C_
P_
P_
P_
P_
P_
P_
P_
P_
P_
P_
P_
P_
P_
P_
Fo
Le
0,080#
0,070#
0,060#
0,050#
ΔHHL#(u.a.)##
0,040#
93
0,030#
0,020#
0,010#
0,000#
Fase#gasosa# Clorofórmio# Água#
Figura 28: H≠L (u.a.) de todas as porfirinas e clorinas estudadas. As colunas em vermelho representam dados simulados em fase gasosa, as verdes
em clorofórmio e as azuis em água.
De acordo com os gráficos das Figuras 27 e 28, primeiramente observa-se que pratica-
mente não há variação do H≠L independente da complexação com metais e/ou adição de
ligantes. Para as porfirinas, observa-se que a adição de ligantes fenı́licos e piridı́nicos faz
com que esse H≠L diminua, mais acentuadamente em clorofórmio. Observa-se também
que a complexação tende a aumentar ligeiramente o H≠L . Não há diferença significativa
em relação a variação das posições dos ligantes cis e trans, porém um pouco considerável
de P 4b1cA para P 4b1cB e de P 4b4cA para P 4b4cB.
Foram analisados também os gráficos de densidade eletrônica dos orbitais de fronteira
desses compostos, que participam das transições referentes a banda Qx , com o intuito
de saber quais partes das moléculas estão contribuindo para as mesmas. A adição de
ligantes e do metal complexado pode acarretar na alteração do rendimento quântico do
fotossensibilizador. Segundo o Modelo de Quatro Orbitais de Gouterman, a porfirina livre
tem simetria D2h , assim como a ftalocianina. É possı́vel observar tal simetria nos gráficos
obtidos dos orbitais que participam das transições das bandas Q, HOMO-1, HOMO,
LUMO e LUMO+1, para a molécula P 4H, representados na Figura 29.
Figura 29: Orbitais de fronteira HOMO-1, HOMO, LUMO e LUMO+1 de porfirina livre
(P 4H), calculados com OLYP/6-31g(d) em THF (IEFPCM).
Figura 30: Orbitais de fronteira HOMO-1, HOMO, LUMO e LUMO+1 da porfirina complexada
com Zn (P Zn 4H), calculados com OLYP/6-31g(d) em THF (IEFPCM).
94
HOMO LUMO
Figura 31: Orbitais de fronteira HOMO e LUMO da porfirina complexada com Al (P Al 4H),
calculados com OLYP/6-31g(d) em THF (IEFPCM).
o Zinco nunca participa dos orbitais HOMO e HOMO-1, porém participam um pouco
dos orbitais LUMO, apesar de a presença de metais não alterar significativamente o
H≠L apresentado na Figura 28. Para todas as moléculas com ligantes, não se observa
Figura 32: Orbitais de fronteira HOMO, LUMO e LUMO+1 de porfirina P 4b1cA, calculados
com OLYP/6-31g(d) em água (IEFPCM).
HOMO LUMO+1
Figura 33: Orbitais de fronteira HOMO e LUMO da porfirina P 4fenil, calculados com OLYP/6-
31g(d) em água (IEFPCM).
95
com metal, de simetria D4h , apresentam os orbitais LUMO distribuı́dos por toda a molécula.
Nos dois casos, sem participação de ligantes. Tanto para os orbitais HOMO quanto para
LUMO, os metais de complexação não participam, porém o H≠L de energia chega a cair
pela metade, em relação as não complexadas.
Figura 34: Orbital de fronteira HOMO para as moléculas F 1oacc0 e F n2acc0, calculados com
OLYP/6-31g(d) em água (IEFPCM).
Já as clorinas, apesar de serem similares as porfirinas, não apresentam a mesma simetria,
sendo C2v , tanto as livres quanto as complexadas, por terem uma dupla de um dos anéis
saturada. Essa simetria causa uma menor simetria nos orbitais e menor degenerêscencia
das transições referentes a banda Q, comparado com a porfirina, como pode-se observar
nas Figuras 35 e 36. Observando os dados de orbitais de fronteira das clorinas, o Zinco não
participa de nenhum dos orbitais analisados, assim como o H≠L não altera com relação a
clorina não complexada.
Figura 35: Orbitais de fronteira HOMO-1, HOMO, LUMO e LUMO+1 de porfirina complexada
com Zn (P Zn 4H), calculados com OLYP/6-31g(d) em água (IEFPCM).
Figura 36: Orbitais de fronteira HOMO-1, HOMO, LUMO e LUMO+1 de porfirina complexada
com Zn (P Zn 4H), calculados com OLYP/6-31g(d) em água (IEFPCM).
96
Experimentalmente, já se sabia que a geração de oxigênio singleto das porfirinas e
clorinas diminui com a complexação por metal por deslocar a densidade eletrônica para
o mesmo. Em porfirinas complexadas com Zinco II, o rendimento de oxigênio singleto
chega a cair quatro vezes em relação a porfirina livre, porém o mesmo não acontece com
as ftalocianinas que ficam mais fotoestáveis e mais geradoras de oxigênio singleto quando
complexadas (com Alumı́nio e Zinco II) [87, 235, 236].
Segundo o Teorema de Koopmans, a energia necessária para retirar um elétron de uma
molécula é conhecida como potencial de ionização e é dado pela Equação 6.1, onde N é o
número de elétrons. Em simulação computacional esse teorema pode ser utilizado dentro
da teoria da DFT e corresponde ao oposto da energia do orbital HOMO. Pelos Gráficos
37 e 38 é possı́vel analisar a energia de HOMO das ftalocianinas, porfirinas e clorinas
estudadas, respectivamente.
97
%
e%
e%
h%
e%
g%
4c
f%
%
%
h3
c2
d2
2e
2e
2e
c2
c0
c0
H%
3%
%
3t%
3_
4H
3%
3%
3%
3%
_b
_4
2c
ch
_a
ac
ac
_c
_c
_c
_c
h
h
Al_
nh
1o
1n
1n
2n
2n
3n
4n
n3
n2
Zn
Zn
Zn
Zn
Zn
Zn
Zn
H%
F4
F_
F_
F_
F_
F_
F_
F_
F_
F_
F_
F_
F_
F_
F_
F_
F_
F_
F_
0,000%
!0,050%
!0,100%
!0,150%
!0,200%
EH%(u.a.)%%
!0,250%
98
!0,300%
!0,350%
!0,400%
!0,450%
!0,500%
Figura 37: Energia do orbital HOMO (u.a.) de todas as ftalocianinas estudadas. As colunas em vermelho representam dados simulados em fase
gasosa, as verdes em clorofórmio e as azuis em água.
%
nil
%
%
h3
AC
AT
e%
H%
H%
A%
A%
A%
B%
B%
%
2fe
4H
%
4H
nil
yn
2c
l an
n%
_4
_4
1c
1c
2c
2c
3c
4c
4c
%
%
4fe
Al_
sca
ud
Si_
Si_
Si_
4H
4H
4b
4b
4b
4b
4b
4b
4b
Zn
Zn
vu
Vis
C_
C_
P_
P_
P_
P_
P_
P_
P_
P_
P_
P_
P_
P_
P_
P_
Fo
Le
!0,050%
!0,150%
EH%(u.a.)%%
!0,250%
!0,350%
99
!0,450%
!0,550%
Fase%gasosa% Clorofórmio% Água%
Figura 38: Energia do orbital HOMO (u.a.) de todas as porfirinas e clorinas estudadas. As colunas em vermelho representam dados simulados em
fase gasosa, as verdes em clorofórmio e as azuis em água.
Analisando os gráficos das Figuras 37 e 38, primeiramente se observa como a energia
de HOMO é diminuı́da quando em fase gasosa, principalmente quando as moléculas estão
carregadas, mostrando a importância de incluir efeitos do solvente nas simulações. Percebe-
se também, para as moléculas carregadas, uma queda nessa energia quando o solvente é o
clorofórmio em relação à água. Para as demais moléculas não há variação considerável na
energia de HOMO em relação aos solventes, mesmo com diferentes metais complexados e
ligantes adicionados. De qualquer forma, não se observa mudança considerável na energia
de HOMO em água em relação às complexações e substituições, porém uma ligeira queda
nas moléculas carregadas em clorofórmio.
Chama-se de afinidade eletrônica a capacidade da molécula de receber um elétron,
que corresponde à energia do orbital LUMO. Pelos Gráficos 39 e 40 é possı́vel analisar a
energia de LUMO das ftalocianinas, porfirinas e clorinas estudadas, respectivamente.
100
c%
e%
e%
h%
e%
g%
f%
%
%
_4
h3
c2
d2
2e
2e
2e
c2
c0
c0
H%
3%
%
t%
4H
%
h3
h3
h3
h3
h3
h3
_b
_4
2c
ch
_a
ac
ac
_c
_c
_c
_c
Al_
nh
1o
1n
1n
2n
2n
3n
4n
n3
n2
Zn
Zn
Zn
Zn
Zn
Zn
Zn
H%
F4
F_
F_
F_
F_
F_
F_
F_
F_
F_
F_
F_
F_
F_
F_
F_
F_
F_
F_
0,000%
!0,050%
!0,100%
!0,150%
!0,200%
EL%(u.a.)%%
101
!0,250%
!0,300%
!0,350%
!0,400%
!0,450%
!0,500%
Fase%gasosa% Clorofórmio% Água%
Figura 39: Energia do orbital LUMO (u.a.) de todas as ftalocianinas estudadas. As colunas em vermelho representam dados simulados em fase
gasosa, as verdes em clorofórmio e as azuis em água.
%
nil
%
%
h3
AC
AT
e%
H%
H%
A%
A%
A%
B%
B%
%
2fe
4H
%
4H
nil
yn
2c
lan
n%
_4
_4
1c
1c
2c
2c
3c
4c
4c
%
%
4fe
Al_
sca
ud
Si_
Si_
Si_
4H
4H
4b
4b
4b
4b
4b
4b
4b
Zn
Zn
vu
Vis
C_
C_
P_
P_
P_
P_
P_
P_
P_
P_
P_
P_
P_
P_
P_
P_
Fo
Le
!0,050%
!0,150%
EL%(u.a.)%%
!0,250%
102
!0,350%
!0,450%
!0,550%
Fase%gasosa% Clorofórmio% Água%
Figura 40: Energia do orbital LUMO (u.a.) de todas as porfirinas e clorinas estudadas. As colunas em vermelho representam dados simulados em
fase gasosa, as verdes em clorofórmio e as azuis em água.
Assim como para a energia de HOMO, observa-se a mesma tendência em relação a
adição de solvente e a constante dielétrica do solvente em questão. As moléculas carregadas
são mais sensı́veis à presença do dielétrico, principalmente para solventes menos polares.
Esse comportamento em relação as energias de HOMO e LUMO era esperado já que pela
análise dos orbitais de fronteira, praticamente não há participação dos ligantes e baixa
participação dos metais complexados.
Nesta seção, apenas as bandas Q serão analisadas, já que a banda Soret se encontra
fora da janela terapêutica de PDT para câncer. Visando também melhor representar
o meio em que a terapia ocorre, os solventes clorofórmio e água foram escolhidos pela
similariedade entre suas constantes dielétrias e os meio interno e externo, respectivamente,
da bicamada lipı́dica. Como escolhido no estudo metodológico apresentado anteriormente,
todos os cálculos de espectros de absorção daqui em diante foram calculados da seguinte
forma: otimização da geometria molecular com DFT/B3LYP/LanL2DZ seguido do cálculo
TD-DFT/OLYP/6-31g(d), e quando com solvente, o modelo IEFPCM foi aplicado. Nas
Tabelas 10 e 11 são apresentados os dados gerados para espectro de absorção em clorofórmio
e água, respectivamente, para todas as ftalocianinas estudadas.
103
Dipolo de transição
Molécula Carga Energia ⁄máx (nm) f Orbitais Contribuição
x y z
1,87 665 0,4274 133-135 0,58431 -0,0582 -0,0013 0,0000
F 4H 0
1,87 662 0,4979 133-134 0,58423 0,0015 -3,2933 0,0000
147-148 -0,1471 3,0604 -0,4500
F Al 4H 0 1,81 686 0,4249 0,68856
147-149 -0,0222 0,4495 3,0639
147-148 -3,2118 -0,0582 0,0000
F Zn 4H 0 1,89 657 0,4774 0,68964
147-149 0,0580 -3,2117 0,0000
1,85 672 0,4570 145-146 0,49316 -3,0554 0,08679 0,1314
F 1nh3 1
1,89 657 0,4913 145-147 0,49352 -0,8822 -3,1367 0,0475
1,86 668 0,4592 153-155 0,58415 3,1002 0,6211 0,3210
F 1nh3 4c 1
1,87 662 0,4905 153-154 0,58364 0,6458 -3,2032 0,0915
1,87 665 0,4711 157-159 0,46679 2,9633 -1,2200 -0,1998
F 2nh3 2
1,87 663 0,4963 157-158 0,46732 -1,4500 -2,9542 -0,075
1,82 679 0,4807 157-158 0,58300 3,1615 0,8513 -0,1808
F 2nh3t 2
1,90 654 0,4984 157-159 0,58234 0,8572 -3,1604 -0,0668
1,86 667 0,5401 169-170 0,58565 3,3391 0,8303 -0,1345
F 3nh3 3
1,88 659 0,4551 169-171 0,58398 0,8187 -3,0296 -0,1690
1,87 664 0,4838 181-183 0,66935 -3,1208 -0,8963 -0,1893
F 4nh3 4
1,87 662 0,5518 181-182 0,68191 -1,0374 3,3063 -0,1078
1,85 672 0,4598 153-154 0,43522 -3,0499 0,9147 0,1661
F nh2ch3 1
1,89 657 0,4914 154-155 0,43540 0,9351 3,1223 -0,0640
1,85 672 0,4638 157-158 0,44553 -3,0575 0,9214 0,2421
F n3ch3 1
1,89 657 0,4933 157-159 0,44567 0,9363 3,1273 -0,0931
1,82 681 0,4237 152-153 0,55588 -1,8012 2,5019 -0,0001
F 1oacc0 0
1,86 667 0,4919 152-154 0,56016 2,5978 2,0116 0,0001
1,86 666 0,4669 170-173 0,58109 3,1927 0,0705 -0,2099
F n2acc0 0
1,87 662 0,4922 170-172 0,58051 0,0858 -3,2730 0,0102
1,76 706 0,4336 215-216 0,58708 -2,8879 1,3141 -0,1205
F Zn ac2e 0
1,80 688 0,5058 215-217 0,59278 -1,3915 -3,0813 0,1734
1,75 707 0,4335 219-220 0,58680 -3,0984 0,6949 -0,1065
F Zn bc2e 0
1,80 689 0,5043 219-221 0,59293 -0,7027 -3,3039 0,1786
1,76 706 0,4398 223-224 0,58722 -3,1485 0,5515 -0,0994
F Zn cd2e 0
2,23 688 0,5045 223-225 0,59281 -0,5385 -3,3335 0,1685
1,76 705 0,4548 236-237 0,58786 -3,1909 0,6083 -0,1000
F Zn c2ef 0
1,80 687 0,5035 236-238 0,59222 -0,6038 -3,3164 0,1639
1,76 705 0,4592 240-241 0,58795 3,2424 -0,3735 -0,0817
F Zn c2eg 0
1,80 687 0,5071 240-242 0,59209 0,3511 3,3645 0,1548
1,76 706 0,4488 250-251 0,58756 3,2275 0,0118 -0,1271
F Zn c2eh 0
1,80 688 0,5045 250-252 0,59286 -0,0604 3,3754 0,1914
Tabela 10: Energia (eV), comprimento de onda ⁄máx (nm), força do oscilador (f ), orbitais
envolvidos e suas contribuições, e dipolo de transição das transições, calculadas via
DFT/OLYP/6-31g(d) em clorofórmio (IEFPCM), responsáveis pelas bandas Q das
ftalocianinas estudadas
104
Dipolo de transição
Molécula Carga Energia ⁄máx (nm) f Orbitais Contribuição
x y z
1,87 662 0,3948 133-135 0,58020 2,9339 0,0041 0,0000
F 4H 0
1,88 661 0,4614 133-134 0,58015 -0,0045 3,1683 0,0000
147-148 -0,4132 2,9794 0,0956
F Al 4H 0 1,81 686 0,3946 0,58101
147-149 -0,0040 0,0962 -2,9828
147-148 3,0974 -0,0910 0,0000
F Zn 4H 0 1,90 652 0,4475 0,58659
147-149 0,0907 3,0974 0,0000
1,96 667 0,4280 145-147 0,57800 2,9619 -0,7772 -0,1283
F 1nh3 1
1,89 658 0,4627 145-146 0,57801 0,7908 3,0645 -0,0452
1,86 665 0,4294 153-155 0,54050 3,0305 -0,3400 0,3206
F 1nh3 4c 1
1,87 663 0,4529 153-154 0,53995 -0,3018 -3,1295 -0,0074
1,87 663 0,4381 157-159 0,50241 -0,0671 3,0918 0,0031
F 2nh3 2
1,87 663 0,4742 157-158 0,50248 3,1973 0,2786 -0,2173
1,85 669 0,4505 157-159 0,56472 -3,0186 -0,8871 0,1724
F 2nh3t 2
1,89 657 0,4639 157-158 0,56301 -0,8851 3,0410 0,0674
1,87 664 0,5087 169-170 0,56699 3,0462 1,3621 -0,0956
F 3nh3 3
1,87 661 0,4318 169-171 0,56513 1,4097 -2,7168 -0,1870
1,87 664 0,4838 181-183 0,66935 -3,1208 -0,8963 -0,1893
F 4nh3 4
1,87 662 0,5518 181-182 0,68191 -1,0374 3,3063 -0,1078
1,86 667 0,4297 153-155 0,57848 2,9500 -0,8382 -0,1619
F nh2ch3 1
1,88 658 0,4598 153-154 0,57831 0,8636 3,0351 -0,0591
1,86 667 0,4353 157-159 0,57691 2,9571 -0,8695 -0,2334
F n3ch3 1
1,89 658 0,4628 157-158 0,57667 0,8884 3,0372 -0,0871
1,82 681 0,3864 152-153 0,55303 -1,6288 2,4530 0,0000
F 1oacc0 0
1,87 665 0,4593 152-154 0,55783 2,5554 1,8766 0,0001
1,87 664 0,4399 170-173 0,58030 3,0956 -0,0024 -0,2020
F n2acc0 0
1,88 661 0,4591 170-172 0,57953 0,0143 -3,1602 0,0155
1,77 702 0,3994 215-216 0,58166 -2,7300 1,3255 -0,1181
F Zn ac2e 0
1,80 686 0,4684 215-217 0,58780 1,3645 -2,9484 0,1638
1,76 704 0,4124 219-220 0,58251 -3,0220 0,6441 -0,1000
F Zn bc2e 0
1,80 687 0,4696 219-221 0,58790 -0,6121 -3,1963 0,1703
1,76 703 0,4197 223-224 0,58329 -3,0743 0,5055 -0,0944
F Zn cd2e 0
1,81 686 0,4703 223-225 0,58812 -0,4641 -3,2219 0,1643
1,77 702 0,4326 236-237 0,58412 -3,1232 0,4803 -0,0938
F Zn c2ef 0
1,81 685 0,4659 236-238 0,58768 -0,4460 -3,2069 0,1611
1,77 701 0,4375 240-241 0,58417 3,1655 -0,2745 -0,0743
F Zn c2eg 0
1,81 685 0,4690 240-242 0,58740 0,2206 3,2408 0,1487
1,77 702 0,4286 250-251 0,58364 3,1458 0,0128 -0,1201
F Zn c2eh 0
1,81 686 0,4682 250-252 0,58788 -0,0910 3,2451 0,1826
Tabela 11: Energia (eV), comprimento de onda ⁄máx (nm), força do oscilador (f ), orbitais
envolvidos e suas contribuições, e dipolo de transição das transições, calculadas
via DFT/OLYP/6-31g(d) em água (IEFPCM), responsáveis pelas bandas Q das
ftalocianinas estudadas
105
Para as porfirinas e clorinas, são apresentados os resultados obtidos para espectro de
absorção em clorofórmio e água nas Tabelas 12 e 13, respectivamente.
Todas as ftalocianinas estudadas tem ⁄máx de absorção dentro da janela terapêutica
proposta para tratamentos de câncer (500-800 nm), apresentando comprimentos de onda
em torno de 660 nm. Apenas o grupo F Zn x2ex, que possui grupo éter e álcool, apresenta
um leve aumento da banda Qx para em torno de 700 nm. Já as porfirinas e clorinas,
possuem todas as suas moléculas dentro da faixa terapêutica quando em água, sendo que o
grupo de porfirinas contendo ligantes piridinicos (grupo P 4bxc) apresentam um aumento
relativamente alto em relação à porfirina livre, e tem seus comprimentos de onda em torno
de 770 nm. No solvente clorofórmio, uma das moléculas, P 4b1cAB, apresenta sua banda
Qx fora da janela terapêutica, em 910 nm. Outros porfirinas contendo ligantes sobem o
⁄máx para em torno de 620 nm. A presença de ı́ons metálicos, para todas as famı́lias, não
mostram alterar o ⁄máx , exceto para Si com fenil como ligante axial, levando a banda Qx
da porfirina para 620 nm.
A variação de ⁄máx praticamente não varia comparando o efeito dos diferentes solventes,
sendo que as ftalocianinas apresentam um aumento em torno de 3 nm da água para o
clorofórmio, sendo que a molécula que mais varia é a F 2nh3t, 10 nm. Já as porfirinas são
mais sensı́veis aos seus ligantes. Em geral apresentam a mesma variação em torno de 3 nm
de aumento da água para o clorofórmio, assim como as ftalocianinas, porém ligantes como
a piridina fazem com que a diferença entre ⁄máx simulado em clorofórmio e água aumente
bastante, como é o caso da P 4b1cB, com um aumento de 175 nm. Ligantes como fenil,
álcool ou ácido carboxı́lico não apresentam variação maior que 3 nm. Já para os ligantes
axiais do Si complexado, apenas o ligante fenil faz com que essa variação aumente 30 nm
da água para o clorofórmio.
A força do oscilador é sempre em torno de 0,4 para as ftalocianinas e relativamente
baixo para as porfirinas e clorinas, exceto quando há piridinas como ligantes (famı́lia
P 4bxc), álcool (Levulan e Foscan) ou ácido carboxı́lico e cetona (Visudine). A complexação
por metais de transição causa uma queda na força do oscilador de todas as moléculas,
independente do ligante axial do Si. Em relação ao solvente, todas as moléculas, pelo
menos em uma das transições, apresentam uma queda na força do oscilador quando em
clorofórmio, exceto F Al 4H, que apresenta um leve aumento em sua única transição e
P Al 4H juntamente com P Zn 4H e P 4fenil que não apresentam variação da força do
oscilador nos diferentes solventes.
Todas as transições com maior contribuição referentes as bandas Q envolvem orbitais
HOMO-1, HOMO, LUMO e LUMO+1. Ora uma das transições é de HOMO para LUMO
e a outra HOMO para LUMO+1 (ftalocianinas e Foscan), ora é de HOMO para LUMO e
a outra HOMO-1 para LUMO (porfirnas e clorinas, exceto Foscan).
Também foi avaliada a variação do comprimento de onda de absorção da banda Qx
versus a energia livre de solvatação na água (Figura 41).
106
Dipolo de transição
Molécula Carga Energia ⁄máx (nm) f Orbitais Contribuição
x y z
2,13 581 0,0031 81-82 0,54044 0,2451 -0,0049 0,0003
P 4H 0
2,27 546 0,0006 81-83 0,49598 -0,0029 -0,1070 0,0000
95-97 0,49130 -0,0027 -0,0581 0,0000
P Al 4H 0 2,26 549 0,0002
95-96 0,49130 0,0580 0,0027 0,0000
95-96 0,49457 0,0156 0,0190 -0,0027
P Zn 4H 0 2,27 546 0,0000
95-97 0,49454 -0,0195 0,0158 0,0000
104-105 0,55646 0,0030 0,1588 -0,0002
P Si 4H 0 2,20 563 0,0014
95-97 0,55632 -0,1576 0,0028 -0,0004
1,99 622 0,0484 161-162 0,56313 0,0005 -0,9954 0,0000
P 4fenil 0
2,12 586 0,0675 161-163 0,53515 1,1407 -0,0004 0,0129
1,68 737 0,2401 165-166 0,48671 2,3841 0,3790 -0,0061
P 4b1cA 1
1,73 718 0,1944 164-166 0,51609 2,0393 -0,6598 0,0183
1,36 910 0,0406 165-166 0,67190 -1,0243 0,4006 -0,0817
P 4b1cB 1
1,61 769 0,0177 164-166 0,68239 -0,6626 -0,0993 -0,0194
1,57 789 0,2653 169-170 0,68008 2,6246 0,0174 0,0077
P 4b2cAC 2
1,63 758 0,0002 168-170 0,60457 0,0254 0,0727 -0,0009
1,60 777 0,3793 169-170 0,59653 3,0929 0,3674 0,0237
P 4b2cAT 2
1,71 724 0,1087 168-170 0,60446 -1,3364 0,8973 0,0140
1,67 744 0,1993 173-174 0,58029 -2,1600 -0,4627 0,0268
P 4b3cA 3
1,71 724 0,0852 173-175 0,54812 0,3410 1,3832 0,0143
1,82 681 0,0352 177-179 0,53021 -0,0075 0,8879 0,0012
P 4b4cA 4
1,85 671 0,1698 177-178 0,59369 1,9362 0,0042 -0,0246
1,87 664 0,0081 177-178 0,61781 0,1232 0,4023 -0,0298
P 4b4cB 4
1,94 639 0,0027 177-179 0,53323 0,1790 -0,1506 -0,0346
2,23 556 0,0000 104-106 0,48970 -0,0066 0,0086 0,0001
P Si 2ch3 0
2,24 554 0,0018 104-105 0,55579 0,1816 0,0198 0,0001
1,99 622 0,0025 136-137 0,63237 0,0000 -0,2274 -0,0051
P Si 2fenil 0
1,99 622 0,0009 136-138 0,64065 -0,0109 -0,0061 0,1348
2,02 613 0,0018 151-152 0,42429 -0,1130 0,1514 -0,0062
Levulan 0
2,15 576 0,0222 150-152 0,38480 0,0729 0,6451 -0,0091
2,00 621 0,0099 185-187 0,49114 0,2292 -0,3860 -0,0340
Visudyne 0
2,11 587 0,0448 185-186 0,51302 -0,9120 0,1732 0,0682
1,99 622 0,0349 177-178 0,57617 0,4889 -0,6891 -0,0009
Foscan 0
2,11 588 0,0400 177-179 0,55248 -0,7369 -0,4881 0,0005
2,17 572 0,0863 82-83 0,59427 -0,0019 1,2776 0,0000
C 4H 0
2,28 544 0,0068 81-83 0,55564 0,3481 0,0009 0,0001
2,21 561 0,1465 96-97 0,63727 -0,0002 1,6445 -0,0003
C Zn 4H 0
2,30 540 0,0050 95-97 0,55729 -0,2983 0,0000 -0,0001
Tabela 12: Energia (eV), comprimento de onda ⁄máx (nm), força do oscilador (f ), orbitais
envolvidos e suas contribuições, e dipolo de transição das transições, calculadas via
DFT/OLYP/6-31g(d) em clorofórmio (IEFPCM), responsáveis pelas bandas Q das
porfirinas estudadas
107
Dipolo de transição
Molécula Carga Energia ⁄máx (nm) f Orbitais Contribuição
x y z
2,15 577 0,0003 81-82 0,53755 0,2368 0,0052 -0,0001
P 4H 0
2,28 543 0,0003 81-83 0,49312 0,0026 0,0781 0,0000
95-96 0,37108 -0,0001 -0,0578 0,0000
P Al 4H 0 2,27 546 0,0002
95-97 0,37107 0,0576 -0,0001 0,0230
95-96 0,48381 0,0085 -0,0038 0,0001
P Zn 4H 0 2,27 546 0,0000
95-97 0,48379 -0,0037 -0,0088 -0,0031
104-105 0,54680 0,0071 0,1297 0,0000
P Si 4H 0 2,22 558 0,0009
104-106 0,54671 0,1289 -0,0070 0,0000
2,00 620 0,0485 161-162 0,56122 -0,0006 0,9948 0,0000
P 4fenil 0
2,12 589 0,0625 161-163 0,53327 -1,0967 0,0004 -0,0017
1,79 694 0,2826 165-166 0,55772 2,5400 0,0828 -0,0063
P 4b1cA 1
1,87 664 0,0780 164-166 0,55980 1,2091 -0,4933 0,0230
1,69 735 0,0533 165-166 0,61238 -1,0406 0,4463 -0,0873
P 4b1cB 1
1,9 653 0,0013 164-166 0,56231 0,1560 0,0485 -0,0152
1,74 712 0,2373 169-170 0,66217 2,3592 0,0156 0,0096
P 4b2cAC 2
1,78 697 0,0131 168-170 0,53036 -0,0030 0,5490 -0,0013
1,71 727 0,3413 169-170 0,59048 2,84 0,2459 0,0196
P 4b2cAT 2
1,83 676 0,0803 168-170 0,58348 -1,1203 0,7305 0,0178
1,76 703 0,1609 173-174 0,53141 1,9272 0,0860 -0,0177
P 4b3cA 3
1,82 680 0,0916 173-175 0,48183 0,8821 1,1280 -0,0151
1,85 670 0,0347 177-179 0,53083 -0,0099 0,8747 0,0012
P 4b4cA 4
1,89 657 0,1500 177-178 0,58301 1,8014 0,0043 -0,0228
1,90 652 0,0085 177-178 0,60483 0,1119 0,4103 -0,0274
P 4b4cB 4
1,98 624 0,0022 177-179 0,51568 0,1581 -0,1386 -0,0346
2,25 551 0,0003 104-106 0,48151 0,0614 0,0290 0,0001
P Si 2ch3 0
2,25 551 0,0013 104-105 0,50956 0,1541 0,0025 0,0001
2,09 592 0,0026 136-137 0,51314 0,0004 0,2255 0,0018
P Si 2fenil 0
2,10 591 0,0009 136-138 0,53140 0,0108 0,0020 -0,1298
2,03 610 0,0015 151-152 0,41362 -0,1340 -0,1144 0,0064
Levulan 0
2,16 573 0,0232 150-152 0,38794 -0,0662 -0,6583 -0,0092
2,00 619 0,0075 185-187 0,49008 -0,2161 -0,3245 0,0312
Visudyne 0
2,12 586 0,0473 185-186 0,53015 -0,9316 0,1960 0,0724
1,99 623 0,0343 177-178 0,58255 0,5410 -0,6412 0,0003
Foscan 0
2,10 592 0,0329 177-179 0,54197 -0,5904 -0,5398 -0,0007
2,18 570 0,0761 82-83 0,59168 -0,0017 1,1950 0,0001
C 4H 0
2,29 542 0,0048 81-83 0,54933 -0,2935 -0,0006 -0,0001
2,22 559 0,1261 96-97 0,63083 -0,0001 1,5230 -0,0002
C Zn 4H 0
2,29 541 0,0050 95-97 0,54863 -0,2985 0,0000 0,0000
Tabela 13: Energia (eV), comprimento de onda ⁄máx (nm), força do oscilador (f ), orbitais
envolvidos e suas contribuições, e dipolo de transição das transições, calculadas
via DFT/OLYP/6-31g(d) em água (IEFPCM), responsáveis pelas bandas Q das
porfirinas estudadas
108
710%
F_Zn_cd2e%
F_Zn_bc2e%
F_Zn_c2eh%
F_Zn_c2eg% F_Zn_c2ef%
690%
F_Al_4H%
F_nh2ch3%
F_1oacc0%
F_Zn_4H%
650%
630%
Foscan% P_4fenil%
Visudyne%
P_4b3cA% P_4b2cAC%
109
610%
P_4b1cB% Levulan%
P_4b2cAT%
P_Si_2fenil%
590% P_4b4cB%
P_4b4cA% P_4H%
570% C_4H%
P_Si_4H% C_Zn_4H%
P_4b1cA% P_Si_2ch3%
550%
P_Al_4H% P_Zn_4H%
530%
H0,700% H0,600% H0,500% H0,400% H0,300% H0,200% H0,100% 0,000%
ΔGsolv"
Figura 41: Gsolv (eV) e comprimento de onda de absorção da banda Qx (nm) de todas as moléculas estudadas.
Da mesma forma que para o volume e o dipolo, comparando o Gsolv com o compri-
mento de onda de absorção da banda Qx , não é possı́vel estabelecer uma tendência. A
adição de ligantes com carga continua acentuada no gráfico e não há um padrão em relação
ao ⁄máx da banda Qx .
110
2
3
4
5
6
7
111
8
9
10
11
12
Carga Metal Areap Areag Vol Dip HOMO LUMO HOMO-LUMO solvat CompOnd
112
Figura 44: Análise HCA para variáveis, considerando nove das doze estudadas, para todo o
conjunto de moléculas.
113
3
Visudyne
F_Zn_c2eh
F_Zn_c2eg
F_Zn_c2ef
2
F_Zn_cd2e
F_Zn_bc2e
F_n2acc_0
1 F_1nh3
F_1oacc_0
Factor2
F_n3ch3 P_4b1cA
F_4H
F_1nh3_4cF_2nh3
114
P_4b2cACP_4b2cAT
F_nh2ch3
0 F_3nh3 F_Al_4H F_Zn_4H
P_4b1cB
P_4fenil
F_2nh3t Foscan
F_4nh3
-1 Levulan
P_4b3cA
P_Si_2fenil C_Zn_4H
P_Zn_4H
P_Al_4H
P_Si_4H
C_4H
P_4b4cA
P_4b4cB P_4H P_Si_2ch3
-2
-4 -2 0 2
Factor1
Figura 45: Análise PCA, gráfico de scores considerando nove das doze variáveis estudadas, para todo o conjunto de moléculas.
Vol
oscila HOMO
0.4
CompOnd
0.2
Dip solvat
Factor2
115
Metal
0.0
-0.2 Carga
LUMO
Figura 46: Análise PCA, gráfico de loadings considerando nove das doze variáveis estudadas, para todo o conjunto de moléculas.
A seguir, retirando essas três variáveis (áreas e diferença entre energia de HOMO e de
LUMO), foram realizadas as análises HCA (Hierarchical Cluster Analysis), apresentadas
nas Figuras 43 e 44, e PCA (Principal Components Analysis), apresentadas nas Figuras
45 e 46. Na HCA, verifica-se que porfirinas e ftalocianinas não são totalmente distinguidas
entre si. Inclusive, os fármacos Levulan e Foscan são agrupados já no inı́cio, embora o
primeiro seja uma porfirina e o segundo, uma clorina, enquanto Visudyne, também porfirina,
encontra-se bastante separado e mais correlacionado a ftalocianinas. A contribuição das
variáveis para a formação dos agrupamentos é dada no gráfico da Figura 44. Volume e
força do oscilador separam a maior parte das ftalocianinas no grupo superior, uma vez
que essas moléculas são as que apresentam maiores valores para essas variáveis, mas é
interessante notar que modificações estruturais como as realizadas as porfirinas P 4b2cA
(CIS e TRANS ) as faz serem agrupadas com as ftalocianinas. A análise PCA confirma essa
distinção observada na análise HCA para os fármacos. Entretanto, as demais amostras
não encontram-se completamente agrupadas. Por esse motivo, foram refeitas as análises
considerando:
• neste mesmo sentido, foi retirada a variável carga, uma vez que o dipolo deve refletir
melhor a distribuição de cargas da molécula como um todo;
Desse modo, as Equações 6.2 e 6.3 as Figuras 47, 48, 49 e 50 apresentam as análises
refeitas considerando apenas cinco variáveis restantes.
116
Figura 47: Análise HCA para amostras de todas as moléculas, considerando cinco variáveis.
117
Figura 48: Análise HCA para variáveis, considerando cinco delas, para todas as moléculas.
118
F_1nh3_4c
F_1nh3
F_nh2ch3
F_n3ch3
P_4b1cA
P_4b2cAC F_3nh3
F_Zn_4H F_Zn_c2eg
C_Zn_4H F_Al_4H
F_1oacc_0 F_Zn_c2ef
F_4H P_4b1cB
F_n2acc_0 F_Zn_bc2e
0 F_Zn_cd2e
Factor2
C_4H
P_Zn_4H P_Al_4H F_Zn_c2eh
P_4H
P_Si_2ch3 Levulan F_2nh3
P_4fenil Foscan
119
Visudyne
F_4nh3
-2
P_4b4cB P_4b4cA
-3 -1 1
Factor1
Figura 49: Análise PCA, gráfico de scores considerando cinco variáveis, para todas as moléculas.
Dip
0.6
solvat
oscila
0.4
Factor2
120
LUMO
0.2
0.0
Vol
-0.2 0.2
Factor1
Figura 50: Análise PCA, gráfico de loadings considerenado cinco variáveis, para todas as moléculas.
Agora, foi criado um grupo com poucos integrantes que consistem das menores porfirinas
e clorinas, incluindo Levulan e Foscan, e um grande grupo com os demais compostos e o
Visudyne. O que também fica evidenciado pela análise de PCA; embora essas análises
sejam bastante preliminares e careçam de correlação com dados de testes biológicos; que
os compostos agrupados com Visudyne sejam promissores para as mesmas aplicações desse
fármaco, assim como os demais compostos com relação à atividade do Levulan e do Foscan.
É importante salientar que, embora um grupo pequeno e bastante similar estrutural-
mente tenha sido agrupado com Levulan e Foscan, foi verificado que algumas ftalocianinas
foram agrupadas bastante próximas ao Visudyne, que é uma porfirina, indicando que
propriedades semelhantes podem ser obtidas para estruturas um pouco diversas.
121
6.3 Dı́meros
Tendo em vista que quando ocorre a dimerização de um composto, este pode passar a
absorver em ⁄máx diferentes do que quando estão em monômeros e este fenômeno pode
alterar a eficiência da PDT [53,87,96], foram avaliadas as contribuições para a conformação
de dı́meros de algumas das moléculas estudadas.
Considerando que as interações existentes na formação de dı́meros são interações
de dispersão, foi com o programa CPMD [237], que permite a inclusão de correções no
funcional BLYP, chamadas DCACP’s, que foram realizados cálculos para este estudo. O
mesmo foi iniciado pela otimização de três dı́meros do tipo empilhado da ftalocianina
livre F 4H: o dı́mero Dim0, em que os monômeros são sobrepostos átomo a átomo, o
dı́mero Dim45, em que um dos monômeros foi rotacionado 45¶ em relação ao Dim0, e o
dı́mero Dim90, em que um dos monômeros foi rotacionado 90¶ também em relação ao
Dim0, fazendo com que os átomos de hidrogênio ligados aos nitrogênios do anel central
não eclipsaram (Figura 51).
As energias de interação (Tabela 14) mostram que os três dı́meros são estáveis em fase
gasosa ( Edim ¡ 0), e o Dim45 é o mais estável entre eles. O grande sistema fi conjudado
deve permitir diversos mı́nimos de energia para os dı́meros, mas o objetivo é avaliar se as
modificações na estrutura do corante desestabilizam a formação do dı́mero. Assim, para os
demais compostos, apenas foram avaliadas as estruturas de partida mais estáveis, ou seja,
em que um dos monômeros está rotacionado 45¶ em relação ao outro monômero, exceto
para três casos. Um deles são os dı́meros da molécula P Al 4H, em que foi necessário o
estudo de estruturas que alteram a posição do ligante axial do Alumı́nio, o Cloro. No
dı́mero DimC, os dı́meros estão contrários de forma que os ligantes axiais voltados para
fora; no DimF, as moléculas estão de frente, com os ligantes axiais voltados para dentro; e
no DimS, os monômeros estão em sequência, com os ligantes axiais voltados para o mesmo
lado. O segundo caso acontece com os dı́meros da molécula P Si 4H, em que um dos
dimeros (Dim45D) tem um dos seus monômeros deslocado no eixo x. E o terceiro caso,
um dos dı́meros da molécula P 4b2cAC está com rotação de 45¶ com os ligantes contendo
nitrogênio sobrepostos e o outro dı́mero está com a mesma rotação porém com os ligantes
contendo nitrogênio não sobrepostos.
Nas Figuras 51 a 58 são apresentadas as estruturas dos dı́meros estudados das moléculas
F 4H, P 4H, P Al 4H, P Zn 4H, P Si 4H, P 4fenil, P 4b2cAC e F 2nh3, respectivamente,
obtidos com a otimização incluindo DCACP’s.
Na Tabela 14 estão a menor distância entre os dois monômeros que compõem o dı́mero
(sempre entre os planos das moléculas), os valores de energia para todas as conformações
de dı́meros, a energia do monômero e a energia de dimerização, ou seja, a diferença dessa
energia para o dobro da energia de um monômero.
Para alguns dı́meros, a energia de dimerização é positiva, indicando que não há formação
do dı́mero em fase gasosa. É importante ressaltar a grande diferença causada pela adição
122
Dim0 Dim45 Dim90
Figura 51: Estruturas dos dı́meros de F 4H utilizados após otimização incluindo DCACP’s.
Figura 53: Estruturas dos dı́meros de P Al 4H utilizados após otimização incluindo DCACP’s.
123
Figura 54: Estrutura do dı́mero de P Zn 4H utilizado após otimização incluindo DCACP’s.
Dim45 Dim45D
Figura 55: Estruturas dos dı́meros de P Si 4H utilizados após otimização incluindo DCACP’s.
Figura 56: Estrutura do dı́mero de P 4fenil utilizado após otimização incluindo DCACP’s.
124
Dim45 Dim45S
Figura 57: Estruturas dos dı́meros de P 4b2cAC utilizados após otimização incluindo DCACP’s.
Figura 58: Estrutura do dı́mero de F 2nh3 utilizado após otimização incluindo DCACP’s.
de ligantes na energia de dimerização nos casos P 4b2cAC e F 2nh3, elevando essa energia,
indicando a não formação de dı́meros. Já para o ligante fenil, em F 4fenil, parece acontecer
o contrário, o ligante facilita a formação do dı́mero em relação a porfirina sem ligantes
P 4H.
Com intuito de avaliar a contribuição do solvente para a estabilização dos dı́meros,
foram realizados cálculos de energia para os dı́meros e monômeros em fase gasosa, água
e clorofórmio, com o modelo de solvatação IEFPCM, como implementado no pacote
computacional Gaussian 09. Para os dı́meros P 4b2cAC e F 2nh3 os cálculos no Gaussian 09
não convergiram. Com esses cálculos foi possı́vel realizar o ciclo termidinâmico apresentado
na Figura 59 e resolver as Equações 6.4 a 6.6.
Para os dı́meros DimC e DimF de P Al 4H, Dim45 de P Zn 4H e os dois dı́meros de
P Si 4H, a inclusão do efeito do solvente desfavoreceu a formação de dı́meros. Paraos
demais dı́meros, essa inclusão indicou facilitação na formação de dı́meros. Porém, nenhum
dos dı́meros apresentou troca de sinal no valor de G0solv . O dı́mero P 4H passou a ter
uma energia de simerização relativamente próxima de zero.
125
Figura 59: Ciclo termodinâmico usado no processo de dimerização.
j
dG = 0 ∆ G0 + Gúf gæsolv (D) ≠ G0solv ≠ 2 Gúf gæsolv (M ) = 0 (6.4)
= Edimero +
Gúsolv ≥ Gúf gæsolv (6.6)
Com a adição do efeito do solvente na energia de dimerização, dos supostos dı́meros que
não indicavam dimerização, apenas o dı́mero 45 de P Si 4H passou a indicar dimerização.
Porém, para quase todos os dı́meros, exceto para DimC e DimF de P Al 4H, para o Dim45D
126
de P Si 4H e para o dı́mero de P 4fenil, a inclusão do efeito do solvente desfavoreceu a
formação de dı́meros. Apenas o dı́mero de P 4fenil, antes com energia de dimerização
positiva, passou a apresentar a energia de dimerização negativa após a inclusão de solvente.
O dı́mero de P 4H passou a ter uma energia de dimerização relativamente próxima de
zero.
Ainda é possı́vel aproveitar os dados para avaliar quais metais complexados ajudam
a formar dı́meros, visto que a complexação contribui para essa formação. No caso, as
energias de Gtotal em solvente mais negativas mostram que a complexação com Alumı́nio é
a menos efetiva para minimizar a formação de dı́meros. Em segundo lugar, o Dim45 com
Silı́cio, seguido do dı́mero com Zinco e por último o Dim45D com Silı́cio.
Além disso, observa-se que a ftalocianina livre tem maior capacidade de formar dı́meros
que a porfirina livre. Provavelmente isso se deve ao sistema de elétrons fi conjugados da
ftalocianina ser maior que o da porfirina, aumentando a interação entre os monômeros.
Outro ponto importante a ser avaliado na formação de dı́meros é se a absorção é afetada
pela agregação, tanto o comprimento de onda de absorção quanto a força do oscilador
dessa transição, já que podem haver interações construtivas ou destrutivas entre os orbitais
moleculares. A Tabela 16 mostra os comprimentos de onda máximo (⁄máx ) de absorção e
força do oscilador das principais transições dos dı́meros estudados e seus monômeros.
127
Fase gasosa Clorofórmio Água
Molécula Conformação Qx Qy Qx Qy Qx Qy
⁄máx (nm) f ⁄máx (nm) f ⁄máx (nm) f ⁄máx (nm) f ⁄máx (nm) f ⁄máx (nm) f
Monômero 644 0,2696 633 0,3341 655 0,4274 662 0,4979 662 0,3948 661 0,4614
Dim0 618 0,3519 604 0,4550 634 0,5834 626 0,7215 633 0,5336 627 0,6691
F 4H
Dim45 631 0,3448 629 0,3868 633 0,5800 630 0,6645 617 0,5309 611 0,6474
Dim90 613 0,3770 612 0,3752 631 0,6347 630 0,6309 632 0,5851 631 0,5808
Monômero 585 0,0010 548 0,0004 581 0,0031 546 0,0006 577 0,0003 543 0,0003
P 4H
Dim45 711 0,0000 692 0,0000 708 0,0000 689 0,0000 704 0,0000 685 0,0000
Monômero 554 0,0001 x x 549 0,0002 x x 546 0,0002 x x
DimC 620 0,0000 x x 609 0,0000 x x 604 0,0001 x x
P Al 4H
DimF 656 0,0011 x x 632 0,0004 x x 623 0,0003 x x
DimS - - x x 1077 0,0001 x x 804 0,0002 x x
Monômero 546 0,0002 x x 546 0,0000 x x 546 0,0000 x x
P Zn 4H
Dim45 599 0,0001 x x 598 0,0002 x x 598 0,0001 x x
Monômero 573 0,0012 x x 563 0,0014 x x 558 0,0009 x x
P Si 4H Dim45 643 0,0000 x x 633 0,0000 x x 643 0,0000 x x
Dim45D 650 0,0001 x x 641 0,0001 x x 650 x x
Monômero 622 0,0244 x x 622 0,0484 x x 620 0,0485 x x
P 4fenil
Dim45 708 0,0001 x x 708 0,0001 x x 706 0,0001 x x
128
Tabela 16: Comprimentos de onda máximo (⁄máx ) de absorção (nm) e força do oscilador das principais transições dos dı́meros estudados e seus
monômeros
É possı́vel observar que todos os dı́meros estão dentro da janela terapêutica da PDT,
apesar de em alguns casos como os dı́meros de P 4H e P Al 4H terem apresentado um
deslocamento de mais de 100 nm. O comprimento de onda de absorção, no caso da
ftalocianina, apresenta um leve deslocamento para ⁄máx menores, maior para os dı́meros
com menor interação. Já para as porfirinas, sempre o dı́mero apresenta comprimento de
onda maior. Esse deslocamento é menor quando há metal complexado, principalmente
para o Zinco. A substituição do ligante fenil, apesar de apresentar um deslocamento
relativamente alto de 86 nm do monômero para o dı́mero, ainda assim é menor do que o
deslocamento apresentado pela a porfirina livre (126 nm). Os dı́meros que apresentam
menor variação com relação ao comprimento de onda do monômero, são os dı́meros da
ftalocianina.
Observa-se uma interação construtiva, aumentando a força do oscilador das transições,
apenas para a ftalocianina, não chega a ser o dobro, porém sempre mais alto. Já para as
porfirinas, muitas das forças do oscilador se anulam completamente, indicando que quando
dimerizados, esses compostos perdem sua eficiência fotodinâmica.
Outra forma de avaliação encontrada foi comparar o quanto o dı́mero se solvata e o
comprimento de onda de absorção. O gráfico com essas informações estão na Figura 60.
129
0,0400%
P_4fenil%&%Dim45%
0,0200%
P_Zn_4H%
P_Al_4H%
P_Si_4H%
P_4fenil%
P_Si_4H%&%Dim45D%
P_4H%
F_4H%
0,0000% P_4H%&%Dim45%
500% 520% 540% 560% 580% 600% 620% 640% 660% 680% 700% 720% 740%
&0,0200%
&0,0400%
ΔΔGsolv'
Água%
Clorofórmio%
&0,0600%
130
Monômero%
P_Zn_4H%&%Dim45%
&0,0800%
P_Si_4H%&%Dim45%
&0,1000%
F_4H%&%Dim90% F_4H%&%Dim0%
&0,1400%
Comprimento'de'onda'(nm)'
Figura 60: Gráfico do Gsolvatação de cada dı́mero pelo seu comprimento de absorção máximo em água e clorofórmio, além dos ⁄máx de absorção
de cada monômero.
Não há mudança significativa no Gsolvatação em relação aos dois solventes, conse-
quentemente, no comprimento de onda também não.
Para os dois cálculos de dinâmica molecular realizados, foram extraı́dos dados dos
200 ns de simulação obtidos, para certificar-se da espessura da membrana, verificando a
distância entre o centro de massa dos grupos fosfato das duas camadas (Figura 61).
Outra propriedade a ser analisada em uma dinâmica molecular de membranas é a
área ocupada por um lipı́deo. Esses dados também foram extraı́dos das duas dinâmicas,
mostrados na Figura 62.
Observando estes dois gráficos, é possı́vel dizer que a dinâmica entrou em equilı́brio
a partir do 100¶ ns. Sendo assim, foram medidos os dois parâmetros anteriores, fazendo
uma média a partir do 100 ns até o fim da simulação. A espessura da membrana no
equilı́brio foi 42,18 Å para a Dinâmica F e 42,75 Å para a Dinâmica P. A área por lipı́deo
51,65 Å2 da Dinâmica F e 50,93 Å2 da a Dinâmica P, condizem razoavelmente bem com
medições experimentais e cálculos computacionais encontrados na literatura, desta mesma
membrana (POPG) ou dessa mixada com outra (Tabela 17).
131
Figura 61: Gráfico da distância (Å) entre o centro de massa dos grupos fosfato das duas
camadas lipı́dicas da membrana pelo tempo de simulação (ns), onde a curva preta
representa Dinâmica F e a curva vermelha, Dinâmica P.
Figura 62: Gráfico da área (Å2 ) ocupada por um lipı́deo, pelo tempo de simulação (ns), onde a
curva preta representa Dinâmica F e a curva vermelha, Dinâmica P.
132
Referência Membrana Modo de obtenção Espessura (Å) Área por lipideo (Å2 )
Cowley, 1978 [238] POPG Experimental 39,0 66,0
DM [1] 39,7 ± 0,6 63,6 ± 1,2
Janosi, 2010 POPC/POPG
DM [2] 39,1 ± 0,6 63,4 ± 0,8
Róg, 2003 [227] POPE/POPG DM 38,7 63,9
Elmore, 2006 [68] POPG DM 42,7 ± 0,6 56,1
Zhao, 2007 [61] POPG DM 43,9 ± 0,2 53,0 ± 0,6
Polyansky, 2007 [239] POPG DM 39,4 56,5
Dickey, 2008 2008 [63] POPG DM 44,5 ± 0,5 54,8 ± 0,4
Tolokh, 2009 [71] POPG DM 36,3 66,3 ± 1,4
Experimental 36,3 66,3
Kucerka, 2012 [69] POPG Experimental 36,5 66,0
DM 36,7 66,0
DM [3] 37,8 57,1
Mao, 2013 [70] POPG
DM [4] 43,2 53,4
133
6.4.1.2 Equilı́brio da molécula
Figura 63: Orientação da molécula F 2nh3 (Dinâmica F), onde as esferas laranjas são os grupos
fosfato e as esferas azuis representam as caudas dos lipideos.
134
Figura 64: Orientação da molécula P 4b2cAT (Dinâmica P), onde as esferas laranjas são os
grupos fosfato e as esferas azuis representam as caudas dos lipideos.
Figura 65: Gráfico da orientação das moléculas F 2nh3 (linha preta) e P 4b2cAT (linha verme-
lha) durante o tempo da dinâmica molecular, sendo o ângulo da orientação feito
entre a normal do plano da molécula e o eixo z.
135
Figura 66: Gráfico da profundidade das moléculas F 2nh3 (linha preta) e P 4b2cAT (linha
vermelha) durante o tempo de simulação, obtido a partir da densidade de massa da
molécula em relação ao eixo z.
tem suas cargas nas suas extremidades, já a porfirina tem suas cargas relativamente mais
ao centro da molécula.
Segundo o gráfico de distribuição de oxigênio e suas espécies reativas apresentado por
Cordeiro e colaboradores [51], a posição que as moléculas ficaram nas duas dinâmicas é
relativamente favorável para a interação do fotossensibilizador com o oxigênio molecular,
sendo a posição da F 2nh3 mais vantajosa que a P 4b2cAT.
Após terminadas as dinâmicas, foram calculados os estados excitados de cada molécula
a partir de um cálculo QM/MMPol.
136
50
40
30
Frequência
20
10
0
1 1,5 2 2,5 3
Energia (eV)
40
30
Frequência
20
10
0
1 1,5 2 2,5 3
Energia (eV)
137
Figura 69: Representação esquemática de uma monocamada construı́da pelo modelo contı́nuo
eletrostático.
138
1,9%
Qx% Qy%
1,89%
1,88%
1,87%
1,86%
1,85%
)20% )15% )10% )5% 0% 5% 10% 15% 20%
Tabela 18: Energia de absorção (eV), comprimento de onda (nm) e dipolo das bandas Qx e Qy ,
bem como a diferença de energia entre elas, geradas pelos métodos QM/MMPol,
com e sem a presença da membrana, e contı́nuo eletrostático
na Tabela 18, estão os valores de absorção dos métodos QM/MMPol e método contı́nuo
eletrostático no ponto de equilı́brio.
Foram utilizados dados experimentais da ftalocianina livre para comparação [240]:
⁄Qx =686nm; ⁄Qy =622nm; ⁄Qx -Qy =64nm. Para a banda Qx , comparando o método
contı́nuo eletrostático com o experimental, a diferença é de 38 nm, enquanto QM/MMPol
tem apenas 2 nm de diferença. Para a banda Qy , essa diferença aumenta nos dois casos:
o contı́nuo eletrostático apresenta 43 nm e o QM/MMPol com membrana apresenta 34
nm de diferença em relação ao dado experimental. Para os dois métodos, a tendência é
subestimar o primeiro estado excitado e superestimar o segundo. Os dois métodos têm a
mesma tendência sobre o dipolo, exceto QM/MMPol na ausência da membrana.
Comparando a diferença entre as bandas Qx e Qy , o método QM/MMPol prevê de forma
mais eficaz essa diferença em relação ao método contı́nuo eletrostático, que praticamente
torna degeneradas as duas transições. Isso acontece pois o método QM/MMPol inclui
as distorções (movimentos similares à batida de asa de uma borboleta) estruturais da
ftalocianina causadas por efeitos térmicos, que influenciam suas propriedades eletrônicas.
Além disso, esses efeitos são também responsáveis pelo alargamento da banda Qx com a
variação da temperatura [79, 240].
Como os resultados obtidos pelo método QM/MMPol, incluindo a membrana e sem a
139
50
40
30
Frequência
20
10
0
1 1,5 2 2,5 3
Energia (eV)
membrana são muito similares, foi realizada uma segunda dinâmica molecular da molécula
apenas em água. Dessa forma, é possı́vel verificar se o afastamento das bandas Qx e
Qy é causada pela presença da membrana durante a dinâmica molecular ou se apenas a
presença de um solvente pode levar a resultados similares. A dinâmica feita apenas em
água foi menos custosa que a dinâmica incluindo a membrana, assim seria possı́vel obter os
resultados desejados mais rapidamente. Foi feito um histograma com os ⁄máx de absorção
obtidos com os últimos 100 ns da dinâmica em água para melhor análise, apresentados na
Figura 71.
Visto que as energias de absorção obtidos a partir da dinâmica em água, para as
bandas Qx e Qy são 1,791 e 1,870, e os ⁄máx 694 e 664, respectivamente, é possı́vel afirmar
que o afastamento das bandas Qx e Qy não são causadas pela presença da membrana.
Portanto, é possı́vel obter resultados para essas bandas realizando uma dinâmica molecular
da molécula apenas em solvente, diminuindo o custo do método.
6.5 Fosforescência
Outra propriedade importante para se conhecer sobre um fotossensibilizador é o seu
comprimento de onda de emissão. No caso das moléculas estudadas neste trabalho, a
emissão ocorre através da fosforescência, ou seja, o elétron no estado excitado singleto,
realiza uma conversão intersistemas (processo CIS da Figura 1) trocando de spin e
140
passando para o estado excitado tripleto. Do estado excitado tripleto, ele decai até o
estado fundamental, caracterizando a fosforescência (processo P da Figura 1), emitindo
uma quantidade energia especı́fica para cada molécula. Durante a PDT, essa energia pode
ser transferida para outras moléculas, como oxigênio triplete.
Sabe-se que a emissão dessas moléculas é bastante influenciada pelo meio em que
estão, o que dificulta a comparação com os dados obtidos computacionalmente [241, 242].
Segundo Melo et al., a molécula do Levulan, em meio hidrofóbico, apresenta o méximo de
emissão em 634 nm [241]. Já o Visudyne, em plasma, apresenta um máximo de emissão
em 692 nm [243], e o Foscan, em meio celular, em 652 nm [242]. A ftalocianina livre,
em tolueno, cuja constante dielétrica é 2,38, apresentou um pico de emissão em 698 nm
segundo Van Mingroot et al. [244]. Para Das [245], em solução aquosa, a porfirina livre
apresentou emissão em 610 nm e complexada com Zn2+ , em 575 nm.
Utilizando a mesma metodologia da membrana descrita por IEFPCM da secção anterior,
foram feitos cálculos de energia do estado tripleto excitado e singleto, mantendo a geometria
do estado tripleto, para todos os pontos ao longo do gradiente dielétrico. A diferença de
energia entre esses dois estados é o que se espera ser a energia de emissão da molécula
( Eemissão ). Nas Tabelas 30, 32, 31 e 33 que se encontram no Apêndice A, estão o
Eemissão para todos os pontos ao longo da membrana descrita por IEFPCM. Nesse
capı́tulo, apenas são apresentados os valores de Eemissão para o ponto de equilı́brio da
molécula na membrana (ponto 0), na Tabela 19.
Comparando os valores experimentais citadas acima com os dados obtidos com as
simulações ( Eemissão ), colunas 1 e 2 da Tabela 19), os seguintes erros são apresentados:
para Foscan, 31%; para o Levulan, 28%; para o Visudyne, 25%; para a F 4H, 43%; para a
P 4H, 30% e para P Zn 4H, 18%. A média desses erros é 29%, portanto foram incluı́das
duas colunas na Tabela 19 retirando 29% dos valores obtidos com as simulações. Além
disso, também são apresentados os valores de variação mádima entre Eemissão ao longo
de toda a membrana.
Analisando a Tabela 19, observa-se que o método mantém a mesma tendência, superes-
timando os valores de Eemissão . Porém, levando em conta o erro relativamente alto em
relação aos dados experimentais e a variação desse erro, talvez seja necessário um estudo
mais aprofundado dos estados excitados dessas moléculas para investigar outros tripletos
que não o primeiro, para conferir os que participam efetivamente da fosforescência.
De toda forma, ainda que esse estudo necessite ser mais aprofundado, esta forma de
simulação do ambiente da membrana se mostrou muito poderosa, pois é possı́vel verificar
os corantes cujas espécies são mais sensı́veis ao meio a um relativamente baixo custo
computacional.
141
Eemissão
Variação Eemissão
no ponto de equilı́brio
Molécula ao longo da membrana de PCM
Simulado real Simulado Corrigido
(eV) (nm) (eV) (nm) (eV) (nm)
F 4H 1,020 1216 1,436 863 0,057 69
F 1nh3 1,016 1220 1,431 866 0,037 43
F 1nh3 4c 1,008 1230 1,428 874 0,033 39
F 2nh3 1,014 1223 1,428 868 0,035 41
F 2nh3t 1,027 1207 1,447 857 0,071 83
F 3nh3 1,019 1217 1,435 864 0,049 60
F 4nh3 1,214 1022 1,709 725 0,067 56
P 4H 1,432 866 2,016 615 0,053 32
P Al 4H 1,683 737 2,371 523 0,199 94
P Zn 4H 1,765 703 2,485 499 0,043 17
P Si 4H 1,621 765 2,283 543 0,100 47
P 4b2cAC 1,207 1027 1,700 729 0,384 400
P 4b2cAT 1,137 1090 1,602 774 0,410 380
P Si 2ch3 1,742 712 2,453 505 0,047 19
P Si 2fenil 1,822 680 2,566 483 0,018 7
Levulan 1,416 876 1,994 622 0,083 51
Visudyne 1,349 919 1,900 653 0,124 85
Foscan 1,311 946 1,847 671 0,130 91
C 4H 1,445 858 2,035 609 0,050 30
Tabela 19: Energia de absorção (eV), comprimento de onda (nm) e dipolo das bandas Qx e Qy ,
bem como a diferença de energia entre elas, geradas pelos métodos QM/MMPol,
com e sem a presença da membrana, e contı́nuo eletrostático
142
7 Conclusões
144
o custo computacional do método QM/MM é alto em relação sua acuracidade.
145
Referências
9 LIAO, M.-S.; WATTS, J. D.; HUANG, M.-J. DFT study of unligated and ligated
manganese(II) porphyrins and phthalocyanines. Inorganic chemistry, v. 44, n. 6, p.
1941–9, mar. 2005. ISSN 0020-1669.
147
24 REN, X.-f. et al. A density functional theory study on photophysical properties of red
light-emitting materials: Meso-substituted porphyrins. Journal of Photochemistry
and Photobiology A: Chemistry, v. 203, n. 2-3, p. 92–99, abr. 2009. ISSN 10106030.
25 REN, X.-F. et al. Density functional theory study on photophysical properties of
the porphyrins derivatives with through-bond energy transfer characters. Organic
Electronics, v. 11, n. 6, p. 979–989, jun. 2010. ISSN 15661199.
26 ROBERTSON, C. A.; EVANS, D. H.; ABRAHAMSE, H. Photodynamic therapy
(PDT): a short review on cellular mechanisms and cancer research applications for
PDT. Journal of photochemistry and photobiology. B, Biology, v. 96, n. 1, p. 1–8, jul.
2009. ISSN 1873-2682.
27 PALMA, M.; CáRDENAS-JIRóN, G. I.; Menéndez Rodrı́guez, M. I. Effect of chlorin
structure on theoretical electronic absorption spectra and on the energy released by
porphyrin-based photosensitizers. The journal of physical chemistry. A, v. 112, n. 51,
p. 13574–83, dez. 2008. ISSN 1520-5215.
28 DOUGHERTY, T. J. Yearly Review - Photodynamic Therapy. Photochemistry and
Photobiology, v. 58, n. 6, p. 895–900, 1993.
29 LU, T.-T. et al. Density functional theory studies of N-protonation of the free base
phthalocyanine. Journal of Molecular Structure: THEOCHEM, v. 860, n. 1-3, p.
141–149, jul. 2008. ISSN 01661280.
30 BATHELT, C. M. et al. Electronic structure of compound I in human isoforms
of cytochrome P450 from QM/MM modeling. Journal of the American Chemical
Society, v. 127, n. 37, p. 12900–8, out. 2005. ISSN 0002-7863.
31 INSTITUTE, N. C. Cancer Trends Progress Report. [S.l.]: www.cancer.gov.
32 DOLMANS, D. E. J. G. J.; FUKUMURA, D.; JAIN, R. K. Photodynamic therapy
for cancer. Nature reviews. Cancer, v. 3, n. 5, p. 380–7, maio 2003. ISSN 1474-175X.
33 RAVANAT, J. L. et al. Singlet oxygen induces oxidation of cellular DNA. The Journal
of biological chemistry, v. 276, n. 8, p. 40601–4, fev. 2001. ISSN 0021-9258.
34 SHARMAN, W. M.; ALLEN, C. M.; LIER, J. E. van. Photodynamic therapeutics:
basic principles and clinical applications. Drug Discovery Today, v. 4, n. 11, p.
507–517, nov. 1999. ISSN 13596446.
35 STARK, G. Functional consequences of oxidative membrane damage. The Journal
of membrane biology, v. 205, n. 1, p. 1–16, maio 2005. ISSN 0022-2631.
36 DINI, L. et al. Rose bengal acetate photodynamic therapy-induced autophagy. Cancer
biology & therapy, v. 10, n. 10, p. 1048–55, nov. 2010. ISSN 1555-8576.
37 KONAN, Y. N.; GURNY, R.; ALLéMANN, E. State of the art in the delivery of pho-
tosensitizers for photodynamic therapy. Journal of Photochemistry and Photobiology
B: Biology, v. 66, n. 2, p. 89–106, mar. 2002. ISSN 10111344.
38 BUYTAERT, E.; DEWAELE, M.; AGOSTINIS, P. Molecular effectors of multiple
cell death pathways initiated by photodynamic therapy. Biochimica et biophysica
acta, v. 1776, n. 1, p. 86–107, set. 2007. ISSN 0006-3002.
148
39 KUBLER, A. C. Photodynamic therapy. Medical Laser Application, v. 20, n. 1, p.
37–45, maio 2005. ISSN 16151615.
41 BROWN, S. B.; BROWN, E. A.; WALKER, I. The present and future role of
photodynamic therapy in cancer treatment. The lancet oncology, v. 5, n. 8, p. 497–
508, ago. 2004. ISSN 1470-2045.
149
53 ALLISON, R. R. et al. Photosensitizers in clinical PDT. Photodiagnosis and Pho-
todynamic Therapy, v. 1, n. 1, p. 27–42, maio 2004. ISSN 15721000.
60 ENGELMAN, D. M. Membranes are more mosaic than fluid. Nature, v. 438, n. 7068,
p. 578–80, dez. 2005. ISSN 1476-4687.
62 MILLS, F. D. et al. The helical structure of surfactant peptide KL4 when bound to
POPC: POPG lipid vesicles. Biochemistry, v. 47, n. 32, p. 8292–8300, 2008. ISSN
00062960.
63 DICKEY, A.; FALLER, R. Examining the Contributions of Lipid Shape and He-
adgroup Charge on Bilayer Behavior. Biophysical Journal, v. 95, p. 2636–2646,
2008.
150
66 CHEN, W. W. et al. Dimers of human —-defensins and their interactions with the
POPG membrane. Molecular Simulation, v. 39, n. 11, p. 849–859, set. 2013. ISSN
0892-7022.
67 MANNA, M.; MUKHOPADHYAY, C. Molecular dynamics simulations of the in-
teractions of kinin peptides with an anionic POPG bilayer. Langmuir : the ACS
journal of surfaces and colloids, v. 27, n. 7, p. 3713–22, abr. 2011. ISSN 1520-5827.
68 ELMORE, D. E. Molecular dynamics simulation of a phosphatidylglycerol membrane.
FEBS letters, v. 580, n. 1, p. 144–8, jan. 2006. ISSN 0014-5793.
69 KUCERKA, N. et al. Scattering density profile model of POPG bilayers as determined
by molecular dynamics simulations and small-angle neutron and X-ray scattering
experiments. The journal of physical chemistry. B, v. 116, n. 1, p. 232–9, jan. 2012.
ISSN 1520-5207.
70 MAO, Y. et al. Binding competition to the POPG lipid bilayer of Ca2+, Mg2+,
Na+, and K+ in different ion mixtures and biological implication. The journal of
physical chemistry. B, v. 117, n. 3, p. 850–8, jan. 2013. ISSN 1520-5207.
71 TOLOKH, I. et al. Binding free energy and counterion release for adsorption of the
antimicrobial peptide lactoferricin B on a POPG membrane. Physical Review E,
v. 80, n. 3, p. 031911, set. 2009. ISSN 1539-3755.
72 CORDEIRO, R. M.; MIOTTO, R.; BAPTISTA, M. S. Photodynamic efficiency
of cationic meso-porphyrins at lipid bilayers: insights from molecular dynamics
simulations. The journal of physical chemistry. B, v. 116, n. 50, p. 14618–27, dez.
2012. ISSN 1520-5207.
73 LIU, Z. et al. Theoretical investigation of the molecular, electronic structures and
vibrational spectra of a series of first transition metal phthalocyanines. Spectrochimica
acta. Part A, Molecular and biomolecular spectroscopy, v. 67, n. 5, p. 1232–46, ago.
2007. ISSN 1386-1425.
74 BEREZIN, K. V.; NECHAEV, V. V. Calculation of Chlorin IR Spectrum by Density
Functional Theory. Chemistry of Natural Compounds, v. 39, n. 6, p. 540–548, nov.
2003. ISSN 0009-3130.
75 SESALAN, B. ; KOCA, A.; GUL, A. Water soluble novel phthalocyanines containing
dodeca-amino groups. Dyes and Pigments, v. 79, n. 3, p. 259–264, dez. 2008. ISSN
01437208.
76 GOKSEL, M.; DURMUS, M.; ATILLA, D. Amino-functionalized water-soluble
zinc phthalocyanines: Synthesis, photophysical, photochemical and protein binding
properties. Journal of Photochemistry and Photobiology A: Chemistry, v. 266, p.
37–46, ago. 2013. ISSN 10106030.
77 MASHAZI, P. et al. Applications of polymerized metal tetra-amino phthalocyanines
towards hydrogen peroxide detection. Journal of Porphyrins and Phthalocyanines,
v. 14, n. 03, p. 252–263, mar. 2010. ISSN 1088-4246.
78 KOBAYASHI, M. et al. Solvent effect on the aggregation of amphiphilic phthalocya-
nines substituted by polyethyleneoxide. Journal of Porphyrins and Phthalocyanines,
v. 16, n. 2, p. 183–191, fev. 2012. ISSN 1088-4246.
151
79 CABRAL, B. J. et al. Free base phthalocyanine: Influence of thermal effects and
dimerization on the electronic absorption spectrum. Chemical Physics Letters, v.
595-596, p. 97–102, mar. 2014. ISSN 00092614.
80 MOISENOVICH, M. M. et al. Novel Photosensitizers Trigger Rapid Death of Ma-
lignant Human Cells and Rodent Tumor Transplants via Lipid Photodamage and
Membrane Permeabilization. Plos One, v. 5, n. 9, 2010.
81 GURINOVICH, G. P.; SEVCHENKO, A. N.; Solov’ev K. N. The Spectroscopy of
the Porphyrins. Soviet Physics Uspekhi, v. 6, n. 1, p. 173–234, 1963.
82 GOUTERMAN, M. Study of the Effects of Substitution on the Absorption Spectra of
Porphin. The Journal of Chemical Physics, v. 30, n. 5, p. 1139, 1959. ISSN 00219606.
83 ALLISON, R. R.; SIBATA, C. H. Oncologic photodynamic therapy photosensitizers:
a clinical review. Photodiagnosis and photodynamic therapy, v. 7, n. 2, p. 61–75, jun.
2010. ISSN 1873-1597.
84 BONNETT, R. Photosensitizers of the porphyrin and phthalocyanine series for
photodynamic therapy. Chemical Society Reviews, v. 24, n. 1, p. 19–33, 1995. ISSN
0306-0012.
85 HASHIMOTO, T. et al. Theoretical Study of the Q and B Bands of Free-Base,
Magnesium, and Zinc Porphyrins, and Their Derivatives. The Journal of Physical
Chemistry A, v. 103, n. 12, p. 1894–1904, mar. 1999. ISSN 1089-5639.
86 CHAUDHURI, T. et al. Environmental effects on luminescence of meso-tetra-2-
chlorophenylporphyrin. Spectrochimica acta. Part A, Molecular and biomolecular
spectroscopy, v. 76, n. 2, p. 230–8, jul. 2010. ISSN 1873-3557.
87 DEROSA, F. S.; BENTLEY, V. L. B. Photodynamic therapy of skin cancers:
Sensitizers, clinical studies and future directives. Pharmaceutical Research, v. 17, p.
1447–1455, 2000.
88 ALMLOF, J. Ab initio calculations on porphin. International Journal of Quantum
Chemistry, v. 8, n. 6, p. 915–924, nov. 1974. ISSN 0020-7608.
89 De Angelis, F. et al. Time-dependent and coupled-perturbed DFT and HF investi-
gations on the absorption spectrum and non-linear optical properties of push–pull
M(II)–porphyrin complexes (M=Zn, Cu, Ni). Chemical Physics Letters, v. 447, n. 1-3,
p. 10–15, out. 2007. ISSN 00092614.
90 GHOSH, A. et al. Effect of five membered versus six membered meso-substituents on
structure and electronic properties of Mg(II) porphyrins: a combined experimental
and theoretical study. Inorganic chemistry, v. 49, n. 18, p. 8287–97, set. 2010. ISSN
1520-510X.
91 MINAEV, B.; GREN, H. A. Theoretical DFT study of phosphorescence from porphy-
rins. Chemical Physics, v. 315, n. 3, p. 215–239, ago. 2005. ISSN 03010104.
92 WANG, L.-F.; MENG, X.-W.; TANG, F.-Q. Density functional theory study of
electronic absorption spectra and intermolecular interactions of porphyrin–borate
complexes. Journal of Molecular Structure: THEOCHEM, v. 956, n. 1-3, p. 26–32,
set. 2010. ISSN 01661280.
152
93 ZHANG, X. X. et al. IR and Raman vibrational assignments for metal-free phtha-
locyanine from density functional B3LYP/6-31G(d) method. Chinese Journal of
Chemistry, v. 22, n. 4, p. 325–332, 2004.
97 ISHII, K.; IWASAKI, M.; KOBAYASHI, N. The Q absorption band of the excited
triplet phthalocyanine dimer. Chemical Physics Letters, v. 436, p. 94–98, 2007. ISSN
00092614.
101 BALU, R.; BYRD, E. F. C.; RICE, B. M. Assessment of dispersion corrected atom
centered pseudopotentials: Application to energetic molecular crystals. Journal of
Physical Chemistry B, v. 115, p. 803–810, 2011. ISSN 15206106.
104 PLIEGO Jr, J. R. Modelos continuos do solvente: fundamentos. Quı́mica Nova, v. 29,
n. 3, p. 535–542, jun. 2006. ISSN 0100-4042.
153
107 MELLO, P. Homem-de et al. The effects of solvation in the theoretical spectra of
cationic dyes. Theoretical Chemistry Accounts, v. 113, n. 5, p. 274–280, maio 2005.
ISSN 1432-881X.
108 ALMERINDO, G. I.; TONDO, D. W.; PLIEGO, J. R. Ionization of Organic Acids
in Dimethyl Sulfoxide Solution: A Theoretical Ab Initio Calculation of the p K a
Using a New Parametrization of the Polarizable Continuum Model. The Journal of
Physical Chemistry A, v. 108, n. 1, p. 166–171, jan. 2004. ISSN 1089-5639.
109 MIERTUS, S.; SCROCCO, E.; TOMASI, J. Electrostatic interaction of a solute with
a continuum. A direct utilizaion of AB initio molecular potentials for the prevision of
solvent effects. Chemical Physics, v. 55, n. 1, p. 117–129, fev. 1981. ISSN 03010104.
110 STRATMANN, R. E.; SCUSERIA, G. E.; FRISCH, M. J. An efficient implementation
of time-dependent density-functional theory for the calculation of excitation energies
of large molecules. The Journal of Chemical Physics, v. 109, n. 19, p. 8218, 1998.
ISSN 00219606.
111 ENGELMANN, F. M. et al. Interaction of cationic meso-porphyrins with liposomes,
mitochondria and erythrocytes. Journal of bioenergetics and biomembranes, v. 39,
n. 2, p. 175–85, abr. 2007. ISSN 0145-479X.
112 AL-ABDUL-WAHID, M. S. et al. A combined NMR and molecular dynamics study
of the transmembrane solubility and diffusion rate profile of dioxygen in lipid bilayers.
Biochemistry, v. 45, n. 35, p. 10719–28, out. 2006. ISSN 0006-2960.
113 SUBCZYNSKI, W. K.; HYDE, J. S. Concentration of oxygen in lipid bilayers using
a spin-label method. Biophysical journal, v. 41, n. 3, p. 283–6, mar. 1983. ISSN
0006-3495.
114 ETHIRAJAN, M. et al. The role of porphyrin chemistry in tumor imaging and
photodynamic therapy. Chemical Society reviews, v. 40, n. 1, p. 340–62, jan. 2011.
ISSN 1460-4744.
115 KEPCZYNSKI, M. et al. Do liposome-binding constants of porphyrins correlate with
their measured and predicted partitioning between octanol and water? Photochemis-
try and photobiology, v. 76, n. 2, p. 127–34, ago. 2002. ISSN 0031-8655.
116 BRONSHTEIN, I. et al. Porphyrin depth in lipid bilayers as determined by iodide and
parallax fluorescence quenching methods and its effect on photosensitizing efficiency.
Biophysical journal, v. 87, n. 2, p. 1155–64, ago. 2004. ISSN 0006-3495.
117 SHOLTO, A.; EHRENBERG, B. Hydrophobicity, topography in membranes and
photosensitization of silicon phthalocyanines with axial ligands of varying lengths.
Photochemical & photobiological sciences : Official journal of the European Photoche-
mistry Association and the European Society for Photobiology, v. 7, n. 3, p. 344–51,
mar. 2008. ISSN 1474-905X.
118 CAMINOS, D. A. et al. Mechanisms of Escherichia coli photodynamic inac-
tivation by an amphiphilic tricationic porphyrin and 5,10,15,20-tetra(4-N,N,N-
trimethylammoniumphenyl) porphyrin. Photochemical & photobiological sciences
: Official journal of the European Photochemistry Association and the European
Society for Photobiology, v. 7, n. 9, p. 1071–8, set. 2008. ISSN 1474-905X.
154
119 ENGELMANN, F. M. et al. Determination of n-octanol/water partition and mem-
brane binding of cationic porphyrins. International journal of pharmaceutics, v. 329,
n. 1-2, p. 12–8, fev. 2007. ISSN 0378-5173.
121 GINER CASARES, J. J. et al. Methylene blue adsorption on a DMPA lipid langmuir
monolayer. Chemphyschem : a European journal of chemical physics and physical
chemistry, v. 11, n. 10, p. 2241–7, jul. 2010. ISSN 1439-7641.
124 ERDTMAN, E. et al. Modelling the behavior of 5-aminolevulinic acid and its alkyl
esters in a lipid bilayer. Chemical Physics Letters, v. 463, n. 1-3, p. 178–182, set.
2008. ISSN 00092614.
125 ERDTMAN, E. et al. Modelling the behavior of 5-aminolevulinic acid and its alkyl
esters in a lipid bilayer. Chemical Physics Letters, v. 463, n. 1-3, p. 178–182, set.
2008. ISSN 00092614.
126 GRATER, F. et al. Protein/Ligand Binding Free Energies Calculated With Quantum
Mechanics/Molecular Mechanics. Journal of Physical Chemistry B, v. 109, n. 20, p.
10474–10483, maio 2005. ISSN 1520-6106.
128 RICCARDI, D.; SCHAEFER, P.; CUI, Q. pKa calculations in solution and proteins
with QM/MM free energy perturbation simulations: a quantitative test of QM/MM
protocols. The journal of physical chemistry. B, v. 109, n. 37, p. 17715–33, set. 2005.
ISSN 1520-6106.
129 BARONE, V.; IMPROTA, R.; REGA, N. Computation of protein pK’s values by
an integrated density functional theory/Polarizable Continuum Model approach.
Theoretical Chemistry Accounts: Theory, Computation, and Modeling (Theoretica
Chimica Acta), v. 111, n. 2-6, p. 237–245, mar. 2004. ISSN 1432-881X.
130 JENSEN, J. H. et al. Prediction and rationalization of protein pKa values using
QM and QM/MM methods. The journal of physical chemistry. A, v. 109, n. 30, p.
6634–43, ago. 2005. ISSN 1089-5639.
131 IOZZI, M. F. et al. A polarizable continuum approach for the study of heterogeneous
dielectric environments. The Journal of chemical physics, v. 124, n. 18, p. 184103,
maio 2006. ISSN 0021-9606.
155
132 BENZI, C. et al. A combined theoretical and experimental approach to determining
order parameters of solutes in liquid crystals from 13C NMR data. The journal of
physical chemistry. B, v. 109, n. 7, p. 2584–90, mar. 2005. ISSN 1520-6106.
133 BENZI, C.; COSSI, M.; BARONE, V. Accurate prediction of electron-paramagnetic-
resonance tensors for spin probes dissolved in liquid crystals. The Journal of chemical
physics, v. 123, n. 19, p. 194909, dez. 2005. ISSN 0021-9606.
134 FREDIANI, L. et al. A polarizable continuum model for molecules at diffuse interfaces.
The Journal of chemical physics, v. 120, n. 8, p. 3893–907, fev. 2004. ISSN 0021-9606.
135 FREDIANI, L.; POMELLI, C. S.; TOMASI, J. n-Alkyl alcohols at the water/vapour
and water/benzene interfaces: a study on phase transfer energies. Physical Chemistry
Chemical Physics, v. 2, n. 21, p. 4876–4883, 2000. ISSN 14639076.
136 HOSHI, H. et al. Medium effects on the molecular electronic structure. I. The formu-
lation of a theory for the estimation of a molecular electronic structure surrounded
by an anisotropic medium. Journal of Chemical Physics, v. 87, n. 2, p. 1107, 1987.
ISSN 00219606.
137 RILEY, K. E.; VONDRáSEK, J.; HOBZA, P. Performance of the DFT-D method,
paired with the PCM implicit solvation model, for the computation of interaction
energies of solvated complexes of biological interest. Physical chemistry chemical
physics : PCCP, v. 9, n. 41, p. 5555–60, nov. 2007. ISSN 1463-9076.
138 PARISIO, G. et al. Polarity-sensitive fluorescent probes in lipid bilayers: bridging
spectroscopic behavior and microenvironment properties. The journal of physical
chemistry. B, v. 115, n. 33, p. 9980–9, ago. 2011. ISSN 1520-5207.
139 STERN, H. A.; FELLER, S. E. Calculation of the dielectric permittivity profile
for a nonuniform system: Application to a lipid bilayer simulation. The Journal of
Chemical Physics, v. 118, n. 7, p. 3401, 2003. ISSN 00219606.
140 TANIZAKI, S.; FEIG, M. A generalized Born formalism for heterogeneous dielectric
environments: application to the implicit modeling of biological membranes. The
Journal of chemical physics, v. 122, n. 12, p. 124706, mar. 2005. ISSN 0021-9606.
141 IM, W.; FEIG, M.; BROOKS III, C. L. An implicit membrane generalized born
theory for the study of structure, stability, and interactions of membrane proteins.
Biophysics Journal, v. 85, p. 2900–18, 2003.
142 KESSEL, A.; MUSAFIA, B.; BEN-TAL, N. Continuum solvent model studies of the
interactions of an anticonvulsant drug with a lipid bilayer. Biophysical journal, v. 80,
n. 6, p. 2536–45, jun. 2001. ISSN 0006-3495.
143 BEN-TAL, N. et al. Free-energy determinants of alpha-helix insertion into lipid
bilayers. Biophysical journal, v. 70, n. 4, p. 1803–12, maio 1996. ISSN 0006-3495.
144 BERNèCHE, S.; NINA, M.; ROUX, B. Molecular dynamics simulation of melittin in
a dimyristoylphosphatidylcholine bilayer membrane. Biophysical journal, v. 75, n. 4,
p. 1603–18, out. 1998. ISSN 0006-3495.
145 BIGGIN, P. C. et al. Simulation studies of alamethicin-bilayer interactions. Biophy-
sical journal, v. 72, n. 2 Pt 1, p. 627–36, mar. 1997. ISSN 0006-3495.
156
146 BLOM, J. G.; PELETIER, M. A. A continuum model of lipid bilayers. European
Journal of Applied Mathematics, v. 15, n. 4, p. 487–508, ago. 2004. ISSN 0956-7925.
147 KESSEL, A.; CAFISO, D. S.; BEN-TAL, N. Continuum solvent model calculations
of alamethicin-membrane interactions: thermodynamic aspects. Biophysical journal,
v. 78, n. 2, p. 571–83, fev. 2000. ISSN 0006-3495.
149 SZABO, A.; OSTLUND, N. S. Modern Quantum Chemistry. 1. ed. Mineola, Nova
York: Dover Publications, 1996.
151 FRISCH, M. J. et al. Gaussian 09. Wallingford CT: Gaussian, Inc., 2009.
153 GODBOUT, N. et al. Optimization of Gaussian-type basis sets for local spin density
functional calculations. Part I. Boron through neon, optimization technique and
validation. Canadian Journal of Chemistry, v. 70, n. 2, p. 560–571, fev. 1992. ISSN
0008-4042.
154 SOSA, C. et al. A local density functional study of the structure and vibratio-
nal frequencies of molecular transition-metal compounds. The Journal of Physical
Chemistry, v. 96, n. 16, p. 6630–6636, ago. 1992. ISSN 0022-3654.
156 KOHN, W.; SHAM, L. J. Self-Consistent Equations Including Exchange and Cor-
relation Effects. Physical Review, v. 140, n. 4A, p. A1133–A1138, nov. 1965. ISSN
0031-899X.
157 MARQUES, M. A. L.; BOTTI, S. Oque é e para que serve a Teoria dos Funcionais
da Densidade. Gazeta da Fı́sica, v. 29, p. 10–15, 2006.
160 SCHAFER, M. et al. Systematic study of parameters influencing the action of Rose
Bengal with visible light on bacterial cells: comparison between the biological effect
and singlet-oxygen production. Photochemistry and photobiology, v. 71, n. 5, p.
514–23, maio 2000. ISSN 0031-8655.
157
161 ZHAO, Y.; TRUHLAR, D. G. Density functionals with broad applicability in
chemistry. Accounts of chemical research, v. 41, n. 2, p. 157–67, fev. 2008. ISSN
0001-4842.
162 PARR, R. G.; YANG, W. Density-functional theory of atoms and molecules. Oxford:
Oxford University Press, 1989.
164 ADAMO, C.; BARONE, V. Toward reliable adiabatic connection models free from
adjustable parameters. Chemical Physics Letters, v. 274, n. 1-3, p. 242–250, ago.
1997. ISSN 00092614.
167 LIN, I.-C. et al. Predicting noncovalent interactions between aromatic biomolecules
with London-dispersion-corrected DFT. The journal of physical chemistry. B, v. 111,
n. 51, p. 14346–54, dez. 2007. ISSN 1520-6106.
168 JAVIER LUQUE, F. et al. Continuum solvation models: Dissecting the free energy
of solvation. Physical Chemistry Chemical Physics, v. 5, n. 18, p. 3827, 2003. ISSN
1463-9076.
169 LI, J. et al. Extension of the platform of applicability of the SM5.42R universal sol-
vation model. Theoretical Chemistry Accounts: Theory, Computation, and Modeling
(Theoretica Chimica Acta), v. 103, n. 1, p. 9–63, nov. 1999. ISSN 1432-881X.
170 LI, J. et al. Universal reaction field model based on ab initio Hartree–Fock theory.
Chemical Physics Letters, v. 288, n. 2-4, p. 293–298, maio 1998. ISSN 00092614.
171 GAO, J.; LUQUE, F. J.; OROZCO, M. Induced dipole moment and atomic charges
based on average electrostatic potentials in aqueous solution. The Journal of Chemical
Physics, v. 98, n. 4, p. 2975, 1993. ISSN 00219606.
172 PLIEGO, J. R.; RIVEROS, J. M. Parametrization of the PCM model for calculating
solvation free energy of anions in dimethyl sulfoxide solutions. Chemical Physics
Letters, v. 355, n. 5-6, p. 543–546, abr. 2002. ISSN 00092614.
173 FU, Y. et al. First-principle predictions of absolute pKa’s of organic acids in dimethyl
sulfoxide solution. Journal of the American Chemical Society, v. 126, n. 3, p. 814–22,
jan. 2004. ISSN 0002-7863.
158
175 LUQUE, F. J. et al. Extension of MST/SCRF method to organic solvents:Ab
initio and semiempirical parametrization for neutral solutes in CCl4. Journal of
Computational Chemistry, v. 17, n. 7, p. 806–820, maio 1996. ISSN 0192-8651.
177 CANCES, E.; MENNUCCI, B.; TOMASI, J. A new integral equation formalism
for the polarizable continuum model: Theoretical background and applications to
isotropic and anisotropic dielectrics. The Journal of Chemical Physics, v. 107, n. 8,
p. 3032, 1997. ISSN 00219606.
178 MENNUCCI, B.; CANCèS, E.; TOMASI, J. Evaluation of Solvent Effects in Isotropic
and Anisotropic Dielectrics and in Ionic Solutions with a Unified Integral Equation
Method: Theoretical Bases, Computational Implementation, and Numerical Applica-
tions. The Journal of Physical Chemistry B, v. 101, n. 49, p. 10506–10517, dez. 1997.
ISSN 1520-6106.
179 TOMASI, J.; MENNUCCI, B.; CAMMI, R. Quantum mechanical continuum solva-
tion models. Chemical reviews, v. 105, n. 8, p. 2999–3093, ago. 2005. ISSN 0009-2665.
180 CAMMI, R.; TOMASI, J. Nonequilibrium solvation theory for the polarizable
continuum model: A new formulation at theSCF level with application to the case
of the frequency-dependent linear electric response function. International Journal
of Quantum Chemistry, v. 56, n. S29, p. 465–474, fev. 1995. ISSN 0020-7608.
185 GOUJON, F. et al. Direct Monte Carlo simulations of the equilibrium properties of
n-pentane liquid–vapor interface. The Journal of Chemical Physics, v. 116, n. 18,
p. 8106, 2002. ISSN 00219606.
159
188 WEINER, S. J. et al. A new force field for molecular mechanical simulation of
nucleic acids and proteins. Journal of the American Chemical Society, v. 106, n. 3, p.
765–784, fev. 1984. ISSN 0002-7863.
195 HESS, B. et al. GROMACS 4: Algorithms for Highly Efficient, Load-Balanced, and
Scalable Molecular Simulation. Journal of Chemical Theory and Computation, v. 4,
n. 3, p. 435–447, mar. 2008. ISSN 1549-9618.
198 DARDEN, T.; YORK, D.; PEDERSEN, L. Particle mesh Ewald: An Nlog(N) method
for Ewald sums in large systems. The Journal of Chemical Physics, v. 98, n. 12, p.
10089, 1993. ISSN 00219606.
200 SKJEVIK, A. A. et al. LIPID11: a modular framework for lipid simulations using
amber. The journal of physical chemistry. B, v. 116, n. 36, p. 11124–36, set. 2012.
ISSN 1520-5207.
201 MACKERELL, A. D. et al. All-atom empirical potential for molecular modeling and
dynamics studies of proteins. The journal of physical chemistry. B, v. 102, n. 18, p.
3586–616, abr. 1998. ISSN 1520-6106.
202 WANG, J. et al. Development and testing of a general amber force field. Journal of
computational chemistry, v. 25, n. 9, p. 1157–74, jul. 2004. ISSN 0192-8651.
160
203 OOSTENBRINK, C. et al. A biomolecular force field based on the free enthalpy
of hydration and solvation: the GROMOS force-field parameter sets 53A5 and
53A6. Journal of computational chemistry, v. 25, n. 13, p. 1656–76, out. 2004. ISSN
0192-8651.
206 FELLER, S. E. et al. Constant pressure molecular dynamics simulation: The Langevin
piston method. The Journal of Chemical Physics, v. 103, n. 11, p. 4613, 1995. ISSN
00219606.
207 MARTYNA, G. J. et al. Explicit reversible integrators for extended systems dynamics.
Molecular Physics, v. 87, n. 5, p. 1117–1157, abr. 1996. ISSN 0026-8976.
211 ÖHRN, A.; KARLSTRöM, G. A theoretical study of the solvent shift to the transition
in formaldehyde with an effective discrete quantum chemical solvent model including
non-electrostatic perturbation. Molecular Physics, v. 104, n. 19, p. 3087–3099, out.
2006. ISSN 0026-8976.
212 LIN; GAO, J. Solvatochromic Shifts of the n fi* Transition of Acetone from Steam
Vapor to Ambient Aqueous Solution: A Combined Configuration Interaction QM/MM
Simulation Study Incorporating Solvent Polarization. Journal of Chemical Theory
and Computation, v. 3, n. 4, p. 1484–1493, jul. 2007. ISSN 1549-9618.
161
Chemical Theory and Computation, v. 8, n. 7, p. 2197–2203, jul. 2012. ISSN 1549-
9618.
216 KAILA, V. R. I.; WIKSTROM, M.; HUMMER, G. Electrostatics, hydration, and
proton transfer dynamics in the membrane domain of respiratory complex I. Procee-
dings of the National Academy of Sciences, v. 111, n. 19, p. 6988–6993, maio 2014.
ISSN 0027-8424.
217 WARSHEL, A.; LEVITT, M. Theoretical Studies of Enzymic Reactions: Dielectric,
Electrostatic and Steric Stabilization of the Carbonium Ion in the Reaction of
Lysozyme. Journal of Molecular Biology, v. 103, p. 227–249, 1976.
218 SENN, H. M.; THIEL, W. QM/MM Methods for Biomolecular Systems. Angewandte
Chemie International Edition, v. 48, n. 7, p. 1198–1229, fev. 2009. ISSN 14337851.
219 FIELD, M. J.; BASH, P. A.; KARPLUS, M. A combined quantum mechanical
and molecular mechanical potential for molecular dynamics simulations. Journal of
Computational Chemistry, v. 11, n. 6, p. 700–733, jul. 1990. ISSN 0192-8651.
220 SINGH, U. C.; KOLLMAN, P. A. A combinedab initio quantum mechanical and
molecular mechanical method for carrying out simulations on complex molecular
systems: Applications to the CH3Cl + Cl? exchange reaction and gas phase
protonation of polyethers. Journal of Computational Chemistry, v. 7, n. 6, p. 718–
730, dez. 1986. ISSN 0192-8651.
221 CURUTCHET, C. et al. Electronic Energy Transfer in Condensed Phase Studied by
a Polarizable QM/MM Model. Journal of Chemical Theory and Computation, v. 5,
n. 7, p. 1838–1848, jul. 2009. ISSN 1549-9618.
222 LIN, I.-C. et al. Library of dispersion-corrected atom-centered potentials for genera-
lized gradient approximation functionals: Elements H, C, N, O, He, Ne, Ar, and Kr.
Physical Review B, v. 75, n. 20, p. 205131, maio 2007. ISSN 1098-0121.
223 CASE, D. A. et al. AMBER 12. São Francisco: University of California, 2012.
224 CASE, D. A. et al. AmberTools13. São Francisco: Universidade da California, 2012.
225 GUNDERTOFTE, K. A comparison of conformational energies calculated by several
molecular mechanics methods. Journal of computational chemistry, v. 17, n. 4, p.
429–449, 1996.
226 HOBZA, P. et al. Performance of empirical potentials (AMBER, CFF95, CVFF,
CHARMM, OPLS, POLTEV), semiempirical quantum chemical methods (AM1,
MNDO/M, PM3), and ab initio Hartree–Fock method for interaction of DNA bases:
Comparison with nonempirical beyond Hartree–Fock res. Journal of computational
chemistry, v. 18, n. 9, p. 1136–1150, 1997.
227 RóG, T.; MURZYN, K.; PASENKIEWICZ-GIERULA, M. Molecular dynamics simu-
lations of charged and neutral lipid bilayers: treatment of electrostatic interactions.
Acta biochimica Polonica, v. 50, n. 3, p. 789–98, jan. 2003. ISSN 0001-527X.
228 CORNELL, W. D. et al. A Second Generation Force Field for the Simulation of
Proteins, Nucleic Acids, and Organic Molecules. Journal of the American Chemical
Society, v. 117, n. 19, p. 5179–5197, maio 1995. ISSN 0002-7863.
162
229 JóJáRT, B.; MARTINEK, T. A. Performance of the general amber force field in
modeling aqueous POPC membrane bilayers. Journal of computational chemistry,
v. 28, n. 12, p. 2051–8, set. 2007. ISSN 0192-8651.
230 THOMPSON, M. a. QM/MMpol: A Consistent Model for Solute/Solvent Polari-
zation. Application to the Aqueous Solvation and Spectroscopy of Formaldehyde,
Acetaldehyde, and Acetone. The Journal of Physical Chemistry, v. 100, p. 14492–
14507, 1996. ISSN 0022-3654.
231 JO, S. et al. CHARMM-GUI: a web-based graphical user interface for CHARMM.
Journal of computational chemistry, v. 29, n. 11, p. 1859–65, ago. 2008. ISSN
1096-987X.
232 JO, S. et al. CHARMM-GUI Membrane Builder for mixed bilayers and its application
to yeast membranes. Biophysical journal, v. 97, n. 1, p. 50–8, jul. 2009. ISSN 1542-
0086.
233 CORMANICH, R. A.; NUNES, C. A.; FREITAS, M. P. Desenhos de estruturas
quı́micas correlacionam-se com propriedades biológicas: MIA-QSAR. Quı́mica Nova,
v. 35, n. 6, p. 1157–1163, 2012. ISSN 0100-4042.
234 HASEL, W.; HENDRICKSON, T. F.; STILL, W. C. A rapid approximation to the
solvent accessible surface areas of atoms. Tetrahedron Comp. Meth., v. 1, p. 103–116,
1988.
235 SETúBAL, C. A. Procura por novos fotossensibilizadores para uso em terapia
fotodinâmica. Setor de Ciências Exatas da Universidade Federal do Paraná, Exatas
da Universidade Federal do Paraná, Curitiba, 2007.
236 STERNBERG, E. D.; DOLPHIN, D.; BRUCKNER, C. Porphyrin-based photosensi-
tizers for use in photodynamic therapy. Tetrahedron, v. 54, n. 17, p. 4151–4202, abr.
1998. ISSN 00404020.
237 PARRINELLO, M.; ANDREONI, W. CPMD - Carr Parrinello Molecular Dynamics.
[S.l.]: http://www.cpmd.org/„ 2000.
238 COWLEY, A. C. et al. Measurement of repulsive forces between charged phospholipid
bilayers. Biochemistry, v. 17, n. 15, p. 3163–3168, jul. 1978. ISSN 0006-2960.
239 POLYANSKY, A. A.; VOLYNSKY, P. E.; EFREMOV, R. G. Computer simulations
of membrane-lytic peptides: perspectives in drug design. Journal of bioinformatics
and computational biology, v. 5, n. 2B, p. 611–26, abr. 2007. ISSN 0219-7200.
240 EDWARDS, L.; GOUTERMAN, M. Porphyrins XV. Vapor absorption spectra and
stability: Phatalocyanines. Journal of Molecular Spectroscopy, v. 33, n. 2, p. 292–310,
fev. 1970. ISSN 00222852.
241 MELO, T. B.; REISAETER, G. The Physicochemical State of Protoporphyrin IX
in aqueous solution investigated by Fluorescence and Light Scattering. Biophysical
Chemistry, v. 25, p. 99–104, 1986.
242 ATIF, M. A fluorescence micro-spectroscopy technique for the study of intracellular
photobleaching of mTHPC. Journal of Optoelectronics and Advanced Materials, v. 16,
n. 7-8, p. 798–803, 2014.
163
243 AQUARON, R. et al. Simple, reliable and fast spectrofluorometric method for
determination of plasma Verteporfin (Visudyne) levels during photodynamic therapy
for choroidal neovascularization. Cellular and molecular biology (Noisy-le-Grand,
France), v. 48, p. 925–930, 2002. ISSN 01455680.
244 VAN MINGROOT, H. et al. The emission at 669 nm of metal free phthalocyanine
in toluene and 1-bromonaphthalene solutions. 1996, 397–402 p.
245 DAS, R. R. Absorption and Emission Spectral Studies on the Dimerization of Free
Porphyrin and its Zinc(II), COPPER(II) and Nickel(II) Derivatives in Aqueous
Solution. Journal of Inorganic Nuclear Chemistry, v. 37, p. 153–157, 1975.
164
A Apêndice
Qx Qy Soret Qy -Soret
Metodologia
⁄ (nm) Erro (%) ⁄ (nm) Erro (%) ⁄ (nm) Erro (%) ⁄ (nm) Erro (%)
B3LYP/6-31g(d) 600,5 12,5 593,5 4,7 328,0 3,5 265,5 6,2
B3PW91/6-31g(d) 597,5 12,9 590,0 5,3 337,0 0,9 253,0 10,6
M06-L/6-31g(d) 609,0 11,2 603,0 3,2 338,5 0,4 264,5 6,5
B3LYP/DGDZVP 610,0 11,1 603,0 3,2 332,0 2,4 271,0 4,2
B3LYP/LanL2DZ 600,0 12,5 598,0 4,0 349,0 -2,6 249,0 12,0
PM6 783,3 -14,2 537,9 13,7 349,5 -2,8 188,4 33,4
Experimental 686,0 - 623,0 - 340,0 - 283,0 -
Qx Soret Qx -Soret
Metodologia
⁄ (nm) Erro (%) ⁄ (nm) Erro (%) ⁄ (nm) Erro (%)
B3LYP/6-31g(d) 593,0 11,1 331,0 3,6 262 19,0
B3PW91/6-31g(d) 587,5 11,9 328,5 4,4 259,0 19,9
M06-L/6-31g(d) 602,0 9,7 370,0 -7,7 232,0 28,3
B3LYP/DGDZVP 601,5 9,8 320 6,8 281,5 13,0
B3LYP/LanL2DZ 596,0 10,6 342,0 0,4 254,0 21,5
PM6 650,0 2,5 403,0 -17,3 247,0 21,9
Experimental 667,0 - 343,5 - 323,5 -
F Zn 4H F 1nh3
F 1nh3 4c F 2nh3
166
F 2nh3t F 3nh3
F 4nh3 F nh2ch3
F n3ch3 F 1oacc0
167
F n2acc0 F Zn ac2e
F Zn bc2e F Zn cd2e
F Zn c2ef F Zn c2eg
168
F Zn c2eh
Figura 72: Estruturas das ftalocianinas otimizadas por DFT/B3LYP/LanL2DZ em fase gasosa.
Esfera cinza = carbono, branca = hidrogênio, azul = nitrogênio, vermelha =
oxigênio, rosa = alumı́nio, roxa = zinco e verde = cloro.
Qx Qy Soret Qx-Soret
Metodologia
⁄ (nm) Erro (%) ⁄ (nm) Erro (%) ⁄ (nm) Erro (%) ⁄ (nm) Erro (%)
B3LYP/6-31g(d) 600,5 12,5 593,5 4,7 328,0 3,5 265,5 6,2
B3PW91/6-31g(d) 597,5 12,9 590,0 5,3 337,0 0,9 253,0 10,6
M06-L/6-31g(d) 609,0 11,2 603,0 3,2 338,5 0,4 264,5 6,5
OLYP/6-31g(d) 626,0 8,7 617,0 1,0 357,0 8,7 260,0 8,1
B3LYP/DGDZVP 610,0 11,1 603,0 3,2 332,0 2,4 271,0 4,2
B3LYP/LanL2DZ 600,0 12,5 598,0 4,0 349,0 -2,6 249,0 12,0
ZINDO 853,0 -24,3 798,0 -28,1 349,0 -2,6 449,0 -58,7
OLYP/6-31g(d) com
644,0 6,1 633,0 -1,6 367,0 -7,9 266,0 6,0
geometria B3LYP/LanL2DZ
Experimental 686,0 - 623,0 - 340,0 - 283,0 -
Tabela 22: Comprimentos de onda ⁄ (nm) experimental (em fase gasosa [29]) e calculado com
diversas metodologias em fase gasosa, referentes as bandas Qx e Soret, diferença
entre Qx e Soret e os erros entre comprimentos de onda experimental e calculado,
para F Zn 4H
Qx Soret Qx-Soret
Metodologia
⁄ (nm) Erro (%) ⁄ (nm) Erro (%) ⁄ (nm) Erro (%)
B3LYP/6-31g(d) 593,0 11,1 331,0 3,6 262 19,0
B3PW91/6-31g(d) 587,5 11,9 328,5 4,4 259,0 19,9
M06-L/6-31g(d) 602,0 9,7 370,0 -7,7 232,0 28,3
OLYP/6-31g(d) 618 7,3 353 -2,8 265 18,1
B3LYP/DGDZVP 601,5 9,8 320 6,8 281,5 13,0
B3LYP/LanL2DZ 596,0 10,6 342,0 0,4 254,0 21,5
ZINDO 882,0 -32,2 297,0 -13,5 585,0 -80,8
OLYP/6-31g(d) com
623 6,6 370,5 -7,9 252,5 21,9
geometria B3LYP/LanL2DZ
Experimental 667,0 - 343,5 - 323,5 -
Tabela 23: Comprimentos de onda ⁄ (nm) experimental (em THF) e calculado com diversas
metodologias (fase gasosa), referentes as bandas Qx e Soret, diferença entre Qx e
Soret e os erros entre comprimentos de onda experimental e calculado, para F Zn 4H
169
P 4H P Al 4H
P Zn 4H P Si 2Cl
P 4fenil P 4b1cA
170
P 4b1cB P 4b2cAC
P 4b2cAT P 4b3cA
P 4b4cA P 4b4cB
171
Molécula Carga Gf g Gag Gsolv
F4H 0 -1667,639 -1667,671 -0,032
F Al 4H 0 -1683,639 -1683,679 -0,040
F Zn 4H 0 -1732,158 -1732,207 -0,049
F 1nh3 1 -1801,884 -1802,015 -0,869
F 1nh3 4c 1 -1880,445 -1880,581 -0,136
F 2nh3 2 -1936,090 -1936,358 -0,268
F 2nh3t 2 -1936,092 -1936,360 -0,268
F 3nh3 3 -2070,260 -2070,706 -0,446
F 4nh3 4 -2204,387 -2205,048 -0,661
F nh2ch3 1 -1880,442 -1880,553 -0,111
F n3ch3 1 -1919,720 -1919,822 -0,102
F 1oacc0 0 -1970,644 -1970,698 -0,054
F n2acc0 0 -2217,848 -2217,907 -0,059
F Zn ac2e 0 -2723,228 [1] [1]
F Zn bc2e 0 -2762,516 -2762,619 -0,103
F Zn cd2e 0 -2837,714 -2837,826 -0,112
F Zn c2ef 0 -3003,448 -3003,564 -0,116
F Zn c2eg 0 -3078,591 -3078,708 -0,117
F Zn c2eh 0 -3183,335 -3183,455 -0,120
P 4H 0 -989,132 -989,158 -0,026
P Al 4H 0 -1005,126 -1005,165 -0,039
P Zn 4H 0 -1053,652 -1053,689 -0,037
P Si 4H 0 -1021,866 -1021,911 -0,044
P 4fenil 0 -1912,885 -1912,914 -0,029
P 4b1cA 1 -1968,588 -1968,670 -0,082
P 4b1cB 1 -1968,584 -1968,657 -0,073
P 4b2cAC 2 -2024,225 -2024,429 -0,205
P 4b2cAT 2 -2024,230 -2024,428 -0,198
P 4b3cA 3 -2079,809 -2080,170 -0,361
P 4b4cA 4 -2135,332 -2135,939 -0,607
P 4b4cB 4 -2135,327 -2135,940 -0,613
P Si 2ch3 0 -1222,107 -1222,139 -0,032
P Si 2fenil 0 -1605,441 -1605,483 -0,041
C 4H 0 -990,326 -990,348 -0,022
C Zn 4H 0 -1054,839 -1054,876 -0,036
Foscan 0 -2213,725 -2213,801 -0,76
Levulan 0 -1837,527 -1837,579 -0,052
Visudyne 0 -2367,011 -2367,054 -0,043
Tabela 24: Energia livre de Gibbs em fase gasosa ( Gf g ) e em água ( Gag ) e diferença entre
elas ( Gsolv ) de todas as moléculas estudadas. [1] Cálculo não convergido.
172
P Si 2ch3 P Si 2fenil
Levulan Visudyne
Figura 73: Estruturas das porfirinas otimizadas por DFT/B3LYP/LanL2DZ em fase gasosa.
Esfera cinza = carbono, branca = hidrogênio, azul = nitrogênio, vermelha =
oxigênio, rosa = alumı́nio, roxa = zinco, verde = cloro e amarelo=silı́cio
173
Foscan C 4H
C Zn 4H
Figura 74: Estruturas das clorinas otimizadas por DFT/B3LYP/LanL2DZ em fase gasosa.
Esfera cinza = carbono, branca = hidrogênio, azul = nitrogênio e roxa = zinco.
174
Fase gasosa Água Clorofórmio
Molécula
EH EL H≠L (u.a.) H≠L (eV) EH EL H≠L (u.a.) H≠L (eV) EH EL H≠L (u.a.) H≠L (eV)
F 4H -0,16569 -0,11448 0,051 1,394 -0,17049 -0,12050 0,050 1,360 -0,16809 -0,11762 0,050 1,373
F Al 4H -0,17127 -0,12169 0,050 1,349 -0,17480 -0,12647 0,048 1,315 -0,17293 -0,12412 0,049 1,328
F Zn 4H -0,16447 -0,11232 0,052 1,419 -0,16612 -0,11409 0,052 1,416 -0,16521 -0,11317 0,052 1,416
F 1nh3 -0,22875 -0,18228 0,046 1,265 -0,17202 -0,12152 0,051 1,374 -0,18071 -0,12950 0,051 1,394
F 1nh3 4c -0,21453 -0,18905 0,025 0,693 -0,17054 -0,12067 0,050 1,357 -0,17668 -0,12609 0,051 1,377
F 2nh3 -0,28979 -0,23905 0,051 1,381 -0,17350 -0,12344 0,050 1,362 -0,19308 -0,14244 0,051 1,378
F 2nh3t -0,28848 -0,24106 0,047 1,290 -0,17339 -0,12295 0,050 1,373 -0,19276 -0,14215 0,051 1,377
F 3nh3 -0,34894 -0,29883 0,050 1,364 -0,17502 -0,12485 0,050 1,365 -0,20520 -0,15474 0,050 1,373
F 4nh3 -0,40676 -0,35500 0,052 1,408 -0,17652 -0,12626 0,050 1,368 -0,21894 -0,16804 0,051 1,385
175
F nh2ch3 -0,22637 -0,17605 0,050 1,369 -0,17207 -0,12161 0,050 1,373 -0,18078 -0,12957 0,051 1,394
F 3nh3 -0,22564 -0,17525 0,050 1,371 -0,17214 -0,12167 0,050 1,373 -0,18082 -0,12962 0,051 1,393
F 1oacc0 -0,16338 -0,11331 0,050 1,363 -0,16918 -0,12033 0,049 1,329 -0,16641 -0,11710 0,049 1,342
F n2acc0 -0,16518 -0,11406 0,051 1,391 -0,16094 -0,12016 0,041 1,110 -0,16263 -0,11731 0,045 1,233
F Zn ac2e -0,15454 -0,10552 0,049 1,334 -0,16078 -0,11208 0,049 1,325 -0,15832 -0,10954 0,049 1,327
F Zn bc2e -0,15391 -0,10492 0,049 1,333 -0,15993 -0,11128 0,049 1,324 -0,15794 -0,10921 0,049 1,326
F Zn cd2e -0,15393 -0,10494 0,049 1,333 -0,16018 -0,11148 0,049 1,325 -0,15803 -0,10928 0,049 1,327
F Zn c2ef -0,15670 -0,10770 0,049 1,333 -0,16056 -0,11181 0,049 1,327 -0,15906 -0,11026 0,049 1,328
F Zn c2eg -0,15762 -0,10858 0,049 1,334 -0,16065 -0,11187 0,049 1,327 -0,15940 -0,11057 0,049 1,329
F Zn c2eh -0,15395 -0,10495 0,049 1,333 -0,16025 -0,11153 0,049 1,326 -0,15815 -0,10938 0,049 1,327
Tabela 25: Energias (u.a.) dos orbitais HOMO ( H ), LUMO ( L) e diferença entre eles ( H≠L ) em u.a. e eV, para todas as ftalocianinas estudadas
Fase gasosa Água Clorofórmio
Molécula
EH EL H≠L (u.a.) H≠L (eV) EH EL H≠L (u.a.) H≠L (eV) EH EL H≠L (u.a.) H≠L (eV)
P 4H -0,16453 -0,09435 0,070 1,910 -0,17056 -0,09903 0,072 1,946 -0,16751 -0,09640 0,071 1,935
P Al 4H -0,17828 -0,10324 0,075 2,042 -0,18184 -0,10493 0,077 2,093 -0,17965 -0,10329 0,076 2,078
P Zn 4H -0,16749 -0,09173 0,076 2,062 -0,16752 -0,09238 0,075 2,045 -0,16726 -0,09182 0,075 2,053
P Si 4H -0,17947 -0,10899 0,070 1,918 -0,18583 -0,11220 0,074 2,004 -0,18229 -0,10974 0,073 1,974
P 4fenil -0,15659 -0,09171 0,065 1,766 -0,16631 -0,10140 0,065 1,766 -0,16236 -0,09744 0,065 1,767
P 4b1cA -0,24947 -0,20885 0,041 1,105 -0,17704 -0,12381 0,053 1,448 -0,19011 -0,14103 0,049 1,336
P 4b1cB -0,24393 -0,21836 0,026 0,696 -0,17410 -0,11632 0,058 1,572 -0,18633 -0,13913 0,047 1,284
P 4b2cAC -0,32756 -0,29173 0,036 0,975 -0,18395 -0,12953 0,054 1,481 -0,21269 -0,16390 0,049 1,328
P 4b2cAT -0,32914 -0,28990 0,039 1,068 -0,18468 -0,13245 0,052 1,421 -0,21229 -0,16369 0,049 1,323
P 4b3cA -0,40114 -0,35878 0,042 1,153 -0,19052 -0,13344 0,057 1,553 -0,23226 -0,17818 0,054 1,472
P 4b4cA -0,48503 -0,42624 0,059 1,600 -0,19734 -0,13589 0,061 1,672 -0,25435 -0,19391 0,060 1,645
176
P 4b4cB -0,48563 -0,42223 0,063 1,725 -0,19527 -0,12944 0,066 1,791 -0,25322 -0,18820 0,065 1,769
P Si 2ch3 -0,16701 -0,09505 0,072 1,958 -0,17700 -0,10365 0,073 1,996 -0,17311 -0,10029 0,073 1,982
P Si 2fenil -0,16449 -0,10186 0,063 1,704 -0,18041 -0,10714 0,073 1,994 -0,17557 -0,10437 0,071 1,937
C 4H -0,15808 -0,09214 0,066 1,794 -0,16388 -0,09761 0,066 1,506 -0,16101 -0,09488 0,066 1,800
C Zn 4H -0,15494 -0,09033 0,065 1,758 -0,15719 -0,09204 0,065 1,803 -0,15611 -0,09124 0,065 1,765
Foscan -0,15318 -0,08761 0,066 1,784 -0,16613 -0,10078 0,065 1,226 -0,16065 -0,09515 0,066 1,782
Levulan -0,15914 -0,09070 0,068 1,862 -0,16951 -0,09946 0,070 1,778 -0,16531 -0,09589 0,069 1,889
Visudyne -0,16383 -0,09645 0,067 1,834 -0,17326 -0,10471 0,068 1,906 -0,17002 -0,10178 0,068 1,857
Tabela 26: Energias (u.a.) dos orbitais HOMO ( H ), LUMO ( L) e diferença entre eles ( H≠L ) em u.a. e eV, para todas as porfirinas e clorinas
estudadas
Dipolo de transição
Molécula Carga Energia ⁄ (nm) f Orbitais Contribuição
x y z
1,93 644 0,2696 133-135 0,56180 2,3901 0,0011 0,0000
F 4H 0
1,96 633 0,3341 133-134 0,56505 -0,0015 2,6381 0,0000
147-148 -0,1189 2,4751 0,0784
F Al 4H 0 1,89 657 0,2843 0,56735
147-149 -0,0032 0,0781 -2,4781
147-148 2,4992 -0,5852 0,0000
F Zn 4H 0 1,96 634 0,3159 0,56474
147-149 -0,5851 -2,4991 0,0000
1,86 665 0,3151 145-147 0,56069 -2,4579 0,9201 0,0857
F 1nh3 1
1,95 635 0,2863 145-148 0,56921 -0,8786 -2,2831 0,0525
1,89 654 0,3177 153-155 0,56250 2,4948 0,7537 0,2333
F 1nh3 4c 1
1,96 632 0,3168 153-156 0,56304 0,7548 -2,4525 0,1016
1,90 654 0,2813 157-158 0,49496 2,3762 -0,6255 -0,1450
F 2nh3 2
1,93 643 0,3102 157-159 0,51102 -0,9140 -2,3935 0,0662
1,78 695 0,3254 157-158 0,56288 2,6327 0,7070 -0,1278
F 2nh3t 2
1,97 628 0,3005 157-159 0,56195 0,6204 -2,4133 -0,0491
1,89 657 0,3530 169-170 0,56738 2.7245 0,4509 -0,1025
F 3nh3 3
1,96 633 0,2814 169-171 0,55693 0,4125 -2,3837 -0,1103
1,93 643 0,3336 181-183 0,64435 1,5393 0,7685 0,1430
F 4nh3 4
1,96 634 0,3973 181-182 0,66195 0,9222 -2,7269 0,0798
1,87 663 0,3161 153-154 0,56416 -2,4613 0,9073 0,1141
F nh2ch3 1
1,95 634 0,2915 153-155 0,56753 -0,8887 -2,3008 0,0660
1,87 662 0,3180 157-158 0,56482 -2,4840 0,8518 0,1754
F n3ch3 1
1,96 634 0,2987 157-159 0,56885 -0,8333 -2,3513 0,0845
1,90 654 0,2816 152-153 0,46994 -2,0499 1,3644 -0,0001
F 1oacc0 0
1,93 641 0,3337 152-154 0,47355 1,4736 2,2077 0,0001
1,92 646 0,3168 171-173 0,54233 2,5900 0,0389 -0,1591
F n2acc0 0
1,96 633 0,3330 171-172 0,54319 0,0466 -2,6331 0,0099
1,82 683 0,2869 215-216 0,56615 -2,3258 1,0146 -0,0932
F Zn ac2e 0
1,87 662 0,3588 215-217 0,57400 -1,1259 -2,5561 0,1406
1,81 684 0,2882 219-220 0,56577 -2,4840 0,5542 -0,0838
F Zn bc2e 0
1,87 663 0,3586 219-221 0,57406 -0,5811 -2,7323 0,1442
1,82 683 0,2949 223-224 0,56670 -2,5353 0,4409 -0,0802
F Zn cd2e 0
1,87 662 0,3613 223-225 0,57416 -0,4402 -2,7677 0,1389
1,82 682 0,3107 236-237 0,67175 -2,5783 0,5713 -0,0840
F Zn c2ef 0
1,87 661 0,3578 236-238 0,67703 -0,5594 -2,7311 0,1304
1,82 682 0,3149 240-241 0,56839 2,6328 -0,3622 -0,0660
F Zn c2eg 0
1,88 661 0,3613 240-242 0,57338 0,3346 2,7809 0,1228
1,82 683 0,3031 250-251 0,56682 2,6075 0,0536 -0,0974
F Zn c2eh 0
1,87 662 0,3589 250-252 0,57416 -0,1129 2,7907 0,1552
Tabela 27: Energia (eV), comprimento de onda ⁄ (nm), força do oscilador (f ), orbitais envolvidos
e suas contribuições, e dipolo das transições, calculadas via DFT/OLYP/6-31g(d)
em fase gasosa, responsáveis pelas bandas Q das ftalocianinas estudadas
177
Dipolo de transição
Molécula Carga Energia ⁄ (nm) f Orbitais Contribuição
x y z
81-83 0,53808
2,12 585 0,0010 -0,1391 0,0025 -0,0001
80-82 -0,42134
P 4H 0
81-82 0,49519
2,26 548 0,0004 0,0022 0,0806 0,0000
80-83 0,44031
95-97 0,49823
0,0004 0,0305 0,0000
94-96 0,45127
P Al 4H 0 2,24 554 0,0001
95-96 0,49822
0,0308 -0,0005 0,0001
94-97 -0,45127
95-97 0,48829
-0,0579 0,0350 -0,0002
94-96 0,46218
P Zn 4H 0 2,27 546 0,0002
95-96 0,48830
-0,0349 -0,0570 -0,0030
94-97 -0,46217
104-105 0,57200
0,0051 0,1528 0,0000
103-106 -0,41013
P Si 4H 0 2,1634 573 0,0012
104-106 0,57191
0,1522 -0,0050 0,0000
103-105 0,41026
161-163 0,55364
1,99 622 0,0244 -0,0004 -0,7069 0,0000
160-162 0,39844
P 4fenil 0
161-162 0,52259
2,12 584 0,0393 -0,8698 0,0002 0,0193
160-163 0,40523
1,39 889 0,0165 164-166 0,67875 -0,2023 -0,6643 0,0115
1,59 778 0,2779 165-166 0,64608 2,6610 -0,1824 0,0092
1,77 700 0,0498 163-166 0,67961 1,0709 0,0242 0,0017
2,14 581 0,0225 165-167 0,50219 -0,2079 -0,6224 -0,0106
P 4b1cA 1
156-166 0,50192
2,23 556 0,0147 -0,4446 -0,2680 -0,0099
165-169 -0,38538
2,31 537 0,0440 165-168 0,47971 0,7997 -0,3712 -0,0090
2,35 528 0,0322 164-169 0,67155 -0,7476 0,0299 -0,0032
0,80 1544 0,0160 165-166 0,67224 -0,8467 0,3039 -0,0575
1,03 1202 0,0258 164-166 0,66785 -1,0099 0,0091 -0,0128
1,73 718 0,1423 165-167 0,50245 -1,8072 0,3141 0,0022
P 4b1cB 1
1,80 688 0,0069 164-167 0,65198 -0,1639 -0,3589 -0,0322
2,11 588 0,0112 165-169 0,51161 0,0200 -0,4640 -0,0257
2,22 560 0,0211 165-168 0,44071 -0,2342 -0,5925 0,0300
1,23 1008 0,1336 169-170 0,66633 2,1049 0,0217 0,0163
1,27 978 0,0034 168-170 0,65429 0,0803 -0,3202 -0,0027
P 4b2cAC 2
1,53 812 0,0028 168-171 0,68880 0,2032 -0,1838 -0,0140
1,55 801 0,1707 169-171 0,55439 0,0194 2,1217 -0,0042
1,34 925 0,1143 169-170 0,52133 1,7090 0,7488 0,0304
1,42 871 0,1609 168-170 0,50861 2,1028 -0,4362 -0,0051
P 4b2cAT 2
1,58 783 0,0000 168-171 0,53698 -0,0081 0,0032 0,0003
1,69 735 0,0005 169-171 0,51805 -0,1035 0,0242 -0,0006
1,26 983 0,0696 173-174 0,65759 1,0619 0,9123 0,0223
1,28 966 0,0020 172-174 0,67139 0,2066 0,1441 -0,0279
P 4b3cA 3
1,33 931 0,0527 173-175 0,64508 0,7283 -1,0410 0,0112
1,36 910 0,0093 172-175 0,68893 0,1777 0,4972 0,0014
1,78 670 0,0203 177-179 0,52290 -0,0117 0,6823 0,0004
1,81 686 0,1220 177-178 0,62729 1,6593 0,0028 -0,0242
P 4b4cA 4
2,02 614 0,2046 176-178 0,43278 0,0267 2,0339 -0,0010
2,04 607 0,0804 176-179 0,55389 -1,2674 0,0212 0,0159
1,80 688 0,0041 177-178 0,65087 0,1170 0,2777 0,0333
1,86 667 0,0017 177-179 0,58656 0,1305 -0,1380 0,0328
P 4b4cB 4
2,00 619 0,0431 176-178 0,59689 0,9341 -0,0775 0,0007
2,03 612 0,0490 176-179 0,55816 -0,0145 -0,9937 0,0005
2,22 558 0,0008 104-105 0,54231 0,1186 0,0219 -0,0003
P Si 2ch3 0
2,25 551 0,0003 104-106 0,50042 -0,0269 0,0655 -0,0001
1,71 719 0,0081 136-137 0,68973 0,4381 0,0008 0,0047
P Si 2fenil 0
1,72 719 0,0005 136-138 0,69035 0,1103 -0,0011 0,0014
2,00 620 0,0180 177-179 0,55212 -0,3231 0,5138 0,0009
Foscan 0
2,13 583 0,0266 177-178 0,53435 -0,6382 -0,3214 -0,0003
2,00 618 0,0014 151-152 0,43784 -0,0798 0,1460 -0,0014
Levulan 0
2,14 579 0,0113 151-153 0,39344 0,0675 0,4587 -0,0062
1,99 624 0,0085 185-187 0,49276 0,1934 -0,3681 -0,0312
Visudyne 0
2,11 587 0,0213 185-186 0,49857 -0,6329 0,1008 0,0446
Tabela 28: Energia (eV), comprimento de onda ⁄ (nm), força do oscilador (f ), orbitais en-
volvidos e suas contribuições, e dipolo de transição das transições, calculadas via
DFT/OLYP/6-31g(d) em fase gasosa (IEFPCM), responsáveis pelas bandas Q das
porfirinas estudadas
178
Dipolo de transição
Molécula Carga Energia ⁄(nm) f Orbitais Contribuição
x y z
82-83 0,57227
2,18 569 0,0478 -0,0013 0,9463 0,0000
81-84 -0,41032
C 4H 0
81-83 0,55423
2,28 543 0,0039 0,2639 0,0007 0,0000
82-84 0,43529
96-87 0,62023
2,24 553 0,0866 -0,0002 1,2555 0,0000
95-98 -0,33557
C Zn 4H 0
95-97 0,53729
2,31 537 0,0018 -0,1789 -0,0001 0,0001
96-98 0,45206
Tabela 29: Energia (eV), comprimento de onda ⁄ (nm), força do oscilador (f ), orbitais en-
volvidos e suas contribuições, e dipolo de transição das transições, calculadas via
DFT/OLYP/6-31g(d) em fase gasosa (IEFPCM), responsáveis pelas bandas Q das
clorinas estudadas
179
Molécula Zint (D12) E trip ≠ sing (eV) E trip ≠ sing (nm)
20 1,008 1230
12 0,993 1248
8 1,008 1230
4 1,013 1223
0 1,027 1207
F 2nh3t
-4 1,029 1205
-8 1,040 1193
-12 1,047 1184
-16 1,050 1181
-20 1,065 1165
20 1,020 1216
12 0,981 1264
8 1,029 1204
4 1,024 1210
0 1,019 1217
F 3nh3
-4 1,020 1215
-8 1,018 1218
-12 1,018 1218
-16 1,024 1211
-20 1,001 1239
20 1,198 1035
12 1,197 1036
8 1,187 1045
4 1,207 1027
0 1,214 1022
F 4nh3
-4 1,208 1026
-8 1,208 1026
-12 1,230 1008
-16 1,240 1000
-20 1,254 989
180
Molécula Zint (D12) E trip ≠ sing (eV) E trip ≠ sing (nm)
20 1,458 851
12 1,457 851
8 1,455 852
4 1,444 858
0 1,432 866
P 4H
-4 1,417 875
-8 1,414 877
-12 1,409 880
-16 1,405 882
-20 1,406 882
20 1,724 719
12 1,724 719
8 1,721 721
4 1,705 727
0 1,683 737
P Al 4H
-4 1,667 744
-8 1,665 745
-12 1,526 813
-16 1,659 747
-20 1,660 747
20 1,801 688
12 1,801 688
8 1,799 689
4 1,791 692
0 1,765 703
P Zn 4H
-4 1,757 706
-8 1,758 705
-12 1,775 698
-16 1,776 698
-20 1,777 698
20 1,679 738
12 1,679 739
8 1,675 740
4 1,654 750
0 1,621 765
P Si 4H
-4 1,599 775
-8 1,594 778
-12 0,583 783
-16 0,580 784
-20 0,579 785
20 1,300 954
12 1,287 963
8 1,268 978
4 1,248 994
0 1,207 1027
P 4b2cAC
-4 1,120 1107
-8 1,038 1194
-12 1,007 1231
-16 0,970 1278
-20 0,916 1354
20 1,379 899
12 1,355 915
8 1,317 942
4 1,254 989
0 1,137 1090
P 4b2cAT
-4 1,083 1144
-8 1,043 1189
-12 0,993 1249
-16 0,969 1279
-20 0,973 1274
181
Molécula Zint (D12) E trip ≠ sing (eV) E trip ≠ sing (nm)
20 1,766 702
12 1,766 702
8 1,763 703
4 1,756 706
0 1,742 712
P Si 2ch3
-4 1,723 720
-8 1,723 720
-12 1,722 720
-16 1,719 721
-20 1,719 721
20 1,822 681
12 1,822 681
8 1,822 681
4 1,822 680
0 1,822 680
P Si 2fenil
-4 1,817 682
-8 1,819 682
-12 1,814 683
-16 1,803 688
-20 1,804 687
20 1,466 845
12 1,466 845
8 1,466 846
4 1,462 848
0 1,416 876
Levulan
-4 1,392 890
-8 1,389 893
-12 1,386 895
-16 1,384 896
-20 1,383 896
20 1,412 878
12 1,412 878
8 1,412 878
4 1,408 880
0 1,349 919
Visudyne
-4 1,353 917
-8 1,351 917
-12 1,288 963
-16 1,347 920
-20 1,345 922
20 1,469 844
12 1,469 844
8 1,469 844
4 1,468 845
0 1,445 858
C 4H
-4 1,422 872
-8 1,419 874
-12 1,419 874
-16 1,420 873
-20 1,426 870
20 1,399 886
12 1,399 886
8 1,397 888
4 1,392 891
0 1,311 946
Foscan
-4 1,283 967
-8 1,262 982
-12 1,277 971
-16 1,284 966
-20 1,269 977
182
B SÚMULA CURRICULAR
DADOS PESSOAIS
184
BETTANIN, F.; PERUSSI, J. R.; SILVA, A. B.; HOMEM-DE-MELLO, P. Comparação de
geometrias de corantes fotossensı́veis geradas por cálculos DFT com diferentes funcionais.
33o Reunião Anual da Sociedade Brasileira de Quı́mica (SBQ), Águas de Lindóia, SP,
Brasil. 2010.
BUCK, S. T. G.; BETTANIN, F.; HOMEM-DE-MELLO, P.; SILVA, A. B.; IMASATO,
H.; PERUSSI, J. R. Implications in the molecular and electronic properties of Xanthene
dyes on their photodynamic efficiency.15o International Congress in Photobiology (ICP),
Düsseldorf, NW, Alemanha. 2009.
BUCK, S. T. G.; BETTANIN, F.; HOMEM-DE-MELLO, P.; SILVA, A. B.; IMASATO,
H.; PERUSSI, J. R. Correlating molecular and electronic properties with photodynamic
efficiency of xantenes dyes. International Workshop Nanophotonic and Biophotonics,
Recife, PE, Brasil. 2009.
BETTANIN, F.; HOMEM-DE-MELLO, P.; PERUSSI, J. R.; SILVA, A. B. Evaluation
of density-functionals to study the emission spectra of cationic and anionic dyes. 11o
International Conference on Advanced Materials (ICAM), Rio de Janeiro, RJ, Brasil. 2009.
TRABALHOS APRESENTADOS EM ESCOLAS
185
2010.
BETTANIN, F.; PERUSSI, J. R.; SILVA, A. B.; HOMEM-DE-MELLO, P. Comparação
entre os funcionais B3LYP, PW91, LSDA para estudo teórico de propriedades eletrônicas
de fotossensibilizadores. 1o Workshop Ciência e Tecnologia Quı́mica, Universidade Federal
do ABC, Santo André, SP, Brasil. 2009.
CAPÍTULOS DE LIVROS PUBLICADOS
186