Anna V Nominé

Thomas Gries

Cédric Noel

Alexandre Nominé

Valentin A Milichko

Thierry Belmonte
email: thierry.belmonte@univ-lorraine.fr

A V Nominé

Synthesis of nanomaterials by electrode erosion using discharges in liquids

Keywords: number: 52.80.Wq Discharge in liquids and solids Spark discharges, Submerged discharges, Time-resolved optical emission spectroscopy, Liquid nitrogen

come

to clarify breakdown mechanisms in these media, [START_REF] Ushakov | Impulse breakdown of liquids[END_REF] an undertaking that is still running. Dischargeassisted electrochemistry inspired biologists who dreamt of creating conditions favourable to the emergence of life thanks to these non-faradic electrochemical reactions. [START_REF] Harada | [END_REF] In the middle of the twentieth century, the so-called "electrode effect", i.e. the glow discharge taking place at either electrode depending on the voltage polarity, originally observed in molten salts, was also observed in aqueous solutions. 8 Hickling and Newns 3 discovered that the chemical performances of the GDE processes were several times the values of Faraday's law of electrolysis. The chemical products were different from those of conventional electrolysis such as the release of large amounts of molecular hydrogen and hydrogen peroxide, as well as the H• and OH• radicals in the liquid phase, among other secondary reactions products. 9,10 Very early, erosion of submerged electrodes was noticed. Joseph Priestly described the phenomenon of material erosion by electric spark in 1878. [START_REF] Webzell | [END_REF] But it is only in 1930 that he patented the idea to exploit it for electrical discharge machining (EDM). [START_REF] Gusseff | Method and apparatus for the electrolytic treatment of metals[END_REF] This process is among the earliest non-traditional manufacturing processes. 13 Material removal by controlled erosion through a series of sparks was actually developed in the USSR during the 1940s. B.R. and N.I. Lazarenko first applied it to a machine for stock removal. 14 Full commercial exploitation started only in the 1960s. The production of particles by EDM was considered as a drawback, making the process drift. Powdery by-products were simply not considered as a potential subject of study.

Nairne was the first to make silver and copper wires explode in 1774. 15 Apparently, the production of ultrafine powders by the exploding wire technique was begun by Abrams, [START_REF] Abrams | Production and Analysis of Radioactive Aerosols[END_REF] who studied radioactive Al, U, and Pu aerosols. In this process, the breakage of the wire in the liquid state leads to the disintegration of the wire in the form of hot microdrops or clusters. [START_REF] Sarkisov | [END_REF] Joncich and Reu in 1964 chose to make explosions in liquid nitrogen to produce metal nitrides particles. [START_REF] Joncich | Synthesis of inorganic binary compounds using exploding wire techniques[END_REF] The wires are most often placed in water. [START_REF] Lebedev | [END_REF]20 This allows one to increase the energy deposition and to use the converging shock wave generated by the exploding wire array to reach high pressure on the axis of the implosion. 21 The production of nanoparticles by this process is still the subject of active research. 20,22,23 Explosions in general imply the presence of a plasma phase and a liquid phase.

By comparing the explosive detonation parameters with the carbon phase diagram and performing thermodynamic calculations, it has been shown empirically and theoretically that free carbon in detonation products of powerful condensed carbon-containing individual explosives with a negative oxygen balance should condense in a diamond or liquid phase. 24 This applies to the synthesis of nanodiamonds discovered first in the USSR in 1963 by Volkov, Danilenko and Elin. 25 The situation changed dramatically with the advent of laser-based processes. Maiman constructed in 1960 the first functional pulsed-laser using a ruby crystal and emitting at 694 nm, 26 enabling the development of pulse laser ablation shortly afterwards by Smith and Turner. 27 Patil et al. 28 introduced the concept of pulsed laser ablation at the solid-liquid interface in 1987. This method, known as Liquid Phase Pulsed Laser Ablation (LP-PLA), Laser Ablation in Liquids (LAL) or Laser Ablation Synthesis in Solution (LASiS), is based on the erosion of a solid target immersed in a solution by a laser beam focused through the liquid onto the solid surface. Progress in laser technology led to ultrashort pulse lasers (i.e. picosecond in 1964, 29,30 and femtosecond in 1981, 31 ), which brought considerable perspectives to ablation in liquids by affecting basically the mechanisms of the erosion of solids. 32 Moreover, accessing to synthesis conditions that are far from thermodynamic equilibrium opened up new possibilities, among which the synthesis of nanomaterials with original properties. 33 The idea of using discharges to create conditions for such a purpose made submerged discharges a possible alternative to lasers as sources of energy deposition. The first works mentioning the synthesis of "fine powders" by discharges in liquids were published by Ishibashi et al. 34 They proposed a spark discharge method, in which a spark discharge takes place at the contact points of pellets dipped in the liquid medium. There is no doubt that nanoparticles were also synthesized in these conditions. Dubovoy et al. 35 proposed a similar approach in 1985 as well as Sato et al. 36 in 1992. The first experiments mentioning the synthesis of nanomaterials were certainly carried out by Ishigami et al. 37 in the USA who wanted to propose a simplified synthesis method of carbon nanotubes by arc discharges. This strategy eliminates nearly all of the complex and expensive machinery associated with conventional nanotube growth techniques. Sano et al. 38 in the UK demonstrated one year later with a similar device the possibility of producing carbon onions. The topic starts taking off in 2004, likely because of the ease of running these processes. The idea of resorting to ultrasound assistance came up quickly 39 as it was studied to improve the material removal rate, to reduce the tool wear rate and to optimise the surface quality of the machined workpieces in EDM. 40,41 The idea of locating the plasma outside the discharge came back in 2005 with an alternating current (AC) plasma source. 42 Two years before, generation of nanomaterials by melting and evaporation of a sacrificial metal electrode in a microplasma was demonstrated for the first time by Shimizu et al. 43 .

Next, other solutions of GDE-type processes were proposed. [44][45][46] Corona discharges usually employed in GDE were replaced by various atmospheric plasmas among which Dielectric Barrier Discharges. 47 The concept was applied to surface functionalization of nanoparticles by Mariotti et al. 48 The formation of bubbles (or their pre-existence) in liquids, submitted to an injection current, was the most studied mechanism responsible for breakdown. The formation of pre-breakdown bubbles was examined in the 60s. 49,[START_REF] Balygin | Electric Strength of Liquid Dielectrics[END_REF] Discharges in bubbles injected on purpose in a liquid became a new process as such when Yamabe et al. proposed this concept in 1993. 51 Their utilization for nanoparticles synthesis was proposed only very recently. Shin et al. [START_REF] Shin | [END_REF] produced carbon nanomaterials by discharges in Ar/O2 gas bubbles injected in 1-hexanol and by Yamada et al. who used a gas-liquid slug flow reactor system. 53 A bubble in a liquid can be considered as the counterpart of a droplet in a gas. Both are submitted to discharge-liquid interactions in the processes we are interested in. This is especially the case of warm dense matter produced for instance by ultra-intense laser excitation of liquid fuel targets, [START_REF] Atzeni | The Physics of Inertial Fusion: Beam-Plasma Interaction, Hydrodynamics, Hot Dense Matter[END_REF] a domain far beyond the scope of this paper. The technological importance of introducing droplets in discharges was first foreseen at the end of the 19 th century by Schoop in Switzerland. He used flames to coat metal surfaces with lead and tin, [START_REF] Ulrich Schoop | Elektrolytischer Wasserzersetzungs-Apparat[END_REF][START_REF] Siegmann | [END_REF] inventing plasma spray. The process was strongly improved in the 50s with the development of plasmatron systems. [START_REF] Heimann | Plasma-spray coating: principles and applications[END_REF] Droplets in cold plasmas were studied much later because studies were for long devoted to the sole interaction between electric fields and droplets in neutral gas, after the discovery by Taylor in 1964 of a condition for stable liquid cone existence. [START_REF] Taylor | [END_REF] Hager and Dovichi studied the behaviour of microscopic liquid droplets near a strong electrostatic field and developed the droplet electrospray. 59 The repulsion was a result of droplet charging via the corona discharge created at the tip of the high-voltage probe. In a related study, Kim and Dunn investigated the formation of progeny droplets ejected from a parent droplet, in the presence of an intense electric field. 60 However, investigation of the behaviour of droplets in cold non-equilibrium discharges is much recent.

The interest of using pulsed discharge to improve combustion in engines required to investigate this topic. 61,62 Ward et al. 63 used an ultrasonic nozzle to atomize acrylic acid monomer into an atmospheric pressure glow discharge (APGD) to deposit polymeric coatings, opening up the way to several other processes alike. [64][65][66] Aggregation in low-pressure activated vapour phases leads to ultra-small nanoparticles that can be collected in liquids with extremely low vapour pressure, like oils or ionic liquids, which enables the production of nanofluids. The vapour phase can be created by any physical source like evaporation 67 or sputtering. 68 The liquid can also be directly decomposed by an electron beam for instance. 69 We shall not forget to mention that liquid droplets in discharges are also of primary importance for climatologists 70 and astrophysicists. 71 For example, the atmosphere of Titan is constantly bombarded by galactic cosmic rays and Saturnian magnetospheric electrons causing the formation of free electrons and primary ions, which are then stabilized by ion cluster formation and charging of aerosols. These charged particles accumulate in drops in cloud regions of the troposphere.

The synthesis of nanoparticles by discharge-liquid interactions involves the presence of all states of matter. The previous history, because it is intentionally short, is biased to select some seminal works as breakthroughs in the topic. The following paper is centred on nanoparticles synthesis. Other correlated aspects can be found with much detail in some reference review papers:

• For discharges in liquids, see Refs. 72-83 • For discharges in contact with liquids, see Refs. 48, 84-96 • For laser in liquids, see Refs. 32,33, 97-104 • For spray in discharges, see Refs. 105-110

• For electrical wire explosion, see Refs. 22, 111 • For sputtering onto liquids, see Refs. 112-117 or in books, see Refs. 6, 81, 118-122. To the best of our knowledge, no review or book is available about nanoparticles synthesis for discharges in bubbles.

II. EXPERIMENTAL APPROACHES

