Scattering of compact oscillons

We study various aspects of the scattering of generalized compact oscillons in the signum-Gordon model in (1+1) dimensions. Using covariance of the model we construct traveling oscillons and study their interactions and the dependence of these interactions on the oscillons’ initial velocities and their relative phases. The scattering processes transform the two incoming oscillons into two outgoing ones and lead to the generation of extra oscillons which appear in the form of jet-like cascades. Such cascades vanish for some values of free parameters and the scattering processes, even though our model is non-integrable, resemble typical scattering processes normally observed for integrable or quasi-integrable models. Occasionally, in the intermediate stage of the process, we have seen the emission of shock waves and we have noticed that, in general, outgoing oscillons have been more involved in their emission than the initial ones i.e. they have a border in the form of curved worldlines. The results of our studies of the scattering of oscillons suggest that the radiation of the signum-Gordon model has a fractal-like nature.


JHEP01(2020)006
do not radiate at all and so they can live forever. Such infinitely long living objects are called "breathers" which distinguishes them from standard oscillons. In this paper we describe our study of another interesting group of oscillons which share many properties of oscillons and breathers. A solution of such a type was first discovered in the signum-Gordon model [24]. The signum-Gordon oscillon is an exact solution which, if not perturbed, behaves like a breather i.e. it can live forever without emitting any radiation. This is a very interesting and extremely rare behaviour for time dependent solutions in nonintegrable field theories. Of course, due to the non-integrability of the model a perturbed signum-Gordon oscillon would emit some radiation. Such radiation often takes the form of emissions of smaller oscillon-like packages. So, this type of an oscillon can be thought of as being a stable (or perhaps metastable) time dependent non-topological solution of a non-integrable model. Moreover, some very special perturbations of such oscillons lead to more general, exact and infinitely long lived oscillons (generalized oscillons). Such oscillons were constructed in [25,26].
The signum-Gordon model [27] is perhaps the simplest example of a wider class of scalar field-theoretic models with non-analytic potentials. A very important and characteristic property of such models is their possession of compactons [28][29][30][31][32][33] and scaling symmetry [34]. This symmetry makes these models relevant in the description of dynamics of fields in other models with approximate scaling symmetry in the limit of small amplitudes [35][36][37][38]. In other words, the signum-Gordon model can be thought of as emerging, in this limit, from models containing non-analytic potentials. This shows that studies of solutions of this model are useful and can have relevance in the description of some aspects of solutions of other models with non-analytic potentials. Of course, due to an often encountered rich structure of minima of such more general models, the field configurations with larger amplitudes could also have some nontrivial topology (kinks, skyrmions etc.) [5,28,35,39] and so their complete dynamics would be essentially different from the dynamics of the signum-Gordon compactons.
In what follows we present some basic notions about the signum-Gordon model. The model is defined by the Lagrangian density L = 1 2 ∂ µ φ∂ µ φ − |φ| (1.1) and its dynamics is described by solutions of the Euler-Lagrange equation The Euler-Lagrange equations contain a term sgn(φ) = ∂ ∂φ |φ| = ±1 and so they do not include the vacuum solution φ = 0. In order to include explicitly the vacuum solution into the set of solutions of (1.2) we require that sgn(0) := 0. The model (1.1) has naturally appeared in the study of the behaviour of scalar fields in the vicinity of minima of V-shaped potentials i.e. potentials whose left and right side derivatives at minima are different from each other. Such models are perfectly well-defined from a physical point of view. Moreover, in some cases they can be seen as field-theoretic limits of certain mechanical models, which certainly admit experimental realizations. In fact, it was such mechanical models that led to scalar field models with non-analytic potentials [28].

JHEP01(2020)006
The physical origin of models with non-analytic potentials is wider than the continuous limit of mechanical models. Quite recently it was reported in [35] that models with Vshaped potentials may be obtained from other physical models when a parametrization associated with the symmetry reduction leads to new field variables that are restricted, i.e they cannot take arbitrarily large (or small) values. In the case of models with a mechanical realisation such a restriction is a priori imposed on the system. The restrictions on values of fields lead to some inconvenience in the description of the dynamics of the system which, in such a case, is governed by both the Euler-Lagrange equations and the extra condition on the time derivative of the field. For instance, the mechanical model, from which the signum-Gordon model originates, has a continuous limit described by the field variable that satisfiesφ ≥ 0. Thus the model possesses the potential with an infinite barrier at φ = 0 i.e.Ṽ (φ) =φ forφ ≥ 0 andṼ (φ) = ∞ forφ < 0. The field must also satisfy the reflection condition ∂ tφ → −∂ tφ atφ = 0. One can avoid such an inconvenient reflection condition by introducing an auxiliary model with a new field φ ∈ (−∞, ∞) and the potential V (φ) = |φ|. 1 This new model is so-called the unfolded model. The dynamics of this auxiliary field can be mapped onto the dynamics of the physical field through the folding transformation, see [27]. Thus the signum-Gordon model and other models of this type can describe behaviour of physical systems with restrictions on the values of scalar fields.
The signum-Gordon model is certainly non-integrable. This conclusion can be drawn from the existence of radiation in numerical simulations of generic initial field configurations. In fact, very little is known about the nature of this radiation. In this paper we describe results of our study of some aspects of the radiation in the signum-Gordon model. We pay particular attention to the exact time-dependent solutions of the model known as compact oscillons 2 [24,25], which are our principal candidates for constituents of the radiation. The signum-Gordon oscillons rely on three principal properties of models with V -shaped potentials: the existence of compact solutions, their scale invariance (exact or approximate) and the lack of linearization of small amplitude oscillations. The existence of compact solutions, like compact oscillons in particular, follows from the fact that models with standard kinetic and gradient terms in the Lagrangian approach vacuum in a quadratic manner if the potential has a V -shaped form close to its minimum [28]. The scale invariance [34] is a straightforward consequence of the form of the field equations. In the case of the signum-Gordon model the scale invariance is exact because sgn(φ) is a scale invariant term. A very important consequence of this fact is the existence the self-similar solutions and oscillons of all scales of energy and length. Thus the perturbed oscillons may lose energy by the emission of smaller (perturbed) oscillons. We will demonstrate in this paper that this is exactly what happens and so that the oscillons are main ingredients of the radiation of the signum-Gordon model. Finally, the absence of linear small amplitude oscillations follows from the fact that sgn(φ) term cannot be linearized at φ = 0. This 1 This potential generates the term sgn(φ) = dV dφ in the field equations which justifies the name of the model. 2 In the literature the infinitely long lived oscillons are more often called breathers than oscillons. Here we keep the name oscillon following the original paper [24].

