if subset is None:
        subset = list( range( 1, len( molecule )+1 ) )