Lorenz Type Behaviors in the Dynamics of Laser Produced Plasma

: An innovative theoretical model is developed on the backbone of a classical Lorenz system. A mathematical representation of a di ﬀ erential Lorenz system is transposed into a fractal space and reduced to an integral form. In such a conjecture, the Lorenz variables will operate simultaneously on two manifolds, generating two transformation groups, one corresponding to the space coordinates transformation and another one to the scale resolution transformation. Since these groups are isomorphs various types isometries become functional. The Lorenz system was further adapted to describe the dynamics of ejected particles as a result of laser matter interaction in a fractal paradigm. The simulations were focused on the dynamics of charged particles, and showcase the presence of current oscillations, a heterogenous velocity distribution and multi-structuring at di ﬀ erent interaction scales. The theoretical predictions were compared with the experimental data acquired with noninvasive diagnostic techniques. The experimental data conﬁrm the multi-structure scenario and the oscillatory behavior predicted by the mathematical model.


Introduction
Laser ablation embodies a series of phenomena with a complex interconnection between them [1]. For the better understanding of this physical process there have been multiple diagnostics techniques implemented which showcased the particle removal process [2], particle dynamics after ejection [3], plasma formation and expansion [4][5][6], etc. Complementary, theoretical aspects of the process have been steadily showcased in the past 30 years [7][8][9]. The difficulty in developing theoretical models for a multi-physics process like laser ablation comes from the combination of different elements from laser physics, solid state physics and plasma physics. Thus, an appropriate theoretical approach should be able to transition different interaction scale (ns or fs for the laser beam, fs or ps for the laser matter interaction and microsecond for the plasma expansion) and still keep a connective link between them.
The development of theoretical models that can accurately describe the ablation process has been a main direction for the understanding the underline phenomena and their interdependencies. A robust model needs to contain as input variables, laser properties, the nature of the material, ejection of the material and formation of a transient plasma. This is a difficult task as the majority of the existent theoretical approaches focus on particular sequences of the laser ablation. Therefore, we found classical theoretical approach that focused on ultra-short laser interaction with dielectric material or metals [10], ns-laser ablation [11], thermal evaporation of material, particle ejection mechanisms and plasma formation or expansion. In the past decade there have been some proposals based on a fractal paradigm [12][13][14]. These attempts were focused on the behavior of the laser produced plasma, mainly on peculiar results like multiple structuring during expansion [6], particle oscillations [5] and temporal and spatial distribution of some plasma parameters [13]. The advantage of representing the laser ablation process in a multifractal space is given by the ease with which we can transition between different interaction scales and different sequences of the process. This was shown for a laser produced plasma on a single element [13] and multi-component target [12]. The core premise of the model is a hydrodynamic description of the plasma plume expansion, without an explicit involvement of the experimental parameters.
In this paper we develop a mathematical model starting from classic Lorenz system and we transpose it in nondifferential (fractal) representation. From such a perspective, Lorenz type variable will depend both on space coordinates and scale resolutions. Consequently, every variable will act as a limit of a family of functions which are non-differentiable for null scale resolutions and differentiable for non-null scale resolutions; every variable will operate simultaneously on two manifolds, one generated by the space coordinates' transformation group and another one generated by the scale resolution transformation group; these two groups are isomorphs, so that various types of isometries can be applied (embeddings, compactizations, etc.). The system will be solved in a fractal space. The final aim of this development is the understanding the dynamics of a laser produced plasma with a multi-component structure, generated on a complex target. The solution allows the simulation of particle velocity distribution across a wide range of scale resolutions; the charge particle spatial distribution, showcasing the formation of a space charge double layer at the interface between plasma structure; and charged current temporal trances at various scale resolutions. The theoretical predictions are compared with experimental data extracted from a laser produced plasma on a chalcopyrite target by means of Intensified Coupled Charged Device (ICCD) fast camera imaging and optical emission spectroscopy. The theoretical simulations are confirmed by the experimental data.

