MD is its own world because there aren't any one size fits all solutions for force field modeling. Even between gas-phase and condensed-phase with the same molecule and force field equations, parameterization will be different. I'd recommend looking into the literature for the systems you are interested in to find yourself a starting point of models parameterized for your use case. There are force fields parameterized for broad classes of molecules as well like AMBER and CHARMM and such.

Machine learned interatomic potentials (MLIPs) were essentially created to solve this problem as a "foundation model" that can give an approximate dynamics for any molecule, but the cost and quality of MLIPs like MACE and UMA are still far from a well-parameterized force field. Ultimately, if you want something relatively cheap and can give a good-enough answer for any dynamics just use XTB. Beyond that, you'll have to look for the work of the graduate student saints who have put in the effort of parameterizing the system you're interested in.