Non-equilibrium reaction and relaxation dynamics in a strongly interacting explicit solvent: F + CD3CN treated with a parallel multi-state EVB model.
Non-equilibrium reaction and relaxation dynamics in a strongly interacting explicit solvent: F + CD3CN treated with a parallel multi-state EVB model.
复制标题
强相互作用显式溶剂中的非平衡反应和弛豫动力学:用并行多态 EVB 模型处理的 F CD3CN。
DOI:
--
复制
发表时间:
2014
影响因子:
4.4
通讯作者:
J. Harvey
中科院分区:
文献类型:
--
作者:
D. Glowacki;A. Orr;J. Harvey
We describe a parallelized linear-scaling computational framework developed to implement arbitrarily large multi-state empirical valence bond (MS-EVB) calculations within CHARMM and TINKER. Forces are obtained using the Hellmann-Feynman relationship, giving continuous gradients, and good energy conservation. Utilizing multi-dimensional Gaussian coupling elements fit to explicitly correlated coupled cluster theory, we built a 64-state MS-EVB model designed to study the F + CD3CN → DF + CD2CN reaction in CD3CN solvent (recently reported in Dunning et al. [Science 347(6221), 530 (2015)]). This approach allows us to build a reactive potential energy surface whose balanced accuracy and efficiency considerably surpass what we could achieve otherwise. We ran molecular dynamics simulations to examine a range of observables which follow in the wake of the reactive event: energy deposition in the nascent reaction products, vibrational relaxation rates of excited DF in CD3CN solvent, equilibrium power spectra of DF in CD3CN, and time dependent spectral shifts associated with relaxation of the nascent DF. Many of our results are in good agreement with time-resolved experimental observations, providing evidence for the accuracy of our MS-EVB framework in treating both the solute and solute/solvent interactions. The simulations provide additional insight into the dynamics at sub-picosecond time scales that are difficult to resolve experimentally. In particular, the simulations show that (immediately following deuterium abstraction) the nascent DF finds itself in a non-equilibrium regime in two different respects: (1) it is highly vibrationally excited, with ∼23 kcal mol(-1) localized in the stretch and (2) its post-reaction solvation environment, in which it is not yet hydrogen-bonded to CD3CN solvent molecules, is intermediate between the non-interacting gas-phase limit and the solution-phase equilibrium limit. Vibrational relaxation of the nascent DF results in a spectral blue shift, while relaxation of the post-reaction solvation environment results in a red shift. These two competing effects mean that the post-reaction relaxation profile is distinct from what is observed when Franck-Condon vibrational excitation of DF occurs within a microsolvation environment initially at equilibrium. Our conclusions, along with the theoretical and parallel software framework presented in this paper, should be more broadly applicable to a range of complex reactive systems.
登录
查看更多内容
影响因子:
1.7
作者:
H. Werner;G. Knizia;F. Manby
通讯作者:
H. Werner;G. Knizia;F. Manby
影响因子:
2.3
作者:
Reed S
通讯作者:
Reed S
影响因子:
15
作者:
Glowacki, David R.;Liang, C. H.;Pilling, Michael J.
通讯作者:
Pilling, Michael J.
DOI:
10.1039/c2sc21267f
发表时间:
2013
期刊:
Chem. Sci.
影响因子:
--
作者:
Abou-Chahine F
通讯作者:
Abou-Chahine F
影响因子:
21.8
作者:
Otto, R.;Brox, J.;Wester, R.
通讯作者:
Wester, R.