Route to Non-Differentiability
Non-linearity and chaoticity are fundamental attributes of a Lorenz type system. Usually, the classical models built to describe the Lorenz type dynamics, are developed on the presumption of differentiability and integrability at all scale resolutions. The successes of these approaches need to be understood sequentially, on domains in which the differentiability and integrability are still respected. The differential and integral mathematical procedures fail when we attempt to describe nondifferentiable dynamics of Lorenz type system.
In order to describe the non-differentiable dynamics of a Lorenz type system and attempt to remain tributary to the well-established differentiable and integral mathematical procedures, it is necessary to introduce the scale resolution both in the variable and differentiable equations describing such dynamics. All the variables are describing a Lorenz type system from a non-differentiable perspective, which depends on both space-time coordinates and the scale resolution. Consequently, instead of operating, for example, with a single variable described by a nondifferentiable function, we will operate only with the approximations of this function obtained by mediating it at various resolution scales. This special function will act as a limit of a family of functions which are non-differentiable for zero scale resolution and differentiable for non-zero scale resolutions.
This way of describing nondifferentiable dynamics of a Lorenz type system implies both the generation of a new geometry of movement and a new class of Lorenz type models (which we will refer to as fractal/multifractal Lorenz systems). According to these geometric studies for which the movement laws, invariants of the space time coordinates' transformation are integrated over the scale resolutions which are invariant of the scale resolution's transformation. If we further admit that these geometric structures are based on the concept of multifractality, then the holographic implementation of movement-seen through scale relativity theory with an arbitrary fractal dimension or through operational procedures in the description of physical systems [14]-becomes essential when describing the dynamics of non-differentiable Lorenz type systems. However, both ways presented above imply the definition of nondifferentiable dynamics of a Lorenz type system on continuous but nondifferentiable curves. Therefore: (i) Any variable used to describe the dynamics of a nondifferential Lorenz type system will be described through multifractal mathematical functions dependent on both the spatial and temporal coordinates, and on the scale resolution. (ii) The laws describing these dynamics are invariant with respect to the spatial coordinates and temporal transformation, and the scale resolution transformation. (iii) The constraints on the Lorenz type system dynamics, described through continuous and differentiable curves of a Euclidian space, are replaced by the dynamics of a system lacking any constraints, and being described by continuous and nondifferentiable curves in a multifractal space. (iv) Between any two points in the multifractal space there is an infinity of curves describing the dynamics of a systems (its geodesics). The indiscernibility between these curves is a natural property of multifractalization through stochasticization; meanwhile, their discernibility is the result of a selection process based on the principle of maximum informational energy [14]. From such a perspective, any Lorenz type system with dynamics described by continuous and differentiable curves has hidden dissipative information (lacks memory). Otherwise, Lorenz type systems described by continuous and nondifferentiable curves have explicit information (presents memory).

Scale Resolutions
Let us consider a multifractal function (representing any of the variables describing the dynamics of a Lorenz type system) f (u), defined in the closed interval u ∈ [a, b]. Let us also consider the set of values for the u variable: We will define f (u, µ) as the broken line connecting the points: Then f (u, µ) becomes a µ-scale approximation.
Let us now consider f (u, µ), the µ-scale approximation of the same multifractal function. Since f (u) is self-similar virtually everywhere, if µ and µ are considered small, the approximations f (u, µ) and f (u, µ) lead to the same result when multifractal phenomena are investigated. Comparing these two situations, there is an infinitesimal increase/decrease dµ of µ that would correspond to an increase/decrease dµ of µ, only in the case of scale contraction or dilatation. In that case, the following relationship is satisfied: Therefore, the rational for the scale µ + dµ and dµ needs to be constant. In these conditions we can consider the infinitesimal transformation of the scale as: by performing such a transformation, in the case of f (u, µ), it results in: Furthermore, if we use stop at the first approximation of the function we get: Meaning, Moreover, we note that, for a fixed arbitrary µ 0 , so that Equation (7) can be rewritten as: In the end we will obtain: which showcases the dilatation/contraction operator: Equation (11) also showcases the fact that the intrinsic variation of the scale resolution is not in µ but on ln µ µ 0 .

