Esta dissertação investigou o comportamento de propriedades termodinâmicas de moléculas aromáticas utilizando o método de dinâmica molecular. Foram determinadas sete propriedades termodinâmicas para o benzeno e naftaleno nas suas fases líquidas, a saber: massa específica, compressibilidade isotérmica, coeficiente de expansão térmica, capacidade térmica a pressão e a volume constante, coeficiente de Joule-Thomson e velocidade do som. Para o benzeno, foi realizado um estudo na região crítica com a determinação do ponto crítico para cada um dos campos de forças utilizados, OPLS e TraPPE. O ponto crítico foi encontrado com auxílio do estimador Mustistate Bennett Acceptance Ratio (MBAR) e com um método de busca implementado no algoritmo R2W. O campo de força TraPPE foi mais preciso para encontrar o ponto crítico do benzeno, assim como também levou aos menores erros quando comparados com dados experimentais. A metodologia utilizada permitiu determinar valores de massa específica com erros inferior 0,0025% em toda fase líquida e, no máximo, 12% para as demais propriedades calculadas. Para o naftaleno, apesar da falta de dados experimentais para uma comparação completa, os resultados estão dentro das ordens de grandeza esperadas e são próximos dos dados experimentais disponíveis. A metodologia proposta é geral, logo pode ser utilizada para determinar propriedades de hidrocarbonetos aromáticos de maior tamanho molecular, como as presentes em piches de petróleo e derivados, para o cálculo acurado de propriedades termodinâmicas desses sistemas que não estão disponíveis tanto teórico como experimentalmente.