Oi Mariana, tudo bem? Esses fatores ótimos dependem muito da base escolhida (tanto os estabelecidos da literatura quanto os que você fez)? Qual base você usou nesse caso?   E cálculos com efeito de solvente, sabe se esses fatores mudam muito?   E para cálculos com métodos de função de onda? Há dados na literatura e/ou você planeja fazer esse tipo de estudo para métodos como MP2 ou CCSD?

Parabéns pelo trabalho Mari! Vocês fizeram algum critério na escolha dos 10 funcionais ou os funcionais utilizados foram selecionados aleatoriamente?

Olá Mariana, bastante interessante o seu trabalho. Tenho algumas questões. 1) Como era de se esperar, a presença de hidrogênio leva a um desvio maior do modelo harmônico devido a alta anarmonicidade de estiramentos C-H e similares. Você não gerou "scaling factor" em função do tipo de vibração ou da região da frequência? Isso seria interessante. 2) Por que você não usou funcionais híbridos duplos como o B2PLYP? Esse funcional apresenta em geral uma precisão maior para geometrias e frequências harmônicas. 3) Como você lidou com as frequências com menos de 100 cm-1 para o cálculo da entalpia? Essas frequências em geral possuem um grande erro numérico. 4) Poderia me dar mais detalhes do cálculo anarmônico? Como você fez para utilizá-lo na correção do harmônico? Pergunto isso porque existem diversos protocolos para isso como VSCF e VPT2.  Além disso, o cálculo anarmônico gera sempre mais bandas, em particular na região acima de 2500 cm-1 (bandas de combinação/overtones).