堆积的胸腺嘧啶体系光物理失活的动力学模拟.pdf

No ve mbe r 物理化学学 W u l i Hu a x u e X u e b a o 1 A ct a 尸 一 C h im S in 2 0 1 1 2 7 1 1 2 5 5 9 2 5 6 4 2 5 5 9 Ar t icl e 豆育升 www wh x b p k u e d u cn 堆积的胸腺嘧啶体系光物理失活的动力学模拟 李伟 袁帅 张文英 李安阳 舒坤贤 重庆邮电大学高性能计算与应用研究所 重庆 4 0 0 0 6 5 唐红 摘要 采用半经典动力学方法模拟了丌 堆积的胸腺嘧啶体系最低激发态的光物理失活过程 设置激光脉冲仅 作用于一个胸腺嘧啶分子 另一胸腺嘧啶分子T 保持基态 模拟发现由于T与T 之间存在丌 士 I 积相互作用 导 致电荷转移 形成T带负电荷 T 带正电荷的激基复合物 由于相邻分子的空间效应阻碍了激发的T分子到达圆 锥相交所必需的强烈扭曲 激基复合物的寿命比单体增长 当分子间距离缩短至0 3 n m后 T分子C 一C 键扭 曲程度最大 此时发生电荷重组 两个胸腺嘧啶分子均恢复电中性 电荷重组诱导T r 分子发生畸变 并在C 5 一 C 扭曲最大时避免相交 体系衰减至基态 T和T 分子均恢复平面构型 关键词 胸腺嘧啶 半经典动力学模拟 电荷转移 激基复合物 失活 中图分类号 O6 4 1 06 4 4 Dy n amics Simul a t ion o f Ph o t op h y si C a I De a ct i v a t i on Pa t h wa y f or St a c k e d Th ym i n e s DOU Yu Sh e n g L I W ei YUAN Sh u a i ZHANG W e n Yi n g L I An Ya n g SHU Ku n Xi a n T ANG Ho n g I n s t it u t e o f Hi g h P e rf o r ma n ce C o m p u t i n g a n d A p p l ica t i o n R e s e a r ch C e n t e r o fN e t wo r k a n d C o m p u t i n C h o n g q i n g U n i v e r s ity o f P o s t s a n d T e l e co mmu n ica t io n s C h o n g q i n g 4 0 0 0 6 5 P R C h i n a Ab s t r a ct A s e mi cl a s s i cal dy n a mi cs s i mu l a t i o n s t u d y wa s u n d e r t a k e n t o d e t e r mi n e t h e p h o t op h y s i c a l d e a ct i v a t i o n o f t h e I o we s t e x ci t e d s t a t e o f t wo s t a c k e d t h y mi n e s On ly o n e t h y mi n e r e f e r r ed t o a s wa s e x c i t e d b y a la s e r p u l s e an d t h e o t h e r mole c u l e r e f e r r e d t o a s T r ema i n e d j n t h e gr o u n d s t a t e Th e s i mu l a t i on r e s u lt s s h o w t h a t ch a r g e t r an s f e r b e t we e n t h e t wo t h y mi n e s b e ca u s e o f a E s t a c k i n g i n t e r a ct i o n l e a d s t o t h e f o r ma t i o n o f an e x c i me r s t a t e wh i ch i n c l u d e s a n e g a t i v e T a n d a p o s i t i v e T Ad di t i o n a l l y l h e s i mu l a t i o n s t u d y i n d i ca t e s t h a t a s t e r i c e f f e ct o f t h e n e i g h b o r i n g b a s e s i n h i b i t s t h e o u t o f p lan e d e f o r ma t i o n wh i c h i S e s s e n t i a I i n a c c e s s i n g t h e c o n i caI i n t e r s e ct i o n b e t we e n t h e l o we s t e l e ct r o n i c e x c i t e d s t a t e a n d t h e g r o u n d s t a t e T h e s t er i c e ff e ct e v e n t u a l l y l e a d s t o a I on g e r ele c t r on i c e x c i t e d s t a t e I i f e t i me f o r t h e t wo s t a ck ed t h y mi n e s Th e s i mu l a t i o n r e s u lt s r e v e al t h a t wh e n t h e in t e r b a s e d i s t an ce i S I e s s t h a n 0 3 n m t h e mo l e c u l e T h a s a r e ma r k a b l e d e f o r ma t i o n a t i t s C5 a n d C6 s i t e s r e s u l t i n g i n c h a r g e r e co mb i n a t i o n Th e c h a r g e r e comb i n a t i on u lt i ma t e l y ma k e s t h e s y s t e m elect r o n i ca l l y n eu t r a I On t h e o t h e r h a n d t h e mo l e c u le T h a s a s t r on g tw i s t ab ou t i t s C5 一 C6 b o n d i n t h e pr o x i mi t y o f t h e a v o i d e d cr o s s in g b y wh i c h t h e s y s t e m d e c a y s t o t h e gr ou n d s t a t e Fi n a l ly t h e tw o t h y mi n e mole c u l e s i n t h e i r gr o u n d s t a t e s r e co v e r t h ei r p lan ar g e o me t r i e s Ke y Wor ds Th y mi n e Se mi c l a s s i cal d y n a mi c s i mu la t i o n Ch a r g e t r a n s f e r Ex c i mer De a c t i v a t i o n Re ce iv e d Ma y2 4 20 1 1 Re vis e d J u l y2 6 2 01 l P u b l is h e do nWe b S e p t e mbe r 5 2 0 1 1 C o r r e s p o n d i n g a u t h o r E ma i l d o u y s cq u p t e d u cn T h e p r o j e ct wa s s u p p o r t e d b y t h e Na t i o n a l Na t u r a l S ci e n ce F o u n d a t i o n o f C h i n a 2 1 0 7 3 2 4 2 Na mr a l S cie n ce F o u n d a t i o n o f C h o n g q i n g C h i n a cs t c2 0 1 1 j j A0 0 0 0 9 a n dp r o j e ct o f t h e S ci e n ce T e ch n o l o g y F o u n d a t i o n o f C h o n g q i n gE d u ca t i o nC o mmi t t e e C h i n a KJ 0 0 5 0 7 国家 自然科学基金 2 1 0 7 3 2 4 2 重庆市 自然科学基金 cs t c2 0 1 l j j A0 0 0 0 9 重庆市教委科学技术项 目 KJ 1 0 0 5 0 7 资助 E d i t o r i a l o f f ice o f A ct a P h y s i co C h i mi ca S i n i ca 25 6 0 Act a P h y S 一 C h i m S i n 2 0 1 1 0 I 2 7 1 引 言 脱氧核糖核酸 D N A 在紫外光 2 0 0 3 0 0 n m 辐 射下容易发生光化反应而造成损伤 其 中一种光化 产物是在 D NA同一条链上相邻的嘧啶碱基 问形成 环丁烷 型嘧啶二聚体 C P D 自从半世纪 以前发现 C P D后 很多实验 及理论 研究都致力于深入 理解 C P D形成的机理过程 胸腺嘧啶二聚是一个超 快光化反应 S ch r e ie F 观察到 d T 在 2 6 6 n m紫外 光辐射后 1 p s 内形成ci s s y n 型环丁烷胸腺 嘧啶二 聚体 T T 表 明该 反应遵循 无能垒 的非绝 热机 理 D NA中 T T 的产率很低 仅 2 一 3 但 晶 体研究表明 1 6 1 7 当两个胸腺嘧啶分子处于其二聚体 构型的位置时 二聚反应的量子产率几乎为 1 因而 S ch r e ie r 等 提 出在 D NA中观察 到二聚反应 的量子 产率低是因为D NA螺旋链中连续碱基对问平均 3 6 的扭转角太大而使反应难以发生 分子动力学 MD1 模拟 探讨 了可二聚的胸腺嘧啶构型 以C 一C C 一C 之间的平均距离 以及 C 一C 一C 一C 扭 转角 大小确定胸腺嘧啶单体二聚反应的可能性 提 出最容易二聚的反应构型为 d 0 3 4 n m 2 7 可二聚 的区间在 0 3 n m d 0 4 5 n m 及0 q 4 7 电子 激 发 单 重 态 是 二 聚 体 形 成 的 主 要 通 道 B o g g io P a s q u a等 用 完 全 活 性 空 间 自洽 场 C AS S C F 方法确定并描述 了圆锥相交 C I 的特 点 推定从最初单重激发态经 C I 点能超快地形成二 聚 体 或发生光物理失活 核酸多聚体 中由于存在 堆积效应 其激发态动力学特征与单一碱基明显不 同 最 显著 的特征 是存在长寿命 的 激基 复合体 e x ci me r e x cip l e x 对 e x cime r 长寿命激 发态 的一种 解释 是 由于激发态 能量对结构变化不敏感 堆积 碱 基体系相对 温和 的振动 导致 了较平缓的激发 态势能面 延滞 了相交点 的出现 使得激发态寿命 延长 由于嘌呤分子的7 g 体系比嘧啶分子更大 因此 嘌呤分子之 间的堆积 力更大 嘌呤 e x ci me r 的寿命 更长 可达 1 0 1 0 0 p s 相对而言 嘧啶 e x cime r 的 寿命较短 仅比单体激发态寿命略长 由于胸腺嘧啶二聚体量子产率很低 表明堆积 胸腺嘧啶体系的激发态 以光物理方式失活为主 受 限于过于简单的反应坐标 量子化学计算难 以顾及 分子整体的运动 无法得到光物理失活 的详细动力 学信息 本文采用半经典动力学方法模拟 了堆积 的 胸腺嘧啶体系最低激发态的光物理失活过程 可以 得到每一个反应 自由度 以及能量随时问的变化 模 拟结果可 以提供堆积 的碱基体系 的长寿命激发态 的形成和失活 的详细动力学信息 为理解这一反应 机理提供一定的理论支持 2 模型与方法 本 文 采 用 半 经 典 电 子 辐 射 离 子 动 力 学 S E R I D 模型模拟激光诱导碱基光化学反应的动力 学计算模型 如所谓半经典方法是指对 电子运动由 含 时S ch r d i n g e r 方程确定 而辐射场和核运动则由 经典力学处理 这种模拟方法考虑激光场对 电子结 构 的作 用 并 且采 用 非 绝 热动 力 学 的近似 方 法 E h r e n f e s t 方法 和紧束缚近似来建立模型 在 该方法 中 电子运动 由含时 薛定谔方程描 述 i f S 日 1 代表原子轨道 的重叠矩阵 激光脉冲辐射场的矢 势 可以通过含 时的P e ie r l s 替代 与电子耦合 X X e x P I 2 式中月 和 Ha X X 是基 函数 a 和 b 在 P e ie r l s 替代前 后的哈密顿矩阵元 代表核坐标 q e 是 电 子电荷 哈密顿矩 阵元 重叠矩阵元采用紧束缚近似 密度泛 D F T B 方法计算 这一近似不需要引入另 外的参量 并且原则上适用于任意强度 的含时外场 矢势 是一个 向量函数 全面地描述激光脉冲 其 中 包括脉冲 时宽 光子频率 脉冲强度和光子极化等 原子核运动的力由广义的E h r e n f e s t 方程来描述 一 吉 吉 蕞 暑 一 a U r a 3 式 中 是 离 子一 离子 排 斥势 O HI O X O S O X 和 a a 是 和 相对于每一个原子位置 的梯 度 所有 的梯度矩 阵元可经过数值求解法确定 原 子位置 由经典的牛顿方程描述 F d t 式中 是原子的原子量 牛顿方程 由V e l o cit y V e 卜 l e t 方法求解 DF T B方法 已经成功地应用于研究多种体系 的基态性质 最近发现 DF T B方法作为 D F T响应方 程 的简化计算 其对激发态能量和结构的计算结果 可与 a b in it io 和含时密度泛函理论 T DD F T 计算的 结果媲美 在某些氨基酸模型体系中 D F T B方法 NO 1 1 豆育升等 堆积的胸腺嘧啶体系光物理失活的动力学模拟 2 5 6 l 可 以描述极化和 电荷转移 效应 对于本文所研 究 的较大的生物分子体系 采用 D F T B方法计算 的电 荷转移激发态可以满足定性描述的需要 近年来该 模型成功地模拟 了几个生物分子光化学反应 包括 偶氮苯及其衍生物的异构化 腺嘌呤无辐射超快 失活 胸腺嘧啶光致二聚反应 及其二聚体 的光致 解离 反应 并描述了堆积的腺嘌呤体系超长寿命 的e x cime r 和b o n d e d e x cime r 3 结果与讨论 3 1 优化基态构型 为 了研究堆积效应对碱基激 发态动力学 的影 响 我 们 从 其 B型 结 构 的 d T 中随 机截 取 两 个 相 邻 胸腺 嘧 啶分 子片 段 忽 略核 糖 和磷 酸 根 由 于考虑 堆积效应 的影 响 在 保持重 叠程度不 变 的 基 础 上 平 移 使分 子 间距离 缩 短 调 整后 的 结构 参数 为D C 一C 5 0 3 6 0 n m D C C 6 0 3 6 4 a m C C C C 2 7 4 见 1 所不 在室温 不引入激光 的条件下 经过 5 0 0 f s 的动 力学模 拟 得 到堆积胸 腺嘧 啶体系 的基 态平 衡构 型 将平衡构型再进行 5 0 0 f s 的模拟 间隔 5 0 f s 取 得 1 0 个构型作为初始构型 初始构型包括 了原子坐 标和速度 加速度等动力学信息 3 2 模拟结果与讨论 根据 e x cime r 的定义 本文 的模拟中激光只作 用于一个胸腺嘧啶分子 选择频率为4 0 e V 最大半 波 宽f F WH M 为 2 5 f s 的激光脉冲照射位 于上方 的 胸腺嘧啶分子 T T分子被激发 位于下方的腺嘌呤 分 子 T 未 受 激 光 照 射 而 保 持 基 态 脉 冲 频 率 与 S E R I D方法优化得到 的基态平衡状态 下最高 占据 分 子轨道 H OMO 与最低未 占据 分子轨道 L U MO 的能差相 吻合 我们对 1 0 个 不 同初始构型进 行共 4 0 0 余次模拟 发现7 g 堆积 的胸腺嘧啶体系除了形成 o c5 C 6 C C 2 7 4 图1 堆积的胸腺嘧啶体系结构 a 与原子编号及参数 b F ig 1 Co n f o r ma t i o n a a n d a t o mic l a b e l in g a n d p a r a me t e r s b f o r t h e s t a ck e d t h y min e s s y s t e m C P D 还存在形成 e x cime r 的光物理失活通道 形成 C P D的研 究 已在他处 发表 本文着重 讨论 e x cime r 的 光 物 理 失活 通 道 此 处选 择 辐 射 通 量 密度 为 6 6 8 5 J 1 T I 的一条轨迹作为代表进行报道 图2描述体系轨道能量和 电子布居随时间的变 化 从 图2 a 看 出 初始 时刻 L UMO能量 为一 2 0 e V H0M0能量 为一 6 5 e V 激光作用 以后 L UM0能量 显著下降 1 0 0 f s 后降为平均一 3 5 e V 并在此值上下 波动 同时 H OMO能量 略上升 在一 6 0 e V上下波 动 7 0 0 f s 后 HO MO能量逐 渐上升 并在 7 6 7 f s 与 L U MO发生能级简并 即避免相交 此后两个轨道 能量恢复至其初 始值 图2 b 显示 了电子跃 迁 的细 节 激光 作用 以后 约 1 2 个 电子从 HO MO激 发至 L UMO 避免相交 即7 4 0 f s 时发生电子从L U MO至 HO MO的跃迁 分 子恢 复 到基 态 激 发态 寿命 为 69 0 f s 图 3 描述 了C 一C 和 C 一C 6 键长 随时间的变 化 初始时刻 C 一C 和 C 5 一C 6 键长为 0 1 3 8 n m 是 典型 的共轭体系 C C双键键长 激光作用 以后 由 于 T分子被激 发后 键被破坏 C 一C 转变为 键 其键 长延长至 0 1 5 4 n m 而 T 分子未被激发 C 一 C 6 键长保持不变 7 0 f S 以后 c 一c 键长稍缩短 变 为平 均 0 1 4 8 a m C 一C 键 长 延长 至 平 均 0 1 4 2 n m 6 6 0 7 4 0 f s 时 C 一C 6 键长从 0 1 4 2 n m延长至 图2 两种堆积的胸腺嘧啶体系的HOMO与L UMO的能 量 a 和电子布居 b 随时间的变化图 F i g 2 Or b i t a l e n e r g y o f t h e HOMO a n d L U MO a a n d e l e ct r o n p o p u l a t io n b o f t wo s t a ck e d t h y mi n e s v a r ia t io n wit h t im e HOMO t he h ig h e s t o ccu p ie d mo l e cul a r o r b it a l LUM 0 t h e l o we s t u n o ccu pie d mo l e cu l a r o r b it a l NO 1 1 豆育升等 堆积 的胸腺嘧啶体系光物理失活 的动力学模拟 2 5 6 3 激发态分子 的氧化性和还 原性都 比其相应 的基 态 分子有显著增加 根据 R e h m We l l e r 方程 A G 9 6 4 8 E D D E A A e 1 一 o 巴 可判 断激发 的胸腺 嘧啶分子和基态分子之 问电子 转移的 自由能变化 式中 G o 是 电子给体 D和 电子 受体 A之 间发生光诱导 电子转移 反应 生成溶剂 分 离的 自由基离子对时的标准 自由能变 E D D 和 A 1 分别为D的氧化电位和 A的还原 电位 为受激 分子 的激发能 E o 为生成溶剂分离 的 自由基离 子 对所需 的静 电作用 能 在 极性溶剂 中其值很小 可 忽略 胸腺嘧啶的氧化 电位 和还原 电位分别为2 1 1 和 一 2 1 8 e V 其最 低 激发 态 的垂 直激 发 能 为 4 8 e V 由此判 断 胸腺嘧啶 e x cime r 间的 电子转移在 热力学上是可行的 图6 是分子间距离与T分子环上的净 电荷关系 图 此处的分子间距离由C 一C 和 C 一C 距离的平 均值表示 原子 电荷 由电子波 函数 向原子轨道 的投 影计算得到 分子 问距离的初始值为 0 3 6 n m 在该 距 离下 7 C 堆积力较强 T分子受激 以后 约 1 2 个 电 子 从 T 分 子 向 T转 移 以填 充 HOMO的 空 穴 7 0 6 0 0 f s 之 间 T分子净 电荷保持一1 2 e 这可 以解 释 该时问段 C 一C 键延长和 C 一C 键缩短 的现象 当电子转移到 T分子后 C 一C 之问 电子密度增大 使得 C 一C 键级增加 键长缩短 反之 电子从 T 分 子迁移后 C 一C 之间 电子密度减小 使得键级 降 低 键 长增长 6 0 0 f s以后 分 子 问距 离缩短 为 0 3 n m 对 比图4 a 发现 此时 C 一C 键扭 曲程度 最大 0 N C C C 1 达 到一 2 8 6 强烈 的分子 问相互作 用可 能改变 分子的氧 化 还 原电位而 导致 电荷重组 电 E C 8 亡 鲤 图6 分子间距离与T分子净电荷的关系 Fig 6 Re l at ion be t we e n int e rba s e dis t a nce a nd ne t ch a r g e o nT 子 从 T 回迁 到 两个 胸 腺 嘧 啶恢 复 电 中性 在 电荷 转 移 阶 段 即 e x cime r 阶段 带 正 电荷 的 T 一 直保 持 平面 构 型 电荷 重 组后 T 分 子发 生 明 显 畸变 其 C 一C d 键 在 7 4 0 f s 扭 曲程 度 最大 此 时0 N C C C 为一l 6 2 比6 0 0 f s 时减小 正是 由于 C 一C 的强烈扭 曲导致体系发生非绝热衰减而返 回基态 本研 究模拟 的多条e x cime r 光物理的失活轨迹 激发态寿命为 6 0 0 1 2 0 0 f s 均 比相 同方法模拟的单 体胸腺嘧啶激发态寿命长 其失活 的方式也十分类 似 均在 T 分子畸变程度最大时发生无辐射跃迁衰 减至基态 所有光 物理 的失活通道 中的e x cime r 均 在避免相交点以前都已经通过 电荷重组恢复中性 与堆积 的胸腺嘧啶体系相 比 腺嘌呤分子的兀 共轭体系体系较大 因而它们 的堆积相互作用较强 这种 更强 的相 互作用 导致 腺嘌 呤腺 嘌 呤 AA e x cime r 稳 定性更高 其寿命远大于胸腺嘧啶一 胸腺嘧 啶 T T e x cime r 而且 A A e x ci me r 之 间的相互作用 也更复杂 失活 的通道也更 多 例如可 以通过 键合 的激基复合体 b o n d e d e x ci me r 中间体 的方式失 活 此外 不 同碱基之间的堆积形成的e x ci p l e x 动力 学性质更加复杂 也更具有挑战性 是我们下一步将 要开展 的研究 4 结论 采用半经典动力学方法对激光诱导下 堆积的 胸腺 嘧啶体系 的最低激发态 的失活过程进行 了模 拟 模拟脉冲频率与DF T B水平计算的 I H o 能 差相当 为4 0 e V F WH M 2 5 f s 设置激光脉冲仅作 用于一个胸腺嘧啶分子 另一胸腺嘧啶分子 T 保 持基态 激光作用 以后 约 1 2 个 电子从HOMO激发 至 L UMO 形成了T带负电荷 T 带正电荷 的激基复 合物 相邻分子的空问效应使得T分子强烈扭 曲受 到阻碍 而这种强烈的扭 曲正是到达圆锥相交所必 需 的 因此激基 复合物 的寿命 比单体长 当分子间 距离缩短至0 3 n m后 T分子 的C 一C 键扭 曲程度 最大 此 时发 生电荷 重组 两个胸腺嘧啶分子恢复 中性 电荷重组诱导T 分子发生畸变 并在C 一C 扭 曲最大时发生避免相交 体系衰减至基态 两个 胸腺嘧啶分子均恢复平面构型 Re f e r e il ce s 1 Be u k e r s R E k e r A P M L o h ma n P H M D N ARe p a i r 2 5 6 4 A ct a P h y s 一 C h i m S i n 2 0 1 1 厂0 1 2 7 f 8 1 1 9 f 2 O 1 2 1 2 2 2 3 2 4 2 5 2 6 2 7 2 8 2 0 0 8 7 5 3 0 M e l n ik ov a V O An a n t ha s wa m y H N M u t at Re s 2 0 0 5 57 1 9 1 Ca d e t J S a g e E Do u k i M u t a t Re s 2 00 5 5 7 1 3 M o u r e t S Ba d o uin C Cha r ve r o n M Fa v ie r A Ca d e t J Do u k i Pr o c Na t l Ae a d S ci U s A 2 0 0 6 1 0 3 3 7 6 5 Do u ki T Re y n a u d An g e l in A Ca d e t J S a g e E Bio ch e mis t r y 2 0 0 3 4 2 9 2 2 1 D u r b e e j B E r i k s s o n L A P h o t o ch e m Ph o t o b i o l 2 0 0 3 7 8 1 5 9 S ch r e ie r J Sch r a d e r E Ko l l e r F O Gil ch P Cr e s po He m fin d e z C E S wa min a t h a n V N Ca r e l l T Zin t h W Ko h l e r B S cie n ce 2 0 0 7 31 5 6 2 5 Sch r e ie r W J Ku b o n J Re g n e r N Am Ch e m S o c 2 0 0 9 1 3 1 5 0 3 8 Z h a n g R B E r i k s s o n L A J P h y s C h e m B 2 0 0 6 1 1 0 7 5 5 6 Du r b e e j B E r i k s s o n L A Ph o t o ch e m P h o t o b i o 1 A 2 0 0 2 1 5 2 9 5 Bo g g io Pa s q u a M Gr o e n h o f G Sch a f e r L V Gmb mfil l e r H Ro b b M A Am Ch e m S o c 2 0 0 7 1 2 9 1 0 9 9 6 Bl a n ca f o r t L M ig a n i A Am Che m S o c 2 0 0 7 1 2 9 1 4 5 4 0 La w Y K Az a d i J Cr e s p o He m n d e z C E Ol mo n E Ko h l e r B Bio ph y s 2 0 0 8 9 4 3 5 9 0 J o h ns o n A W ie s t O P Ch e m B 2 0 0 7 1 1 1 1 4 3 9 8 M cCul l a g h M Ha r ih a r a n M Le wis F D M a r k o v it s i D Do u k i S ch a t z G C 尸 Ch e m B 2 0 1 0 1 1 4 5 21 5 E is in g e r J La m o l a A Bio ch e m Bio ph y s Re s Co mmu n 1 9 6 7 2 8 5 58 Eis in g e r J Sh ul ma n R G Pr o c Na t 1 Aca d S ci A 1 9 67 5 8 8 95 Eis in g e r J Gu 6 r o n M S h u l ma n R G Ya ma n e Pr o c Na t 1 Aca d S c A 1 9 6 6 5 5 1 01 5 Bir k s J B Na t u r e 1 9 6 7 21 4 1 1 8 7 T a k a y a T S u C d e La Ha r pe K Cr e s po He m find e z C E Ko h l e r B Pr o c Na t 1 Ac a d S e i A 2 0 0 8 1 0 5 1 02 8 5 Co n t i I Al t o e P S t e n t a M Ga r a v e l l i M Or l a nd i G P C h e m C h e m P h y s 2 0 1 0 1 2 5 0 1 6 Cr e s po He m n d e z C Ko h l e r B Nat u r e 2 0 0 5 4 3 6 1 1 41 S ch wa l b N K T e mp s F S cie n ce 2 0 0 8 3 2 2 2 4 3 Kwo k W M M a C Ph il l ip s D L Am C he m S o c 2 0 0 6 1 2 8 1 l 8 9 4 Cr e s p o He m md e z C E Ko h l e r B Ph y s Ch e m B 2 0 04 1 0 8 1 1 1 8 2 Kwo k M M a C P h il l ip s D L Ch e m B 2 0 0 9 1 1 3 1 1 5 2 7 Ho l m A I S Nie l s e n L M Ko h l e r B Ho ffma n n S V Nie l s e n S B Ph y s Ch e m Ch e m Ph y s 2 01 0 1 2 3 4 2 6 Co h e n B La r s o n M H Ko h l e r B C h e m P 2 0 0 8 3 5 0 1 6 5 f 2 9 3 0 3 1 3 2 3 3 1 3 4 3 5 3 6 3 7 f 3 8 f 3 9 f 4 0 1 4 1 4 2 1 f 43 1 f 4 4 f 45 1 4 6 4 7 4 8 1 f 4 9 5 O 5 1 5 2 Do u 55 T o r r a l v a B Al l e n R M o d Opt ics 2 0 0 3 5 0 2 61 5 Do u 5 5 T o r r a l v a B Al l e n R Ch e m Ph y s L e t t 1 9 9 8 3 7 8 3 2 3 Gr a f M Vo g l P Re v B 1 99 5 5 1 4 9 Ha u g k M Eis ne r J F r a u e n h e im T S e if e r t G S t e i n b e r g M P h y s S t a t u s S o l i d iB 2 0 0 0 2 1 7 4 7 3 F r a ue n h e im S e if e r t G El s t n e r M Nie h a u s A KS h l e r C Am k r e u t z M S t e r n b e r g M Ha j n a l Z Di C a r l o A S u h a i S J P h y s C o n d e n s Ma t e r 2 0 0 2 1 4 3 0 1 5 W a n k o M Ga r a v e l l i M Be r n a r d i F Nie ha u s T A F r a u e n h e im T El s m e r M Ch e m Ph y s 2 0 0 4 1 2 0 1 6 7 4 Zh e ng G Lu n d b e r g M J a k o ws k i J Vr e ve n Fr is ch M J Mo r o k u ma K I n t Q