Non-Differentiable Lorenz Type Systems
Let us now consider the classic Lorenz system [15,16] described in a non-dimensional coordinate system by the differentiable equations: In Equation (12) the variables are obviously differentiable and can be integrated. According to the paradigm presented above, when we consider the non-differentiable and non-integrable variables, the system of Equation (12) becomes a multifractal Lorenz type system. For this to occur, the variable needs to be dependent both on the spatial and temporal coordinates, and on the scale resolution. Moreover, the affine parameter T which characterizes the trajectories in the space phase (X, Y and Z) needs to be dependent on the resolution scale. Using such a hypothesis, the variables X, Y and Z will operate simultaneously on two manifolds, one generated by the space-time coordinates' transformation group and another one generated by the scale resolution transformation group. Since these two group are isomorphs, or in an extremely restrictive case self-morphs, embedding type isometries, self-embedding isometries, compactization type isometries, etc., become functional.
A possible embedding scenario can be performed through scaling: with, In Equations (13) and (14), ε can functionally be dependent on only one fractal dimension, D F , following the relationship ε = ε(D F ). In that case, using Equation (12) we will describe the dynamics of a mono-fractal Lorenz type system. We can find another case of ε which can be dependent on a singularity spectrum f α following the relationship ε = ε( f α ). In that case, Equation (12) describes the dynamics of a multifractal Lorenz type system.
The scaling Equation (13) specifies the fact that, simultaneously with the contraction in the space phase attributed to a mono-fractal or multifractal dynamics of a Lorenz type system, the dilation of time takes place. As such, the classical Lorenz system (12) through (13) and (14) becomes a multifractal Lorenz type system described by the following system: either through ε = ε(D F ) or ε = ε( f α ).
In this system the dynamics variables X, Y and Z are multi fractal functions, as they depend not only on the space-time coordinates but also on the scale resolution. Therefore, instead of working with X, Y and Z variables, which are described by non-differentiable functions, we will operate only with approximations of these functions obtained by their averaging at various scale resolutions given by Equation (14). Any of the variables X, Y or Z will behave like limits of a family of functions, which are non-differentiable on null scale resolutions ( ε → 0 , and µ = µ 0 ) and differentiable on nonzero scale resolutions (ε 0). In these conditions the movement laws given by Equation (15), regardless of any coordinate transformations, can be integrated with the scale laws, regardless of any scale resolution transformations.
In such context, a further in-depth study of the systems dynamics needs to be done for the case in which the scale resolution is null, thus for the case in which Equation (15) defines a non-differential system. Let us consider that in Equation (15) the scale resolution in null, thus the restriction ε → 0 is satisfied. Then Equation (15) takes the following form:

Motion Integration
The system in Equation (16) presents some interesting features, such as the phase space volume associated to the dynamics of the system being conserved, since it satisfies the relationship: This means that through the restriction ε → 0 , describing the transition from a differentiable Lorenz type system to a non-differentiable Lorenz type system, the system evolves from a dissipative one to a non-dissipative one. The non-differentiable Lorenz type system admits two integrals of the motions. The first is: This relation is obtained by multiplying the first equation from (16) with Y and further substituting the term XY through the third equation of Equation (16) and finally integrating the result. The second movement's integral is: Relation (19) is achieved by multiplying the second equation from Equation (16) with Y, substituting the terms XY and XYZ with the third equation from Equation (16), and finally integrating the result. Let us further note that the dynamics introduced through a non-differential Lorenz type system are described analytically through elliptic functions. In order to showcase this, let us consider the square of Equation (16) using Equations (18) and (19), written as: It reveals the elliptic integral: dX with Admitting that the P(X) polynom can be written as: where u 2 1 and u 2 2 have the following expressions: The integral (21) becomes: dX Furthermore, by using the substitutions: the integral (24) takes the Legendre form: The solution to that integral is the Jacobi elliptic function sn, or in its original coordinates: From here, through Equation (18) we will get the expression for Z: while by using Equation (19) we can obtain the expression for Y: Now, the expressions for X, Y and Z can be simplified even more if we take into account the relationship: cn is a Jacobi elliptical function, s is the modulus of the elliptical functions cn and sn and s is the complementary modulus s 2 = 1 − s 2 a and ∆T = T − T 0 .
In the following we will attempt to implement the solution gained by Lorenz system to a well-known system covering the ejection and dynamics of particles.

