def lammps_pos_text_charged(structure, charges, species):
    """Create LAMMPS position text for a system of charged particles."""
    pos_text = '\n'
    for count, (pos, chrg, spec) in enumerate(zip(structure.positions,
        charges, species)):
        pos_text += '%i %i %f %f %f %f \n' % (count + 1, spec, chrg, pos[0],
            pos[1], pos[2])
    return pos_text
