Dinámica molecular


Dinámica molecular

La Dinámica Molecular (DM) es una técnica de simulación en la que se permite que átomos y moléculas interactúen por un período, permitiendo una visualización del movimiento de las partículas. Originalmente fue concebida dentro de la física teórica, aunque hoy en día se utiliza sobre todo en biofísica y ciencia de materiales. Su campo de aplicación va desde superficies catalíticas hasta sistemas biológicos como las proteínas. Si bien los experimentos de cristalografía de rayos X permiten tomar "fotografías estáticas" y la técnica de RMN nos da indicios del movimiento molecular, ningún experimento es capaz de acceder a todas las escalas de tiempo involucradas. Resulta tentador, aunque no es enteramente correcto, describir a la DM como un "microscopio virtual" con alta resolución espacial y temporal.

En general, los sistemas moleculares son complejos y consisten de un gran número de partículas, por lo cual sería imposible encontrar sus propiedades de forma analítica. Para evitar este problema, la DM utiliza métodos numéricos. La DM representa un punto intermedio entre los experimentos y la teoría. Puede ser entendida como un experimento en la computadora.

Sabemos que la materia está constituida de partículas en movimiento e interacción al menos desde la época de Boltzmann en el siglo XIX. Pero muchos aún se imaginan a las moléculas como los modelos estáticos de un museo. Richard Feynman dijo en 1963 que "todo lo que hacen los seres vivos puede ser entendido a través de los saltos y contorsiones de los átomos."[1] Una de las contribuciones más importantes de la dinámica molecular es crear conciencia de que el DNA y las proteínas son máquinas en movimiento. Se le utiliza para explorar la relación entre estructura, movimiento y función.

La dinámica molecular es un campo multidisciplinario. Sus leyes y teorías provienen de las Matemáticas, Física y Química. Emplea algoritmos de las Ciencias de la Computación y Teoría de la información. Permite entender a los materiales y las moléculas no cómo entidades rígidas, sino como cuerpos animados. También se le ha llamado "estadística mecánica numérica" o "la visión de Laplace de la mecánica Newtoniana", en el sentido de predecir el futuro al animar las fuerzas de la naturaleza.

Para utilizar esta técnica de forma correcta, es importante entender las aproximaciones utilizadas y evitar caer en el error conceptual de que estamos simulando el comportamiento real y exacto de un sistema molecular. La integración de las ecuaciones de movimiento están mal condicionadas, lo cual genera errores numéricos acumulativos, que pueden ser minimizados seleccionando apropiadamente los algoritmos, pero no eliminados del todo. Por otro lado, las interacciones entre las partículas se modelan con un campo de fuerza aproximado, que puede o no ser adecuado dependiendo del problema que queremos resolver. De cualquier forma, la dinámica molecular nos permite explorar su comportamiento representativo en el espacio fásico.

Ejemplo de de una simulación de un sistema simple por el método de dinámica molecular: deposición de un Átomo de Cu en una superficie de Cu (001). Cada círculo ilustra la posición de un átomo; note que las verdaderas interacciones atómicas usadas en simulación son más complejas que las bidimensionales mostradas en la figura.

En la DM, hay que balancear el costo computacional y la fiabilidad en los resultados. En la DM clásica se utilizan las Ecuaciones de Newton, cuyo costo computacional es mucho menor que el de las de la mecánica cuántica. Es por ello que muchas propiedades que pueden resultar de interés, como la formación o ruptura de enlaces no puedan ser estudiadas mediante este método ya que no contempla estados excitados o reactividad.

Existen métodos híbridos denominados QM/MM (Quantum Mechanics/Molecular Mechanics) en los que un centro reactivo es tratado de modo cuántico mientras que el ambiente que lo rodea se trata de modo clásico. El desafío en este tipo de métodos resulta en la definición de manera precisa de la interacción entre los dos formas de describir el sistema...

El resultado de una simulación de dinámica molecular son las posiciones X y velocidades V de cada átomo de la molécula, para cada instante en el tiempo discretizado. A esto se le llama trayectoria.

Contenido

Principios físicos

Colectividad microcanónica (NVE)

La forma más simple de dinámica molecular ocurre en el colectividad microcanónica. En él, el sistema está aislado: su volumen no se altera (V) y no intercambia masa (N) ni energía (E) con el entorno. Para un sistema de N partículas con coordenadas X y velocidades V, se puede plantear el siguiente par de ecuaciones diferenciales de primer orden

F(X) = -\nabla U(X)=M\dot{V}(t)
\dot{X}(t)=V(t)

La función de energía potencial U(X) son las atracciones y repulsiones que sienten los átomos entre sí debido a los enlaces químicos, interacciones electrostáticas, van der Waals, etc. de las moléculas. A U(X) también se le conoce como campo potencial de fuerza y es una función de las coordenadas de las partículas X. Normalmente proviene de cálculos de química cuántica y/o experimentos espectroscópicos. Sin embargo, el campo de fuerza generalmente tiene una forma funcional que lo hace pertenecer a la mecánica clásica.

