**Are you an EPFL student looking for a semester project?**

Work with us on data science and visualisation projects, and deploy your project as an app on top of GraphSearch.

Person# Marcin Rafal Buchowiecki

This page is automatically generated and may contain information that is not correct, complete, up-to-date, or relevant to your search query. The same applies to every other page on this website. Please make sure to verify the information with EPFL's official sources.

Related units

Loading

Courses taught by this person

Loading

Related research domains

Loading

Related publications

Loading

People doing similar research

Loading

Related research domains (1)

Related publications (3)

Courses taught by this person

Monte Carlo method

Monte Carlo methods, or Monte Carlo experiments, are a broad class of computational algorithms that rely on repeated random sampling to obtain numerical results. The underlying concept is to use rando

No results

Loading

Loading

Loading

Related units (1)

People doing similar research (127)

Marcin Rafal Buchowiecki, Jing Huang, Jiri Vanicek

The primary H/D kinetic isotope effect on the intramolecular proton transfer in malonaldehyde is determined from quantum instanton path integral Monte Carlo simulations on a fully dimensional and validated potential energy surface for temperatures between 250 and 1500 K. Our calculations, based on thermodynamic integration with respect to the mass of the transferring particle are significantly accelerated by the direct evaluation of the kinetic isotope effect instead of computing it as a ratio of two rate constants. At room temperature, the KIE from the present simulations is 5$\pm$4. The KIE is found to vary considerably as a function of temperature and the low-T behaviour is dominated by the fact that the free energy derivative in the reactant state increases more rapidly than at the transition state. Detailed analysis of the various contributions to the quantum rate constant together with estimates for rates from conventional transition state theory and from periodic orbit theory suggest that the KIE in malonaldehyde is dominated by zero point energy effects and that tunneling plays a minor role at room temperature

Marcin Rafal Buchowiecki, Jiri Vanicek

A general method for the direct evaluation of the temperature dependence of the quantum-mechanical reaction rate constant in many-dimensional systems is described. The method is based on the quantum instanton approximation for the rate constant, thermodynamic integration with respect to the inverse temperature, and the path integral Monte Carlo evaluation. It can describe deviations from the Arrhenius law due to the coupling of rotations and vibrations, zero-point energy, tunneling, corner-cutting, and other nuclear quantum effects. The method is tested on the Eckart barrier and the full-dimensional H+H2→H2+H reaction. In the temperature range from 300 to 1500 K, the error of the present method remains within 13% despite the very large deviations from the Arrhenius law. The direct approach makes the calculations much more efficient, and the efficiency is increased even further (by up to two orders of magnitude in the studied reactions) by using optimal estimators for reactant and transition state thermal energies. Which of the estimators is optimal, however, depends on the system and the strength of constraint in a constrained simulation.

Marcin Rafal Buchowiecki, Jiri Vanicek

The Feynman path integral approach for computing equilibrium isotope effects and isotope fractionation corrects the approximations made in standard methods, although at significantly increased computational cost. We describe an accelerated path integral approach based on three ingredients: the fourth- order Takahashi-Imada factorization of the path integral, thermodynamic integration with respect to mass, and centroid virial estimators for relevant free energy derivatives. While the frst ingredient speeds up convergence to the quantum limit, the second and third improve statistical convergence. The combined method is applied to compute the equilibrium constants for isotope exchange reactions H2+D=H+HD and H2+D2=2HD.

2013