Publications

2013

• Absolute atomic oxygen and nitrogen densities in radio-frequency driven atmospheric pressure cold plasmas: Synchrotron vacuum ultra-violet high-resolution Fourier-transform absorption measurements
  
  • Niemi K.
  • O'Connell D.
  • de Oliveira N.
  • Joyeux D.
  • Nahon L.
  • Booth Jean-Paul
  • Gans T.
  Applied Physics Letters, American Institute of Physics, 2013, 103 (3), pp.034102. Reactive atomic species play a key role in emerging cold atmospheric pressure plasma applications, in particular, in plasma medicine. Absolute densities of atomic oxygen and atomic nitrogen were measured in a radio-frequency driven non-equilibrium plasma operated at atmospheric pressure using vacuum ultra-violet (VUV) absorption spectroscopy. The experiment was conducted on the DESIRS synchrotron beamline using a unique VUV Fourier-transform spectrometer. Measurements were carried out in plasmas operated in helium with air-like N2/O2 (4:1) admixtures. A maximum in the O-atom concentration of (9.1 ± 0.7)×1020 m−3 was found at admixtures of 0.35 vol. %, while the N-atom concentration exhibits a maximum of (5.7 ± 0.4)×1019 m−3 at 0.1 vol. %. (10.1063/1.4813817)
  DOI : 10.1063/1.4813817
• Design and Testing of a Surface Switch for the Dynamic Load Current Multiplier on the SPHINX Microsecond LTD
  
  • Maysonnave Thomas
  • Bayol Frédéric
  • Demol Gauthier
  • d'Almeida Thierry
  • Morell Alain
  • Lassalle Francis
  • Grunenwald Julien
  • Chuvatin Alexandre S.
  • Pecastaing Laurent
  • de Ferron Antoine Silvestre
  IEEE Transactions on Plasma Science, Institute of Electrical and Electronics Engineers, 2013, 41 (10), pp.2593 - 2599. SPHINX is a microsecond linear transformer driver located at Atomic Energy Comission (CEA) Gramat (France), which can deliver a current pulse of 6 MA within 800 ns in a Z -pinch load. Using the concept of the dynamic load current multiplier (DLCM), which was proposed by Chuvatin, we expect to increase the load current above 6 MA, while decreasing its rise time to ~ 300 ns. The DLCM developed by the CEA Gramat and International Technologies for High Pulsed Power (ITHPP) is a compact system made up of concentric electrodes (autotransformer), a dynamic flux extruder (cylindrical wire array), a vacuum convolute (eight post-hole rods), and a closing switch (compact vacuum surface switch). The latter is a key component of the system, which is used to prevent the current from flowing into the load until the inductance builds up due to the implosion of the wire array. This paper presents the design and testing of the DLCM surface switch, resulting from both electrostatic simulations and experiments on the SPHINX generator. These studies, carried out either with or without load (open circuit), were valuable for a first experimental evaluation of the DLCM scheme in a microsecond regime and provided detailed information on the surface switch behavior. (10.1109/TPS.2013.2245148)
  DOI : 10.1109/TPS.2013.2245148
• Redistribution of high energy alpha particles due to sawteeth with partial reconnection
  
  • Farengo R.
  • Ferrari H.E.
  • García-Martínez P.L.
  • Firpo Marie-Christine
  • Ettoumi Wahb
  • Lifschitz A.F.
  Nuclear Fusion, IOP Publishing, 2013, 53, pp.043012. The redistribution of high energy alpha particles due to internal kink modes is studied in plasmas with ITER-like parameters. The exact particle trajectories in the total fields, equilibrium plus perturbation, are calculated. The equilibrium magnetic field is obtained by analytically solving the GradShafranov equation and the perturbed electric and magnetic fields are reconstructed using ideal MHD and the experimental information about the displacement eigenfunction. The (1, 1), (2, 2) and (2, 1) modes are included and the effect of changing their amplitude and frequency is determined. The results show that if the conditions are similar to those reported in Igochine et al (2007 Nucl. Fusion 47 23), the peak density of counter-passing particles decreases between 25% and 40% (depending on the energy); the peak of the trapped particles density shifts outwards by approximately 10% of the minor radius and the total on axis density decreases by more than 25%. This redistribution occurs inside the q = 1 surface. The addition of a (2, 1) mode, which can produce the stochastization of the magnetic field, significantly increases particle redistribution and allows particles to spread beyond the q = 1 surface. Different groups of particles (co-passing, counter-passing, trapped) respond differently to the perturbations. (10.1088/0029-5515/53/4/043012)
  DOI : 10.1088/0029-5515/53/4/043012
