%%’(
!!!!!!!!物理化学学报!!
!!!!!!!!!!!!!!-./+01234!51$678$9#$!
D89B7:G
分子模拟研究气体在室温离子液体中的溶解度@
吴晓萍
摘要
刘志平汪文川
!北京化工大学化学工程学院分子与材料模拟研究室$教育部纳米材料重点实验室$北京%&&&*?
多种咪唑类离子液体中的溶解度$包括*种侧链长度的阳离子和’种阴离子#首先考察了*+>个离子对组成的体系中溶质分子插入次数对计算结果的影响$在此基础上计算了不同温度下气体在%
度较实验值偏高$精度与文献模拟结果相当;并预测了较高温度下CD*气体在I7.0.JIKL)J和I7.0.JIML>J中的溶解度#计算结果也正确地反映了+种气体在I7.0.JIML>J中溶解度实验值的相对大小#另外考察了常温下几种气体在不同室温离子液体中的溶解度$模拟结果表明气体在含有较长碳链和双
关键词:室温离子液体$分子动力学模拟$粒子插入法$气体溶解度中图分类号:D>)’#%
在化工生产过程中广泛使用的挥发性有机溶剂
是造成大气污染的主要原因$因此如何寻找新的“绿
$以减少挥发性有机溶剂的使用量$一直是色溶剂”
绿色化学的重要课题#室温离子液体即是近年来非
I%
常被看好的一类“绿色溶剂”#它们是在环境温度下呈液态的一类有机盐$由于其蒸汽压几乎检测不到$被认为是一种理想的挥发性有机溶剂的替代物#
化工分离过程是离子液体十分有前景的一个应用领域$如天然气的脱硫脱碳I)J、超临界CD*和离子液体构成的反应分离耦合体系I+
本文在作者先前建立的咪唑类离子液体的力
E*、D*、CH)场I%>J基础上$用分子模拟方法预测了FG、
及CD*等+种气体在%
;/3B102.4:Q/N12BGB54BR54/9:$I7.0.JIML>J
%模拟方法介绍
%&%分子力场及模拟细节
I7.0.JIKL)J及I7.0.JIML>J的全原子力场参数取自作者的前期工作I%>J#而%
对于溶质气体分子$氩气由T:
@
*&&+
国家自然科学基
!
刘志平等:分子模拟研究气体在室温离子液体中的溶解度
表!气体分子的力场参数D,[L6$\
==1?
!!^,3b7/=@2Y4A/0=2
b7Y4Y—4!0/=?2X0/0:2YX0/002
!!—!XX—XYXY—X
・Y_,7]6&)R$**‘+H%**‘GaH
:$:#OO:!!:#=OB::!!:#OB=0$:#10BO
:#1A=::#10:9:#09AA:#11=::#1:1::#09B9:#1:11
:#??1?:#0B=:#:A0
:#=:?:
:#0??1
:#=:?:
:#A:=O
:#=0=:
:#011?:#OO?1
:#==O0
式中&I为I
个粒子和其余%个粒子的相互作用能8〈6.N&U
Q-M
/0B2
的超额化学势!6.0可表示为
〈$6.N&U
!6.0TU&I#L+%
&0R
其中〈$6.N&U
温等压系综平均
#
释超额化学势!6.808而根据溶液热力学8气体的亨利常数(4可由下式计算:
#!6.80
(4&
分子动力学&’()模拟采用’(*+,’-.程序/012#模拟过程采用!
文献/$O2采用的体系仅包括$0@个离子对8这对于其研究的性质8如密度、内能和径向分布函数而言已足够8而对本文研究的气体溶解度8模拟体系都包括了0BO个离子对
粒子插入法是用于测定纯流体或混合物样本化学势的一种经典方法P通常称之为Q-M
考虑一含!个粒子的体系8体积为$
〈6.N&U
&=R
&1R
其中#$代表溶剂的数密度,时的气体溶解度!0
#结果与讨论
#
插入次数的确定
根据上述Q-M
可能穷尽体系构象空间的所有“空位”8即能量较低的位置8才能得到可靠的结果8因此足够的插入次数非常必要,8标准偏差为%,8
表#溶质分子的插入次数对计算结果的影响
D,[L60J;;6FV
c*3V6H!!X0…/[H-H2/I\A2!!X0…/[H-H2/I\A2YX0…/[H-H2/I\A2YX0…/[H-H2/I\A2(4:46+7*%3F
#‘C0?@1010?@101
?::::V-H63!‘GaH
U=
0:::::V-H63!‘GaH
(4‘’>,00:0B=O#O?#0
(4‘’>,0000O:9#:==#1
))6,
!!!!!!!!!!!!!!!
$%&’()
表-气体在室温离子液体#$%&%’#()*’和#$%&%’#+),’
中亨利常数计算结果
H=2&G9IG@JK$LM%@LN=@NL%[email protected]=@/
!!!!!!120.0317583=N/.OOGJG@NNG0PGJ=N
=IF>7=
QJ)*8)C*(9)C*B()()
R(*D(*BCD)*(D((8*996
;(
(BC9(996*(BC9(996*
((,(*)9(([email protected]
))B)(B)8(
D8C6)(B
)’6(’9D’D
)DB(,)(C,
8’8B’()C
:I6
:;(
验值1),!))3及文献模拟值1)6!)*3的比较’图中实线表示实验值#虚线表示计算值#点划线及点线表示文献模拟值’从图中可以看出#:;(气体在120.0314563中溶解度的计算值与实验值吻合非常好#本文的计算结果明显优于7
本文模拟结果显示:;(气体在120.0317583中的溶解度高于120.0314563#但:=/G@=等1(83的实验结果表明#:;(气体在120.0314563和120.0317583中的溶解度实际相差无几’由于两种离子液体具有同一阳离子#模拟计算得到的溶解度存在偏差则来自45E75E8、6二者力场的差异’
力场开发中电荷的分配存在一定任意性#在进行力场开发时#由于缺乏能量方面的实验数据进行验证而难以判断各种力场电荷取值的合理性1)83’本文对:;(气体溶解度的计算结果表明#先前力场中75E8各原子所分配的电荷值可能不尽合理#在今后的力场开发中#可以考虑增加气体的溶解度数据作为调节参数的依据’
最大偏差已经小于*+’在综合考虑了计算精度和速度的平衡之后#本文以后的计算中对每个溶剂构象都进行了(,,,,,次溶质分子的插入尝试’!
首先#用-./%0方法预测了几种气体在120.0314563和120.0317583中不同温度下的亨利常数#结果如表9所示’从表中可以看出#气体溶质的亨利常数基本上随着温度的增加而增大#也就是说气体在120.0314563和120.0317583中的溶解度基本上随温度升高而减小’从实验结果可知#该变化趋势与:;(气体在室温离子液体中溶解度的变化趋势是一致的1),!))3’表明仅通过分子间的物理作用可正确地反映:;(气体与120.0314563或120.0317583之间的相互作用’
!
图)和图(表示不同温度下:;(气体在120.0314563和120.0317583中的溶解度计算值与实
图
.不同温度下/0
!气体在#$%&%’#()*’中的溶解度)&1
!!!=&::4>4;994%?4>45
=
图!不同温度下/0!气体在#$%&%’#+),’中的溶解度)&1
!!!4;994%?4>45
!!!=
!
刘志平等:分子模拟研究气体在室温离子液体中的溶解度
$$9$
对比+,-.-/+@19/和+
上碳原子的数目会增加56’在离子液体中的溶解度3因此采用AB’!
值得注意的是4前文述及G表&?模拟表明456’
在+,-.-/+@19/的溶解度明显小于其在+,-.-/+012/的溶解度4而表9显示4其在+
图
#!$%&时气体在’()*)+’,-.+中溶解度计算值和
实验值的比较
-*/
6*)3294*1:9:5;*)5:494!$%&
H!IPQ.RS
#结论
H$I通过对特定体系溶质分子插入次数的比较4确定了合适的模拟条件:体系大小为’J2个离子对4对$%%%个溶剂分子构象进行插入时4每个构象的插入次数应不少于’%%%%%次3
=’I模拟结果正确反映了56’气体在+,-.-/+@19/和+,-.-/+012/中溶解度随温度的变化趋势4表明仅通过分子间的物理作用可以正确地反映56’气体与室温离子液体之间的相互作用3对于+,-.-/+@19/体系4本文模拟结果明显优于0!KLD等+$9/的结果4而对于+,-.-/+012/体系4本文结果和CDM.NN等+>J/的模拟结果非常相近4均高于实验值4而0!KLD等+>9/的模拟结果则较实验值偏低#
H&IO.K
H9I通过对常温下几种气体在不同离子液体+,-.-/+@19/、+,-.-/+AB’!/、+
!
图&为’()*时几种气体在+,-.-/+012/中溶解度计算值与实验值+$&/的比较3由于!’的溶解度超出检测范围4因而图中所示为实验值的上限#由图&可见4虽然本文模拟得到的气体溶解度普遍较实验值为高4但正确地反映了各气体在+,-.-/+012/中溶解度大小的相对顺序:即56’75897:;76’7!’3如前节所述4这种偏差的来源是01
表9所示为’()*下气体在几种室温离子液体中超额化学势及亨利常数的模拟结果3
前文=’#’#>节?已经述及模拟得到’()*下56’气体在+,-.-/+@19/中溶解度与实验结果符合非常好;在表9中456’气体在+,-.-/+AB’!/中亨利常数的模拟值为’32C0D4也接近其实验值&3EC0D+>&/3其在+
对比+,-.-/+AB’!/和+,-.-/+@19/456’在前者中的亨利常数为在后者中的9%F4亦即其在前者中的溶解度比后者要大’倍以上3其它气体也存在类似关系4即常温下气体更易溶于+,-.-/+AB’!/3
离子和具有较长侧链的咪唑型离子液体#
表?!$%&下气体在几种室温离子液体中超额化学势及亨利常数计算结果
AD,YU9Z.-LYDP.
]DR:;!’6’58956’
+,-.-/+@19/・!^T_-
!!!23$9!!!(3&E!!E3%%!!239)
!8^C0D$J2JE’’’%$E(232
+,-.-/+AB’!/・!^T_-
!!!J3’E!!)3’J!!J3E’!!939’
!8^C0DEJ’9))(J&’32
+
・!^T_-
!!!93E2!!!E39%!!!93(2!!!&3%E!
!8^C0D2E$(2E&&9$3J
+
!!93EE!!E39&!!J3%&!!&3$(
!8^C0D2&$)JE%&&$3E
))4(
!
)(R46K9\L
!!!!!!!!!!!!!!!
$%&’()
&-aF-/7b0%&+,F&127-EF&1=-%
&+
01C-(5OR1+2+7A+>?2117,222/)L6)6O313#X’S’;
)K>-F#A’B’;YF1
51;/)(L9\
B1/F1#Z’Z’Y’E-H&’(-=&5,-I?2117,51;/)9>%;+7#X’P’:’;
)K\LR
)\0-,3+&7#Z’;V-N[+2#Y’H&’(-=&5,-=&5,-H&’(-?5*7*,5:/K(9Z&+N1’;
2111,*;/4R
(50-81(5OR1+2+7A+>?5***,5:;/(L
01c‘+&&#T’O’;J-21/%*-.+7#X’E-F,-=&5,-()X%2C+’;
G8
((T+7=2-C
5**
(R>8FI12=7+.#Z’B’;>11[7%
6K6
(4V-/%E#U’X’=&5,-H&’(-?5*
(6Z&&+
(!(!(45
(K:1/+’;O313#X’S’;0%22%‘#J’?’;U2+
X’Q’;01C-
O+//%#>2-?()*+,-./KLRO+//%#1A5B>#5A+>2(?011-,2/RKR
VF#0’;S%F#W’;0-#>2’(+(C8D>’?2112,A31%#T’U’;34/)69
018=%
=&5,-G8
U&1’Z’;Y1#1A5?5***89**/(\
]F#A’;U2+’Z’;51:/(4R9
Y%&I2+8#X’T’;O+//%/HA8D1
V+&=%
S+-E#V’F/L5M-=&5,-?N/#-OD-?2111,9*/V177+27,3+-/#B’;R99(
JFE1#T’;^-1#X’;01F2+2#]’E-=&5,-O/L-P>#>?)5S1E;7#Z’B’;
2119,7;/94K
))YF77%#>?
2119,7;/4\5
)(U1&=F7#*’D’;:F&I+2=7%F%#Y’;T+B1%&-#T’V’
E-H&’(-=&5,-I?2117,51;/9()
01C-’;
2112,51
)4T+7,31E;7#X’;B
!
V@#^-1%!B-
>?@#A3-!B-
VZP]#V+
!P+K+(+8/8UB825
=>’?$@&?J3+7%&FI-&-=-+7%__-.+C17+7#:d(#P(#d(#Z21
I8V-/%E=+7=!;12=-,&+-
A
]177%&FI-&-=8
M
*+,+-.+/:012,34#(556;*+.-7+/:018(9#(556’J3+B2%N+,=OF;;%2=+/I8POQ:!(5(RK5)5
:%22+7;%?@#A3-!B-
%%’(
!!!!!!!!物理化学学报!!
!!!!!!!!!!!!!!-./+01234!51$678$9#$!
D89B7:G
分子模拟研究气体在室温离子液体中的溶解度@
吴晓萍
摘要
刘志平汪文川
!北京化工大学化学工程学院分子与材料模拟研究室$教育部纳米材料重点实验室$北京%&&&*?
多种咪唑类离子液体中的溶解度$包括*种侧链长度的阳离子和’种阴离子#首先考察了*+>个离子对组成的体系中溶质分子插入次数对计算结果的影响$在此基础上计算了不同温度下气体在%
度较实验值偏高$精度与文献模拟结果相当;并预测了较高温度下CD*气体在I7.0.JIKL)J和I7.0.JIML>J中的溶解度#计算结果也正确地反映了+种气体在I7.0.JIML>J中溶解度实验值的相对大小#另外考察了常温下几种气体在不同室温离子液体中的溶解度$模拟结果表明气体在含有较长碳链和双
关键词:室温离子液体$分子动力学模拟$粒子插入法$气体溶解度中图分类号:D>)’#%
在化工生产过程中广泛使用的挥发性有机溶剂
是造成大气污染的主要原因$因此如何寻找新的“绿
$以减少挥发性有机溶剂的使用量$一直是色溶剂”
绿色化学的重要课题#室温离子液体即是近年来非
I%
常被看好的一类“绿色溶剂”#它们是在环境温度下呈液态的一类有机盐$由于其蒸汽压几乎检测不到$被认为是一种理想的挥发性有机溶剂的替代物#
化工分离过程是离子液体十分有前景的一个应用领域$如天然气的脱硫脱碳I)J、超临界CD*和离子液体构成的反应分离耦合体系I+
本文在作者先前建立的咪唑类离子液体的力
E*、D*、CH)场I%>J基础上$用分子模拟方法预测了FG、
及CD*等+种气体在%
;/3B102.4:Q/N12BGB54BR54/9:$I7.0.JIML>J
%模拟方法介绍
%&%分子力场及模拟细节
I7.0.JIKL)J及I7.0.JIML>J的全原子力场参数取自作者的前期工作I%>J#而%
对于溶质气体分子$氩气由T:
@
*&&+
国家自然科学基
!
刘志平等:分子模拟研究气体在室温离子液体中的溶解度
表!气体分子的力场参数D,[L6$\
==1?
!!^,3b7/=@2Y4A/0=2
b7Y4Y—4!0/=?2X0/0:2YX0/002
!!—!XX—XYXY—X
・Y_,7]6&)R$**‘+H%**‘GaH
:$:#OO:!!:#=OB::!!:#OB=0$:#10BO
:#1A=::#10:9:#09AA:#11=::#1:1::#09B9:#1:11
:#??1?:#0B=:#:A0
:#=:?:
:#0??1
:#=:?:
:#A:=O
:#=0=:
:#011?:#OO?1
:#==O0
式中&I为I
个粒子和其余%个粒子的相互作用能8〈6.N&U
Q-M
/0B2
的超额化学势!6.0可表示为
〈$6.N&U
!6.0TU&I#L+%
&0R
其中〈$6.N&U
温等压系综平均
#
释超额化学势!6.808而根据溶液热力学8气体的亨利常数(4可由下式计算:
#!6.80
(4&
分子动力学&’()模拟采用’(*+,’-.程序/012#模拟过程采用!
文献/$O2采用的体系仅包括$0@个离子对8这对于其研究的性质8如密度、内能和径向分布函数而言已足够8而对本文研究的气体溶解度8模拟体系都包括了0BO个离子对
粒子插入法是用于测定纯流体或混合物样本化学势的一种经典方法P通常称之为Q-M
考虑一含!个粒子的体系8体积为$
〈6.N&U
&=R
&1R
其中#$代表溶剂的数密度,时的气体溶解度!0
#结果与讨论
#
插入次数的确定
根据上述Q-M
可能穷尽体系构象空间的所有“空位”8即能量较低的位置8才能得到可靠的结果8因此足够的插入次数非常必要,8标准偏差为%,8
表#溶质分子的插入次数对计算结果的影响
D,[L60J;;6FV
c*3V6H!!X0…/[H-H2/I\A2!!X0…/[H-H2/I\A2YX0…/[H-H2/I\A2YX0…/[H-H2/I\A2(4:46+7*%3F
#‘C0?@1010?@101
?::::V-H63!‘GaH
U=
0:::::V-H63!‘GaH
(4‘’>,00:0B=O#O?#0
(4‘’>,0000O:9#:==#1
))6,
!!!!!!!!!!!!!!!
$%&’()
表-气体在室温离子液体#$%&%’#()*’和#$%&%’#+),’
中亨利常数计算结果
H=2&G9IG@JK$LM%@LN=@NL%[email protected]=@/
!!!!!!120.0317583=N/.OOGJG@NNG0PGJ=N
=IF>7=
QJ)*8)C*(9)C*B()()
R(*D(*BCD)*(D((8*996
;(
(BC9(996*(BC9(996*
((,(*)9(([email protected]
))B)(B)8(
D8C6)(B
)’6(’9D’D
)DB(,)(C,
8’8B’()C
:I6
:;(
验值1),!))3及文献模拟值1)6!)*3的比较’图中实线表示实验值#虚线表示计算值#点划线及点线表示文献模拟值’从图中可以看出#:;(气体在120.0314563中溶解度的计算值与实验值吻合非常好#本文的计算结果明显优于7
本文模拟结果显示:;(气体在120.0317583中的溶解度高于120.0314563#但:=/G@=等1(83的实验结果表明#:;(气体在120.0314563和120.0317583中的溶解度实际相差无几’由于两种离子液体具有同一阳离子#模拟计算得到的溶解度存在偏差则来自45E75E8、6二者力场的差异’
力场开发中电荷的分配存在一定任意性#在进行力场开发时#由于缺乏能量方面的实验数据进行验证而难以判断各种力场电荷取值的合理性1)83’本文对:;(气体溶解度的计算结果表明#先前力场中75E8各原子所分配的电荷值可能不尽合理#在今后的力场开发中#可以考虑增加气体的溶解度数据作为调节参数的依据’
最大偏差已经小于*+’在综合考虑了计算精度和速度的平衡之后#本文以后的计算中对每个溶剂构象都进行了(,,,,,次溶质分子的插入尝试’!
首先#用-./%0方法预测了几种气体在120.0314563和120.0317583中不同温度下的亨利常数#结果如表9所示’从表中可以看出#气体溶质的亨利常数基本上随着温度的增加而增大#也就是说气体在120.0314563和120.0317583中的溶解度基本上随温度升高而减小’从实验结果可知#该变化趋势与:;(气体在室温离子液体中溶解度的变化趋势是一致的1),!))3’表明仅通过分子间的物理作用可正确地反映:;(气体与120.0314563或120.0317583之间的相互作用’
!
图)和图(表示不同温度下:;(气体在120.0314563和120.0317583中的溶解度计算值与实
图
.不同温度下/0
!气体在#$%&%’#()*’中的溶解度)&1
!!!=&::4>4;994%?4>45
=
图!不同温度下/0!气体在#$%&%’#+),’中的溶解度)&1
!!!4;994%?4>45
!!!=
!
刘志平等:分子模拟研究气体在室温离子液体中的溶解度
$$9$
对比+,-.-/+@19/和+
上碳原子的数目会增加56’在离子液体中的溶解度3因此采用AB’!
值得注意的是4前文述及G表&?模拟表明456’
在+,-.-/+@19/的溶解度明显小于其在+,-.-/+012/的溶解度4而表9显示4其在+
图
#!$%&时气体在’()*)+’,-.+中溶解度计算值和
实验值的比较
-*/
6*)3294*1:9:5;*)5:494!$%&
H!IPQ.RS
#结论
H$I通过对特定体系溶质分子插入次数的比较4确定了合适的模拟条件:体系大小为’J2个离子对4对$%%%个溶剂分子构象进行插入时4每个构象的插入次数应不少于’%%%%%次3
=’I模拟结果正确反映了56’气体在+,-.-/+@19/和+,-.-/+012/中溶解度随温度的变化趋势4表明仅通过分子间的物理作用可以正确地反映56’气体与室温离子液体之间的相互作用3对于+,-.-/+@19/体系4本文模拟结果明显优于0!KLD等+$9/的结果4而对于+,-.-/+012/体系4本文结果和CDM.NN等+>J/的模拟结果非常相近4均高于实验值4而0!KLD等+>9/的模拟结果则较实验值偏低#
H&IO.K
H9I通过对常温下几种气体在不同离子液体+,-.-/+@19/、+,-.-/+AB’!/、+
!
图&为’()*时几种气体在+,-.-/+012/中溶解度计算值与实验值+$&/的比较3由于!’的溶解度超出检测范围4因而图中所示为实验值的上限#由图&可见4虽然本文模拟得到的气体溶解度普遍较实验值为高4但正确地反映了各气体在+,-.-/+012/中溶解度大小的相对顺序:即56’75897:;76’7!’3如前节所述4这种偏差的来源是01
表9所示为’()*下气体在几种室温离子液体中超额化学势及亨利常数的模拟结果3
前文=’#’#>节?已经述及模拟得到’()*下56’气体在+,-.-/+@19/中溶解度与实验结果符合非常好;在表9中456’气体在+,-.-/+AB’!/中亨利常数的模拟值为’32C0D4也接近其实验值&3EC0D+>&/3其在+
对比+,-.-/+AB’!/和+,-.-/+@19/456’在前者中的亨利常数为在后者中的9%F4亦即其在前者中的溶解度比后者要大’倍以上3其它气体也存在类似关系4即常温下气体更易溶于+,-.-/+AB’!/3
离子和具有较长侧链的咪唑型离子液体#
表?!$%&下气体在几种室温离子液体中超额化学势及亨利常数计算结果
AD,YU9Z.-LYDP.
]DR:;!’6’58956’
+,-.-/+@19/・!^T_-
!!!23$9!!!(3&E!!E3%%!!239)
!8^C0D$J2JE’’’%$E(232
+,-.-/+AB’!/・!^T_-
!!!J3’E!!)3’J!!J3E’!!939’
!8^C0DEJ’9))(J&’32
+
・!^T_-
!!!93E2!!!E39%!!!93(2!!!&3%E!
!8^C0D2E$(2E&&9$3J
+
!!93EE!!E39&!!J3%&!!&3$(
!8^C0D2&$)JE%&&$3E
))4(
!
)(R46K9\L
!!!!!!!!!!!!!!!
$%&’()
&-aF-/7b0%&+,F&127-EF&1=-%
&+
01C-(5OR1+2+7A+>?2117,222/)L6)6O313#X’S’;
)K>-F#A’B’;YF1
51;/)(L9\
B1/F1#Z’Z’Y’E-H&’(-=&5,-I?2117,51;/)9>%;+7#X’P’:’;
)K\LR
)\0-,3+&7#Z’;V-N[+2#Y’H&’(-=&5,-=&5,-H&’(-?5*7*,5:/K(9Z&+N1’;
2111,*;/4R
(50-81(5OR1+2+7A+>?5***,5:;/(L
01c‘+&&#T’O’;J-21/%*-.+7#X’E-F,-=&5,-()X%2C+’;
G8
((T+7=2-C
5**
(R>8FI12=7+.#Z’B’;>11[7%
6K6
(4V-/%E#U’X’=&5,-H&’(-?5*
(6Z&&+
(!(!(45
(K:1/+’;O313#X’S’;0%22%‘#J’?’;U2+
X’Q’;01C-
O+//%#>2-?()*+,-./KLRO+//%#1A5B>#5A+>2(?011-,2/RKR
VF#0’;S%F#W’;0-#>2’(+(C8D>’?2112,A31%#T’U’;34/)69
018=%
=&5,-G8
U&1’Z’;Y1#1A5?5***89**/(\
]F#A’;U2+’Z’;51:/(4R9
Y%&I2+8#X’T’;O+//%/HA8D1
V+&=%
S+-E#V’F/L5M-=&5,-?N/#-OD-?2111,9*/V177+27,3+-/#B’;R99(
JFE1#T’;^-1#X’;01F2+2#]’E-=&5,-O/L-P>#>?)5S1E;7#Z’B’;
2119,7;/94K
))YF77%#>?
2119,7;/4\5
)(U1&=F7#*’D’;:F&I+2=7%F%#Y’;T+B1%&-#T’V’
E-H&’(-=&5,-I?2117,51;/9()
01C-’;
2112,51
)4T+7,31E;7#X’;B
!
V@#^-1%!B-
>?@#A3-!B-
VZP]#V+
!P+K+(+8/8UB825
=>’?$@&?J3+7%&FI-&-=-+7%__-.+C17+7#:d(#P(#d(#Z21
I8V-/%E=+7=!;12=-,&+-
A
]177%&FI-&-=8
M
*+,+-.+/:012,34#(556;*+.-7+/:018(9#(556’J3+B2%N+,=OF;;%2=+/I8POQ:!(5(RK5)5
:%22+7;%?@#A3-!B-