Se realizaron simulaciones de dinámica molecular utilizando el modelo de energía de adhesión (AE) para investigar la morfología de crecimiento del trinitroresorcínato de plomo 2,4,6 (LTNR, estifnato de plomo) en vacío y en diferentes disolventes. Se calcularon las energías de adsorción de LTNR en los planos cristalinos (001), (110), (011), (020), (111), (200) y (201). Mientras tanto, se predijo la morfología cristalina en disolventes de etanol, tolueno, diclorometano, acetona, dimetilsulfóxido (DMSO) y agua a 298 K calculando las energías de interacción entre los disolventes y los planos cristalinos. Los resultados calculados muestran que la morfología de los cristales de LTNR en diferentes disolventes es significativamente diferente. En tolueno, las morfologías cristalinas de LTNR son planas, mientras que en disolventes puros de etanol, acetona y DMSO, el número de planos cristalinos aumenta y el grosor del cristal es mayor. En el agua, LTNR tiende a formar cristales tabulares, lo que es similar a los resultados experimentales. Tanto el análisis de la función de distribución radial (RDF) como el desplazamiento cuadrático medio (MSD) revelan que los enlaces de hidrógeno dominan las interacciones entre LTNR y las moléculas de disolvente. Las moléculas de disolvente con coeficientes de difusión más altos exhiben mayores tendencias de desorción de las superficies cristalinas, lo que puede reducir sus efectos inhibitorios en planos cristalográficos específicos. Sin embargo, no existe una correlación directa entre los coeficientes de difusión del disolvente y las tasas de crecimiento de los planos cristalinos, lo que sugiere que la cinética de adhesión en la superficie o las barreras de energía interfacial juegan un papel más crítico en el crecimiento cristalino.