JHEP01(2020)006
implies that, independently of their size, the signum-Gordon oscillons are fundamentally non-linear field configurations. The existence of many exact solutions of the signum-Gordon model in (1+1) dimensions follows from the fact that equation (1.2) reduces to a non-homogeneous wave equation on segments of the x axis, where the sign of the scalar field is fixed, and so, on each segment, it has the general solution of the form: where F and G are arbitrary functions. The main point here is that this reduction is local i.e. the size and the localization of the segments of constant sign varies with time. This is a direct manifestation of a non-linear character of the model. Solutions like (1.3) are called partial solutions. The exact global solution of the model is given by the explicitly known set of properly patched partial solutions which must hold for arbitrary times. Determining such a closed set of partial solutions usually generates some small technical difficulties. They correspond to the unpleasant side of finding solutions of models with V -shaped potentails.
The main aim of this paper is to describe the results of our study of interactions between two oscillons. Our motivation for such a study follows from the fact that the previous numerical simulations of models with V -shaped potentials [37] have indeed observed collisions of oscillon-like structures in emerging radiation. The study of systems containing colliding oscillons may throw new light on the nature of the radiation of the signum-Gordon model. In particular, we are interested in the scattering processes involving only two oscillons. Unfortunately, even such a 'simple' scattering process is too complicated for a purely analytical investigation. For this reason we have used the numerical analysis as our principal tool and have restricted our attention to initial field configurations which possess certain symmetries. We hope that our results will find applications not only for models with a single minimum potential but also for models with multi V -shaped minima that can support the existence of compact kinks etc.
The paper is organized as follows. In section 2 we present some facts about generalized oscillons with vanishing total linear momentum. Then, making use of the Lorentz covariance of the model, we construct traveling oscillons. In section 3 we study scattering of two oscillons. Due to the compactness of the oscillons we construct the initial configurations of such oscillons by considering simple superpositions of non-overlapping oscillons. Further on in this section we also discuss shock waves and oscillons with non uniformly moving borders. The last section presents a short summary of our results and some comments.

Generalized oscillons
The generalized oscillons are exact compact solutions of the signum-Gordon model which are distinguished by the fact that their borders move periodically from left to right and back again. The first particular solutions of this class, characterized by constant velocity of the motion, were reported in [25] and, due to this nature of the motion, were called swaying oscillons. A further generalization of the swaying oscillons to arbitrary periodic JHEP01(2020)006 motions of the borders was discussed in [26]. In order to simplify our initial discussions we first consider the field configurations describing oscilllons swaying with a constant velocity v. The more general oscillons are discussed later in section 3.5 where we also comment on the results of scattering of compact oscillons.
Before we go further let us establish the terminology for different types of oscillons considered in this paper. The oscillons which are a priori exact solutions (like initial field configurations) will be called the exact oscillons if their borders do not move periodically or, if they do move, the generalized exact oscillons. The oscillons produced in the process of the scattering are numerical solutions of the signum-Gordon equation, hence, they are not a priori exact. They will be called quasi-oscillons in the case of their strong similarity to the exact oscillons and perturbed oscillons if such a similarity is only approximate i.e. when they are approximately periodic, emit radiation etc.

Oscillons at rest
In our work here we are particularly interested in scattering of generalised exact oscillons. Such oscillons move uniformly (modulo a periodic motion given by v) in a certain reference frame S. We shall refer to this frame as to the laboratory reference frame. Moreover, we preserve symbols t and x for coordinates exclusively in this reference frame. On the other hand, the reference frame of the oscillon is denoted by S and called the rest frame of the oscillon. The precise meaning of expression "rest frame" in the case of generalised exact oscillons is the following one: it is the inertial reference frame in which the total linear momentum of the oscillon vanishes. The space time coordinates in S are denoted by t and x .
The basic oscillon has period T = 1. Smaller and bigger oscillons, which differ by their period, can be obtained from the basic one by the scaling transformation which exploits the scaling symmetry of the signum-Gordon equation, see section 3.1.1. The Minkowski diagram presented in figure 1 shows regions of validity of the partial solutions φ k (t , x ), k = {C, L 1 , L 2 , L 3 , R 1 , R 2 , R 3 } that together describe the motion of the oscillon. These partial solutions are given by quadratic polynomials in variables t and x in the rest frame of the oscillon. Physically relevant parts of the polynomial solutions are restricted to some intervals of t and x . In order to get a periodic solution for any t one has to replace t by a periodic function of t . Below we discuss this process in detail.
We start our discussion with a set of partial solutions which are valid only in the interval t ∈ 0, 1 2 . These solutions have been given in [25] and they are denoted by the letter ϕ. Here, we present them in a new notation. In fact, there are seven partial solutions in the interval t ∈ 0, 1 2 . Among them, four solutions are essentially different, namely 2) JHEP01(2020)006 Figure 1. The world sheet of the generalized exact oscillon seen in its own rest frame S . Only in the interval t ∈ 0, 1 2 φ k = ϕ k holds; other partial solutions are given by (2.9). The axes x of the laboratory reference frame seen in frame S have an inclination with respect to axis x . Here we present three different cases of inclination corresponding to the velocity u = −V of the laboratory frame S that moves to the left on the diagram S . The angle of this inclination between axes x and x is given by arctan(V ).
The other three partial solutions ϕ R j (t , x ; v), j = 1, 2, 3 can be obtained from those shown above by performing the transformation

JHEP01(2020)006
Note that the solution ϕ C (t, x; v) is invariant under this transformation. Note also that all solutions ϕ k (t , x ; v) are negative-valued in their domains. Each solution ϕ k (t , x ; v) is valid only in a specific region of the Minkowski diagram. For this reason we define a few region step functions Π k (t , x ; v) which are equal to unity only in the region in which the given partial solution holds and vanish outside this region: where θ(z) is the unit step function θ(z) = 0 for z < 0 and θ(z) = 1 for z ≥ 0. One can The parabolic functions like (2.1)-(2.4) are not periodic whereas the oscillon is a periodic solution. In order to give formulas which are valid for any t , we define two periodic functions τ (z) := 1 π arcsin(| sin(πz)|), (2.7) σ(z) := sgn(sin(2πz)), (2.8) where τ (z) maps any t onto the intarval [0, 1 2 ] and σ(z) = ±1 agrees with classical derivative of τ (z). The functions (2.7) and (2.8) allow us to construct the periodic solutions. They are presented in figure 2. The function σ(z) is needed to describe the changes of the sign of the partial solutions at t = n 2 . Any partial solution (for arbitrary t ) can be written in terms of basic solutions ϕ k (t , x ; v), where k = {C, L 1 , · · · , R 3 }, that are valid only on the interval t ∈ 0, 1 2 . The partial solutions presented in the Minkowski diagram in figure 1 have the form (2.9) The functions (2.7) and (2.8) allow us to introduce a more compact notation than the one in [37]. Note that here the functions φ − L 1 /R 1 and φ + L 3 /R 3 have been absorbed into the definitions of functions φ L 1 /R 1 . Similarly, φ − L 3 /R 3 and φ + L 1 /R 1 have been absorbed into the definitions of functions φ L 3 /R 3 . The total solution is given by a continuous function which is a sum of non overlapping partial solutions (2.9). The derivative of the oscillon solution with respect to time is given by  where σ 2 (z) = 1. Note that all the derivatives of the region step functions Π α can be ignored because the sum of partial solutions is a continuous function so there is no reason to expect delta functions at the matching points.
In figure 3 we present three snapshots of the generalized exact oscillon solution φ(t , x ; v) and its time derivative ∂ t φ(t , x ; v) at t = 0.1, t = 0.4 and t = 0.75. Solutions presented in figures (a) and (d) consists of (from left to right)

