A new method to calculate broadband dielectric spectra of solvents from molecular dynamics simulations demonstrated with polarizable force fields.

Simulating the dielectric spectra of solvents requires the nuanced definition of inter- and intra-molecular forces. Non-polarizable force fields, while thoroughly benchmarked for dielectric applications, do not capture all the spectral features of solvents, such as water. Conversely, polarizable force fields have been largely untested in the context of dielectric spectroscopy but include charge and dipole fluctuations that contribute to intermolecular interactions. We benchmark non-polarizable force fields and the polarizable force fields AMOEBA03 and HIPPO for liquid water and find that the polarizable force fields can capture all the experimentally observed spectral features with varying degrees of accuracy. However, the non-polarizable force fields miss at least one peak. To diagnose this deficiency, we decompose the liquid water spectra from polarizable force fields at multiple temperatures into static and induced dipole contributions and find that the peak originates from induced dipole contributions. Broadening our inquiry to other solvents parameterized with the AMOEBA09 force field, we demonstrate good agreement between the experimental and simulated dielectric spectra of methanol and formamide. To produce these spectra, we develop a new computational approach to calculate the dielectric spectrum via the fluctuation dissipation theorem. This method minimizes the error in both the low and high frequency portions of the spectrum, improving the overall accuracy of the simulated spectrum and broadening the computed frequency range.

CryoEM reveals the complex self-assembly of a chemically driven disulfide hydrogel.

Inspired by the adaptability of biological materials, a variety of synthetic, chemically driven self-assembly processes have been developed that result in the transient formation of supramolecular structures. These structures form through two simultaneous reactions, forward and backward, which generate and consume a molecule that undergoes self-assembly. The dynamics of these assembly processes have been shown to differ from conventional thermodynamically stable molecular assemblies. However, the evolution of nanoscale morphologies in chemically driven self-assembly and how they compare to conventional assemblies has not been resolved. Here, we use a chemically driven redox system to separately carry out the forward and backward reactions. We analyze the forward and backward reactions both sequentially and synchronously with time-resolved cryogenic transmission electron microscopy (cryoEM). Quantitative image analysis shows that the synchronous process is more complex and heterogeneous than the sequential process. Our key finding is that a thermodynamically unstable stacked nanorod phase, briefly observed in the backward reaction, is sustained for ∼6 hours in the synchronous process. Kinetic Monte Carlo modeling show that the synchronous process is driven by multiple cycles of assembly and disassembly. The collective data suggest that chemically driven self-assembly can create sustained morphologies not seen in thermodynamically stable assemblies by kinetically stabilizing transient intermediates. This finding provides plausible design principles to develop and optimize supramolecular materials with novel properties.