• Structure of nonlocality of plasma turbulence
  
  • Gürcan Özgür D.
  • Vermare Laure
  • Hennequin Pascale
  • Berionni Vincent
  • Diamond P.H.
  • Dif-Pradalier Guilhem
  • Garbet X.
  • Ghendrih Philippe
  • Grandgirard Virginie
  • Mcdevitt C.J.
  • Morel Pierre
  • Sarazin Y.
  • Storelli A.
  • Bourdelle C.
  • Tore Supra Team
  Nuclear Fusion, IOP Publishing, 2013, 53, pp.073029. Various indications on the weakly nonlocal character of turbulent plasma transport both from experimental fluctuation measurements from Tore Supra and observations from the full-f, flux-driven gyrokinetic code GYSELA are reported. A simple Fisher equation model of this weakly nonlocal dynamics can be formulated in terms of an evolution equation for the turbulent entropy density, which contains the basic phenomenon of radial turbulence spreading in addition to avalanche-like dynamics via coupling to profile modulations. A derivation of this model, which contains the so-called beach effect, a diffusive and convective flux components for the flux of turbulence intensity, in addition to linear group propagation is given, starting from the drift-kinetic equation. The proposed model has the form of a transport equation for turbulence intensity, and may be considered as an addition to transport modelling. The kinetic fluxes given, can be computed using model closures, or local gyrokinetics. The model is also used in a particular setup that represents the near edge region as a relatively stable zone between the core and edge region where the energy injection is locally more substantial. It is observed that with constant, physical coefficients, the model gives a convincing qualitative profile of fluctuation intensity when the turbulence is coming from the core region with either a group velocity or a convective flux. (10.1088/0029-5515/53/7/073029)
  DOI : 10.1088/0029-5515/53/7/073029
• Localizing Transport Barriers in Degenerate 3/2 D.O.F. Hamiltonian Systems with Applications to Magnetic Confinement Fusion
  
  • Constantinescu D.
  • Firpo Marie-Christine
  International journal of bifurcation and chaos in applied sciences and engineering, World Scientific Publishing, 2013, 23, pp.1350034. The existence of transport barriers in 3/2 degrees of freedom degenerate Hamiltonian systems is studied using the associated stroboscopic maps. For small enough amplitudes of the perturbations, a transport barrier (formed by infinitely many invariant rotational circles) is proven to exist and to form in the degenerate annulus. These results are applied to Hamiltonian models which describe some magnetic configurations in tokamak plasmas. Read More: http://www.worldscientific.com/doi/abs/10.1142/S021812741350034X (10.1142/S021812741350034X)
  DOI : 10.1142/S021812741350034X
• Particle-in-cell simulation of an electronegative plasma under direct current bias studied in a large range of electronegativity
  
  • Oudini N.
  • Raimbault Jean-Luc
  • Chabert Pascal
  • Meige A.
  • Aanesland Ane
  Physics of Plasmas, American Institute of Physics, 2013, 20, pp.043501. A one-dimensional electronegative plasma situated between two symmetrical parallel electrodes under DC bias is studied by Particle-In-Cell simulation with Monte Carlo Collisions. By varying the electronegativity &#945; &#8801; n&#8722;/ne from the limit of electron-ion plasmas (negative ion free) to ion-ion plasmas (electron free), the sheaths formation, the negative ion flux flowing towards the electrodes, and the particle velocities at the sheath edges are investigated. Depending on &#945;, it is shown that the electronegative plasma behavior can be described by four regimes. In the lowest regime of &#945;, i.e., &#945; < 50, negative ions are confined by two positive sheaths within the plasma, while in the higher regimes of &#945;, a negative sheath is formed and the negative ion flux can be extracted from the bulk plasma. In the two intermediate regimes of &#945;, i.e., 50 < &#945; < 10⁵, both the electron and the negative ion fluxes are involved in the neutralization of the positive ions flux that leaves the plasma. In particular, we show that the velocity of the negative ions entering the negative sheath is affected by the presence of the electrons, and is not given by the modified Bohm velocity generally accepted for electronegative plasmas. For extremely high electronegativity, i.e., &#945; > 10⁵, the presence of electrons in the plasma is marginal and the electronegative plasma can be considered as an ion-ion plasma (electron free). (10.1063/1.4798501)
  DOI : 10.1063/1.4798501
