Méthodes déterministes et stochastiques pour la simulation moléculaire, Deterministic and stochastic methods for molecular simulation

Kimiya Minoukadeh

Description

Sous la direction de Eric Cances
Thèse soutenue le 24 novembre 2010: Paris Est
La simulation moléculaire est un outil indispensable pour comprendre le comportement de systèmes complexes pour lesquels les expériences s'avèrent coûteuses ou irréalisables à l'heure actuelle. Cette thèse est dédiée aux aspects méthodologiques de la simulation moléculaire et comprend deux volets. Le premier volet porte sur la recherche de chemins de réaction et de points col d'une surface d'énergie potentielle. Nous proposons, dans le chaptire 3, une amélioration d'une des méthodes de cette classe, appelée 'Activation Relaxation Technique(ART). Nous donnons également une preuve de convergence pour un algorithme prototype. Le deuxieme volet porte sur le calcul d'énergie libre pour les transitions caractérisées par une coordonnée de réaction. Nous nous plaçons dans le cadre d'une méthode d'échantillonnage d'importance adaptative, appelée 'Adaptive Biasing Force' (ABF). Ce volet comprend en soi deux sous-parties. La première partie (chapitre 5) s'attache à montrer l'applicabilité à un système biomoléculaire, d'une nouvelle mise en oeuvre parallèle d'ABF, nommée 'multiple-walker ABF' (MW-ABF), consistant à utiliser plusieurs répliques. Cette mise en oeuvre s'est avérée utile pour surmonter des problèmes liés à un mauvais choix de coordonnée de réaction. Nous confirmons ensuite ces résultats numériques en étudiant la convergence théorique d'un algorithme d'ABF adapté. Le chapitre 6 comprend une étude de convergence en temps long utilisant les méthodes d'entropie relative et les inégalités de Sobolev logarithmiques
-Recherche d'états de transition
-Calcul d'énergie libre
-Inégalités de Sobolev logarithmiques
Molecular simulation is an essential tool in understanding complex chemical and biochemical processes as real-life experiments prove increasingly costly or infeasible in practice . This thesis is devoted to methodological aspects of molecular simulation, with a particular focus on computing transition paths and their associated free energy profiles. The first part is dedicated to computational methods for reaction path and transition state searches on a potential energy surface. In Chapter 3 we propose an improvement to a widely-used transition state search method, the Activation Relaxation Technique (ART). We also present a local convergence study of a prototypical algorithm. The second part is dedicated to free energy computations. We focus in particular on an adaptive importance sampling technique, the Adaptive Biasing Force (ABF) method. The first contribution to this field, presented in Chapter 5, consists in showing the applicability to a large molecular system of a new parallel implementation, named multiple-walker ABF (MW-ABF). Numerical experiments demonstrated the robustness of MW-ABF against artefacts arising due to poorly chosen or oversimplified reaction coordinates. These numerical findings inspired a new study of the longtime convergence of the ABF method, as presented in Chapter 6. By studying a slightly modified model, we back our numerical results by showing a faster theoretical rate of convergence of ABF than was previously shown
-Transition state search
-Free energy computations
-Logarithmic Sobolev inequalities
Source: http://www.theses.fr/2010PEST1034/document

Informations

• Publié par: Thesee
• Nombre de lectures: 28
• Langue: Français
• Poids de l'ouvrage: 2 Mo

Extrait

