运用密度泛函活性理论的信息论方法研究苯并富烯衍生物的芳香性

物州　化学学报　Ａｃｔａ　Ｐｈ　’ｓ．一Ｃｈｉｍ．Ｓｉｎ．２０１８，３４（６），６３９—６４９　６３９　［Ａｒｔｉｃｌｅ］　ｄｏｉ：１０．３８６６／ＰＫＵ．ＷＨＸＢ２０ｌ　７ｌ０２３１　Ｉｌ，Ｈ　饵　ｗｈｘｂ．ｐｋｕ．ｅｄｕ．ＣＩＩ　Ａｒｏｍａｔｉｃｉｔｙ　Ｓｔｕｄｙ　ｏｆ　Ｂｅｎｚｅｎｅ－Ｆｕｓｅｄ　Ｆｕｌｖｅｎｅ　Ｄｅｒｉｖａｔｉｖｅｓ　Ｕｓｉｎｇ　ｔｈｅ　Ｉｎｆｏｒｍａｔｉｏｎ－Ｔｈｅｏｒｅｔｉｃ　Ａｐｐｒｏａｃｈ　ｉｎ　Ｄｅｎｓｉｔｙ　Ｆｕｎｃｔｉｏｎａｌ　Ｒｅａｃｔｉｖｉｔｙ　Ｔｈｅｏｒｙ　Ｙｕ　Ｄｏｎｇｈａｉ１，ＲＯＮＧ　Ｃｈｕｎｙｉｎｇ１·，ＬＵ　Ｔｉａｎ２，ＤＥ　ＰＲＯＦＴ　Ｆｒａｎｋ３，ＬＩＵ　Ｓｈｕｂｉｎ４，　，Ｋｅｙ　Ｌａｂｏｒａｔｏｒｙ　ｏｔ＂Ｃｈｅｍｉｃａｌ　Ｂｉｏｌｏｇｙ　ａｎｄ　Ｔｒａｄｉｔｉｏｎａｌ　Ｃｈｉｎｅｓｅ　Ｍｅｄｉｃｉｎｅ　Ｒｅｓｅａｒｃｈ　ｆＭｉｎｉｓｔｒｙ　ｏｌ、Ｅｄｕｃａｔｉｏｎ　ｏｆＣｈｉｎａ）．Ｃｏｌｌｅｇｅ　０ｒ　Ｃｈｅｍｉｓｔｒｙ　ａｎｄ　Ｃｈｅｒｏｉｃａｌ　Ｅｎｇｉｎｅｅｒｉｎｇ．Ｈｕｎａｒｌ　ＮＯｒｔＩｌａｌ　Ｕｎｉｖｅｒｓｉｔｙ，Ｃｈａｎｇｓｈａ　４　１　００８　１　Ｐ＿Ｒ．（、ｈｉｎａ．　Ｂｅｉｉｌｉｎｇ　Ｋｅｉｎ　Ｒｅｓｅａｒｃｈ　Ｃｅｎｔｅｒ　ｆｏｒ　Ｎａｔｕｒａｌ　Ｓｃｉｅｎｃｅｓ，Ｂｅｉｉｉｎ　ｌ　０００２２　Ｐ　Ｒ．Ｃｈｉｎａ．　、Ｒｅｓｅａｒｃｈ　Ｇｒｏｕｐ　ｏｒＧｅｎｅｒａ】Ｃｈｅｍｉｓｔ￣。ｖ（ＡＬＧＣ）．Ｖｒｉｊｅ　Ｕｎｉｖｅｒｓｉｔｅｉｔ　Ｂｒｕｓｓｅｌ（ＶＵＢ），Ｐｌｅｉｎｌａａｎ　２　ｌ（１５０　Ｂｒｕｓｓｅｌｓ，Ｂｅｌｇｉｕｎ１．　０　Ｒｅｓｅａｒｃｈ　Ｃｏｍｐｕｔｉｎｇ　Ｃｅｎｔｅｒ，Ｕｎｉｖｅｒｓｉｔｙ　ｏｆＮｏｒｔｈ　Ｃａｒｏｌｉｎａ．Ｃｈａｐｅｌ　ＨⅢ．ＮＣ　２７５９９—３４２０，ＵＳＡ．　Ａｂｓｔｒａｃｔ：　Ａｌｔｈｏｕｇｈ　ａ　Ｉａｒｇｅ　ｖａｒｉｅｔｙ　ｏｆ　ａｒｏｍａｔｉｃ　ｓｙｓｔｅｍｓ　ｈａｖｅ　ｂｅｅｎ　ｕｎｖｅｉｌｅｄ　ｉｎ　ｔｈｅ　ｌｉｔｅｒａｔｕｒｅ，ｊｕｓｔｉｆｙｉｎｇ　ｔｈｅｉｒ　ｏｒｉｇｉｎ　ｏｆ　ｓｔａｂｉｌｉｔｙ　ａｎｄ　ｕｎｄｅｒｓｔａｎｄｉｎｇ　ｔｈｅｉｒ　ｎａｔｕｒｅ　ｏｆ　ａｒｏｍａｔｉｃｉｔｙ　ｉｓ　ｓｔｉｌｌ　ａｎ　ｕｎａｃｃｏｍｐｌｉｓｈｅｄ　ｔａｓｋ．Ｉｎ　ｔｈｉｓ　ｗｏｒｋ．ｕｓｉｎｇ　ｔｏｏｌｓ　ｒｅｃｅｎｔｌｙ　ｄｅｖｅｌｏｐｅｄ　ｂｙ　ｕｓ　ｗｉｔｈｉｎ　ｔｈｅ　ｄｅｎｓｉｔｙ　ｆｕｎｃｔｉｏｎａＩ　ｒｅａｃｔｉｖｉｔｙ　ｔｈｅｏｗ　ｆｒａｍｅｗｏｒｋ．　ｗｈｅｒｅ　ｗｅ　ｅｍｐｌｏｙ　ｓｉｍｐｌｅ　ｄｅｎｓｉｔｙ　ｆｕｎｃｔｉｏｎａｌｓ　ｔｏ　ｑｕａｎｔｉｆｙ　ｍｏｌｅｃｕｌａｒ　ｓｔｒｕｃｔｕｒａｌ　—　ａｎｄ　ｒｅａｃｔｉｖｉｔｙ　ｐｒｏｐｅｒｔｉｅｓ．ｗｅ　ｅｘａｍｉｎｅ　ｔｈｅ　ａｒｏｍａｔｉｃｉｔｙ　ｃｏｎｃｅｐｔ　ｆｒ０ｍ　ａ　ｄｉｆｅｒｅｎｔ　ｐｅｒｓｐｅｃｔｉｖｅ．Ｕｓｉｎｇ　ｓｉｘ　ｑｕａｎｔｉｔｉｅｓ　ｆｒＯｍ　ｔｈｅ　ｉｎｆｏｒｍａｔｉｏｎ—ｔｈｅｏｒｅｔｉｃ　ａｐｐｒｏａｃｈ．　ｎａｍｅｌｙ．１ｈｅ　Ｓｈａｎｎｏｎ　ｅｎｔｒｏｐｙ，Ｆｉｓｈｅｒ　ｊｎｆｏｒｍａｔｉｏｎ．Ｇｈｏｓｈ—Ｂｅｒｋｏｗｉｔｚ—Ｐａｒｒ　ｅｎｔｒｏｐｙ，Ｏｎｉｃｅｓｃｕ　ｉｎｆｏｒｍａｔｉｏｎ　ｅｎｅｒｇｙ，ｉｎｆｏｒｍａｔｉｏｎ　ｇａｉｎ，ａｎｄ　ｒｅｌａｔｉｖｅ　Ｒｅｎｙｉ　ｅｎｔｒｏｐｙ，ａｎｄ　ｆｏｕｒ　ａｒｏｍａｔｉｃｉｔｙ　ｄｅｓｃｒｉｐｔｏｒｓ，ｎａｍｅｌｙ，ｔｈｅ　ａｒｏｍａｔｉｃ　ｓｔａｂｉｌｉｚａｔｉｏｎ　ｅｎｅｒｇｙ（ＡＳＥ）ｉｎｄｅｘ，ｔｈｅ　ｈａｒｍｏｎｉｃ　ｏｓｃｉｌｌａｔｏｒ　ｍｏｄｅＩ　ｏｆ　ａｒｏｍａｔｉｃｉｔｙ（ＨＯＭＡ）ｉｎｄｅｘ，ｌｈｅ　ａｒｏｍａｔｉｃ　ｆｌｕｃｔｕａｔｉｏｎ（ＦＬＵ）ｉｎｄｅｘ，ａｎｄ　ｔｈｅ　ｎｕｃｌｅｕｓ—ｉｎｄｅｐｅｎｄｅｎｔ　ｃｈｅｍｉｃａｌ　ｓｈｉｆｔ　（ＮｌＣＳ）ｉｎｄｅｘ．ｗｅ　ｓｙｓｔｅｍａｔｉｃａｌｌｙ　ｅｘａｍｉｎｅｄ　ｔｈｅ　ｃｏｒｒｅｌａｔｉｏｎｓ　ｂｅｔｗｅｅｎ　ｓｕｂｓｔｉｔｕｔｅｄ　ｆｕｌｖｅｎｅ　ｄｅｒｉｖａｔｉｖｅｓ　ｆｕｓｅｄ　ｗｉｔｈ　ｏｎｅ，ｔｗｏ，　ａｎｄ　ｔｈ　ｒｅｅ　ｂｅｎｚｅｎｅ　ｒｉｎｇｓ．Ａｍｏｎｇ　ｔｈｅ　１　４　ｂｅｎｚｏｆｕｌｖｅｎｅ　ｄｅｒｉｖａｔｉｖｅｓ　ｓｔｕｄｉｅｄ　ｉｎ　ｔｈｉｓ　ｗｏｒｋ．ｔｈｅｒｅ　ｗｅｒｅ　ｓｅｖｅｎ　ｓｉｎｇＩｅ—ｆｕｓｅｄ．ｆｏｕｒ　ｄｏｕｂｌｅ．ｆｕｓｅｄ．ａｎｄ　ｔｈ　ｒｅｅ　ｔｒｉｐｌｅ—ｆｕｓｅｄ　ｂｅｎｚｏｆｕｌｖｅｎｅ　ｄｅｒｉｖａｔｉｖｅｓ．Ｏｕｒ　ｒｅｓｕｌｔｓ　ｓｈｏｗ　ｔｈａｌ　ｔｈｅ　ａｒｏｍａｔｉｃｉｔｙ　ｉｎｄｅｘｅｓ　ａｒｅ　ｏｆｔｅｎ　ｗｅｌＩ　ｃｏｒｒｅｌａｔｅｄ　ｗｉｔｈ　ｏｎｅ　ａｎｏｔｈｅｒ　Ｔｈｅ　ｓａｍｅ　ｉｓ　ｔｒｕｅ　ｆｏｒ　ｉｎｆｏｒｍａｔｉｏｎ—ｔｈｅｏｒｅｔｉｃ　ｑｕａｎｔｉｔｉｅｓ．Ｍｏｒｅｏｖｅｒ．ｔｈｅｓｅ　ｃｏｒｒｅｌａｔｉｏｎｓ　ａｒｅ　ｖａｌｉｄ　ａｃｒｏｓｓ　ａｌＩ　ｓｅｒｉｅｓ　ｏｆ　ｂｅｎｚｏｆｕｌｖｅｎｅ　ｄｅｒｉｖａｔｉｖｅｓ　ｗｉｔｈ　ｄｉｆｅｒｅｎｔ　ｒｉｎｇ　ｓｔｒｕｃｔｕｒｅｓ．Ｔｈｅ　ｃｒｏｓｓ　ｃｏｒｒｅｌａｔｉｏｎｓ　ｂｅｔｗｅｅｎ　ｉｎｆｏｒｍａｔｉｏｎ—ｔｈｅｏｒｅｔｉｃ　ｑｕａｎｔｉｔｉｅｓ　ａｎｄ　ａｒｏｍａｔｉｃｉｔｙ　ｉｎｄｅｘｅｓ　ｗｅｒｅ　ｕｓｕａｌｌｙ　ｓｔｒｏｎｇ．Ｈｏｗｅｖｅｒ，ｔｗｏ　ｃｏｍｐｌｅｔｅｌｙ　ｏｐｐｏｓｉｔｅ　ｐａｔｔｅｒｎｓ　ｗｅｒｅ　ｏｂｓｅｍｅｄ：ａｓ　ａ　ｃｏｎｓｅｑｕｅｎｃｅ．ｔｈｅｓｅ　ｃｏｒｒｅｌａｔｉｏｎｓ　ａｒｅ　ｎｏｔ　ｖａｌｉｄ　ａｃｒｏｓｓ　ａｌｌ　ｓｅｒｉｅｓ　ｏｆ　ｂｅｎｚＯｆｕｌｖｅｎｅ　ｄｅｒｉｖａｔｉｖｅｓ．Ｔｈｅ　ｎａｔｕｒｅ　ｏｆ　ｔｈｅｓｅ　ｃｏｒｒｅｌａｔｉｏｎｓ　ｄｅｐｅｎｄｓ　ｏｎ　ｔｈｅ　ｎａｔｕｒｅ　ｏｆ　ｔｈｅ　ｒｉｎｇ　ｓｔｒｕｃｔｕｒｅ．Ｔｈｅ　ｔｗｏ　ｇｒｏｕｐｓ　ｏｆ　ｓｙｓｔｅｍｓ，ｅａｃｈ　ｏｂｅｙｉｎｇ　ｔｈｅ　ｓａｍｅ　ｃｒｏｓｓ—ｃｏｒｒｅｌａｔｉｏｎ　ｐａＲｅｒｎｓ．ｈａｖｅ　ａ　ｔｏｔａｌ　ｎｕｍｂｅｒ　ｏｆ　４ｎ　４－２　ａｎｄ　４ｎ丌ｅｌｅｃｔｒｏｎｓ．ｒｅｓｐｅｃｔｉｖｅｌｙ．ｗｈｉｃｈ　ａｒｅ　ｉｎ　ａｇｒｅｅｍｅｎｔ　ｗｉｔｈ　Ｈ０ｃｋｅｌ　ｓ　ｒｕｌｅ　ｏｆ　ａｒｏｍａｔｉｃｉｔｙ　ａｎｄ　ａｎｔｉａｒｏｍａｔｉｃｉｔｙ．Ｃｏｍｐａｒｅｄ　ｗｉｔｈ　ｔｈｅ　ｒｅｓｕｌｔｓ　ｏｂｔａｉｎｅｄ　ｆｏｒ　ｓｙｓｔｅｍｓ　ｗｉｔｈｏｕｔ　ａ　ｂｅｎｚｅｎｅ　ｆｕｓｅｄ　ｒｉｎｇ，ｔｈｅ　ｃｏｒｒｅｌａｔｉｏｎ　ｐａｔｔｅｒｎｓ　ｏｆ　ｔｈｅｓｅ　ｑｕａｎｔｉｔｉｅｓ　ｗｅｒｅ　ａｌｗａｙｓ　ｆｏｕｎｄ　ｔｏ　ｂｅ　ｔｈｅ　ｓａｍｅ，ｂｏｔｈ　ｗｉｔｈ　ａｎｄ　ｗｉｔｈｏｕｔ　ｆｕｓｅｄ　ｂｅｎｚｅｎｅ　ｒｉｎｇｓ　Ｔｈｉｓ　ｓｕｇｇｅｓｔｓ　ｔｈａｔ．ｄｅｓｐｉｔｅ　ｂｅｎｚｅｎｅ’ｓ　ａｒｏｍａｔｉｃｉｔｙ．ｉｔｓ　ｆｕｓｉｏｎ　ｗｉｔｈ　ａ　ｆｕｌｖｅｎｅ　ｍｏｉｅｔｙ　ｄｏｅｓ　ｎｏｌ　ｍｏｄｉｆｙ　ｔｈｅ　ａｒｏｍａｔｉｃｉｔｙ　ａｎｄ　ａｎｔｉａｒｏｍａｔｉｃｉｔｙ　ｏｆ　ｔｈｅ　ｆｕｌｖｅｎｅ　ｒｉｎｇ　Ｔｈｅｓｅ　ｒｅｓｕｌｔｓ　ｃｏｎｆｉｒｍ　ｔｈａｔ　ｔｈｅ　ｆｕｓｉｏｎ　ｏｆ　ｂｅｎｚｅｎｅ　ｒｉｎｇｓ　ｗｉｔｈ　ａ　ｆｕｌｖｅｎｅ　ｍｏｉｅｔｙ　ｈａｓ　ｎｏ　ｉｎｆｌｕｅｎｃｅ　ｏｎ　ｔｈｅ　ａｒｏｍａｔｉｃ　ｎａｔｕｒｅ　ｏｆ　ｔｈｅ　ｆｕｌｖｅｎｅ　ｍｏｉｅｔｙ．Ｔｈｕｓ．ｔｈｅ　ａｒｏｍａｔｉｃｉｔｙ　ａｎｄ　ａｎｔｉａｒｏｍａｔｉｃｉｔｙ　ｏｆ　ｂｅｎｚｅｎｅ—ｆｕｓｅｄ　ｆｕｌｖｅｎｅ　ｄｅｒｉｖａｔｉｖｅｓ　ａｒｅ　ｓｏｌｅｌｙ　ｄｅｔｅｒｍｉｎｅｄ　ｂｙ　ｔｈｅ　ｆｕｌｖｅｎｅ　ｍｏｉｅｔｙ．Ｔｈｅｓｅ　ｒｅｓｕｌｔｓ　ｓｈｏｕｌｄ　ｐｒｏｖｉｄｅ　ａ　ｎｅｗ　ｕｎｄｅｒｓｔａｎｄｉｎｇ　ｏｆ　ｔｈｅ　ｏｒｉｑｉｎ　ａｎｄ　ｎａｔｕｒｅ　ｏｆ　ａｒｏｍａｔｉｃｉｔｙ　ａｎｄ　ａｎｔｉａｒｏｍａｔｉｃｉｔｙ．　Ｒｅｃｅｉｘｅｄ：Ａｕｇｕｓｔ　３Ｉ　！（）Ｉ　７：Ｒｅｖｉｓｅｄ：Ｏｃｔｏｂｅｒ１　８　２０Ｉ　７：Ａｃｃｅｐｔｅｄ：Ｏｃｔｏｂｅｒ１８　２０Ｉ７：Ｐｕｂｌｉｓｈｅｄ　ｏｎＩｉｌｌｅ：Ｏｃｔｏｂｅｒ　２３，２（）Ｉ　７　‘Ｃｏｎ’ｅｓｐｏｎｄｉｎｇ　ａｕｔｈｏｒｓ　Ｅｍａｉｌ：ｒ（】ｎｇｃｈｕｎｙｊｎｇ＠ａｌｉｙＬｍ　ＣＯＩｌｌ；Ｔｅ１．　８６—７３　：ｌ－８８８７２５３３（Ｃ　Ｙ　Ｒ）．Ｅｍａｉｌ：ｆｄｅｐｒｏｆ（ａ）ｖｕｂ　ａｃ　ｂｅ；Ｔｅｌ：１　３２（０）２６２９３３　１　０　Ｄ　Ｉ】ｌ１，　ＥｍａｉＩ：ｓｈｕｂｉｎ（ａ　ｅｍａｉＩ　Ｉ．ｉｌｌＣ　ｅｄｕ：Ｔｅｌ：＋１．９Ｉ９—９６２—４０３２　ｒＳ　Ｂ　Ｌ¨１　ＣＹＲ　ａｎｄ　ＳＢＬ　ａｃｋｎｏｗｌｅｄｇｅ　ｓｕｐｐｏｒｔ　ｆｒｏｍ　ｔｈｅ　ＮａｔｉｏｎａＩ　Ｎａｔｕｒａｌ　Ｓｃｉｅｎｃｅ　Ｉ：ｏｕｎｄａｔｉｏｎ　ｏｆ、Ｃｈｉ　ｒｌａ　ｆ２１　５０３０７６ｌ　ａｎｄ｝Ｉｕｕａｎ　Ｐｒｏｖｉｎｃｉａ１　Ｎａｔｕｒａｌ　Ｓｃｉｃｎｃｃ　ｌ￣＇ｏｔｌｎｄａｔｉｏｌｌ　ｏｆ　Ｃｈｉｎａ（２｛）１　７，１．Ｊ３２（）１）Ｉ）ＨＹ　ａｃｋｎｏｗｌｅｄｇｅｓ　ｔｈｅ　ｓｕｐｐｏｒｔ　ｌ＂ｔｏｍ　ｔｈｅ　ｔｔｕｒｌａｌｌ　ＰｒｏｖｉｎｃｉａＩ　ｌｎｎｏｖａｔｉｏｕ　Ｆｏｕｎｄａｔｊｏｎ　ｌｂｒ　Ｐｏｓｔｇｒａｄｕａｔｅｓ（ＣＸ２Ｉ１ｌ　７Ｂ１　７９）　ａｎｄ【’ｈｉｎａ　Ｓｃｈｏｌａｒｓｈｉｐ　Ｃｏｍｌｃｉｌ（２０Ｉ　７０６７２００Ｉ　５）．ＦＤＰ　ａｃｋｎｏｗｌｅｄｇｅｓ　ｔｈｅ　Ｒｅｓｅａｒｃｈ　Ｆｏｕｎｄａｔｉｏｎ　Ｆｌａｎｄｅｒｓ（ＦＷＯ）　ｒ　ＣＯｌｌｔｉｎＵＯｋｌＳ　ｓｕｐｐｏｒｔ　ｔｏ　ｈｉｓ　ｇｒｏｕｐ　ａｌ１【】ｔｈｅ　Ｖｒｉｉｅ　Ｕｎｉｖｅｒｓｉｔｃｉｔ　Ｂｒｕｓｓｅｌ　ｆｏｒ　ｓｕｐｐｏｒｔ，ａｍｏｌｌｇ　ｏｔｈｅｒ　ｔｈｒｏｕｇｈ　ａ　Ｓｔｒａｔｅｇｉｃ　Ｒｅｓｅａｒｃｈ　Ｐｒｏｇｒａｍ　ａｗａｒｄｅｄ　ｔｏ　ｈｉｓ　ｇｒｏｕｐ　Ｆｉｎａｌｌｙ．ＦＤＰ　ａｃｋｎｏｗｌｅｄｇｅｓ　ｔｈｅ　Ｆｒａｎｃｑｕｉ　ｆｏｕｎｄａｔｉｏｎ　ｆｂｒ　ａ　ｐｏｓｉｔｉｏｎ　ａｓ　Ｆｒａｎｃｑｕｉ　Ｒｅｓｅａｒｃｈ　Ｐｒｏｆｅｓｓｏｒ　Ｅｄｉｔｏｒｉａ１　ｏｌｔ５ｃｅ　ｏｆ　Ａｅｔａ　Ｐｈｙｓｉｃｏ—Ｃｈｉｍｉｃａ　Ｓｉｎｉｃａ　６４０　Ａｃｔａ　Ｐｈｙｓｉｃｏ—Ｃｈｉｍｉｃａ　Ｓｉｎｉｃａ　Ｖｏ１．３４　Ｋｅｙ　Ｗｏｒｄｓ：Ａｒｏｍａｔｉｃｉｔｙ；Ａｎｔｉａｒｏｍａｔｉｃｉｔｙ；Ｂｅｎｚｏｆｕｌｖｅｎｅ；Ｉｎｆｏｒｍａｔｉｏｎ—ｔｈｅｏｒｅｔｉｃ　ａｐｐｒｏａｃｈ；Ｄｅｎｓｉｔｙ　ｆｕｎｃｔｉｏｎａ　ｒｅａｃｔｉｖｉｔｙ　ｔｈｅｏｒｙ　运用密度泛函活性理论的信息论方法研究苯并富烯衍生物的芳香性　余东海’，荣春英’，，卢天２，ＤＥ　ＰＲＯＦＴ　Ｆｒａｎｋ３，，刘述斌４，　ｉ湖南师范大学化学化工学院教育部化学生物和中药研究重点实验室，长沙４１００８１　２北京科音自然科学研究中心，北京１０００２２　３布鲁塞尔自由大学普通化学研究组，布鲁塞尔普莱园赖２号，１０５０，比利时　北卡罗莱纳大学大学教堂山校区研究计算中心，北卡２７５９９．３４２０，美国　摘要：文献中已有越来越多的芳香性体系被发现，同时也有越来越多的芳香性指标被提出来，但是如何解释芳香化合物　稳定性的起源以及理解芳香性的本质仍然是当今理论化学中一个悬而未决的难题。运用我们新近提出的密度泛函活性理　论信息论方法，不久前我们曾对一系列富烯衍生物进行了系统研究并得到了一个全新的认识。本文进一步探讨苯并富烯　衍生物的芳香性行为，目的在于考察一个或多个苯环与富烯连接之后其芳香性发生变化的情况。运用香农熵，费舍尔信　息，Ｇｈｏｓｈ．Ｂｅｒｋｏｗｉｔｚ．Ｐａｒｒ熵，Ｏｎｉｃｅｓｃｕ信息能，信息增益，以及相对Ｒ￣ｎｙｉ熵六个信息量，和四种芳香指标，ＡＳＥ，　ＨＯＭＡ，ＦＬＵ￣ＩＩＮＩＣＳ，我们系统地研究了信息量和芳香性指标在单、双、三苯并富烯衍生物中的相关性。我们发现，不　管是否有苯环与富烯相连，芳香指标和信息量的交叉相关性都是一样的。这表明，虽然苯环本身具有芳香性，但苯环与　富烯相连并不能改变富烯的芳香性与反芳香性本质。苯并富烯衍生物与富烯衍生物的芳香性和反芳香性一致。苯并富烯　衍生物的芳香性和反芳香性完全取决于富烯本身的芳香性和反芳香性。这些结果为认识和理解复杂体系芳香性和反芳香　性起源和本质将提供有益的启示。　关键词：芳香性：反芳香性；苯并富烯；信息论方法；密度泛函活性理论　中图分类号：Ｏ６４１　１　Ｉｎｔｒｏｄｕｃｔｉｏｎ　ｉｎｄｅｐｃｎｄｅｎｔ　ｃｈｅｍｉｃａｌ　ｓｈｉｔｆｓ），ｅｔｃ．１，２８　３６　Ｔｈｅｓｅ　ｄｅｓｃｒｉｐｔｏｒｓ　ａｒｅ　ｏｎｌｙ　ｉｎｄｉｃａｔｉｖｅ　ｏｆ　ｏｎｅ　ｏｒ　ｆｅｗ　ｆａｃｅｔｓ　ｏｆ　ｔｈｅ　ｃｏｍｐｌｉｃａｔｅｄ　ｎａｔｕｒｅ　ｏｆ　ｒｏｍａｔｉａｃｉｔｙ．Ｎｏ　ｓｉｎｇｌｅ　ｉｎｄｅｘ　ｉｓ　ａｂｌｅ　ｔｏ　ａｃｃｏｕｎｔ　ｆｏｒ　ａｌｌ　ｋｉｎｄｓ　Ａｒｏｍａｔｉｃｉｔｙ　ｉｓ　ｏｎｅ　ｏｆ　ｔｈｅ　ｍｏｓｔ　ｗｉｄｅｌｙ　ｕｓｅｄ　ｃｈｅｍｉｃａｌ　ｃｏｎｃｅｐｔｓ　１．Ｉｔ　ｒｅｆｅｒｓ　ｔｏ　ｔｈｅ　ｅｘｔｒａ　ｓｔａｂｉｌｉｔｙ　ａｎｄ　ｓｕｂｓｉｄｅｄ　ｒｅａｃｔｉｖｉｔｙ　ｏｆ　ａ　ｐｌａｎａｒ　ｃｙｃｌｉｃ　ｓｔｒｕｃｔｕｒｅ　ｃｏｍｐａｒｅｄ　ｔｏ　ｏｔｈｅｒ　ｃｏｎｆｏｒｍａｔｉｏｎｓ　ｏｆ　ｔｈｅ　ｓａｍｅ　ｍｏｌｅｃｕｌａｒ　ｆｏｒｍｕｌａ　２－５．Ｉｔ　ｉｓ　ｏｒｉｇｉｎａｔｅｄ　ｆｒｏｍ　ｔｈｅ　ａｄｄｉｔｉｏｎａｌ　ｄｅｌｏｃａｌｉｚａｔｉｏｎ　ｏｆ　ｅｌｅｃｔｒｏｎｓ，ｅｉｔｈｅｒ瓦　ｏｒ　，ｌｅａｄｉｎｇ　ｏｆ　ａｒｏｍａｔｉｃｉｔｙ　ｍｅｎｔｉｏｎｅｄ　ａｂｏｖｅ　８，３７．　Ｒｅｃｅｎｔｌｙ，ｗｅ　ｔａｃｋｌｅｄ　ｔｈｉｓ　ｍａｔｔｅｒ　ｆｒｏｍ　ａ　ｃｏｍｐｌｅｔｅｌｙ　ｄｉｆｆｅｒｅｎｔ　ｐｅｒｓｐｅｃｔｉｖｅ　３８．Ｕｓｉｎｇ　ｔｈｅ　ｑｕａｎｔｉｔｉｅｓ　ｗｅ　ｎｅｗｌｙ　ｄｅｖｅｌｏｐｅｄ　ｗｉｔｈｉｎ　ｔｏ　ｔｈｅ　ｒｅｄｉｓｔｒｉｂｕｔｉｏｎ　ｏｆｔｈｅ　ｅｌｅｃｔｒｏｎ　ｄｅｎｓｉｔｙ　ｏｆｔｈｅ　ｅｎｔｉｒｅ　ｓｙｓｔｅｍ　６－９．　ｔｈｅ　ｆｒａｍｅｗｏｒｋ　ｏｆ　ｄｅｎｓｉｔｙ　ｆｕｎｃｔｉｏｎａｌ　ｒｅａｃｔｉｖｉｔｙ　ｔｈｅｏｒｙ（ＤＦＲＴ）３９，４０，　ｗｉｔｈ　ｓｕｂｓｔｉｔｕｔｅｄ　ｆｕｌｖｅｎｅ　ｄｅｒｉｖａｔｉｖｅｓ　ａｓ　ｔｈｅ　ｍｏｄｅｌ　ｓｙｓｔｅｍ，ｗｅ　Ｄｅｓｐｉｔｅ　ｏｆ　ｉｔｓ　ｇｅｎｅｒａｌ　ａｃｃｅｐｔａｎｃｅ　ａｎｄ　ｗｉｄｅｓｐｒｅａｄ　ｕｓｅ，ｔｈｅｒｅ　ａｒｅ　ｓｔｉｌｌ　ｌｏｔｓ　ｏｆ　ｃｏｎｔｒｏｖｅｒｓｉｅｓ　ａｂｏｕｔ　ａｒｏｍａｔｉｃｉｔｙ　ａｎｄ　ｉｔｓ　ｃｏｕｎｔｅｒｐａｒｔ．　ａｎｔｉａｒｏｍａｔｉｃｉｔｙ，ｉｎ　ｔｈｅ　ｌｉｔｅｒａｔｕｒｅ　９－１２．Ｔｈｉｓ　ｉｓ　ｍａｉｎｌｙ　ｃａｕｓｅｄ　ｂｙ　ｔｈｅ　ｆａｃｔ　ｔｈａｔ　ｔｈｅｒｅ　ａｒｅ　ｍａｎｙ　ｄｉｆｆｅｒｅｎｔ　ｃａｔｅｇｏｒｉｅｓ　ｏｆ　ａｒｏｍａｔｉｃｉｔｙ　ｄｉｓｃｏｖｅｒｅｄ　ｉｎ　ｔｈｅ　ｌｉｔｅｒａｔｕｒｅ，ｓｕｃｈ　ａｓ　Ｈｆｉｃｋｅｌ，Ｍ６ｂｉｕｓ，ｅｘｃｉｔｅｄ　ｓｔａｔｅ，ｓｐｈｅｒｉｃａｌ，ｃｕｂｉｃ，ｏｃｔａｈｅｄｒａｌ，ｍｅｔａｌｌｉｃ，ｃｈｅｌａｔｅｄ，ａｎｄ　ｓｏ　ｃｏｍｐａｒｅｄ　ｃｈａｎｇｉｎｇ　ｐａｔｔｅｎｓｒ　ｏｆ　ｉｎｆｏｒｍａｔｉｏｎ—ｔｈｅｏｒｅｔｉｃ（ＩＴ）　ｑｕａｎｔｉｔｉｅｓ　ｓｕｃｈ　ａｓ　Ｓｈａｎｎｏｎ　ｅｎｔｒｏｐｙ，Ｆｉｓｈｅｒ　ｉｎｆｏｒｍａｔｉｏｎ，　Ｇｈｏｓｈ—Ｂｅｒｋｏｗｉｔｚ·Ｐａｒｒ　ｅｎｔｒｏｐｙ，ｉｎｆｏｒｍａｔｉｏｎ　ｇａｉｎ，Ｏｎｉｃｅｓｃｕ　ｉｎｆｏｒｍａｔｉｏｎ　ｅｎｅｒｇｙ，ａｎｄ　ｒｅｌａｔｉｖｅ　Ｒ６ｎｙｉ　ｅｎｔｒｏｐｙ　ｏｎ　ｔｈｅ　ｒｉｎｇ　ｃａｒｂｏｎ　ａｔｏｍｓ　ｏｆ　ｔｈｅｓｅ　ｓｙｓｔｅｍｓ　ｗｉｔｈ　ｆｏｕｒ　ｒｅｐｒｅｓｅｎｔａｔｉｖｅ　ｏ１３　．Ａｎｄ　ａｔ　ｔｈｅ　ｓａｍｅ　ｔｉｍｅ　ｔｈｅｒｅ　ａｒｅ　ｄｏｚｅｎｓ　ｏｆ　ａｒｏｍａｔｉｃｉｔｙ　ｒｏｍａｔｉａｃｉｔｙ　ｉｎｄｅｘｅｓ，ＦＬＵ，ＨＯＭＡ，ＡＳＥ，ａｎｄ　ＮＩＣＳ　３８．Ｏｕｒ　ｒｅｓｕｌｔｓ　ｕｎｖｅｉｌｅｄ　ｔｗｏ　ｏｐｐｏｓｉｔｅ　ｐａｔｔｅｒｎｓ　ｆｏｒ　ｃｒｏｓｓ　ｃｏｒｒｅｌａｔｉｏｎｓ　ｂｅｔｗｅｅｎ　ＩＴ　ｑｕａｎｔｉｔｉｅｓ　ａｎｄ　ａｒｏｍａｔｉｃｉｔｙ　ｉｎｄｅｘｅｓ　ｆｏｒ　ｄｉｆｆｅｒｅｎｔ　ｉｎｄｅｘｅｓ　ｐｒｏｐｏｓｅｄ，ｗｈｏｓｅ　ａｃｒｏｎｙｍｓ　ｉｎｃｌｕｄｅ　ＡＳＥ／ＩＳＥ　（ａｒｏｍａｔｉｃ／ｉｓｏｍｅｒｉｃ　ｓｔａｂｉｌｉｚａｔｉｏｎ　ｅｎｅｒｇｙ），　ＲＥ（ｒｅｓｏｎａｎｃｅ　ｅｎｅｒｇｙ），Ａｉ（Ｊｕｌｇ　ｃｏｎｃｅｐｔ），ＨＯＭＡ（ｈａｒｍｏｎｉｃ　ｏｓｃｉｌｌａｔｏｒ　ｍｏｄｅｌ　ｒｉｎｇ　ｓｔｒｕｃｔｕｒｅｓ，　ａｎｄ　ｔｈｅｓｅ　ｒｉｎｇ—ｓｔｒｕｃｔｕｒｅ—ｄｅｐｅｎｄｅｎｔ　ｃｏｒｒｅｌａｔｉｏｎｓ　ａｒｅ　ｉｎ　ｇｏｏｄ　ａｇｒｅｅｍｅｎｔ　ｗｉｔｈ　ＨＯｃｋｅｌ’ｓ　４ｎ＋２　ｒｕｌｅ　ｏｆ　ａｒｏｍａｔｉｃｉｔｙ　ａｎｄ　４ｎ　ｒｕｌｅ　ｏｆ　ａｎｔｉａｒｏｍａｔｉｃｉｔｙ　ｌ，４２．　Ｉｎ　ｔｈｉｓ　ｗｏｒｋ，ｗｅ　ｃｏｎｔｉｎｕｅ　ｔｈｅ　ｓｔｕｄｙ　ｕｓｉｎｇ　ａ　ｎｅｗ　ｍｏｄｅｌ　ｓｙｓｔｅｍ．Ｗｅ　ｆｕｓｅ　ｏｎｅ　ｏｒ　ｍｏｒｅ　ｂｅｎｚｅｎｅ　ｒｉｎｇｓ．ｗｈｉｃｈ　ａｒｅ　ｋｎｏｗｎ　ｔｏ　ｂｅ　ａｒｏｍａｔｉｃ．ｗｉｔｈ　ｆｕｌｖｅｎｅ　ｄｅｒｉｖａｔｉｖｅｓ．Ｔｈｅ　ｐｕｒｐｏｓｅ　ｏｆ　ｔｈｅ　ｐｒｅｓｅｎｔ　ｓｔｕｄｙ　ｉｓ　ｔｏ　ｅｘａｍｉｎｅ　ｗｈｅｔｈｅｒ　ｏｒ　ｎｏｔ　ｔｈｅ　ｃｈａｎｇｉｎｇ　ｐａｔｔｅｒｎｓ　ｏｆ　ｃｒｏｓｓ　ｃｏｒｒｅｌａｔｉｏｎｓ　ｄｉｓｃｏｖｅｒｅｄ　ｐｒｅｖｉｏｕｓｌｙ　ａｒｅ　ｓｔｉｌ１　ｉｎ　ｏｆ　ａｒｏｍａｔｉｃｉｔｙ），Ｊｕｇ　ｉｎｄｅｘ，Ｂｉｒｄ　ｉｎｄｅｘ，ＰＤＩ（ｐａｒａ—ｒｅｌａｔｅｄ　ｄｅｌｏｃａｌｉｚａｔｉｏｎ　ｉｎｄｅｘ）．ＦＬＵ（ａｒｏｍａｔｉｃ　ｆｌｕｃｔｕａｔｉｏｎ　ｉｎｄｅｘ１，ＭＣＩ　（ｍｕｌｔｉｃｅｎｔｅｒ　ｉｎｄｉｃｅｓ），　ｎｇ（ｍｕｌｔｉｃｅｎｔｅｒ　ｉｎｄｅｘ），，ＮＧ（ｎｏｒｍａｌｉｚｅｄ　吨），／ｙＢ（ｎｏｒｍａｌｉｚｅｄ　ＭＣＩ），ＥＬＦ　（ｅｌｅｃｔｒｏｎ　ｌｏｃａｌｉｚａｔｉｏｎ　ｆｕｎｃｔｉｏｎ），ＡＴＩ（ａｖｅｒａｇｅ　ｔｗｏ—ｃｅｎｔｅｒ　ｉｎｄｉｃｅｓ），ＰＬＲ（ｐａｒａ　ｌｉｎｅａｒ　ｒｅｓｐｏｎｓｅ　ｉｎｄｅｘ），ＡＩＣＤ（ａｎｉｓｏｔｒｏｐｙ　ｏｆ　ｔｈｅ　ｉｎｄｕｃｅｄ　ｃｕｒｒｅｎｔ　ｄｅｎｓｉｔｙ），ＲＲＣＳ（ｄｅｎｓｉｔｙ　ｏｆ　ｒｉｎｇ　ｃｒｉｔｉｃａｌ　ｐｏｉｎｔ），ＮＩＣＳ（ｎｕｃｌｅｕｓ　Ｎｏ．６　ｄｏｉ：ｌ　０．３　８６６／ＰＫＵ．ＷＨＸＢ２０　１　７　１　０２３　１　６４１　ｐｌａｃｅ．Ｗｅ　ａｒｅ　ａｌｓｏ　ｉｎｔｅｒｅｓｔｅｄ　ｔｏ　ｓｅｅ　ｉｆ　ｔｈｅ　ａｒｏｍａｔｉｃｉｔｙ　ｏｆ　ｆｕｓｅｄ　ｂｅｎｚｅｎｅ　ｒｉｎｇｓ　ｉｓ　ｉｍｐａｃｔｅｄ　ｉｎ　ａｎｙ　ｍａｎｎｅｒ　ｂｙ　ｔｈｅ　ｆｕｌｖｅｎｅ　（Ａ，Ｂ）＝４∑　Ｓｉ，ｊ（Ａ）Ｓｉ，ｊ（Ｂ），ｗｈｅｒｅ　Ｓｉｊ（Ａ）ｉｓ　ｔｈｅ　ｏｖｅｒｌａｐ　ｏｆｔｈｅ　ｍｏｌｅｃｕｌａｒ　ｏｒｂｉｔａｌｓ　ａｎｄ，ｗｉｔｈｉｎ　ｔｈｅ　ｂａｓｉｎ　ｏｆ　ａｔｏｍ　Ａ．Ｔｈｅ　ｓｍａｌｌｅｒ　ｔｈｅ　ＦＬＵ　ｖａｌｕｅ，ｔｈｅ　ｓｔｒｏｎｇｅｒ　ｔｈｅ　ａｒｏｍａｔｉｃｉｔｙ　３５．　Ｔｈｅ　ｌａｓｔ　ｏｎｅ　ｅｍｐｌｏｙｅｄ　ｉｎ　ｏｕｒ　ｓｔｕｄｙ　ｉｓ　ｔｈｅ　ｎｕｃｌｅｕｓ．　ｄｅｒｉｖａｔｉｖｅｓ．Ｒｅｓｕｌｔｓ　ｏｂｔａｉｎｅｄ　ｆｒｏｍ　ｔｈｉｓ　ｓｔｕｄｙ　ｓｈｏｕｌｄ　ｐｒｏｖｉｄｅ　ｕｒｔｆｈｅｒ　ｕｎｄｅｒｓｔａｎｄｉｎｇ　ａｂｏｕｔ　ｔｈｅ　ｎａｔｕｒｅ　ａｎｄ　ｏｒｉｇｉｎ　ｏｆ　ａｒｏｍａｔｉｃｉｔｙ　ａｎｄ　ａｎｔｉａｒｏｍａｔｉｃｉｔｙ　ｆｏｒ　ｔｈｅ　ｓｙｓｔｅｍｓ　ｓｔｕｄｉｅｄ　１０，１９，４３，４４．　ｉｎｄｅｐｅｎｄｅｎｔ　ｃｈｅｍｉｃａｌ　ｓｈｉｔ（ｆＮＩＣＳ）ｄｅｒｉｖｅｄ　ｆｒｏｍ　ｔｈｅ　ｅ　ｂｃｔ　ｏｆ　ａｒｏｍａｔｉｃ　ｒｉｎｇ　ｃｕｒｒｅｎｔ．Ｉｔ　ｉｓ　ｆｏｕｎｄ　ｔｈａｔ　ａ　ｄｉａｔｒｏｐｉｃ（ｄｉａｍａｇｎｅｔｉｃ）　２　ＴｈｅｏｒｅｔｉｃａＩ　ｆｒａｍｅｗｏｒｋ　Ｔｈｅ　ｏｒｉｇｉｎａ１　ｄｅｆｉｎｉｔｉｏｎ　ｏｆ　ａｒｏｍａｔｉｃｉｔｙ　ｆｅａｔｕｒｅｓ　ｔｈｅ　ｅｘｔｒａ　ｒｉｎｇ　ｃｕｒｒｅｎｔ　ｉｓ　ａｓｓｏｃｉａｔｅｄ　ｗｉｔｈ　ａｒｏｍａｔｉｃｉｔｙ，ｗｈｅｒｅａｓ　ａ　ｐａｒａｔｒｏｐｉｃ　ｒｐａｒａｍａｇｎｅｔｉｃ）ｒｉｎｇ　ｃｕｒｒｅｎｔ　ｓｉｎａｌｇｓ　ａｎｔｉａｒｏｍａｔｉｃｉｔｙ．１１１ｉｓ　ｄｉｆｆｅｒｅｎｃｅ　ｉｎ　ｒｉｎｇ　ｃｕｒｒｅｎｔｓ　ｇｅｎｅｒａｔｅｓ　ｎｏｔｉｃｅａｂｌｅ　ｄｉｆｆｅｒｅｎｃｅｓ　ｉｎ　ＮＭＲ　ｃｈｅｍｉｃａｌ　ｓｈｉｆｔｓ．ａｎｄ　ｔｈｕｓ　ｃａｎ　ｂｅ　ｕｓｅｄ　ｔｏ　ｑｕａｎｔｉｆｙ　ｓｔａｂｉｌｉｔｙ　ｏｆ　ａ　ｐｌｎａａｒ　ｃｙｃｌｉｃ　ｓｔｒｕｃｔｕｒｅ　ｗｉｔｈ　ａ　ｃｏｎｊｕｇａｔｅｄ　ｓｙｓｔｅｍ　ｏｆ　４，２＋２　ｅｌｅｃｔｒｏｎｓ　４１．Ｔｏ　ｑｕａｎｔｉｆｙ　ａｒｏｍａｔｉｃｉｔｙ，ａ　ｎｕｍｂｅｒ　ｏｆ　ｄｅｓｃｒｉｐｔｏｒｓ　ｈａｖｅ　ｂｅｅｎ　ｐｒｏｐｏｓｅｄ　ｕｓｉｎｇ　ｄｉｆｅｒｅｎｔ　ｐｒｏｐｅｒｔｉｅｓ　ａｓ　ｔｈｅ　ｃｏｎｓｅｑｕｅｎｃｅｓ　ｏｆ　ａｒｏｍａｔｉｃｉｔｙ．Ｗｅ　ｓｅｌｅｃｔ　ｆｏｕｒ　ｒｅｐｒｅｓｅｎｔａｔｉｖｅ　ｄｅｓｃｒｉｐｔｏｒｓ　ｆｒｏｍ　ｔｈｅ　ｖｉｅｗｐｏｉｎｔ　ｏｆ　ｔｈｅ　ｅｎｅｒｇｅｔｉｃｓ（ＡＳＥ）　，　一，　ｅｌｅｃｔｒｏｎ　ｄｅｌｏｃａｌｉｚａｔｉｏｎ（ＦＬＵ）６，３５，ｇｅｏｍｅｔｒｉｃ（ＨＯＭＡ）３３　ａｎｄ　ｍａｇｎｅｔｉｃ（ＮＩＣＳ）ｃｒｉｔｅｒｉａ　４５，ｉｎ　ｔｈｉｓ　ｓｕｔｄｙ．ＡＳＥ（ａｒｏｍａｔｉｃ　ｓｔａｂｉｌｉｚａｔｉｏｎ　ｅｎｅｒｇｉｅｓ）ｉｓ　ａｎ　ｅｎｅｒｇｅｔｉｃｓ　ｄｅｓｃｒｉｐｔｏｒ　ｄｅｆｉｎｅｄ　ａｓ　ｔｈｅ　ｔｏｔａｌ　ｅｎｅｒｇｙ　ｄｉｆｆｅｒｅｎｃｅ　ｏｆ　ａｎ　ｉｓｏｄｅｓｍｉｃ　ｏｒ　ｈｏｍｏｄｅｓｍｏｔｉｃ　ｒｅａｃｔｉｏｎ　ｏｒ　ｉｓｏｍｅｒｉｚａｔｉｏｎ　ｂｅｔｗｅｅｎ　ａｎ　ａｒｏｍａｔｉｃ　ｓｔｒｕｃｔｕｒｅ　ａｎｄ　ｉｔｓ　ｏｔｈｅｒ　ｓｔｒｕｃｔｕｒｅｓ　ｗｉｔｈ　ｔｈｅ　Ｓ８ｉＴＩｅ　ｃｈｅｍｉｃａｌ　ｆｏｒｍｕｌａ．Ｔａｋｅ　ｂｅｎｚｅｎｅ　ｏｒ　ｔｏｌｕｅｎｅ　ａｓ　ａｎ　ｅｘａｍｐｌｅ，ｔｈｅ　ｓｉｍｐｌｅｓｔ　ｉｓｏｄｅｓｍｉｃ，　ｈｏｍｏｄｅｓｍｏｔｉｃ，ａｎｄ　ｉｓｏｍｅｒｉｚａｔｉｏｎ　ｒｅａｃｔｉｏｎｓ　ｒｅｐｒｅｓｅｎｔｉｎｇ　ＡＳＥ　ａｒｅ．ｒｅｓｐｅｃｔｉｖｅｌｙ　，５，４６，４７．　３Ｈ２Ｃ＝ＣＨＥ＋３Ｈ３Ｃ—ＣＨ３＝＝＞Ｃ６Ｈ６＋６ＣＨ４　ｆ１）　３ｃｉｓ—Ｈ２Ｃ＝ＣＨ．ＣＨ＝ＣＨ２＝＝＞Ｃ６Ｈ６＋３ＨＥＣ＝ＣＨＥ　（２）　Ａ　ｐｏｓｉｔｉｖｅ　ｖａｌｕｅ　ｏｆ　ＡＳＥ　ｄｅｎｏｔｅｓ　ｔｈａｔ　ｂｅｎｚｅｎｅ　１ｓ　ｍｏｒｅ　ｓｔａｂｌｅ　ｔｈａｎ　ｉｔｓ　ｒｅｆｅｒｅｎｃｅ，ｅｔｈｙｌｅｎｅ　ｏｒ　１，３－ｂｕｔａｄｉｅｎｅ，ａｎｄ　ｔｈｕｓ　ｔｈｅ　ｅｘｉｓｔｅｎｃｅ　ｏｆ　ａｒｏｍａｔｉｃｉｔｙ　ｓｈｏｕｌｄ　ｂｅ　ｉｎ　ｐｌａｃｅ，ｗｈｅｒｅａｓ　ａ　ｎｅｇａｔｉｖｅ　ＡＳＥ　ｖａｌｕｅ　ｓｕｇｇｅｓｔｓ　ｔｈａｔ　ｔｈｅ　ｓｙｓｔｅｍ　ｉｓ　ｏｆ　ａｎｔｉａｒｏｍａｔｉｃ　ｎａｔｕｒｅ　４８．　Ｔｈｅ　ＨＯＭＡ（ｈａｒｍｏｎｉｃ　ｏｓｃｉｌｌａｔｏｒ　ｍｏｄｅｌ　ｏｆ　ａｒｏｍａｔｉｃｉｔｙ）　ｉｎｄｅｘ　ｉｓ　ａ　ｇｅｏｍｅｔｒｉｃａｌ　ｍｅａｓｕｒｅｍｅｎｔ　ｏｆ　ｔｈｅ　ｅｑｕａｌｉｚａｔｉｏｎ　ｏｆ　ｃｈｅｍｉｃａｌ　ｂｏｎｄｓ　ｏｎ　ａｎ　ａｒｏｍａｔｉｃ　ｒｉｎｇ．Ｉｔｓ　ｄｅｆｉｎｉｔｉｏｎ　ｉｓ　ｔｈｅ　ｆｏｌｌｏｗｉｎｇ　ＨＯＭＡ＝卜　∑（‰一Ｒｉ）。＝卜［ａ（Ｒｏｐｔ－Ｒａ　）　＋　∑（尼　一Ｒｉ）　］　ｆ４）　，ｚ　胛　、　ｗｈｅｒｅ／，／ｉｓ　ｔｈｅ　ｂｏｎｄ　ｎｕｍｂｅｒ　ｏｆｔｈｅ　ｃｏｎｓｉｄｅｒｅｄ　ｒｉｎｇ，Ｒａｖ　ｉｓ　ｔｈｅ　ａｖｅｒａｇｅ　ｂｏｎｄ　ｌｅｎｇｔｈ，Ｒｉ　ｉｓ　ｈｔｅ　ｂｏｎｄ　ｌｅｎｇｔｈ（ａｌｌ　ｉｎ　Ａ，１　Ａ＝０．１　ｎｍ），ｏｆ　ｅａｃｈ　ｂｏｎｄ　ｏｎ　ｔｈｅ　ｒｉｎｇ，ａｎｄ　ｉｓ　ｔｈｅ　ｎｏｒｍａｌｉｚａｔｉｏｎ　ｃｏｎｓｔａｎｔ（ｆｏｒ　Ｃ—Ｃ　ｂｏｎｄｓ　＝２５７．７、．Ａ　ＨｏＭＡ　ｖａｌｕｅ　ｏｆ＜　０，＝０，ａｎｄ＝１　ｒｅｐｒｅｓｅｎｔｓ　ａｎｔｉｒａｏｍａｔｉｃ，ｎｏｎａｒｏｍａｔｉｃ，ａｎｄ　ｐｅｒｆｅｃｔ　ａｒｏｍａｔｉｃ　ｓｙｓｔｅｍｓ，ｒｅｓｐｅｃｔｉｖｅｌｙ　３３，４９　Ａｓ　ａｎ　ｅｌｅｃｔｒｏｎｉｃ　ｄｅｓｃｒｉｐｔｏｒ，ｔｈｅ　ａｒｏｍａｔｉｃ　ｆｌｕｃｔｕａｔｉｏｎ　ｉｎｄｅｘ　（ＦＬＵ）ｄｅｓｃｒｉｂｅｓ　ｔｈｅ　ｆｌｕｃｔｕａｔｉｏｎ　ｏｆ　ｅｌｅｃｔｒｏｎｉｃ　ｃｈａｒｇｅ　ｂｅｔｗｅｅｎ　ａｄｊａｃｅｎｔ　ａｔｏｍｓ　ｉｎ　ａ　ｇｉｖｅｎ　ｒｉｎｇ．Ｉｔｓ　ｄｅｆｉｎｉｔｉｏｎ　ｉｓ　ａｓ　ｆｏｌｌｏｗｓ，　Ｆｌｕ（Ａ　－－￣Ｂ）Ｂ叫Ｊ］　［　ｃ　（Ａ，Ｂ）　２［Ⅳ（Ａ）一　（Ａ）】　：ｊ１　（Ａ　Ｂ）＞Ｆ１　（Ｂ　，ａｎ—１　ＦＩｕｆＡ—Ｂ）＜Ｆ１ｕ（Ｂ－－＊Ａ），　ｄ　ａｒｏｍａｔｉｃｉｔｙ　ａｎｄ　ａｎｔｉａｒｏｍａｔｉｃｉｔｙ．Ｉｎ　ｆｏｒｍｕｌａｔｉｏｎ．ＮＩＣＳ　ａｔ　ｔｈｅ　ｃｈｏｓｅｎｐｏｉｎｔＲＮｌｏｃａｔｅｄ　ａｔ　ｏｒｏｎｔｏｐ　ｏｆａ　ｒｉｎｇ　ｃａｎｂｅ　ｄｅｓｃｒｉｂｅｄ　ａｓ　ｔｈｅ　ｓｕｍ　ｏｆ　ｐａｒｔｉａｌ　ｃｈｅｍｉｃａｌ　ｓｈｉｔｆｓ　ａｒｉｓｉｎｇ　ｆｒｏｍ　ｏｃｃｕｐｉｅｄ　ｍｏｌｅｃｕｌａｒ　ｏｒｂｉｔａｌｓ　（ｏ　一昙击　（　０　ＪＩ　Ｊ笮　　㈤　／　（６）　ｗｈｅｒｅ　ｔｈｅ　ｆｉｒｓｔ　ａｎｄ　ｓｅｃｏｎｄ　ｔｅｒｍｓ　ｏｎ　ｔｈｅ　ｒｉｇｈｔ—ｈａｎｄ　ｓｉｄｅ　ｏｆ　Ｅｑ．　（１）ａｒｅ　ｔｈｅ　ｄｉａｍａｇｎｅｔｉｃ　ａｎｄ　ｐａｒａｍａｇｎｅｔｉｃ　ｃｏｎｔｒｉｂｕｔｉｏｎｓ，　ｒｅｓｐｅｃｔｉｖｅｌｙ，ＬＮ＝ＦＮ×　ｒｅｆｅｒｓ　ｔｏ　ｔｈｅ　ａｎｇｕｌａｒ　ｍｏｍｅｎｔｕｍ　ｏｐｅｒａｔｏｒ．ｒＮ：　—ＲＮ，ａｎｄ，ｉｓ　ｔｈｅ　ｅｌｅｃｔｒｏｎ　ｐｏｓｉｔｉｏｎ．Ａ　ｍｏｒｅ　ｎｅｇａｔｉｖｅ　ｖａｌｕｅ　ｏｆ　ｏ　ｉｓ　ａｎ　ｉｎｄｉｃａｔｉｏｎ　ｏｆ　ａ　ｓｔｒｏｎｇｅｒ　ａｒｏｍａｔｉｃｉｔｙ，　ｗｈｅｒｅａｓ　ａ　ｐｏｓｉｔｉｖｅ　ｏ－ｖａｌｕｅ　ｓｕｇｇｅｓｔｓ　ｔｈａｔ　ｔｈｅ　ｒｉｎｇ　ｉｓ　ａｎｔｉａｒｏｍａｔｉｃ　４，４５，５０－５３．　Ｏｕｒ　ｒｅｃｅｎｔ　ｔｏｏｌ　ｔｏ　ｄｅａｌ　ｗｉｔｈ　ａｒｏｍａｔｉｃｉｔｙ　ａｎｄ　ａｎｔｉａｒｏｍａｔｉｃｉｔｙ　ｉｓ　ｔｈｅ　ｉｎｆｏｒｍａｔｉｏｎ－ｔｈｅｏｒｅｔｉｃ　ａｐｐｒｏａｃｈ　ｆｒｏｍ　ＤＦＲＴ　３９　ｗｈｅｒｅ　ｗｅ　ｍａｋｅ　ｕｓｅ　ｏｆ　ｓｉｍｐｌｅ　ｄｅｎｓｉｔｙ　ｆｕｎｃｔｉｏｎａｌｓ　ｔｏ　ａｐｐｒｅｃｉａｔｅ　ｐｈｙｓｉｏｃｈｅｍｉｃａｌ　ｐｒｏｐｅｒｔｉｅｓ　ｌｉｋｅ　ａｒｏｍａｔｉｃｉｔｙ，ｓｔｅｒｉｃ　ｅｆｆｅｃｔ，　ｅｌｅｃｔｒｏｐｈｉｌｉｃｉｔｙ，　ｎｕｃｌｅｏｐｈｉｌｉｃｉｔｙ，　ｎａｄ　ｓｔｅｒｅｏｓｅｌｅｃｔｉｖｉｔｙ．　Ｑｕａｎｔｉｔｉｅｓ　ｃｏｎｓｉｄｅｒｅｄ　ｉｎｃｌｕｄｅ　Ｓｈａｎｎｏｎ　ｅｎｔｒｏｐｙ　ｓｓ　５４　＝一ｆ　ｐ（ｒ）ｌｎｐ（ｒ）ｄｒ　（７）　Ｆｉｓｈｅｒ　ｉｎｆｏｒｍａｔｉｏｎ／ｅ　５４　＝　（８）　Ｏｎｉｃｅｓｃｕ　ｉｎｆｏｒｍａｔｉｏｎ　ｅｎｅｒｇｙ　ｏｆ　ｏｒｄｅｒ　５５　＝　ｐ　（９）　Ｇｈｏｓｈ—Ｂｅｒｋｏｗｉｔｚ—Ｐａｒ（ＧＢＰ）ｅｎｔｒｏｐｙ　ＳＧＢＰ　＝　卜　ｉｎｆｏｒｍａｔｉｏｎ　ｇａｉｎ（ａｌｓｏ　ｃａｌｌｅｄ　Ｋｕｌｌｂａｃｋ－Ｌｅｉｂｌｅｒ　ｄｉｖｅｒｇｅｎｃｅ）　ＩＧ　３９　，Ｇ＝　）ｌｎ　出　…）　ａｎｄ　ｒｅｌａｔｉｖｅ　Ｒｒｎｙｉ　ｅｎｔｒｏｐｙ　Ｒｒ　ｏｆ　ｏｒｄｅｒ　ｎ　５５，　＝　－ｎ［　ｄ，１　Ｉｎ　ａｂｏｖｅ　ｆｏｒｍｕｌａｓ，ｐ（ｒ）ｉｓ　ｔｈｅ　ｅｌｅｃｔｒｏｎ　ｄｅｎｓｉｔｙ，Ｖｐ（ｒ）ｉｓ　ｔｈｅ　ｄｅｎｓｉｔｙ　ｇｒａｄｉｅｎｔ，ｐｏ（ｒ）ｉｓ　ｔｈｅ　ｒｅｆｅｒｅｎｃｅ　ｓｔａｔｅ　ｄｅｎｓｉｔｙ　ｓａｔｉｓｌｙｉｎｇ　ｔｈｅ　ｓａｍｅ　ｎｏｒｍａｌｉｚａｔｉｏｎ　ｃｏｎｄｉｔｉｏｎ　ａｓ　ｐ（　，　，；Ｐ）ｉｓ　ｔｈｅ　ｋｉｎｅｔｉｃ　６４２　Ａｃｔａ　Ｐｈｙｓｉｃｏ—Ｃｈｉｍｉｃａ　Ｓｉｎｉｃａ　Ｖ０１．３４　＝＝＼　Ｒ　３Ｍ　４Ｍ一　４Ｍ　５Ｍ　６Ｍ＋　６Ｍ一　７Ｍ　ｂｅｎ３Ｍ　一ｂｅｎ４Ｍ—　＿ｂｅｎ４Ｍ＋　一ｂｅｎ５Ｍ　一ｂｅｎ６Ｍ＋　一ｂｅｎ６Ｍ—　一ｂｅｎ７Ｍ　一Ｒ　Ｒ　Ｒ　Ｒ　ｄｉｂｅｎ　５Ｍ　ｄｉｂｅｎ６Ｍ＋＿ｄｉｂｅｎ６Ｍ—＿ｄｉｂｅｎ７Ｍ　一ｔｒｉｂｅｎ５Ｍ　一ｔｒｉｂｅｎ６Ｍ＋一ｔｒｉｂｅｎ６Ｍ—　＿Ｓｃｈｅｍｅ　１　Ｓｕｂｓｔｉｔｕｔｅｄ　ｆｕｌｖｅｎｅ　ａｎｄ　ｂｅｎｚｅｎｅ　ｆｕｓｅｄ　ｆｕｌｖｅｎｅ　ｄｅｒｉｖａｔｉｖｅｓ　ｅｘａｍｉｎｅｄ　ｉｎ　ｔｈｉｓ　ｗｏｒｋ　ｅｎｅｒｇｙ　ｄｅｎｓｉｔｙ，ｔｒＦ（ｒ；Ｐ）　ｉｓ　ｔｈｅ　Ｔｈｏｍａｓ—Ｆｅｒｍｉ　ｋｉｎｅｔｉｃ　ｅｎｅｒｇｙ　ｄｅｎｓｉｔｙ，ａｎｄ　ｃ　ａｎｄ　ｋ　ａｒｅ　ｃｏｎｓｔａｎｔｓ．Ｔｈｅｓｅ　ｄｉｆｆｅｒｅｎｔ　ｔｈｅｏｒｅｔｉｃａｌ　ｑｕａｎｔｉｔｉｅｓ　ｈａｖｅ　ｂｅｅｎ　ｐｒｅｖｉｏｕｓｌｙ　ａｐｐｌｉｅｄ　ｔｏ　ａｐｐｒｅｃｉａｔｅ　ｎｕｍｅｒｏｕｓ　ｍｏｌｅｃｕｌａｒ　ｐｒｏｃｅｓｓｅｓ　ａｎｄ　ｐｈｙｓｉｏｃｈｅｍｉｃａｌ　ｐｒｏｐｅｒｔｉｅｓ．Ｉｎ　ｔｈｉｓ　ｗｏｒｋ，ｗｅ　ａｐｐｌｙ　ｔｈｅｍ　ｔｏ　ａ　ｎｅｗ　ｍｏｄｅｌ　ｓｙｓｔｅｍ　ａｎｄ　ｅｘａｍｉｎｅ　ｔｈｅｉｒ　ｃｏｒｒｅｌａｔｉｏｎｓ　ｗｉｔｈ　ａｒｏｍａｔｉｃｉｔｙ　ｄｅｓｃｒｉｐｔｏｒｓ　ｌｉｋｅ　ＮＩＣＳ，ＦＬＵ，　Ｈ０ＭＡ　ａｎｄ　ＡＳＥ．　ｓａｍｅ　１ｅｖｅ１　ｏｆ　ｔｈｅｏｒｙ．Ｔｏ　ｐｅｒｆｏｒｍ　ｔｈｅ　ａｔｏｍｉｃ　ｐａｒｔｉｔｉｏｎ．Ｂｅｃｋｅ’ｓ　ｕｚｚｙ　ａｔｏｍ　ａｐｐｒｏａｃｈ　６５，Ｂａｄｅｒ’ｆｓ　ｚｅｒｏ—ｆｌｕｘ　ａｔｏｍｓ　ｉｎ．ｍｏｌｅｃｕｌｅｓ　ｃｒｉｔｅｒｉｏｎ　６６ａｎｄ　Ｈｉｒｓｈｆｅｌｄ’ｓ　ｓｔｏｃｋｈｏｌｄｅｒ　ａｐｐｒｏａｃｈ　ｏｒｅ　ｐｏｓｓｉｂｌｅ　６７．　，Ａｓ　ｈａｖｅ　ｂｅｅｎ　ｄｅｍｏｎｓｔｒａｔｅｄ　ｅａｒｌｉｅｒ，ｔｈｅｓｅ　ｔｈｒｅｅ　ａｐｐｒｏａｃｈｅｓ　ｙｉｅｌｄ　ｑｕａｌｉｔａｔｉｖｅｌｙ　ｓｉｍｉｌｒ　ｒｅｓｕｌｔａｓ　５５，６８，６９Ｉｎ　ｔｈｉｓ　ｗｏｒｋ，ｗｅ　．ｃｈｏｏｓｅ　ｔｈｅ　Ｈｉｒｓｈｆｅｌｄ’ｓ　ｓｔｏｃｋｈｏｌｄｅｒ　ａｐｐｒｏａｃｈ　ｔｏ　ｐａｒｔｉｔｉｏｎ　ａｔｏｍｓ　ｉｎ　ｍｏｌｅｃｕｌｅｓ　ｔｏ　ｏｂｔａｉｎ　ａｔｏｍｉｃ　ｖａｌｕｅｓ　ｏｆ　ａｂｏｖｅ　ＩＴ　ｑｕａｎｔｉｔｉｅｓ　ｏｎ　ｃａｒｂｏｎｓ　ａｔｏｍｓ　ｉｎ　ｔｈｅ　ａｒｏｍａｔｉｃ／ａｎｔｉａｒｏｍａｔｉｃ　ｒｉｎｇ．Ｉｎ　ｔｈｉｓ　ｓｔｕｄｙ．　３　Ｃｏｍｐｕｔａｔｉｏｎａｌ　ｄｅｔａｉｌｓ　Ｓｃｈｅｍｅ　１　ｓｈｏｗｓ　ｔｈｅ　ｍｏｄｅｌ　ｓｙｓｔｅｍｓ　ｉｎｖｅｓｔｉｇａｔｅｄ　ｉｎ　ｔｈｉｓ　ｗｅ　ｃｈｏｓｅ　ＮＩＣＳ（１），ｗｈｉｃｈ　ｍｅａｎｓ　ｔｈａｔ　ｔｈｅ　ｃｈｏｓｅｎ　ｐｏｉｎｔ　ＲＮ　ｉＳ　１　Ａ　ａｂｏｖｅ　ｔｈｅ　ｃｅｎｔｅｒ　ｏｆ　ｔｈｅ　ａｒｏｍａｔｉｃ　ｒｉｎｇ　３０，３１．Ｔｈｅ　Ｎ］ＣＳ（１、　ａｒｏｍａｔｉｃｉｔｙ　ｉｎｄｅｘ　ｗａｓ　ｃａｌｃｕｌａｔｅｄ　ｕｓｉｎｇ　ｔｈｅ　ｐｒｏｃｅｄｕｒｅ　ｆｒｏｍ　ｔｈｅ　ｌｉｔｅｒａｔｕｒｅ　７，１４，５６．Ｔｈｅ　ＡＳＥ　ｉｎｄｅｘ　ｆｏｒ　ｆｕｓｅｄ　ｆｕｌｖｅｎｅ　ｄｅｒｉｖａｔｉｖｅｓ　ｆｏｒ　３ＭＲ一７ＭＲ　ｓｐｅｃｉｅｓ　ｉｎ　Ｓｃｈｅｍｅ　ｌ　ｗａｓ　ｃｏｍｐｕｔｅｄ　ａｓ　ｔｈｅ　ｔｏｔａｌ　ｅｎｅｒｇｙ　ｄｉｆｆｅｒｅｎｃｅ　ｆｒｏｍ　ｔｈｅ　ｉｓｏｍｅｒｉｚａｔｉｏｎ　ｒｅａｃｔｉｏｎ　ｓｈｏｗｎ　ｉｎ　Ｓｃｈｅｍｅ　２　ｗｉｔｈ　ｔｈｅ　ｓａｍｅ　１ｅｖｅｌ　ｏｆ　ｔｈｅｏｒｙ．Ｔｈｅ　ａｒｏｍａｔｉｃ　ｉｎｄｅｘ　０ｆ　ＦＬＵ　ａｎｄ　ＨＯＭＡ　ｗｅｒｅ　ｃａｌｃｕｌａｔｅｄ　ｕｓｉｎｇ　ＭｕｌｔｉＷＦＮ　３．３．９　ｗｏｒｋ．Ｔｒｉ（３ＭＲ），ｔｅｔｒａ（４ＭＲ＋ａｎｄ　４ＭＲ—ｗｉｔｈ　ｏｎｅ　ｐｏｓｉｔｉｖｅ　ａｎｄ　ｎｅｇａｔｉｖｅ　ｃｈａｒｇｅ，ｒｅｓｐｅｃｔｉｖｅｌｙ），ｐｅｎｔａ（５ＭＲ），ｈｅｘａ（６ＭＲ＋　ａｎｄ　６ＭＲ一，ｗｉｔｈ　ｏｎｅ　ｐｏｓｉｔｉｖｅ　ａｎｄ　ｎｅｇａｔｉｖｅ　ｃｈａｒｇｅ，　ｒｅｓｐｅｃｔｉｖｅｌｙ），ａｎｄ　ｈｅｐｔａ一（７ＭＲ）ｆｕｌｖｅｎｅ　ｄｅｒｉｖａｔｉｖｅｓ　ｗｅｒｅ　ｗｈａｔ　ｗｅ　ｃｏｎｓｉｄｅｒｅｄ　ｂｅｆｏｒｅ，ｓｅｒｖｉｎｇ　ａｓ　ｔｈｅ　ｒｅｆｅｒｅｎｃｅ　ａｎｄ　ｔｈｕｓ　ｔｈｅ　ｓｔａｒｔｉｎｇ　ｐｏｉｎｔ　ｏｆ　ｔｈｉｓ　ｗｏｒｋ．Ｔｈｅ　ｎｅｗ　ｓｙｓｔｅｍｓ　ａｒｅ　ｆｕｓｅｄ　ｆｕｌｖｅｎｅ　ｄｅｒｉｖａｔｉｖｅｓ　ｗｉｔｈ　ｏｎｅ，ｔｗｏ，ａｎｄ　ｔｈｒｅｅ　ｂｅｎｚｅｎｅ　ｒｉｎｇｓ，ｓｈｏｗｎ　ｉｎ　Ｓｃｈｅｍｅ　１．ｄｅｎｏｔｅｄ　ｂｙ　ｂｅｎｎＭ，ｄｉｒｉｂｅｎｎＭ，　ｂｅｎｎＭ，ａｎｄ　ｔｐｒｏｇｒａｍ　ａｔ　ＨＦ／６—３　１１＋Ｇ（　ｌｅｖｅｌ　ｏｆｔｈｅｏｒｙ　６１，７０￣Ｕｎｉｔｓ　ｆｏｒ　ＡＳＥ　ｉｓ　ｋＪ·ｍｏｌ一，ｆｏｒ　ＮＩＣＳ　ｉｓ　１０．ａｎｄ　ｆｏｒ　ＩＴ　ｑｕａｎｔｉｔｉｅｓ　ａｔｏｍｉｃ　ｕｎｉｔｓ．　ｒｅｓｐｅｃｔｉｖｅｌｙ．ｗｉｔｈ　＝３—７．Ｔｈｅ　ｒｉｎｇｓ　ｗｉｔｈｉｎ　ｔｈｅ　ｍｏｄｅｌｓ　ａｒｅ　ｍａｒｋｅｄ　ｂｙ　Ａ，Ｂ，Ｃ　ａｎｄ　Ｄ（ｓｅｅ　Ｓｃｈｅｍｅ　１），ｗｉｔｈ　Ｒｉｎｇ　Ｂ　ｓｔａｎｄｉｎｇ　ｆｏｒ　ｔｈｅ　ｆｕｌｖｅｎｅ　ｒｉｎｇ．１１１ｅ　ｓｕｂｓｔｉｔｕｔｉｎｇ　ｇｒｏｕｐｓ　Ｒ　ｗａｓ　ｃｈｏｓｅｎ　ｆｒｏｍ　ｔｈｅ　ｆｏｌｌｏｗｉｎｇ　ｇｒｏｕｐｓ，Ｒ＝Ｈ，ＣＨ３，ＣＣＨ，ＣＭｅ３，　４　Ｒｅｓｕｌｔｓ　ａｎｄ　ｄＩｓＣｕＳＳｉＯｎ　Ａｓ　ｉｌｌｕｓｔｒａｔｉｖｅ　ｅｘａｍｐｌｅｓ．Ｔａｂｌｅ　１　ｓｈｏｗｓ　ｔｈｅ　ｎｕｍｅｒｉｃａｌ　ｒｅｓｕｌｔｓ　ｏｆ　ｔｈｒｅｅ　ａｒｏｍａｔｉｃｉｔｙ　ｄｅｓｃｒｉｐｔｏｒｓ，ＦＬＵ，ＨＯＭＡ，ａｎｄ　ＮＩＣＳ（１），　ｏｒ　ｔｆｈｒｅｅ　ｂｅｎｚｅｎｅ　ｆｕｓｅｄ　ｆｕｌｖｅｎｅ　ｄｅｒｉｖａｔｉｖｅｓ　ｗｉｔｈ　２５　ｓｕｂｓｔｉｔｕｔｉｎｇ　ｇｒｏｕｐｓ．Ａｓ　ｃａｎ　ｂｅ　ｓｅｅｎ　ｆｒｏｍ　ｔｈｅ　ｔａｂｌｅ，ｓａｍｅ　ａｓ　ｗｈａｔ　ｗｅ　ｈａｖｅ　ＣＮ，Ｃ０ＮＨ２，ＣＯＣＨ３，ＣＦ３，ＣＨ２，ＣＣ一，ＣＯＯ一，Ｆ，Ｂ（ＯＨ）２，ＯＨ，　０ＣＨ３，０，ＮＨ２，Ｎ０２，ＮＯ，ＮＭｅ２，ＮＨ一，ＮＨ；，ＮＮ　，ａｎｄ　ＳｉＭｅ３　５６，”．ｗｉｔｈ　ｔｈｅ　ｃｒｉｔｅｒｉｏｎ　ｔｈａｔ　ｔｈｅ　ｏｐｔｉｍｉｚｅｄ　ｓｔｍｃｍｒｅ　ｏｂｓｅｒｖｅｄ　ｉｎ　ｆｕｌｖｅｎｅ　ｄｅｒｉｖａｔｉｖｅｓ；ｆｉ１　ｔｈｅｓｅ　ｓｐｅｃｉｅｓ　ｃｏｕｌｄ　ｂｅ　ｂｏｔｈ　ａｒｏｍａｔｉｃ　ａｎｄ　ａｎｔ　ａｒｏｍａｔｉｃ，ａｃｃｏｒｄｉｎｇ　ｔｏ　ｔｈｅ　ｄｅｆｉｎｉｔｉｏｎ　ｏｆ　ｅａｃｈ　ｓｈｏｕｌｄ　ｂｅ　ｐｌａｎａｒ．Ｍｏｄｅｌ　ｔｒｉｂｅｎ７Ｍ　ｉｓ　ｎｏｔ　ｐｌａｎａｒ　ｓｏ　ｉｔ　ｗｉｌｌ　ｎｏｔ　ｂｅ　ｃｏｎｓｉｄｅｒｅｄ　ｈｅｒｅ．　Ａｌｌ　ｃａｌｃｕｌａｔｉｏｎｓ　ｗｅｒｅ　ｐｅｒｆｏｒｍｅｄ　ａｔ　ｔｈｅ　ＤＦＴ　．　Ｂ３ＬＹＰ／６—３　１　１　Ｇ（ｄ，ｐ、　一　ｌｅｖｅｌ　ｏｆ　ｔｈｅｏｒｙ　ｕｓｉｎｇ　Ｇａｕｓｓｉａｎ　０９　ｐａｃｋａｇｅ　ｖｅｒｓｉｏｎ　Ｅ０１　６２　ｗｉｔｈ　ｔｈｅ　ｔｉｇｈｔ　ＳＣＦ　ｃｏｎｖｅｒｇｅｎｃｅ　ｃｒｉｔｅｒｉｏｎ　ａｎｄ　ｕｌｔｒａｆｉｎｅ　ｉｎｔｅｇｒａｔｉｏｎ　ｇｒｉｄｓ　６３．Ａ　ｓｉｎｇｌｅ—ｐｏｉｎｔ　ｒｅｑｕｅｎｃｙ　ｃａｌｃｕｌａｔｉｏｎ　ｗａｓ　ｆｏｌｌｆｏｗｅｄ　ｔｏ　ｅｎｓｕｒｅ　ｔｈａｔ　ｔｈｅ　ｆｉｎａｌ　ｓｔｒｕｃｔｕｒｅ　ｏｂｔａｉｎｅｄ　ｈａｓ　ｎｏ　ｉｍａｇｉｎａｒｙ　ｆｒｅｑｕｅｎｃｙ．Ｔｈｅ　ＭｕｌｔｉＷＦＮ　３．３．９　ｐｒｏｇｒａｍ　ｄｅｖｅｌｏｐｅｄ　ｂｙ　ｏｎｅ　ｏｆ　ｔｈｅ　ｐｒｅｓｅｎｔ　ａｕｔｈｏｒｓ　ｗａｓ　ｕｓｅｄ　ｔｏ　ｃａｌｃｕｌａｔｅ　ａｌｌ　ｉｎｆｏｒｍａｔｉｏｎ—ｔｈｅｏｒｅｔｉｃ　ｑｕａｎｔｉｔｉｅｓ　ｉｎｔｒｏｄｕｃｅｄ　ａｂｏｖｅ　ｗｉｔｈ　ｔｈｅ　ｃｈｅｃｋ　ｐｏｉｎｔ　ｆｉｌｅ　ｇｅｎｅｒａｔｅｄ　ｆｒｏｍ　ｔｈｅ　ａｂｏｖｅ　Ｇａｕｓｓｉａｎ　ｃａｌｃｕｌａｔｉｏｎｓ　ａｓ　ｔｈｅ　ｉｎｐｕｔ　ｆｉｌｅ　６４．Ｔｏ　ｏｂｔａｉｎ　ｔｈｅ　ｅｌｅｃｔｒｏｎ　ｄｅｎｓｉｔｙ　ｆｕｒ　ｔｈｅ　ｉｓｏｌａｔｅｄ　ｓｔａｔｅ，　ｗｅ　ｅｍｐｌｏｙｅｄ　ｔｈｅ　ｓｐｈｅｒｉｃａｌｌｙ．ａｖｅｒａｇｅｄ　ｅｌｅｃｔｒｏｎ　ｄｅｎｓｉｔｙ　ｏｆ　ｔｈｅ　ｎｅｕ￣ａｌ　ａｔｏｍ　ａｔ　ｔｈｅ　Ｒ一　Ｒ　Ｒ　Ｓｃｈｅｍｅ　２　Ｔｈｅ　ｉｓｏｍｅｒｉｚａｔｉｏｎ　ｒｅａｃｔｉｏｎｓ　ｕｓｅｄ　ｔｏ　ｃａｌｃｕｌａｔｅ　ｔｈｅ　ＡＳＥ　ｉｎｄｅｘ　ｆｏｒ　ｓｉｎｇｌｙ　ｆｕｓｅｄ　ｂｅｎｚｏｆｕｌｖｅｎｅｓ．　Ｎｏ．６　ｄｏｉ：１０．３８６６／ＰＫＵ．ＷＨＸＢ２０１７１０２３　１　６４３　ＴａｂＩｅ　Ｎｕｍｅｒｉｃａｌ　ｒｅｓｕｌｔｓ　ｏｆ　ＦＬＵ，ＨＯＭＡ，ａｎｄ　ＮＩＣＳ（１）ｆｏｒ　Ｒｉｎｇ　Ｂ　ｏｆ　ｂｅｎ＿５Ｍ，ｄｉｂｅｎ＿５Ｍ，ａｎｄ　ｔｒｉｂｅｎ＿ｂｅｎ　ａｓ　ｉｌｌｕｓｔｒａｔｉｖｅ　ｅｘａｍｐｌｅｓ．　ｏｆ　ｔｈｅｓｅ　ａｒｏｍａｔｉｃｉｔｙ　ｄｅｓｃｒｉｐｔｏｒｓ；（ｉｉ）ａｒｏｍａｔｉｃｉｔｙ　ａｎｄ　ａｎｔｉａｒｏｍａｔｉｃｉｔｙ　ｄｅｐｅｎｄ　ｏｎ　ｔｈｅ　ｎａｔｕｒｅ　ｏｆ　ｔｈｅ　ｓｕｂｓｔｉｔｕｔｉｎｇ　ｇｒｏｕｐ；　ｂｅ　ｎｏｔｉｃｅｄ：（ｉ）Ｔｈｅｓｅ　ｃｏｒｒｅｌａｔｉｏｎ　ｐａｔｔｅｒｎｓ　ａｒｅ　ｔｈｅ　ｓａｍｅ　ａｓ　ｔｈｏｓｅ　ｏｆ　ｕｌｆｖｅｎｅ　ｄｅｒｉｖａｔｉｖｅｓ；ａｎｄ（ｉｉ）ｔｈｅｓｅ　ｐａｔｔｅｒｎｓ　ａｒｅ　ｐｒｅｓｅｒｖｅｄ　ｗｉｔｈ　ｔｈｅ　ｃｈａｎｇｅｄ　ｎｕｍｂｅｒ　ｏｆ　ｆｕｓｅｄ　ｂｅｎｚｅｎｅ　ｒｉｎｇｓ，ｗｉｔｈ　ｏｎｅ　ｂｅｎｚｅｎｅ，　ｔｏ　ｔｗｏ，ａｎｄ　ｔｗｏ　ｔｈｒｅｅ　ｆｕｓｅｄ　ｂｅｎｚｅｎｅ　ｔｉｎｇｓ．Ｔｈｅｓｅ　ｒｅｓｕｌｔｓ　ｓｕｇｇｅｓｔ　ｔｈａｔ　ｎｏ　ｍａｔｔｅｒ　ｈｏｗ　ｍａｎｙ　ｂｅｎｚｅｎｅ　ｔｉｎｇｓ　ａｒｅ　ｍｅｒｇｅｄ　ｗｉｔｈ　ｔｈｅ　ｆｕｌｖｅｎｅ　ｒｉｎｇ，ｔｈｅ　ｃｏｒｒｅｌａｔｉｏｎ　ｐａｔｔｅｍ　ｏｆ　ａｒｏｍａｔｉｃｉｔｙ　ｄｅｓｃｒｉｐｔｏｒｓ　ｉｓ　ａｌｗａｙｓ　ｔｈｅ　ｓａｍｅ．Ｔｈｅ　ｓａｍｅ　ｔｅｎｄｅｎｃｙ　ｉｓ　ｔｒｕｅ　ｆｏｒ　ＩＴ　ｑｕａｎｔｉｔｉｅｓ，ａｓ　ｓｈｏｗｎ　ｉｎ　Ｔａｂｌｅ　４　ｆｏｒ　ｂｅｎ　５Ｍ．ｄｉｂｅｎ　５Ｍ．ａｎｄ　ｔｒｉｂｅｎ　５Ｍ　ｓｙｓｔｅｍｓ　ａｓ　ａｎ　ｅｘａｍｐｌｅ．ｆｒｏｍ　ｗｈｉｃｈ　ｔｈｅ　ｅｘａｃｔｌｙ　ｓａｍｅ　ｃｏｎｃｌｕｓｉｏｎｓ　ｃａｎ　ｂｅ　ｄｒａｗｎ．　Ｆｏｒ　ｃｒｏｓｓ－ｃｏｒｒｅｌａｔｉｏｎｓ　ｂｅｔｗｅｅｎ　ａｒｏｍａｔｉｃｉｔｙ　ｉｎｄｉｃｅｓ　ａｎｄ　ＩＴ　ｑｕａｎｔｉｔｉｅｓ，ｏｕｒ　ｐｒｅｖｉｏｕｓ　ｒｅｓｕｌｔｓ　ｈａｖｅ　ｓｈｏｗｎ　ｔｈａｔ　ｔｈｅｉｒ　ｐａｔｔｅｒｎｓ　ａｒｅ　ｄｅｐｅｎｄｅｎｔ　ｏｆ　ｔｈｅ　ｆｕｌｖｅｎｅ　ｒｉｎｇ　ｓｉｚｅ．Ａｓ　ｓｈｏｗｎ　ｉｎ　Ｔａｂｌｅ　５．　ｔｈｉｓ　ｓａｍｅ　ｃｏｎｃｌｕｓｉｏｎ　ｉｓ　ｓｔｉｌｌ　ｖａｌｉｄ　ｆｏｒ　ｂｅｎｚｅｎｅ　ｆｕｓｅｄ　ｆｕｌｖｅｎｅ　ｄｅｒｉｖａｔｉｖｅｓ．Ｆｏｒ　ｅｘａｍｐｌｅ，ｗｅ　ｆｉｎｄ　ｔｈａｔ　ｂｅｎ　５Ｍ　ａｎｄ　ｂｅｎ　６Ｍ—　ｈａｖｅ　ｏｐｐｏｓｉｔｅ　ｃｈａｎｇｉｎｇ　ｐａｔｔｅｒｎｓ　ｂｅｔｗｅｅｎ　ａｒｏｍａｔｉｃｉｔｙ　ｉｎｄｉｃｅｓ　ｎｄ　ＩＴ　ｑｕａｎｔｉａｔｉｅｓ．Ｉｎ　ａｄｄｉｔｉｏｎ，ｔｈｉｓ　ｏｐｐｏｓｉｔｅ　ｃｈａｎｇｅ　ｐａｔｔｅｒｎ　ｐｅｒｓｉｓｔｓ　ｗｉｔｈ　ｔｈｅ　ｉｎｃｒｅａｓｅｄ　ｎｕｍｂｅｒ　ｏｆ　ｕｓｅｄ　ｆｂｅｎｚｅｎｅ　ｒｉｎｇｓ．Ｔｈａｔ　ｉｓ，ｆｒｏｍ　ｂｅｎ　５Ｍ，ｔｏ　ｄｉｂｅｎ　５Ｍ，ａｎｄ　ｔｏ　ｔｒｉｂｅｎ　５Ｍ，ｔｈｅ　ｓｉｇｎ　ｏｆ　ｔｈｅ　ｃｏｒｒｅｌａｔｉｏｎ　ｃｏｅｆｉｃｉｆｅｎｔ　ｄｏｅｓ　ｎｏｔ　ｃｈａｎｇｅ．Ａｇａｉｎ，ｔｈｅｓｅ　ｒｅｓｕｌｔｓ　ｉｎｄｉｃａｔｅ　ｔｈａｔ　ｔｈｅ　ｃｏｒｒｅｌａｔｉｏｎ　ｐａｔｔｅｒｎ　ｏｆ　ｂｅｎｚｅｎｅ　ｆｕｓｅｄ　ｆｕｌｖｅｎｅ　ｄｅｒｉｖａｔｉｖｅｓ　ｉｓ　ｄｉｃｔａｔｅｄ　ｔｈｅ　ｓｕｂｓｔｉｔｕｔｅｄ　ｆｕｌｖｅｎｅｓ．　ａｎｄ（ｉｉｉ）ｄｉｆｅｒｅｎｔ　ａｒｏｍａｔｉｃｉｔｙ　ｄｅｓｃｒｉｐｔｏｒｓ　ｏｆｔｅｎ　ｙｉｅｌｄ　ｄｉｆｆｅｒｅｎｔ　ｃｏｎｃｌｕｓｉｏｎｓ　ｏｎ　ｗｈｅｔｈｅｒ　ａ　ｓｐｅｃｉｅｓ　ｉｓ　ａｒｏｍａｔｉｃ　ｏｒ　ａｎｔｉａｒｏｍａｔｉｃ．　Ｔａｂｌｅ　２　ｅｘｈｉｂｉｔｓ　ｔｈｅ　ａｖｅｒａｇｅ　ｖａｌｕｅ　ｏｆ　ｒｉｎｇ　ｃａｒｂｏｎ　ａｔｏｍｓ　ｏｆ　ｅｉｇｈｔ　ｉｎｆｏｒｍａｔｉｏｎ—ｔｈｅｏｒｅｔｉｃ　ｑｕａｎｔｉｔｉｅｓ，Ｓｈａｎｎｏｎ　ｅｎｔｒｏｐｙ　Ｓｓ，　Ｆｉｓｈｅｒ　ｉｎｆｏｒｍａｔｉｏｎ／Ｆ，Ｇｈｏｓｔ—Ｂｅｒｋｏｗｉｔｚ—Ｐａｒｒ　ｅｎｔｒｏｐＹ＆ＢＰ，　ｉｎｆｏｒｍａｔｉｏｎ　ｇａｉｎ／ｃ，Ｏｎｉｃｅｓｃｕ　ｉｎｆｏｒｍａｔｉｏｎ　ｅｎｅｒｇｙ　ｏｆ　ｏｒｄｅｒｓ　２　ａｎｄ　３，Ｅ２　ａｎｄ　Ｅ３，ａｎｄ　ｒｅｌａｔｉｖｅ　Ｒｄｎｙｉ　ｅｎｔｒｏｐｙ　ｏｆ　ｏｒｄｅｒｓ　２　ａｎｄ　３，　ｒ２　ａｎｄ　Ｒｒ３．ｆｏｒ　Ｒｉｎｇ　Ｂ　ｏｆ　ｔｈｅ　ｄｉｂｅｎ　５Ｍ　ｓｙｓｔｅｍ　ａｓ　ｉｌｌｕｓｔｒａｔｉｖｅ　ｅｘａｍｐｌｅｓ．Ｃｏｍｐａｒｉｎｇ　ｗｉｔｈ　ｔｈｅ　ｄａｔａ　ｆｒｏｍ　ｏｕｒ　ｐｒｅｖｉｏｕｓ　ｓｔｕｄｙ　ｆｏｒ　ｓｕｂｓｔｉｕｔｅｄ　ｆｕｌｔｖｅｎｅ　ｄｅｒｉｖａｔｉｖｅｓ，ｗｅ　ｆｏｕｎｄ　ｔｈａｔ　ｅａｃｈ　ｏｆ　ｔｈｅｓｅ　ＩＴ　ｑｕａｎｔｉｔｉｅｓ　ｉｓ　ｓｌｉｇｈｔｌｙ　ｆｌｕｃｔｕａｔｉｎｇ　ａｒｏｕｎｄ　ｉｔｓ　ｃｈａｒａｃｔｅｒｉｓｔｉｃ　（ａｖｅｒａｇｅ）ｖａｌｕｅ，ｂｕｔ，ａｓ　ｗｉｌｌ　ｂｅ　ｓｈｏｗｎ　ｂｅｌｏｗ，ｉｔ　ｉｓ　ｔｈｅ　ｐａｔｔｅｒｎ　ｏｆ　ｔｈｅｉｒ　ｆｌｕｃｔｕａｔｉｏｎ　ｔｈａｔ　ｗｉｌ１　ｒｅｆｌｅｃｔ　ｔｈｅ　ｎａｔｕｒｅ　ｏｆ　ａｒｏｍａｔｉｃｉｔｙ　ａｎｄ　ａｎｔｉａｒｏｍａｔｉｃｉｔｙ　ｏｆ　ｔｈｅｓｅ　ｓｙｓｔｅｍｓ．　Ｆｏｌｌｏｗｉｎｇ　ｏｕｒ　ｐｒｅｖｉｏｕｓ　ｗｏｒｋ　３８，６９，７１　７６，Ｔａｂｌｅｓ　３－５　ｄｉｓｐｌａｙ　ｔｈｅ　ｃｏｒｒｅｌａｔｉｏｎ　ｃｏｅ伍ｃｉｅｎｔｓ　ａｍｏｎｇ　ｔｈｅ　ａｒｏｍａｔｉｃｉｔｙ　ａｎｄ　ＩＴ　ｑｕａｎｔｉｔｉｅｓ　ｓｔｕｄｉｅｄ　ｉｎ　ｔｈｉｓ　ｗｏｒｋ．Ｔａｂｌｅ　３　ｓｈｏｗｓ　ｔｈｅ　ｉｎｔｅｒ－ｃｏｒｒｅｌａｔｉｏｎ　ｂｅｔｗｅｅｎ　ｅａｃｈ　ｐａｉｒ　ｏｆ　ａｒｏｍａｔｉｃｉｔｙ　ｄｅｓｃｒｉｐｔｏｒｓ　ｏｒ　ｆＲｉｎｇ　Ｂ　ｏｆ　ｎｉｎｅ　ｂｅｎｚｅｎｅ　ｆｕｓｅｄ　ｆｕｌｖｅｎｅ　ｓｙｓｔｅｍｓ，ｆｒｏｍ　ｗｈｉｃｈ　ｗｅ　ｃａｎ　ｓｅｅ　ｔｈａｔ　ｔｈｅ　ｃｏｒｒｅｌａｔｉｏｎ　ｐａｔｔｅｍ　ｆｉ．ｅ．．ｐｏｓｉｔｉｖｅ　ａｎｄ　ｎｅｇａｔｉｖｅ　ｓｉｇｎｓ　ｏｆ　ｔｈｅ　ｃｏｒｒｅｌａｔｉｏｎ　ｃｏｅｆｉｃｉｆｅｎｔ）ｉｓ　ａｌｗａｙｓ　ｔｈｅ　Ｔｏ　ｖｉｓｕａｌｉｚｅ　ｔｈｅ　ａｂｏｖｅ　ｒｅｓｕｌｔｓ．ｕｓｉｎｇ　ｔｈｅ　ｃｏｒｒｅｌａｔｉｏｎ　ｂｅｔｗｅｅｎ　ｓａｍｅ．Ｔｈａｔ　ｉｓ，ＦＬＵ　ｖｓ　ＮＩＣＳ（１）ｉｓ　ｐｏｓｉｔｉｖｅｌｙ　ｃｏｒｒｅｌａｔｅｄ　ｆｗｉｔｈ　ａ　ｐｏｓｉｔｉｖｅ　ｃｏｒｒｅｌａｔｉｏｎ　ｃｏｅｆｉｃｉｆｅｎｔ）．ｗｈｅｒｅａｓ　ＦＬＵ　ｖｓ　ＨＯＭＡ　ａｎｄ　Ｈ０ＭＡ　ｖｓ　ＮＩＣＳ（１１　ｒｅ　ａｎｅｇａｔｉｖｅｌｙ　ｃｏｒｒｅｌａｔｅｄ．Ｔｗｏ　ｐｏｉｎｔｓ　ａｒｅ　ｔｏ　ＮＩＣＳ（１）ａｎｄ　ＳＧＢＰ　ａｓ　ａｎ　ｅｘａｍｐｌｅ，Ｆｉｇ．１　ｅｘｈｉｂｉｔｓ　ｔｈｅ　ｓｔｒｏｎｇ　ｌｉｎｅａｒ　ｒｅｌａｔｉｏｎｓｈｉｐｓ　ｏｎ　ｔｈｅ　ｆｕｌｖｅｎｅ　ｒｉｎｇ　ｏｆ　０—３　ｂｅｎｚｅｎｅ　ｆｕｓｅｄ　ｓｙｓｔｅｍｓ．Ｔｌ１ｅ　ａｖｅｒａｇｅ　ｖａｌｕｅ　ｏｆ　ａｌｌ　ｃａｒｂｏｎ　ａｔｏｍｓ　ｏｎ　ｔｈｅ　ｆｕｌｖｅｎｅ　Ａｃｔａ　Ｐｈｙｓｉｃｏ—Ｃｈｉｍｉｃａ　Ｓｉｎｉｃａ　Ｖ＿ｏ１．３４　Ｔａｂｌｅ　２　Ｎｕｍｅｒｉｃａｌ　ｖａｌｕｅｓ　ｏｆ　ｉｎｆｏｒｍａｔｉｏｎ．ｔｈｅｏｒｅｔｉｃ　ｑｕａｎｔｉｔｉｅｓ　ｆｏｒ　Ｒｉｎｇ　Ｂ　ｏｆ　ｄｉｂｅｎ５Ｍ　ａｓ　ｉｌｌｕｓｔｒａｔｉｖｅ　ｅｘａｍｐｌｅｓ·　＿Ｔａｂｌｅ　３　Ｔｈｅ　ｃｏｒｒｅｌａｔｉｏｎ　ｃｏｅｆｉｃｉｆｅｎｔ　Ｒ　ｂｅｔｗｅｅｎ　ａｒｏｍａｔｉｃｉｔｙ　ｉｎｄｅｘｅｓ　ＦＬＵ，ＨＯＭＡ，ａｎｄ　ＮＩＣＳ（１）ｆｏｒ　ｉｎｇ　ＲＢ　ｏｆ　ａ　ｆｅｗ　ｂｅｎｚｅｎｅ　ｆｕｓｅｄ　ｆｕｌｖｅｎｅ　ｓｙｓｔｅｍｓ·　ｒｉｎｇ　ｗａｓ　ｅｍｐｌｏｙｅｄ　ｉｎ　ｔｈｅ　Ｆｉｇｕｒｅ　３８，７７．Ａ　ｄｏｔｔｅｄ　ｌｉｎｅ　ｉｓ　ｕｓｅｄ　ｉｎ　ｃｌｅａｒ　ｔｈａｔ　ｔｈｅ　ｃｏｒｒｅｌａｔｉｏｎｓ　ａｒｅ　ｓｔｒｏｎｇｅｒ　ｏｎ　Ｒｉｎｇ　Ａ　ｔｈａｎ　Ｒｉｎｇ　Ｂ，　ｅｖｅｎ　ｔｈｏｕｇｈ　ｔｈｅ　ｑｕａｌｉｔａｔｉｖｅ　ｔｒｅｎｄ　ｉＳ　ｅｖｉｄｅｎｔ．Ｗｉｔｈ　ｔｈｅ　ｄｉｆｆｅｒｅｎｃｅ　ｉｎ　ｔｈａｔ．ｔｈｅ　ｃｏｒｒｅｌａｔｉｏｎｓ　ｉｎ　Ｆｉｇ．２　ａｒｅ　ｆｏｒ　ｏｎｅ　ｏｆ　ｔｈｅ　ｆｕｓｅｄ　ｔｈｅ　Ｆｉｇｕｒｅｔｏ　ｓｅｐａｒａｔｅ　ａｒｏｍａｔｉｃ（Ｎ￣ｃｓ０）＜０）ａｎｄ　ａｎｔｉａｒｏｍａｔｉｃ　（Ｎ￣ＣＳ（１）＞０）ｒｅｇｉｏｎｓ　３０，３１．Ａｓ　ｃａｎ　ｂｅ　ｓｅｅｎ　ｆｒｏｍ　ｔｈｅ　ｆｉｇｕｒｅ，　ｔｈｅｒｅ　ａｒｅ　ｔｗｏ　ｄｉｆｆｅｒｅｎｔ　ｃａｔｅｇｏｒｉｅｓ　ｏｆ　１ｉｎｅａｒ　ｃｏｒｒｅｌａｔｉｏｎｓ，ｏｎｅ　ｂｅｎｚｅｎｅ　ｒｉｎｇｓ，ｎｏｔ　ｆｏｒ　ｔｈｅ　ｆｕｌｖｅｎｅ　ｒｉｎｇ．Ｃｏｍｐａｒｅｄ　ｗｉｔｈ　ｔｈｅｉｒ　ｃｏｕｎｔｅｒｐａｒｔ　ｉｎ　Ｆｉｇ．１，ｉｔ　ｂｅｃｏｍｅｓ　ａｐｐａｒｅｎｔ　ｔｈａｔ　ｔｈｅ　ｐａｔｔｅｒｎ　ｏｆ　ｃｈａｎｇｅｓ　ｗｉｔｈ　ｒｅｓｐｅｃｔ　ｔｏ　ｔｈｅ　ｃｈａｎｇｅ　ｏｆ　ｓｕｂｓｔｉｔｕｔｉｎｇ　ｇｒｏｕｐｓ　ｗｉｔｈ　ａ　ｐｏｓｉｔｉｖｅ　ｓｌｏｐｅ　ａｎｄ　ｔｈｅ　ｏｔｈｅｒ　ｗｉｔｈ　ａ　ｎｅｇａｔｉｖｅ　ｓｌｏｐｅ，ａｎｄ　ｔｈｅｓｅ　ｐａｔｔｅｒｎｓ　ｄｏ　ｎｏｔ　ｃｈａｎｇｅ　ｗｉｔｈ　ｒｅｓｐｅｃｔ　ｔｏ　ｔｈｅ　ｎｕｍｂｅｒ　ｏｆ　ｔｈｅ　ｆｕｓｅｄ　ｂｅｎｚｅｎｅｓ　７８，７９．Ｎｏｔｉｃｅ　ｔｈａｔ，ａｓ　ｉｎｄｉｃａｔｅｄ　ｅａｒｌｉｅｒ，ｔｈｅ　ｔｏｔａｌ　ｎｕｍｂｅｒ　ｏｆ　ｅｌｅｃｔｒｏｎｓ　ｉｎ　ｔｈｅｓｅ　ｔｗｏ　ｃａｔｅｇｏｒｉｅｓ　ｏｆ　ｃｏｒｒｅｌａｔｉｏｎｓ　ｒｅ　ｉｎ　ｃｏｎｓｅｎｓｕｓ　ｗｉａｔｈ　Ｈｆｉｃｋｅｌ’Ｓ　４ｎ＋２　ｒｕｌｅ　ｏｆａｒｏｍａｔｉｃｉｔｙ　ａｎｄ　４ｎ　ｒｕｌｅ　ｏｆａｎｔｉａｒｏｍａｔｉｃｉｔｙ　４１．　Ｆｉｇ．２　ｄｉｓｐｌａｙｓ　ｔｈｅ　ｓａｍｅ　ｃｏｒｒｅｌａｔｉｏｎ　ｐａｔｔｅｍｓ　ｂｅｔｗｅｅｎ　ＳＧＢＰ　ｂｅｔｗｅｅｎ　ＮＩＣＳ（１）ａｎｄ　ＳｏＢＰ　ｉｓ　ｔｈｅ　ｓａｍｅ，ｃｏｎｆｉｒｍｉｎｇ　ｗｈａｔ　ｗｅ　ｃｏｎｃｌｕｄｅｄ　ｅａｒｌｉｅｒ　ｔｈａｔ　ｉｔ　ｉＳ　ｔｈｅ　ｆｕｌｖｅｎｅ　ｒｉｎｇ　ｔｈａｔ　ｄｉｃｔａｔｅｓ　ｔｈｅ　ｃｈａｎｇｉｎｇ　ｐａｔｔｅｒｎ　ｏｆ　ｔｈｅ　ｂｅｎｚｏｆｕｌｖｅｎｅ　ｓｙｓｔｅｍｓ　３８．Ｗｅ　ｋｎｏｗ　ｔｈａｔ　ｂｅｎｚｅｎｅ　ｉｔｓｅｌｆ　ｉＳ　ａｒｏｍａｔｉｃ，ｂｕｔ　ｗｈｅｎ　ｆｕｓｅｄ　ｗｉｔｈ　ｆｕｌｖｅｎｅ，　ｓｏｍｅｔｉｍｅ　ｉｔ　ｂｅｃｏｍｅｓ　ａｎｔｉａｒｏｍａｔｉｃ，ａｃｃｏｒｄｉｎｇ　ｔｏ　ｔｈｅ　ｃｒｉｔｅｒｉｏｎ　ｏｆ　ｎｄ　ａＮＩＣＳ（１）ｏｆ　ｒｉｎｇ　Ａ　ｆｏｒ　ｔｈｒｅｅ　ｂｅｎｚｏｆｕｌｖｅｎｅｓ　ａｓ　Ｆｉｇ．１．Ｉｔ　ｉｓ　ＮＩＣＳ（１）２９　３１，４５．　ＮＯ．６　ｄｏｉ：１　０．３　８６６／ＰＫＵ．ＷＨＸＢ２０　１　７　１　０２３　１　６４５　Ｔａｂｌｅ　５　Ｔｈｅ　ｃｏｒｒｅｌａｔｉｏｎ　ｃｏｅｆｉｃｉｆｅｎｔｓ　ｂｅｔｗｅｅｎ　ａｒｏｍａｔｉｃ　ｉｎｄｅｘｅｓ　ａｎｄ　ＩＴ　ｑｕａｎｔｉｉｔｅｓ　ｆｏｒ　Ｒｉｎｇ　Ｂ　ｏｆ　ａ　ｆｅｗ　ｂｅｎｚｅｎｅ　ｆｕｓｅｄ　ｆｕｌｖｅｎｅ　ｄｅｒｉｖａｔｉｖｅｓ．　６４６　Ａｃｔａ　Ｐｈｙｓｉｃｏ—Ｃｈｉｍｉｃａ　Ｓｉｎｉｃａ　Ｖｌ０１．３４　３ｇ　６　§　ｍｃｓ（１１＝０　。　。　。｛　３　Ｍｃ５ｒ工Ｊ　０　ａ㈣ａｆ　ｍａ　Ｎｉｃｓｏ］。　Ｆｉｇ．Ｉ　Ｓｔｒｏｎｇ　ｌｉｎｅａｒ　ｃｏｒｒｅｌａｔｉｏｎｓ　ｂｅｔｗｅｅｎ　ＳＧａｒ　ａｎｄ　ＮＩＣＳ（ｔ）ｆｏｒ（ａ）ｓｕｂｓｔｉｔｕｔｅｄ　ｆｏｌｖｅｎｅ　ｄｅｒｉｖａｔｉｖｅｓ　ａｎｄ　ｆｏｒ　Ｒｉｎｇ　Ｂ　ｏｆ（ｂ）ｓｉｎｇｌｙ　ｆｕｓｅｄ　ｂｅｎｚｏ　ｆｕｌｖｅｎｅｓ；（ｃ）ｄｏｕｂｌｙ　ｆｕｓｅｄ　ｂｅｎｚｏｆｕｈ＇ｅｎｅｓ；ａｎｄ（ｄ）ｔｒｉｐｌｙ　ｆｕｓｅｄ　ｂｅｎｚｏｆｕｌｖｅｎｅｓ·　Ｔｈｅ　ｄｏｔｔｅｄｌｉｎｅｉｓＮＩＣＳｆＩ）＝０　０　缱｛＆ｒ　ｆ　｝“￣ｉ！：ＡＳＥ＝０　３８　７　．２　１　Ｎ托ｓｌ蛆　１５　．６　Ｉ　３Ｎ｜ｃｓｔｌ｝　Ｆｉｇ．２　Ｓｔｒｏｎｇ　ｌｉｎｅａｒ　ｃｏｒｒｅｌａｔｉｏｎｓ　ｂｅｔｗｅｅｎ＆　ＢＰ　ａｎｄ　ＮＩＣＳ（Ｉ）ｏｆ　Ｒｉｎｇ　Ａ　ｆｏｒ　ｆｏｒ（ａ）ｓｉｎｇｌｙ　ｆｕｓｅｄ　ｂｅｎｚｏ　ｆｕｌｖｅｎｅｓ；（ｂ）ｄｏｕｂｌｙ　ｆｕｓｅｄ　ｂｅｎｚｏｆｕｌｖｅｎｅｓ；　ａｎｄ（ｅ）ｔｒｉｐｌｙ　ｆｕｓｅｄ　ｂｅｎｚｏｆｕｌｖｅｎｅｓ．　Ｔｈｅ　ｄｏｔｔｅｄ　ｌｉｎｅ　ｉｓ　ＮＩＣＳ（１１＝０　Ｆｉｇ．３　ＳＷｏｎｇ　ｌｉｎｅａｒ　ｃｏｒｒｅｌａｔｉｏｎｓ　ｂｅｔｗｅｅｎ＆：ＢＰ　ａｎｄ　ＡＳＥ　ｆｏｒ（ａ）　Ｔｏ　ｖｅｒｉｒｙ　ｏｕｒ　ａｂｏｖｅ　ｆｉｎｄｉｎｇｓ，Ｆｉｇ．３　ｓｈｏｗｓ　ｔｈｅ　ＡＳＥ　ｓｕｂｓｔｉｔｕｔｅｄ　ｆａｌｖｅｎｅ　ｄｅｒｉｖａｔｉｖｅｓ　ｏｆ　ｄｉｆｆｅｒｅｎｔ　ｒｉｎｇ　ｓｉｚｅｓ；ａｎｄ（ｂ）ｓｉｎｇｌｙ　ｆｕｓｅｄ　ｂｅｎｚｏｆｕＩｖｅｎｅｓ．　Ｔｈｅ　ｄｏｔｔｅｄｌＩｎｃｉｓＡＳＥ＝０　ｃｏｒｒｅｌａｔｉｏｎ　ｒｅｓｕｌｔｓ　ｗｉｔｈ　ｏｎｅ　ｏｆ　ＩＴ　ｑｕａｎｔｉｔｉｅｓ，ＳｃＢＰ，ｆｏｒ　ｓｐｅｃｉｅｓ　（ａ）ｗｉｔｈｏｕｔ　ａｎｄ（ｂ）ｗｉｔｈ　ｏｎｅ　ｂｅｎｚｅｎｅ　ｆｕｓｉｏｎ．Ａｇａｉｎ，ａ　ｄｏｔｔｅｄ　ｌｉｎｅ　ｉＳ　ｕｓｅｄ　ｔｏ　ｓｅｐａｒａｔｅ　ａｒｏｍａｔｉｃ　ａｎｄ　ａｎｔｉａｒｏｍａｔｉｃ　ｒｅｇｉｏｎｓ．　Ｃｏｍｐａｒｉｎｇ　ｔｈｅ　ｐｌｏｔｓ　ｆｏｒ　ｔｈｅ　ｔｗｏ　ｓｙｓｔｅｍｓ．ｗｉｔｈ　ａｎｄ　ｗｉｔｈｏｕｔ　ｐｏｓｉｔｉｖｅｌｙ　ｓｌｏｐｅｄ　ｌｉｎｅｓ　ａｎｄ　ｔｗｏ　ｎｅｇａｔｉｖｅｌｙ　ｓｌｏｐｅｄ　ｌｉｎｅｓ，　ｒｅｓｐｅｃｔｉｖｅｌｙ．Ｐｕｔ　ｔｏｇｅｔｈｅｒ，ｏｕｒ　ｒｅｓｕｌｔｓ　ｄｅｍｏｎｓｔｒａｔｅ　ｔｈａｔ　ｂｅｎｚｅｎｅ　ｆｕｓｉｏｎ　ｄｏｅｓ　ｎｏｔ　ｃｈａｎｇｅ　ｔｈｅ　ｃｏｒｒｅｌａｔｉｏｎ　ｐａｔｔｅｒｎｓ　ａｍｏｎｇ　ｂｅｎｚｅｎｅ　ｆｕｓｉｏｎ，ｉｔ　ｃａｎ　ｂｅ　ｒｅａｄｉｌｙ　ｓｅｅｎ　ｔｈａｔ　ｔｈｅ　ｃｈａｎｇｉｎｇ　ｐａｔｔｅｒｎｓ　ｏｆ　ＡＳＥ　ｗｉｔｈ　ｒｅｓｐｅｃｔ　ｔｏ　ＳＧＢＰ　ａｒｅ　ｔｈｅ　ｓａｍｅ，ｅａｃｈ　ｗｉｔｈ　ｔｗｏ　ＮＯ．６　ｄｏｉ：ｌ０．３８６６／ＰＫＵ．ＷＨＸＢ２０１７１０２３１　６４７　ａｒｏｍａｔｉｃｉｔｙ　ｉｎｄｉｃｅｓ，ａｍｏｎｇ　ＩＴ　ｑｕａｎｔｉｔｉｅｓ，ａｓ　ｗｅｌｌ　ａｓ　ｃｒｏｓｓ　ｃｏｒｒｅｌａｔｉｏｎｓ　ｂｅｔｗｅｅｎ　ｔｈｅｓｅ　ｔｗｏ　ｃａｔｅｇｏｒｉｅｓ　ｏｆ　ｑｕａｎｔｉｔｉｅｓ．Ｔｈｅ　ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ　ｏｆ也ｅｓｅ　ｃｏｒｒｅｌａｔｉｏｎｓ　ｆｏｒ　ｔｈｅ　ｂｅｎｚｅｎｅ　ｆｕｓｅｄ　ｆｕｌｖｅｎｅ　ｄｅｒｉｖａｔｉｖｅｓ　ｉｓ　ｄｉｃｔａｔｅｄ　ｂｙ　ｔｈｅ　ｆｕｌｖｅｎｅ　ｍｏｉｅｔｙ．　；Ｏｓｕｋａ，Ａ．；Ｋｉｍ，Ｄ．Ｎａｔ．Ｃｈｅｍ．２００９，，，１１３．　ｆ１Ｏ）　Ｙｏｏｎ，Ｚ．Ｓ．ｄｏｉ：１０．１０３８／ｎｃｈｅｍ．１７２　ｃｋｓｏｎ，Ｃ．Ｋ．；Ｚａｋｈａｒｏｖ，Ｌ．Ｎ．：Ｈａｌｅｙ，Ｍ．Ｍ．　ＡｍＣｈｅｍ．（１１）　Ｆｒｅｄｅｒｉ　．Ｓｏｃ．２０１６，１３８，１６８２７．ｄｏｉ：１０．１０２１￣ａｃｓ．６ｂｌ１３９７　５　Ｃｏｎｃｌｕｓｉｏｎｓ　Ａｓ　ａ　ｃｏｎｔｉｎｕａｔｉｏｎ　ｏｆ　ｏｕｒ　ｅｎｄｅａｖｏｒ　ｔｏ　ａｐｐｒｅｃｉａｔｅ　ｍｏｌｅｃｕｌｒ　ａｄｙｒｅｖ，Ａ．Ｉ．；Ｐｏｐｏｖ，Ｉ．Ａ．；Ｓｔａｒｉｋｏｖａ，Ａ．Ａ．；Ｓｔｅｇｌｅｎｋｏ，Ｄ．Ｖ　（１２）　ＢｏｌＣｈｅｍ．一　２０１７．ｄｏｉ：１０．１００２／ｃｈｅｍ．２０１７０２０３５　ａｒｏｍａｔｉｃｉｔｙ　ａｎｄ　ａｎｔｉａｒｏｍａｔｉｃｉｔｙ　ｕｓｉｎｇ　ｉｎｆｏｒｍａｔｉｏｎ．ｔｈｅｏｒｅｔｉｃ　ｑｕａｎｔｉｔｉｅｓ　ｆｒｏｍ　ｄｅｎｓｉｔｙ　ｆｕｎｃｔｉｏｎａｌ　ｒｅａｃｔｉｖｉｔｙ　ｔｈｅｏｒｙ，ｉｎ　ｔｈｉｓ　ｗｏｒｋ，ｗｅ　ｈａｖｅ　ｓｙｓｔｅｍａｔｉｃａｌｌｙ　ｅｘａｍｉｎｅｄ　ｓｕｂｓｔｉｔｕｔｅｄ　ｆｕｌｖｅｎｅ　ｄｅｒｉｖａｔｉｖｅｓ　ｆｕｓｅｄ　ｗｉｔｈ　ｏｎｅ，ｔｗｏ，ａｎｄ　ｔｈｒｅｅ　ｂｅｎｚｅｎｅ　ｒｉｎｇｓ．Ｗｅ　；Ａｎｇｌａｄａ，Ｊ．Ｍ．；Ｔｏｒｒｅｎｔ—Ｓｕｃａｒｒａｔ，Ｍ．　Ｏｒｇ．Ｃｈｅｍ．　（１３）　Ｍａｒｃｏｓ，Ｅ．２０１４，７９，５０３６．ｄｏｉ：１０．１０２１￣ｏ５００５６９ｐ　ｅｆｌ，Ｍ．；Ｌａｔｏｓ－Ｇｒａ２ｙｆｉｓｋｉ，Ｌ．；Ｓｐｒｕｔｔａ，Ｎ．；Ｃｈｗａｌｉｓｚ，Ｐ；　（１４）　Ｓｔ￣ｐｉｈａｖｅ　ｃｏｎｓｉｄｅｒｅｄ　ｔｈｅ　ｃｏｒｒｅｌａｔｉｏｎｓ，ｆｏｒ　ｔｈｅｓｅ　ｓｙｓｔｅｍｓ，ａｍｏｎｇ　ｆｏｕｒ　ａｒｏｍａｔｉｃｉｔｙ　ｄｅｓｃｒｉｐｔｏｒｓ，ＡＳＥ，ＦＬＵ，ＨＯＭＡ，ａｎｄ　ＮＩＣＳ，ａｎｄ　ｓｉｘ　ｉｎｆｏｒｍａｔｉｏｎ—ｔｈｅｏｒｅｔｉｃ　ｑｕａｎｔｉｔｉｅｓ，Ｓｈａｎｎｏｎ　ｅｎｔｒｏｐｙ，Ｆｉｓｈｅｒ　ｉｎｆｏｒｍａｔｉｏｎ，Ｇｈｏｓｈ·Ｂｅｒｋｏｗｉｔｚ—Ｐａｒｒ　ｅｎｔｒｏｐｙ，ｉｎｆｏｒｍａｔｉｏｎ　ｇａｉｎ，　Ｏｎｉｃｅｓｃｕ　ｉｎｆｏｒｍａｔｉｏｎ　ｅｎｅｒｇｙ，ａｎｄ　ｒｅｌａｔｉｖｅ　Ｒｅｎｙｉ　ｅｎｔｒｏｐｙ，ａｓ　ｗｅｌｌ　ａｓ　ｃｒｏｓｓ　ｃｏｒｒｅｌａｔｉｏｎｓ　ｂｅｔｗｅｅｎ　ｔｈｅｓｅ　ｔｗｏ　ｃａｔｅｇｏｒｉｅｓ　ｏｆ　ｐｒｏｐｅｒｔｉｅｓ．Ｏｕｒ　ｒｅｓｕｌｔｓ　ｓｈｏｗ　ｔｈａｔ　ｗｉｔｈ　ａｎｄ　ｗｉｔｈｏｕｔ　ｆｕｓｅｄ　ｂｅｎｚｅｎｅ　ｒｉｎｇｓ，ａｂｏｖｅ　ｃｏｒｒｅｌａｔｉｏｎ　ｐａｔｔｅｍｓ　ａｒｅ　ｔｈｅ　ｓａｍｅ，　ｓｕｇｇｅｓｔｉｎｇ　ｔｈａｔ，ｅｖｅｎ　ｔｈｏｕｇｈ　ｗｅ　ｋｎｏｗ　ｔｈａｔ　ｂｅｎｚｅｎｅ　ｉｔｓｅｌｆ　ｉｓ　ｒａｏｍａｔｉｃ，ｉｔｓ　ｆｕｓｉｏｎ　ｗｉｔｈ　ｔｈｅ　ｆｕｌｖｅｎｅ　ｍｏｉｅｔｙ　ｄｏｅｓ　ｎｏｔ　ｃｈａｎｇｅ　ｔｈｅ　ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ　ｏｆ　ｔｈｅ　ａｒｏｍａｔｉｃ　ａｎｄ　ａｎｔｉａｒｏｍａｔｉｃ　ｎａｔｕｒｅ　ｏｆ　ｔｈｅ　ｆｕｌｖｅｎｅ　ｒｉｎｇ．　Ｉｎ　ｏｔｈｅｒ　ｗｏｒｄｓ，　ａｒｏｍａｔｉｃｉｔｙ　ａｎｄ　ａｎｔｉｒａｏｍａｔｉｃｉｔｙ　ｏｆ　ｂｅｎｚｅｎｅ　ｆｕｓｅｄ　ｆｕｌｖｅｎｅ　ｄｅｒｉｖａｔｉｖｅｓ　ａｒｅ　ｓｏｌｅｌｙ　ｄｅｔｅｒｍｉｎｅｄ　ｂｙ　ｔｈｅ　ｆｕｌｖｅｎｅ　ｍｏｉｅｔｙ．Ｔｈｅｓｅ　ｒｅｓｕｌｔｓ　ｓｈｏｕｌｄ　ｓｈｅｄ　ｎｅｗ　ｌｉｇｈｔ　ｏｎ　ｏｕｒ　ｕｎｄｅｒｓｔｎａｄｉｎｇ　ａｂｏｕｔ　ｔｈｅ　ｏｒｉｇｉｎ　ａｎｄ　ｎａｔｕｒｅ　ｏｆ　ａｒｏｍａｔｉｃｉｔｙ　ｎａｄ　ｎａｔｉａｒｏｍａｔｉｃｉｔｙ　ｏｆｒ　ｔｈｉｓ　ａｎｄ　ｏｔｈｅｒ　ｓｙｓｔｅｒｎｓ．　Ｒｅｆｅｒｅｎｃｅｓ　（１）Ｓｏｌｈ，Ｍ．Ｆｒｏｎｔ．Ｃｈｅｍ．２０１７，５，１．ｄｏｉ：１０．３３８９／ｆｃｈｅｍ．２０１７．０００２２　（２）Ｈｅｒｎｄｏｎ，Ｗ　Ｃ．；Ｍｉｌｌｓ，Ｎ．Ｓ．　Ｏｒｇ．Ｃｈｅｍ．２００５，７０，８４９２．　ｄｏｉ：１０．１０２１￣ｏ０５１２８９ｂ　（３）Ｓｋ０ｖ’Ａ．Ｂ．；Ｂｒｏｍａｎ，Ｓ．Ｌ．；Ｇｅｒｔｓｅｎ，Ａ．Ｓ．；Ｅｌｍ，Ｊ．；Ｊｅｖｒｉｃ，Ｍ．；　Ｃａｃｃｉ　ｎｉ，Ｍ．；Ｋａｄｚｉｏｌａ，Ａ．；Ｍｉｋｋｅｌｓｅｎ，Ｋ．Ｖ；Ｎｉｅｌｓｅｎ，Ｍ．Ｂ．　Ｃｈｅｍ．一　２０１６，２２，ｌ４５６７．ｄｏｉ：１０．１００２／ｃｈｅｍ．２０１６０ｌ１９０　（４）Ａｉｈａｒａ，Ｊ．一Ｉ．　Ａｍ．Ｃｈｅｍ．Ｓｏｃ．２００６，１２８，２８７３．　ｄｏｉ：１０　１０２１￣ａ０５６４３０ｃ　（５）Ｃｙｒａｆｉｓｋｉ，Ｍ．Ｋ．Ｃｈｅｍ．Ｒｅｖ．２００５，１０５，３７７３．　ｄｏｉ：１０．１０２ｌ／ｃｒ０３００８４５　（６）Ｆｅｉｘａｓ，Ｆ；Ｍａｔｉｔｏ，Ｅ．；Ｐｏａｔｅｒ，Ｊ．；Ｓｏｌｈ，Ｍ．Ｃｈｅｍ．Ｓｏｃ．Ｒｅｖ．２０１５，　４４．６４３４．ｄｏｉ：１０．１０３９／Ｃ５ＣＳＯ００６６Ａ　（７）Ｓｕｍｉｔａ，Ａ．；Ｑａｓｏｎｏｏ，Ｍ．；Ｂｏｂｌａｋ，Ｋ．Ｊ．；Ｏｈｗａｄａ，Ｔ．；Ｋｌｕｍｐｐ，Ｄ．　Ａ．　ｅｍ．一Ｅｕｒ　２０１７，２３，２５６６．ｄｏｉ：１０．１００２／ｃｈｅｍ．２０１６０６０３６　（８）Ｓｚｃｚｅｐａｎｉｋ，Ｄ．ｗ；Ｓｏｌａ，Ｍ．；Ａｎｄｒｚｅｊａｋ，Ｍ．；Ｐａｗｅ￣ｅｋ，Ｂ．；　Ｄｏｍｉｎｉｋｏｗｓｋａ，Ｊ．；Ｋｕｋｕｔｋａ，Ｍ．；Ｄｙｄｕｃｈ，Ｋ．；Ｋｒｙｇｏｗｓｋｉ，　Ｍ．；　Ｓｚａｔｙｌｏｗｉｃｚ，Ｈ．Ｊ＝Ｃｏｍｐｕｔ．Ｃｈｅｍ　２０１７，３８，１６４０．　ｄｏｉ：１０．１００２￣ｅｅ．２４８０５　（９）Ｌｉｕ，Ｊ．Ｚ．；Ｍａ，Ｊ．；Ｚｈａｎｇ，Ｋ．；Ｒａｖａｔ，Ｐ．；Ｍａｃｈａｔａ，Ｐ．；Ａｖｄｏｓｈｅｎｋｏ，　Ｓ．；Ｈｅｎｎｅｒｓｄｏｒ￡Ｆ；Ｋｏｍｂｅｒ，Ｈ．；Ｐｉｓｕｌａ，Ｗ．；Ｗｅｉｇａｎｄ，Ｊ．Ｊ．；ｅｔ　ａ１．　Ａｍ．Ｃｈｅｍ．Ｓｏｃ．２０１７，１３９，７５１３．ｄｏｉ：１０．１０２１￣ａｃｓ．７ｂ０１６１９　Ｓｚｔｅｒｅｎｂｅｒｇ，Ｌ．Ａｎｇｅｗ．Ｃｈｅｍ．２００７，Ｉ１９，８０１５．　ｄｏｉ：１０．１００２／ａｎｇｅ．２００７００５５５　（１５１　Ｓｚｙｓｚｋｏ，Ｂ．；Ｓｐｒｕｔｔａ，Ｎ．；Ｃｈｗａｌｉｓｚ，Ｐ；Ｓｔｑｐｉｅｆｌ，Ｍ．；　Ｌａｔｏｓ—Ｇｒａ２ｙｆｉｓｋｉ，Ｌ　Ｃｈｅｍ．－Ｅｕｒ．　２０１４，２０，１９８５．　ｄｏｉ：１０．１００２／ｃｈｅｍ．２Ｏ１３０３６７６　（１６）　Ｍｏｂｉｕｓ，Ｋ．；Ｐｌａｔｏ，Ｍ．；Ｋｌｉｈｍ，Ｇ．；Ｌａｕｒｉｃｈ，Ｃ．；Ｓａｖｉｔｓｋｙ，Ａ．；　Ｌｕｂｉｔｚ，　；Ｓｚｙｓｚｋｏ，Ｂ．；Ｓｔｅｐｉｅｎ，Ｍ．；Ｌａｔｏｓ—Ｇｒａｚｙｎｓｋｉ，Ｌ．　Ｃｈｅｍ．Ｃｈｅｍ．Ｐｈｙｓ　２０１５，ｉ　７，６６４４．ｄｏｉ：１０．１０３９／Ｃ４ＣＰ０５７４５Ｇ　（１７）　Ｍ６ｂｉｕｓ，Ｋ．；Ｓａｖｉｔｓｋｙ，Ａ．；Ｌｕｂｉｔｚ，　；Ｐｌａｔｏ，Ｍ．Ａｐｐ１．Ｍａｇｎ．Ｒｅｓｏｎ．　２０１６．４７，７５７．ｄｏｉ：１０．１００７／ｓ００７２３—０１６－０７８９—１　（１８）　Ｍａｒｃｏｓ，Ｅ．；Ａｎｇｌａｄａ，Ｊ．Ｍ．；Ｔｏｒｒｅｎｔ－Ｓｕｃａｒｒａｔ，Ｍ．　Ｐｈｙｓ．Ｃｈｅｍ．Ｃ　２０１２，１１６，２４３５８．ｄｏｉ：１０．１０２１￣ｐ３０８３６１２　ｆ１９）　Ｒｏｓｅｎｂｅｒｇ，Ｍ．；Ｄａｈｌｓｔｒａｎｄ，Ｃ．；Ｋｉｌｓｆｉ，Ｋ．；Ｏｔｔｏｓｓｏｎ，Ｈ．Ｃｈｅｍ　Ｒｅｖ．　２０１４，１１４，５３７９．ｄｏｉ：１０．１０２１／ｃｒ３００４７１ｖ　（２Ｏ）　Ｏｔｔｏｓｓｏｎ，Ｈ．Ｎａｔ．Ｃｈｅｍ．２０１２，４，９６９．ｄｏｉ：１０．１０３８／ｎｃｈｅｍ．１５１８　（２１）　Ｌｉｅｗ，Ｊ．Ｙ＿；Ｂｒｏｗｎ，Ｊ．Ｊ．；Ｍｏｏｒｅ，Ｅ．Ｇ．；Ｓｃｈｗａｌｂｅ，Ｍ．Ｃｈｅｍ．－Ｅｕｒ　２０１６，２２，１６１７８．ｄｏｉ：１０．１００２／ｃｈｅｍ．２０１６０２１８９　（２２　Ｎａｓｋａｒ，Ｓ．；Ｄａｓ，Ｍ．ＡＣＳ　Ｏｍｅｇａ　２０１７，２，１７９５．　ｄｏｉ：１０．１０２１／ａｃｓｏｍｅｇａ．７ｂ００２７８　（２３）　ｖａｎ　Ｌｅｅｕｗｅｎ，Ｔ；Ｐｏｌ，Ｊ．；Ｒｏｋｅ，Ｄ．；Ｗｅｚｅｎｂｅｒｇ，Ｓ．Ｊ．；Ｆｅｆｉｎｇａ，Ｂ．Ｌ．　Ｏｒｇ．Ｌｅｔｔ．２０１７，１９，１４０２．ｄｏｉ：１０．１０２１／ａｃｓ．ｏｒｇｌｅａ．７ｂ００３１７　（２４）　Ｆｅｉｘａｓ，Ｆ；Ｍａｔｉｔｏ，Ｅ．；Ｐｏａｔｅｒ，Ｊ．；Ｓｏｌｈ，Ｍ．ＷＩＲＥｓ：Ｃｏｍｐｕｔ．Ｍｏ１．　Ｓｃｉ．２０１３。３．１０５．ｄｏｉ：１０．１００２／ｗｃｍｓ．１１１５　（２５）　ＢＯｈｌ，Ｍ．；Ｈｉｒｓｃｈ，Ａ．Ｃｈｅｍ．尺　２００１，１０１，１１５３．　ｄｏｉ：１０．１０２１／ｃｒ９９０３３２ｑ　（２６）　Ｅｌ　Ｂａｋｏｕｒｉ，Ｏ．；Ｄｕｒａｎ，Ｍ．；Ｐｏａｔｅｒ，Ｊ．；Ｆｅｉｘａｓ，Ｆ；Ｓｏｌａ，Ｍ．　．　Ｃｈｅｍ　Ｃｈｅｍ．Ｐｈｙｓ．２０１６，１８，１１７００．ｄｏｉ：１０．１０３９／Ｃ５ＣＰ０７０１１Ｂ　（２７）　Ｚｈｕ，Ｃ．Ｑ．；Ｙａｎｇ，ｃ．ｘ．；Ｗａｎｇ，ＹＨ．；Ｌｉｎ，Ｇ．；Ｙａｎｇ，ＹＨ．；Ｗａｎｇ，　Ｘ．Ｙ；Ｚｈｕ，Ｊ．；Ｃｈｅｎ，Ｘ．Ｙ；Ｌｕ，Ｘ．；Ｌｉｕ，Ｇ．；ｅｔａ１．Ｓｃｉ　Ａｄｖ．２０１６，２，　ｅ１６０１Ｏ３１．ｄｏｉ：１０．１１２６／ｓｃｉａｄｖ．１６０１０３１　（２８１　Ａｙｕｂ，Ｒ．；Ｂａｋｏｕｒｉ，Ｏ．Ｅ．；Ｊｏｍｅｒ，Ｋ．；Ｓｏｌａ，Ｍ．；Ｏｔｔｏｓｓｏｎ，Ｈ．　Ｏｒｇ．　Ｃｈｅｍ．２０１７，８２，６３２７．ｄｏｉ：１０．１０２１／ａｃｓ．ｊｏｃ．７ｂ００９０６　（２９）　Ｓｕｚｕｋｉ，Ｓ．；Ｍｏｒｉｔａ，Ｙ；Ｆｕｋｕｉ，Ｋ．；Ｓａｔｏ，Ｋ．；Ｓｈｉｏｍｉ，Ｄ．；Ｔａｋｕｉ，Ｔ；　Ｎａｋａｓｕｊｉ，Ｋ．　Ａｍ．Ｃｈｅｍ．Ｓｏｃ．２００６，１２８，２５３０．　ｄｏｉ：１０．１０２１／ｊａ０５８３８７ｚ　（３０）　Ｃｈｅｎ，Ｚ．Ｆ；Ｗａｒｍｅｒｅ，Ｃ．Ｓ．；Ｃｏｒｍｉｎｂｏｅｕｆ，Ｃ．；Ｐｕｃｈｔａ，Ｒ．；ｖｏｎ　Ｒａｇｕ６　Ｓｃｈｌｅｙｅｒ，Ｐ　Ｃｈｅｍ．　Ｐｖ　２００５，１０５，３８４２．　６４８　Ａｃｔａ　Ｐｈｙｓｉｃｏ—Ｃｈｉｍｉｃａ　Ｓｉｎｉｃａ　Ｖ０１．３４　ｄｏｌ：１０．１０２ｌ／ｃｒ０３００８８＋　（３１）Ｓｃｈｌｅｙｅｒ，Ｐ　ｖ　Ｒ．；ＭａｅｒｋｅＬ　Ｃ．；Ｄｒａｎｓｆｅｌｄ，Ａ．；Ｊｉａｏ，Ｈ．；ｖａｉｌ　Ｅｉｋｅｍａ　Ｈｏｍｍｅｓ，Ｎ．Ｊ．Ｒ．　Ａｍ．Ｃｈｅｍ　Ｓｏｃ　１９９６，１１８，６３１７．　ｄｏｉ：１０．１０２１Ｏａ９６０５８２ｄ　（３２）Ｆａｌｌａｈ—Ｂａｇｈｅｒ－Ｓｈａｉｄａｅｉ，Ｈ．；Ｗａｎｎｅｒｅ，Ｃ．Ｓ．；Ｃｏｒｍｉｎｂｏｅｕｆ，Ｃ．；　Ｐｕｃｈｔａ，Ｒ．；Ｓｃｈｌｅｙｅｒ，Ｐ　ｖ　Ｒ．Ｏｒｇ　Ｌｅｔｔ．２００６，８，８６３．　ｄｏｉ：１０．１Ｏ２１／０１０５２９５４６　（３３）Ｋｒｙｇｏｗｓｋｉ，　Ｍ．；Ｓｚａｔｙｌｏｗｉｃｚ，Ｈ．；Ｓｔａｓｙｕｋ，Ｏ．Ａ．；　Ｄｏｍｉｎｉｋｏｗｓｋａ，Ｊ．；Ｐａｌｕｓｉａｋ，Ｍ．Ｃｈｅｍ．Ｒｅｖ．２０１４，１１４，６３８３．　ｄｏｉ：ｌ０．１０２１／ｃｒ４００２５２ｈ　（３４）Ｐｏａｔｅｒ，Ｊ．；Ｆｒａｄｅｒａ，Ｘ；Ｄｕｒａｎ，Ｍ．；Ｓｏｌｈ，Ｍ．Ｃｈｅｍ－Ｅｕｒ．　２００３，９，　４００．ｄｏｉ：１０　１　００２／ｃｈｅｍ．２００３９００４　１　（３５）Ｍａｔｉｔｏ，Ｅ．；Ｄｕｒａｎ，Ｍ．；Ｓｏｌｈ，Ｍ．　Ｃｈｅｍ．ＰＪ　．２００５，１２２，０１４１０９．　ｄｏｉ：１Ｏ．１０６３／１　１８２４８９５　（３６）Ｈｏｎｇ，Ｙ；Ｏｈ，Ｊ．；Ｓｕｎｇ，Ｙ　Ｍ．；Ｔａｎａｋａ，Ｙ；Ｏｓｕｋａ，Ａ．；Ｋｉｍ，Ｄ．　ＡｎｇｅｗＣｈｅｍ．１ｎｔ．Ｅｄ．２０１７，５６，２９３２．　ｄｏｉ：１０　１００２／ｎａｉｅ＿２Ｏ１６ｌ１４３ｌ　（３７）Ｆｅｉｘａｓ，Ｆ．；Ｍａｔｉｔｏ，Ｅ．；Ｐｏａｔｅｒ，Ｊ．；Ｓｏｌ／ｔ，Ｍ．　Ｃｏｍｐｕｔ．Ｃｈｅｍ　２００８，　２９．１５４３．ｄｏｉ：１０．１００２￣ｃｃ．２０９１４　（３８）Ｙｕ，Ｄ．Ｈ．；Ｒｏｎｇ，Ｃ．Ｙ；Ｌｕ，　；Ｃｈａｔｔａｒａｊ，ＲＫ．；Ｄｅ　Ｐｒｏｆｌ，Ｆ；Ｌｉｕ，Ｓ　Ｂ．Ｐｈｙｓ．Ｃｈｅｍ　Ｃｈｅｍ　Ｐｈｙｓ．２０１７，１９，、８６３５．　ｄｏｉ：１０　１０３９／Ｃ７ＣＰ０３５４４Ｆ　（３９）Ｌｉｕ，Ｓ．Ｂ．Ａｃｔａ　一Ｃｈｉｍ　Ｓｉｎ　２０１６，３２，９８．　ｄｏｉ：１０．３８６６／ＰＫＵ．ＷＨＸＢ２０１５ｌ０３０２　（４０）Ｌｉｕ，Ｓ．Ｂ．Ａｃｔａ　Ｐ厅　－Ｃｈｉｍ．Ｓｉｎ．２００９，２５，５９０．　ｄｏｉ：１　０．３８６６／ＰＫＵ．ＷＨＸＢ２００９０３３２　（４１）Ｚｈａｏ，Ｌ．Ｌ．；Ｇｒａｎｄｅ—Ａｚｔａｔｚｉ，Ｒ．；Ｆｏｒｏｕｔｎａ·Ｎｅｊａｄ，Ｃ．；Ｕｇａｌｄｅ，Ｊ．　Ｍ．；Ｆｒｅｎｋｉｎｇ，Ｇ．ＣｈｅｍｉｓｔｒｙＳｅｌｅｃｔ　２０１７，２，８６３．　ｄｏｉ：１０．１００２／ｓｌｃｔ　２０１６０２０８０　（４２）Ｃｉｅｓｉｅｌｓｋｉ，Ａ．；Ｋｒｙｇｏｗｓｋｉ，Ｔ　Ｍ；Ｃｙｒａｎｓｋｉ，Ｍ．Ｋ．；Ｂａｌａｂａｎ，Ａ．　Ｐｈｙｓ．Ｃｈｅｍ　Ｃｈｅｍ．Ｐｈｙｓ．２０１１，１３　３７３７．　ｄｏｊ：１Ｏ　ｌＯ３９／ＣＯＣＰＯ１４４６Ｊ　（４３）Ｂｅｍａｓｃｏｎｉ，Ｃ．Ｆ；Ｗｅｎｚｅｌ，Ｅ　Ｊ．　Ｏｒｇ　Ｃｈｅｍ．２０１０，７５，８４２２．　ｄｏｉ：１０．１０２１￣ｏ１０１７１９ｚ　（４４）Ｍｏ，Ｙ　Ｒ．；ＳｃｈｌｅｙｅＬ　Ｐ　ｖ．Ｒ．Ｃｈｅｍ．－Ｅｕｒ　２００６，１２，２００９．　ｄｏｉ：１Ｏ．１Ｏ０２／ｃｈｅｍ．２００５００３７６　（４５）Ｇｅｒｓｈｏｎｉ－Ｐｏｒａｎｎｅ，Ｒ．；Ｓｔａｎｇｅｒ，Ａ．Ｃｈｅｍ．Ｓｏｃ．Ｒ　２０１５，４４，６５９７　ｄｏｉ：１Ｏ．１０３９／Ｃ５ＣＳＯ０ｌ１４Ｅ　（４６）Ｄｅ　Ｏｌｉｖｅｉｒａ，Ｂ．Ｇ．Ｐｈｙｓ　Ｃｈｅｍ．Ｃｈｅｍ．Ｐｈｙｓ　２０１３，１５，３７．　ｄｏｉ：ｌ０．１０３９／Ｃ２ＣＰ４ｌ７４９Ａ　（４７）Ｚｈｕ，Ｊ．；Ａｎ，Ｋ．；Ｓｃｈｌｅｙｅｒ，Ｅ　ｖ．Ｒ．Ｏｒｇ．Ｌｅｔｔ．２０１３，１５，２４４２．　ｄｏｉ：ｌ０．１０２１／ｏ１４００９０８ｚ　（４８）　Ｃｙｒａｆｉｓｋｉ，Ｍ．Ｋ．；Ｋｒｙｇｏｗｓｋｉ，　Ｍ．；Ｋａ￣ｉｔｚｋ　Ａ　Ｒ．；Ｓｃｈｌｅｙｅｒ，Ｐ　Ｖ．Ｒ．　Ｏｒｇ．Ｃｈｅｍ　２００２，６７，１３３３．ｄｏｉ：ｌＯ．１０２１￣ｏ０１６２５５ｓ　（４９）Ｋｒｙｇｏｗｓｋｉ，　Ｍ．；Ｃｙｒａｆｌｓｋｉ，Ｍ．Ｋ．Ｃｈｅｍ．Ｒｅｖ．２００１，１０１，１３８５．　ｄｏｉ：】０．】Ｏ２】／ｃｒ９９０３２６ｕ　（５０）Ａｌｖａｒｅｚ—Ｔｈｏｎ，Ｌ．；Ｍａｍｍｉｎｏ，Ｌ．１ｎｔ．　Ｑｕａｎｔｕｍ　Ｃｈｅｍ．２０１７，１１　７，　ｅ２５３８２，ｄｏｉ：１０．１００２／ｑｕａ．２５３８２　（５　１）Ｋａｒａｄａｋｏｖ，Ｐ　Ｂ．；Ｈｅａｍｓｈａｗ，Ｐ．；Ｈｏｍｅｒ，Ｋ　Ｅ　Ｏｒｇ．Ｃｈｅｍ　２０１６，８ｌ，１１３４６．ｄｏｉ：１Ｏ．１０２１／ａｃｓ．ｊｏｃ．６ｂ０２４６０　（５２）Ｂｏｌｌ，Ｍ；Ｈｉｌｋｅｒ，　Ａ．；Ｓａｌｏｍｏｎ，Ｇ．；Ｏｍｒａｎ，Ａ．；Ｎｅｓｐｏｌｏ，Ｊ．；　Ｐｏｌｌｅｔ，Ｌ．；Ｂｌｏｃｈ，１．；Ｇｒｏｓｓ，Ｃ．Ｓｃｉｅｎｃｅ　２０１６，３５３，１２５７．　ｄｏｉ：１０．１１２６／ｓｃｉｅｎｃｅ．ａａｇ１６３５　（５３）Ｓｕｎｄｈｏｌｍ，Ｄ．；Ｆｌｉｅｇｌ，Ｈ．；Ｂｅｒｇｅｒ，Ｒ．Ｊ．Ｆ　ＷＩＲＥｓ．＂Ｃｏｍｐｕｔ．Ｍｏｌ　Ｓｃｉ．２０１６，６，６３９．ｄｏｉ：１０．１００２／ｗｃｍｓ．１２７０　（５４　Ｌｉｕ，Ｓ．Ｂ．Ｊ．Ｃｈｅｍ　ＰｈｙＳ．２００７，１２６，１９１１０７．　ｄｏｉ：１０．１０６３／１．２７４１２４４　（５５）Ｌｉｕ，Ｓ．Ｂ．；Ｒｏｎｇ，Ｃ．　；Ｗｕ，Ｚ．Ｍ．；Ｌｕ，Ｔ．Ａｃｔａ　Ｐ　—Ｃｈｉｍ　Ｓｉｎ．　２０１５，２０５７．ｄｏｉ：１０．３８６６／ＰＫＵ．ＷＨＸＢ２０１５０９１８３　（５６）Ｎｏｏｒｉｚａｄｅｈ，Ｓ．；Ｓｈａｋｅｒｚａｄｅｈ，Ｅ．Ｃｏｍｐｕｔ　Ｔｈｅｏｒ　Ｃｈｅｍ．２０１１，９６４，　１４１．ｄｏｉ：１０．１０１６￣．ｃｏｍｐｔｃ．２０１０　１２．０１２　（５７）Ｍ６１１ｅｒｓｔｅｄｔ，Ｈ．；Ｐｉｑｕｅｒａｓ，Ｍ．Ｃ；Ｃｒｅｓｐｏ，Ｒ．；Ｏｔｔｏｓｓｏｎ，Ｈ．　Ａｍ．　Ｃｈｅｍ．Ｓｏｃ．２００４，１２６，１３９３８．ｄｏｉ：１０．１０２１￣ａ０４５７２９ｃ　（５８）Ｂｅｃｋｅ，Ａ．Ｄ．　Ｃｈｅｍ．Ｐｈｙｓ．１９９３，９８，５６４８　ｄｏｉ：ｏｒｇ／１０．１０６３／１．４６４９１３　（５９）Ｌｅｅ，Ｃ．；Ｙａｎｇ，ｗ；Ｐａｒｒ，Ｒ．Ｇ　ＪＤ『　Ｒｅｖ．Ｂ１９８８，３７，７８５．　ｄｏｉ：１０．１　１０３／ＰｈｙｓＲｅｖＢ．３７．７８５　（６０）Ｆｒａｎｃｌ，Ｍ．Ｍ．；Ｐｉｅｔｒｏ，ｗ　Ｊ．；Ｈｅｈｒｅ，Ｗ．Ｊ．；Ｂｉｎｋｌｅｙ，Ｊ．Ｓ．；Ｇｏｒｄｏｎ，　Ｍ．Ｓ．；ＤｅＦｒｅｅｓ，Ｄ．Ｊ．；Ｐｏｐｌｅ，Ｊ．Ａ．　Ｃｈｅｍ．｜Ｄ　１９８２，７７，３６５４．　ｄｏｉ：１０　ｌ０６３／１　４４４２６７　（６１）Ｒａｓｓｏｌｏｖ，Ｖ　Ａ．；Ｒａｔｎｅｒ，Ｍ．Ａ．；Ｐｏｐｌｅ，Ｊ．Ａ．；Ｒｅｄｆｅｍ，Ｐ　Ｃ．；Ｃｕｒｔｉｓｓ，　Ｌ．Ａ．　Ｃｏｍｐｕｔ　Ｃｈｅｍ．２００１．２２．９７６．ｄｏｉ：１０．１００２￣ｃｃ．１０５８　（６２）Ｆｒｉｓｃｈ，Ｍ．Ｊ；Ｔｒｕｃｋｓ，Ｇ．Ｗ；Ｓｃｈｌｅｇｅｌ，Ｈ．Ｂ．；Ｓｃｕｓｅｒｉａ，Ｇ．Ｅ．；　Ｒｏｂｂ，Ｍ．Ａ．；Ｃｈｅｅｓｅｍａｎ，Ｊ．Ｒ．；Ｓｃａｌｍａｎｉ，Ｇ．；Ｂａｒｏｎｅ，Ｖ；　Ｍｅｎｎｕｃｃｉ，Ｂ．；Ｐｅｔｅｒｓｓｏｎ，Ｇ　Ａ．；ｅｔ　ａｌ　Ｇａｕｓｓｉａｎ　０９，Ｒｅｖｉｓｉｏｎ　Ｅ．Ｏ１：　Ｇａｕｓｓｉａｎ　Ｉｎｃ．：Ｗａｌｌｉｎｇｆｏｒｄ，Ｃ　ＵＳＡ，２０１３．　（６３）Ｋｕｄｉｎ，Ｋ．Ｎ．；Ｓｃｕｓｅｒｉａ，Ｇ．Ｅ．；Ｃａｎｃ６ｓ，Ｅ．　Ｃｈｅｍ．Ｐ　．２００２，　』，　８２５５．ｄｏｉ：１０．１０６３／１．１４７０１９５　（６４）Ｌｕ，　；Ｃｈｅｎ，Ｆｗ　Ｃｏｍｐｕｔ．Ｃｈｅｍ．２０１２，３３，５８０．　ｄｏｉ：１０．１００２￣ｃｃ．２２８８５　（６５）Ｂｅｃｋｅ，Ａ．Ｄ．　Ｃｈｅｍ．　１９８８，８８，２５４７．ｄｏｉ：　１０．１０６３／１　４５４０３３　（６６）Ｒｉｃｈａｒｄ，Ｆ；Ｂａｄｅｒ，Ｒ．Ａｔｏｍｓ　ｉｎ　Ｍｏｌｅｃｕｌｅｓ：ａ　Ｑｕａｎｔｕｍ　Ｔｈｅｏｒｙ；　Ｏｘｆｏｒｄ　Ｕｎｉｖｅｒｓｉｔｙ　Ｐｒｅｓｓ，Ｏｘｆｏｒｄ，ＵＫ，１９９０．　（６７）Ｈｉｒｓｈｆｅｌｄ，Ｆ．Ｌ．Ｔｈｅｏｒ　Ｃｈｉｍ　Ａｃｔａ　１９７７，４４，１２９．　ｄｏｉ：１０，ｌ００７／ｂｆ００５４９０９６　（６８）Ｚｈｏｕ，Ｘ．　；Ｒｏｎｇ，Ｃ．５Ｌ；Ｌｕ，Ｔ；Ｌｉｕ，Ｓ．Ｂ　Ａｃｔａ尸协　－Ｃｈｉｍ　Ｓｉｎ　２０１４．３Ｏ．２０５５．ｄｏｉ：１Ｏ．３８６６／ＰＫＵ．ＷＨＸＢ２０１４０９１９３　ｔ６９）Ｌｉｕ，Ｓ．Ｂ．Ｊ　Ｐｈｙｓ　Ｃｈｅｍ　Ａ　２０１５，１１９，３１０７．　ｄｏｉ：１０．１０２１／ａｃｓ．ｊｐｃａ．５ｂ００４４３　（７０）Ｋｏｈｎ，ｗ；Ｓｈａｍ，Ｌ．Ｊ．　Ｒｅｖ．１９６５，１４０，ＡＩ　１３３．　ｄｏｉ：１０．１１０３／ＰｈｙｓＲｅｖ．１４０．Ａ１１３３　ＮＯ．６　ｄｏｉ：１０．３８６６／ＰＫＵ．ＷＨＸＢ２０１７１０２３１　６４９　：Ｒｏｎｇ，Ｃ．Ｙ；Ｙｉｎｇ，Ｄ．Ｌ．；Ｌｉｕ，Ｓ．Ｂ．　Ｔｈｅｏｎ　Ｃｏｍｐｕｔ．　（７１）　Ｚｈａｏ，Ｄ．Ｂ．Ｃｈｅｍ．２０１３，１２，１３５００３４．ｄｏｉ：１０．１１４２￣０２１９６３３６１３５００３４ｘ　；Ｗｕ，Ｚ．Ｍ．；Ｒｏｎｇ，Ｃ．Ｙ；Ｌｕ，　；Ｈｕａｎｇ，Ｙ；Ｌｉｕ，Ｓ．Ｂ．　（７２）　Ｗｕ，Ｗ．Ｊ．１４９６．ｄｏｉ：ｌ０．１０３９／（２６ＣＰ０６３７６Ｄ　（７６）Ｚｈｏｕ，Ｘ．Ｙ；Ｒｏｎｇ，Ｃ．Ｙ；Ｌｕ，工；Ｚｈｏｕ，Ｅ　Ｐ．；Ｌｉｕ，Ｓ．Ｂ．　．　Ｃｈｅｍ．Ａ　２０１６，１２０，３６３４．ｄｏｉ：１０．１０２１／ａｃｓ．ｊｐｃａ．６ｂ０１１９７　（７７）Ｎｏｏｒｉｚａｄｅｈ，Ｓ．；Ｓｈａｋｅｒｚａｄｅｈ，Ｅ．Ｐｈｙｓ．Ｃｈｅｍ　Ｃｈｅｍ．Ｐｈｙｓ　２０１０，　１２，４７４２．ｄｏｉ：１０．１０３９／Ｂ９１６５０９Ｆ　Ｐｈｙｓ．Ｃｈｅｍ．Ａ　２０１５．１１９，８２１６．ｄｏｉ：１０．１０２ｌ／ａｃｓｄｐｃａ．５１３０４３０９　；Ｒｏｎｇ，Ｃ．Ｙ；Ｌｕ，　；Ａｙｅｒｓ，Ｅ、Ⅳ．；Ｌｉｕ，Ｓ．Ｂ．　（７３）　Ｗｕ，Ｚ．Ｍ．．　Ｃｈｅｍ．Ｃｈｅｍ．Ｐｈｙｓ．２０１５，　７，２７０５２．ｄｏｉ：１０．１０３９／Ｃ５ＣＰ０４４４２Ａ　ｆ７８）Ｓｎｙｄｅ￣Ｇ．Ｊ．　Ｐｈｙｓ．Ｃｈｅｍ．Ａ　２０１２，１１６，５２７２．　（７４）　Ｈｕａｎｇ，Ｙ；Ｒｏｎｇ，Ｃ．Ｙ；Ｚｈａｎｇ，Ｒ．Ｑ．；Ｌｉｕ，Ｓ．Ｂ．　Ｍｏ１．Ｍｏｄｅ１．　２０１６，２３，３．ｄｏｉ：１０．１００７／ｓ００８９４－０１６—３１７５一ｘ　ｕ，Ｓ．Ｂ．；Ｒｏｎｇ，Ｃ．Ｙ；Ｌｕ，　Ｐｈｙｓ．Ｃｈｅｍ．Ｃｈｅｍ　ｆ７５１　Ｌｉ．２０１７，１９，　ｄｏｉ：１０．１０２１￣ｐ３０４０１５ｋ　（７９）Ｈａｌｔｏｎ，Ｂ．Ｅｕｒ．　Ｏｒｇ，Ｃｈｅｍ．２００５，２００５，３３９１．　ｄｏｉ：１０．１００２／ｅｊＯＣ．２００５００２３１