• Normal regime of the weak-current mode of an rf capacitive discharge
  
  • Lisovskiy V. A.
  • Yegorenkov V. D.
  • Artushenko E.
  • Booth Jean-Paul
  • Martins S.
  • Landry K.
  • Douai D.
  • Cassagne V.
  Plasma Sources Science and Technology, IOP Publishing, 2013, 22 (1), pp.015018. This paper studies the normal and abnormal regimes of a weak-current rf discharge in ammonia, nitrogen, hydrogen and N 2 O for the rf electric field frequencies of 13.56 and 27.12 MHz. We reveal that only the abnormal regime of burning is observed at low pressures when the current growth is accompanied by an rf voltage increase while the surface of the electrodes is completely covered with the discharge. The normal regime occurs at higher gas pressures when the current growth is due to the increase in the surface area occupied by the discharge on the electrodes. The discharge burns in the abnormal mode after the surface of the electrodes is completely covered with the discharge. We demonstrate that the normal current density is directly proportional to the gas pressure and it depends approximately on the square of the rf electric field frequency. We present an analytical model for two limiting cases: constant free path length and constant mobility of positive ions furnishing a satisfactory description of the experimental data. (10.1088/0963-0252/22/1/015018)
  DOI : 10.1088/0963-0252/22/1/015018
• Cluster Observations of Ion-Scale Magnetic Structures and their Coupling with Whistler Waves During the August 17th 2003 Substorm Event
  
  • Tenerani Anna
  • Le Contel Olivier
  • Califano F.
  • Robert Patrick
  • Fontaine Dominique
  • Cornilleau-Wehrlin Nicole
  • Sauvaud J.-A.
  Journal of Geophysical Research Space Physics, American Geophysical Union/Wiley, 2013, 118, pp.6072-6089. We provide evidence of the simultaneous occurrence of large-amplitude, quasi-parallel whistler mode waves and ion-scale magnetic structures, which have been observed by the Cluster spacecraft in the plasma sheet at 17 Earth radii, during a substorm event. It is shown that the magnetic structures are characterized by both a magnetic field strength minimum and a density hump and that they propagate in a direction quasi-perpendicular to the average magnetic field. The observed whistler mode waves are efficiently ducted by the inhomogeneity associated with such ion-scale magnetic structures. The large amplitude of the confined whistler waves suggests that electron precipitations could be enhanced locally via strong pitch angle scattering. Furthermore, electron distribution functions indicate that a strong parallel heating of electrons occurs within these ion-scale structures. This study provides new insights on the possible multiscale coupling of plasma dynamics during the substorm expansion, on the basis of the whistler mode wave trapping by coherent ion-scale structures. (10.1002/jgra.50562)
  DOI : 10.1002/jgra.50562
• Non thermal Pasma regeneration of Acetone adsorbed TiO<SUB>2</SUB> Surface
  
  • Sivachandiran Loganathan
  • Thévenet Frédéric
  • Rousseau Antoine
  Plasma Chemistry and Plasma Processing, Springer Verlag, 2013. Improvement of indoor air quality regarding volatile organic compounds (VOCs) requires the development of innovative oxidation processes. This paper investigates the coupling of a metal oxide sorbent with non-thermal plasma (NTP) in an especially designed reactor. TiO2 was selected as model sorbent and acetone was used as model VOC. The analyses of gas phase species at the reactor downstream have been performed using FTIR spectroscopy. In a first step, acetone adsorption on TiO2 surface under dry air was characterized in terms of total amount adsorbed, as well as reversibly and irreversibly adsorbed fractions. Obtained results were compared and discussed with literature in terms of acetone reactive adsorption on TiO2 surface. Mesityloxide was proposed as the major compound in the irreversibly adsorbed fraction. In a second time, acetone saturated TiO2 surface was exposed to NTP surface discharge. Irrespectively of the injected power, <30 % of the initially adsorbed acetone has been recovered as CO, CO2 and desorbed acetone. Finally, thermal desorptions have been performed. They evidenced that (1) NTP treatment modifies the nature of the adsorbed organic species, (2) mineralization rate is considerably improved. Based on desorbed species temporal profile analysis, carboxylates and more especially formates are suggested as major adsorbed species after NTP treatment (Pinj > 0.2 W). This hypothesis has been evaluated and confirmed. This paper finally evidenced that NTP can be used as an efficient pretreatment technique to promote the mineralization of adsorbed acetone for further thermal treatment (10.1007/s11090-013-9463-7)
  DOI : 10.1007/s11090-013-9463-7
