Sensitivity analysis and optimization of reaction rate

The chemical reaction rate from reactant to product depends on the geometry of
potential energy surface (PES) as well as the temperature. We consider a design problem of how to choose the best PES from a given family of smooth potential functions in order to maximize (or minimize) the reaction rate for a given chemical reaction. By utilizing the transition-path theory, we relate reaction rate to committor functions which solves boundary-value elliptic problems, and
perform the sensitivity analysis of the underlying elliptic equations via adjoint approach. We derive the derivative of the reaction rate with respect to the potential function. The shape derivative with respect to the domains defining reactant and product is also investigated. The numerical optimization method based on the gradient is applied for two simple numerical examples to demonstrate the feasibility of our approach.