Travelling oscillons
The signum-Gordon equation is invariant under the Lorentz transformations. Thus, traveling compact oscillons can be obtained from the non traveling ones by an appropriate Lorentz transformation. In particular, the oscillon with a non vanishing linear momentum is obtained from the generalized exact oscillon by a Lorentz boost. In what follows we assume that the laboratory reference JHEP01(2020)006 frame S moves with velocity u = ∓V with respect to the rest frame of the oscillon S . Thus, the oscillon has velocity u = ±V in the laboratory frame S. The field configuration that describes the traveling oscillon is a function of coordinates t and x and is given by performing the transformations in (2.9) and obtaining the partial solutions The derivatives of the fields with respect to time t are given by the expression where d dz τ (z) = σ(z) and σ 2 (z) = 1 at open supports of the partial solutions. Hence, the travelling oscillon in S is a solution of (2.11) given by a sum of non overlapping partial solutions (2.13), namely, Note, that the axis x i.e the line t = 0 in S is not a simultaneity line in S. The scalar field φ(t , x ; v) vanishes at the horizontal lines t = n 2 , n = 0, ±1, ±2, . . . which are shown in figure 1. It shows that the oscillon seen in the laboratory reference frame S has some isolated traveling zeros. Such isolated zeros are given by points of intersection of lines parallel to the axis x (given by t = const) with the lines t = n 2 . The number of isolated zeros and the composition of the oscillon (types of partial solutions seen in S at given instant of time t) depend on the value of the velocity u = ±V which the oscillon has in S. The axes x and x form an angle arctan(V ) so for V < 1 2+v there is only one point of intersection of straight lines parallel to x with the line t = n 2 . For 1 2+v < V < 1, a second isolated zero arises at some value of the time interval.
According to the diagram presented in figure 1, the initial (at t = 0) configuration of the field in S, obtained by the Lorentz boost of the oscillon given in S , consists of different sets of partial solutions corresponding to different values of the velocity V . Let us consider the oscillon that moves with the velocity u = +V in the laboratory reference frame. The oscillon configuration ψ at t = 0 consists of (from left to right)   The compactness of exact oscillons allows the construction of some multi-oscillon configurations which are exact solutions of the signum-Gordon equation. The only condition to satisfy is the non-overlapping of the supports of individual oscillons. A generic initial configuration {Ψ(x), Ψ t (x)} containing two travelling oscillons is given by the superposition of non-overlapping travelling oscillons obtained from generalized exact oscillons φ(t, x; v 1 ) and φ(t, x; v 2 ) by the transformations which are symmetries of the signum-Gordon equation: • Poincaré transformatons in (1+1) dimensions: boosts, spatial and temporal translations, spatial and temporal reflections, • sign flipping of the field φ → −φ.
Two individual oscillons, i = 1, 2, are obtained by transformations where u i are boost velocities, t 0i are temporal translations, x 0i are spatial translations, λ i stand for scale parameters and ε i = ±1 and i = ±1 for reflections. Let us fix the sign of and so the sign of v can be absorbed into the combination of spatial reflexion and translation. The initial configuration {Ψ(x), Ψ t (x)} is now given by where the non-overlapping of their supports restricts the values of admissible spatial translations. It turns out that some of left-over parameters can be omitted without any loss of generality. The only relevant parameters of the initial configurations involving two oscillons are those which are not equivalent i.e. they cannot be related by a symmetry transformation. Thus the set of relevant free parameters contains relative scale, relative velocity etc. In some cases it will be more convenient fix one of them and so only study the dependence on the last one.
In figure 5 we present an example of the initial configuration containing two oscillons which move in opposite directions with velocities u 1 = 0.5 and u 2 = −0.7. This configuration was obtained by taking the generalized exact oscillons with v 1 = 0.3 and v 2 = 0.34 and λ 1 = λ 2 = 1. We have set t 0 = 0 and ε i = 1 as well as i = 1. The oscillons were shifted in space taking x 01 = 0 and x 02 = −1.3.

Symmetric configurations
A numerical study of the scattering process shows that the process depends on many parameters like v 1 , v 2 , the relative velocity of the oscillons, their initial distance, time shift and reflections. In order to simplify the set of parameters we have decided to restrict our considerations to symmetric and anti-symmetric initial configurations which certainly reduces the number of free parameters. It has turned out that even with this restriction we have been left with a sufficiently rich set of physical systems. We think that more general configurations are even interesting; however, their systematic study would have required much more work so we have decided to put main effort on symmetric configurations.  In order to get a symmetric Ψ (s) or an antisymmetric Ψ (a) configuration we can take a single exact oscillon parametrized by v and perform a sequence of symmetry transformations which leads to ψ(t + t 0 , x + x 0 ; v, u)| t=0 , where the boost velocity is chosen to be u = +V with V ≥ 0. The second oscillon can by obtained from this result by the spatial reflection x → −x and, optionally, by sign flipping. Naturally, x 0 must be chosen in a way that the support of ψ lies on negative semiaxis x. In such a case the support of the oscillon and its mirror image do not overlap. The resultant initial symmetric and antisymmetric configurations are given by Note that an alternative symmetric (antisymmetric) initial configuration can also be obtained by taking u = −V and shifting the resulting oscillon to the right i.e. choosing x 0 < 0. The second oscillon is obtained by the mirror reflection x → −x (and possibly the sign flip). An example of such an initial configuration is shown in figure 7 (shadowed regions). The configuration shown in figure 6 is marked in figure 7 by curves without shadowing.

Phase of the oscillon
Unlike for the case of compact kinks, the scattering process of two oscillons depends on the initial distance between them. This observation follows from the fact that the shape of the JHEP01(2020)006 oscillons changes with time and the outcome of the scattering process depends strongly 3 on the shapes of oscillons at the moment when their supports begin to collide. Thus the phase of the oscillon is another relevant parameter which must be taken into account in the analysis of the scattering of oscillons. Two traveling oscillons that differ exclusively by the value of spatial translation are said to have the same phase. The phase of the oscillon, in its own rest frame S , is the number α ∈ [0, 1), where the lowest value of α = 0 represents an oscillon configuration at t = 0 and the upper limiting value α = 1 = T rest (for λ = 1) is given by the period of the oscillon. The period of the oscillon in the laboratory reference frame S, in which the oscillon has a certain velocity V , is given by The phase of the oscillon in the laboratory reference frame S can be chosen again as the number α ∈ [0, 1), whose upper limit is given by T lab /γ. Note, that two oscillons with the same phase in two different inertial reference frames describe different field configurations. Below we describe in more detail the choice of the initial symmetric (antisymmetric) configurations containing two generalized exact oscillons.
The uniform motion of the oscillons from t = 0 to the moment of the collision results in a variation of their individual phases or the variation of their common phase in the case of symmetric (antisymmetric) initial configurations. The variation of the phase depends on the initial distance between support of two oscillons. Since the oscillons do not interact until they supports begin to overlap one can eliminate this initial distance without any loss of generality. This can be done choosing properly the value of spatial translation. The condition that oscillons begin to collide at t = 0 (their supports touch each other) makes the parameter of spatial shift a function of the phase i.e. x 0 = x 0 (α).
In order to set up the phase of oscillation at t = 0 one can make use of the translational symmetry t → t + t 0 of the signum-Gordon equation. Due to the periodicity of the solution the parameter t 0 can be chosen as t 0 = αγ. A sequence of generalized exact oscillons with different phases α is plotted in figure 8a. The configurations α = 0 and α = 1 differ exclusively by a spatial translation which implies that they have equal phases. In figure 8b we plot a worldsheet of the generalized exact oscillon in the laboratory reference frame in which it moves with the velocity u = +V . The endpoints of the oscillon are marked bỹ x(α) (the left one) and x(α) (the right one). It is pretty clear from this diagram that the length of the oscillon in the laboratory reference frame is a periodic function of the time when v = 0 and it is equal to 1/γ (for λ = 1) when v = 0. In order to get an initial configuration for the scattering process we translate the oscillon to the left by a distance x 0 (α). The second oscillon is obtained by the mirror reflection of the first one. Hence the initial configurations that are subject of considerations in this paper are those given by (3.3) with t 0 = αγ and x 0 = x 0 (α). In figure 9 we plot the initial profiles of the field Ψ (s) (x) with different phases.