Plasma Modelling
The fractal representation of phenomena like laser produced plasma has a simplified elegance, as it confines all the complex intricate behaviors in a handful of parameters. This is a great advantage when attempting to simulate a wide range of behaviors for multiple external conditions, spatial and temporal coordinates, etc. However, the interpretation of the obtained results requires a direct correspondence between the compact fractal parameters and real, measurable plasma parameters. Given the description of our three normalized functions (X, Y and Z) they define complex functions that will be used further to describe formation and evolution of the ejected particles at different scale resolutions. Therefore, in our paradigm: X will define the particle distributions, Y will define the charged particle current as a function and Z will define the charge density fluctuations induced by some unbalances of the transient electrical field during expansion, and can be associated with the ratio between the kinetic and thermal energy of the ejected particles. Let us note that we will operate with the previous functions and parameters as normalized quantities. Since the Laser Produced Plasmas (LPP) dynamics are described through continuous and non-differentiable curves (fractal curves with different degrees of fractality), we had to operate in this situation with what we call the singularity spectrum. If in the system the dynamics are characterized by only one fractal dimension then we are dealing with mono fractal dynamics. If for a system (which is the case for laser produced plasmas) we deal with simultaneous dynamics with various fractal dimensions, then the role of the singularity spectra is not only to showcase the variation domain of these dimensions but also the characteristic classes associated by defining the strange attractors.
In Figure 1a,b we have represented the particle velocity distribution at various scale resolutions and at various moments in time. We can see from the two representations that the evolution of the plasma in time, implies the presence of families of particles defined by specific scale resolutions and velocity distributions. With an increase of the scale resolution we see that the fractal image of the laser produced plasma showcases the appearance of a one particle distribution centered on a relatively high velocity, at low resolutions scale and a consistent increase number of distributions centered around lower velocities. The presence of multiple structures in laser produced plasma has been previously attested and reported as a direct result of the multiple ejection mechanism and plasma interaction with the background gas. We note that the presence of multiple ejection mechanisms [17][18][19][20] (Coulomb explosion, explosive boiling, phase explosion, etc.) will lead to the presence, within the plasma volume, of particles or plasma structure defined by different fractalization and scale resolutions. The simulations also showcase a transition from quasi mono-energetic particles ejected through a single ablation mechanism to an amalgam of particles with different energetic distributions specific to each ejection mechanism involved.
In Figure 1b we represented the particle distribution with the scale resolution for various moments in time. As the plasma evolves we see the presence of multiple distribution sites centered across low (s < 0.3) and high (0.4 < s < 1) resolution scales. This suggests that the first ejected particles (∆T = 40) are defined only by one scale resolution, characteristic of the Coulomb explosion mechanism. For longer period of time, when the thermal mechanism is dominant, we see multiple distributions centered on higher values of the scale resolutions. Therefore, our simulations capture the ejection of a series of particles with different kinetic energies and fractalizations. As such, we can now define a scale resolution for each of the ejection mechanisms (~0.2 for Coulomb Explosion,~0.6 for Explosive Boiling and~1 for the removal of complex structures) ensuring a wide covering of our theoretical model, transcending the Coulomb temporal scale, up to the explosive boiling ejection scenario and the ejection of clusters and more complex structures. image of the laser produced plasma showcases the appearance of a one particle distribution centered on a relatively high velocity, at low resolutions scale and a consistent increase number of distributions centered around lower velocities. The presence of multiple structures in laser produced plasma has been previously attested and reported as a direct result of the multiple ejection mechanism and plasma interaction with the background gas. We note that the presence of multiple ejection mechanisms [17][18][19][20] (Coulomb explosion, explosive boiling, phase explosion, etc.) will lead to the presence, within the plasma volume, of particles or plasma structure defined by different fractalization and scale resolutions. The simulations also showcase a transition from quasi monoenergetic particles ejected through a single ablation mechanism to an amalgam of particles with different energetic distributions specific to each ejection mechanism involved. In Figure 1b we represented the particle distribution with the scale resolution for various moments in time. As the plasma evolves we see the presence of multiple distribution sites centered across low (s < 0.3) and high (0.4 < s < 1) resolution scales. This suggests that the first ejected particles (ΔT = 40) are defined only by one scale resolution, characteristic of the Coulomb explosion mechanism. For longer period of time, when the thermal mechanism is dominant, we see multiple distributions centered on higher values of the scale resolutions. Therefore, our simulations capture the ejection of a series of particles with different kinetic energies and fractalizations. As such, we can now define a scale resolution for each of the ejection mechanisms (~ 0.2 for Coulomb Explosion, ~0.6 for Explosive Boiling and ~1 for the removal of complex structures) ensuring a wide covering of our theoretical model, transcending the Coulomb temporal scale, up to the explosive boiling ejection scenario and the ejection of clusters and more complex structures. In Figure 2a we have represented the charge density fluctuations induced by some unbalances with respect to the plasma volume ( ). We can see that for short moments of time we observed a separation of the charges. The obtained distribution in the fractal space resembles the double layer distribution in the case of classical plasma physics [18]. As the time evolves, the double-layer-like distribution moves towards higher distances, and for longer moments of time a secondary one appears. This behavior is in line with the transient double layer scenario published by Bulgakova [20]. The charge separation that occurs during the Coulomb explosion will act as a driving force behind the expansion of the particles. The double layer formed during the particle expansion will lead to the acceleration of ions and the deceleration of the electron. This exchange will induce some current oscillations [21] with the frequency depending of the nature of the material, and thus on its fractality. The formation of the second plasma structure, expanding with a lower velocity and being In Figure 2a we have represented the charge density fluctuations induced by some unbalances with respect to the plasma volume (ξ). We can see that for short moments of time we observed a separation of the charges. The obtained distribution in the fractal space resembles the double layer distribution in the case of classical plasma physics [18]. As the time evolves, the double-layer-like distribution moves towards higher distances, and for longer moments of time a secondary one appears. This behavior is in line with the transient double layer scenario published by Bulgakova [20]. The charge separation that occurs during the Coulomb explosion will act as a driving force behind the expansion of the particles. The double layer formed during the particle expansion will lead to the acceleration of ions and the deceleration of the electron. This exchange will induce some current oscillations [21] with the frequency depending of the nature of the material, and thus on its fractality. The formation of the second plasma structure, expanding with a lower velocity and being described as having a different temperature [22], will also lead to the formation of a double layer between the fast structure and slow structure [20]. This phenomenon was showcased by us in Figure 2a, where we can see the appearance of a second signature of the double layer in the charged particle distribution at a later evolution time. The intensity of the function represented in Figure 2a increases with the increase of time and we can see that at longer expansion time, where the formation of a second plasma structure becomes visible, another more intense double layer forms. This result was expected, as the multiple structures formed within the plasma volume suffered a spatial and temporal expansion. In order to maintain the composed global shape and to compensate the loses, particle density, with an increase in the separation between the two structures, of the electrical field defining the double layers increase their amplitudes. In Figure 2b, we have represented the particle current evolution in time for various scale resolutions. We noticed that the particles present an oscillatory dynamic, as was previously reported by [23] and experimentally proven in [24] and [25]. The frequency changes with the scale resolution, thus we deduce that each component of the plasma might present a different oscillatory behavior dictated by the characteristics of the double layer formed in the area separating the structures. with the increase of time and we can see that at longer expansion time, where the formation of a second plasma structure becomes visible, another more intense double layer forms. This result was expected, as the multiple structures formed within the plasma volume suffered a spatial and temporal expansion. In order to maintain the composed global shape and to compensate the loses, particle density, with an increase in the separation between the two structures, of the electrical field defining the double layers increase their amplitudes. In Figure 2b, we have represented the particle current evolution in time for various scale resolutions. We noticed that the particles present an oscillatory dynamic, as was previously reported by [23] and experimentally proven in [24] and [25]. The frequency changes with the scale resolution, thus we deduce that each component of the plasma might present a different oscillatory behavior dictated by the characteristics of the double layer formed in the area separating the structures.

