  The process e+e- -> gamma gamma (gamma) is studied using data collected by
the OPAL detector at LEP between the years 1997 and 2000. The data set
corresponds to an integrated luminosity of 672.3pb-1 at centre-of-mass energies
lying between 181 GeV and 209 GeV. Total and differential cross- sections are
determined and found to be good agreement with the predictions of QED. Fits to
the observed angular distributions are used to set limits on parameters from
several models of physics beyond the Standard Model such as cut-off parameters,
contact interactions of the type e+e-gamma gamma, gravity in extra spatial
dimensions and excited electrons. In events with three photons in the final
state the mass spectrum of photon pairs is investigated. No narrow resonance
X-> gamma gamma is found and limits are placed on the product of the Xgamma
production cross-section and branching ratio.
