Miranda
Miranda
1 Miranda’s Proposal
Los potenciales quı́micos perturbados para la especie donadora (Do) de electrones, el cual corresponde al nu-
cleófilo, y aceptora (Ac) de de electrones, el cual corresponde al electrófilo, tienen las siguientes expresiones:
1
Numerador = (γµDo IDo + ADo )2 (1 + γµAc )2
+ 2(γµDo IDo + ADo )(IAc + γµAc AAc )(1 + γµDo )(1 + γµAc ) (12)
+ (IAc + γµAc AAc )2 (1 + γµDo )2 .
Expresión final:
− (γµDo IDo + ADo )2 (1 + γµAc )2 + 2(γµDo IDo + ADo )(IAc + γµAc AAc )(1 + γµDo )(1 + γµAc ) + (IAc + γµAc AAc )2 (1 + γµDo )2
ωU F =
2ζηDo (IDo − ADo )(1 + γµDo )2 (1 + γµAc )2
(13)
Expresión final simplificada:
2
− (γµDo IDo + ADo )(1 + γµAc ) + (IAc + γµAc AAc )(1 + γµDo )
ωU F = (14)
2ζηDo (IDo − ADo )(1 + γµDo )2 (1 + γµAc )2
Simplificando el denominador
Reescribiendo la expresión
2
− (γµDo IDo + ADo )(1 + γµAc ) + (IAc + γµAc AAc )(1 + γµDo )
ωU F = (18)
4(IDo − ADo )(1 + γµDo )(1 + γµAc )2
Simplificando el numerador
" #2
2 γµDo IDo + ADo IAc + γµAc AAc
(γµDo IDo γµAc ) γµAc AAc )(1 γµDo ) (1+γµDo )2 (1+γµAc )2
+ ADo )(1 + + (IAc + + = +
1 + γµDo 1 + γµAc
(19)
Se llega a la siguiente expresión final simplificada:
2
Do Ac
γµ IDo +ADo IAc +γµ AAc
(1 + γµDo ) 1+γµDo + 1+γµAc
−
ωU F = (20)
4(IDo − ADo )
Si se hace la siguiente consideración:
2
Se simplifica el numerador
2
Numerador = [(γµ IDo + ADo )(1 + γµ ) + (IAc + γµ AAc )(1 + γµ )]
2
= [(1 + γµ )(γµ IDo + ADo + IAc + γµ AAc )] (24)
2 2
= (1 + γµ ) (γµ IDo + ADo + IAc + γµ AAc ) .
Se simplifica el denominador
Denominador = 2ζηDo (IDo − ADo )(1 + γµ )2 (1 + γµ )2 = 2ζηDo (IDo − ADo )(1 + γµ )4 (25)
h i2
γµ (IDo +AAc )+(ADo +IAc )
(1 + γµ ) 1+γµ
−
ωU F = (34)
4(IDo − ADo )
O, simplificando aún más:
2
− [γµ (IDo + AAc ) + (ADo + IAc )]
ωU F = (35)
4(1 + γµ )(IDo − ADo )
Desarrollando los cuadrados,
− γµ2 (IDo + AAc )2 + 2γµ (IDo + AAc )(ADo + IAc ) + (ADo + IAc )2
ωU F = (36)
4(1 + γµ )(IDo − ADo )
− γµ2 IDo
2
+ 2γµ2 IDo AAc + γµ2 A2Ac + 2γµ IDo ADo + 2γµ IDo IAc + 2γµ AAc ADo + 2γµ AAc IAc + A2Do + 2ADo IAc + IAc
2
ωU F =
4(1 + γµ )(IDo − ADo )
(37)
3
− γµ2 IDo
2
+ 2γµ2 IDo AAc + γµ2 A2Ac + 2γµ IDo ADo + 2γµ IDo IAc + 2γµ AAc ADo + 2γµ AAc IAc + A2Do + 2ADo IAc + IAc
2
ωU F =
4(1 + γµ )(IDo − ADo )
(38)
Multiplicamos ambos lados por el denominador para eliminarlo:
− 2 2 2 2 2 2 2
ωU F ·4(1+γµ )(IDo −ADo ) = γµ IDo +2γµ IDo AAc +γµ AAc +2γµ IDo ADo +2γµ IDo IAc +2γµ AAc ADo +2γµ AAc IAc +ADo +2ADo IAc +IAc
(39)
Expandimos el lado izquierdo:
− − −
4ωU F (1 + γµ )(IDo − ADo ) = 4ωU F (IDo − ADo ) + 4ωU F γµ (IDo − ADo ) (40)
Agrupamos los términos que contienen γµ :
− − 2 2 2 2 2
4ωU F (IDo −ADo )+4ωU F γµ (IDo −ADo ) = γµ (IDo +2IDo AAc +AAc )+2γµ (IDo ADo +IDo IAc +AAc ADo +AAc IAc )+ADo +2ADo IAc +IAc
(41)
Reorganizamos la ecuación para agrupar los términos en γµ :
−
γµ2 (IDo
2
+ 2IDo AAc + A2Ac ) + γµ 2(IDo ADo + IDo IAc + AAc ADo + AAc IAc ) − 4ωU
F (IDo − ADo ) + (42)
−
2 2
ADo + 2ADo IAc + IAc − 4ωU F (IDo − ADo ) = 0 (43)
(44)
2
a = IDo + 2IDo AAc + A2Ac (46)
−
b = 2(IDo ADo + IDo IAc + AAc ADo + AAc IAc ) − 4ωU F (IDo − ADo ) (47)
−
c = A2Do + 2ADo IAc + IAc
2
− 4ωU F (IDo − ADo ) (48)
−
2
4ωU F (IDo − ADo ) − 2IDo ADo − 2IDo IAc − 2AAc ADo − 2AAc IAc − (50)
2 2 − 2 2
4 −(IDo + 2IDo AAc + AAc ) 4ωU F (IDo − ADo ) − ADo − 2ADo IAc − IAc (51)
−
2
α = 4ωU F (IDo − ADo ) − 2IDo ADo − 2IDo IAc − 2AAc ADo − 2AAc IAc − (52)
2 2 − 2 2
4 −(IDo + 2IDo AAc + AAc ) 4ωU F (IDo − ADo ) − ADo − 2ADo IAc − IAc (53)
− √
−4ωU F (IDo − ADo ) − 2IDo ADo − 2IDo IAc − 2AAc ADo − 2AAc IAc α
γµ = 2 2 ± (54)
−2(IDo + 2IDo AAc + AAc ) −2(IDo + 2IDo AAc + A2Ac )
2
4
2 Other approach minimization
Se tiene que las expresiones perturbadas para el potencial quı́mico y la dureza quı́mica son
(µ−
Do )
−
2
ωU F =− (57)
2ηDo
Introduciendo las expresiones para los parámetros perturbados
2
γ Do IDo +ADo
− µ 1+γ Do
µ
−
ωU F = (58)
2ζηDo (IDo − ADo )
Se tiene que
!2 !2
γµDo IDo + ADo γµDo IDo + ADo
− = (59)
1 + γµDo 1 + γµDo
!2
γµDo IDo + ADo (γµDo IDo + ADo )2
= (60)
1 + γµDo (1 + γµDo )2
Do
(γµ IDo +ADo )2
Do )2
(1+γµ
−
ωU F = (61)
2ζηDo (IDo − ADo )
ζηDo = 1 (63)
5
2.2 Otra consideración
Si consideramos que γ y η son diferentes a 1, se tiene que
6
3.1 Definición de g(γµDo )
Definimos g(γµDo ) como:
d u u′ v − uv ′
= (81)
dx v v2
donde:
d
u′ (γµDo ) = (γIDo )2 + 2γIDo ADo + A2Do
dγ
2
= 2γIDo + 2IDo ADo
d
v ′ (γµDo ) = 2(IDo − ADo ) (1 + γ)2
dγ
= 2(IDo − ADo ) · 2(1 + γ)
= 4(IDo − ADo )(1 + γ)
dg 2
(2γIDo + 2IDo ADo )v(γ) − u(γ)v ′ (γ)
= (82)
dγ v(γ)2
Donde:
7
4.2 Paso 2: Reescribir el denominador
El denominador es:
v(γ)2 = 4(IDo − ADo )2 (1 + γ)4 (83)
Que podemos expresar como:
2
dg (2γIDo + 2IDo ADo ) · 2(IDo − ADo )(1 + γ)2
=
dγ [2(IDo − ADo )(1 + γ)2 ]2
(γ IDo + 2γIDo ADo + A2Do ) · 4(IDo − ADo )(1 + γ)
2 2
−
[2(IDo − ADo )(1 + γ)2 ]2
Desarrollando el numerador:
2
(2γIDo + 2IDo ADo )(1 + γ) − 2(γ 2 IDo
2
+ 2γIDo ADo + A2Do )
2
= 2γIDo + 2γ 2 IDo
2
+ 2IDo ADo + 2γIDo ADo
− 2γ 2 IDo
2
− 4γIDo ADo − 2A2Do
2
= 2γIDo + 2IDo ADo + 2γIDo ADo − 4γIDo ADo − 2A2Do
2
= 2γIDo + 2IDo ADo − 2γIDo ADo − 2A2Do
8
4.7 Expresión final de la derivada
Combinando todos los términos:
2
dg 2γIDo + 2IDo ADo − 2γIDo ADo − 2A2Do
= (88)
dγ 2(IDo − ADo )(1 + γ)3
Factorizando el numerador:
dg IDo γ + IDo +A
2
Do
= 3
(90)
dγ (1 + γ)
df
γµDo (k+1) = γµDo (k) − α · (94)
dγµDo Do =γ Do (k)
γµ µ
7 Condición de parada
El algoritmo se detiene cuando:
9
(γµDo IDo + ADo )2
f (γ, η) = − NM ayr = 0 (98)
2ζη (IDo − ADo )(1 + γµDo )2
Do
(γIDo + ADo )2
= NM ayr (100)
2ζ(IDo − ADo )(1 + γ)2
Multiplicamos ambos lados por el denominador:
γ 2 IDo
2
+ 2γIDo ADo + A2Do = 2ζ(IDo − ADo )NM ayr (1 + 2γ + γ 2 )
= 2ζ(IDo − ADo )NM ayr
+ 4ζ(IDo − ADo )NM ayr γ
+ 2ζ(IDo − ADo )NM ayr γ 2 (102)
γ 2 IDo
2
+ 2γIDo ADo + A2Do − 2ζ(IDo − ADo )NM ayr (1 + 2γ + γ 2 ) = 0 (103)
Agrupamos términos semejantes:
γ 2 IDo
2
− 2ζ(IDo − ADo )NM ayr
+γ (2IDo ADo − 4ζ(IDo − ADo )NM ayr )
+ A2Do − 2ζ(IDo − ADo )NM ayr = 0
(104)
2
a = IDo − 2ζηDo (IDo − ADo )NM ayr (105)
b= 2IDo ADo − 4ζηDo (IDo − ADo )NM ayr (106)
c= A2Do − 2ζηDo (IDo − ADo )NM ayr (107)
r
b2 − 4ac
γ = −b ± (108)
2a
La solución para γ es:
q
(2IDo ADo − 4ζ(IDo − ADo )NM ayr )2 − 4 IDo
2
− (2IDo ADo − 4ζ(IDo − ADo )NM ayr ) ± − 2ζ(IDo − ADo )NM ayr A2Do − 2ζ(IDo − ADo )NM ayr
γ=
2 (I 2 − 2ζ(IDo − ADo )NM ayr )
(109)
Si η es igual a 1:
q
(2IDo ADo − 4(IDo − ADo )NM ayr )2 − 4 IDo
2
− (2IDo ADo − 4(IDo − ADo )NM ayr ) ± − 2(IDo − ADo )NM ayr A2Do − 2(IDo − ADo )NM ayr
γ= 2
2 IDo − 2(IDo − ADo )NM ayr
(110)
10
8 Normalización de las NM ayr
NMayr : [NMayr-Min , NMayr-Max ] (111)
La ecuación general de normalización es:
0 = aNMayr-Min + b (113)
1 = aNMayr-Max + b (114)
De la condición para el mı́nimo (113) se despeja b:
b = −aNMayr-Min (115)
NMayr NMayr-Min
Nscaled = − (120)
NMayr-Max − NMayr-Min NMayr-Max − NMayr-Min
NMayr − NMayr-Min
Nscaled = (121)
NMayr-Max − NMayr-Min
11