Laboratorio di Simulazione Atomistica e Fluidodinamica

Laboratorio di Simulazione
Atomistica e Fluidodinamica
dinamica molecolare classica
Mauro Chinappi
Center for Life Nano Science@Sapienza
Istituto Italiano di Tecnologia
… per lavorare meglio
●
Scaricare i software NAMD e VMD sui
vostri computer
Provare VMD e NAMD seguendo i tutorial presenti sul
sito (con VMD è facile, con NAMD potrebbe essere un
po' complesso). Basta anche solo l'inizio di un tutorial,
serve solo sapere se funzionano.
●
Mauro Chinappi
Dinamica molecolare classica (in pillole)
Ogni singolo atomo è modellato come un punto
materiale di massa m e carica q
Su ogni atomo agisce una forza (conservativa)
dipendente dalle posizioni degli altri atomi
Posizioni e velocità degli atomi sono ottenute
integrando le equazioni di Hamilton
Inoltre
Possibilità accoppiamento sistema con bagno termico
o pistone
●
Condizioni periodiche al bordo
●
Possibilità di aggiungere forse esterne
(anche non conservative)
●
Mauro Chinappi
Struttura generale di un codice MD
Configurazione
iniziale
Modello
(potenziale di
interazione)
Istruzioni
varie
Mauro Chinappi
Osservabili
A(t) = A(p(t),q(t))
Integrazione
equazioni moto
Configurazioni
(p(t), q(t))
Altre info
Integrazione: termostati e barostati
Eq. Hamilton integrate
con algoritmo di Verlet
NVE
(Hamilton)
NVT
(termostato)
NPT
(barostato)
Mauro Chinappi
Modifica
evoluzione
velocità
2
mv
T =∑i =1
3 N kb
N
Modifica
N
N
1
1
evoluzione
P= NkT −
r f
∑
∑
posizioni
ij ij
V
3 kT
i=1 j=i1
(e velocità)
[
]
Struttura generale di un codice MD
Configurazione
iniziale
Modello
(potenziale di
interazione)
Istruzioni
varie
Mauro Chinappi
Osservabili
A(t) = A(p(t),q(t))
Integrazione
equazioni moto
Configurazioni
(p(t), q(t))
Altre info
Potenziali di interazione
Forze incluse nei modelli di interazione per MD classica
●
Dispersione (dipolo-indotto dipolo-indotto)
●
Volume escluso
●
Elettrostatica
●
Legami covalenti
Mauro Chinappi
}
Lennard-Jones
Periodicità, raggio di cut-off e minima immagine
Mauro Chinappi
Potenziali di interazione
Forze incluse nei modelli di interazione atomistica
per MD classica
●
Dispersione (dipolo-indotto dipolo-indotto)
●
Volume escluso
●
Elettrostatica
●
Legami covalenti
Mauro Chinappi
Problema
periodicità
(somme di Ewald)
Potenziali di interazione
Forze incluse nei modelli di interazione per MD classica
●
Dispersione (dipolo-indotto dipolo-indotto)
●
Volume escluso
●
Elettrostatica
●
Legami covalenti
Vincoli o potenziali
armonici?
Mauro Chinappi
Potenziali di interazione
LJ + elettrostatica
Mauro Chinappi
Qualche domanda (… a cui dovreste saper rispondere)
Quanti atomi ha (più o meno) una piccola proteina
globulare ?
●
I famosi legami idrofobi, dove sono nella MD?
●
Quante molecole di acqua ci sono in un cubo di 4
nanometri di lato?
●
Fornire degli argomenti (convincenti) del fatto nei liquidi
esistono forze attrattive tra le molecole
●
Mauro Chinappi
Suggerimenti
Scirvere un paio di programmini FORTRAN
(programma_1, scrivo dei dati in un file, programma_2 li leggo
e calcolo la media etc.)
●
Usare gnuplot per visualizzare i dati generati dal programma_1
●
Tutorial di VMD (è anche abbastanza divertente)
●
Ripasso unità di misura (Angstrom, picosecondo) ed
“equivalenze” (Angstrom/ps = ??? m/s)
●
Mauro Chinappi
Lab. di simulazione atomistica e fluidodinamica
?
Mauro Chinappi
NAMD e VMD
NAnoscale Molecular Dynamics
●
University of Illinois at Urbana-Champaign
●
http://www.ks.uiuc.edu/Research/namd/
●
Altri codici per MD classica
●
GROMACS
LAMMPS
●DL_POLY (DL_PROTEINS e simili)
●AMBER
●
Mauro Chinappi
Struttura generale di un codice MD
Configuarazione
iniziale
Modello
(potenziale di
interazione)
Istruzioni
varie
Mauro Chinappi
Osservabili
A(p,q)
Integrazione
Configurazioni
(p(t), q(t))
Altre info
NAMD I/O (base)
Standard
output
File .pdb
File .psf
(e par)
File .conf
Mauro Chinappi
Integrazione
(NVE,
Vari NVT,
Vari NPT,
etc )
.dcd
Altre info
.pdb (configurazione iniziale)
http://www.rcsb.org/
ATOM 1 N LYS A 1 ­14.189 28.577 49.986 1.00 59.37 N
ATOM 2 CA LYS A 1 ­14.427 27.969 48.643 1.00 60.26 C
ATOM 3 C LYS A 1 ­14.388 26.456 48.756 1.00 60.10 C
ATOM 4 O LYS A 1 ­14.139 25.752 47.779 1.00 60.59 O
ATOM 5 CB LYS A 1 ­13.396 28.459 47.636 1.00 60.34 C
Mauro Chinappi
.psf (configurazione iniziale)
6682 !NATOM
1 U 1 MET N NH3 ­0.300000 14.0070 0
2 U 1 MET HT1 HC 0.330000 1.0080 0
3 U 1 MET HT2 HC 0.330000 1.0080 0
4 U 1 MET HT3 HC 0.330000 1.0080 0
5 U 1 MET CA CT1 0.210000 12.0110 0
6 U 1 MET HA HB 0.100000 1.0080 0
...
NOTA
●
●
Per generare un psf a partire da un pdb è necessario
usare uno strumento esterno (psfgen) che è incluso in
VMD.
In generale la creazione dei file con la configurazione
iniziale e del campo di forza del sistema è un aspetto
non banale.
Mauro Chinappi
.conf
Informazioni minime
Percorso file pdb e psf (configurazione iniziale e potenziale
di interazione)
●
Nomi file output
●
Raggio di cutoff e altri parametri per Lennard Jones
●
Parametri per integrazione elettrostatica
●
Altre info
Termostati/ barostati
●
Minimizzazione
●
Etc etc ...
●
Mauro Chinappi
Suggerimenti
Scaricare i software NAMD e VMD sui
vostri computer
●
Provare VMD e NAMD seguendo i tutorial
presenti sul sito (con VMD è facile, con
NAMD potrebbe essere un po'
complesso). Basta anche solo l'inizio di un
tutorial, serve solo sapere se funzionano.
●
http://www.ks.uiuc.edu/Training/Tutorials/
●
Mauro Chinappi
Lab. di simulazione atomistica e fluidodinamica
?
Mauro Chinappi