We investigate when shrinking or steady Ricci solitons admit certain classes of harmonic functions, in connection to their structure at infinity (number of ends). If time allows, we will also discuss integral bounds of some quantities involving the curvature, which appear in the classification of locally conformally flat shrinking Ricci solitons. This is joint work with Natasa Sesum.