姜黄素与邻苯二酚共晶的太赫兹光谱.pdf

姜黄素与邻苯二酚共晶的太赫兹光谱*郑转平刘榆杭赵帅宇蒋杰伟卢乐(西安邮电大学电子工程学院,西安710121)(2023年 5月 6 日收到;2023年 6月 21 日收到修改稿)姜黄素是一种具有抗炎、抗氧化及抗癌等作用的常用药物,但是其在水中的溶解度较低.近年来,药物共晶是一种增强水溶性有限药物溶解度和溶解性的有效方法.基于此,本文利用太赫兹时域光谱技术(THz-TDS)研究了姜黄素与邻苯二酚共晶的太赫兹光谱.首先测试了姜黄素、邻苯二酚、二者物理混合及其共晶在 0.53.5THz 的实验谱.实验数据显示共晶获得的位于 3.31THz 处的吸收峰明显区别于原料物质,表明THz-TDS 可以有效鉴别姜黄素、邻苯二酚及其共晶体;基于密度泛函理论,对姜黄素与邻苯二酚共晶体可能存在的 4 种理论晶型进行结构优化与频谱模拟,其中理论晶型模拟结果与实验数据比较符合.研究发现,共晶体是通过姜黄素的羰基 C10=O3 与邻苯二酚的羟基 O61H55 形成氢键而成,共晶的 THz 特征吸收峰是在氢键带动下由共晶中两分子的官能团共同作用产生的.关键词：姜黄素,共晶,太赫兹波,吸收谱PACS：32.30.Bv,31.15.E,33.20.TpDOI:10.7498/aps.72.202307391引言药物共晶是药物活性成分(API)与共晶形成物(CCF)在氢键和其他非共价键的作用下结合而成的晶体.它是一种重要的多组分复合物,在同一晶格中压缩至少两个由非共价键连接的分子.研究表明,药物共晶的结构在一定程度上会对药物分子的理化性质产生影响,可以显著改变药物活性成分的物理化学性质,如溶解度、溶解速率、稳定性、渗透性、可食性和生物利用度等,使得其在药物开发中有巨大的潜在价值13.现阶段,研究共晶的分子结构排布及其微观的作用模式对药物应用领域具有重要的现实意义46.传统研究中,多采用X衍射、红外与拉曼光谱仪来检测共晶,但需要色谱法等辅助表征晶体结构7,8.而振动光谱技术源于分子振动能级的跃迁,并伴有转动能级改变,被广泛用于物质的表征辨别中9.太赫兹波作为转动光谱,其波长介于微波与红外波之间,其波段范围位于 0.110THz,由于光子能量与大部分药物分子的转动振动能级之间跃迁的能量相当,能够敏锐地探测到分子低频振动区域内的偶极跃迁,并对分子间弱相互作用力(氢键、范德瓦耳斯力)及分子官能团振动形成“指纹光谱”,受到了越来越多的关注10.近年来,由于太赫兹波较强的穿透性和光子能量低等特点,基于太赫兹波的太赫兹时域光谱技术(THz-TDS),已经被越来越多地应用于药物共晶的检测表征与结构分析中11,12.如国内 Cai 等12利用太赫兹时域光谱技术在室温下研究了阿昔洛韦、富马酸、阿昔洛韦与富马酸共晶体的 THz 光谱,发现共晶体与原料物质的吸收峰完全不同;Wang 等13研究了药物分子异烟酰胺与硝苯地平共晶体的两种晶型,测定了共晶在 85325K 的 THz 光谱,结果显示晶型和在 THz 区域有“指纹光谱”.国外Margaret等14研究了药物分子马来酸依那普利*国家自然科学基金(批准号:62276210)和陕西省自然科学基础研究计划(批准号:2022JM-380)资助的课题.通信作者.E-mail:2023中国物理学会ChinesePhysicalSocietyhttp:/物理学报ActaPhys.Sin.Vol.72,No.17(2023)173201173201-1共晶的两种多形态的 THz 光谱,并结合理论模拟分析了晶型的稳定性;Luczynska 等15基于多种光谱技术研究了溴苯甲酸与 2,6-二甲基吡嗪共晶的结构,揭示了共晶体系中形成了中等强度的 OHN 氢键.姜黄素(curcumin,CUR)是一种天然的多酚类化合物,是姜黄的主要活性成分之一,具有抗炎、抗氧化、抗肿瘤及保护肝脏和肾脏等多种药理活性,且对人体的毒副作用很小,非常适合于药物医学方面的应用1618.但由于其性质不稳定,在碱性条件下主要以酮-稀醇的形式存在,在酸性中性的条件下主要以 b-二酮的形式存在,难溶于水,在水中的溶解度仅为 7.8mg/L19.其生物利用率低、稳定性差等理化性质极大的限制了 CUR 在药物方面的应用发展20.在不改变 CUR 药物活性成分的前提下改善 CUR 的理化性质,已经成为近年来 CUR 研究者们主要待解决问题.已有研究表明,CUR 共晶体在水中溶解度高于 CUR 本身21.但是目前有关 CUR 共晶的研究主要集中在共晶体的制备以及对共晶溶解度、溶解速率、稳定性等理化性质22.而对 CUR 共晶体的分子结构及光谱表征的研究一直是有待解决的问题.有机化合物邻苯二酚(CTL)是一种重要的化工中间体,其生产方便成本低廉.邓彩赟等23利用物理研磨法制备了 CUR 与 CTL 共晶,采用溴化钾压片法测试并表征了 4000400cm1波数的红外光谱数据,根据实验数据,推测出 CUR 与 CTL共晶中存在OHO 氢键.但是,对于氢键的具体连接位置以及形成的共晶分子中 CUR 与 CTL的分子个数比例尚不清楚.因此本工作将以 CTL作为姜黄素共晶合成物,采用 THz-TDS 对 CUR、CTL、CUR-CTL 共晶、CUR-CTL 物理混合进行THz 光谱测试;通过原料物质与共晶体的 X 衍射光谱验证,发现基于 THz 特征吸收峰,THz-TDS可以实现对共晶体的定性表征;构建共晶可能存在的 4 种晶型结构,基于 DFT 理论,对其进行结构优化和光谱计算,通过与实验光谱结果进行匹配,推测共晶体的晶型结构,并对共晶体的吸收峰来源进行振动归属.2实验与理论 2.1 样品制备实验所用的 CUR(纯度98%,CAS:458-37-7)和 CTL(纯度99%,CAS:120-80-9)均购买于上海阿拉丁生化科技股份有限公司,后续使用均未进一步纯化.在共晶的制作过程中,首先将 CUR与 CTL 分别研磨成粉末状;接着将 CUR 和 CTL按照物质的量比 11 投料,称量 CUR(0.921g,2.5mmol)和 CTL(0.275g,2.5mmol)加入研钵,室温研磨 30min 得到 CUR 和 CTL 共晶粉末,在研磨过程中固体粉末从原来 CUR 的黄色逐渐变成橘红色;最后称量共晶粉末 100mg,在 12MPa 的压力下制成测试样品片.在物理混合制备过程中,分别取 CUR 和 CTL 适量加入研钵,在室温下充分研磨,称量研磨好的 CUR(50mg)和 CTL(50mg)加入研钵搅拌均匀,混合样品呈黄色,在 12MPa的压力下制成测试样品片.2.2 实验装置X 衍射谱测试采用 SmartLab9kW(日本理学 RIGAKU).其中仪器最大输出功率为 9kW,额定电压为 60kV,工作电流为 200mA,扫描角度2q 范围为 590.51C49C52C57H50C53C56H59H58H48C61O60O55H54H14H1O7C8C28C29H9C24C25H26C27H30C31H32C33H10C3O34C36H35H11C4O37C38H39C40H41C42C6O46C44C45H12C13C47H20C21H23H22H5O19H8H16C17H2O15H(a)(b)图1CUR(a)与 CTL(b)的分子结构Fig.1.MolecularStructureofCUR(a)andCTL(b).物理学报ActaPhys.Sin.Vol.72,No.17(2023)173201173201-2THz 光谱测试采用莱仪特太赫兹(天津)科技有限公司的光谱仪(Advantest,TAS7500SU).实验使用的是透射式太赫兹时域光谱系统,其中飞秒激光器产生波长为 8001550nm 的超短激光脉冲,分别采用光电导天线和光电采样方法产生和探测 THz 脉冲,系统测试光谱范围为 0.56.0THz,光谱分辨率 7.8GHz,测试环境为室温(295K).为了减少空气中的水分对太赫兹波的吸收,常将THz-TDS 系统罩在充满干燥氮气的箱体内,以减小空气中水分对 THz 波的吸收和影响,实验中相对湿度低于 2%.2.3 理论计算理论模拟采用原子轨道线性组合方法(linearcombinationofatomicorbitals,LCAO)24.它是量子化学中用于求解分子轨域的一种方法,通过对原子轨域进行线性叠加来构造分子轨域,继而更有效地描述电子结构、结合及反应性质.具体的计算方法是 B3LYP25,用到的基组为 6-311+G(d,p),收敛标准为 tight,计算结果没有使用标度因子修正.图 1 为 CUR 与 CTL 的分子结构图,本文在构建CUR-CTL 共晶可能存在的结构体系基础上,采用量化密度泛函理论(densityfunctionaltheory,DFT)对共晶体进行理论计算26.DFT 由于包含了相关能效应,因此能够比 Hartree-Fock 方法更精确地描述研究对象,但由于相关泛函的精确形式未知,因此采用B3LYP 杂化泛函,在6-311+G(d,p)基组下对共晶体进行结构优化和频率计算.通过计算结果与实验结果进行比对,推测出 CUR-CTL共晶体的可能晶型.3实验结果与讨论 3.1 CUR,CTL,CUR-CTL 共晶及其物理混合的 THz 实验谱图 2 是 CUR,CTL,CUR-CTL 共晶与物理混合的太赫兹频域吸收光谱.由图可见,CUR 获得6 个 THz 吸收峰,分别位于 0.67,1.03,2.05,2.30,2.60,3.07THz,其中 0.67THz 和 1.03THz 属于弱吸收峰;CTL 获得了 3 个 THz 吸收峰,分别位于 1.95,2.49,3.41THz,其中 1.95THz 属于强吸收峰;CUR-CTL 物理混合物获得了 4 个吸收峰,分别位于0.67,1.97,2.49,3.06THz,其中1.97THz和2.49THz 与CTL 的THz 吸收峰匹配,0.67THz和 3.06THz 与 CUR 的吸收峰匹配;CUR-CTL 共晶获得了 3 个吸收峰,分别位于 1.03,1.95THz和 3.31THz,其中 1.03THz 与姜黄素吸收峰匹配,1.95THz 与邻苯二酚的THz 吸收峰匹配,3.31THz是产生的新THz 吸收峰.由实验结果可见,CUR-CTL 物理混合在0.67THz 的吸收峰主要来自于 CUR 的 0.67THz,1.97THz 吸收峰主要来自于 CTL 在 1.95THz 和CUR 在 2.05THz 的叠加,2.49THz 吸收峰主要来自于CTL 在2.49THz 和CUR 在2.30THz 的叠加,3.06THz 的吸收峰主要来自于CUR 的3.07THz,物理混合的 THz 吸收峰均来自于 CUR 与 CTL,并未产生新的吸收峰,说明二者的物理混合中 CUR与 CTL 各自的分子结构不变,分子间不存在相互作用力.在 CUR-CTL 的共晶体中,出现了姜黄素的 弱 吸 收 峰 1.03THz,邻 苯 二 酚 的 强 吸 收 峰1.95THz,以及新的吸收峰 3.31THz,说明共晶中姜黄素和邻苯二酚之间产生了新的分子三维结构,分子间产生新的弱相互作用力.这种作用力的存在使得姜黄素的弱吸收峰(1.03THz)出现在共晶中,邻苯二酚 1.95THz 的强吸收峰在共晶中的吸收幅度减弱.此结果表明 CUR-CTL 共晶区别于原料物质,晶体的物象结构相比于原料物质已经发生了变化.说明基于THz-TDS 技术可以对 CUR-CTL共晶进行表征.3.2 CUR,CTL,CUR-CTL 共晶的 XDR图谱图 3 是 CUR,CTL 及 CUR-CTL 共晶的 X 衍0.51.01.52.02.53.03.5Absorption/arb.unitsFrequency/THz3.313.061.972.491.952.493.412.052.302.603.071.031.031.950.670.67CUR-CTLPhysical mixtureCTLCUR图2CUR,CTL,CUR-CTL 物理混合、CUR-CTL 共晶的 THz 谱Fig.2.THzspectraofCUR,CTL,CUR-CTLphysicalmix-ing,andCUR-CTLco-crystal.物理学报ActaPhys.Sin.Vol.72,No.17(2023)173201173201-3射实验谱.由图3 可见,CUR 在7.88,13.94,15.56,16.34,17.47,26.16,27.02,28.14和 28.96的特征峰在共晶中消失,CTL在 16.22,20.02,23.6,24.52,24.86,32.38,32.78,33.58,33.96,34.42的特征峰在共晶中消失,而 CUR-CTL 共晶在11.22,13.6,20.86和 21.56出现了新的特征峰20.此结果说明 CUR-CTL 共晶体不是两种原料物CUR 和 CTL 的简单混合,而是形成了自己独有的衍射峰,表明 CUR-CTL 共晶相比于原料物质其物相结构发生了变化,确认了共晶样品的生成.此实验结果验证了基于 THz-TDS 可以辨别检测CUR-CTL 共晶体的准确性.4理论模拟及结果分析 4.1 共晶的结构模型根据 CUR 与 CTL 的氢键供受体位置以及氢键合成规律,本文模拟了 CUR-CTL 共晶体的 4 种可能晶型,如图 4 所示.红色阴影表示 CUR 与CTL 的氢键连接结构.其中,理论晶型(图 4(a)是 CUR 的羰基 C11=O4,C10=O3 分别与 CTL的羟基 O61H55,O60H54 通过双氢键连接而成,理论晶型(图 4(b)是 CUR 的羰基 C11=O4与 CTL 的羟基 O61H55 通过单氢键作用而成,理论晶型(图 4(c)是 CUR 的羰基 C10=O3 与5101520253035CTLCURCUR-CTLIntensity/counts2/(O)图3CUR,CTL,CUR-CTL 共晶的 X 衍射实验谱Fig.3.X-ray diffraction experimental spectra of CUR,CTL,andCUR-CTLco-crystal.19H17H16C18H5O14H1O8C7C28C25H29H9C26C27H30C31H32C33H10C34C35H11C36H37C38H39C40H43H45H41C42C44C12C13C46C47H2O15H21H20C23H22H3O4O55H54H61O60O49C40C53C59H50C51C52C58H57H56H19H43H42C45H44C12C2O15H13C46C41C39C37C47H6O20C22H23H38H11C34C35H10C32C30C31H3O4O55H61O60O54H48C49C50C56H53C59H52C51C58H57H9C26C27H24C25H7C1O14H8C28C29H5O16C18H17H(a)(b)15H41C13C46C47H5O20C12C2O44C45H42C22H21H39C43H23H4O11C38H37C34C35H10C32C30C31H60O3O61O54H51C61O48C49C50C53C56H57H52C58H59H27H1O9C28C26C25H29H8C7C24C14H5O16C19H18H17H15H2O46C6O12C44C43H41C13C20C11C45H42C23H22H39C37C38H34C35H47H10H3O4O32H30H31H27H25H9C26C24C7C16H17H18H52C1O14H8C28C29H5O19H69H55H61O49C48C53C59H10O54H58H57H51C75O50C56H65C72H73H71H64C70H63C74O68H62C67C66C(c)(d)图4CUR-CTL 共晶分子结构(a)晶型;(b)晶型;(c)晶型;(d)晶型Fig.4.MolecularstructuresofCUR-CTLco-crystal:(a)FormI;(b)formII;(c)formIII;(d)form.物理学报ActaPhys.Sin.Vol.72,No.17(2023)173201173201-4CTL 的羟基 O61H55 通过单氢键作用而成,理论晶型(图4(d)是CUR 的羰基C11=O4,C10=O3 分别与两个 CTL 分子的羟基 O75H69,O61H55 通过双氢键连接而成.4.2 太赫兹光谱理论分析图 5 是 CUR-CTL 共晶体的理论模拟结果与实验谱.其中,理论数据由洛伦兹函数展宽,半高宽为 5cm1.由图 5 可见,晶型,的模拟数据与实验结果相差较大,晶型除了 2.06THz 与实验的 1.95THz 匹配较好以外,其余的模拟结果与实验谱相差较大.晶型在研究范围内获得了 5 个光学模式,分别位于 1.03,1.95,2.47,2.71THz 和3.34THz,其中 1.03,1.95THz 和 3.34THz 分别与实验吸收峰 1.03,1.95THz 和 3.31THz 相匹配,理论数据与实验峰匹配较好.但是,晶型中获得的 2.47THz 和 2.71THz 的光学模式并未对应实验峰.可能的原因一是理论计算中采用密度泛函理论 B3LYP 方法计算的频率是谐振频率,并没有考虑非谐振效应的影响,且选择的基组也会在一定程度上影响分子体系能量计算的准确性,产生高估或低估的现象,从而出现实验与理论频谱数据的偏差;二是理论结构优化是在绝对 0K 下进行的,而实验测试在常温下(295K).Otsuka等27研究得出,随着温度的升高,样品 THz 谱的分辨率随之下降,容易导致较小的吸收峰被淹没无法分辨.由理论数据发现,晶型的分子结构能量在 4 种晶型中最低,分子结构最稳定,由此推断 CUR-CTL 共晶体系为理论晶型,CUR 的羰基(C10=O3)与CTL 的羟基(O61H55)通过氢键(O61H55O3C10)连接形成新的晶体.表 1 为共晶形成前后分子的键长键角.如表 1所示,氢键C10=O3H55O61 的形成使得C10=O3(CUR)的键长由0.12984nm 缩短到0.12131nm,O61H55(CTL)的键长由 0.08141nm 伸长到了0.09698nm,使得键角C34C10O3 由 117.96300增至 124.79870,键角C49O61H55 由 111.14600增至 116.49998,其中 C61H55 的键长变化与邓彩赟等23的研究成果一致.由此可见,氢键的生成改变了 CUR 与 CTL 分子内部的结构参数.表1姜黄素与邻苯二酚键长和键角Table1.Bondlengthandbondangleofcurcuminandcatechol.Bondlength/nmCUR28CTL29Co-crystalC10=O30.129840.12131O61H550.081410.09698Bondangle/()CURCTLCo-crystalC34C10O3117.93600124.79870C49O61H55111.14600116.49998表 2 为 CUR-CTL 共晶体实验数据与理论晶型的计算结果.基于理论与实验数据匹配,CUR-CTL 共晶体在 1.03THz(图 6(a)处的振动模式来自于 CUR 分子中碳链 C34C10C32 的摆动,该摆动带动苯环R2 面外摆动,并通过氢键C10=O30H55O61 带动 CTL 分子的面内摆动;1.95THz(图 6(b)处的振动模式主要来源于 CUR 分子中-C16H3与苯环 R2 的面外扭动,同时通过氢键C10=O30H55O61 带动 CTL 分子轻微的面外摆动;3.34THz(图 6(c)处的模式主要是来源于氢键 C10=O30H55O61 的伸缩和 CUR 碳链C30C32C10 的摆动,同时带动苯环 R2 与 R1面外摆动.通过对共晶吸收峰来源的归纳,发现 1.03THz吸收峰是在氢键的带动下由邻苯二酚面内摆动与姜黄素的碳链、苯环 R2 共同作用产生,其中邻苯二酚作用模式比较明显,所以使得姜黄素的弱吸收峰在共晶中重现;1.95THz 的吸收峰中,相较于姜黄素苯环 R2 和甲基 CH3的较明显扭动,邻苯二酚的作用模式明显较弱,所以使得邻苯二酚的强吸收峰 1.95THz 在共晶中吸收幅度减弱;3.34THz 属于共晶体分子间作用产生的新吸收峰,主要是氢键和碳链作用带动苯环引起的.总体来说,共晶中的3 个吸收峰的来源都不是单个分子的作用,而是在Frequency/THz0.51.01.52.02.53.03.53.341.032.472.711.782.543.063.392.562.872.533.041.952.06CUR-CTL experimentalForm IVForm IIIForm IIForm IAbsorption/arb.units3.311.951.03图5CUR-CTL 共晶的实验与模拟 THz 图谱Fig.5.Experimental and simulated THz spectra of CUR-CTLco-crystal.物理学报ActaPhys.Sin.Vol.72,No.17(2023)173201173201-5氢键带动下由共晶中两分子的官能团共同作用产生的,氢键等弱相互作用力的存在不仅改变了两个分子的结构参数,而且重新建立了新的分子间作用力,继而影响了共晶的作用模式,使得 CUR-CTL共晶在太赫兹波段呈现出与原料试剂不同的THz 吸收峰.5结论本文采用 THz-TDS 测量了 CUR,CTL 及其物理混合和共晶的太赫兹光谱.实验数据表明,CUR,CTL 及其物理混合在测试范围内分别获得6 个、3 个及 4 个 THz 吸收峰;共晶获得 3 个吸收峰,分别位于 1.03,1.95THz 和 3.31THz,此吸收峰有别于其原料物质及其物理混合样品,说明共晶相比于原料物质形成了新的物象结构.此外,借助X 衍射实验谱验证了 THz-TDS 可以很好完成共晶的表征.基于密度泛函理论,模拟了 CUR-CTL共晶可能存在的 4 种晶型结构,研究结果表明共晶体系是由单个姜黄素和邻苯二酚分子通过单氢键连接而成.最后,基于 THz 实验吸收峰与理论数据的良好匹配效果,分析归纳了共晶 THz 吸收峰的来源,发现其 THz 吸收峰不是来源于单个分子的振动,而是通过氢键 C10=O30H55O61 由两个分子的官能团共同振动引起的.本研究将为CUR-CTL 共晶的表征、共晶分子的结构体系及低频作用模式的确认提供有效的参考数据.参考文献 KulkarniA,SheteS,HolV,BachhavR2020AJPS131041PrabhakarP,GiridharS2019Indian J.Pharm.Educ.535632NovenaL.M,AthimoolamS,AnithaR,BahadurSA2021J.Mol.Struct.124920223AkshitaJ,MansiP,NeelimaD,KunalC,Maninder2022J.Pharm.11127884George C P,Thorat S H,Shaligram P S,Suresha P R,5表2CUR-CTL 共晶体实验和理论计算数据Table2.TheresultsofexperimentalandcalculateddataofCUR-CTLco-crystal.Frequency/THzVibrationalmodeExpt.Theo.1.031.03BendingofC34C10C32,Out-of-planeoscillationofR2andR2,InplaneoscillationofCTL1.951.95Outofplanetorsionof-C16H3andR2,OutofplanerockingoscillationofCTL3.313.34ThestretchingvibrationofhydrogenbondC10=O30H55O61,SwingofC30C32C10,Out-of-planeoscillationofR2andR123H13C21H47H6O20C46C12C2O15H41C37C38H45H42C44C4O39C43H11C40H34C10C3O35H36H32C(a)61O55H54H50O49C48C50C51C52C53C59H1O58H57H33H56H30C9C26C24C7C8C28C27H25H14H5O16C19H18H17H29HR3R2R1(b)23H13C21H6O20C46C41C42C44C12C39C37C11C34C10C32C30C9C26C24C7C8C28C16C15H45H47H43H38H40H35H36H55H54H59H58H57H56H33H50C51C52C53C48C49C60O51O4O3O1O5O2O29H17H18H19H25H14HR3R1R227H(c)19H17H16C5O14H1O8C7C28C29H24C25H26C9C27H30C31H32C33H35H10C3O4O54H60O61O55H48C49C50C56H53C52C51C57H34C59H58H11C37C38H36H40H39C41C45H42C43H44C12C13C46C47H2O6O15H20C21H22H23HR1R2R3R1图6CUR-CTL 共晶在不同峰位的振动模式(a)1.03THz;(b)1.95THz;(c)3.34THzFig.6.VibrationalmodesofCUR-CTLco-crystalatdiffer-entpeaks:(a)1.03THz;(b)1.95THz;(c)3.34THz.物理学报ActaPhys.Sin.Vol.72,No.17(2023)173201173201-6GonnadeRG2020Crystengcomm226137SakamotoN,TsunoN,KoyamaR,GatoK,TitapiwatanakunV,TakatoriK,FukamiT2021Chem.Pharm.Burl.6919956LiuLX,LiuM,ZhangY,FengY,WuL,ZhangL,ZhangYJ,LiuYL,ZouDY2022J.Mol.Struct.125013187YongD,HongXF,QiZ,HuiLZ,ZhiH2016Spectrochim.Acta.A1535808WanM,FangJY,XueJD,LiuJJ,QinJY,HongZ,LiJS,DY2022Int.J.Mol.Sci.2385509DuanTC,YanSj,ZhaoY,SunTY,LiYG,ZhuZ2021Acta Phys.Sin.70248702(inChinese)段铜川,闫韶健,赵妍,孙庭钰,李阳梅,朱智2021物理学报7024870210ZhengZP,LiuYH,ZengF,ZhaoSY,ZhuLP2023ActaPhys.Sin.72083201(inChinese)郑转平,刘榆杭,曾方,赵帅宇,朱礼鹏2023物理学报7208320111CaiQ,XueJD,WangQ,DuY2017Spectrochim.Acta.A1862912WangPF,ZhaoJT,ZhangYM,ZhuZJ,LiuLY,ZhaoHG,YangXC,YangXN,SunXH,HeMG2022Int.J.Pharmaceut.62012175913Margaret D,Mizuki M,Kei H,Timothy K 2020 J.Phys.Chem.A124979314Luczynska,Katarzyna,Druzbicki,Kacper,Lyczko,Krzysztof,Dobrowolski2015J.Phys.Chem.B119685215Adrieli S H,Matheus G L,Radharani B 2022 Ind.Crop.16Prod.177114501YaoTM,SrinivasJW2022Food Hydrocolloids12610746617SaffarionpourS,DiosadyLL2022Curr.Opin.Food.Sci.4315518SakinehM,AliHA,PouyaM2020Chromatographia83129319ChengGL,DengGY,JiangCJ2018Chin.J.Mod.DrugApp.355(inChinese)程桂林,邓彩赟,蒋成君2018中国现代应用药学35520RibasMM,SakataGS,SantosAE,MagroCD,AguiarGP,VladimirML2019J.Supercrit.Fluid.15210456421Indumathi S,Jenna S M,Sohrab R,Sameer D V 2018 J.Chem.Data.63365222DengCY,JiangCJ2018J.Zhejiang Univ.Sci.B3016(inChinese)邓彩赟,蒋成君2018浙江科技学院学报301623RoothanCC1951Rev.Mod.Phys.236924StephensPJ,DevlinFJ,ChabalowskiCF1994J.Phys.Chem.981162325ChenT,YuLX,TangZQ,LiZ,HuFG2022Chem.Phys.56211167626OtsukaYT,ItoA,TakeuchiM,SasakiT,TanakaHJ2020J.Brug Deliv.Sci.Tec.5610121527HjorthHT,JanK,ArvidM,BertilS,ZnzellCR,EricJB1982Acta Chem.Sca.3647528WunderlichH,MootzD1971Acta Cry.Sect.B27168429物理学报ActaPhys.Sin.Vol.72,No.17(2023)173201173201-7Terahertz spectra of curcumin and catechol co-crystals*ZhengZhuan-PingLiuYu-HangZhaoShuai-YuJiangJie-WeiLuLe(School of Electronic Engineering,Xian University of Posts&Telecommunications,Xian 710121,China)(Received6May2023;revisedmanuscriptreceived21June2023)AbstractCurcumin(CUR)isacommonlyusedpharmaceuticalwithanti-inflammatory,antioxidantandanti-cancereffects,but its solubility in water is relatively low.In recent years,pharmaceutical co-crystal has been aneffectivemethodofenhancingthesolubilityoflimitedwater-solublepharmaceuticals.Basedonthis,terahertztime-domainspectroscopy(THz-TDS)isusedtostudytheTHzspectraofcurcumin-catecholco-crystal.Firstly,theexperimentalspectraofcurcumin,catechol(CTL),theirphysicalmixtureandtheirco-crystalaremeasuredinarangeof0.53.5THz,respectively.TheexperimentaldatashowthatCURobtainssixTHzabsorptionpeaks,whileCTLpossessthreeTHzabsorptionpeaks,thephysicalmixtureobtainsfourabsorptionpeaks,andtheirCUR-CTLco-crystalobtainsthreeabsorptionpeaks.TheseresultsindicatethatTHz-TDScaneffectivelyidentifycurcumin,catecholandtheirco-crystals.Thefactthattheabsorptionpeakat3.31THzobtainedinco-crystal is entirely different from those of raw materials,implying that new weak interactional forces aregeneratedbetweenCURmoleculeandCTLmolecule,theco-crystalformsanewthree-dimensionalstructurecomparedwiththeirrawmaterials.TheseresultsarealsoverifiedbyX-raydiffractionspectraofrawmaterialandtheirCo-crystal.Moreover,fourpossibletheoreticalformsofcurcumin-catecholco-crystalareoptimizedandsimulatedbyusingdensityfunctionaltheory(DFT).Thecalculatedresultsindicatethatthedataofco-crystal form III are in good agreement with the experimental spectrum,and the simulation effectivelyreconstructs the experimental spectrum.So it can be inferred that the co-crystal is formed through thehydrogen bond between the carbonyl C10=O3 of CUR and the hydroxy O61-H55 of CTL.In addition,dependingonthegoodmatchbetweenexperimentaldataandtheoreticalresults,itisfoundthatthethreeabsorptionpeaksintheco-crystaldonotoriginfromtheactionofasinglemolecule,butthejointactionofthefunctional groups of the two molecules under the driving by the hydrogen bond.The existence of weakinteractionforces,suchasthehydrogenbond,notonlychangesthestructuralparametersofthetwomolecules,butalsoreestablishesanewintermolecularforce,whichthenaffectstheinteractionalmotionso