Simulating Hydrogen-Bond Structure and Dynamics in Glassy Solids Composed of Imidazole Oligomers

详细信息    查看全文

• 作者：Jacob A. Harvey ; Scott M. Auerbach
• 刊名：Journal of Physical Chemistry B
• 出版年：2014
• 出版时间：July 10, 2014
• 年：2014
• 卷：118
• 期：27
• 页码：7609-7617
• 全文大小：495K
• ISSN：1520-5207

文摘

We simulated structural and dynamical properties of imidazoles tethered to aliphatic backbones to determine how chain length influences the competition between extended hydrogen-bond networks and imidazole reorientation dynamics. We performed molecular dynamics simulations on hypothetical solids using the GAFF Amber force field over the temperature range 300鈥?00 K, for chain lengths varying from monomers to pentamers. We investigated the effect of heterogeneity by simulating monodisperse and polydisperse solids with the same average chain length. We computed hydrogen-bond cluster sizes and percolation ratios; orientational order parameters associated with imidazole rings, tethering linkers, and backbones; and orientational correlation functions for imidazole rings. We found the surprising result that chain-length heterogeneity negligibly affects system density at standard pressure, the fraction of percolating hydrogen-bonded clusters, and the order parameters for backbone, linker, and imidazole ring. Decreasing oligomer chain length from pentamers to shorter chains decreases the tendency to form percolating hydrogen-bond networks while dramatically decreasing imidazole ring reorientation times from a broad range of 100鈥?00 ps for pentamers down to 20 ps for monomers, hence quantifying the competition between hydrogen-bond cluster size and reorientation rate. The computed orientational order parameters suggest the following hierarchy of structural excitations: imidazole ring reorientation in the range 400鈥?00 K, linker motion around 500鈥?00 K, and backbone relaxation at 600鈥?00 K in this model. The question remains for this class of systems which of these motions is crucial for facile proton transport.