There are two main ways to generate a direct current (DC) discharge in a liquid to synthesise nanoparticles. A low-voltage generator (typically below 1 kV) can be used, which requires to put the two electrodes in contact, the discharge being created when the gap between the electrodes increases, either by material erosion or by mechanical separation. The current 123,124 or the voltage 125 are, most often, used as closed-loop feedback signals to keep the inter-electrode gap distance, measured with a detector, as constant as possible. For instance, Bera et al. 126 developed an optoelectronic system made of three main components to maintain constant a preselected distance between the electrodes (Fig. 1): a photosensor for an optical emission diagnosis, a feedback loop which is composed of analog electronics and a computing unit, and a servo-unit for axis translation of the anode. Resorting to high-voltage discharges (typically beyond 1 kV) enables breakdown at short gaps (i.e. from 100 µm to 1 mm typically, depending to the strength of the dielectric liquid). If one electrode is insulating, AC or radiofrequency (RF) excitation can also be used. 72 Experimentally, reactor designs are many (Fig. 2), even though it is relatively simple to set up a process, as only a power generator, two electrodes and a vessel containing the liquid are required. Indeed, the power source can be used to generate a DC, AC (low frequency), radio frequency, microwave, unipolar pulsed discharge or bipolar pulsed discharge. The electrodes can be a pin, a plate or even materials granules. The shape of the vessel can be a cylinder or a cone to prevent the drift of particles from the discharge zone. The process can be assisted by ultrasound. 151,152 One electrode can vibrate to facilitate electrical breakdown and mix the treated solution. [START_REF] Hamdan | [END_REF] (c) supercritical or cryogenic fluids, 148,155 (d) multi-pin-toplate, [START_REF] Vanraes | Electrical discharge in water treatment technology for micropollutant decomposition. Plasma science and technology-progress in physical states and chemical reactions[END_REF] (e) pulsed arc with vibrating electrode, 156 (f) brush-to-plate, [START_REF] Vanraes | Electrical discharge in water treatment technology for micropollutant decomposition. Plasma science and technology-progress in physical states and chemical reactions[END_REF] (g) microwave discharge with quarter-wavelength antenna, 75 (h) microwave discharge with coaxial line, 75 (i) diaphragm discharge, 157 (j) spark discharges between metal chunks, 158 (k) funnel-type reactor, 159 (l) capillary discharge. 160 Discharges created within liquids are often close to the equilibrium, 161

III. A SIMPLIFIED SKETCH OF EVENTS

The fastest voltage rise time used to generate discharges in liquids is around 150 picoseconds. 162 For the synthesis of nanoparticles, it is always beyond 2 nanoseconds (i.e. within the so-called short time scale) and can extend to the second time scale. This means that typical phenomena found in target ablation by femtosecond and even picosecond lasers (leading to the so-called ultra-short interactions) are beyond the scope of this work.

The synthesis of nanoparticles by electrode erosion using discharges in liquids can be described by the following simplified chronological sketch of events (Fig. 3). 163 Before plasma ignition, charges are injected into the liquid during the pre-breakdown phase. Either electrons are injected at the cathode or impurities (and even the liquid itself if it is aprotic) can be split into ion pairs of opposite charges at the anode. Charge injection drives the electrohydrodynamic movement of the liquid, possibly leading to turbulent flows. Once charges are injected in the liquid, breakdown occurs when the liquid density becomes low enough to enable charge multiplications. When the discharge ignites, if the dissipated power is high enough, a shockwave can be emitted. The current flows back and forth across the electrode-liquid interfaces. Both electrodes are submitted to these current oscillations that can be described by a damped RLC circuit. The discharge volume extends towards the ground electrode under the electric field, forming an ionized channel that crosses the inter-electrode gap, typically in several nanoseconds for sub-millimetric gap distances, a bit like streamers in dielectric barrier discharges. When the channel reaches the ground electrode, the voltage drops suddenly, transforming the discharge channel into an arc discharge. The electrodes, on contact with the hightemperature arc discharge, start melting after a few hundred nanoseconds, releasing a metallic vapour in the discharge. Once the metallic vapour is emitted, nanoparticles are synthesized by condensation. Other erosion processes may also be involved, depending on the conditions, leading to different size distributions. Produced particles are next transferred into the liquid phase. Once the discharge phase stops, the gas volume it occupied expands and collapses several times within the millisecond time scale, leading to bubble oscillations. An energy balance is useful at this stage. 164 Basically, the largest part of the energy delivered by the source is used to create the plasma (~95%). Electrode erosion consumes only about one percent of the total energy (~1%), whereas the bubble dissipates a few percent (~5 %). The amount of energy spent in the shockwave is almost negligible (<0.1%). Of course, these figures are affected by the type of discharge used, the inter-electrode gap distance, the nature of the liquid, etc. However, we readily understand that the overall yield of nanoparticles production is a key issue in this kind of process, even though large quantities of nanoparticles can be produced because of high production rates (typically ~100 mg h -1 and up to 10 g h -1 ). 89

IV. ELECTRODE EROSION

Electrodes can be eroded by different mechanisms. This produces impacts with complex shapes. As mentioned previously, some recently identified ablation mechanisms in laser processes with ultra-short time scale (from a few femtoseconds to a few picoseconds) will not be considered. 165,166

A. Erosion mechanisms

At short pulses, thermal processes must be considered, which is not true for ultra-short processes. The melting of the electrodes by the discharge is the most important mechanism that leads to particles synthesis. The modelling of the process is generally attuned to those developed for the description of arc-surface interaction in domains like welding or circuit-breakers. [167][168][169] However, several other mechanisms are possible, depending on experimental conditions. Whatever the conditions, the mean energy of ions in these high-pressure discharges is too low to promote any sputtering effect as it is too often claimed.

Melting and vaporization (evaporation and sublimation)

The discharge pressure at ignition is typically beyond hundreds of bars. It is still around a few tens of bars 200-250 ns after ignition, at the end of a nanosecond-pulsed discharge. 170 This means that the energy of most particles in the ionized gas is certainly too low for sputtering. The melting of metal electrodes is not due to the Joule effect either, the resistivity of the material being too low until dissolved gases start forming bubbles, which affects the material conductivity. Following Hamdan et al. 163 , electrodes melt essentially by irradiation from the blackbody emission of the high-temperature discharge. The maximum temperature reaches a fraction of eV (typically around 5000 K), discharges being close to local thermodynamic equilibrium. As radiation is emitted similarly on either side of the discharge channel, erosion spots are similar on either electrode. The time to melting is of the order of 100 ns for a micrometric erosion spot on an aluminium electrode.

Local melting of the electrode leads to the formation of a liquid well. The bottom of the well is in contact with a mushy region if the electrode is an alloy and not a simple metal. The interface between the surface of the liquid well and the discharge is not clearly described for several reasons: the time evolution of the pressure is not well known, the existence of a plasma sheath is contentious as it should be extremely thin, the surface is likely to be unstable, etc.

The heated electrode emits a vapour either by evaporation from the molten electrode or by sublimation from the solid electrode, and this vapour condensates in colder parts where nanoparticles form. This mechanism together with nucleation (described hereinafter) is responsible for the synthesis of the smallest particles generated by discharge in liquids. Vaporization at the surface of the molten pool can be estimated as follows. 171 The rate of the atomic flux (m -2 s -1 ) leaving the surface during normal evaporation and due to a pressure gradient is given by the Knudsen-Langmuir equation:

( ) P P RT M m s p - = ∞ • π α 2 , ( 1
)
M is the molar mass of the evaporating molecule or atom. The factor α (<1 and depends on temperature) is called the efficiency coefficient. R is the ideal gas constant. P is the vapour phase pressure and ∞ s P the saturation pressure at the temperature T of the liquid surface. These quantities are related by the Rankine form of the Clausius-Clapeyron equation, assuming the latent heat of vaporization Lv constant:

                - - = ∞ b v s T T R L P P 1 1 exp 0 , (2)
