All simulations discussed herein were performed using LAMMPS. All simulations were conducted using previously established methods and methodological choices. All simulations used periodic boundary conditions, and a cutoff of 15.0 Å for Lennard-Jones interactions. The SDK model lacks electrostatics. The equations of motion for the NVE simulations were evolved using the velocity-Verlet formulation. Isothermal conditions were realized using four-member Nosé-Hoover chain thermostats 27 with 1-ps coupling constants. All NPT n-C 720 H 1442 simulations were conducted using isotropic barostats except the infimer NPT simulations and the NPT 270 K nucleation simulations, which used anisotropic barostats. Pressure was regulated using 4-member barostat chains with 10-ps coupling constants while particle positions, momenta, and simulation cell sizes were evolved according to the scheme of Shinoda and coworkers. A time step of 5.0 fs was used for all simulations except: the infimer simulations at 270 K and 300 K, the nucleation simulations at 270 K and the NVE nucleation simulations. These simulations used a time step of 2.5 fs, which yielded a drift rate of less than 0.05%/100 million time steps in system total energy during tests at 300 K. For the nucleation simulations, system properties and configurations were saved every 10,000 and 40,000 time steps, respectively.