Workshops and programs
 January 25th27th – Piecewise Deterministic Markov Processes and Sampling – ERC is cofunding this event (http://wikimath.univmlv.fr/pdmp/doku.php/events:pdmp_marne_2017)
 Maths/industry workshop on numerical simulation in the pharmaceutical industry (Friday, April 28th, at Institut Henri Poincaré, Paris)
 IHP Trimester Stochastic Dynamics Out of Equilibrium: a program on the study of nonequilibrium dynamics in statistical physics was held in Spring 2017 at Institut Henri Poincaré in Paris. Several events emphasized the importance of developing dedicated numerical methods for more efficiently simulating nonequilibrium systems, in particular introductory lectures during the preschool at CIRM and a workshop on numerical methods for nonequilibrium systems. See here for the general program page. The ERC grant MSMath is cofunding this event, which is coorganized by Gabriel Stoltz.
 IPAM Long Program on “Complex HighDimensional Energy Landscapes”, September 11th – December 15th 2017. Tony Lelièvre is one of the organizers with Cecilia Clementi, Graeme Henkelman, Richard Hennig, Mitchell Luskin, Noa Marom, Petr Plechac and Christof Schuette.
 ICTS program on “Large deviation theory in statistical physics: Recent advances and future challenges”, August 14th – October 13th 2017. Tony Lelièvre is one of the organizers with Arvind Ayyer, Frank den Hollander, Abhishek Dhar, Juan P. Garrahan, Christopher Jarzynski, Manjunath Krishnapur, Sanjib Sabhapandit and Hugo Touchette.
Working Group

January 10th – Spectral theory for random Poincaré maps (M. Baudel).
 January 30th Fitting nonbonded parameters from and for condensed phase simulations (F. Hedin).
 The quality of atomistic simulations depends decisively on the accuracy of the underlying energy function (force field). Of particular importance for condensedphase properties are nonbonded interactions, including the electrostatic and LennardJones
terms. Permanent atomic multipoles (MTPs) are an extension to common pointcharge (PC) representations in atomistic simulations. MTPs are commonly determined from and fitted to an ab initio Electrostatic Potential (ESP), and LennardJones (LJ) parameters are obtained from comparison of experimental and computed observables using molecular dynamics (MD) simulations. For this a set of thermodynamic observables such as density, heat of vaporization, and hydration free energy is chosen, to which the parametrization is fitted.
A toolkit for optimizing nonbonded interactions for atomistic force fields of different qualities will be presented. It supports fitting of standard PCbased force fields and more physically motivated multipolar (MTP) force fields. Validation was done with a set of 20 molecules ranging from Nmethylacetamide and benzene to halogenated benzenes, phenols, anilines, and pyridines. ( http://pubs.acs.org/doi/abs/10.1021/acs.jcim.6b00280 )  Investigating the stability of solvated Hæmoglobin Tetramers using Molecular Dynamics (Florent Hedin ).
Hæmoglobin is a metalloprotein (containing iron), in charge of oxygen transport in red blood cells of animals. It transports oxygen from the lungs to the rest of the body where it releases the oxygen for cell use. Human hæmoglobin is a tetrameric protein consisting of two α and two β subunits.
Each subunit contains a heme group at the center to which molecular oxygen or other ligands bind: the ligandbound (oxy) state is identified as Rstate or 2DN2, and the ligandfree (deoxy) state as T state or 2DN3. The distance between the two terminal Histidines of the β chains is characteristic of each state: this distance fluctuates between 10 – 15 Å for the compact oxy 2DN3 state, and 30 – 35 Å for the more extended deoxy 2DN2 state. Although the deoxy state should be thermodynamically stable and remain in an extended state (as it is the case experimentally), several previous MD studies 1 found that the tetramers adopt a compact structure, similar to the oxy one (i.e. a T ↔ R transition was found).The stability of the tetramer in water boxes of increasing size is investigated, over several hundreds of nanoseconds, using the CHARMM ForceField and the GROMACS software. A coarsegrained definition of the water density is used for analyzing instantaneous fluctuations of th density at the proteinsolvent interface.
 The quality of atomistic simulations depends decisively on the accuracy of the underlying energy function (force field). Of particular importance for condensedphase properties are nonbonded interactions, including the electrostatic and LennardJones
 February 9th – On the coarsegraining of Hamiltonian dynamics (V. Ehrlacher)
 February 27th
 Letif Mones, de Cambridge et Warwick, Gaussian process regression in potential of mean force calculations.
 Stefan Chmiela (Institute of Software Engineering and Theoretical Computer Science, Berlin)
 Résumé : Using conservation of energy – a fundamental property of closed classical and quantum mechanical systems – we develop an efficient gradientdomain machine learning (GDML) approach to construct accurate molecular force fields using a restricted number of samples from ab initio molecular dynamics (AIMD) trajectories. The GDML implementation is able to reproduce global potentialenergy surfaces of intermediatesize molecules with an accuracy of 0.3 kcal/mol for energies and 1 kcal/mol/A for atomic forces using only 1000 conformational geometries for training. We demonstrate this accuracy for AIMD trajectories of molecules, including benzene, toluene, naphthalene, ethanol, uracil, and aspirin. The challenge of constructing conservative force fields is accomplished in our work by learning in a Hilbert space of vectorvalued functions that obey the law of energy conservation. The GDML approach enables quantitative molecular dynamics simulations for molecules at a fraction of cost of explicit AIMD calculations, thereby allowing the construction of efficient force fields with the accuracy and transferability of highlevel ab initio methods.
 Résumé : Using conservation of energy – a fundamental property of closed classical and quantum mechanical systems – we develop an efficient gradientdomain machine learning (GDML) approach to construct accurate molecular force fields using a restricted number of samples from ab initio molecular dynamics (AIMD) trajectories. The GDML implementation is able to reproduce global potentialenergy surfaces of intermediatesize molecules with an accuracy of 0.3 kcal/mol for energies and 1 kcal/mol/A for atomic forces using only 1000 conformational geometries for training. We demonstrate this accuracy for AIMD trajectories of molecules, including benzene, toluene, naphthalene, ethanol, uracil, and aspirin. The challenge of constructing conservative force fields is accomplished in our work by learning in a Hilbert space of vectorvalued functions that obey the law of energy conservation. The GDML approach enables quantitative molecular dynamics simulations for molecules at a fraction of cost of explicit AIMD calculations, thereby allowing the construction of efficient force fields with the accuracy and transferability of highlevel ab initio methods.
 Letif Mones, de Cambridge et Warwick, Gaussian process regression in potential of mean force calculations.
 March 14th – About “Grandes déviations à l’équilibre pour des systèmes invariants par translation” (J. Reygner).
 April 5th – About “Coarsegraining limits in FokkerPlanck equations: a variational approach inspired by large deviations ” (U. Sharma)
 May 11, 02:00 pm, at CERMICS (room F102). Laura Joana Silva Lopes will present : Méthodes pour obtenir le chemin réactif préférentiel en dynamique moléculaire.
 May 19th, 11:00 am, at CERMICS (B211 seminar room). Frédéric Lavancier (Laboratoire Jean Leray, Nantes), will present : processus ponctuels determinantaux
 May 23th, 10;00 am, at CERMICS (B211 seminar room). Erwan Scornet (CMAP) will present: Walk on random forests.
Abstract: The recent and ongoing digital world expansion now allows anyone to have access to a tremendous amount of information. However collecting data is not an end in itself and thus techniques must be designed to gain indepth knowledge from these large data bases. This has led to a growing interest for statistics, as a tool to find patterns in complex data structures, and particularly for turnkey algorithms which do not require specific skills from the user.
Such algorithms are quite often designed based on a hunch without any theoretical guarantee. Indeed, the overlay of several simple steps (as in random forests or neural networks) makes the analysis more arduous. Nonetheless, the theory is vital to give assurance on how algorithms operate thus preventing their outputs to be misunderstood. Among the most basic statistical properties is the consistency which states that predictions are asymptotically accurate when the number of observations increases. In this talk, I will present a first result on Breiman’s forests consistency and show how it sheds some lights on its good performance in a sparse regression setting.

October 5th, 10;00 am, at CERMICS (B211 seminar room). Grégoire Ferré (CERMICS) will present: On the discretization of FeynmanKac semigroups
Abstract : In this presentation, I will present a framework for the numerical analysis in the discretization time step of FeynmanKac semigroups. I will first explain why these semigroups naturally appear in several fields, in particular large deviation theory and Diffusion Monte Carlo (DMC). Then, I will present elements of analysis of the error on the invariant measure that is done when discretizing a stochastic differential equation. Finally, I will show how we extended these results to the case of FeynmanKac semigroups, and why this analysis is interesting in general to discretize non probabilityconserving processes, reviewing some theoretical and practical open problems suggested by this work.. This analysis is supported by relevant numerical applications.
 November 9th, 10am at CERMICS . Laurent Michel (Université de Nice).
 November 23rd, 10am at CERMICS. N. Nuksen.
 December 14th, 10am, A. Schlichting.