def METHOD_NAME(self):
    bond_length = self.system.box_l[0] * \
        np.sqrt(3) / 4. / (self.diamond_params['MPC'] + 1)
    bond = espressomd.interactions.HarmonicBond(
        k=1.5, r_0=bond_length, r_cut=3)
    self.system.bonded_inter.add(bond)
    espressomd.polymer.setup_diamond_polymer(system=self.system,
                                             bond=bond,
                                             **self.diamond_params)
    self.system.time_step = 0.1
    self.system.cell_system.skin = bond_length / 4
    self.node_parts = self.system.part.select(
        type=self.diamond_params['type_nodes'])
    self.charged_mono = self.system.part.select(
        type=self.diamond_params['type_cM'])
    self.noncharged_mono = self.system.part.select(
        type=self.diamond_params['type_nM'])