How can i select all P atoms within 10 Angstrom of another P atom?

[Ibrahim Mohamed]

for each P atom in a lipid leaflet i want to select other P atoms within 10 Angstrom. how can i do this?

I have tried this:

local_p_atoms = u.select_atoms('(around 10 name P) and name P')
for f in u.trajectory:
print(local_p_atoms.names)

but it outputs and empty array

Thanks