WebMO Taller 26 08 2019
WebMO Taller 26 08 2019
WebMO Taller 26 08 2019
QUÍMICA
TUTORIAL – TALLER QUÍMICA COMPUTACIONAL
WEBMO
Dr. Aldo F. Combariza
___________________________________________________________________________________
Introducción a WebMO
WebMO es una interface que permite hacer cálculos de química computacional, utilizando diferentes
herramientas, en linea.
Ingresa a WebMO utilizando una cuenta de “Invitado” (guest) utilizando la siguiente información:
URL: http://webmo.net/
Ingrese usando el enlace “Working Demo” ubicado en el menú izquierdo de la página de WebMO.
Se abrirá la siguiente página:
Ingrese al enlace “WebMO Demo Server” (Resaltado en la figura), utilizando los siguientes datos:
Username: guest
Password: guest
Una vez verificada la información por parte del servidor, aparecerá la página del “Administrador de
Trabajos de WebMO” (WebMO Job Manager).
Al entrar al servidor, se muestra una lista de los trabajos de simulación que están ya sea en ejecución
o terminados.
2. Crear estructuras de compuestos químicos en WebMO
Ubica el cursor sobre la pestaña “New Job”, en la parte superior izquierda y haz click en el enlace
“Create New Job”. Esto te llevará al editor de moléculas de WebMO, se mostrará la página “Build
Molecule”.
Como ejemplo, crea la molécula de agua H2O.
● Selecciona el modo “Build”, donde podrás adicionar átomos y enlaces, en la parte superior
del menú, o directamente desde el menú gráficode la izquierda, .
● Una vez en modo “Build”, aparecerá una línea de texto en la parte inferior que dice “Build
Mode - (click to add atom, click and drag = add atom & bond, letter = change atom)”, como se
muestra en la figura. Las posibles acciones en el modo “Build” se encuentran en la barra de la
izquierda.
● Puedes seleccionar un átomo directamente en el menú despleglable en la parte superior, o
puedes usar el teclado para escribirel símbolo del átomo e.g. “h” para el hidrógeno, “o” para el
oxígeno etc, o también puedes usar la tabla pperiódica que se despliega a través de la
herramienta .
● Con un click en el área de trabajo, se ubica el átomo seleccionado en la ventana. Haciendo
click entre un átomo y deslizando el cursor hasta un segundo átomo permite que se cree un
enlace; enlaces dobles y triples se crean de la misma manera, seleccionando nuevamente la
misma pareja de átomos.
●
●
●
●
●
● Si cometes un error, use “control-z” o “Edit - Undo”. Si no puedes deshacer la acción, o si
has hecho un desastre monumental, noo te procupes, simplemente use “control-n” or “File –
New”, para iniciar de cero.
● Una vez cmpletada la molécula, use “Clean-Up – Comprehensive” o el botón “Clean-Up”,
3. Una vez terminada la molécula de agua, experimenta con el editor de moléculas creando algunas
moléculas.
El Ciclohexano debería ser bastante sencillo de crear. Crea el anillo de átomos de carbonoy usa el
botón para adicionar los hidrógenos necesarios. Utiliza el botón de Mecánica Molecular, ,
para optimizar la estructura. Con que geometría se optimiza el ciclohexano?
Usa el botón de rotación, , para girar tu molécula, el botón para trasladarla y el botón
para magnificar/reducir. También puedes ajustar las longitudes de enlace, ángulos y ángulos diedros
utilizando usando el botón, , para seleccionar los átomos de interés y luego el botón de abajo,
ajuste de longitudes de enlace, , ajuste de ángulos de enlace, , o ajuste de diedros, , los
cuales aparecen dependiendo del número de átomos que selecciones.
Rota, translade y magnifica o reduce la escala de la molécula, de la manera que consideres necesario,
para obtener una orientación adecuada. Una vezestés conforme con la imagen, oprime ALT+Print
Screen en tu PC para capturar la pantalla y pegarla en un programa de edición de gráficos, editala y
copiala en este documento, tal como se indica.
Para crear la molecula de bromoeteno, CH2CHBr, ubica un átomo de carbono y haz click con el botón
derecho del ratón para seleccionar la hibridación “sp2”, oprime sobre el carbono y arrastra para ubicar
el segundo átomo de carbono, y oprime y arrastre desde el primer carbono al segundo para crear el
doble enlace. Ahora adiciona el átomo de Br (usa la tabla periódica, ya que si escribes “B” se
seleccionará el átomo de Boro) y los átomos de hidrógeno restantes, optimiza, rota, magnifica y copia
la imagenen la siguiente página. Continúa creando las moléculas de ácido acético y benceno.
Ahora intenta crear un dímero de agua (H2O)2. Usa mecánica molecular para optimizar esta
estructuray cópiala en la siguiente página.
Crees que el Campo de Fuerza de mecánica molecular utilizado incluye un término para describir el
enlace de hidrógeno?
Inserta las figuras de las moléculas creadas aquí:
Taller 2 and 3: Cálculos a nivel Semi-Empirico, HF, DFT y MP2
1. Optimización de la geometría
● Para efectos de comparación, iniciaremos con la molécula de agua. Inicia un nuevo trabajo
y crea la molécula de H2O.
● Oprime en “Choose engine” en el menú de la izquierda, u oprime en la imagen de la flecha,
. Selecciona “Gaussian” como la máquina de cálculo para ejecutar un trabajo de estructura
electrónica a nivel “ab initio and semi-empirical” utilizando el programa Gaussian.
● Ahora oprime sobre la flecha: . Esto nos lleva a la página “Configure Gaussian Job
Options”:
Coloca tu nombre o tus iniciales en el campo “Job Name”, seguido de H2O para identificar tu trabajo y
poder revisarlo una vez termine.
Inicialmente vamos a correr un cálculo HF/6-31G(d). Primero necesitamos optimizar la estructura del
agua hasta alcanzar su mínimo nivel de energía y luego calcularemos sus frecuencias vibracionales
(para verificar que tenemos un mínimo global y para observar sus modos vibracionales) Para hacerlo:
Para visualizar los resultados oprime sobre el nombre del trabajo (Taller-WebMO-H2O, en este caso).
Esto abre la ventana “View Job”. Para el ejemplo del agua se tiene:
Oprime “Raw Output”, en el menú de la izquierda para observar la salida del programa en una
ventana nueva. Revisa el texto para ver los valores calculados.
También puedes exportar la geometría de la molécula a un nuevo archivo (utilizando una variedad de
formatos), con el botón “Export Molecule”. Por ejemplo, el formato de coordenadas cartesianas
“Gaussian Cartesian” te da los valores de las posiciones de los átomos en Å.
Se requiere la descripción de la geometría de las estructuras optimizadas. Para hacer esto, oprime
sobre el botón, , en el menú de herramientas.
Oprimiendo sobre un átomo y luego manteniendo la tecla “Shift” y oprimiendo sobre otro átomo ,
tendrá como efecto el mostrar la distancia entre esos átomos (en Å), por ejemplo, la distancia entre
los átomos de H en la molécula de H2O es de 1.508 Å, tal como aparece en la parte inferior de la
pantalla de la ventana “Molecule Viewer”.
Oprimiendo sobre un átomo y manteniendo despues la tecla “Shift” y oprimiendo sobre otros dos
átomos te da el valor del ángulo (°), oprimiendo sobre un átomo adicional te dará el ángulo diedro (°),
que es el ángulo del cuarto átomo con respecto al plano formado por los otros tres átomos. Por
supuesto, la molécula de agua solo posee tres átomos y por tanto no existe un ángulo diedro.
Si te desplazas hacia abajo en la página “View Job” podrás observar un resumen de los cálculos
elaborados. Para el ejemplo actual, HF/6-31+G(d), para la molécula de agua tenemos:
Oprimiendo en el icono de la lupa, , frente a “IR Spectrum” en la sección “Vibrational Modes”, se
puede generar el espectro de IR que se predice para el H2O. Por ejemplo,
Los modos vibracionales pueden animarse, utilizando el icono, , ubicado frente a cada una de las
frecuencias calculadas.
HF
3-21G
6-31G(d)
6-311+G(d,p)
B3LYP
3-21G
6-31G(d)
6-311+G(d,p)
MP2
3-21G
6-31G(d)
6-311+G(d,p)
Selecciona los cuatro trabajos en el Job Manager y oprime el botón Spreadsheet (Hoja de Cálculo)
para obtener una comparación inmediata de los cuatro cálculos. La energía de Solvatación es el
cambio de energía cuando entre la molécula en fase gaseosa y la molécula en solución, i.e., Esolvation =
Esolution − Egas.
Elabora una tabla especificando el tipo de molécula, energía en fase gaseosa, energía en solución y
energía de solvatación.
Cual de los isómeros es mas estable en fase gaseosa? Cual isómero es mas estable en agua? Comente
acerca de las razones para que exista una diferencia en las energías de solvatación para los isómeros.
2. Vinyl alcohol conformers.
Build vinyl alcohol-0°, vinyl alcohol-180°, and acetaldehyde. When building these molecules, adjust
dihedral angles to result in planar, staggered conformations. Perform a Geometry Optimization,
Hartree-Fock 3-21G calculation adding the line “%Mem=16MB” as the first line of the input file to
reduce memory usage and speed up the calculation.
Vinyl alcohol transition state: Build an approximate transition state for the isomerization of vinyl
alcohol-0° and vinyl alcohol-180° as follows. View a geometry optimized vinyl alcohol job that
successfully ran, and choose New Job Using This Geometry. Adjust the H-O-C=C dihedral angle to 90°.
Perform a “Transition State Optimization”, Hartree-Fock 3-21G calculation.
Make a table with columns for conformation, energy (Hartree), relative energy (kJ/mol). The relative
energy should be calculated relative to the global minimum. Note that 1 Hartree = 2625.5 kJ/mol.
Include both stable conformations and the transition state in the table.
4. Silane elimination reaction transition state.
In this exercise, the transition state for the reaction SiH4 → SiH2 + H2 will be found. Here we will
“interpolate” a “QST2” guess for the transition state. We will also use the semi-empirical PM3 method
for increased speed, but at the cost of less accuracy. Start by building SiH4 and set up a PM3
optimization as a template for this calculation. Before submitting your job, preview and generate the
input file and edit the displayed text into the following (being very careful to include the blank lines):
#N PM3/3-21G Opt=QST2
01
Si
H 1 B1
H 1 B2 2 A1
H 1 B3 2 A2 3 D1
H 1 B4 2 A3 3 D2
A3 96.849005
A2 84.283836
A1 110.00006
D2 -99.87710
D1 -81.78112
B4 1.9999999
B3 1.9999997
B2 1.4799991
B1 1.4799999
01
Si
H 1 B1
H 1 B2 2 A1
H 1 B3 2 A2 3 D1
H 1 B4 2 A3 3 D2
A3 109.37081
A2 109.30689
A1 110.00006
D2 120.13988
D1 -120.0214
B4 1.4800000
B3 1.4800000
B2 1.4799991
B1 1.4799999
Describe how the transition state search algorithm for this exercise differs from the exercise 2. When
would you use each algorithm?
Workshop4: Solvation and/or Excited States
Solvation
1. Let’s revisit the dipole moment of water. Set up a B3LYP/6-31+G(d,p) Geometry Optimization job
for water. Look at the preview file and edit the command line to include a polarizable background
(PCM) by clicking on the “Advanced Options” tab on the “Configure Gaussian Job” page and selecting
“water” as a solvent. Then go to preview and generate the input file used. You could edit the
command line (the one starting #N) in the input file directly by adding SCRF=(PCM,Solvent=Water).
We also want to add POP=MK to the command line to perform a charge calculation.
What happens with the water dipole moment? From classical simulation we would expect a higher
dipole moment about 2.25 Debye compared to 1.85 Debye in the gas phase.
2. Switch from water to H2S by changing the letter O to S in the input file. Run with geometry
optimization first in vacuum and then with SCRF=(PCM, Solvent=dichloromethane). When you add
PCM you should chose a standard solvent with a similar dielectric constant as H2S (ie somewhere
between 8 and 9 could be reasonable choices, dichloromethane is 8.93). Comments? Does the
geometry change compared to that of the water molecule? What are the bond lengths and the bond
angle for H2S compared to H2O. The experimental gas phase dipole moment is 0.97 Debye, while one
seems to need about 2.17 Debye in the liquid phase to be able to fit the experimental properties.
3. Run a methanol B3LYP/6-31+G(d,p) Geometry Optimization job in vacuum with POP=MK and/or
POP=NBO. The charge distribution on the three CH3 hydrogens may not become completely
symmetric. You may rotate around the O-C bond and see if you can fix this. Compare the fractional
charges to those of a GROMACS/GROMOS united atom model (the CH3 group is here treated as one
united atom) with the charges (H:+0.398 , O:-0.574 and CH3:+0.176 ) and with an all atom model
Pereira et.al., J. Phys. Chem. A, 105, 1909 -1925, (2001) which gives H:+0.4214, O:-0.5943,
C:-0.1820,3*(H:0.1183)
Check the stability of your the results if you include PCM with solvent=methanol and POP=MK and/or
POP=NBO and compare them to the figures above; the dielectric constant of methanol is around 30.
4. The smallest amino acid glycine, NH2CH2COOH, has a zwitterionic isomer, in which the carboxylic
acid proton is transferred to the amino group. Build “neutral” glycine and perform a Geometry
Optimization calculation. Using the resulting geometry, perform a Molecular Energy calculation at the
same level of theory but using the Advanced Options to select water as the solvent. Repeat the gas
phase Geometry Optimization and aqueous Molecular Energy calculations for the zwitterionic
structure.
Select the four jobs in Job Manager and click the Spreadsheet button to obtain an immediate
comparison among the four calculations. Solvation energy is the energy change upon solvation, i.e.,
Esolvation = Esolution − Egas. Make a table with columns for molecule, gas phase energy, solution energy,
and solvation energy.
Which isomer is more stable in the gas phase? Which isomer is more stable in water? Comment on
the reason for the difference in solvation energies between the two isomers.
Excited States
The default excited state method in WebMO is CI singles (which is not all that good), nevertheless…
1. UV-Vis Spectra of Conjugated Aldehydes
Build and perform a standard (ground state) “Geometry Optimization” calculations for the following
molecules, for example, at HF/6-31G. Using the resulting geometry, then perform “UV-Vis Spectra”
calculations for each molecule.
● 2-propenal (Abs Max 209 nm)
● 2-butenal (Abs Max 221 nm)
● 2,4-pentadienal (Abs Max 251 nm)
The first calculated excited state (A″) corresponds to a forbidden n→π* transition and is therefore
extremely weak. The second excited state (A′) corresponds to an allowed π→π* transition and is
strongly absorbing. Compare the computed absorption maximum to the experimentally observed
maximum in a table with columns for molecule, computed maximum, experimental maximum,
and difference. Comment on any trends that are observed.
Which vibrational modes undergo the largest change upon electronic excitation? Offer an explanation
for your result, noting that S1 is pyramidal and that the S0→S1 electronic transition is an n→π*
transition.