Determination of the endpoints x(α) andx(α) of a oscillon
We start with the determination of the function x 0 (α) which describes the position of theright endpoint of the oscillon. To obtain it is sufficient to restrict considerations to v ≥ 0 because any non-travelling oscillon satisfies relation φ(t, x; −v) = φ(t, 1 − x). So we shall consider here the boosts in two directions, namely u = ±V , where 0 ≤ V < 1. In all formulae containing "±" the upper sign corresponds to u = +V and the lower one The simultaneity line t = αγ in the laboratory reference frame S is described by the equation in the reference frame of the oscillon S . In coordinates in S the right border of the oscillon is a worldline given by where τ (t ) is a function defined in (2.7). Eliminating x from (3.5) and (3.6) we find that JHEP01(2020)006 where is a straight line and τ (t ) is a saw-shape function plotted in figure 10. The plot shows two cases which differ by the sign of the boost velocity u.
Note that due to the periodicity of the oscillon the parameter α is restricted to the inteval 0 ≤ α < 1. This can be seen from figure 8b (case u = +V ), where the simultaneity lines t = 0 and t = γ, parallel to the axis x, cross the left border of the oscillon at (t , x ) = (0, 0) and (t , x ) = (1, 0) and therefore two oscillons seen in the laboratory frame S at the instants of time t 0 = 0 and t 0 = γ have the same shape. The same is true for u = −V case.
Two straight lines y(t , α) with α = 0 and α = 1, plotted in figure 10a, cross the axis t at t = −V and t = 1 − V . Similarly, in the case u = −V plotted figure 10b, they cross the axis t at t = V for α = 0 and t = 1 + V for α = 1. In order to cover the whole range of velocities V ∈ [0, 1) we shall consider the equation (3.7) on two intervals; namely, on the interval −1 < t < 1 for u = +V and on the interval 0 < t < 2 for u = −V . The saw-shape function τ (t ) is given by the pieces of straight lines

JHEP01(2020)006
There are two different cases dependent on the absolute value of the boost velocity V = |u|. They are separated by the critical case for which the velocity has value The straight line y(t , α) crosses maxima of τ (t ) i.e. τ (t max ) = 1 2 at t max = − 1 2 , 1 2 for u = +V and at t max = 1 2 , 1 in the case u = −V . It crosses the minimum τ (t min ) = 0 at t min = 0 for u = +V and at t min = 1 for u = −V . The corresponding values of the parameter α are denoted by α (±) l for lower maximum, α (±) u for upper maximum and α (±) 0 for the minimum. The signs "±" correspond to u = ±V . They have the form The upper letter stands for the upper sign (+) and the lower one for (−).
When the boost velocity is less than the critical value, V < V c , then the straight line y(t , α = 0) crosses the saw-shape chain τ (t ) at a certain point given by t belonging to the interval B D and y(t , α = 1) crosses the chain τ (t ) at the point with t that belongs to F . Note that the intervals B D and D F are covered only partially by α. On the other hand, for V > V c (the case sketched in figure 10) the relevant intervals are A C , B D and C E . In the interval A C the parameter α belongs to 0 ≤ α ≤ α In this case the intervals A C and C E are covered only partially by α. The solution of the equation (3.7) is given by where the values of (a, b) are determined by the velocity V and the phase α. A position of the right endpoint of the oscillon, obtained from the spacetime interval, takes the value where x , introduced in (3.6), is a function of α given by x (α) = 1 + v at (α) + b . After some manipulations one gets

JHEP01(2020)006
This results shows that x 0 is a linear function of α in the intervals in which (a, b) remain constant functions of α.
Next we determinate the functionx 0 (α) which describes the position of the left endpoint of the oscillon. The left border of the oscillon is described by the worldline x = vτ (t ) in the rest frame of the oscillon. In order to get the function t (α) one has to solve the equation Taking τ (t ) = a t + b we find that t (α) is given by The coefficients (a , b ) correspond to 0 ≤ t ≤ 1 2 and 1 2 ≤ t < 1 and they are given by We denote the solutions of the equation c . They are given by . Finally, taking similar steps for the case of the 'right endpoint' we find

Remarks on the initial configurations
Having determined the expressions for x 0 (α) andx 0 (α) we can now construct arbitrary initial configurations containg generalized exact oscillons which begin to collide at t = 0. In the case of oscillons with |u| < v an additional caution is necessary. In figure 11 we present the worldsheets of two generalized exact oscillons which form a symmetric configuration. The initial field configuration for the scattering process is given by  t (x)) cannot be obtained by approaching two travelling oscillons. A requirement that worldsheets of colliding oscillons have no intersections for t < 0 excludes such field configurations.
Although we have paid the main attention only to the symmetric (antisymmetric) initial configurations it is quite clear that nonsymmetric configurations with vanishing JHEP01(2020)006 total momentum can be obtained by taking colliding oscillon with different phases. An example of two oscillons with phases α L and α R is given by and shown in figure 12a. This field configuration is obtained by shifting one of the oscillons by distance x 0 (α L ) (placing it to the left of x = 0) and the second one by x 0 (α R ) and reflection x → −x (placing it to the right of x = 0). Another possibility is shown in figure 12b. The first oscillon, with u = −V , is shifted byx 0 (α R ) (which places it to the right of x = 0) and the second oscillon, with u = +V , is shifted by x 0 (α L ) (which places it to the left of x = 0). In this case It is quite clear from diagrams in figure 12 that there exist a set of phases for which the worldsheets of oscillons overlap for t < 0.

