def METHOD_NAME(self, initial_molecule):
    name = 'scan_density'
    self.children[name] = AMSMDScanDensityJob(
        name=name, 
        nsteps=self.nsteps[name],
        settings=self.settings,
        scan_density_upper=self.scan_density_upper, 
        molecule = initial_molecule,
        **self.kwargs,
    )
    return self.children[name]