@@ -11,3 +11,7 @@ Parameters leading to the plots displayed in the report are listed in captions.

## Data analysis

Integration data are saved in .npy files thanks to the save function of Numpy. Then they are loaded into a analysis notebooks in which plots are implemented. **NB:** for the heavy integrations required by full filamentation data is not saved, because of the size of field and charge matrices. Analysis is done just after integration in the notebook and plots are saved.

## Filamentation Loop

In order to integrate the filamentation equation for various input powers, the script /codes/production/CN_classes/filamentation_loop.py can be run. A list of entire multiples of the critical power can be chosen and filamentation is calculated with those values. Resulting matrices are not stored but plots of maximum field intensity and charge density are saved.

where $K$ is the number of photons absorbed at a time, $K =7$ for air, and $\beta^{(K)}$ accounts for the efficiency of this absorption.

where $K$ is the number of photons absorbed at a time, $K =7$ for air and the wavelength considered hereafter, and $\beta^{(K)}$ accounts for the efficiency of this absorption.

The second defocusing effects is due to the \textbf{reaction of charges} in the plasma to the excitation of the electric field. Considering a Drude model for electrons, their answer to a periodic excitation creates a delay in their \textbf{polarisation}. This answer reads

\begin{equation*}

...

...

@@ -346,7 +346,9 @@ The different Crank-Nicholson solvers, according to the effects implemented, are

Each of those classes are accompanied by a Jupyter Notebook named \texttt{*\_script.ipynb}, which instantiates the solver and runs the integration. Then, for linear integrations, data are saved in a \texttt{.npy} file and analysed in another Jupyter Notebook at \texttt{/codes/analyse}. For non linear integrations, results are analysed directly in the integration notebook.

An independant script, \texttt{filamentation\_loop.py}, stored with solving classes, can be run to calculate filamentation for a list of input powers. In this case, resulting data is not returned and cannot be accessed, but plots of maximum field intensity, along propagation axis, and of charge density created by the pulse, along propagation axis, are saved.

The parameters used to create the results displayed in this report are listed in the captions of figures, and some options are proposed in the notebooks, to test the influence of the various effects involved.