La trayectoria de las partículas es discreta en el tiempo. Normalmente se elige un paso de tiempo suficientemente pequeño (p.ej. 1 femtosegundo) para evitar errores numéricos de discretización. Para cada paso de tiempo, se integra la posición X y velocidad V con un método simpléctico como la integración de Verlet. Dadas las posiciones iniciales (p.ej. la estructura de rayos X de una proteína) y las velocidades iniciales (p.ej. aleatorias y Gaussianas), es posible calcular todas las posiciones y velocidades en el futuro.

La colectividad NVE es difícil de realizar experimentalmente. Por ello, generalmente se utiliza la colectividad NVT o NPT en las simulaciones. Sin embargo, las simulaciones NVE son importantes para verificar que un campo de fuerza combinado con un algoritmo de integración conserva la energía del sistema.

Ensamble canónico (NVT)

En el ensamble canónico, el volumen no se altera (V) y no se intercambia masa (N). La temperatura (T) se mantiene alrededor de la media deseada. La temperatura instantánea del sistema no es constante, sino sólo su promedio. En la colectividad NVT, la energía de los procesos endotérmicos y exotérmicos (transiciones conformacionales entre estados con diferentes energías potenciales) es intercambiada con un termostato.

Existe una gran variedad de termostatos para añadir o remover energía para mantener la temperatura promedio constante. Incluyen la reescalación de velocidades, dinámica Langevin, termostato Nosé-Hoover y el termostato Berendsen. No es fácil obtener una distribución de microestados correspondiente a la colectividad canónica, ya que ello depende del tamaño del sistema, selección del termostato, parámetros del termostato, paso de tiempo e integrador utilizados.

Software

En la actualidad existen varios programas capaces de llevar a cabo las simulaciones de DM, y a su vez existen varios campos de fuerza, que en general pueden utilizarse con diversos programas. CHARMM, AMBER, NAMD y GROMACS son algunos de los paquetes más populares.

Referencias

  1. Feynman, Richard (1963). Lectures on Physics. 

Véase también

Enlaces externos


Wikimedia foundation. 2010.

Mira otros diccionarios:

  • Dinámica Química — Saltar a navegación, búsqueda La Dinámica Química es una parte de la Química Física que estudia las interacciones entre las partículas que conforman los diferentes sustancias que intervienen en una reacción química y también los estudios sobre… …   Wikipedia Español

  • Modelado molecular — El Modelado molecular es un término general que engloba métodos teóricos y técnicas computacionales para modelar o imitar el comportamiento de moléculas. Las técnicas son utilizadas en los campos de la Química computacional, Biología… …   Wikipedia Español

  • Hamiltoniano molecular — Este artículo o sección necesita ser wikificado con un formato acorde a las convenciones de estilo. Por favor, edítalo para que las cumpla. Mientras tanto, no elimines este aviso puesto el 14 de septiembre de 2011. También puedes ayudar… …   Wikipedia Español

  • Imagen molecular — Las Imágenes Moleculares son la formación de imágenes a nivel molecular que se origina del campo de la radio farmacología debido a la necesidad de un mejor entendimiento de los caminos (o trayectorias) moleculares fundamentales en los organismos… …   Wikipedia Español

  • Wikiproyecto:Biología celular y molecular — Atajo PR:BCMPR:BCM Wikiproyecto: Biología Celular y Molecular …   Wikipedia Español

  • Máquina molecular — Una máquina molecular, o nanomáquina,[1] se ha definido como un discreto número de componentes moleculares que han sido diseñados para llevar a cabo movimientos mecánicos (de salida), en respuesta a determinados estímulos (entrada).[2] A menudo… …   Wikipedia Español

  • Física molecular — compuestos formados por moleculas.. La Física Molecular estudia problemas relacionados con la estructura atómica de la materia y su interacción con el medio. Por ejemplo, se tratan problemas como dinámica de reacciones, dispersión, interacciones… …   Wikipedia Español

  • Movimiento (física) — Para otros usos de este término, véase Movimiento. El movimiento es un cambio de posición respecto del tiempo. En mecánica, el movimiento es un fenómeno físico que se define como todo cambio de posición en el espacio que experimentan los cuerpos… …   Wikipedia Español

  • Química computacional — Saltar a navegación, búsqueda La química computacional es una rama de la química que utiliza computadores para ayudar a resolver problemas químicos. Utiliza los resultados de la química teórica, incorporados en algún software para calcular las… …   Wikipedia Español

  • Folding@home — Saltar a navegación, búsqueda El programa mostrando una simulación 3D Folding@home (en inglés se puede traducir como: «plegando en casa») es un proyecto de computación distribuida diseñado para realizar simulaciones por ordenador de plegamiento… …   Wikipedia Español