Related Experiment Video
Updated: May 2, 2026

Multiscale Sampling of a Heterogeneous Water/Metal Catalyst Interface using Density Functional Theory and Force-Field Molecular Dynamics
Published on: April 12, 2019
Improving the efficiency of Monte Carlo surface hopping calculations
1Department of Chemistry, Tulane University , New Orleans, Louisiana 70118, United States.
This study enhances surface hopping semiclassical calculations using Monte Carlo methods. Averaging trajectory sections and employing Sobol sequences improve numerical efficiency and accuracy with fewer samples.
Area of Science:
- Quantum chemistry
- Computational physics
- Chemical dynamics
Background:
- Semiclassical methods are crucial for simulating quantum dynamics.
- Surface hopping is a common technique to model non-adiabatic transitions.
- Improving the efficiency and accuracy of these methods is an ongoing challenge.
Purpose of the Study:
- To investigate a Monte Carlo-based surface hopping method for semiclassical calculations.
- To enhance the numerical efficiency of the surface hopping approach.
- To assess the impact of trajectory sampling techniques on accuracy.
Main Methods:
- A surface hopping method incorporating a Monte Carlo procedure for hop decisions.
- Averaging over multiple copies of trajectory sections within interaction regions.
- Utilizing Sobol sequences for selecting initial trajectory momentum.
Main Results:
- Numerical efficiency is improved by averaging trajectory sections.
- Accurate results are achievable with reduced trajectory sample sizes.
- Sobol sequences demonstrate effectiveness in trajectory initializations.
Conclusions:
- The proposed Monte Carlo surface hopping method offers improved efficiency.
- Averaging trajectory segments and using Sobol sequences are effective strategies.
- This approach enables accurate semiclassical calculations with smaller computational cost.
Related Concept Videos
Mechanistic Models: Compartment Models in Algorithms for Numerical Problem Solving
In individual population analyses, different algorithms are employed, such as Cauchy's method, which uses a...
Maxwell-Boltzmann Distribution: Problem Solving
This distribution function f(v) is defined by saying that the expected number N (v1,v2) of particles with speeds between v1 and v2 is given by
Energy Conservation and Bernoulli's Equation
All the terms in the equation have the dimension of energy per unit volume. The kinetic energy per unit volume is called the kinetic energy density, and the potential energy per unit volume is...
Numerical Calculations
The solution to a problem is obtained using different methods. While manually solving algebraic symbols is one of the most common methods, the graphical method is often preferred. Computers...
Reaction Mechanisms: Rate-limiting Step Approximation
Methods of Medium Optimization