Experimental Confirmation
In order to verify if our theoretical assumptions can have experimental correspondents we investigated by means of ICCD fast camera imaging and space-and-time-resolved optical emission spectroscopy of a plasma generated by an ns laser beam on a Chalcopyrite mineral target. The experiments were performed in fixed external conditions (laser fluence 5 J/cm 2 and background pressure of 10 -2 Torr). The choice of the relatively simple mineral comes from its composition, having elements with different physical properties (S, Cu and Fe) which will allow a better showcase of phenomena like: particle separation, ionic oscillations and plume splitting. Further details on the experimental set-up can be found in [6,12].
In Figure 3a we have represented the spatial distribution of spectral region 420-430 nm after a time delay of 650 ns. This region is significative for our laser produced plasmas as it contains all the elements composing the target (Cu I-II, Fe I-II and S II). We noticed that the atomic emission lines can be found only for short distance while the ionic lines start and end emission at considerably longer distances. This result confirms other findings from literature and underlines the separation of the plasma components into fast and slow [26], respecting the ejection mechanism behind each species. In Figure 3b we have plotted the spatial distribution of representative emission lines for all the atomic and ionic species. Due to the difference in expansion velocities between ionic (S II-15.4 km/s, Fe II-14.3 km/s and Cu II-11.4 km/s) and atomic species (Fe I-6.2 km/s and Cu I-4.4 km/s), the ions can be seen expanding at longer distances with respect to the target surface. The slightly increased velocity of S ions compared with Fe or Cu ones is a direct consequence of the acceleration in the double layer generated through Coulomb expansion [27]. We also notice periodic fluctuations (oscillations) on the spatial distribution of all the ionic species. The period of these oscillation are of approximately 900

