Fundamentos de Mecánica Estadística y...
Transcript of Fundamentos de Mecánica Estadística y...
![Page 1: Fundamentos de Mecánica Estadística y Simulacionesm3g.iqm.unicamp.br/main/didatico/celfi/clase1.pdfBola 1 Bola 2 Bola 3 1 1 1 1 1 2 1 1 3 Nuevas reglas: A cada balde corresponde](https://reader033.fdocumento.com/reader033/viewer/2022053117/609ac52d2c2a07326d61fe0d/html5/thumbnails/1.jpg)
Fundamentos de Mecánica Estadística y Simulaciones
Leandro MartínezInstituto de Química - Universidad de Campinas (UNICAMP), Brasil
http://leandro.iqm.unicamp.br
![Page 2: Fundamentos de Mecánica Estadística y Simulacionesm3g.iqm.unicamp.br/main/didatico/celfi/clase1.pdfBola 1 Bola 2 Bola 3 1 1 1 1 1 2 1 1 3 Nuevas reglas: A cada balde corresponde](https://reader033.fdocumento.com/reader033/viewer/2022053117/609ac52d2c2a07326d61fe0d/html5/thumbnails/2.jpg)
Bola 1 Bola 2 Bola 3
1 1 1 1 1 2 1 1 3 1 2 1 1 2 2 1 2 3 1 3 1 1 3 2 1 3 3 2 1 1 2 1 2 2 1 3 2 2 1 2 2 2 2 2 3 2 3 1 2 3 2 2 3 3 3 1 1 3 1 2 3 1 3 3 2 1 3 2 2 3 2 3 3 3 1 3 3 2 3 3 3
3 x 3 x 3 = 27 resultados posibles
- Todos son igualmente probables.
![Page 3: Fundamentos de Mecánica Estadística y Simulacionesm3g.iqm.unicamp.br/main/didatico/celfi/clase1.pdfBola 1 Bola 2 Bola 3 1 1 1 1 1 2 1 1 3 Nuevas reglas: A cada balde corresponde](https://reader033.fdocumento.com/reader033/viewer/2022053117/609ac52d2c2a07326d61fe0d/html5/thumbnails/3.jpg)
Bola 1 Bola 2 Bola 3
1 1 1 1 1 2 1 1 3 1 2 1 1 2 2 1 2 3 1 3 1 1 3 2 1 3 3 2 1 1 2 1 2 2 1 3 2 2 1 2 2 2 2 2 3 2 3 1 2 3 2 2 3 3 3 1 1 3 1 2 3 1 3 3 2 1 3 2 2 3 2 3 3 3 1 3 3 2 3 3 3
Entre los resultados posibles:
El número “1” aparece: 27 veces
El número “2” aparece: 27 veces
El número “3” aparece: 27 veces
La probabilidad de encontrar una bola en particular en el balde 1 es de un tercio (~33%)
![Page 4: Fundamentos de Mecánica Estadística y Simulacionesm3g.iqm.unicamp.br/main/didatico/celfi/clase1.pdfBola 1 Bola 2 Bola 3 1 1 1 1 1 2 1 1 3 Nuevas reglas: A cada balde corresponde](https://reader033.fdocumento.com/reader033/viewer/2022053117/609ac52d2c2a07326d61fe0d/html5/thumbnails/4.jpg)
Nuevas reglas: A cada balde corresponde una “energia”.
“Energia del balde 1” = 1“Energia del balde 2” = 2“Energia del balde 3” = 3
Solo son admisibles sorteos en que la energía total es 5
Bola 1 Bola 2 Bola 3
1 1 1 1 1 2 1 1 3 1 2 1 1 2 2 1 2 3 1 3 1 1 3 2 1 3 3 2 1 1 2 1 2 2 1 3 2 2 1 2 2 2 2 2 3 2 3 1 2 3 2 2 3 3 3 1 1 3 1 2 3 1 3 3 2 1 3 2 2 3 2 3 3 3 1 3 3 2 3 3 3
Hay, ahora, 6 resultados admisibles.
Sorteo con restricciones
![Page 5: Fundamentos de Mecánica Estadística y Simulacionesm3g.iqm.unicamp.br/main/didatico/celfi/clase1.pdfBola 1 Bola 2 Bola 3 1 1 1 1 1 2 1 1 3 Nuevas reglas: A cada balde corresponde](https://reader033.fdocumento.com/reader033/viewer/2022053117/609ac52d2c2a07326d61fe0d/html5/thumbnails/5.jpg)
Entre los resultados admisibles:
El número “1” aparece: 9 veces
El número “2” aparece: 6 veces
El número “3” aparece: 3 veces
La probabilidad de que una bola esté en el balde 1 es 9/18 = 50%
Bola 1 Bola 2 Bola 3
1 1 1 1 1 2 1 1 3 1 2 1 1 2 2 1 2 3 1 3 1 1 3 2 1 3 3 2 1 1 2 1 2 2 1 3 2 2 1 2 2 2 2 2 3 2 3 1 2 3 2 2 3 3 3 1 1 3 1 2 3 1 3 3 2 1 3 2 2 3 2 3 3 3 1 3 3 2 3 3 3
Sorteo con restricciones
![Page 6: Fundamentos de Mecánica Estadística y Simulacionesm3g.iqm.unicamp.br/main/didatico/celfi/clase1.pdfBola 1 Bola 2 Bola 3 1 1 1 1 1 2 1 1 3 Nuevas reglas: A cada balde corresponde](https://reader033.fdocumento.com/reader033/viewer/2022053117/609ac52d2c2a07326d61fe0d/html5/thumbnails/6.jpg)
Conclusión:
Si se pone como restricción en el sorteo una condición sobre lasuma de una propiedad de posibles resultados, losresultados con mayores valores de esa propiedad se hacen relativamente menos probables.
Sin restricciones Con restricciones
33%33% 33%33% 33%33% 50%
33%17%
![Page 7: Fundamentos de Mecánica Estadística y Simulacionesm3g.iqm.unicamp.br/main/didatico/celfi/clase1.pdfBola 1 Bola 2 Bola 3 1 1 1 1 1 2 1 1 3 Nuevas reglas: A cada balde corresponde](https://reader033.fdocumento.com/reader033/viewer/2022053117/609ac52d2c2a07326d61fe0d/html5/thumbnails/7.jpg)
Gas (ideal) en un recipiente cerrado, con paredes conductoras de calor
Fluctuaciones en la Temperatura = Fluctuaciones en la Energia (cinética)
En cada instante del tiempo, las moléculas del gas tienen un conjunto distintode posiciones y velocidades: microestados
A cada microestado corresponde una energía
PREGUNTA: Con que probabilidad se encuentra cada microestado?
(en particular, esto respondería cual es la distribución de temperaturas)
![Page 8: Fundamentos de Mecánica Estadística y Simulacionesm3g.iqm.unicamp.br/main/didatico/celfi/clase1.pdfBola 1 Bola 2 Bola 3 1 1 1 1 1 2 1 1 3 Nuevas reglas: A cada balde corresponde](https://reader033.fdocumento.com/reader033/viewer/2022053117/609ac52d2c2a07326d61fe0d/html5/thumbnails/8.jpg)
Estrategia:
Construcción de un ENSEMBLE (conjunto) de sistemas, que forman como un todo un sistema aislado
La energia TOTAL es constante.
Los sistemas intercambian calor.
HIPOTESIS:
Los sistemas se distribuyen en losmicroestados posibles (posicionesy velocidades) ALEATORIAMENTE
...pero...
LA ENERGIA TOTAL ES CONSTANTE
![Page 9: Fundamentos de Mecánica Estadística y Simulacionesm3g.iqm.unicamp.br/main/didatico/celfi/clase1.pdfBola 1 Bola 2 Bola 3 1 1 1 1 1 2 1 1 3 Nuevas reglas: A cada balde corresponde](https://reader033.fdocumento.com/reader033/viewer/2022053117/609ac52d2c2a07326d61fe0d/html5/thumbnails/9.jpg)
Estrategia:
Construcción de un ENSEMBLE (conjunto) de sistemas, que forman como un todo un sistema aislado
Sorteo aleatório de sistemas en sus posibles microestados,donde cada microestado tiene una energia distinta, y hayuna restricción sobre la energia total:
RESULTADO: Microestados con energías menoresson mas probables que microestados con energíasmayores.
![Page 10: Fundamentos de Mecánica Estadística y Simulacionesm3g.iqm.unicamp.br/main/didatico/celfi/clase1.pdfBola 1 Bola 2 Bola 3 1 1 1 1 1 2 1 1 3 Nuevas reglas: A cada balde corresponde](https://reader033.fdocumento.com/reader033/viewer/2022053117/609ac52d2c2a07326d61fe0d/html5/thumbnails/10.jpg)
Estrategia:
Construcción de un ENSEMBLE (conjunto) de sistemas, que forman como un todo un sistema aislado
RESULTADO: Microestados con energías menoresson mas probables que microestados con energíasmayores.
Resultado analítico: (muchos sistemas,
muchos microestados)
Donde E es la energía del microestado.
![Page 11: Fundamentos de Mecánica Estadística y Simulacionesm3g.iqm.unicamp.br/main/didatico/celfi/clase1.pdfBola 1 Bola 2 Bola 3 1 1 1 1 1 2 1 1 3 Nuevas reglas: A cada balde corresponde](https://reader033.fdocumento.com/reader033/viewer/2022053117/609ac52d2c2a07326d61fe0d/html5/thumbnails/11.jpg)
Estrategia:
Construcción de un ENSEMBLE (conjunto) de sistemas, que forman como un todo un sistema aislado
Energía media del ensemble:
![Page 12: Fundamentos de Mecánica Estadística y Simulacionesm3g.iqm.unicamp.br/main/didatico/celfi/clase1.pdfBola 1 Bola 2 Bola 3 1 1 1 1 1 2 1 1 3 Nuevas reglas: A cada balde corresponde](https://reader033.fdocumento.com/reader033/viewer/2022053117/609ac52d2c2a07326d61fe0d/html5/thumbnails/12.jpg)
Energía media:
Cuenta alternativa:
![Page 13: Fundamentos de Mecánica Estadística y Simulacionesm3g.iqm.unicamp.br/main/didatico/celfi/clase1.pdfBola 1 Bola 2 Bola 3 1 1 1 1 1 2 1 1 3 Nuevas reglas: A cada balde corresponde](https://reader033.fdocumento.com/reader033/viewer/2022053117/609ac52d2c2a07326d61fe0d/html5/thumbnails/13.jpg)
Propriedades termodinámicas:
Postulado: La energía media del ensemble es la energía interna termodinámica.
... otras asociaciones adecuadas .. .
“Función de partición”Cantidad fundamental:
![Page 14: Fundamentos de Mecánica Estadística y Simulacionesm3g.iqm.unicamp.br/main/didatico/celfi/clase1.pdfBola 1 Bola 2 Bola 3 1 1 1 1 1 2 1 1 3 Nuevas reglas: A cada balde corresponde](https://reader033.fdocumento.com/reader033/viewer/2022053117/609ac52d2c2a07326d61fe0d/html5/thumbnails/14.jpg)
Termodinámica estadística vs. realidad
Las propriedades medias NO resultan de ningún sorteo aleatorio, pero de las propriedades mecánicas del sistema.
Hipótesis: Las propriedades calculadas como medias de los sorteos en el ensemble corresponden a las propriedades termodinámicas
vs.
ESTA HIPÓTESIS SE VERIFICA EXPERIMENTALMENTE
![Page 15: Fundamentos de Mecánica Estadística y Simulacionesm3g.iqm.unicamp.br/main/didatico/celfi/clase1.pdfBola 1 Bola 2 Bola 3 1 1 1 1 1 2 1 1 3 Nuevas reglas: A cada balde corresponde](https://reader033.fdocumento.com/reader033/viewer/2022053117/609ac52d2c2a07326d61fe0d/html5/thumbnails/15.jpg)
Termodinámica estadística vs. realidad
Hipótesis ergódica
Tiempo?
Hay que hacer alguna suposición sobre la manera como el sistema evoluciona en el tiempo...
Dada una propriedad A, el sistema pasa un tiempo en microestadoscon valor A = A' proporcional al número de microestados donde A = A'
o
El sistema pasa por todos los microestados antes de volver al microestado inicial
(implica que todos los microestados son igualmente probables)
![Page 16: Fundamentos de Mecánica Estadística y Simulacionesm3g.iqm.unicamp.br/main/didatico/celfi/clase1.pdfBola 1 Bola 2 Bola 3 1 1 1 1 1 2 1 1 3 Nuevas reglas: A cada balde corresponde](https://reader033.fdocumento.com/reader033/viewer/2022053117/609ac52d2c2a07326d61fe0d/html5/thumbnails/16.jpg)
Termodinámica estadística vs. realidad
Hipótesis ergódica
Tiempo?
Las medias de propriedades calculadas sobre el ensemblecorresponden a las medias calculadas sobre el tiempo si la trayectoriaes suficientemente larga.
- Observación (simulación) en el tiempo permite el estudio de propriedades usando la mecánica estadística
- Esta hipótesis puede ser puesta a prueba haciendo solamente simulaciones (ensembles vs. dinámica).
![Page 17: Fundamentos de Mecánica Estadística y Simulacionesm3g.iqm.unicamp.br/main/didatico/celfi/clase1.pdfBola 1 Bola 2 Bola 3 1 1 1 1 1 2 1 1 3 Nuevas reglas: A cada balde corresponde](https://reader033.fdocumento.com/reader033/viewer/2022053117/609ac52d2c2a07326d61fe0d/html5/thumbnails/17.jpg)
Resumen
Pregunta: Como calcular propriedades de sistemas macroscópicos a partir de sus propriedades microscópicas?
Estrategia: Construir conjuntos (ensembles) de sistemas equivalentes, y suponer que los sistemas se distribuyen aleatoriamente entre sus posibles estados microscópicos, satisfaciendo condiciones necesarias del conjunto.
Hipótesis: Las propriedades calculadas a partir del resultado de este sorteo corresponden a propriedades termodinámicas.
Hipótesis: Las medias calculadas a partir del sorteo corresponden a medias temporales: dinámica
![Page 18: Fundamentos de Mecánica Estadística y Simulacionesm3g.iqm.unicamp.br/main/didatico/celfi/clase1.pdfBola 1 Bola 2 Bola 3 1 1 1 1 1 2 1 1 3 Nuevas reglas: A cada balde corresponde](https://reader033.fdocumento.com/reader033/viewer/2022053117/609ac52d2c2a07326d61fe0d/html5/thumbnails/18.jpg)
Tipos de ENSEMBLESLos sistemas deben satisfacer las propriedades termodinámicas de interés
Ensemble NVT : “Canónico”
Ensemble μVT: “Grande canónico”
Ensemble NpT“Isotérmico-isobárico”
Ensemble NVE“Microcanónico”
Energía constante.
T variablep variable
![Page 19: Fundamentos de Mecánica Estadística y Simulacionesm3g.iqm.unicamp.br/main/didatico/celfi/clase1.pdfBola 1 Bola 2 Bola 3 1 1 1 1 1 2 1 1 3 Nuevas reglas: A cada balde corresponde](https://reader033.fdocumento.com/reader033/viewer/2022053117/609ac52d2c2a07326d61fe0d/html5/thumbnails/19.jpg)
Simulaciones
![Page 20: Fundamentos de Mecánica Estadística y Simulacionesm3g.iqm.unicamp.br/main/didatico/celfi/clase1.pdfBola 1 Bola 2 Bola 3 1 1 1 1 1 2 1 1 3 Nuevas reglas: A cada balde corresponde](https://reader033.fdocumento.com/reader033/viewer/2022053117/609ac52d2c2a07326d61fe0d/html5/thumbnails/20.jpg)
SIMULACIONES
Medias sobre conjuntos Medias temporales
Propriedades termodinámicas
ergodicidad
Muestreo aleatorio de microestados(configuraciones)
Muestreo de configuracionesa lo largo del tiempo
![Page 21: Fundamentos de Mecánica Estadística y Simulacionesm3g.iqm.unicamp.br/main/didatico/celfi/clase1.pdfBola 1 Bola 2 Bola 3 1 1 1 1 1 2 1 1 3 Nuevas reglas: A cada balde corresponde](https://reader033.fdocumento.com/reader033/viewer/2022053117/609ac52d2c2a07326d61fe0d/html5/thumbnails/21.jpg)
SIMULACIONES
Medias sobre conjuntos Medias temporales
Propriedades termodinámicas
ergodicidad
Muestreo aleatorio de microestados(configuraciones)
Muestreo de configuracionesa lo largo del tiempo
Simulaciones estocásticas Dinámica ... Molecular
![Page 22: Fundamentos de Mecánica Estadística y Simulacionesm3g.iqm.unicamp.br/main/didatico/celfi/clase1.pdfBola 1 Bola 2 Bola 3 1 1 1 1 1 2 1 1 3 Nuevas reglas: A cada balde corresponde](https://reader033.fdocumento.com/reader033/viewer/2022053117/609ac52d2c2a07326d61fe0d/html5/thumbnails/22.jpg)
SIMULACIONES DE MONTE-CARLO (ESTOCÁSTICAS)
Objetivo: Calcular el valor de una propriedad cualquiera.
(por ejemplo: la energía interna)
Ejemplo: Líquido contenido en recipiente con volumen y temperatura constantes.
![Page 23: Fundamentos de Mecánica Estadística y Simulacionesm3g.iqm.unicamp.br/main/didatico/celfi/clase1.pdfBola 1 Bola 2 Bola 3 1 1 1 1 1 2 1 1 3 Nuevas reglas: A cada balde corresponde](https://reader033.fdocumento.com/reader033/viewer/2022053117/609ac52d2c2a07326d61fe0d/html5/thumbnails/23.jpg)
SIMULACIONES DE MONTE-CARLO (ESTOCÁSTICAS)
E = V + K
E
E
E
E
Como hacer un sorteo de sobrelos microestados que respete laprobabilidad relativa de cada región?
Para T definido, K=constante:
![Page 24: Fundamentos de Mecánica Estadística y Simulacionesm3g.iqm.unicamp.br/main/didatico/celfi/clase1.pdfBola 1 Bola 2 Bola 3 1 1 1 1 1 2 1 1 3 Nuevas reglas: A cada balde corresponde](https://reader033.fdocumento.com/reader033/viewer/2022053117/609ac52d2c2a07326d61fe0d/html5/thumbnails/24.jpg)
SIMULACIONES DE MONTE-CARLO (ESTOCÁSTICAS)
Como hacer un sorteo de sobrelos microestados que respete laprobabilidad relativa de cada región?
![Page 25: Fundamentos de Mecánica Estadística y Simulacionesm3g.iqm.unicamp.br/main/didatico/celfi/clase1.pdfBola 1 Bola 2 Bola 3 1 1 1 1 1 2 1 1 3 Nuevas reglas: A cada balde corresponde](https://reader033.fdocumento.com/reader033/viewer/2022053117/609ac52d2c2a07326d61fe0d/html5/thumbnails/25.jpg)
SIMULACIONES DE MONTE-CARLO
Como hacer un sorteo de sobrelos microestados que respete laprobabilidad relativa de cada región?
Este proceso genera poblaciones de i y j con las proporciones correctas
- Empiezo en un dado i
- Genero un nuevo punto j
- Si Pj > Pi acepto j (con probabilidad 1)
- Si Pj < Pi acepto j con probabilidad Pj/Pi
Definir probabilidades de transición consistentes
Como la tasa de aceptación es 1 en una dirección, y en la otra, la tasa relativa es correcta.
![Page 26: Fundamentos de Mecánica Estadística y Simulacionesm3g.iqm.unicamp.br/main/didatico/celfi/clase1.pdfBola 1 Bola 2 Bola 3 1 1 1 1 1 2 1 1 3 Nuevas reglas: A cada balde corresponde](https://reader033.fdocumento.com/reader033/viewer/2022053117/609ac52d2c2a07326d61fe0d/html5/thumbnails/26.jpg)
- Empiezo en un dado i
- Genero un nuevo punto j
- Si Vj > Vi acepto j con probabilidad Pj/Pi
- Si Vj < Vi acepto j (con probabilidad 1)
SIMULACIONES DE MONTE-CARLO
Si Pj < Pi oPj/Pi > rand[0,1]
Si no, vuelvo a la estructura anterior
= Crear una nueva estructura(generalmente perturbando la anterior)
Si Pj < Pi oPj/Pi > rand[0,1]
![Page 27: Fundamentos de Mecánica Estadística y Simulacionesm3g.iqm.unicamp.br/main/didatico/celfi/clase1.pdfBola 1 Bola 2 Bola 3 1 1 1 1 1 2 1 1 3 Nuevas reglas: A cada balde corresponde](https://reader033.fdocumento.com/reader033/viewer/2022053117/609ac52d2c2a07326d61fe0d/html5/thumbnails/27.jpg)
SIMULACIONES DE MONTE-CARLO
Conjunto de estructuras que son un muestreo correctamente ponderado del conjunto total
- Cálculo de propiedades que no son explícitamente dependientes del tiempo.
- Si la hipótesis ergódica es correcta, las propriedades calculadas corresponden a las medias temporales sobre la evolución dinámica del sistema
![Page 28: Fundamentos de Mecánica Estadística y Simulacionesm3g.iqm.unicamp.br/main/didatico/celfi/clase1.pdfBola 1 Bola 2 Bola 3 1 1 1 1 1 2 1 1 3 Nuevas reglas: A cada balde corresponde](https://reader033.fdocumento.com/reader033/viewer/2022053117/609ac52d2c2a07326d61fe0d/html5/thumbnails/28.jpg)
SIMULACIONES
Medias sobre conjuntos Medias temporales
Propriedades termodinámicas
ergodicidad
Muestreo aleatorio de microestados(configuraciones)
Muestreo de configuracionesa lo largo del tiempo
Simulaciones estocásticas Dinámica ... Molecular
![Page 29: Fundamentos de Mecánica Estadística y Simulacionesm3g.iqm.unicamp.br/main/didatico/celfi/clase1.pdfBola 1 Bola 2 Bola 3 1 1 1 1 1 2 1 1 3 Nuevas reglas: A cada balde corresponde](https://reader033.fdocumento.com/reader033/viewer/2022053117/609ac52d2c2a07326d61fe0d/html5/thumbnails/29.jpg)
1 - Cómo los objetos interactúan:
2 -Posiciones y velocidades iniciales: 30 de Junio de 2009, 18h 37m 54s ...
3 - Integración de las ecuaciones de movimiento:
SIMULACIONES DE DINÁMICA MOLECULAR
![Page 30: Fundamentos de Mecánica Estadística y Simulacionesm3g.iqm.unicamp.br/main/didatico/celfi/clase1.pdfBola 1 Bola 2 Bola 3 1 1 1 1 1 2 1 1 3 Nuevas reglas: A cada balde corresponde](https://reader033.fdocumento.com/reader033/viewer/2022053117/609ac52d2c2a07326d61fe0d/html5/thumbnails/30.jpg)
SIMULACIONES DE DINÁMICA MOLECULAR
Interacciones intra e inter-moleculares
Propiedades medias
Condiciones termodinámicas
Escalas de tiempo (pico, nano, mili-segundos)
: del orden de fs (10-12 s) - 1 o 2 fs para simulacionesatomísticas
En el intervalo los átomos se mueven con
aceleración constante: aproximación
![Page 31: Fundamentos de Mecánica Estadística y Simulacionesm3g.iqm.unicamp.br/main/didatico/celfi/clase1.pdfBola 1 Bola 2 Bola 3 1 1 1 1 1 2 1 1 3 Nuevas reglas: A cada balde corresponde](https://reader033.fdocumento.com/reader033/viewer/2022053117/609ac52d2c2a07326d61fe0d/html5/thumbnails/31.jpg)
SIMULACIONES DE DINÁMICA MOLECULAR
Sin aproximaciones, las ecuaciones de Newton conservan la energía:
Segunda lei de Newton:
cqd
![Page 32: Fundamentos de Mecánica Estadística y Simulacionesm3g.iqm.unicamp.br/main/didatico/celfi/clase1.pdfBola 1 Bola 2 Bola 3 1 1 1 1 1 2 1 1 3 Nuevas reglas: A cada balde corresponde](https://reader033.fdocumento.com/reader033/viewer/2022053117/609ac52d2c2a07326d61fe0d/html5/thumbnails/32.jpg)
SIMULACIONES DE DINÁMICA MOLECULAR
En general, NO QUEREMOS energia constante...
Queremos Temperatura constante
Ex. Monte-Carlo:
Como mantener la TEMPERATURA constante, si la integración de las equacionesconserva la ENERGIA TOTAL?
TERMOSTATOS
Isocinético
Berendsen
Noose-Hover
Langevin
Rescalonamiento de las velocidadesNo generan ensembles correctos
Modificación de las ecuaciones de movimientoGeneran ensembles correctos
![Page 33: Fundamentos de Mecánica Estadística y Simulacionesm3g.iqm.unicamp.br/main/didatico/celfi/clase1.pdfBola 1 Bola 2 Bola 3 1 1 1 1 1 2 1 1 3 Nuevas reglas: A cada balde corresponde](https://reader033.fdocumento.com/reader033/viewer/2022053117/609ac52d2c2a07326d61fe0d/html5/thumbnails/33.jpg)
SIMULACIONES DE DINÁMICA MOLECULAR
TERMOSTATOS
Isocinético:
A cada tanto, reescalonar las velocidades de modoque correspondan a la temperatura deseada.
O sea, multiplicar todas las velocidades por hace con que la energía cinética mediacorresponda a la temperatura deseada.
Usando:
![Page 34: Fundamentos de Mecánica Estadística y Simulacionesm3g.iqm.unicamp.br/main/didatico/celfi/clase1.pdfBola 1 Bola 2 Bola 3 1 1 1 1 1 2 1 1 3 Nuevas reglas: A cada balde corresponde](https://reader033.fdocumento.com/reader033/viewer/2022053117/609ac52d2c2a07326d61fe0d/html5/thumbnails/34.jpg)
SIMULACIONES DE DINÁMICA MOLECULAR
TERMOSTATOS
Isocinético:
A cada tanto, reescalonar las velocidades de modoque correspondan a la temperatura deseada.
![Page 35: Fundamentos de Mecánica Estadística y Simulacionesm3g.iqm.unicamp.br/main/didatico/celfi/clase1.pdfBola 1 Bola 2 Bola 3 1 1 1 1 1 2 1 1 3 Nuevas reglas: A cada balde corresponde](https://reader033.fdocumento.com/reader033/viewer/2022053117/609ac52d2c2a07326d61fe0d/html5/thumbnails/35.jpg)
SIMULACIONES DE DINÁMICA MOLECULAR
TERMOSTATOS
Berendsen:
Instantáneamente lleva de a
Lleva de a lentamente, dependiendo del valor de
: paso de tiempo de la integración numérica: paso de tiempo de la integración numérica
: constante de acoplamiento del termostato
![Page 36: Fundamentos de Mecánica Estadística y Simulacionesm3g.iqm.unicamp.br/main/didatico/celfi/clase1.pdfBola 1 Bola 2 Bola 3 1 1 1 1 1 2 1 1 3 Nuevas reglas: A cada balde corresponde](https://reader033.fdocumento.com/reader033/viewer/2022053117/609ac52d2c2a07326d61fe0d/html5/thumbnails/36.jpg)
SIMULACIONES DE DINÁMICA MOLECULAR
TERMOSTATOS
Berendsen:
Dinámicas que no generan formalmente un ensemble canónico
Importante para sistemas con pocas partículas: Grandes fluctuaciones de temperatura
Menos importante para sistemas con muchas partículas: Pequeñas fluctuaciones de temperatura
Mejor NO usar estos termostatos en simulaciones productivas.
![Page 37: Fundamentos de Mecánica Estadística y Simulacionesm3g.iqm.unicamp.br/main/didatico/celfi/clase1.pdfBola 1 Bola 2 Bola 3 1 1 1 1 1 2 1 1 3 Nuevas reglas: A cada balde corresponde](https://reader033.fdocumento.com/reader033/viewer/2022053117/609ac52d2c2a07326d61fe0d/html5/thumbnails/37.jpg)
SIMULACIONES DE DINÁMICA MOLECULAR
TERMOSTATOS QUE MUESTREAN BIEN
partícula real
partícula virtual
![Page 38: Fundamentos de Mecánica Estadística y Simulacionesm3g.iqm.unicamp.br/main/didatico/celfi/clase1.pdfBola 1 Bola 2 Bola 3 1 1 1 1 1 2 1 1 3 Nuevas reglas: A cada balde corresponde](https://reader033.fdocumento.com/reader033/viewer/2022053117/609ac52d2c2a07326d61fe0d/html5/thumbnails/38.jpg)
SIMULACIONES DE DINÁMICA MOLECULAR
TERMOSTATOS QUE MUESTREAN BIEN
Termostato (o dinámica de) Langevin:
La partícula real está imersa enun baño de partículas virtuales, muchomenores
- Movimiento Browniano
Componentes: 1. Las partículas son desaceleradas por la fricción 2. Las partículas sufren choques estocásticos
Fricciónproporcional a la velocidad
Choques estocásticos de un bañocon temperatura
Variable aleatória (gaussiana)con média 0 y desvio standard 1
= coeficiente de fricción
![Page 39: Fundamentos de Mecánica Estadística y Simulacionesm3g.iqm.unicamp.br/main/didatico/celfi/clase1.pdfBola 1 Bola 2 Bola 3 1 1 1 1 1 2 1 1 3 Nuevas reglas: A cada balde corresponde](https://reader033.fdocumento.com/reader033/viewer/2022053117/609ac52d2c2a07326d61fe0d/html5/thumbnails/39.jpg)
SIMULACIONES DE DINÁMICA MOLECULAR
TERMOSTATOS QUE MUESTREAN BIEN
Termostato (o dinámica de) Langevin:
La partícula real está imersa enun baño de partículas virtuales, muchomenores
- Movimiento Browniano
tiene constante.
Implementación:
Cambios instantáneos de velocidad
Fricción
![Page 40: Fundamentos de Mecánica Estadística y Simulacionesm3g.iqm.unicamp.br/main/didatico/celfi/clase1.pdfBola 1 Bola 2 Bola 3 1 1 1 1 1 2 1 1 3 Nuevas reglas: A cada balde corresponde](https://reader033.fdocumento.com/reader033/viewer/2022053117/609ac52d2c2a07326d61fe0d/html5/thumbnails/40.jpg)
Fin Clase 1
![Page 41: Fundamentos de Mecánica Estadística y Simulacionesm3g.iqm.unicamp.br/main/didatico/celfi/clase1.pdfBola 1 Bola 2 Bola 3 1 1 1 1 1 2 1 1 3 Nuevas reglas: A cada balde corresponde](https://reader033.fdocumento.com/reader033/viewer/2022053117/609ac52d2c2a07326d61fe0d/html5/thumbnails/41.jpg)
SIMULACIONES DE DINÁMICA MOLECULAR
TERMOSTATOS QUE MUESTREAN BIEN
Nosé-Hover:
Movimiento:
Se propagan lasposiciones para losátomos reales y para s
Interación conátomos “virtuales”
![Page 42: Fundamentos de Mecánica Estadística y Simulacionesm3g.iqm.unicamp.br/main/didatico/celfi/clase1.pdfBola 1 Bola 2 Bola 3 1 1 1 1 1 2 1 1 3 Nuevas reglas: A cada balde corresponde](https://reader033.fdocumento.com/reader033/viewer/2022053117/609ac52d2c2a07326d61fe0d/html5/thumbnails/42.jpg)
SIMULACIONES DE DINÁMICA MOLECULAR
TERMOSTATOS QUE MUESTREAN BIEN
Nosé-Hover:
Se propagan las posiciones para los átomos reales y para s
La energía total U se conserva para el sistema completo.
Los átomos reales están en contacto térmico con elsistema virtual. (Hay transferéncia de energia entre losgrados de libertad de los átomos reales y los de s).
La energía de los átomos reales oscila alrededor de lala energía interna termodinámica a la temperatura deseada.
EL SISTEMA COMPLETO ES MICROCANÓNICO, PERO EL SISTEMA REAL ES CANÓNICO(energía constante) (temperatura)
![Page 43: Fundamentos de Mecánica Estadística y Simulacionesm3g.iqm.unicamp.br/main/didatico/celfi/clase1.pdfBola 1 Bola 2 Bola 3 1 1 1 1 1 2 1 1 3 Nuevas reglas: A cada balde corresponde](https://reader033.fdocumento.com/reader033/viewer/2022053117/609ac52d2c2a07326d61fe0d/html5/thumbnails/43.jpg)
SIMULACIONES DE DINÁMICA MOLECULAR
TERMOSTATOS QUE MUESTREAN BIEN*
*Las estructuras, como una simulación de MC
Pero la dinámica delsistema puede no sermás realista.
Propiedades dinámicasgeneralmente se obtienen de simulacionesNVE.
![Page 44: Fundamentos de Mecánica Estadística y Simulacionesm3g.iqm.unicamp.br/main/didatico/celfi/clase1.pdfBola 1 Bola 2 Bola 3 1 1 1 1 1 2 1 1 3 Nuevas reglas: A cada balde corresponde](https://reader033.fdocumento.com/reader033/viewer/2022053117/609ac52d2c2a07326d61fe0d/html5/thumbnails/44.jpg)
SIMULACIONES DE DINÁMICA MOLECULAR
TERMOSTATOS
Isocinético:
A cada tanto, reescalonar las velocidades de modoque correspondan a la temperatura deseada.
No sabemos, al principio,cual es la energía totala la cual corresponderáuna energía cinéticamedia adecuada.
Energía total =energía interna termodinámica