We give two double copy prescriptions which construct asymptotically flat
solutions in gravity from asymptotically flat gauge fields. The first
prescription applies to radiative fields, which are non-linear vacuum solutions
determined by characteristic data at null infinity. For any two such radiative
gauge fields (linear or non-linear), the characteristic data of a radiative
metric, dilaton and axion is constructed by a simple `squaring' procedure,
giving a classical double copy at the level of radiation fields. We demonstrate
the procedure with several examples where the characteristic data can be
explicitly integrated; for linear fields this also sheds light on the
twistorial description of Weyl double copy. Our second prescription applies to
all asymptotically flat fields at the level of their asymptotic equations of
motion: we give a map between any solution of the asymptotic Maxwell equations
and any solution of the asymptotic Einstein equations at null infinity. This
also extends to the asymptotic charges and their duals, preserves the soft and
hard sectors between gauge theory and gravity, and is related to the usual
notion of double copy in scattering amplitudes.