Experimental Confirmation
In order to verify if our theoretical assumptions can have experimental correspondents we investigated by means of ICCD fast camera imaging and space-and-time-resolved optical emission spectroscopy of a plasma generated by an ns laser beam on a Chalcopyrite mineral target. The experiments were performed in fixed external conditions (laser fluence 5 J/cm 2 and background pressure of 10 −2 Torr). The choice of the relatively simple mineral comes from its composition, having elements with different physical properties (S, Cu and Fe) which will allow a better showcase of phenomena like: particle separation, ionic oscillations and plume splitting. Further details on the experimental set-up can be found in [6,12].
In Figure 3a we have represented the spatial distribution of spectral region 420-430 nm after a time delay of 650 ns. This region is significative for our laser produced plasmas as it contains all the elements composing the target (Cu I-II, Fe I-II and S II). We noticed that the atomic emission lines can be found only for short distance while the ionic lines start and end emission at considerably longer distances. This result confirms other findings from literature and underlines the separation of the plasma components into fast and slow [26], respecting the ejection mechanism behind each species. In Figure 3b we have plotted the spatial distribution of representative emission lines for all the atomic and ionic species. Due to the difference in expansion velocities between ionic (S II-15.4 km/s, Fe II-14.3 km/s and Cu II-11.4 km/s) and atomic species (Fe I-6.2 km/s and Cu I-4.4 km/s), the ions can be seen expanding at longer distances with respect to the target surface. The slightly increased velocity of S ions compared with Fe or Cu ones is a direct consequence of the acceleration in the double layer generated through Coulomb expansion [27]. We also notice periodic fluctuations (oscillations) on the spatial distribution of all the ionic species. The period of these oscillation are of approximately 900 kHz in good agreement with other reports of ionic oscillations in laser produced plasmas [21,24,28]. The first attempts for the comprehension of this "peculiar" behavior was based on the formation of single or multiple double-layers in the very vicinity of the target. This picture was the main focus to a long series of papers reporting on charge separation in laser-produced plasma, mainly from the 1980s [29,30]. Eliezer and Hora [23] gathered, in a very comprehensive manner, the state of the art regarding the double and multiple layers in laser-produced plasmas. One of the remarkable results reported are experimental proofs with double-layer electric fields of 105-106 V/cm and widths of 10-100 Debye lengths. In the past few years, three other theoretical approaches where proposed. One based on the fractal model developed as the interaction between two fractal structures [5,31], and their corresponding interface (generally, this interface delineates the double layer), with the second ones being on differential physics [28]: a collisional model based on the plasma ion frequency and electron-ion collision rate in the context of the Lieberman's model for plasma immersion ion implantation, and finally, one based on the AC Josephson effect. So, at this point there is no real consensus for the real mechanism behind the oscillatory behavior but intense theoretical and experimental work is undergoing to shed some light on it. regarding the double and multiple layers in laser-produced plasmas. One of the remarkable results reported are experimental proofs with double-layer electric fields of 105-106 V/cm and widths of 10-100 Debye lengths. In the past few years, three other theoretical approaches where proposed. One based on the fractal model developed as the interaction between two fractal structures [5,31], and their corresponding interface (generally, this interface delineates the double layer), with the second ones being on differential physics [28]: a collisional model based on the plasma ion frequency and electron-ion collision rate in the context of the Lieberman's model for plasma immersion ion implantation, and finally, one based on the AC Josephson effect. So, at this point there is no real consensus for the real mechanism behind the oscillatory behavior but intense theoretical and experimental work is undergoing to shed some light on it. To verify how the distribution of the compositional species within the plasma volume would affect the overall behavior of the laser produced plasma, we recorded the overall emission of the plasma at various moments in time. In Figure 4a we show a representative image of the plasma at a time-delay of 650 ns. If we perform cross section across the main expansion axis (centered around the orthogonal direction on the target in the impact point), three areas are more visible as the time delay increases. The first area seen at a larger distance corresponds to an ion rich area (confirmed by the data presented in Figure 3b) expanding with 16 km/s; the second one corresponds to an atomic rich area (~6 km/s); and the emission maxima corresponds to the maxima of the atomic emission. The last area has a small intensity and it is in the proximity of the target. It corresponds to the emission from atomic species colliding with clusters and microdroplets ejected from the target. This emission was not seen in the spectra's resolved measurements but remains a trademark of a strong thermal mechanism and the ejection of clusters and nanoparticles [32]. The difference in expansion velocity is a signature of the fundamental ejection mechanism behind each structure. The first structure contains mainly ions and it is induced during Coulomb explosion. The second structure contains mainly atomic species ejected by the thermal mechanism which require a longer incubation time (up to a few ns), thus displaying a lower velocity. The third structure mainly contains nanostructures, or clusters ejected directly from the target and expanding with a significantly lower velocity (generally found in the order of hundreds of m/s). To verify how the distribution of the compositional species within the plasma volume would affect the overall behavior of the laser produced plasma, we recorded the overall emission of the plasma at various moments in time. In Figure 4a we show a representative image of the plasma at a time-delay of 650 ns. If we perform cross section across the main expansion axis (centered around the orthogonal direction on the target in the impact point), three areas are more visible as the time delay increases. The first area seen at a larger distance corresponds to an ion rich area (confirmed by the data presented in Figure 3b) expanding with 16 km/s; the second one corresponds to an atomic rich area (~6 km/s); and the emission maxima corresponds to the maxima of the atomic emission. The last area has a small intensity and it is in the proximity of the target. It corresponds to the emission from atomic species colliding with clusters and microdroplets ejected from the target. This emission was not seen in the spectra's resolved measurements but remains a trademark of a strong thermal mechanism and the ejection of clusters and nanoparticles [32]. The difference in expansion velocity is a signature of the fundamental ejection mechanism behind each structure. The first structure contains mainly ions and it is induced during Coulomb explosion. The second structure contains mainly atomic species ejected by the thermal mechanism which require a longer incubation time (up to a few ns), thus displaying a lower velocity. The third structure mainly contains nanostructures, or clusters ejected directly from the target and expanding with a significantly lower velocity (generally found in the order of hundreds of m/s). The experimental data extracted from a plasma generated on chalcopyrite mineral by irradiation with an ns laser beam confirms the theoretical projections presented in Section 3. The temporal separation the plasma components based on their inner properties accurately reflects the separation based on the resolution scale of each type of ablation mechanism and the presence of multiple distribution on the scale resolution representation. One of the most important results is the prediction and confirmation of ionic oscillations by using invasive techniques. The oscillatory behavior coupled with heterogenic dynamics of different ions is well in line with the image depicted by other groups [20,33] and by our group in past papers [5,13,31]. In this paper the oscillatory behavior and the double layer-like distribution appear as natural solutions to the initial paradigm which translated the initial complex Lorenz system from laser to target to plasma. As such, we were able to build a robust theoretical model that can contains the set of parameters for the laser, target and the projection of the particle dynamic after ejection as solutions to the initial system.

