Alcohols are interesting organic compounds that are often used as reactants or solvents in chemical processes, but they also play various roles in many biological processes. As such, they are often the subject of various researches. In this work, we focused on hexan-1-ol, hexan-2-ol, hexan-3-ol, and 3,3-dimethylbutanol. Our goal was to investigate their properties at the molecular level using the computer simulations of molecular dynamics and a selected model based on the TraPPE–UA force field. We started the computer simulations with the energy minimization of the system, followed by the equilibration in the NVT ensemble, the equilibration in the NPT ensemble, and finished with the production step that provided the basic simulation results. By comparing the simulation results for some thermodynamic properties of the simulated system with the experimental results of alcohol from the literature, we investigated how well the chosen model describes the studied system. We found that the TraPPE–UA force field describes the studied primary alcohols slightly better than the secondary alcohols. Since we were also interested in the structural properties of the studied alcohols, we calculated the radial and spatial distribution functions from the simulation results and based on them the average number of hydrogen bonds per hydroxyl group of the alcohol. It turned out that the molecules of primary alcohols can form the most hydrogen bonds on average. Interestingly, the distribution of intramolecular distances from the C1 to the C4 atom in different model alcohols also showed that the rotation of the tert-butyl group around the bond between the C2 and C3 atoms in the model is apparently very limited.
|
