Diffusion Monte Carlo Study of Para -Diiodobenzene Polymorphism Revisited

DSpace/Manakin Repository

Diffusion Monte Carlo Study of Para -Diiodobenzene Polymorphism Revisited

Citable link to this page


Title: Diffusion Monte Carlo Study of Para -Diiodobenzene Polymorphism Revisited
Author: Hongo, Kenta; Watson, Mark A.; Iitaka, Toshiaki; Maezono, Ryo; Aspuru-Guzik, Alan

Note: Order does not necessarily reflect citation order of authors.

Citation: Hongo, Kenta, Mark A. Watson, Toshiaki Iitaka, Alán Aspuru-Guzik, and Ryo Maezono. 2015. “ Diffusion Monte Carlo Study of Para -Diiodobenzene Polymorphism Revisited .” Journal of Chemical Theory and Computation 11 (3) (March 10): 907–917. doi:10.1021/ct500401p.
Full Text & Related Files:
Abstract: We revisit our investigation of the diffusion Monte Carlo (DMC) simulation of p-DIB molecular crystal polymorphism. [J. Phys. Chem. Lett. 2010, 1, 1789-1794] We perform, for the first time, a rigorous study of finite-size effects and choice of nodal surface on the prediction of polymorph stability in molecular crystals using fixed-node DMC. Our calculations are the largest which are currently feasible using the resources of the K computer and provide insights into the formidable challenge of predicting such properties from first principles. In particular, we show that finite-size effects can influence the trial nodal surface of a small (1×1×1) simulation cell considerably. We therefore repeated our DMC simulations with a 1×3×3 simulation cell, which is the largest such calculation to date. We used a DFT nodal surface generated with the PBE functional and we accumulated statistical samples with ∼6.4×105 core-hours for each polymorph. Our final results predict a polymorph stability consistent with experiment, but indicate that results in our previous paper were somewhat fortuitous. We analyze the finite-size errors using model periodic Coulomb (MPC) interactions and kinetic energy corrections, according to the CCMH scheme of Chiesa, Ceperley, Martin, and Holzmann. We investigate the dependence of the finite-size errors on different aspect ratios of the simulation cell (k-mesh convergence) in order to understand how to choose an appropriate ratio for the DMC calculations. Even in the most expensive simulations currently possible, we show that the finite size errors in the DMC total energies are far larger than the energy difference between the two polymorphs, although error cancellation means that the polymorph prediction is accurate. Finally, we found that the T-move scheme is essential for these massive DMC simulations in order to circumvent population explosions and large time-step biases.
Published Version: doi:10.1021/ct500401p
Terms of Use: This article is made available under the terms and conditions applicable to Open Access Policy Articles, as set forth at http://nrs.harvard.edu/urn-3:HUL.InstRepos:dash.current.terms-of-use#OAP
Citable link to this page: http://nrs.harvard.edu/urn-3:HUL.InstRepos:24820079
Downloads of this work:

Show full Dublin Core record

This item appears in the following Collection(s)


Search DASH

Advanced Search