Conclusions
A non-differential Lorenz system was built by projecting a differential Lorenz system on a fractal space. Simulations are performed for a wide range of scale resolutions showcasing the appearance of multiple distribution centered on different velocities attributed to the various plasma formed through different removal mechanisms. Current oscillations were also predicted as a result of the appearance of multiple double layers during expansion.
The theoretical simulations were confronted with experimental data extracted by means of ICCD fast camera imaging and space and time resolved optical emission spectroscopy of a complex plasma generated by ns-laser ablation on a chalcopyrite sample. Space and time resolved measurement revealed an oscillating behavior seen in the emission of the Cu, Fe and S ions. The ions were found to expand with various velocities specific to each species present in the plasma. The ICCD imaging revealed the split into two structures (fast and slow) expanding with different velocities. The values were found consistent with the ones of the individual species seen through spectrally resolved measurements. The experimental data is in good agreement with the major predictions made by the theoretical model based on the Lorenz system.  The experimental data extracted from a plasma generated on chalcopyrite mineral by irradiation with an ns laser beam confirms the theoretical projections presented in Section 3. The temporal separation the plasma components based on their inner properties accurately reflects the separation based on the resolution scale of each type of ablation mechanism and the presence of multiple distribution on the scale resolution representation. One of the most important results is the prediction and confirmation of ionic oscillations by using invasive techniques. The oscillatory behavior coupled with heterogenic dynamics of different ions is well in line with the image depicted by other groups [20,33] and by our group in past papers [5,13,31]. In this paper the oscillatory behavior and the double layer-like distribution appear as natural solutions to the initial paradigm which translated the initial complex Lorenz system from laser to target to plasma. As such, we were able to build a robust theoretical model that can contains the set of parameters for the laser, target and the projection of the particle dynamic after ejection as solutions to the initial system.

Conclusions
A non-differential Lorenz system was built by projecting a differential Lorenz system on a fractal space. Simulations are performed for a wide range of scale resolutions showcasing the appearance of multiple distribution centered on different velocities attributed to the various plasma formed through different removal mechanisms. Current oscillations were also predicted as a result of the appearance of multiple double layers during expansion.
The theoretical simulations were confronted with experimental data extracted by means of ICCD fast camera imaging and space and time resolved optical emission spectroscopy of a complex plasma generated by ns-laser ablation on a chalcopyrite sample. Space and time resolved measurement revealed an oscillating behavior seen in the emission of the Cu, Fe and S ions. The ions were found to expand with various velocities specific to each species present in the plasma. The ICCD imaging revealed the split into two structures (fast and slow) expanding with different velocities. The values were found consistent with the ones of the individual species seen through spectrally resolved measurements. The experimental data is in good agreement with the major predictions made by the theoretical model based on the Lorenz system.