from openff.toolkit import Molecule, ForceField, Quantity
import openmm


topology = Molecule.from_smiles("CCO").to_topology()
topology.box_vectors = Quantity([3, 3, 3], "nanometer")  # same behavior with and without this

system = ForceField("openff-2.1.0.offxml").create_openmm_system(topology)

for force in system.getForces():
    if isinstance(force, openmm.NonbondedForce):
        print(f"{openmm.__version__=}")
        print(f"{force.getUseSwitchingFunction()}")
        print(f"{force.getSwitchingDistance()}")
        print(f"{force.getNonbondedMethod()}")
        print(f"{force.getCutoffDistance()}")
        continue