axsk / sqracore.jl

Approximation of the generator of a Diffusion process by pointwise potential evaluations.

SqraCore.jl

Approximation of the generator of a Diffusion process by pointwise potential evaluations.

Theory

Using the Square-Root approximation (SQRA) as in [6]:

The generetor of the Koopman operator of the diffusion in a potential $U$: $dX = -\nabla U dt + \sigma dB_t$ on a Voronoi tesselation of the state space with centers $x_i$ is approximated by the matrix $Q$:

$$Q_{ij} = C_{ij} \sqrt\frac{\pi_j}{\pi_i}$$ for $i\neq j$ and $Q_{ii} = -\sum_{j\neq i} Q_{ij}$ with $$C_{ij} = \beta^{-1}\frac{A_{ij}}{h_{ij} V_i}$$ where $A_{ij}$ and $V_{ij}$ are the (common boundary) areas and volumes of the cells, $h_{ij}$ the distance between the centers of cells $i,j$.

It ows its name to the fact that it is proportional to the square-root of the fractions of the stationary distribution $\pi_i$ at the cell centers, $$\pi_i = \exp (-\beta U(x_i)).$$

References

• [1] Weber (2010). "A subspace approach to Molecular Markov State Models via a New Infinitesimal Generator"
• [2] Lie, Fackeldey, M. Weber (2013)
• [4] Donati, Heida, Keller, Weber (2018). "Estimation of the infinitesimal generator by square-root approximation."
• [3] Donati, Weber, Keller (2021). "Markov models from the square root approximation of the Fokker–Planck equation: calculating the grid-dependent flux."
• [5] Donati, Weber, Keller (2022). "A review of Girsanov reweighting and of square root approximation for building molecular Markov state models."
• [6] Schütte, Klus, Hartmann (2022). "Overcoming the Timescale Barrier in Molecular Dynamics: Transfer Operators, Variational Principles, and Machine Learning."

sqracore.jl's Issues

TagBot trigger issue

This issue is used to trigger TagBot; feel free to unsubscribe.

If you haven't already, you should update your TagBot.yml to include issue comment triggers.
Please see this post on Discourse for instructions and more details.

If you'd like for me to do this for you, comment TagBot fix on this issue.
I'll open a PR within a few hours, please be patient!

Scale Q with inv(beta)