Il y en a un certain nombre, mais aucun n'est simple.
Si tu veux du simple, et que tu as un algo de résolution de déterminants (ce que tu es censé savoir faire) il ne te faut que quelques lignes pour implémenter ta méthode de Huckel.
La "computational chemistry" n'est pas un domaine simple, et loin d'être abordable en prépa.
Il faut le plus possible éviter les softs qui utilisent des méthodes que tu ne connais pas, et qui sont à des années lumières du programme de prépa, que ce soit dans les calculs menés, les approximations réalisées, ou sur le domaine numérique.
Un soft comme Hulis cité ci-dessus, de même qu'un bon solver "simple" et célèbre EHT (ci-dessous) n'utilisent pas des notions vues en prépa, que tu ne sauras pas justifier, expliquer ni présenter à l'oral.
https://comp.chem.umn.edu/eht/
Typiquement, la charge d'un atome est implémentée en prenant la théorie de Mulliken, et les dérivées de l'énergie avec Hellmann-Feynman. (Et c'est du Fortran 77, donc pas un language vu en prépa. )
On est très loin du calcul de prépa, que tu dois être capable d'implémenter sur n'importe quel outi de calcul basique à la Matlab, ou directement en Python comme c'est ton programme.
Quoi qu'il en soit, il est indispensable de n'utiliser un soft qu'en respectant ses citations et en le comprenant dans son intégralité, y compris à l'échelle du calcul qui est mené dedans.
En gros, si tu n'as pas lu une très longue liste de références et de cours sur le sujet, c'est inadéquat.
(typiquement : la moitié des refs citées dans
https://dasher.wustl.edu/chem478/ )
Enfin, pour le principe de le citer,
Avogadro reste un excellent outil
https://jcheminf.biomedcentral.com/arti ... -2946-4-17
https://avogadro.cc/
Il est scriptable en python et permet de faire beaucoup plus que du Huckel, mais donne, entre autres, les orbitales moléculaires du système qu'on dessine dessus.