Antisymmetric configurations
In this section we present some numerical results for the scattering of two oscillons which form an antisymmetric initial configuration of the signum-Gordon field. We have chosen antisymmetric initial data to be discussed here first because the result of the scattering of oscillons in such a case is not as complex as that for symmetric configurations. The fact that initial configuration is antisymmetric implies that  in regions x < 0 and x > 0 effectively splits into two independent problems containing an initial oscillon and the boundary condition (3.13).
In our numerical study we have evolved an antisymmetric configuration without assuming the condition χ = 0 at x = 0. Looking at results presented in figure 13 we clearly see that the condition (3.13) is satisfied. This effect manifests itself in the presence of vertical white segments at diagrams which are located at x = 0. Another important observation is the absence of radiation in the central region of the diagrams independently on the value of the initial speed of oscillons. The sources of radiation generated in this process are irregular borders of two outgoing oscillons, see figure 13a.
One can note that irregularities of the border are more likely to appear for small velocities V of colliding oscillons than for the larger ones. Moreover, in spite of being irregular the outgoing oscillons radiate significantly less than it would be expected for strongly perturbed oscillons. In fact, the outgoing oscillons with irregular borders belong to a more general class of oscillons. We discuss this subject in more deatil in section 3.5.
Another interesting question is the dependence of the scattering process on the phase α that characterizes the initial configuration. In figure 14 we present the results of the scattering at V = 0.8 for four values of the phase. The figures show that the evolution of the initial configurations does not depend much on the value of the parameter α.
Another parameter that the incoming oscillons depend on is the speed v of the oscillon border in its rest frame. Our numerical studies have shows that irregularities of borders of two outgoing oscillons increase with increasing of the value of parameter v which determines the swaying motion of the incoming oscillons. In figure 15 we present the results of the scattering of an anti-symmetric configuration of oscillons with speed V = 0.8, phase α = 0 and speeds of the border v = 0.47 and v = 0.82. The figures show that as the borders of outgoing oscillons became more and more irregular the number and size of the oscillons emitted from these borders increases.
It is quite notable that the radiation generated in the process of the scattering of antisymmetric configurations is emitted only from the border irregularities of the outgoing JHEP01(2020)006               A typical situation with v = 0 is shown in figure 16. Two incoming exact oscillons form the initial configuration at t = 0. This configuration is adopted as an input for our numerical simulations.
We have performed many numerical simulations for distinct velocities and phases of scattered oscillons. In figure 17 we plot the scattering process of two exact oscillons with JHEP01(2020)006 v = 0 and α = 0 for six values of boost velocity V . All initial configurations contain two exact oscillons whose supports touch each other at t = 0. Numerics shows that the result of collision strongly depends on the velocity V of colliding oscillons. All six diagrams contain two main oscillon-like objects (they perhaps can be identified with perturbed oscillons) that emerge shortly after collision. They move with the velocity which is almost equal to velocity of colliding oscillons. The emerging oscillon-like objects obtained in the process of the scattering are significantly less regular when initial velocities of oscillons are small. For higher velocities the main outgoing oscillons are quite regular. Such field configuration are very close to the exact (generalized) oscillons. In the rest of the paper we call them quasi-oscillons. A fundamental difference between diagrams is visible in their central region where radiation appears in the form of jets. Looking more carefully at the radiation we see that it contains certain structures that strongly resemble oscillons. They look a bit like perturbed oscillons. Note also that the emerging oscillons obtained for small velocities V emit oscillon-like objects directly from their irregular border. A presence of radiation in the scattering process reflects the non-integrable character of the signum-Gordon model.
In figure 18 we present the numerical results of scattering of two oscillons with V = 0.93, v = 0 for four initial phases α = 0, α = 0.25, α = 0.84 and α = 0.93. The figure demonstrates that the initial phase of the colliding oscillons is indeed a relevant scattering parameter. The form of the jets in each subfigure of figure 18 is clearly different. This difference clearly shows that the result of the scattering process is very sensitive to the value of the phase α.
Looking at figures 17 and 18 we see that in the first stage of the scattering the interacting oscillons exist on support that shrinks from its initial size 2L to a certain minimal size 2L min where L = √ 1 − V 2 is the size of the oscillon in the laboratory reference frame. This process takes some time t s . For t > t s we observe emerging of two main oscillons and the appearance of waves of energy identified with radiation. In order to evaluate the time t s we assume that the left (right) border of the oscillon, which initially moves with velocity u = V (u = −V ), moves freely with the velocity that the oscillon has in the laboratory reference frame until it hits the future light cone of the event (0, 0); see figure 16. This assumption allows to determine the event P s with the coordinates (t s , x s ) in the laboratory JHEP01(2020)006 14) The formula (3.14) is valid exclusively for v = 0. Due to the symmetry of the initial configuration we see that the minimal support size of the scattering oscillons can be estimated by the expression ∆x min = 2 1−V 1+V . In figure 19 we show numerical data (dots) and analytical curves (solid lines) representing the characteristic time of the scattering t s and the minimum size of the oscillon configuration ∆x min at t s .     and ∆x min = 2t s .

High velocities -formation of shock waves
Next we discuss in more detail some of our numerical results. First we note that for small velocities the numerical solution is very irregular. Looking at figures 17a-17c we clearly see a formation of a strongly perturbed oscillon centred at x = 0. This perturbed oscillon is certainly unstable and it radiates out smaller oscillon-like objects. The situation changes for scatterings at high velocities as then the numerical solution has a more regular pattern.
In particular, when velocity of incident exact oscillons is close to unity we observe another JHEP01(2020)006 interesting solution. An example of such a solution is presented in figure 17e. The solution presented there was obtained for the scattering of two exact oscillons with initial phases α = 0, speeds V = 0.93 and no swaying motion i.e. for v = 0. A very characteristic feature of this numerical solution is a presence of regular waves that are localized in a diamond-like shape region on the Minkowski diagram. Such waves emerge shortly after the collision and eventually decay into a sequence of oscillon-like structures. It turns out that the nature of these waves can be understood in terms of the so-called shock waves that are exact solutions of the signum-Gordon model reported in [34]. A shock wave is a particular solution of the signum-Gordon model with two discontinuities that propagate with the speed of light. Here, however, we do not observe such wavefronts due to the presence of two oscillons that move with subluminal velocities. The collapse of the wave suggests that our numerical solution is only an approximation to the exact shock wave which would exist for arbitrary times t > 0.
We can check the hypothesis as to the nature of our wave solution comparing its zeros to the zeros of the analytical shock wave. According to [34] a shock wave solution belongs to the class of the signum-Gordon solutions described by φ(t, x) = θ(−z)W (z), where z = 1 4 (x 2 − t 2 ). The function W (z) obeys the ordinary equation zW (z) + W (z) = sgn(W (z)) and it consists of infinitely many partial solutions W k (z), k ∈ Z matched at points z = −a k . Each partial solution satisfies the equation zW k (z) + W k (z) = (−1) k and the conditions W k (−a k ) = 0 = W k+1 (−a k ) and W k (−a k ) = W k+1 (−a k ). So, it can be written in the form Note that α 1 = 1. Furthermore, it follows from (3.16) that α k+1 is determined by a k and b k . Solving numerically the second equation of (3.16) one gets y k+1 and so a k+1 and b k+1 can be determined too. Thus we note that the zeros of the field φ(t, x) are localized on the hyperbolas Outside of this region they break and produce oscillon-like structures.

Vanishing of the radiation
In our numerical studies we have also spotted a very interesting case. Namely, some symmetric initial configurations evolve in such a way that the resulting field contains only the main quasi-oscillons i.e. the amount of radiation generated in this process is virtually insignificant. This phenomenon was observed in the high speed range (V approximately above 0.7) i.e. when two main outgoing quasi-oscillons had very regular form. We have also found that, for a given velocity V , there are two phases α for which the radiation is absent. Figure 22 shows two examples of the scattering processes containing initial JHEP01(2020)006    . In particular, this shows that the incoming zeros x (0) (t), shown in figure 23a, and x (1) (t), shown in figure 23b, lie on the same straight line. Moreover, it can also be checked that both initial configurations characterized by phases α and α + 1 2 differ only by the sign i.e. ψ → −ψ when α → α + 1 2 , see the insertion plots in figure 23b.
The case presented in figure 22 is not unique. Our numerical results suggest that in the regime of high velocities V one can fine-tune the initial parameters V and α (and also v, see the next section) in such a way that the outgoing oscillons are not accompanied by almost any radiation. The absence of radiation demonstrates that the outgoing quasi-oscillon has virtually the same energy as the incoming exact oscillon. In this case the numerical quasioscillon is very similar to the exact one. In the most common situations there is certain difference of energies and this difference is explained by the release of a very large number of small oscillon-like structures. In figure 24 we plot the energy radiated out by the system as a fraction of its initial energy. Two dark regions in the upper part of the diagram correspond to the choices of the initial parameters (α, V ) such that this fraction is very close to zero. Looking at this picture we see that the initial configurations that minimize the radiation correspond to parameters α and V that lie, approximately, on straight lines. For V < 0.7 the outgoing oscillons are quite irregular; see for instance figure 17a. Usually, some radiation can be seen in the vicinity of outgoing irregular perturbed oscillons and this radiation is emitted from these irregularities. In consequence, the determination of the energy of outgoing oscillons is less reliable for small velocities.

