Computational Physics (physics.comp-ph)

  • PDF
    The identification of meaningful reaction coordinates plays a key role in the study of complex molecular systems whose essential dynamics is characterized by rare or slow transition events. In a recent publication, the authors identified a condition under which such reaction coordinates exist - the existence of a so-called transition manifold - and proposed a numerical method for their point-wise computation that relies on short bursts of MD simulations. This article represents an extension of the method towards practical applicability in computational chemistry. It describes an alternative computational scheme that instead relies on more commonly available types of simulation data, such as single long molecular trajectories, or the push-forward of arbitrary canonically-distributed point clouds. It is based on a Galerkin approximation of the transition manifold reaction coordinates, that can be tuned to individual requirements by the choice of the Galerkin ansatz functions. Moreover, we propose a ready-to-implement variant of the new scheme, that computes data-fitted, mesh-free ansatz functions directly from the available simulation data. The efficacy of the new method is demonstrated on a realistic peptide system.
  • PDF
    In diffusion Monte Carlo (DMC) methods, the nodes (or zeroes) of the trial wave function dictate the magnitude of the fixed-node (FN) error. Within standard DMC implementations, they emanate from short multideterminant expansions, \textitstochastically optimized in the presence of a Jastrow factor. Here, following a recent proposal, we follow an alternative route by considering the nodes of selected configuration interaction (sCI) expansions built with the CIPSI (Configuration Interaction using a Perturbative Selection made Iteratively) algorithm. In contrast to standard implementations, these nodes can be \textitsystematically and \textitdeterministically improved by increasing the size of the sCI expansion. The present methodology is used to investigate the properties of the transition metal sulfide molecule FeS. This apparently simple molecule has been shown to be particularly challenging for electronic structure theory methods due to the proximity of two low-energy quintet electronic states of different spatial symmetry. In particular, we show that, at the triple-zeta basis set level, all sCI results --- including those extrapolated at the full CI (FCI) limit --- disagree with experiment, yielding an electronic ground state of $^{5}\Sigma^+$ symmetry. Performing FN-DMC simulation with sCI nodes, we show that the correct $^{5}\Delta$ ground state is obtained if sufficiently large expansions are used. Moreover, we show that one can systematically get accurate potential energy surfaces and reproduce the experimental dissociation energy as well as other spectroscopic constants.

Recent comments

Chris Ferrie May 01 2017 12:52 UTC

I've posted a public referee report on this paper here:

Salvatore Mandrà Nov 21 2014 13:50 UTC

The manuscript has been widely revised to focus the reader's attention on the proposed method and its application in presence of local disorder.


Salvatore, Gian Giacomo and Alán

Salvatore Mandrà Aug 01 2014 19:11 UTC

Thanks Dr. Hastings for your comment. It is true that the transverse field Ising model does not satisfy the requirements to apply our method with an exponential reduction. Indeed, the opposite would be quite impressive since the random Ising model is a NP-Hard problem and we ourselves would be prett

Matt Hastings Aug 01 2014 16:43 UTC

The "quite general" conditions seem not to include the transverse field Ising model, the subject of most of the intensive numerical work previously. Incidentally, the terms "Lanczos" and "Krylov subspace" might be helpful.

Steve Flammia Jan 10 2014 08:01 UTC

This paper is an important first step towards making numerical and computational studies more easily checkable and reproducible. Very interesting stuff.

Noon van der Silk Jan 10 2014 16:26 UTC

Interestingly, another thematically-similar paper popped up today - (posting the link here as I assume not many people are subbed to cs.DL)