• Investigation of NO and NO<SUB>2</SUB> adsorption mechanisms on TiO<SUB>2</SUB> at room temperature
  
  • Sivachandiran Loganathan
  • Thévenet Frédéric
  • Gravelat Paul
  • Rousseau Antoine
  Applied Catalysis B: Environmental, Elsevier, 2013, 142-143, pp.196-204. In this study, NO and NO2 adsorption mechanisms on TiO2 at room temperature have been investigated separately. Atmospheric pressure gas phase Fourier Transform Infrared Spectroscopy (FTIR) coupled to adsorption and temperature programmed desorption (TPD) experiments have been used to characterize gas phase and adsorbed phase species. TiO2 coated glass beads packed-bed reactor has been designed and used for NO and NO2 adsorption and storage capacity under dark and gas-flowing conditions at room temperature. For NO adsorption, experimental data suggests that, at room temperature, NO shows no significant adsorption on TiO2. On the other hand NO2 adsorbs in a reactive way by evolving NO in the gas phase. We quantitatively evidenced that, the ratio between consumed NO2, desorbed NO2 by TPD after adsorption and produced NO during NO2 adsorption is 3:2:1, by a qualitative and quantitative analysis performed downstream the reactor. Based on the quantitative analysis, a new NO2 adsorption mechanism on TiO2 at room temperature has been proposed, and validated for various NO2 inlet concentrations. The proposed mechanism is valid on the investigated NO2 inlet concentration range. In addition, it was found that, NO formation time, during NO2 adsorption, is rather controlled by TiO2 surface coverage than NO2 inlet concentration; similarly, adsorption time i.e. surface coverage, significantly modifies the nature of adsorbed species on TiO2 surface. (10.1016/j.apcatb.2013.04.073)
  DOI : 10.1016/j.apcatb.2013.04.073
• Modes of underwater discharge propagation in a series of nanosecond successive pulses
  
  • Marinov Ilya
  • Guaitella Olivier
  • Rousseau Antoine
  • Starikovskaia Svetlana
  Journal of Physics D: Applied Physics, IOP Publishing, 2013, 46, pp.464013. Initiation and development of nanosecond electric discharges in deionized water was studied by time-resolved shadowgraphy and emission imaging. Pulses of 410 kV in amplitude, 5 ns rise time and 30 ns duration were applied to a high-voltage (HV) pin-like electrode with the diameter of the pin tip about 2&#956;m. Two distinct bush-like and tree-like modes of the discharge propagation were observed in the same conditions already at the rising slope of HV pulse. Both the modes are supersonic. The probability of finding one of them depends on applied voltage. Bush-mode propagates at 4 kms&#8722;1 at the initial stage; the deposited energy is less than 0.05 mJ. For the tree-like mode, propagation velocity, number of branches and deposited energy show well-prononced dependence on applied voltage. It is suggested that the bush-like discharge ignites in the gaseous cavity, whereas the tree-like discharge propagates in bulk liquid. The role of electrostriciton in the initiation and propagation of two modes is discussed, and comparison with microsecond discharges is given. Transition from bush-like discharge to tree-like mode may occur in successive positive pulse coming 500 ns after initial pulse. (10.1088/0022-3727/46/46/464013)
  DOI : 10.1088/0022-3727/46/46/464013
