Desenvolvemos uma metodologia que nos permite estudar processos Raman, a princípio, de primeira ordem, através de uma abordagem puramente ab initio com o uso da teoria do funcional da densidade (DFT). Em nossos cálculos, a intensidade Raman é obtida através da derivada da susceptibilidade elétrica de um sistema com respeito aos seus modos normais de vibração, fazendo uso do bem conhecido código computacional SIESTA. Nossos resultados preliminares mostram que esta metodologia pode ser aplicada ao estudo de processos de primeira ordem em mono, bi e tricamandas de grafeno, além de outros sistemas de interesse, desde que sejam respeitadas certas condições.