Describe the bug
The documentation specifies the inputs of Topology.chemical_environment_matches specifies that the aromatic_model can be set, but this value isn't used in the code.
Suggested Solution
Pass aromaticity_model and propagate through to the toolkit methods find_smarts_matches in the method FrozenMolecule.chemical_environment_matches()