Abstract: Exclusive dimuon production in ultraperipheral collisions (UPC), resulting from photon-photon interactions in the strong electromagnetic fields of colliding high-energy lead nuclei, PbPb(gamma gamma) -> mu(+) mu(-) (Pb-(*Pb-)(()*())), is studied using L-int = 0.48 nb(-1) of root S-NN = 5.02 TeV lead-lead collision data at the LHC with the ATLAS detector. Dimuon pairs are measured in the fiducial region p(T,mu) > 4 GeV, vertical bar eta(mu)vertical bar < 2.4, invariant m(mu mu) > 10 GeV, and p(T,mu mu) <2 GeV. The primary background from single-dissociative processes is extracted from the data using a template fitting technique. Differential cross sections are presented as a function of m(mu mu), absolute pair rapidity (vertical bar y(mu mu)vertical bar), scattering angle in the dimuon rest frame (vertical bar cos v*(mu mu)vertical bar), and the colliding photon energies. The total cross section of the UPC gamma gamma -> mu(+) mu(-) process in the fiducial volume is measured to be sigma(mu mu)(fid) = 34.1 +0.3(stat.)+0.7(syst.) μb. Generally good agreement is found with calculations from STARlight, which incorporate the leading-order Breit-Wheeler process with no final-state effects, albeit differences between the measurements and theoretical expectations are observed. In particular, the measured cross sections at larger vertical bar y(mu mu)vertical bar are found to be about 10-20% larger in data than in the calculations, suggesting the presence of larger fluxes of photons in the initial state. Modification of the dimuon cross sections in the presence of forward and/or backward neutron production is also studied and is found to be associated with a harder incoming photon spectrum, consistent with expectations.