Denis
on
F
p
des
our
T
obtenir
cteur
le
cteur
titre
Examinateurs
de
olis)
Docteur
P
en
INRIA
Sciences
-
de
ec
l
y
'Université
P
P
-
aris-Est
-
Sp
Univ
écialité
Austin
:
Erw
Ma
ou
théma
Dir
tiques
Eric
École
des
Do
P
ctorale
Co-dir
:
T
MSTIC
-
Présen
on
tée
ec
et
Chris
souten
ersité
ue
ala
par
(Sophia-An
Kimiy
Henkelman
a
y
MINOUKADEH
at
Métho
ersité
des
aris-Dauphine
déterministes
an
et
a
sto
-
c
(Rennes)
hastiques
e
p
:
our
Cancès
la
École
sim
P
ulation
ts
moléculaire
arisT
Soutenance
h
le
e
24
:
no
on
v
Lelièvre
em
École
bre
P
2010
ts
dev
arisT
an
h
t
:
le
Chipot
jury
Univ
comp
Nancy
osé
T
de
y
R
INRIA
app
tip
orteurs
Graeme
:
-
Jean
ersit
Dolbea
of
ul
exas
t
Univ
-
Thèse
tel-00597694, version 1 - 1 Jun 2011tel-00597694, version 1 - 1 Jun 2011particular
is
del,
déterministes
an
et
al
stochastiques
ol
pour
is
la
e
simula
re
tion
n
moléculaire
methods
Résumé
exp
La
fo
sim
on
ulation
e
mol
The
éculaire
e
est
y
un
arising
outil
metho
indisp
rate
ensable
Molecular
p
bio
our
costly
comprendre
ects
le
asso
comp
rea
ortemen
e.
t
tran-
de
(AR
s
ergenc
ystèmes
n-
complexes
tanc
p
e
our
consists
lesquels
a
les
ts
exp
or
é
the
riences
By
s'a
results
v
w
èren
olecular
t
an
coûteuses
complex
ou
pro
irréalisables
pro
à
alistically
l'
metho
heure
ulation,
actuelle.
transition
Cette
proles.
thèse
computational
est
transition
dédiée
tial
aux
w
asp
t
ects
h
métho
Relaxation
dologiques
also
de
lo
la
a
sim
dedicated
ulation
e
moléculaire
e
et
hnique,
comprend
(ABF)
deux
to
v
ted
ol
wing
ets.
molecular
Le
tation,
premie
(MW-ABF).
r
of
v
o
olet
reaction
p
a
orte
ergence
s
ted
ur
tly
la
k
rec
wing
herc
ergence
he
stochastic
de
or
c
tion
hem
ulation
ins
tial
de
n
réac
hemical
t
hem-
i
as
on
n
et
e
de
r
p
This
oin
oted
ts
cal
col
molecular
d'une
a
surface
on
d'éne
and
rgie
free
p
rst
oten
dicated
tielle.
ds
Nous
path
prop
searc
osons,
p
dans
sur
le
Chapter
c
prop
hapitre
v
3
a
,
state
une
d,
amélioration
ctiv
d'une
ec
des
W
métho
re
des
t
de
con
cette
study
classe,
ypical
app
part
elée
free
A
computations.
ctivation
cus
R
an
e-
o
laxation
sampling
T
A
e
F
chnique
d.
(AR
con
T).
eld,
N
e
ous
Chapter
donnons
s
égalemen
applic
t
a
une
ystem
preuv
parallel
e
m
de
er
con
exp
v
the
ergence
against
p
to
o
c
u
v
r
ordinates.
un
ndings
algorithme
study
protot
con
yp
the
e.
as
Le
Chapter
deuxième
a
v
died
olet
e
p
n
orte
y
sur
faster
le
con
calcul
ABF
d'énergie
previously
l
wn.
ibre
f
p
m
our
simula
les
Abstract
tran-
sim
sitions
is
caractérisées
essen
par
to
une
i
co
understanding
ordonnée
c
de
and
réaction.
c
Nous
ical
nous
cesses
plaçons
real-life
dans
erime
le
ts
cadre
v
d'une
increasingly
métho
or
de
e
d'
infeasible.
éc
thesis
han
dev
tillonnage
to
d'im
dologi
p
asp
ortance
of
adaptativ
sim
e,
with
app
particular
elée
cus
A
computing
daptive
paths
Biasing
their
F
ciated
or
energy
c
The
e
part
(ABF).
de
Ce
to
v
metho
olet
for
comprend
ction
deux
and
sous-parties.
state
La
hes
première
a
partie
oten
(c
energy
hapitre
fac
5)
In
s'attac
3
he
e
à
ose
mon
impro
trer
emen
l
to
'applicabilité
widely-used
à
sition
un
searc
système
metho
biomoléculaire,
th
d
A
'
ation
une
T
nouv
hnique
elle
T).
mise
e
en
p
o
s
euvre
n
parall
a
èle
cal
d'ABF,
v
nommée
e
multiple-walker
of
ABF
protot
(MW-ABF),
algorithm.
consistan
second
t
is
à
to
utiliser
e
plusieurs
ergy
ré
W
pliques.
fo
Cette
in
mise
on
en
adaptiv
o
imp
euvre
r
s'est
e
a
tec
v
the
érée
daptiv
utile
Biasing
p
orce
our
metho
surmon
Th
ter
rst
des
tribution
problème
this
s
p
liés
s
à
n
un
in
mauv
5,
ais
in
c
ho
hoix
the
de
abilit
co
to
ordonnée
large
de
s
réaction.
of
Nous
new
conrmons
implemen
ensuite
named
ces
ultiple-w
résultats
k
n
ABF
umériques
Numerical
en
erimen
étu-
demonstrated
dian
robustness
t
MW-ABF
la
artefacts
c
due
on
p
v
orly
ergence
hosen
théorique
o
d'un
ersimplied
algorithme
co
d'ABF
These
adapté
umerical
dans
inspired
le
new
c
of
hapitre
longtime
6.
v
L'étude
of
de
ABF
c
d,
on
presen
v
in
ergence
6.
en
studying
temps
sligh
long
mo
utilise
mo
des
w
métho
bac
des
our
d'en
umerical
tropie
b
relativ
sho
e
a
et
theoretical
des
of
inégalités
v
de
of
Sob
than
olev
as
logarithmiques.
sho
Deterministic
and
Méthodes
tel-00597694, version 1 - 1 Jun 2011tel-00597694, version 1 - 1 Jun 2011on
10
our
Publica
Willaime.
tions
oc
d
9
ans
F.
des
metho
revues
J.
interna
d
tionales
.
1.
Chem.
T.
(2010).
Lelièvre,
K.
K.
of
Minouk
tr
adeh.
gy
L
130,
ong-time
1.
c
ho
onver
d'éner
genc
et
e
010).
of
Comput.
an

A
E.
daptive
M.-C.
Biasing
adeh,
F
impr
or
AR
c
for
e
p
metho
otential
d
es
:
Ph
the
(2009).
bi-chann
ng
e
Minouk
l
Une
c
adaptative
ase
c
.
libr
Sou
F
m
n
i
ach
s
J.
au
Theory
journal
6(4),
Arc
08
h.
1017
Rational
3.
Mec
Cancès,
h.
Legoll,
Anal.
Marinica,
(2010).
Minouk
Preprin
F.
t:
Some
hal-004
ovements
77302.
the
2.
T
K.
d
Minouk
nding
adeh,
ansition
C.
athways
Chip
p
ot,
ener
T.
surfac