Anastasios Tsourtis

Anastasios Tsourtis, Vagelis Harmandaris, Dimitrios Tsagkarogiannis

We present a systematic coarse-graining (CG) strategy for many particle molecular systems based on cluster expansion techniques. We construct a hierarchy of coarse-grained Hamiltonians with interaction potentials consisting of two, three and higher body interactions. The accuracy of the derived cluster expansion based on interatomic potentials is examined over a range of various temperatures and densities and compared to direct computation of pair potential of mean force. The comparison of the coarse-grained simulations is done on the basis of the structural properties, against the detailed all-atom data. We give specific examples for methane and ethane molecules in which the coarse-grained variable is the center of mass of the molecule. We investigate different temperature and density regimes, and we examine differences between the methane and ethane systems. Results show that the cluster expansion formalism can be used in order to provide accurate effective pair and three-body CG potentials at high $T$ and low $ρ$ regimes. In the liquid regime the three-body effective CG potentials give a small improvement, over the typical pair CG ones; however in order to get significantly better results one needs to consider even higher order terms.

Anastasios Tsourtis, Yannis Pantazis, Ioannis Tsamardinos

Inferring the driving equations of a dynamical system from population or time-course data is important in several scientific fields such as biochemistry, epidemiology, financial mathematics and many others. Despite the existence of algorithms that learn the dynamics from trajectorial measurements there are few attempts to infer the dynamical system straight from population data. In this work, we deduce and then computationally estimate the Fokker-Planck equation which describes the evolution of the population's probability density, based on stochastic differential equations. Then, following the USDL approach, we project the Fokker-Planck equation to a proper set of test functions, transforming it into a linear system of equations. Finally, we apply sparse inference methods to solve the latter system and thus induce the driving forces of the dynamical system. Our approach is illustrated in both synthetic and real data including non-linear, multimodal stochastic differential equations, biochemical reaction networks as well as mass cytometry biological measurements.