P0 is the ambient pressure and Tb is the equilibrium boiling temperature at the ambient pressure (the normal boiling temperature). Lv is either the heat of sublimation or the heat of evaporation. If Lv cannot be considered as constant, other expressions (like Dupré or Riedel's formulas) are required.

Vaporization can also occur because of a concentration gradient. 172 The diffusive vaporization rate d m

• is expressed in terms of a phase change at the surface and the subsequent transport of the vaporized species to the bulk gas phase through the mass transfer boundary layer surrounding the pool. The vaporization rate is then defined as:

b s d RT MP K m ∞ • = , ( 3
)
where K is the mass transfer coefficient given by similarity laws.

Spallation

What is referred to as "spallation" corresponds to material ejection driven by relaxation of the stress induced by the treatment. 98,173 Basically, the spallation of a liquid phase in contact with its mother solid phase is due to a drift in the phase diagram where the following thermodynamic pathway crosses the liquid branch of a binodal domain and enters a metastable region. This leads to the nucleation, growth, and coalescence of voids. This induces the formation of a transient foamy structure of interconnected liquid regions, and eventual separation (or spallation) of a thin liquid layer from the bulk of the target.

Material spallation occurs at the liquid-solid boundary, where the spall strength is lower than that in the solid phase. However, the tensile wave is strong enough to cause the cut-off in the liquid phase, according to the applied criteria of nucleation (Fig. 4a).

In a liquid -the situation under vacuum being different -, the particles produced by this mechanism exhibit distribution sizes that are around 10 nm. 165 The foamy structure coarsens with time and eventually decomposes into individual droplets on the time scale of nanoseconds. The top liquid layer loses stability because of Rayleigh-Taylor instability and decomposes into large droplets, estimated to have diameters from hundreds of nanometers to tens of micrometers. 174

Phase explosion, also called explosive boiling

What is referred to as "phase explosion" or "explosive boiling" is a process where a superheated surface region, expanding rapidly, undergoes a fast decomposition into a mixture of vapour and liquid droplets (Fig. 4b). Thermodynamically, it means that solid matter is rapidly superheated up to the thermodynamic critical temperature, at which the spinodal decomposition in vapour and liquid phase in the irradiated material occurs by homogeneous nucleation.

This is typically what can be observed when a very clean glass of water (i.e. without surface defects to promote nucleation of steam bubbles) is heated in a microwave oven without boiling. The superheated metastable liquid can undergo an explosive liquid-vapor phase transition when a massive homogeneous nucleation of vapor bubbles starts, by simply touching the surface of the liquid, for instance. This phenomenon appears at much higher energy than spallation and occurs within the nanosecond time scale. 175 As for spallation, on the time scale of nanoseconds, vapor, small nanoparticles (~10 nm), and large liquid droplets are produced.

Droplet emission due to surface temperature gradients

Non-homogeneous distribution of surface temperature is caused by radial temperature gradients in the discharge. It can also be amplified if the discharge channel is not oriented perpendicularly to the liquid well. Because of this non-homogeneous distribution, two driving forces responsible for expulsion of the molten liquid are generated: one emanating from the recoil pressure effect caused by evaporation, and the other one from the Marangoni effect. Large (beyond 100 nm in diameters and up to hundreds of µm)

droplets are produced by both mechanisms (sometimes called hydrodynamic splashing 176 ) and emitted radially (Fig. 4c).

Recoil pressure gradient

Expulsion of large liquid droplets can be driven at sufficiently high power by the vapour recoil. This phenomenon induces a pressure that splashes the melt and residual recoil pressure redistributes the metallic melt. 100 The melt expulsion by the recoil pressure is caused by the spatial variation of the normal stress exerted by vapour pressure within the hot spot. The gradient of vapour pressure makes the melt flow outwards, from the high-pressure spot centre to the low-pressure periphery. 177

Marangoni's effect

The Marangoni or thermocapillary effect creates a strong tangential stress on a surface submitted to large gradients. The liquid movement within the molten pool is affected by this boundary condition that is due to the dependence of the surface tension on temperature. The surface tension of liquid pure metals decreases with increasing temperature 178 and then, Marangoni's force acts in the opposite direction of the surface temperature gradient. As the temperature usually peaks at the centre of the well, the liquid flows from the centre to the spot edge, which contributes to the expulsion of the melt. At high heating rates and small thicknesses of the molten layer, Marangoni's stresses can lead to the rupture of the molten layer, 179 melt spattering, and formation of droplets. 180

Bursting bubbles

The emission of droplets from the molten well is possible via bursting of bubbles formed in the liquid.

Bubbles within the molten well can be created from boiling of the liquid metal or from absorption of discharge gases. 181 Blisters are thus formed on the surface of the well. 182 A liquid film, covering the surface of the bubble, bursts and releases very fine droplets. 183 The subsequent collapse of the bubble forms a rising jet of liquid that pinches, according to the Rayleigh-Plateau instability, and produces a spray of micron-scale droplets within tens of microseconds. (Fig. 4d). The effect of the electric field on the ejection of droplets from bubbles bursting at the liquid surface was considered by Holgate and

Coppins. 184

B. Shapes of eroded areas

The different shapes of eroded areas, produced by the mechanisms described previously, are depicted in

C. Limited erosion

It is possible to strongly limit electrode erosion by using materials with high melting points (usually tungsten) and low current density. 134,135 When nanoparticles are synthesized from a liquid precursor, this point is critical to keep the as-produced nanofluid with high purity.

V. THERMODYNAMICS

When the size of a particle shrinks, its thermodynamic properties changes. The Gibbs free energy of nanoparticles, made of one element, is determined by the bulk free energy and the surface free energy for both the solid and liquid phase. 188 Surface atoms are dominant for the size effect on the thermodynamic properties of nanoparticles. On the basis of the Gibbs free energy, thermodynamic properties of nanoparticles, such as melting temperature, molar heat of fusion, molar entropy of fusion, and temperature dependences of entropy and specific heat capacity can be determined.

Thermodynamics of alloy nanoparticles is still under debate, as recent results showed the existence of multiple local energy minima submitted to local thermal fluctuations. 189 The way phases separate during ultrafast cooling is difficult to predict. In simple systems like, for instance, two-element alloys undergoing eutectic decomposition with a wide miscibility gap, several situations can be encountered:

• The boundaries of phase compositions as well as transition temperatures are changed. It is the case of the Sn-In system (Fig. 5a). 190

• The decomposition mechanism changes. In the Cu-Ag system, phase separation at the nanoscale occurs above 5 nm particle diameters in the 15 and 30 at.% Ag composition range. 191 The separation into Cu-rich and Ag-rich domains takes place by spinodal decomposition. Phase separation no longer occurs at nucleation sites but through the whole particle at once (Fig. 5b).

• The decomposition follows a specific path that leads to metastable phases that cannot be described by the two previous approaches. For instance, supercooled Au-Si nanoparticles exhibit a unique metastable phase δ1 which grows heteroepitaxially to Au. 192 Mater. 29, 1131 (1993). Copyright 1993 Elsevier. (Fig. 1 in Ref. 193) Then, resorting to metastable phase diagrams is a first approach, which cannot be considered as highly predictable. As an illustration, there are more than 20 metastable Au-Si crystalline phases that have been identified to date.

Eutectic and spinodal decompositions are not easy to distinguish in nanostructures. The Cahn-Hilliard approach predicts the critical wavelength λc over which fluctuations of compositions occurs in spinodal decomposition. This approach is based on the following equation:

4 4 " 0 2 2 2 ~x f D x D t i ij i ij i ∂ ∂ - ∂ ∂ = ∂ ∂ ρ κ ρ ρ , ( 4
)
where ρi is the volume density of atoms i in a binary mixture made of atoms i and j, ij D ~ is the interdiffusion coefficient, x is the spatial coordinate, t is the time, " 0 f is the second derivative by composition of the Helmholz free energy of a unit volume of homogeneous material (See Fig. 6a).

κ is the gradient energy coefficient. When the system is phase-separated into high and low composition regions, kT

2 ξ κ =
. ξ is a length that scales with the interface width. It depends on the crystal lattice. For an fcc structure, ( )

T T T a c c - = 2 ξ : 1 2 1 1 2 -         - = c T T ka κ . ( 5
)
Tc is the maximum temperature of the miscibility gap. a is the lattice parameter. It is possible to show that the following solution:

( ) ( )

      = - x t A t i λ π ρ ρ 2 sin (6)
satisfies Eq.( 4), which leads by substitution to the so-called amplification factor R:

( )

2 " 0 2 2 1 ln h f h D dt t A d R ij         + - = = κ , (7)
Here, λ

π 2 = h
and A is the amplitude of concentration fluctuations with a wavelength of λ. ρ is the average composition. The sign of R is given by the evolution of the amplitude of a concentration fluctuation with a given wavelength: if it increases, then R > 0. The critical wavelength λc corresponds to the condition R = 0. It reads:

" 0 2 8 f c κπ λ - = . ( 8
)
In Fig. 6b, λc is plotted versus composition. 191 The width of the miscibility gap is about 93% and the minimum of λc is about 6 nm at 30 at.% of Ag. With increasing Ag content, λc increases more and more rapidly and tends to infinity close to the spinodal composition ( Consequently, we deduce from Fig. 6b that, below 5 nm typically, the nanoparticle cannot adopt the modulated structure of spinodal decomposition. This is observed experimentally: particles smaller than 5 nm in diameter grow as a solid solution whatever the composition (Fig. 6c). This shows that spinodal decomposition is both size and composition dependent. It also demonstrates that non-equilibrium alloys exist at sufficiently small scale.

0 " 0 = f ).

VI. PARTICLES AGGREGATION AND STABILITY OF DISPERSIONS

Proper conditions are sometimes encountered for particles to aggregate, which depends on physicochemical processes at stake. Then, particles may settle, as they get bigger, on the bottom of the vessel in which they are contained. Because water is often used as primary solvent, the specific case of oxidative processes in liquid has to be specifically emphasized. Indeed, the stability of solutions is a key point in the synthesis of nanofluids. It is usually estimated by resorting to the zeta potential that derives of the electric double layer that forms around particles.

A. Aggregation of nanoparticles

Aggregation of nanoparticles, which leads to the formation of highly porous structures, is of primary importance as it generates multi-modal distributions that can compromise the application of synthesis processes of nanoparticles. 194 Irreversible aggregation of nanoparticles is well understood. The fractal formation in colloidal aggregation and the general underlying mechanism have been extensively investigated. They are well known and they can be appropriately applied on the nanoscale. Conversely, reversible aggregation still deserves to be deepened and it is only touched upon in this paper.

Irreversible network formation

Irreversible aggregation will be described briefly since many reference works provide comprehensive pictures of this phenomenon. [195][196][197][198] The formation of colloidal structure and the irreversible nanoparticle aggregation are well described by the diffusion-limited cluster aggregation (DLCA) model, 199,200 a more accurate version of the diffusionlimited aggregation (DLA) model. 201 The latter considers the existence of a primary particle acting as a seed for the formation of a cluster, to which other individual particles diffuse and attach. The former describes how all particles bind to each other, and diffuse as clusters that can also attach together when in contact.

To account for the compactness of agglomerated structure, often estimated less dense than observed, the reaction-limited cluster aggregation (RLCA) model was proposed. 202,203 It relies on the idea that approaching particles aggregate with a probability less than one. This probability is defined by the distribution of situations where the strength of the repulsion between the particles is inferior to the strength of the attraction force at shorter distances, which enables aggregation. The repulsion force is given by the Derjaguin-Landau-Verwey-Overbeek (DLVO) theory. 204,205 It accounts for van der Waals interactions and electrostatic forces arising between the double layers around the particles. This supplementary degree of freedom favours non-binding collisions between clusters, which increases their interpenetration, and produces denser structure once binding occurs. This approach was used for instance by Ziashahabi et al. 206 to study the formation mechanism of bead-chain-like ZnO nanostructures prepared via DC arc discharge in liquid.

In colloidal solutions, an appropriate change of the conditions destabilizes the nanoparticle suspension, and induces the aggregation of particles. The particles interact differently because their surface state changed in the new situation. 207 Forces exerted on a single particle are modified and the former balance is affected in such a way that aggregation is enabled. This results in chemical bonding, non-covalent bonding being conventionally restricted to self-assembly.

Reversible network formation

Aggregation of nanoparticles is expected to proceed reversibly, and the general macroscopic laws describing the irreversible network formation process are not directly useable on the nanoscale. Reversible network formation is a complex emerging topic, an example of which is as follows.

It is possible to account for the many-body van der Waals interactions on the nanoscale with so-called "many-body dispersion (MBD)" models. 208,209 Van der Waals forces between polarizable non-metallic nanostructures can be understood by collective interactions between wavelike charge density fluctuations. A simple summation over pairwise interactions between instantaneous particle-or fragment-like dipolar fluctuations is not as meaningful. The collective wavelike fluctuations are responsible for the emergence of nontrivial modifications of the power laws that govern non-covalent interactions.

Practically, there are several ways in which the modeling of the formation of nanoparticle networks is approached. A practical review on reversible network formation is available in Ref. 207). Nonetheless, a theoretical framework accounting for the reversibility of the clustering process is still pending.