Oscillons with v = 0
In this section we present some of our results on the scattering of generalized exact oscillons i.e. oscillons which depend on v -the parameter controlling the swaying motion of the JHEP01(2020)006 oscillon endpoints. In figure 25 we show the wordsheets of two such oscillons. The configuration at t = 0 was taken as the initial data for our numerical simulation. In similarity to the case v = 0 we can find the expressions for the characteristic time of the scattering t s and the minimum support size ∆x min = 2|x s | by solving the equations x = −L(α) + wt and x = −t. We have found where L(α) ≡ x 0 (α) − x 0 (α) is the size of the oscillon at t = 0. x 0 (α) is given by (3.10) and x 0 (α) by (3.11). An interesting question now arises of how the replacement of v = 0 by v = 0 modifies the results of the scattering processes. For instance, we can take a v = 0 generalization of the symmetric initial configuration with V = 0.93 and α = 0.414. In figure 26a we present the result of the scattering for v = 0.02. We see that even such small value of the parameter v is sufficient for the appearance of shock waves which further transform into a cascade of oscillons. This demonstrates that the scattering process is quite sensitive to the value of the parameter v. In order to minimize this radiation one can adjust properly the parameter α. We have found that in this specific case the radiation vanishes for α = 0.420 and α = 0.918, see figure 26b and 26d. For higher values of v there was much more radiation emitted during the process of scattering.
In figure 27 and figure 28 we present the cases of v = 0.2 and v = 0.7. In both cases we have found two values of the phase α that minimize the radiation.
Looking carefully at the initial configurations we see that there is a significant difference between the case v = 0 and v = 0. In figure 29 we present the initial field configurations that minimize the radiation. Comparing configuration Ψ (s) (x; v, V, α + ∆α), where ∆α is taken from the numerical simulations, with Ψ (s) (x; v, V, α) we see that the configuration with α+∆α is not equal to the negative of the configuration with α. The difference between JHEP01(2020)006             In figure 30 we have plotted the fraction of the initial energy carried by the radiation. The figure was produced for v = 0.45. The dark regions, corresponding to very low values of the radiated energy, are less regular when compared with figure 24 for the case of vanishing v.

Non-symmetric configurations
The simplest non-symmetric configurations are described by oscillons that differ only by their phases. In order to probe this class of scatterings, we use initial conditions given      by (3.12) taken at time t = 0. In figures 31a-31d we plot four cases of the time evolution for these processes. Interestingly, there are all non-symmetric configurations for which no radiation is present. Figures 31b and 31d show two of these cases.
In order to have a clearer picture of the amount of input energy converted into the radiation, we present a density plot much like the ones presented in figures 24 and 30, except that in this case we have fixed the values V = 0.93 and v = 0, and we vary the parameters α L and α R (now independent of each other). This plot is shown in figure 32.
The case α L = α R corresponds tothe symmetric configurations of the initial conditions, and is marked in the plot as a black dashed line. This line passes through two minima in the emitted radiation, as is expected according from the results discussed in section 3.3. Along this line, the plot has a coincident set of values with the plot in figure 24 along the line given by V = 0.93. Also, as mentioned in section 3.3.3, for null swaying speed of the oscillon endpoints (v = 0), a shift of 1 2 in the phase of a given oscillon produces the same oscillon with a sign change in the value of its field and its time derivative. For this reason, the relation α R = α L ± 1 2 corresponds to the anti-symmetric initial field configurations. This relation is marked in the figure as the two dot-dashed white lines. In agreement with the results from section 3.2, these lines lie on top of the dark regions that correspond to no-radiation zones.
The plot can be seen as periodic both in α L and in α R (although the phase is originally defined as a value between 0 and 1, the field configuration is precisely the same for both these values), leaving it with a toroidal topology. So, the two black strips marked by the white dot-dashed lines form a belt around the torus and become a single continuous region.
We note, still, a second strip of no-radiation forming a 90 degree angle with these lines. One could expect the eventual vanishing of radiation in non-symmetric initial scattering configurations e.g. figures 31b and 31d. Yet, the regularity of this region (it is a straight line) is quite remarkable and requires further considerations.
We present a similar plot with α L fixed in which we varied α R and V . This plot is presented in figure 33 and, along the line given by V = 0.93, its values coincide with those of figure 32 for α L = 0.
We note that the strip along α R = 1 2 corresponds to the anti-symmetric configurations and so it shows little or no radiation around it. Below the value V 0.7 there appears a large irregular region devoid of radiation. From the standpoint of numerical stability of our 3.0% 4.0% 5.0% Figure 33. Fraction of the total energy of the initial configuration carried out by the radiation as a function of the input speed V and phase of right input oscillon α R while holding fixed the right input oscillon's phase α R = 0.
methods for the measurement of radiation in regions of low V , as mentioned in section 3.6, this region could well be a numerical artefact. In a brief investigation of this hypothesis, though, we have found this void to be an accurate description and the region, indeed, represents a zone of no-radiation interactions for very low energy input (low boosts). So the value V 0.7 is critical in the sense that, below it, the entire region for this particular set of values in the parameter space, generates initial conditions to the scattering that produces no radiation, and the region, itself, does not seem to have a very well defined shape.
One possible explanation of this void is that there is, indeed, some generation of the radiation in the scattering process but this radiation happens to be absorbed by the outgoing oscillons, which in turn makes them become more perturbed. This hypothesis can be checked by considering the energy of each individual outgoing oscillon in this region. In figures 34a and 34b we plot the energy balance of each outgoing oscillon compared to its initial energy. This balance is calculated by where E L is the energy of the left outgoing oscillon (which entered the interaction region from the right), E R is the energy of the right outgoing oscillon and E the total input energy. Note that the plot in figure 33 shows the value of the amount of energy lost to the radiation, which is given by so that the total radiation is just the sum of the energy balances of individual scattered oscillons.   From our figures (figures 34a and 34b) we see that each oscillon looses/gains a considerable amount of energy in the process. This amount, in some casees, is close to up to 15% of each oscillon's incident energy, yet the total radiation generated is no larger than about 5.5% of the total incident energy. In figure 35 we present the case for V = 0.5, α L = 0, α R = 0.642, v L = 0 = v R , which corresponds to a configuration that produces no radiation (as it is located within the void region) and, at the same time, has a relatively large energy transfer between the two interacting oscillons (the energy balance plots of these regions show a large energy change of both oscillons). Note that the left input oscillon (outgoing towards the right) is larger than the one outgoing to the left. Most such configurations seem to reproduce this behaviour, and the channel through which this energy is transferred, in all cases, seems to be related to the scale of the outgoing oscillons.

Oscillons with accelerating borders
Looking at figure 17 and figure 13 we note that the outgoing oscillons are significantly different from the exact generalized oscillon with uniformly moving endpoints. The main difference is in the form of worldlines describing the borders of the oscillons which take the form of continuous curves rather than zig-zag piecewise straight lines. Moreover, the curves are surprisingly regular in shape and this suggests that the outgoing objects are JHEP01(2020)006 not just simple perturbations of the generalized exact oscillons. Recently, in [26] further generalization of the signum-Gordon oscillins has been proposed. This generalization leads to the emergence of oscillons with borders that are described by arbitrary time-like curves.
Here we present a construction of oscillons with curvilinear borders and produce some plots of such oscillons after applying to them Lorentz boosts.

