rdf = RadialDistribution(md_trajectory,
                         cutoff_radius=20.0*Angstrom,
                         pair_selection=[[0], None])