Basically, in certain circumstances axion-like particles (ALPs) substantially enhance the photon survival probability $P_{ \gamma \to \gamma} ({ \cal E})$ of a beam emitted by a far-away source through the mechanism of photon-ALP oscillation (${ \cal E}$ denotes the energy). But in order for this to work, an external magnetic field ${ \bf B}$ must be present. In several cases ${ \bf B}$ is modeled as a domain-like network with `sharp edges': all domains have the same size $L_{ \rm dom}$ (set by the ${ \bf B}$ coherence length) and the same strength B, but the direction of ${ \bf B}$ changes randomly and abruptly from one domain to the next. It is obviously a highly mathematical idealization wherein the components of ${ \bf B}$ are discontinuous across the edges (whence the name sharp edges). It is therefore highly desirable to go a step further, and to find out what happens when the edges are smoothed out, namely when the abrupt change of ${ \bf B}$ is replaced by a smooth one. Moreover, this step becomes compelling when the photon-ALP oscillation length $l_{ \rm osc}$ turns out to be comparable to -- or smaller than -- $L_{ \rm dom}$, because in this case the photon survival probability $P_{ \gamma \to \gamma} ({ \cal E})$ critically depends on the domain shape. In the present paper we propose a smoothed out version of the previous domain-like structure of ${ \bf B}$ which incorporates the above changes, and we work out its implications. Even in the present case we are able to solve analytically and exactly the photon/ALP beam propagation equation inside a single smoothed-out domain, thereby evaluating the corresponding $P_{ \gamma \to \gamma} ({ \cal E})$ exactly. Our results is of particular importance in view of the new generation of gamma-ray detectors, since in such a situation $l_{ \rm osc} \lesssim L_{ \rm dom}$ occurs just above the TeV scale.
Categories: astro-ph.HE hep-ph