General properties
For an oscillon with period T the complete solution can be constructed from the restriction of the solution ϕ(t, x) to the interval 0 ≤ t ≤ T /2. Similarly to the oscillons already known, we can get periodic and localized solutions imposing the conditions where f (x) is a continuous function such that f (x) ≤ 0 for all 0 ≤ x ≤ T and f (0) = f (T ) = 0. As a consequence, we assume that the solution ϕ(t, x) is negative for 0 ≤ t ≤ T /2. The oscillon solution is localized in the sense that it is nonzero only in the region between two time-like curves γ L and γ R -the borders of the oscillon, see figure 36. The right border is a displacement by T of the left curve, so that the size of the oscillon remains constant. For t ∈ [0, T /2] the borders move to the right by ∆ and, for t ∈ [T /2, T ], the borders move in the opposite direction -returning to the original position.
The conditions (3.20) and (3.21) are satisfied if we take the following ansatz for the non-zero part of the solution: where F (x) satisfies: For less general oscillons, we can demand that the solutions approach zero smoothly at the left border γ L i.e. the partial solution (3.22) for −t < x < t gives where points belonging to the curve γ L have coordinates related by x = x(t). Expressing condition (3.25) in light-cone coordinates where we have used the fact that (y + , y − ) are related and dy + dy − = − 1 2 dx dt . Since γ L is a time-like curve dy + dy − < 1 2 and so condition (3.27) is equivalent to ∂ + ϕ(y + , y − )| γ L = 0. This leads to F (y + ) = − 1 4 (y + − y − ) for the first expression in (3.22). Taking into account expression (3.24) we finally get f (y + ) = − 1 2 (y + − y − ), (3.28) where y − = g(y + ) is a function of y + representing the left border of the oscillon (worldline γ L ). Similarly, demanding that ∂ x ϕ(t, x)| γ R = 0 at the right border of the oscillon γ R we get a condition which, when written in terms of y − , takes the form ∂ − ϕ(y + , y − )| γ R = 0. Then, from the last expression in (3.22) we get Here y + = h(y − ) is a function of y − representing the right border of the oscillon (worldline γ R ). Note that for points on γ L we have 0 ≤ y + ≤ T /2 + ∆ and for points on γ R we have Using this formalism, all oscillon solutions limited by a pair of identical time-like curves can be constructed by plugging the a priori given trajectories of the border in (3.28) and (3.29). The trajectories should have a form that one can to describe them explicitly as y − = g(y + ) for γ L and y + = h(y − ) for γ R . The only remaining problem is to integrate the resulting expressions to get F (x) -and consequently ϕ(t, x).

Example
The first interesting example we have found is an oscillon with borders having a constant acceleration a in the instantaneous rest frame of the border. In what follows we will use units in which c = 1. In the reference frame of the oscillon, in which the border has acceleration γ −3 a, the trajectory describing the motion of such borders takes the form where v 0 is the velocity of the border at t = 0, γ 0 = (1−v 2 0 ) −1/2 and x 0 = 0 (x 0 = T ) for γ L (γ R ). Note that if the oscillon has no extra motion then its rest frame is just the laboratory reference frame. Using the light-cone coordinates y ± = x ± t we put the expression (3.30) into the form y − = g(y + ), where x 0 = 0, and y + = h(y − ) and where x 0 = T . Then plugging these expressions into (3.28) and (3.29) we get our explicit expression for f (x): Once we know f (x), it is possible to integrate it and get partial solutions ϕ k (t, x) with k ∈ {C, L 1 , L 2 , L 3 , R 1 , R 2 , R 3 }, each one valid in a specific subset of the region between γ L and γ R . Such solutions are given by:

JHEP01(2020)006
where AB = a −2 . Similarly to the oscillons previously presented and discussed we can relate the solutions ϕ R i (t, x) to the solutions ϕ L i (t, x) through the transformations: a → −a. (3.40) Note that the last two transformations are equivalent to the transformations: Once again we can write the complete solution with the help of the step functions. We have where ∆ = x(T /2) − x 0 and so is given by The periodicity of the solution can now be taken into account by involving the generalized forms of the functions τ (z) and σ(z) for an arbitrary period T : σ(z) = sgn sin 2πz T . (3.51)
Comparing the generalized exact oscilons presented in figure 37 with numerical solutions presented in figure 17 and figure 13 we conclude that there is certain similarity between outgoing oscillon-like object produced in the scattering process and the exact oscillons with non-uniformly moving endpoints. This suggests that the process of the scattering of two exact oscillons can lead to the production of field configurations corresponding to perturbed generalized oscillons with non-uniformly moving endpoints. In some cases these perturbations can be almost zero and when this happens one can call the outgoing oscillon-like field configurations quasi-oscillons with non uniformly moving endpoints. The transformation of oscillons from one class into the oscillons belonging to another class during the scattering process is an open problem which requires further investigations.

Fractal nature of the radiation
In this section we discuss the interesting possibility that the radiation generated in the process of the scattering of oscillons has a fractal-like nature. There are two facts that suggest this. Firstly, many numerical simulations show that the radiation of the signum-Gordon model is dominated by oscillating structures that look like travelling oscillons.
Such oscillon-like structures are generated as radiation during the scattering of oscillonlike objects. Alternatively, they are emitted from strongly perturbed oscillons. In fact, a production of small-size oscillons during the evolution of perturbed oscillons has been conjectured in ref. [24]. Our work presented in this paper and also in [37] suggests that this conjecture is true. Secondly, the oscillon-like field configurations exist at arbitrarily small scales. Although the numerical approach does not allow for the arbitrary good resolution we know that the existence of exact oscillons with any size is guaranteed by the dilation symmetry of the signum-Gordon equation (1.2). This symmetry implies that for any real number λ > 0 and any solution of the signum-Gordon equation φ (1) (t, x) the function is also a solution of (1.2). Looking at energy of solutions we see that it scales according to     process is a source of new smaller oscillon-like objects. Certainly, we do not expect that oscillon-like objects seen in our numerical simulations are exact oscillons. On the other hand, many of them are surprisingly regular and sufficiently stable. They all possess characteristics necessary to be called quasi-oscillons. Less regular oscillon-like structures "decay" emitting smaller and more regular oscillating objects. The emission of smaller oscillons is a physical mechanism allowing strongly perturbed oscillons to get rid of a surplus of their energy. Summarizing, we can say that interaction between individual oscillon-like objects (constituents of radiation) produces more and more such objects during their evolution. Finally, we note that the relation (3.53) is quite general and it allows us to take as a solution φ (1) (t, x) not only a single oscillon but the whole diagram. Certainly, there is no qualitative difference between scattering processes involving two oscillons with λ = 1 and a scattering of smaller oscillons with λ 1. In principle, the whole diagram (like figure 27a JHEP01(2020)006   and 27c) could repeat itself at any length scales. Such repetitions of structures involving oscillons at all length scales in the spacetime diagram suggest a possible fractal-like nature of the radiation. This statement still has a status of a conjecture and it certainly deserves further investigation. Below we present only some preliminary results of a numerical study which reinforces this idea. In order to check our hypothesis we have perform high resolution simulations of the scattering processes and then looked at the spacetime diagrams representing the result. In figure 38a we plot the two main outcoming oscillons and the radiation in the central region between these oscillons. Looking in more detail at the region inside the rectangle in figure 38a which we replot in figure 38b we see that there exist a huge number of smaller oscillons invisible in the previous picture. Choosing another rectangular region of figure 38b which we replot in figure 38c we see that, again, it contains many oscillating structures. This result supports our idea of the fractal-like nature of the radiation of the signum-Gordon model.

