Estoy tratando de hacer un nuevo campo de fuerza en AMBER para un nuevo compuesto que incluye boro (que no existe en GAFF). Para ello, he calculado la frecuencia y he encontrado las constantes de fuerza de los enlaces, pero no he podido encontrar las constantes de fuerza para ángulos .
Por ejemplo, he probado $\ce{H2O}$ : He calculado $\ce{H2O}$ frecuencias en Gaussian y obtuve las constantes de fuerza del archivo de registro. En el archivo de registro, el $\ce{O-H}$ La constante de fuerza de estiramiento fue $\pu{9.1465 mDyne/\require{mediawiki-texvc}\AA{}}$ y he convertido esto en aproximadamente $\pu{1315.65 kcal mol^{-1} \AA{}^{-2}}$ . En AMBER, $v = 0.5k(l-l_{o})$ así que lo he multiplicado por 0,5 y he encontrado
ow-hw 657.825que está muy cerca del campo de fuerza GAFF2:
ow-hw 612.98que se puede encontrar en línea 87 . Sin embargo, no pude encontrar una manera de calcular las constantes de fuerza de Gaussian para los ángulos. En GAFF2 para $\ce{H2O}$ los valores son
hw-ow-hw 43.276 104.520 AMBER 1 TIP3P_waterHe buscado y he encontrado muchas respuestas como este Pero en esta respuesta, para el enfoque del Seminario no sé qué valor poner en la ecuación de la matriz. ¿Podría proporcionar un ejemplo en un archivo de registro gaussiano?
4 votos
Hola a todos, creo que esto no debería contar como tarea. No sólo porque es poco probable que se pida como tarea, sino porque la traducción de los cálculos químicos cuánticos a campos de fuerza simples es muy deseable.