nueva página del texto (beta)
Español (pdf)
Artículo en XML
Referencias del artículo
Traducción automática
Enviar artículo por email
Citado por SciELO 
Similares en
SciELO 
Permalink
Introducción
En las últimas décadas, los métodos de simulación molecular de la química cuántica y la dinámica molecular han sido empleados para estimar las propiedades fisicoquímicas de diferentes sistemas químicos con el fin de complementar y explicar los comportamientos observados en ciertos procesos experimentales (Mammino, 2006; Ng, 2017). Para llevar a cabo estos tipos de estudios computacionales es necesario el uso de computadoras de alto rendimiento que permitan realizar las simulaciones de los sistemas moleculares con un gran número de átomos y moléculas (Lui et al., 2016). Existen diferentes programas comerciales como por ejemplo GAUSSIAN (Foresman y Frisch, 2015) y Q-chem (Epifanovsky et al., 2021), y de libre acceso como GAMESS-US (M. S. Gordon y Schmidt, 2005; M. W. Schmidt et al., 1993), ORCA (Neese, 2012; Neese, 2018, Neese, 2022) y GROMACS (Lindahl, Hess, y Van Der Spoel, 2001; Berendsen, van der Spoel, y van Drunen, 1995), los cuales han sido empleados para explorar el comportamiento a nivel atomístico de diferentes sistemas moleculares en fase gaseosa y fase condensada (Edeleva, Van Steenberge, Sabbe, y D’hooge, 2021; Schwartz, 2002). Particularmente, el uso de estos programas involucra un largo proceso de entrenamiento de los usuarios para poder realizar las simulaciones moleculares.
La implementación de los métodos de simulación molecular no está incluida en los pensum de las carreras de química y dichos métodos son enseñados exclusivamente en los últimos semestres de las licenciaturas en química en Latinoamérica. Sin embargo, el uso de estas herramientas computacionales para la enseñanza de la química puede permitir a los estudiantes comprender e interpretar diferentes fenómenos y procesos químicos con el uso de un computador personal o de escritorio (Rodríguez-Becerra et al., 2020; Simpson et al., 2021; Quezada et al.,2017; Bedolla y Bermúdez, 2009; Parra et al., 2021; Saldívar-González et al.,2017). De hecho, a nivel de estructura molecular, con las herramientas computacionales es posible estimar propiedades espectroscópicas de moléculas y hacer la elucidación de sus estructuras moleculares para ser comparadas con las determinadas mediante técnicas espectroscópicas usadas para caracterizar los compuestos obtenidos en síntesis experimental (Snyder y Kucukkal, 2021). Adicionalmente, la inversión necesaria para equipar un laboratorio de docencia con instrumentos espectroscópicos para medir espectros Infrarrojos (IR), Raman y Resonancia Magnética Nuclear (RMN) es muy elevada y la química computacional es una alternativa para la enseñanza y compresión de estas técnicas en la elucidación de estructuras moleculares a nivel universitario. Basado en estos antecedentes, en este trabajo hemos realizado un procedimiento usando cálculos mecánico-cuánticos para la determinación de los espectros IR, Raman y RMN de un grupo compuestos orgánicos (Acetanilida, Acetaldehído, Acetato de Etilo, Ácido Acético, Ácido Benzoico, Ácido Fórmico, Anilina, Ciclohexanona, Etanol, Fenol, Formaldehido) que son usados comúnmente en el laboratorio experimental de química orgánica para demostrar el alcance de las simulaciones computacionales en la determinación de propiedades espectroscópicas de moléculas orgánicas. Para ello, los softwares libres Avogadro-1.2.0 (Hanwell et al., 2012) y Orca-5.0.3 (Neese, 2018) fueron utilizados para este procedimiento. Este procedimiento fue diseñado con la finalidad de incrementar el conocimiento de los estudiantes de la carrera de química en la elucidación de estructuras moleculares de compuestos orgánicos. En la explicación y descripción del procedimiento se utilizó la molécula de acetanilida y los datos de los espectros IR, Raman y RMN del resto de las moléculas orgánicas utilizadas fueron reportadas (ver información de soporte). Finalmente, este manuscrito representa una guía general para estimar los espectros IR, Raman y RMN teóricos de especies orgánicas mediante el uso de software libre.
Metodología empleada y detalles de las simulaciones
Software
Los cálculos mecánico-cuánticos se realizaron con el software Orca-5.0.3 (Neese, 2022), disponible de manera gratuita en la página web. Para realizar la visualización y construcción de los sistemas moleculares se utilizó el programa Avogrado-1.2.0 (Hanwell et al., 2012), el cual puede ser obtenido de la página web.
Construcción de los modelos moleculares
El programa Avogradro-1.2.0 (Hanwell et al., 2012) programa que permite manipular la posición de los átomos en el espacio y optimizar las estructuras moleculares antes de ejecutar los cálculos mecano-cuánticos. Los compuestos usados para la predicción de los espectros fueron construidos y optimizados con la opción del campo de fuerza MMFF94 (Halgren, 1996). Luego, las geometrías moleculares son almacenadas en formato XYZ (Vílchez, 2020). Este proceso fue logrado con el siguiente procedimiento usando una instalación del software en español:
a) En el programa Avogadro-1.2.0., seleccionar la herramienta de construcción que esta simbolizada con un lápiz en el menú del programa. En la casilla de la herramienta, seleccionar el átomo de carbono y desmarcar la opción de ajuste automático de átomos de hidrógeno.
b) Construir el esqueleto de la molécula.
c) Ir al menú “Construir” y seleccionar “Añadir hidrógenos”. Esto permite colocar el número de hidrógeno necesarios en los átomos.
d) Usar la herramienta de optimización automática, con la opción “Gradiente Conjugado” para optimizar la molécula por mecánica molecular usando el campo de fuerza MMFF94.
e) Ir al menú “Archivo” y guardar la molécula en el formato XYZ.
Optimización de las geometrías moleculares usando el programa Orca-5.0.3
El programa Orca-5.0.3 permite realizar cálculos para localizar mínimos de energía sobre la superficie de energía potencial y predecir la estructura en equilibrio del sistema molecular en estudio (Neese, 2022). Para este estudio fue utilizado el funcional B3LYP y el conjunto base 6-31+G(d,p) para la obtención de las propiedades estructurales de las moléculas en fase gaseosa (Lodewyk et al., 2011). Este funcional B3LYP con el conjunto base 6-31+G(d,p) ha sido empleado eficientemente en la optimización de moléculas orgánicas y la estimación de propiedades espectroscópicas con muy buenos resultados (Alecu et al., 2010; Zvereva et al., 2011; Tirado-Rives y Jorgensen, 2008; Kesharwani et al., 2015). Moléculas en fase condensada y efectos de solvatación no fueron considerados en estos cálculos ab-initio. El procedimiento fue realizado usando la molécula de acetanilida (Figura S1 en la información de soporte). En la optimización de esta molécula, los pasos empleados fueron los siguientes:
a) Importar la estructura de la molécula almacenada en el formato xyz en el programa Avogadro-1.2.0. El usuario debe ir al menú “Archivo”, seleccionar “Abrir” y buscar el directorio donde se guardó la estructura construida.
b) Ir al menú “Extensiones” y seleccionar Orca para generar la entrada con las instrucciones para desarrollar el cálculo. Hacer click con el mouse en la opción “Generar entrada de Orca”.
c) Modificar el archivo de entrada con cualquier editor de texto agregando la siguiente línea de instrucción: B3LYP 6-31+G(d,p) OPT DEFGRID2 TIGHTSCF D3BJ PAL8. El archivo de entrada se muestra en la sección S1 de la información de soporte.
d) Hacer clic en la opción “Generar…” y guardar el archivo “acetanilide-opt.inp” en una carpeta del ordenador de trabajo.
e) En el sistema operativo Windows, si Orca está instalado y agregado al PATH es suficiente con abrir una terminal en la carpeta donde esté el archivo y escribir la siguiente instrucción:
$ C:/Orca/orca.exe acetanilide-opt.inp > acetanilide-opt.out
f) Una vez finalizado el cálculo, verificar que la estructura molecular haya alcanzado la convergencia buscando la palabra “HURRAY” en el archivo de salida “acetanilide-opt.out” utilizando cualquier editor de texto.
Obtención de los espectros Infrarrojos y Raman
En esta parte, se muestra como pueden ser obtenidos los espectros infrarrojos (IR) y Raman de las moléculas orgánicas. La molécula de Acetanilida optimizada con el procedimiento descrito en la sección 2.3 es utilizada como ejemplo modelo para la estimación de los espectros. En estudios teóricos se ha observado que los cálculos mecánico-cuánticos que involucran la teoría del funcional de la densidad y distintos conjuntos bases han permitido predecir de forma aceptable las señales de los espectros infrarrojos de diferentes tipos de moléculas (Alecu et al., 2010; Kesharwani et al., 2015). En particular, el funcional B3LYP con la base 6-31G+(d,p) ha sido evaluado ampliamente para predecir las señales espectroscópicas de diferentes compuestos orgánicos de forma aceptable (Katari et al., 2017; Merrick et al., 2007). Con el módulo NumFreq del programa Orca-5.0.3 se pueden determinar los modos vibracionales teóricos para las estructuras optimizadas con la finalidad de generar los espectros Infrarrojo y Raman (Neese, 2022). De igual manera, el mismo funcional y conjunto base en fase gaseosa fue utilizado para los cálculos de frecuencias vibracionales de la acetanilida. Los sobretonos y combinación de bandas no fueron incluidos con este método de cálculos (Sandorfy et al., 2007; Beć et al., 2022; Mondragón, 2020). R47este caso, el procedimiento utilizado es el siguiente:
a) Importar la estructura de la molécula optimizada en el programa Avogadro. El usuario debe ir al menú “Archivo”, seleccionar “Abrir” y buscar el archivo en el formato *.xyz que contiene las coordenadas de la estructura molecular en la carpeta correspondiente.
b) Ir al menú “Extensiones” y seleccionar en Orca “Generar entrada de Orca”. Un ejemplo de la entrada es mostrado en la sección S2 de la información de soporte.
c) Hacer clic en la opción “Generar…” y guardar el archivo con el nombre “acetanilide-ir-raman.inp” en una carpeta del computador donde se realiza el trabajo. La línea “%elprop Polar 1 end” especifica la generación del espectro Raman con Orca-5.0.3.
d) En el sistema operativo Windows, abrir una terminal en la carpeta y usar la siguiente instrucción:
$ C:/Orca/orca.exe acetanilide-ir-raman.inp > acetanilide-ir-raman.out
e) Una vez finalizado el cálculo, en el archivo “acetanilide-ir-raman.out” se tiene una Tabla que muestra en resumen las señales IR y Raman de la molécula utilizada (ver sección S3 en la información de soporte).
f) Para extraer los datos necesarios para realizar las gráficas de los espectros IR y Raman del archivo “acetanilide-ir-raman.out” se puede utilizar el módulo “orca_mapspc” del programa Orca-5.0.3. En este caso, el usuario debe ejecutar las siguientes instrucciones en un terminal de Windows:
Para el espectro IR: $ orca_mapspc acetanilide-ir-raman.out ir -x0400 -x14000 -w20
Para el espectro Raman: $ orca_mapspc acetanilide-ir-raman.out raman -x00 -x14000 -w20
La ejecución de esta instrucción genera 2 archivos en la carpeta de trabajo: uno con la extensión .dat que contiene los datos para graficar y otro con la extensión *.stk que contiene las transiciones individuales (número de onda e intensidad). Usando el archivo con extensión .dat se construye la gráfica del espectro usando los programas QtGrace (Winter, 2021) o Gnuplot (Janert, 2016).
g) Para la generación del espectro Raman se utiliza la referencia de Lui et al. (Liu et al., 2020) (ver información de soporte).
Obtención de los desplazamientos químicos de H 1 y C 13 para la obtención de los espectros de RMN
En esta parte, se muestra el procedimiento para obtener los espectros de RMN para compuestos orgánicos usando los desplazamientos químicos obtenidos con el programa ORCA-5.0.3 (Neese, 2022) mediante el funcional B3LYP y el conjunto base 6-311+G(2d,p). Aquí, el conjunto base 6-311+G(2d,p) ha permitido obtener valores aceptables de los acoplamientos químicos de los compuestos orgánicos en fase gaseosa para la construcción de los espectros RMN de H1 y C13 (Lodewyk et al., 2011). En este sentido, para calcular las constantes de acoplamiento químico y generar el espectro RMN del compuesto se siguen las siguientes instrucciones:
a) Importar la estructura de la molécula optimizada en el programa Avogadro. El usuario debe ir al menú “Archivo”, seleccionar “Abrir” y buscar el directorio donde se almaceno la estructura.
b) Ir al menú “Extensiones” y seleccionar en Orca “Generar entrada de Orca”.
c) En esta parte aparece el módulo para crear la entrada. ¡Modificar la entrada agregando la línea “! B3LYP 6-311+G(2d,p) AutoAux DEFGRID2 TIGHTSCF D3BJ” en la parte superior del archivo para el cálculo de las constantes de acoplamientos químicos para la acetanilida (ver sección S4 en la información de soporte):
d) Hacer clic en la opción “Generar…” y guardar el archivo con el nombre “acetanilide-spin.inp” en una carpeta del computador.
e) Ejecutar el cálculo en el programa Orca-5.0.3. Una vez finalizado el cálculo, el archivo “acetanilide-spin.out” contiene una Tabla resumen con las constantes de acoplamiento químico de la acetanilida en la sección “SUMMARY OF ISOTROPIC COUPLING CONSTANTS (Hz)” (ver sección S4 en la información de soporte).
f) Para generar el espectro RMN se necesita las constantes de acoplamiento químico de la acetanilida y los valores promedio de los desplazamientos químicos del tetrametilsilano (TMS, ver Figura S2 en la información de soporte) para los núcleos de hidrógeno y carbono a un nivel B3LYP/6-311+G(2d,p). Para ello, la entrada es modificada para realizar el cálculo y la construcción del espectro RMN (ver sección S4 en la información de soporte).
g) Hacer clic en la opción “Generar…” y guardar el archivo con el nombre “acetanilide-nmr.inp” en la misma carpeta donde se encuentra “acetanilide-spin.inp”.
h) Ejecutar el cálculo con el programa Orca-5.0.3. El archivo “acetanilide-nmr.out” contiene los valores necesarios para graficar el espectro RMN en la sección “NMR Spectrum” (ver sección S4 en la información de soporte).
Discusión de Resultados obtenidos
Geometría en equilibrio
Inicialmente, los valores de las longitudes de enlaces y los ángulos formados entre enlaces de las moléculas fueron comparados con los valores experimentales reportados con la finalidad de evaluar el alcance del cálculo mecánico-cuántico en el nivel B3LYP/6-31+G(d,p) para predecir las propiedades estructurales de las moléculas usadas en esta propuesta (ver Tabla 1 y 2). En esta parte, los valores de las longitudes y los ángulos formados entre los átomos estimados con este nivel de teoría fueron comparados con los datos obtenidos mediante difracción de rayos X (Brown y Corbridge, 1954) para el mismo compuesto con la intención de explorar el alcance del método DFT B3LYP/6-31+G(d,p) para predecir las propiedades estructurales de los compuestos utilizados.
Como se puede observar en la Tablas S1 y S2 (ver información de soporte), el cálculo mecánico-cuántico en el nivel de teoría B3LYP/6-31+G(d,p) permite predecir con exactitud las propiedades estructurales de los compuestos orgánicos como la acetanilida que fue usada como molécula modelo en este manual. Aquí, el valor del RMSE obtenido para las longitudes de enlaces tiene un valor menor en comparación con el RMSE obtenido para los ángulos entre enlaces. En este caso, el mayor error fue obtenido en la predicción de los ángulos de enlaces. Adicionalmente, para las propiedades estructurales del resto de compuestos orgánicos evaluados fue observado que los valores de RMSE de estas moléculas muestran la misma tendencia para el parámetro RMSE obtenido para la acetanilida. En este sentido, el cálculo mecánico-cuántico en el nivel de teoría B3LYP/6-31+G(d,p) permite predecir con exactitud las propiedades estructurales de las moléculas exploradas en este estudio. Al igual que el método DFT B3LYP/6-31+G(d,p), otros métodos DFT han sido evaluado para determinar su capacidad para predecir propiedades estructurales y energéticas con muy buenos resultados (Brémond et al., 2016), lo cual indica que el método empleado en este procedimiento no es el único capaz de estimar las propiedades estructurales de las moléculas.
Espectros IR y Raman
En la Figura S3(a) de la información de soporte, el espectro infrarrojo experimental (línea punteada en negro) es graficado en conjunto con el espectro infrarrojo obtenido por el cálculo mecánico-cuántico (línea roja). Claramente se observa como las intensidades de las señales en las ordenadas son menos profunda en comparación con los valores experimentales reportados. De igual manera, se puede observar como la posición de las señales calculadas están ligeramente desplazadas en el eje x con respecto a sus valores experimentales, lo cual indica la necesidad de usar un factor de escalamiento en la obtención de las frecuencias vibracionales de los espectros de infrarrojo. En este caso, los valores calculados para las frecuencias vibracionales fueron corregidas usando el factor de escalamiento obtenido por la gráfica que se muestra en la Figura S3(b) en la información de soporte. El factor de escalamiento fue obtenido por regresión lineal ajustando las frecuencias vibracionales calculadas con sus valores experimentales. Este factor permite corregir las desviaciones que tengas las señales estimadas con respecto a las señales experimentales. Por lo general, se utiliza un conjunto de moléculas diferentes (n>50) para obtener un único valor de escalamiento para el conjunto base y funcional utilizado en los cálculos mecánico-cuánticos. Según el ajuste lineal realizado en este estudio, el factor de escalamiento obtenido tiene un valor promedio de 0.9679 para el nivel de teoría B3LYP/6-31+G(d,p) en fase gaseosa (Merrick, et al., 2007). Estimar un buen factor de escalamiento es importante para los estudiantes que desean hacer estimaciones de espectros infrarrojos mediante simulaciones mecánico-cuánticos. La Tabla 3 muestra una comparación entre los valores de frecuencia vibracional para cada uno de los grupos más representativos presentes en la acetanilida. Estas frecuencias vibracionales corresponden a los enlaces N-H, el doble enlace C=C, el enlace C-H presente en el anillo aromático y el enlace C=O del grupo carbonilo.
Como se puede observar en la Tabla S3 (ver información de soporte), el cálculo mecánico-cuántico en el nivel de teoría B3LYP/6-31+G(d,p) predice de forma aceptable las frecuencias vibracionales de los grupos más representativos de la molécula de acetanilida, lo cual indica que este método tiene una buena capacidad para determinar las señales correspondientes a los grupos más representativos de una molécula para el espectro infrarrojo. Esto se muestra en la información de soporte que contiene los valores de las frecuencias vibracionales para 10 moléculas orgánicas adicionales.
Por otra parte, con las frecuencias vibracionales calculadas con el método DFT en el nivel B3LYP/6-31+G(d,p) se determinaron las señales correspondientes al espectro Raman de la acetanilida. Este espectro Raman, estimado con los cálculos computacionales, fue comparado con el espectro experimental (Saito, et al., 2022) en fase sólida de la misma molécula de acetanilida (Figura S4 en la información de soporte). En la Figura S4, se puede observar que las señales teóricas difieren de las reportadas experimentalmente. Las señales del Raman fueron determinadas directamente de las señales obtenidas para el espectro infrarrojo de la molécula en fase gaseosa. En este caso, las intensidades de las actividades son determinadas a partir de las frecuencias vibracionales del infrarrojo usando la ecuación 1 a una longitud de onda y temperatura definida por el usuario (ver sección 2.4). Además, el espectro Raman experimental está sujeto a diferentes efectos que computacionalmente no son tomados en cuenta. Experimentalmente, estos efectos corresponden a la excitación del láser, la resonancia y absorciones electrónicas (Colthup, 2012). De igual manera, se puede incluir el efecto del estado de agregación del compuesto (Alajtal et al., 2010; Tang y Albrecht, 1968). Debido a esto, el espectro Raman obtenido por el método mecánico-cuántico tiene una considerable desviación con respecto al espectro Raman experimental. Para mejorar la predicción de las intensidades de las señales es necesario hacer ajustes a las variables involucradas en la ecuación 1 (ver información de soporte). También, se puede ejecutar un procedimiento para la estimación de un factor de escalamiento si se conocen los valores experimentales de las intensidades de las señales de los compuestos utilizados para correlacionar las intensidades de las señales estimadas con las señales experimentales.
Desplazamientos químicos de RMN de H 1 y C 13
La RMN es una de las técnicas experimentales que permite obtener información muy importante sobre la estructura molecular de un compuesto orgánico. Por esta razón es importante dar a conocer la capacidad que tiene la química computacional para obtener dichos espectros. Los espectros de RMN de hidrógeno y RMN de carbonos de la acetanilida son reportados en las Figuras S5 y S6, respectivamente (ver información de soporte). Como se observa en la Figura S5, el espectro RMN de H1 calculado para la acetanilida muestra el número de señales que corresponden con el número de hidrogeno presentes en la molécula. La acetanilida presenta en su estructura un total de 9 hidrógenos. Sin embargo, en el espectro calculado aparecen un total 5 señales, lo cual es debido al efecto de acoplamiento spin-spin presente en los hidrógenos de carbonos vecinos (Laurella, 2017). Los hidrógenos del grupo acetilo corresponden a los hidrógenos H16, H17 y H18 y los hidrógenos H13, H11, H15, H12 y H14 corresponden al anillo aromático. Al comparar las señales calculadas para el espectro RMN de H1 con los valores experimentales (Saito, et al., 2022) se obtuvo un RMSE equivalente a 0.56. Esto demuestra que el funcional B3LYP con el conjunto base 6-31+G(d,p) permite predecir las señales de H1 con una muy buena correlación (ver Tabla S4 en la información de soporte). Para esta molécula, la mayor diferencia fue encontrada para el hidrógeno del grupo -NH, lo cual se puede atribuir a la posible formación de dímeros de la molécula en fase gaseosa (Lodewyk et al., 2011). Además, los efectos asociados a las interacciones intermoleculares no fueron tomadas en cuenta para el funcional utilizado (Kim y Jordan, 1994; Stephens, et al., 1994). En la Figura S6 en la información de soporte se observa el espectro RMN de C13 calculado en el nivel de teoría B3LYP/6-31+G(d,p) para la acetanilida. Aquí, se muestra que cada señal corresponde al número de carbonos presentes en la molécula. Además, al comparar con las señales experimentales (ver Tabla S5 en la información de soporte) existe una mayor desviación con la señal asociada al carbono del anillo aromático unido al nitrógeno (C1-NH). Los protones presentes en los átomos unidos al átomo de C1 puede ocasionar el incremento en el desdoblamiento de la señal de C13 de este átomo. En este caso, se observa que el átomo de nitrógeno adsorbe densidad electrónica del entorno que rodea al átomo de hidrógeno en el grupo -NH, lo cual ocasiona una mayor desprotección del este átomo de hidrógeno que se traduce en un aumento en el desplazamiento químico del átomo de carbono C1 (Carey, 2007).
Finalmente, la Tabla S5 muestra que la estimación de los desplazamientos químicos para el espectro de RMN de C13 tiene un RMSE igual a 4.81. Este valor es mayor en comparación con el RMSE obtenido para el espectro de RMN de H1, lo cual indica una mayor desviación en la estimación de las señales para el espectro de RMN de C13 utilizando los cálculos mecánico cuántico en el nivel de teoría B3LYP/6-31+G(d,p). Sin embargo, la predicción de las señales es aceptables en ambos casos. “Adicionalmente, este procedimiento fue empleado sobre 10 moléculas orgánicas adicionales y las señales espectroscópicas fueron estimadas de forma aceptable en este estudio para el laboratorio. Los datos de los espectros IR, Raman y RMN del resto de las moléculas orgánicas utilizadas fueron reportadas (ver información de soporte), en donde se muestran los valores de RMSE de cada una de las moléculas estudiadas.”
Limitaciones y problemas
Este procedimiento fue desarrollado utilizando moléculas que están formadas por muy pocos átomos (número de átomos < 25) y para un conjunto base intermedio que garantiza que pueda ser llevado a cabo en computadoras de escritorio con procesadores con una capacidad media como los procesadores tipo Intel i3, i5 e i7. Esto permite obtener resultados para las moléculas empleadas el mismo día de la ejecución de los cálculos. De igual manera, el uso del método DFT B3LYP/6-31+G(d,p) permite obtener buenas predicciones para los diferentes espectros evaluados en estas moléculas con pocos átomos, lo cual ha sido reportado en estudios previos (Merrick et al., 2007; Brémond, et al., 2016). Sin embargo, si se desea realizar este procedimiento en moléculas con un gran número de átomos es recomendable utilizar un conjunto base más grande y otros métodos post-Hartree Fock que reduzcan al mínimo el efecto de la base en la correlación electrónica. Sistemas moleculares con un gran número de átomos requieren de mayor poder de cálculo computacional para evaluar las propiedades espectroscópicas de una molécula en particular. Finalmente, para desarrollar este procedimiento es recomendable tener una información experimental amplia de las diferentes señales espectroscópicas de las moléculas orgánicas que se requieran estudiar con la finalidad de hacer la validación del procedimiento empleado para la predicción de las propiedades espectroscópicas de interés.
Conclusiones
Los cálculos mecánico-cuánticos en el nivel de teoría B3LYP/6-31+G(d,p) han permitido predecir de manera aceptable las señales de los espectros IR, Raman y RMN de H1 y C13 para compuestos orgánicos que son utilizados en el laboratorio de química orgánica. Aquí, el programa ab-initio Orca-5.0.3 permite desarrollar los cálculos mecánico-cuánticos a partir de un archivo de entrada que puede ser construido de una manera muy simple. Esto representa una alternativa para la enseñanza de las técnicas espectroscópicas a los estudiantes universitarios en el laboratorio de química orgánica. Las propiedades estructurales calculadas para la acetanilida tales como la longitud y ángulo de enlace se obtuvieron con gran exactitud al comparar con la data experimental. Las señales asociadas a los espectros Infrarrojos presentan menor desviación en comparación con el espectro Raman, lo cual es debido al uso del factor de escalamiento para ajustar los valores obtenidos de los cálculos mecánico-cuánticos en el espectro infrarrojo. En cambio, las señales del espectro Raman están sujetas a los términos presentes en la ecuación 1 (ver información de soporte) para la estimación de las señales (longitud de onda y temperatura). Del mismo modo, las variaciones en las señales de los espectros infrarrojos y Raman pueden ser atribuidas a la realización de los cálculos mecánico-cuánticos en el vacío para las moléculas seleccionadas, en donde no se toman en cuenta los efectos de solvatación de la molécula, formación de dímeros, la resonancia y absorciones electrónicas (Skyner, et al., 2015). Por su parte, los valores de las señales de los espectros RMN de H1 y C13 fueron obtenidos con una muy buena estimación basado en los valores de RMSE reportados. En esta parte, las desviaciones más importantes para los desplazamientos químicos fueron obtenidas para los átomos unidos al nitrógeno de la molécula de acetanilida. En este sentido, La química computacional es una herramienta fundamental para la enseñanza de la química que permite el uso de una simple computadora para explorar propiedades de diferentes moléculas (Alvarez, et al., 2019).