B. Oxidation in liquids

Discharges are often created in liquids with different oxidation capabilities. Notoriously, water is more oxidizing than ethanol, and ethanol than ethylene glycol. 210 In general, alcohols, glycols, and hydrazine can all serve as reductant source. The choice of the liquid is also known to affect the shape of the nanoparticles. 86,145,211 Discharges in liquids make these effects even more complicated to explain than other processes as liquids are transformed by discharges. Indeed, by-products resulting from liquid-discharge interactions are likely to play a role on the nucleation and growth of nanoparticles.

In the case of water, where air gases dissolve, chemistry is extremely complicated. It includes as main species H2O, O2, N2, OH•, H•, O•, N•, H2O2, NOx, O3, HNO2, HNO3, NO•, ONOOH, and ions. 212 Species responsible for oxidation are likely to be several, but the way oxidation occurs is still unclear.

Simple description explaining how ethanol molecules scavenge OH radicals, and generate reducing species (H and H2): 213,214 OH + C2H5OH  H2O+ C2H4OH (a)

eaq -+ C2H5OH  H + C2H5O - (b) H + C2H5OH  H2+ C2H4OH (c)
is put forward to explain how non-oxidized metal nanoparticles can be produced by discharges in ethanol. Similar reasoning can be applied to ethylene glycol as metallic nanoparticles can be produced as well (see e.g. Refs. 215 and 216).

More complex mechanisms are certainly at stake, as in the case of carbonaceous liquids or liquids containing carbonaceous additives like surfactants, where nanoparticles are often coated by a few layers of carbon. 86 Lee et al. 217 showed by dissolving WCl6 in ethanol that they could synthesize spherical tungsten nanoparticles. By adding an anionic surfactant (SDS) at 30% and 50% (the surfactant/WCl6 molar ratio), spherical particles are no longer agglomerated. A weak fraction of triangular and hexagonal-shaped nanoparticles are also found. The presence of carbon is not studied but is very likely.

Generally, the balance between a too large quantity of surfactant, leading to carbon contamination and a too small quantity, leading to no effect, can only be determined by a trial and error procedure.

C. Stability of dispersions

Because of the importance of nanofluids in many industrial applications, it is essential to ensure a high stability of the dispersion in order to maintain the fluid properties constant in time. The physicochemistry of nanoparticles in liquids is a complex topic, and the examples of gold and silver will be taken to illustrate this aspect. Next, the concept of electrical double layer will be introduced to explain how the zeta-potential is defined, a useful quantity to evaluate the stability of dispersions.

Physicochemistry: examples of Ag and Au

The synthesis of gold nanoparticles by femtosecond laser ablation in aqueous solutions produces colloidal Au nanoparticles that are characterized by their surface hydroxylation. 218 Au-O compounds, followed by a proton loss to give surface Au-O -, results in the negative charging of the nanoparticles.

This increase in the net surface charge of particles limits their coalescence, due to electrostatic repulsion.

If n-propylamine is added to the solution, amine groups react with the nanoparticle surface. This is accompanied by the reduction of the particle size, which indicates that functionalization of Au nanoparticles is likely achieved during their formation by laser ablation.

On the contrary, for silver 219 (but also for iron 220 for instance), surface oxidation affects the stability of the dispersion. Nanoparticles synthesized in solutions can be stabilized by adding sodium dodecyl sulfate (SDS) as a surfactant. The surfactant coverage and the charge state on the nanoparticle surface are closely related to the concentration of SDS, which finally determines the stability of the nanoparticles in the solutions. The nanoparticles tend to be aggregated when the coverage is less than unity, while they are very stable when the surface is covered with a double layer of the surfactant molecules. High stability can be achieved as well by using isopropanol, in absence of additional capping agents. The liquid undergoes a self-inhibited free radical chain mechanism, generating the organic stabilizing agent that covers particles and stabilize them. 221 The stability of the suspension can also be affected by the evolution of the solvent itself that may age or react with chemicals from the gas environment it is in contact with. Another possibility showed recently is water splitting during pulsed laser ablation of metal targets in water. 222,223

Electric double layer (EDL)

The electrical double layer on the surface of a nanoparticle is based on the Gouy-Chapman-Stern model explaining how a surface charge, created by dissociative ionization of surface molecules, is counterbalanced by ions in the electrolyte. Basically, the double layer is made of the Stern (inner) layer where the ions are strongly bound and the Gouy-Chapman diffuse (outer) layer where the ions are loosely bound. Within the diffuse layer, there is a virtual boundary below which ions and particles behave as a whole. This boundary, placed at the zeta potential, is known as the surface hydrodynamic shear or slipping plane. Consequently, when a particle moves in an electrolyte, ions below the slipping plane move along, contrary to ions beyond the plane.

Because a nanoparticle is surrounded by the electrical double layer that sizes one Debye length, its movement can be affected by electric fields. 224 It turns out that the EDL can be polarized under the influence of external fields. This plays an important role in the accurate evaluation of the electrophoretic motion of particles. 225 There are three main polarisation mechanisms: EDL, volume and field polarisations.

The EDL polarisability of charged species can be calculated by the Maxwell-Wagner theory from the permittivities and conductivities of the species and the medium. For a colloid, the concept of surface conductivity must be introduced and its polarizability is given by the Maxwell-Wagner-O'Konski theory. 226,227 EDL polarization cannot be neglected when the Debye length is comparable to the colloidal size, the zeta potential is moderately high, and/or the applied frequency of the external electric field is of the same order as the rate of ionic diffusion over the colloidal size.

In volume polarisation (Fig. 7), the positive and negative ions in the liquid are migrating under the influence of the external field. As these ions cannot penetrate the colloid, the volume of the particle acts as a physical obstacle, which hinders the migration. The subsequent accumulation of charges of opposite signs on either side of the particle induces a dipole moment anti-parallel to the external electric field. Field polarisation is due to the Coulombic attraction/repulsion between a charged colloidal particle and the clouds of surrounding ions. Counter-ions/co-ions tend to accumulate in the tail of the particle as they circumvent the particle, resulting in a dipole moment proportional to the external electric field. This effect is observed for highly charged particles with high zeta potential, the induced dipole moment being then proportional to the external electric field.

VII. PARTICLES' LIFE

Particles' life starts by nucleation, nuclei becomes seeds, and seeds nanocrystals. What is meant by seeds is an intermediate object between a nucleus and a nanocrystal, in which no structure fluctuation (i.e.

changes of bond-orientational order) 228 is possible. From their formation to their transfer into the liquid phase, nanoparticles are submitted to complex environments where huge gradients prevail. Even though discharges in liquids are often close to equilibrium, nanoparticles can adopt complex shapes.

A. Shape control

The shape and structure of Au or Ag nanocrystals strongly affect their spectroscopic properties (like localized surface plasmon resonance, 229 surface-enhanced Raman scattering, 230 Surface Enhanced Fluorescence 231 , plasmon enhanced optically stimulated luminescence, 232 etc.), enabling Enhanced Spectroscopies, reviewed for example in the following contributions. [233][234][235]

From atoms to nuclei

Classical nucleation theory stipulates that the formation of nuclei in supersaturated homogeneous solution is governed by the balance between the bulk and surface energy of the new phase.

        ∆ -         - = T k G T k E A J B ex B a nucl exp exp , (9)
The first exponent (-Ea / kBT) is related to the kinetic barriers with an overall activation energy Ea, while the second exponent (-∆Gex / kBT) represents the thermodynamic barrier. The parameter A is a preexponential factor that depends on the properties of the investigated material. It is possible to predict nucleation rates from this equation for any material at a given level of supersaturation.

However, as discussed by Gebauer et al. 236 , values calculated accordingly can differ by orders of magnitude from experimentally measured data. Non-classical nucleation theories are then needed but are beyond the scope of this work.

From nuclei to seeds

The shape of a single crystal at equilibrium in an inert gas at low temperature (rigorously 0 K) or vacuum can be determined by the Wulff construction. For a fcc metal, possible shapes correspond to polyhedrons (truncated octahedrons). However, the final shape adopted by a fcc nanocrystal can differ a lot from Wulff's shapes. This deviation can be attributed to several reasons: 210 -nucleation and/or growth are far from equilibrium, -surface energies are modified by a capping agent, impurity, or solvent, -twin defects affect nucleation and growth and lead to new shapes such as decahedron and icosahedron with a total free energy lower than that of Wulff's polyhedrons, -the synthesis temperature is high.

There are three main ways to control the shape of nanoparticles.

Interfacial energy

Discharges in liquids are often near thermodynamic equilibrium. The shape of nanoparticles is then controlled by the interfacial free energy, γ. It can be defined as the energy needed for creating a unit area of new surface:

s b P T n N S G ερ γ 2 1 , , =       ∂ ∂ = , ( 10
)
where G is the free energy and S the surface area. For a given surface, Nb is the number of broken bonds, ε the bond strength, and ρs the atom surface density. For an fcc structure with a lattice constant of a, one finds easily that:

2 110 2 100 2 111 36 . 3 4 24 . 4 a a a ε γ ε γ ε γ × = < × = < × = , (11)
Then, single-crystal seed takes an octahedral or tetrahedral shape in order to maximize the expression of {111} facets and minimize the total surface energy. These inequalities can be modified by the introduction of twin defects. The strain energy caused by twin defects greatly increases as the seed grows in size. This critical dependence on size is also sensitive to reaction kinetics, which offers a possibility of control.

Kinetics limitation

If the reaction that forms the atoms building the particle is slow, nuclei and seeds form through random hexagonal close packing (rhcp), together with the inclusion of stacking faults. 237 This limitation by a kinetically controlled step makes seeds adopt shapes that depart from those favoured by thermodynamics. Inclusion of stacking faults and/or twin planes can thus lead to the formation of platelike seeds.

The key is to ensure an extremely low concentration of metal atoms, so the nuclei will not be able to grow auto-catalytically into polyhedral structures. Instead, the atoms will add to the edges of a planar cluster to generate a plate-like seed. This was observed (Refs. 130 and 202) by forming plate-like lead particles using discharges in liquid nitrogen (Fig. 8). is also what Saito et al. 79 observed in similar experiments (Fig. 9).

