Espectros Raman en Gaussian

Question

Quiero simular espectros Raman de benceno en Gaussian 09.

Espectros calculados para un benceno de anillo simple ( $\ce{C6H6}$ ) en Gaussian difiere alrededor de $\pu{40 cm^{-1}}$ para cada pico, en comparación con las bases de datos Raman experimentales.

¿Cómo puedo resolver este problema? Archivo gaussiano creado con GaussView6: MP2 config

DFT config:

Parámetros: (para optimización)

Answer 1

Gracias por adjuntar su archivo de entrada. El primer punto de partida para mejorar estos cálculos es la elección del conjunto de bases. 6-31G es un conjunto de bases bastante pequeño y es probable que sea una fuente de error significativa aquí. Si te gusta utilizar la base de tipo Pople, entonces sería mejor que utilizaras algo como 6-31++G* para que incluyas las funciones de polarización y difusa. Con tu base actual sospecho que te estás perdiendo una polarización importante de los átomos al distorsionarse.

Si realmente quieres una respuesta muy precisa te recomendaría el siguiente esquema. Optimiza tu geometría a nivel de DFT y luego usa esa geometría como punto de partida para una optimización y cálculo de frecuencias a un nivel superior de teoría. Por ejemplo, utilice el conjunto de bases cc-pvtz (un conjunto de bases coherente con la correlación que, en general, es más eficiente que el tipo Pople) y CCSD(T). Esto proporcionará un espectro raman en fase gaseosa muy preciso.

Edita: El único parámetro adicional específico de Gaussian que debería necesitar es el comando Geom=Checkpoint para la optimización CCSD(T) y la ejecución de frecuencias. Esto utilizará su archivo de punto de control existente como un punto de partida geometría. Las otras palabras clave (base y método) son exactamente las que deberían ser cc-pvdz y ccsd(t). Para más información sobre las palabras clave de base recomiendo el manual del usuario.