Conclusions
In this paper we have reported our results on the scattering of compact oscillons in the signum-Gordon model in one spatial dimension. We have looked at two qualitatively distinct initial configurations -symmetric and anti-symmetric one. In both cases the initial configurations consisted of exact compact solutions. Due to the compactness of oscillons JHEP01(2020)006 there was no problem with their overlapping at t = 0. In fact we also evolved oscillons whose supports touched each other but did not overlap at t = 0. A time dependence of the shape of oscillons was responsible for the existence of an additional scattering parameter which we called the phase of the oscillon. This phase was a important quantity and the properties of the scattering process depended very strongly on it.
Looking at the results of the scattering of oscillons we have found that there was a significant qualitative difference between symmetric and anti-symmetric initial configurations. The emission of radiation for anti-symmetric configurations was restricted to situations where outgoing oscillons had irregular borders. Such irregular borders act as sources of radiation which has the form of showers of smaller oscillons sent out from the borders. The central region of the Minkowski diagrams just after emergence of outgoing oscillons was free of radiation. On the other hand, symmetric configurations produced much more radiation than the anti-symmetric ones. In this case the radiation was emitted mainly in the central part of the Minkowski diagram where structures similar to exact shock wave solutions of the signum-Gordon model were formed. These waves were not stable and, eventually, they decaysed into cascades of oscillons. We have spotted that there were special values of the phases of colliding oscillons for which there was almost no radiation. We suspect that this fact was associated with the absence of the shockwave-like structures between outgoing oscillons. The relation between the collapse (decay) of the shockwave-like solution and the appearance of a cascade of oscillons has been found to be a very interesting subject and it requires more thorough analysis than we could carry out in the present paper. We hope to report more on this subject in near future.
Comparing incoming oscillons with outgoing ones we have spotted that, in general, the later ones belong to a wider class of oscillons. This class is characterized by a nonuniform motion of the border of the oscillon in its own rest frame. In our numerical study many of the outgoing oscillons had borders described by a segment of the worldline curve whereas for incoming oscillons these borders were segments of straight lines. Thus the collision transformed the compact oscillons of a very special class into more general compact oscillons.
We have also looked at the radiation of the signum-Gordon model and have found that it possessed what looked like a self-similar structure. Since the model has the scalling symmetry one can show that an exact compact oscillon can have arbitrarily small support and energy. Our numerical studies have shown that small quasi-oscillons were emitted from perturbed oscillons or appeared in the scattering processes of two oscillons-like structures. Since, in general, they were also perturbed the process of emission repeated itself (in principle infinitely many times). This mechanism of emission of oscillons from perturbed objects and the fact that oscillons existed at arbitrarily small scales suggests to us the emergence of dynamical fractals i.e. of the fractal structures in the spacetime diagrams.
Final remarks.
1. Our investigations of the scattering processes have been based primarily on the numerical integration of the signum-Gordon equation. The complexity of this process excludes any analytical approach to this problem. We have made many attempts to JHEP01(2020)006 calculate analytically the evolution of the initial profile containing two exact oscillons. Unfortunately, even before the emergence of main oscillons we have encountered technical difficulties in the construction of partial solutions. Moreover, the number and localization of any partial solution depends on the initial data. In contradiction to the standard analytical non-linear models the perturbative approach cannot be used in the case of signum-Gordon model because of the non analytical character of the potential V = |φ| at φ = 0. Hence, the small perturbations of the vacuum solutions are always nonlinear.
2. One can get some analytical results considering the decay of shock waves in a cascade of oscillons. This subject has been recently studied and the results have been reported in [40].
3. Of course, we can also think of the comparisons of our results with those obtained in other models, such as the Sine-Gordon model of λφ 4 model. However, these models are basically very different as they do not possess compact solitons. Of course, their solitons are exponentially localised and some studies of such solitons have also been performed in much detail. The most comparable studies involved looking at the properties of Sine-Gordon kinks on scattering on various obstructions (potential wells or barriers) and the effects of the obstructions on the properties of the basic solitonic structures. The obstructions generated the emission of kink-antikink pairs, either in the form of breathers or as invidual pairs. And for perturbed models, which still had solitonic solutions, one had emission of long lived breathers (basically oscillons) or annihilations. An interested reader can look at papers [41][42][43][44] and references therein.

A Comments on the numerics
The numeric results presented in this paper have been generated by the use of the standard 4 th -order Runge-Kutta method, integrating the system via the discrete timesteps ∆t. The second-order equation in time has been decomposed into a coupled system of two first-order equations φ(x, t i ) = φ(x, t i−1 ) + ∆t ψ(x, t i−1 ) (A.1) where t i = n ∆t (n = 1, 2, · · · ) and ψ(x, t) = ∂φ ∂t . The spatial dimension was made discrete over N sites of width ∆x, so it had width L = N ∆x.
In most of our simulations we have used a spatial resolution of N = 2 15 (for L = 6 this corresponds to ∆x 1.8 × 10 −4 ). There were two exceptions to this. The first one corresponded to the case of the fractal, which in order to generate and capture the small scale details we had to us N = 2 20 (in this case, the simulation space length was L = 12, leading to ∆x 1.14 × 10 −5 ). The second exception corresponded to the generation of the diagrams of describing the percentage of energy lost to radiation (e.g. figures 24, 30) and JHEP01(2020)006 the balance of energy after the interaction (figures 34a and 34b). Since the value of each pixel of these images is computed based on one entire simulation, in order to speed up the computations (in the case of figure 32, we have used 350 2 = 122500 simulations) we had performed lower resolution simulations, with N = 2 12 .
The timestep value, in all our simulations, was given by the relation a = ∆t ∆x and we have used in all our simulations a = 0.1.
Finally, we would like to point out to the reader the need of some caution while dealing numerically with the very small scale structures appearing within the fractals. This is further explained in section B.

B Caveats
We would like to add a few comments on the difficulties associated with the numerical integration of the signum-Gordon equation. The main difficulty has origin in the fact that the radiation contains perturbed oscillons of arbitrarily small sizes. Certainly, oscillons smaller than the size of numerical domain cannot be seen in the simulation. An obvious solution involves increasing the number of points in the grid. We have run many simulations changing the number of points and comparing the results. Such tests have shown that very small oscillons are very sensitive to the number of points. In some cases, changing the number of points by a factor of two resulted in the appearance or even disappearance of some tiny structures whereas bigger structures remained stable under this procedure. A similar problem was spotted in simulations of a special type of self-similar solutions with an infinite number of zeros on a finite segment. In that case the numerical solution and the analytical one diverged after a very short time (the numerical solution was very unstable), and the increase of the simulation resolution resulted in only very small in stability. On the other hand, our numerical simulations of the exact oscillons did not lead to visible instability within intervals of time corresponding to many oscillon periods. Also, the simulations of exact shock waves were very consistent with the analytical solutions. Thus, in the regions dominated by radiation (or the special self-similar solutions) the solutions of the model were found to be sensitive to the initial conditions. In this sense, the signum-Gordon model shares some properties with chaotic systems. This property is, for instance, one of the main difficulties in generating of high resolution fractals. However, although some results were quite sensitive to the details of the details of the numerical procedures, most of them were not and we strongly believe in their validity.