Xia et al. 210 give a more accurate view of oxidative etching. The distribution of single-crystal versus twinned seeds can be further modified by oxidative etching, in which zero-valent metal atoms are oxidized back to ions. If a ligand for the metal ion is also present in the same solution, a combination of the ligand and O2 can result in a powerful etchant for both the nuclei and seeds.

Consequently, we can readily understand how discharges in colloidal dispersions contribute to the shaping of nanoparticles. As the discharge runs, production of oxidizing species (among which H2O2)

progressively modifies the liquid composition by reacting with ligands and by forming etching agents.

However, these mechanisms have not been clarified yet in the case of discharges in liquids.

From seeds to nanocrystals

The condensation mechanism for the growth of nuclei, i.e., deposition of atoms or molecules on the surface of growing particles, the condensation rate can be approximated by the collision rate of monomers with the nucleus: 239

( ) 2 0 n th cond r r v v + = ρ , ( 12
)
The average thermal velocity of the monomers in the gas is given by

( ) 2 / 1 8 M T k v B th π =
, where ρ is the monomer density, r0 and rn are the hard sphere equivalent radii of a monomer and a cluster. In this expression for the condensation rate, the cross section is taken as the surface area of a sphere of radius rn and the sticking probability upon collision is equal to 1.

Consequently, if this probability is not 1 but dependent on facet orientation, the shape of a crystal can change if there is a mix of different facets on the surface. For example, when metal atoms add to the {100} faces of a nanocube, they migrate to the edges of the face resulting in the elongation of the {111} facets. Progressively, the cube transforms into a cuboctahedron and eventually an octahedron. This shape evolution process is known as overgrowth.

The introduction of a capping agent can alter the growth rates of the different facets, thus dramatically altering the final nanocrystal shape. In discharges in liquids, sodium citrate, 131,138 reverse micelle solution (sodium bis (2-ethylhexyl) sulfosuccinate in dodecane), 240 sodium dodecyl sulfonate, 135 cetyltrimethylammonium bromide (CTAB), 241 among many other capping agents, are used, from recipes developed for chemical processes.

It is important to mention that surfactants and capping agents are very difficult to decompose completely, which results in poor catalytic activities. Discharges in liquids are often presented as efficient methods to synthesize nanoparticles without any surfactant.

B. Other phenomena affecting size

Nanoparticles within the discharge can melt, be vaporized or coalesce, leading to changes in size distributions.

Melting

The melting point of a spherically symmetric nanoparticle decreases as the radius decreases, which is only observed at very low radius (typically below 10 nm). 242 This effect can be described down to about 1 nm for gold 242 or nickel 243 by using a continuum model (Fig. 10). Melting point depression is often accounted for by the generalized Gibbs-Thomson implicit relation:

( ) a l s l s sl m m m m m p m m m p p T T T T T c T T L -         - = +         - +         ∆ +         - ρ ρ κ ρ σ 1 1 2 ln 1 * * * , ( 13
)
Lm is the latent heat [J kg -1 ], Tm the temperature at which the phase change occurs

[K], * m T the bulk phase change temperature [K], cp is the specific heat [J kg -1 K -1 ], ∆cp = cl -cs, σ the surface tension [N m -1 ]
and κ the mean curvature. ρ is the density [kg m -3 ]. p is the pressure [Pa]. Subscripts a, s and l indicate ambient, solid and liquid, respectively.

Surface atoms are less strongly bound to the cluster than bulk atom. Then, melting proceeds by separation of surface atoms from the bulk. This separation is paid for with the latent heat. With a sufficiently large cluster, the energy required is relatively constant since each surface molecule is affected by the same quantity of bulk molecules. However, as the cluster decreases in size, surface molecules are less and less influenced by the bulk, and less and less energy is required for separation.

The change in the ratio of surface to bulk energy may also lead to a structural transition and a reduction in surface tension.

Vaporization

Once the vapour is emitted from the electrode, it condensates farther in colder regions into droplets that are still prone to vaporization along their path. Vaporization of a nanodroplet is affected by the particle size. The equilibrium pressure at the surface of the droplet (Eq. 2) must be corrected. 244,245 The correction is:

        × = ∞ r RT M P P s r s 1 2 exp ρ σ , ( 14
)
r s P is the saturation pressure for a droplet of radius r. σ the surface energy and ρ the density.

Coalescence

Coalescence affects the size distribution of nanoparticles by merging primary particles, i.e. those produced by condensation, into a larger single nanoparticle (see Fig. 11). This is a very complex phenomenon between nanocrystals that can be decomposed into elementary steps: 246 1) Coagulation refers to the attraction and contact between primary particles. In the discharge, it is due to in-flight collisions after which particles remains stuck together. This phenomenon is often referred to as Smoluchowski ripening. 247 2) Rigid body reorientation occurs once the primary particles are close enough to start interacting. An interface forms that is usually misoriented, because of the randomness in the coagulation geometry.

Thus, a number of interface layers are shorn. Eventually, the particles rotate as rigid bodies, driven by mutual torques, maximizing the contact area. A more coherent interface with a larger area is then produced.

3) Formation of defects results from the very fast (few ps) reorientation stage. If primary particles cannot be perfectly aligned because of their geometry, which is the most probable situation, a twin boundarycommonly accompanied by the appearance of misfit dislocations -forms, as rotation is hindered by the newly formed bonds.

4)

Heat release results from free-surface annihilation due to bond formation. As the surface/volume ratio decreases, the part of the surface energy remaining after the creation of the interface transforms into thermal energy. 248 5) Temporary melting of interface has been proposed by Grammatikopoulos et al. 249 If the heat release essentially contributes to the temperature rise of the child particle (case of an adiabatic system), it can be sufficient to reorganize it at the atomic level. A wave that originates at the interface between the primary particles, propagates along the body of the child particle, resulting in its full crystallization.

Concurrently, a temporary disordering of the atoms at the interface and a "softening" of their bonds is observed.

6)

Neck growth is all the more fast as local curvature is high, which speeds up atomic surface diffusion.

A neck is formed on either side of the interface when primary particles coagulate. However, the neck only reaches substantial dimensions if it is promoted by the temporary melting of the interface that enhances the transport of atoms from the core of primary particles to the neck.

7) Plastic deformation. Even after neck formation, the system can still behave within the elastic limit (~0.8 nm after Averback and Zhu 250 ) as a coupled-oscillator (with possible rotation of the two primary particles). Thermal vibration makes each particle pulsate around the initial contact locus. Plastic deformation is then possible by relaxation of accumulated internal shear stress, possibly assisted by the heat, which enables the glide of misfit dislocations and their annihilation at the surface. Protrusions can thus be formed.

8) Consolidation.

The neck thickens at a growth rate that slows down, because the difference in curvature between the neck and the primary particles becomes smaller. The two-particle system gets stiffer, which hinders its pulsation, and makes it behave as a single particle. Consolidation takes about 100-200 ps. Free-surface annihilation becomes negligible, and so the heating it induces. Oftentimes, experimentally grown nanoparticles are trapped in a metastable steady state due to quenching upon deposition. Coalescence is often considered as achieved at this point in theoretical approaches. 9) Slow ageing, however, extends the coalescence process and favours a "sphericization" of the system, as referred to by Lewis et al. 251 , a process that lasts for hundreds of ns up to a few μs. This means that the ratio of the surface area of the particle tends toward the surface area of a sphere of the same volume as the particle. Here, the neck acts as a sink for diffusing atoms.

Oxidation in air

Except for gold, nanoparticles produced by submerged discharges are prone to be oxidized either in air, when they are floating on or taken out of a non-oxidizing liquid. They can also be oxidized on purpose but this aspect is deemed to be beyond the scope of this article. Oxidation of sub-micrometric objects has been studied thoroughly. 252 Only oxidation in dry air at room temperature is described here. This process changes the size of the pristine nanoparticles because of the change of molar volume between the metal and its oxides.

Basically, when the particle is less than 1 µm in diameter, electrical neutrality, assumed in Wagner's theory for thick films, 253 within the oxide layer is not satisfied, and when less than 20 nm, the Nernst-Einstein relationship is no longer appropriate. Thus, theories of thin film growth must consider atom jumps in the presence of large electric fields and the possibility of large space charges.

Cabrera and Mott's theory assumes that electrons can freely pass from the metal to ionize adsorbed oxygen atoms or molecules at the oxide/gas interface, so that the electron electrochemical potential (Fermi level) is equal in the metal and the adsorbed gas layer. 254 The general picture is based on the injection of defects into the oxide at one of the film interfaces (with gas or metal). This process is assumed to be rate-controlling. Using the Kröger-Vink notation, and taking the example of chromium oxidation, one notices that, depending on the temperature and O2 partial pressure, either triply charged chromium vacancies are injected at the gas/oxide interface and diffuse inward:

• + + → h V O O Cr x O 6 2 3 2 3 ' ' ' 2
(p-type Cr2-xO3) -gas/oxide interface or triply charged interstitial chromium atoms are created at the metal/oxide interface and diffuse outward: (n-type Cr2O3-x) -metal/oxide interface In both cases, the accumulation of vacancies at the metal/oxide interface induces the development of Kirkendall porosity. Such a porosity is not observed if the oxide growth is driven by two other types of mechanisms, observed with other metals: oxygen atoms, injected at the gas/oxide interface, diffuse inward (p-type MO1+x) and oxygen vacancies, injected at the metal/oxide interface, diffuse outward (ntype M1+xO). This explains why oxidation of some nanoparticles leads to core-shell structures and/or hollow particles.

For very thin (i.e. where tunnelling can no longer be neglected) films, Fromhold showed that the assumption of electronic equilibrium is reasonable. [START_REF] Fromhold | Theory of metal oxidation[END_REF] In the case of nanoparticles, quantum confinement effects should also intervenes, but to best of our knowledge, no study has been done on this topic yet.

Finally, it is also important to mention that reduction of nanoparticles after oxidation does not lead to a complete recovery of the initial metallic state. For instance, after oxidation in O2, ZnCu clusters (< 10 nm) transform into a polycrystalline cluster consisting of separate CuO and ZnO nanocrystals.