• Frequency dependence of the electrical asymmetry effect in dual-frequency capacitively coupled discharges
  
  • Booth Jean-Paul
  • Lafleur Trevor
  Applied Physics Letters, American Institute of Physics, 2013, 102, pp.154104. Using experimental measurements complemented with particle-in-cell (PIC) simulations, we demonstrate a reduction in the electrical asymmetry effect in capacitively coupled discharges for low excitation frequencies (0.8&#8201;MHz&#8201;<&#8201;f<10&#8201;MHz). These results confirm PIC simulation predictions made by [I. Korolov, Z. Donkó, U. Czarnetzki, and J. Schulze, J. Phys. D: Appl. Phys. 45, 465205 (2012)], where this effect was investigated for dual-frequency discharges in which the driving frequencies are phase-shifted harmonics. Because the reduction in electrical asymmetry is a sensitive function of the secondary electron emission coefficient, we are able to non-invasively estimate this coefficient as 0.035 for argon ions incident on aluminium electrodes. (10.1063/1.4802241)
  DOI : 10.1063/1.4802241
• Sahraoui et al. Reply:
  
  • Sahraoui Fouad
  • Robert Patrick
  • Goldstein M. L.
  • Khotyaintsev Y. V.
  Physical Review Letters, American Physical Society, 2013, 111, pp.149002. A Reply to the Comment by O. Alexandrova, S. D. Bale, and C. Lacombe.
• First demonstration of an asymmetric kinetic equilibrium for a thin current sheet
  
  • Aunai Nicolas
  • Belmont Gérard
  • Smets Roch
  Physics of Plasmas, American Institute of Physics, 2013, 20 (11), pp.110702. The modeling of steady state collisionless asymmetric tangential current layers is a challenging and poorly understood problem. For decades now, this difficulty has been limiting numerical models to approximate equilibria built with locally Maxwellian current layers and theoretical analyses to the very restricted Harris equilibrium. We show how the use of any distribution functions depending only on local macroscopic quantities results in a strong alteration of the current layer internal structure, which converges toward an unpredictable quasi-steady state with emission of ion scale perturbations. This transient can be explained in terms of ion kinetic and electron fluid physics. We demonstrate, for the first time, the validity of an asymmetric kinetic equilibrium model as well as its usability as an initial condition of hybrid kinetic simulations. This offers broad perspectives for the current sheet modeling, for which the early phase of instabilities can be studied within the kinetic formalism. (10.1063/1.4833679)
  DOI : 10.1063/1.4833679
• Magnetic field and dynamic pressure ULF fluctuations in coronal-mass-ejection-driven sheath regions
  
  • Kilpua E. K. J.
  • Hietala H.
  • Koskinen H. E. J.
  • Fontaine Dominique
  • Turc Lucile
  Annales Geophysicae, European Geosciences Union, 2013, 31, pp.1559-1567. Compressed sheath regions form ahead of interplanetary coronal mass ejections (ICMEs) that are sufficiently faster than the preceding solar wind. The turbulent sheath regions are important drivers of magnetospheric activity, but due to their complex internal structure, relatively little is known on the distribution of the magnetic field and plasma variations in them. In this paper we investigate ultra low frequency (ULF) fluctuations in the interplanetary magnetic field (IMF) and in dynamic pressure (P<SUB>dyn</SUB>) using a superposed epoch analysis of 41 sheath regions observed during solar cycle 23. We find strongest fluctuation power near the shock and in the vicinity of the ICME leading edge. The IMF and P<SUB>dyn</SUB> ULF power have different profiles within the sheath; the former is enhanced in the leading part of the sheath, while the latter is increased in the trailing part of the sheath. We also find that the ICME properties affect the level and distribution of the ULF power in sheath regions. For example, sheath regions associated with strong or fast ICMEs, or those that are crossed at intermediate distances from the center, have strongest ULF power and large variation in the power throughout the sheath region. The weaker or slower ICMEs, or those that are crossed centrally, have in general considerably weaker ULF power with relatively smooth profiles. The strong and abrupt decrease of the IMF ULF power at the ICME leading edge could be used to distinguish the ICME from the preceding sheath plasma. (10.5194/angeo-31-1559-2013)
  DOI : 10.5194/angeo-31-1559-2013
• Enhanced Magnetic Compressibility and Isotropic Scale Invariance at Sub-ion Larmor Scales in Solar Wind Turbulence
  
  • Kiyani K. H.
  • Chapman S. C.
  • Sahraoui Fouad
  • Hnat B.
  • Fauvarque O.
  • Khotyaintsev Y. V.
  The Astrophysical Journal, American Astronomical Society, 2013, 763 (1), pp.10. The anisotropic nature of solar wind magnetic turbulence fluctuations is investigated scale by scale using high cadence in situ magnetic field measurements from the Cluster and ACE spacecraft missions. The data span five decades in scales from the inertial range to the electron Larmor radius. In contrast to the inertial range, there is a successive increase toward isotropy between parallel and transverse power at scales below the ion Larmor radius, with isotropy being achieved at the electron Larmor radius. In the context of wave-mediated theories of turbulence, we show that this enhancement in magnetic fluctuations parallel to the local mean background field is qualitatively consistent with the magnetic compressibility signature of kinetic Alfvén wave solutions of the linearized Vlasov equation. More generally, we discuss how these results may arise naturally due to the prominent role of the Hall term at sub-ion Larmor scales. Furthermore, computing higher-order statistics, we show that the full statistical signature of the fluctuations at scales below the ion Larmor radius is that of a single isotropic globally scale-invariant process distinct from the anisotropic statistics of the inertial range. (10.1088/0004-637X/763/1/10)
  DOI : 10.1088/0004-637X/763/1/10
• Electron acceleration to relativistic energies at a strong quasi-parallel shock wave
  
  • Masters A.
  • Stawarz L.
  • Fujimoto M.
  • Schwartz S. J.
  • Sergis N.
  • Thomsen M. F.
  • Retinò Alessandro
  • Hasegawa H.
  • Zieger B.
  • Lewis G. R.
  • Coates A. J.
  • Canu Patrick
  • Dougherty M. K.
  Nature Physics, Nature Publishing Group [2005-....], 2013, 9, pp.164-167. Electrons can be accelerated to ultrarelativistic energies at strong (high Mach number) collisionless shock waves that form when stellar debris rapidly expands after a supernova. Collisionless shock waves also form in the flow of particles from the Sun (the solar wind), and extensive spacecraft observations have established that electron acceleration at these shocks is effectively absent whenever the upstream magnetic field is roughly parallel to the shock-surface normal (quasi-parallel conditions). However, it is unclear whether this magnetic dependence of electron acceleration also applies to the far stronger shocks around young supernova remnants, where local magnetic conditions are poorly understood. Here we present Cassini spacecraft observations of an unusually strong solar system shock wave (Saturn's bow shock) where significant local electron acceleration has been confirmed under quasi-parallel magnetic conditions for the first time, contradicting the established magnetic dependence of electron acceleration at solar system shocks. Furthermore, the acceleration led to electrons at relativistic energies (about megaelectronvolt), comparable to the highest energies ever attributed to shock acceleration in the solar wind. These observations suggest that at high Mach numbers, such as those of young supernova remnant shocks, quasi-parallel shocks become considerably more effective electron accelerators. (10.1038/nphys2541)
  DOI : 10.1038/nphys2541
• Radio frequency current-voltage probe for impedance and power measurements in multi-frequency unmatched loads
  
  • Lafleur Trevor
  • Delattre Pierre-Alexandre
  • Johnson E.V.
  • Booth Jean-Paul
  • Diné Sébastien
  Review of Scientific Instruments, American Institute of Physics, 2013, 84, pp.015001. A broad-band, inline current-voltage probe, with a characteristic impedance of 50 &#937;, is presented for the measurement of voltage and current waveforms, impedance, and power in rf systems. The probe, which uses capacitive and inductive sensors to determine the voltage and current, respectively, can be used for the measurement of single or multi-frequency signals into both matched and unmatched loads, over a frequency range of about 1100 MHz. The probe calibration and impedance/power measurement technique are described in detail, and the calibrated probe results are compared with those obtained from a vector network analyzer and other commercial power meters. Use of the probe is demonstrated with the measurement of power into an unmatched capacitively coupled plasma excited by multi-frequency tailored voltage waveforms. (10.1063/1.4773540)
  DOI : 10.1063/1.4773540
• Collisionless Plasmas in Astrophysics
  
  • Belmont Gérard
  • Grappin Roland
  • Mottez Fabrice
  • Pantellini F.
  • Pelletier Guy
  , 2013 (ISBN: 978-3-527-41074-3).
• First Realistic validation of Plasma Air treatment with QCL-OFCEAS detection of formaldehyde below ppb range
  
  • Guaitella Olivier
  • Barakat Christelle
  • Fasci Eugenio
  • Rousseau Antoine
  • Gorrotxategi Carbajo P.
  • Ventrillard I.
  • Carras M.
  • Maisons G.
  • Romanini D.
  , 2013.