The simulation force field details and methodology are the same as we previously reported. Simulation compositions are reported in Table . Simulations were conducted with the GROMACS 2016.2 software package. Force field details for DMDBPMA and dodecane are reported in Ref. 61. Initial configurations were created with Packmol, followed by energy minimization using a steepest descent algorithm. Each simulation was equilibrated for 5 ns in the NPT ensemble, followed by 20 ns of equilibration in the NVT ensemble and 50 ns of production in the NVT ensemble. In the NPT ensemble, pressure was set to 1 bar with the Berendsen barostat 76 using a 2 ps coupling time and temperature was set to 300 K using a 0.2 ps coupling time with the velocity rescale thermostat. In the NVT ensemble, temperature was set to 300 K with the Nosé-Hoover thermostat 77 with a 0.2 ps coupling time. Dynamics were conducted using a 2 fs time step using the leap-frog Verlet integrator, with hydrogen-containing bonds constrained using the LINCS algorithm. A 15 Å cutoff was used for Lennard-Jones and short-range electrostatic interactions, with particle-mesh Ewald summation used for long-range electrostatics. Trajectories were sampled for analysis at 100 ps intervals.