Subsequent reduction in H2 converts CuO into Cu with ZnO nanocrystal covering its surface. Then, H2 dissociates onto metallic Cu to form H atoms that partially reduce ZnO into CuZn. 256 These theories are based on point defects but oxidation is also affected by the presence of structural defects that play a tremendous role at low temperature, as volume diffusion is strongly limited because of high-energy barriers of jump processes. These barriers can be lowered by stress to a certain extent, which is also an important effect to take into account in order to predict correctly how oxidation occurs in these specific conditions. Finally, in air, not only oxygen but also water can intervene as oxidation agents. The reader is referred to a number of seminal papers illustrating these phenomena in the case of different elements: e.g. oxidation at room temperature of reduced iron oxide nanoparticles by Wang et al. 257 , air oxidation of monodisperse cobalt nanoparticles by Varon et al. 258 , the first stage of air oxidation by Chen et al. 259 or the evolution of the oxide thickness and oxidation rate of Si nanoparticles in air and at room temperature by Yang et al. 260 The effect of the spherical geometry on metal oxidation driven by mobile charged point defect was also studied by Fromhold. 261 The electric potential developed across the oxide is the same as for planar geometry, but the time, for complete oxidation for spherical particles of a given diameter, can be as much as a factor of 3 shorter than that for planar samples of the same thickness.

C. More complex structures

Discharges in liquids can be used to synthesize nanoparticles made of, at least, two different materials that can be organized in different structures depicted in Fig. 12. Both materials are made immiscible, either because thermodynamics imposes it (Fig. 12a to e), or because they are synthesized successively (Figs. 12a,b or c).

Bi-particles (Fig. 12a) have been synthesized by different groups: e.g. Pt-Pd, 262 Ni-Cu, 241 Ag-Pt, 263 Ag-Cu, 264 etc. This kind of structures results from coalescence of particles that are synthesized either simultaneously (from electrodes made of different materials) or separately (in space and/or time).

Tsukanov et al. 265 showed by molecular dynamics that immiscible nanoparticles can interpenetrate if they collide with sufficiently high relative velocities.

Core-shell structures (Fig. 12b) are usually obtained by chained treatments where seed nanoparticles of one material are coated by another material in a second step. The best way to coat nanoparticles already present in a liquid is to decompose a precursor that deposits onto the nanoparticles surface. This is easy to achieve by discharges in contact with liquids for instance.

Conversely, discharges in liquids rather form attached nanoparticles in sequenced treatment. The expansion of the discharge repels nanoparticles added in the liquid from emitted species by the electrodes. Then, one-step synthesis is only possible by mixing intimately the two elements in a crucible: e.g. CuO@Ta2O5. 123 For a carbon shell, one-step synthesis is also possible with a carbonaceous liquid, like benzene or hexane: e.g. (Co, Ni, Fe)@C. [266][267][268] Then, using a liquid precursor (like HAuCl4) in a aqueous colloidal dispersion containing nanoparticles, say Pt, should give Pt@Au nanoparticles in one step. But to the best of our knowledge, no experiment of that kind has ever been done yet.

Pulsed laser ablation also provides similar capabilities. For instance, the introduction of ad hoc solvents lead to the creation of carbon shell or structures. 269 Even more interestingly, the possibility of using CO2 sub-products in solution with NaOH was showed to be efficient for the production of gold-carbon nanocomposites. 270,271 Kabbara et al. 272 showed that it is also possible in two steps, by wrapping one type of particle (Cu) into an ultra-thin foil made of another one (Zn) synthesized beforehand, this mechanical way of forming core-shell nanoparticles requiring no deposition. Decorated nanoparticles (Fig. 12c) are usually obtained by depositing one phase of small size onto another phase of large size. Usually, this is a two-step process. This has been achieved in one step by discharges in liquids for the Si-Sn system. 273 The mechanism relies on grain detachment of a sintered Si-Sn electrode, where Si and Sn are well separated, leading to tin nanoparticles attached to silicon grains.

Nanoparticles with multiple domains (Fig. 12d) are either made of one phase inlaid into another or made of one phase segregated at the grain boundaries of another. These particles, as Janus particles (Fig. 12e) -a special type of nanoparticles with surfaces exhibiting two or more distinct physical properties -have only been synthesized by PLAL but not by discharges in liquids.

VIII. PARTICLES' SIZE DISTRIBUTION

As discharges in liquids most often produce several size distributions, it is essential to describe how to measure them, which requires to master the techniques used in order to get convergent results. It is also possible to change size distributions by fragmenting nanoparticles.

A. Size distributions

Ashkarran et al. 130 (Table 1) obtained almost monodisperse (i.e. with a geometric standard deviation of a log-normal distribution < 1.25) gold colloidal dispersion (mean size: 8 ± 3 nm) by creating discharges between two titanium electrodes in a HAuCl4 solution. This process behaves as an electrochemical process, similarly to discharges in contact with liquids. Submerged discharges in liquids where nanoparticles are produced by electrode erosion produce multimodal size distributions (Fig. 13). 145,274,275 At most, three distributions coexist. Large particles (i.e. > 100 nm and up to tens of µm) are only a few and spread all over the substrate. Nanoparticles of intermediate size (i.e. from 20 to 100 nm) are found either as necklace-like or fractal agglomerated structures or as isolated particles. Nanoparticles of small size (i.e. from 2 to 10 nm) are primary particles.

However, this approach is not satisfactory as the various populations are not associated with their mechanisms of formation. From the list of erosion mechanisms established hereinabove, it is easy to understand that the origin of the largest population of nanoparticles is due to ejection of droplets from the liquid well but there are several possible mechanisms that may be involved. The origin of the smallest distribution can only be due to the rapid condensation of the vapour phase. On the contrary, the origin of the intermediate-size population can be multiple. The smallest nanoparticles can continue growing if they stay within the discharge region, reacting with the emitted vapour that is not homogeneously distributed within the interelectrode gap. They can also grow by coarsening with other particles and then, the distribution of the particles themselves also matters. Consequently, we could expect 4, 5 or even more contributions in a multimodal size distribution. This statement is corroborated by several works where trimodal distributions are proposed (see for instance Ref. 128 and in Ref. 150 Table 1) but with very different mode positions. The issue lies in the possibility for distinguishing accurately all the modes in a distribution of nanoparticles. In LASiS, nanosecond pulsed laser ablation is known to give more or less narrow monomodal distributions that are described by a log-normal distribution, ablation being principally driven by thermal processes. On the other hand, pico and femtosecond ablation may lead to multimodal distributions, due to the multiple ablation mechanisms. In this latter case, size distributions can also be sometimes by described by a log-normal distribution. Therefore, caution must be taken to interpret these size distributions.

Logically, the existence of several erosion mechanisms must lead to an identical number of distributions.

However, one reads commonly that two populations are to be distinguished: small and large. In recent studies based on molecular dynamics simulation, Shih et al. 175,176 investigate the origin of bimodal nanoparticle size distributions produced by laser ablation in liquids. They establish that "three distinct nanoparticle generation mechanisms operating at different stages of the ablation process and in different parts of the emerging cavitation bubble" exist. Two of them produce large nanoparticles (spallation and explosive boiling) and one small nanoparticles (nucleation at the very front of the emerging cavitation bubble). They conclude that "the coexistence of the three distinct mechanisms of the nanoparticle formation at the initial stage of the ablation process can be related to the broad nanoparticle size distributions commonly observed in nanosecond PLAL experiments". To the best of our knowledge, no attempt was made to fit with a multimodal approach, size distributions of nanoparticles produced by ns laser pulses that can be described by a monomodal log-normal function.

The confusion also exists because only the smallest nanoparticles are collected for analysis. For instance, it is possible to proceed as follows. 276 In a discharge-based process, when heavy particles settle at the bottom of the vessel that contains the synthesized colloid suspension, this latter can be cleared from its sediments by pouring carefully the upper part of the solution into another beaker. Then, small quantities of the new solution containing only its lightest particles can be collected with a pipette for further characterization.

There are several techniques that can be used to evaluate the mean diameter of a set of well-defined nanospheres: DLS (dynamic light scattering), AFM (atomic force microscopy) in tapping mode, UV-Visible spectroscopy and two-colour SPR (Surface Plasmon Resonance) spectroscopy. Several papers 277,278 showed that estimates of the mean diameter of well-controlled nanospheres synthesized by wet chemistry and evaluated by the previously-cited techniques are very close, demonstrating the convergence of these methods for almost ideal spheres. On the other hand, for other types of nanoobjects, caution must be taken to interpret results produced by these different diagnostics, as discussed by Tomaszewska et al. 279 For example, differences between TEM and UV-Vis spectroscopy often arise due to:

-possible agglomeration of nanoparticles during deposition on the grid, especially in colloidal dispersions in water; -non-spherical shape of nanoparticles; -a non-representative sampling of nanoparticles for TEM analysis; -nanoparticles are oxidized, which affects the effective dielectric constant to be used in the Mie fit required by UV-Vis spectroscopy.

The existence of an intermediate size distribution is not clearly understood yet, likely because it probably results from distinct phenomena:

-the duration of the nucleation process may strongly vary if the vapour is emitted at the very beginning of the process, when the pressure is high, or much later, when it is reduced; -it is also true for the vapours emitted by the cathode and by the anode; -coalescence of particles certainly affects the size distribution, and smear out transitions between modes, etc.

To conclude about size distribution, it is important to mention that even though it is a useful concept to describe the population of nanoparticles produced by a given process, it cannot be used usually to correlate observations on size modes with underlying growth mechanisms.

B. Fragmentation of nanoparticles

Until now, fragmentation of nanoparticles has not been observed with electric discharges, but this process might exist as well for sufficiently fast electric fields (picosecond-pulsed discharges). Indeed, Giammanco et al. 280 fragmented Au nanoparticles by interaction with the second and the third harmonics of a Nd:YAG picosecond laser. However, laser and discharge treatments are already associated in that purpose. 281 The photothermal mechanism responsible for nanoparticle fragmentation relies on the way temperature rise in response to the laser excitation. At moderate fluences, heating, melting, and evaporation of irradiated nanoparticles enables size reduction of primary nanoparticles, creating a new population of smaller ones by nucleation and growth from the evaporated atoms. At higher fluences, nanoparticles can be either fully evaporated or fragmented by explosive decomposition into vapor and liquid droplets when superheated up to the limit of thermodynamic stability of the molten material. 282 According to Huang and Zhigilei,283 the detailed analysis of the nanoparticle fragmentation mechanisms reveals two distinct pathways of the formation of the fragmentation products: (1) the direct generation within ∼100-200 ps of relatively large nanoparticles with diameters ranging from ∼2.5 to 4 nm and (2) the much more gradual growth of smaller nanoparticles formed by agglomeration and coalescence of atomic clusters on the timescale of tens of nanoseconds. These two pathways are responsible for the early appearance of a bimodal mass-weighted particle size distribution.

