Ab initio stochastic optimization of conformational and many-body degrees of freedom

Moderate to large size molecules in solution have complex energy surfaces due to intramolecular (conformational) and intermolecular (many-body) interactions. The first principles Monte Carlo (FPMC) method, previously shown to effectively locate minimum-energy structures for systems with only many-body complexity, has been extended to address conformational flexibility by adding three new Monte Carlo move types. The primary advantage of the FPMC method is the ability to efficiently locate minimum energy structures of molecules with conformational flexibility in the presence of explicit solvent molecules using highly accurate quantum chemical calculations. The additions to FPMC were validated by studying conformers of glycerol, glyceraldehyde, and a large humic acid monomer unit. The structure of glyceraldehyde in the presence of one and two water molecules was also explored to demonstrate the power of FPMC to study systems with both conformational and many-body degrees of freedom.