IX. CROSSING INTERFACES

A. The discharge-liquid barrier

The discharge-liquid interface is very difficult to probe but interesting results could be obtained thanks to original experimental works. Sano et al. 284 recorded under microgravity conditions the development of an arc discharge submerged in water. They notice, from images of the reactor falling freely in a column under vacuum for about 4 s, that a transparent gas film was generated in the gas phase zone adjacent to the gas-liquid interface. The authors concluded that chemical reactions involving water were possible inside the bubble and not only at the water-gas interface.

The behaviour of nanoparticles in discharge bubbles is much better understood thanks to recent results obtained with synchrotron facilities on pulsed laser ablation in liquid. [285][286][287] It has been shown, with SAXS experiments, that the bubble interface is a quite impervious boundary for the primary and secondary particles. In PLAL, they cross the interface barrier by jet-like emission after the collapse of the first rebound. These two statements must also apply for submerged discharges created between two electrodes. Indeed, bubbles oscillations produce many different kinds of instabilities like Rayleigh-Taylor, Kelvin-Helmholtz, Landau-Darrieus or Birkhoff, leading to jet-like emission (Fig. 14). 288 Basically, the discharge-liquid barrier is rather nanoparticles tight. The transfer from one phase to the other is limited. This is also the case of the liquid-air barrier.

B. The liquid-air barrier and safety aspects

Sano et al. 38,289 synthesized carbon nanomaterials by erosion of graphite electrodes submerged in water.

They noticed that nanoonions (average diameter : 25-30 nm) were majorly present as floating powder on the water surface, the other nanocarbons being rather found located on the bottom of the vessel.

Even though the liquid evaporates, this mechanism does not transport any nanoparticles in air. 290 This is why the dewetting of a colloidal dispersion concentrates nanoparticles at the moving front of the liquid fingers that form upon demixing of the two phases. 291 However, evaporation is affected by the presence of nanoparticles. 292 As long as nanoparticles are left in the liquid, nanofluids in containers can be considered as safe.

However, they also face public concerns about their safety in use as they can be spilled by accident. The lack of legal framework for their handling and recycling is a key issue. 293 It would be prudent to pursue green designs by choosing nontoxic or biodegradable nanoparticles. 294

X. POST-TREATMENT WITH DISCHARGES AND ELECTRIC FIELD

Discharges in liquids or in contact with liquids can be used for appropriate post-treatments of colloidal dispersions. Sometimes, one-dimensional objects are also synthesized together with nanoparticles as their assembly can be driven by electric fields. Then, applying electric fileds with proper strengths in a post-treatment is a way to form new one-dimensional objects.

A. Post-treatment with discharges

Surface functionalization by atmospheric plasmas in contact with colloids

Because of the difficulty to control the complex chemistry induced by the presence of additives in processes using submerged discharges, it is possible to resort to a plasma post-treatment, as developed by Mariotti et al. 48 It relies on the utilization of plasma jets in contact with liquids. In this configuration, electrons are injected into the liquid solution from the plasma phase. More efficiently than radicals formed in the gas phase and ultra-violet radiation, electrons are supposed to mainly contribute to a specific liquid chemistry for the surface engineering of pre-synthesized nanomaterials. This post-processing results in considerable improvements in the stability of a colloidal dispersion, especially in the case of silicon nanocrystals in ethanol or water, which could not be observed with standard electrochemistry, 295 and carbon. 286 Stable Si-ncs can be created with different surface functional groups. This physicochemical change affects the optical properties of the colloid, its photoluminescence emission varying from 590 nm up to 650 nm just by changing the oxygen-bonding configuration at the surface. 297

Post-treatments by discharges in liquids

Discharges in liquids can be used to activate physicochemical processes, for example by enhancing the reactivity of metal particles. High-voltage electrical discharges in liquids 74 are capable of:

-modifying the kinetics of reaction between particles; -decreasing the energy of formation of new phases and compounds; -increasing surface and bulk diffusion of reacting compounds; -modifying the morphology and the dimensional characteristics; -changing physical and mechanical properties of treated materials.

If these effects have been demonstrated, they are not very well documented and still to be investigated in detail. Nonetheless, these empirical post-treatments demonstrate clearly that discharges in liquids are also very promising for controlled modifications of nanoparticles and not only for their synthesis.

B. Role of electric fields on particles

A particle moving in a liquid environment can be submitted to many contributions like gravity, van der Waals or multi-body interactions, electrothermal flows, Brownian motion, fluid dynamics, etc. In this work, because light is shed on discharges produced between submerged electrodes, a special attention is paid to the electrokinetic motion of particles lying in an electrolyte exposed to an externally applied electric field. Electro-osmosis exerts a long-range interaction and does not scale with the particle radius, contrary to dielectrophoresis (see next paragraph) that results from short-range interactions and scales with the particle volume. Both forces can be combined to trap and manipulate nanoparticles. 298 The principal drawback with DC electro-osmosis is the need for high DC voltages (up to 1 kV across a 1 cm gap typically to reach acceptable flow velocities (around 100 μm s -1 ). 299 This can lead to unwanted electrochemical reactions, electrode dissolution, gas bubble production and/or hydrodynamic instability.

Electro-osmosis

Resorting to AC electro-osmosis with voltages of only a few volts by using planar periodic arrays of interdigitated, asymmetric electrodes can overcome these drawbacks.

It was shown by Gierhart et al. 300 that particle motion by AC electro-osmosis greatly influences the morphology of resultant nanoparticle assemblies. 301 With frequencies ranging typically from 0.1 to 10 kHz, this phenomenon is likely to be considered in AC-excited submerged discharges (see Table 1).

Electrophoresis (EP)

The motion of a submicron particle is influenced by the cloud of ions near the charged or polarised particle induced by the external electric field.

The particle can be initially charged and possess an equilibrium surface charge, which need not be uniformly distributed over its surface. The interaction between the electric field and the surface charge corresponds to the linear EP problem, where the dynamic reaction experienced by the particle and the induced mobility are linear in the external electric field. For a polarisable uncharged particle, the external electric field produces an additional space charge distribution in the electrolyte which decays away from the particle. The corresponding induced-charge electrophoresis (ICEP) problem implies that the force and torque experienced by the particle are nonlinear (quadratic) in the electric field.

Dielectrophoresis (DEP)

Dielectrophoresis can be utilized to assemble nanoparticles into wires. Alternating current fields allow assembly of particles without the interference of direct current effects like electro-osmotic and electrochemical effects. The dielectrophoretic force results from the interaction of particle dipoles in a non-uniform DC or AC field. It is conservative and oriented along the field gradient:

( ) 2 3 2 rms m cons dep E a F ∇ ⋅ ℜ = ω ε π (15)
a is the particle radius,. ℜ denotes the real part of the Clausius-Mossotti factor, which is an effective polarisability of the particle in the medium, and it depends on ω, the angular field frequency. It is expressed by the complex relative permittivities of the particle p ε ~ and of the medium m ε ~: (16) ℜ is bounded between -1/2 and 1 and varies with the frequency of the applied field and the complex permittivity of the medium.

( ) ( )         + - ℜ = ℜ m p m p e ε ω ε ε ω ε 2 ~~
( ) ( ) ( ) ( ) ( ) ( )( ) 1 2 2 0 2 2 2 1 1 2 3 2 - ∞ ∞ + ℜ - ℜ + ℜ = + + - +         + - = ℜ MW MW m p MW m p p m m p m p τ ω τ ω σ σ τ σ ε σ ε ε ω ε ε ω ε (17) At high frequency, when m p m p MW ε ε σ σ τ ω 2 2 1 + + = >> -
where σm denotes the electric conductivity, 302

( ) ( )         + - = ℜ → ℜ ∞ m p m p ε ω ε ε ω ε 2 . At low frequency, ( ) ( )         + - = ℜ → ℜ m p m p σ ω σ σ ω σ 2 0
Positive dielectrophoresis occurs when ℜ> 0 ; the force is in the direction of increasing field strength and the particles are attracted by the electrode edges. The opposite behavior is observed for negative values.

In the interelectrode gap, nanoparticles experience a large electric field that induces a dipole in each of them if they are not conductors. In the case of most metallic particles for instance, ℜ =1 because the conduction band electrons redistribute continually to annul the internal electric field. The field frequency needed to keep charges separated on each side of a metallic nanoparticle should excess tens of THz. On the contrary, dielectric nanoparticles with diameters between 10 nm and 100 nm can be trap by an ac field in the MHz range. 303 However, in 2001, Hermanson et al. 304 succeeded in synthesizing gold wires with millimeter length in an AC electric field from a suspension of gold nanoparticles by dielectrophoresis. They suggested some collective effect or complex electrohydrodynamic interactions to account for direct chaining of gold nanoparticles. Barsotti et al. 305 considered in addition to the dipole contribution to the DEP force, effects of higherorder polarisation moments and forces. 306 These higher-order corrections to the DEP force cannot be neglected for conductive particles in AC fields and introduce a non-conservative term proportional to the curl of the vector electromechanical potential n W and the imaginary part of Clausius-Mossotti factor.

Dipolophoresis (DPP)

In the general case where spatially non-uniform electric fields are encountered, both DEP and ICEP contribute to the particle movement. The combined non-linear electrokinetic effect is known as dipolophoresis. 307 DEP and ICEP mechanisms act in different directions and their combination can either be additive or subtractive. Under appropriate circumstances like a constant gradient field, the sum of the two contributions can be null and particle mobility turned to zero.

XI. OUTLOOK AND NEW APPLICATIONS

In this pedagogical review, emphasis is put on specificities of discharges in liquids in nanoparticles' synthesis. As colloids and nanofluids are generated by these processes, application related to these biphasic systems is extremely large. For instance, they can be used:

-to replace water and ethylene glycol as conventional coolants in an automotive car radiator; 308 -to produce hydrogen via their use as working fluid in solar parabolic trough collectors; 309 -to recover oil more efficiently in rock pores by depositing a wedge film making rock surfaces more dewettable; 310 -to decrease by several orders of magnitude bacterial activity including Gram-negative, Grampositive and spore bacteria; 311 -to serve as a new reaction media in the preparation of chemicals, like amides via direct amidation of aliphatic carboxylic acids; 312 -to enable sensing like, for instance, those magnetically polarizable nanoemulsions used for imaging internal defects in materials; 313 etc.

There is a number of challenges that deserve to be tackled: -Recently, nanoalloys, i.e. single solid solution, made of immiscible elements have been synthesized by PLAL. 314,315 Producing them by discharges in liquids has been showed to be possible too. 316 Mechanisms leading to nanoalloying of immiscible elements are still to be clarified.

These new results open up the way to the study of optical properties of metallic nanoalloys. 315,317 Engineering of nanoalloys for optical properties consists in controlling defects, element distributions, oxidation and surface functionalization processes.

Fig. 1 :

 1 Fig. 1: Block diagram showing the concept and major components of the optoelectronic feedback control of the electrode feeding system that is implemented in the present ADS process. Reproduced from D. Bera, E. Brinley, S. C. Kuiry, M. McCutchen, S. Seal, H. Heinrich and B. Kabes, Rev. Sci. Instrum. 76, 033903 (2005), with the permission of AIP Publishing (Fig. 3 in Ref. 126).

153

 153

Fig. 2 :

 2 Fig. 2: Selected examples of submerged discharge processes. a) pin-to-pin spark discharge with ultrasound assistance, 152 (b) pin-to-plate, 154 (c) supercritical or cryogenic fluids, 148,155 (d) multi-pin-to-

  which takes only a few nanoseconds to reach. Orders of magnitude of key parameters of this kind of discharges are: -Electron density ne = 10 18 -10 19 cm -3 ; -Initial pressure P0 = ~100-1000 bars; -Initial temperatures Tr ~ Tv ~ Te < 1 eV; -Current I ~ 1-10 A.

Fig. 3 :

 3 Fig. 3: Simplified sketch of events occurring in submerged discharge processes.

Fig. 4 .

 4 Fig. 4.Of course, these mechanisms can combine over sufficiently large time scales, and give more complex structures. However, most of the eroded areas exhibit these shapes that are the most commonly encountered. Roughly, the size of a crater is a linear function of the square root of the electrical charge dissipated during the discharge process.185

Fig. 4 :

 4 Fig. 4: Typical shapes of impacts and related underlying mechanisms leading to the production of large particles. Pictures in a) and b) are reproduced with permission from Phys. Procedia 86, 66 (2017). Copyright 2017 Elsevier. (Fig. 1 in Ref. 186) and from Mater. Res. Bull. 48, 3268 (2013). Copyright 2013 Elsevier. (Fig. 4 in Ref. 187). Other pictures are current authors'credits. All spot sizes are around 10-20 µm in diameter.

Fig. 5 :

 5 Fig. 5: Examples of metastable phase diagrams. a) Full range metastable phase diagram of Sn/In nanoparticles comparing with the equilibrium phase diagram. Reproduced with permission from Acta Mater. 123, 82 (2017). Copyright 2017 Elsevier. (Fig. 11 in Ref. 191). b) Calculated chemical and coherent spinodal boundaries in the Ag-Cu system. Reproduced with permission from Script. Metall.

Fig. 6 :

 6 Fig. 6: (a) Homogeneous part of the Gibbs free energy f0 and (b) critical wavelength λc in spinodal decomposition versus composition in Cu-Ag system (T = 873 K, Tc = 1053 K, a = 0.4 nm). In (b) full dots mark experimental critical size, the open circle marks the estimated lower value of the critical size. (c) HRTEM images of Cu-Ag nanoparticles in 30 at.% Ag particles and corresponding FFT diffraction patterns. Areas rich in silver or copper are denoted by Ag or Cu. As they contain also Cu or Ag, the

Fig. 7 :

 7 Fig. 7: Illustration of the volume polarisation effect created by an accumulation of negative ions on one side of the particle with negative surface charge, while positive ions are depleted in this region. The induced dipole moment points from the front side toward the back side will be formed and is anti-parallel to the external electric field. Reproduced with permission from Molec. Phys. 113, 2511 (2015). Copyrights Francis & Taylor. (Fig. 12 in Ref. 225).

Fig. 8 :

 8 Fig. 8: Left: Diffraction pattern of a PbO2 single crystal superimposed on the TEM bright field micrograph. Oxidation occurred after air exposure but growth corresponds to Pb in liquid nitrogen. The arrows show directions that are perpendicular to the growth of the hexagon faces. The zone axis is [001]. Reproduced with permission from Particuology 40, 152 (2018). Copyright 2018 Elsevier. (Fig. 7 in Ref. 148). Middle and right: Schematic representation of the radial growth of the lead nanosheet. Reproducedwith permission from J. Appl. Cryst.52, 304 (2019).Copyright 2019 IUCr. (Fig. 8 in Ref. 238)

Fig. 9 :

 9 Fig. 9: HRTEM images of gold nanoparticles synthesized by discharges in reverse micelle solution. The synthesized gold nanoparticles have different crystal features like fcc single-crystalline particles, multiply twinned particles, and incomplete multiply twinned particles (single-nanotwinned fcc configuration). Reproduced with permission from Jap. J. Appl. Phys. 57, 0102A4 (2018). Copyright 2018 The Japan Society of Applied Physics -IOP Science. (Fig. 7 in Ref. 79).

Fig. 10 :

 10 Fig. 10: Experimental data and Gibbs-Thomson curves established from Eq. (1) showing the size dependence of the melting temperature of gold 242 and nickel 243 nanoparticles in an (Y -1/X) graph. The behaviour is almost linear beyond 2 nm in this graph and follows a reciprocal law.

Fig. 11 :

 11 Fig. 11: Parameters affecting nanoparticle coalescence. (a) Effect of temperature and relative orientation on the coalescence of 2 crystalline Pd nanoparticles, 3 nm in diameter. (b) Effect of temperature and size on the coalescence of 2 crystalline Pd nanoparticles of different sizes (5 and 3 nm in diameter), with touching (100)-type facets. Note that (b) is not to scale with (a). (For clarity, temperatures are represented by shades of colours). P. Grammatikopoulos, C. Cassidy, V. Singh and M. Sowwan, Sci. Rep. 4, 5779 (2014); licensed under a Creative Commons Attribution (CC BY) license. (Fig. 2 in Ref. 249).

Fig. 12 :

 12 Fig. 12: Possible structures associated with bi-particles. a) bi-particles (more or less interpenetrated); b) core-shell; c) decorated; d) multiple domains; e) Janus.

Fig. 13 :

 13 Fig. 13: Size distributions of particles produced by discharges in liquid. Three distributions are commonly found. The two lowest distributions can overlap, leading to two main distributions, as it is described in some other processes like pulsed laser ablation in liquid. Reproduced with permission from Mater. Chem. Phys. 142, 199 (2013). Copyright Elsevier. (Fig. 4 in Ref. 274).

Fig. 14 :

 14 Fig. 14: Example of jet-like emission after the collapse of the first rebound of a bubble formed by nanosecond-pulsed discharge in heptane. Jets fragment into groups of aligned bubbles forming threads that point toward the center of the interelectrode gap.

  Electro-osmosis is responsible for the flow driven by the Coulomb force induced by an electric field on the net charge in the electrical double layer at the electrodes. Upon experiencing an electrical potential gradient, the mobile ions in the Gouy-Chapman layer move past the electrode surface and drag molecules of the solution, putting the fluid in motion. Nanoparticles, whatever their charge, are carried by the flow of the liquid and not because of their own electrical double layer. The electro-osmotic flow is evaluated by solving the Navier-Stokes continuity and momentum equations coupled with Laplace's equation for the external electric field and Poisson's equation for the potential within the electric double layer.

  of the square of the root-mean-square electrical field. Then, the strength of the force exerted on a particle depends on 2 rms E ∇ and also on ℜ. The frequency dependence of ℜ is given by:

  Nanoparticles -MPs: Microparticles -NWs: Nanowires -NSs: Nanosheets -NCs: Nanocubes -NRs: Nanorods ** 0, 1+1 or 2 means no electrode, 1 electrode made of a given material and 1 electrode made of another material or 2 electrodes made of the same materials.

Table 1

 1 (see at the end of the manuscript) presents selected examples of discharges in liquids for the synthesis of nanoparticles of simple metals. It is useful to compare conditions and have practical values of parameters. For carbon materials, a thorough review is available.

Table 1 :

 1 Submerged discharges in liquids used to synthesize nano-objects.

DATA AVAILABLITY STATEMENTS FOR AUTHOR USE

Data available on request from the authors

-Chemistry driven by discharges in liquids is a very complex topic, still in its infancy. Discharges in carbonaceous liquids for instance are responsible for the deposit of a carbon layer onto nanoparticles, but they also produce unusual chemical by-products like polyynes. 318 This is also true for the synthesis of uncommon nitrogen polymer 319 or diamondoids. 320 Chemistry driven by discharges in liquids is extremely complex and capable of producing rare byproducts. The possibility of a better control of chemical reactions in the liquid phase is an important stake. It is a major challenge for discharges in contact with liquids 321 but also for discharges in liquids.

-Multimodal size distributions are still an issue for large scale production of well-controlled nanoparticles but huge progress made recently on erosion mechanisms gives promising outlook to achieve unimodal monodisperse size distributions.

The possibility of controlling size distributions in a one-step process is not an easy task. If unimodal distributions have been produced by several groups (see Table 1), none of them achieved to get monodisperse distributions. Discharges in liquids containing precursors are more suited to this goal than processes based on electrodes erosion.

-Control of defects to control properties.

Engineering of defects (twin defects, vacancies, edges, doping, etc.) in nanoparticles is a key topic to better control optical 322 and catalytic properties. 323 Controlling defects in nanoparticles by discharges in liquids remains hitherto a dare to take. This is could be a way to develop low-cost catalysts based on non-rare chemical elements or to tune non-linear effects in nanoparticles for optics.

-Control of surface physicochemical states.

Controlling the grafting of selected chemical functional groups on nanoparticles remains a challenge. It is essential anyway to set the stability of nanoparticles in colloidal dispersion. It also control their activity, not only catalytically but also optically. It may even affect the chemical homogeneity of very small nanoparticles. As seen before, discharge-enhanced post-treatments can be used to functionalize nanoparticles, but resorting to these processes cannot be generalized to any situation.

ACKNOWLEDMENTS

The authors acknowledge the French PIA (programme d'investissements d'avenir) project Lorraine Université d'Excellence (Ref. ANR-15-IDEX-04-LUE) for financial support.