Second Harmonic Scattering of Molecular Aggregates

A general model is developed to describe the polarization-resolved second harmonic scattering (SHS) response from a liquid solution of molecular aggregates. In particular, the molecular spatial order is introduced to consider the coherent contribution, also known as the retarded contribution, besides the incoherent contribution. The model is based on the description of a liquid suspension of molecular dyes represented by point-like nonlinear dipoles, locally excited by the fundamental field and radiating at the harmonic frequency. It is shown that for a non-centrosymmetrical spatial arrangement of the nonlinear dipoles, the SHS response is very similar to the purely incoherent response, and is of electric dipole origin. However, for centrosymmetrical or close to centrosymmetrical spatial arrangements of the nonlinear dipoles, the near cancellation of the incoherent contribution due to the inversion symmetry rule allows the observation of the coherent contribution of the SHS response, also known as the electric quadrupole contribution. This model is illustrated with experimental data obtained for aqueous solutions of the dye Crystal Violet (CV) in the presence of sodium dodecyl sulfate (SDS) and mixed water-methanol solutions of the dye 4-(4– dihexadecylaminostyryl)-N-methylpyridinium iodide (DiA), a cationic amphiphilic probe molecule with a strong first hyper-polarizability; both CV and DiA form molecular aggregates in these conditions. The quantitative determination of a retardation parameter opens a window into the spatial arrangements of the dyes in the aggregates, despite the small nanoscale dimensions of the latter.


Introduction
Molecular aggregates are ubiquitous in nature. They may be found on many occasions, especially in the field of soft matter, where structures like micelles, liposomes, or vesicles are often encountered, and can be used as nanoscale probes for nonlinear optics alongside other nanoparticles [1,2]. The organization of the molecules at the microscopic level in these aggregates defines the properties or the function of these nanostructures. Hence, developing new tools and techniques to investigate how molecules arranges at the nanoscale nanostructures is always welcome, but is often hindered by the dimensions of the structures, often of sizes much smaller than the wavelength of light. Hence, one reverts to methods where the wavelength of the probe is of the order of the size of the aggregates, using X-ray or neutron scattering, for instance [3,4]. As a result, linear optical techniques are of limited use because of the diffraction limit rule. Nonlinear optical techniques, however, may offer new strategies in these studies, especially for even order techniques, like second harmonic generation (SHG), obeying the inversion symmetry cancellation rule. This rule states that SHG is forbidden in a medium possessing inversion symmetry within the electric dipole approximation [5]. Hence, when the nonlinear optical sources are strongly correlated, the coherent response from the system is markedly different from that of an assembly of non-correlated sources. Furthermore, when going beyond the electric dipole approximation, a direct view is opened on the molecular organization in the system [6][7][8]. To investigate this question more thoroughly, the spatial organization in small molecular aggregates has been probed by the frequency doubling conversion process known as second harmonic scattering (SHS). The latter terminology emphasizes that both an incoherent and a coherent response are observed as opposed to hyper-Rayleigh scattering (HRS), where only the incoherent response is observed. The first observation of this process has been reported by R.W. Terhune et al. and discussed by P.D. Maker, S.J. Cyvin et al. and R. Besohn et al. [9][10][11][12]. It has then been extensively used in the past to characterize the first hyper-polarizability of molecular compounds, supporting the extensive molecular engineering design of new organic compounds able to act as molecular probes [13][14][15]. The investigation of elaborated structures, like micelles, liposomes, or even metallic nanoparticles, has also been performed [16][17][18][19]. In these systems, the question of the phase retardation in the SHS response must be raised. In particular, in the case of metallic gold nanoparticles, where the diameter of the system can be tuned without changing other properties, this question has turned into the determination of the role of retardation in the electromagnetic fields to correctly describe the polarization-resolved response at the harmonic frequency [20,21].
In this work, we first present a theoretical foundation of the theory of second harmonic scattering of molecular aggregates. This is an extension of the standard theory developed for second harmonic scattering from molecules, also known as the theory of hyper-Rayleigh scattering, to encompass the case of molecular aggregates. In a second step, we illustrate this theory with two specific cases of molecular aggregates, the two cases differing by the internal molecular organization of the aggregate. In particular, we show unambiguously how the experimental data provide a quantitative insight into this aggregate molecular organization.
Hence, a general model for the SHS response of an assembly of a nonlinear dipole is first developed in order to investigate the impact of molecular organization of small molecular aggregates on the polarization-resolved SHS intensity. As a direct consequence, the weight of the dipolar and quadrupolar contributions is determined through a retardation parameter. This parameter describes the retardation effects, namely the spatial dependence of the electromagnetic field over the volume occupied by the aggregate, introduced in order to fully describe the SHS intensity response from aggregates. It is shown that this model encompasses the previous approaches restricted to the electric dipole approximation and used for the analysis of liquid solutions of well-dispersed, non-interacting molecular compounds, which now account for molecular correlations [22,23]. Experimental data are provided to illustrate the model. These data are obtained from aqueous solutions of Crystal Violet (CV), a non-centrosymmetric octupolar dye, in the presence of sodium dodecyl sulfate (SDS) and a mixed water-methanol solution of 4-(4-dihexadecylaminostyryl)-Nmethylpyridinium iodide (DiA), a cationic amphiphilic probe molecule with a strong first hyper-polarizability.

Experimental
Optics: The optical set-up has already been described elsewhere [24]. Briefly, to perform the SHS measurements, a femtosecond Ti-sapphire oscillator laser providing pulses with 150 fs duration with a repetition rate of 76 MHz was used. The fundamental wavelength was set at 800 nm. The average power was set at about 500 mW. The incident beam was focused into the sample cell. The latter consisted into a standard spectrophotometric cell with quartz windows. The focusing lens was standard consisting in a low numerical aperture microscope objective. The incident beam was linearly polarized and its polarization angle γ defined with a half wave plate. For vertically polarized light, γ = 0 (also noted v), whereas for horizontally polarized light, γ = π/2 (also noted h). The SHS intensity was collected at a right angle with a 5 cm focal length fused silica lens and sent to a monochromator coupled to a cooled photomultiplier tube. The detection was working in the gated photon counting regime. Owing to the low light level, the beam was chopped to remove noise with periods when the beam is blocked. The scattered harmonic light polarization state was selected with an analyzer. For the analyzer set vertically, the angle of polarization Γ = 0 (also noted V) was selected for vertical output polarization, whereas the angle Γ = π/2 (also noted H) was selected for horizontal output polarization. Color filters were used to remove any unwanted light along the beam path. The laser power was continuously monitored to account for intensity fluctuations.

Theory
SHS intensity for molecular aggregates: The electric field amplitude of light scattered by a molecule i located at position → r i at the harmonic frequency Ω = 2ω, where ω is the fundamental frequency in the direction , the expression of which is given by [25]: where K (Ω) = n(Ω)Ω/c is the harmonic wave vector modulus, with n(Ω) being the optical index of the medium at the harmonic frequency, andn = → r /r being the unit vector in the direction of collection. The expression of the nonlinear dipole induced at the harmonic frequency is related to the molecular first hyper-polarizability ↔ β m (i) through: where → E ( → r , ω) is the fundamental wave electric field, withê = cos γX + sin γŶ being its polarization vector defined with the polarization angle γ in the plane perpendicular to the propagation OZ direction. The frame transformation tensor ↔ T (r i ), accounting for the transformation of the molecular first hyper-polarizability tensor from the molecular frame to the laboratory frame, is defined with the standard Euler angles according to Figure 1 below.
For an assembly of correlated molecules, the total field amplitude → E ( → r , Ω) is the sum of the individual amplitudes given in Equation (1). The m molecules of the aggregate are assumed to be fixed relative to each other in space, hence, in position and orientation. Here, it is therefore simply assumed that the time scale of molecular motion in the aggregate is much longer than the time scale of the nonlinear interaction with the electromagnetic wave. With multiple scattering at both frequencies, and local field factors and molecular interactions neglected, then:

of 13
Symmetry 2021, 13, x FOR PEER REVIEW 4 of 14 Figure 1. Illustration of the three reference frames required, the molecular one (blue, smallcase labels), the aggregate one (red, A subscripted smallcase labels), and the laboratory one (black, uppercase labels).
For an assembly of correlated molecules, the total field amplitude is the sum of the individual amplitudes given in Equation (1). The m molecules of the aggregate are assumed to be fixed relative to each other in space, hence, in position and orientation. Here, it is therefore simply assumed that the time scale of molecular motion in the aggregate is much longer than the time scale of the nonlinear interaction with the electromagnetic wave. With multiple scattering at both frequencies, and local field factors and molecular interactions neglected, then: This superposition is built as a coherent superposition of the amplitudes scattered by single molecules. For a large number of molecules present in the aggregate, that is, if m is large, the discrete sum may be replaced by a continuous integration. In a liquid suspension, however, the molecules are distributed into many aggregates, and the latter are randomly distributed in the sampled volume in position and orientation. Hence, for N molecules distributed within n aggregates containing m molecules, one has m n N   , and the total intensity scattered at the harmonic frequency is the incoherent superposition of the intensities scattered by each molecular aggregate.
This intensity is therefore: where the brackets stand for the spatial averaging of the orientation of the aggregate, since it is dispersed in a liquid suspension. This description is the standard theoretical foundation for hyper-Rayleigh scattering.
In the developments of the above foundation for molecular aggregates, it is now important to introduce three distinguished reference frames: the laboratory frame (OXYZ), the molecular frame (Oxyz), and the aggregate frame (OxAyAzA) (see Figure 1). Indeed, the Figure 1. Illustration of the three reference frames required, the molecular one (blue, smallcase labels), the aggregate one (red, A subscripted smallcase labels), and the laboratory one (black, uppercase labels).
This superposition is built as a coherent superposition of the amplitudes scattered by single molecules. For a large number of molecules present in the aggregate, that is, if m is large, the discrete sum may be replaced by a continuous integration. In a liquid suspension, however, the molecules are distributed into many aggregates, and the latter are randomly distributed in the sampled volume in position and orientation. Hence, for N molecules distributed within n aggregates containing m molecules, one has N = n × m, and the total intensity scattered at the harmonic frequency is the incoherent superposition of the intensities scattered by each molecular aggregate.
This intensity is therefore: where the brackets stand for the spatial averaging of the orientation of the aggregate, since it is dispersed in a liquid suspension. This description is the standard theoretical foundation for hyper-Rayleigh scattering. In the developments of the above foundation for molecular aggregates, it is now important to introduce three distinguished reference frames: the laboratory frame (OXYZ), the molecular frame (Oxyz), and the aggregate frame (Ox A y A z A ) (see Figure 1). Indeed, the molecules possess a fixed position and orientation within the aggregate, but the aggregate is randomly distributed in the liquid suspension in position and orientation. Hence, in Equation (4), the averaging procedure is performed on the orientation of the aggregate. Owing to the intermediate frame, the frame transformation from the molecular to the laboratory frame can be performed into two successive steps for a single molecule (see Supplementary File Part A), namely: with the introduction of the frame transformation tensor from the aggregate to the laboratory frame ↔ T A and the molecular first hyper-polarizability tensor ↔ β A (i) of molecule i expressed in the aggregate reference frame. In order to get a better insight into the different contributions to the scattered intensity, Equation (3) is rewritten using the first order approximation: for simplicity and later discussion. Hence, keeping only the term in 1/r, Equation (3) yields: with the introduction of the molecular volume density ρ( → r ) within the aggregate and the wave vector mismatch ∆ → k , the expression of which is simply: with the introduction of the fundamental wavelength λ. If the direction of collection is in the forward direction, namelyn =Ẑ, then the usual phase matching condition of second harmonic generation (SHG) is recovered. Because the condition r << λ is also assumed, i.e., the aggregates are small before the fundamental wavelength, the expansion of the exponential factor of the integral in Equation (7) yields: It is therefore convenient to split the harmonic field amplitude into two terms, namely: The first term → p e f f ,L (Ω) in Equation (10) is the equivalent of Equation (3) at the level of the aggregate instead of the molecule. Hence, the integrated induced dipole of the aggregate → p e f f ,L (Ω) is simply: where the dipolar first hyper-polarizability of the aggregate is given by: i.e., the superposition of the first hyper-polarizabilities of the molecules present in the aggregate. This superposition takes into account both the position and the orientation of the molecules. The second term ↔ q e f f ,L (Ω) of Equation (10) is the first correction to the non-vanishing spatial extension of the aggregate. Therefore, ↔ q e f f ,L (Ω) is the integrated induced quadrupole of the aggregate, and is given by: with the quadrupolar first hyper-polarizability ↔ γ e f f ,A of the aggregate defined through: Symmetry 2021, 13, 206 6 of 13 Note that the vector → r has to be also transformed into the aggregate reference frame.
In the molecular frame, it is simply given by ↔ γ e f f ,A . Finally, it is observed that, from Equation (10), the second contribution scales with the parameter r /λ. Since it has been set beforehand that this parameter is much smaller than unity, this contribution is usually neglected for molecular aggregates, since their average size is much smaller than the fundamental wavelength. Typical aggregates have diameters of the order of a few nanometers, whereas the fundamental wavelength is often in the range of hundreds of nanometers. However, depending on the molecular organization of the aggregate, the induced dipole → p e f f ,L (Ω) may be very small or even vanish. In this case, the quadrupolar contribution may not be negligible any longer, and must be accounted for.
SHS scattered intensity for non-centrosymmetrical aggregates: In the particular case of a non-centrosymmetrical spatial orientational distribution of the molecules in the aggregates, and I = X for Γ = V and I = Z for Γ = H. In particular, the depolarization ratio D = c V /a V is now an indication of the symmetry of the aggregate, and no longer of the molecule. Furthermore, the total intensity scales linearly with the number of aggregates present in the liquid suspension. Note that the resulting apparent first hyper-polarizability of the aggregates depends now on the orientation of all molecules constituting the aggregate, and can be very large, since it is the superposition of the first hyper-polarizabilities of the m molecules contained in the aggregate. For instance, in a perfect alignment of all molecules, the aggregate first hyperpolarizability is indeed as large as m times that of a single molecule, and hence, the intensity scattered at the harmonic intensity scales with the product nm 2 . A typical graph of the vertically polarized SHS intensity collected in the case of a noncentrosymmetric aggregate with an arbitrary depolarization ratio of D = 0.4 is provided in Figure 2 below. SHS scattered intensity for centrosymmetrical aggregates: The second case occurs when the organization of the molecules within the aggregate possesses a center of inversion. This is the organization found in spherical micelles, for instance, or corresponding to a random orientation. As expected for this centrosymmetrical structure, the first contribution vanishes, namely The scattered amplitude at the harmonic frequency is therefore reduced to the second contribution only, which may be recast as: Finally, the intensity again takes the expressions given by Equations (15a) and (15b), but this time:   SHS scattered intensity for centrosymmetrical aggregates: The second case occurs when the organization of the molecules within the aggregate possesses a center of inversion. This is the organization found in spherical micelles, for instance, or corresponding to a random orientation. As expected for this centrosymmetrical structure, the first contribution vanishes, namely → p e f f ,L (Ω) = 0. The scattered amplitude at the harmonic frequency is therefore reduced to the second contribution only, which may be recast as: Finally, the intensity again takes the expressions given by Equations (15a)  with the constant: (20) and the effective parameters: Γ Γ e f f ,L,I JK = n(Ω)γ e f f ,L,YI JK + n(ω)γ e f f ,L,ZI JK (21) with I = X for Γ = V and I = Z for Γ = H. The case of spherical micelles of radius a leads to simple expressions for the intensity coefficients of Equations (19a-c). With a uniform distribution of the molecular dye at the aggregate surface, the molecular density is simply: if one assumes that the nonlinear optical chromophore is point-like. The volume integration for the effective quadrupolar first hyper-polarizability elements leads to a dependence of these elements with the third power of the radius of the sphere. Furthermore, the coefficients a V and c V for the vertical polarization vanish, whereas b V does not, and is proportional to b V ∝ ∆k 2 m 2 a 2 (see Figure 3). For the horizontal polarization, one gets b H = 2a H = 2c H . Furthermore, for an experimental geometry where the collection is performed along the opposite direction of the OY axis of the laboratory, Equation (21)  The case of spherical micelles of radius a leads to simple expressions for the intensity coefficients of Equations (19a-c). With a uniform distribution of the molecular dye at the aggregate surface, the molecular density is simply: (22) if one assumes that the nonlinear optical chromophore is point-like. The volume integration for the effective quadrupolar first hyper-polarizability elements leads to a dependence of these elements with the third power of the radius of the sphere. Furthermore, the  The perfect centrosymmetry rule is not expected to be strictly fulfilled by aggregates However, if this symmetry is close to perfect, contributions from the first contribution ) ( , Ω L eff p  will also contribute simultaneously to the total intensity scattered at the harmonic frequency without overwhelming the second one. Therefore, the resulting polarization-resolved plot may take a profile somewhere in between the two contributions  The perfect centrosymmetry rule is not expected to be strictly fulfilled by aggregates. However, if this symmetry is close to perfect, contributions from the first contribution → p e f f ,L (Ω) will also contribute simultaneously to the total intensity scattered at the harmonic frequency without overwhelming the second one. Therefore, the resulting polarizationresolved plot may take a profile somewhere in between the two contributions presented in Figures 2 and 3. Note also that other orientational distributions may be found, like the azimuthal one, instead of the radial one.
SHS of Molecular Aggregates: To illustrate the problem, we now discuss the case of aggregates constituted by Crystal Violet and DiA, two compound-forming aggregates with different properties depending on the conditions, see below.
Experiments are initially performed for a simple rod-like molecule, 4-(4-dihexadecylaminostyryl)-N-methylpyridinium (DiA) (see Figure 4), the first hyper-polarizability tensor of which is dominated by a single element. This latter dominating element is β m,zzz (i) where the molecular Oz axis is taken along the axis of the charge transfer. Since DiA is rod-like, this molecular axis is nearly oriented along the axis direction of DiA. Since DiA is also amphiphilic due to the two long alkyl chains, it is expected to form micelle-like aggregates, with the Oz axis oriented towards the center of the aggregate. The orientation distribution within the molecular aggregate is therefore expected to closely resemble that of spherical micelles. DiA dissolves in methanol but not in water, and therefore aggregation was induced by varying the volume fraction of water into methanol at a fixed DiA concentration of 12.5 µM [25]. A typical vertically polarized resolved SHS intensity for DiA dispersed in pure methanol is given in Figure 5. This graph allows to immediately identify a V as I V (γ = 0), c V as I V (γ = π/2), and b V as 4I V (γ = π/4) − a V − c V , although an adjustment procedure with Equation (15a) is usually performed, and the condition b V = a V + c V is obeyed.
Symmetry 2021, 13, x FOR PEER REVIEW presented in Figures 2 and 3. Note also that other orientational distributions found, like the azimuthal one, instead of the radial one.
SHS of Molecular Aggregates: To illustrate the problem, we now discuss th aggregates constituted by Crystal Violet and DiA, two compound-forming ag with different properties depending on the conditions, see below.
Experiments are initially performed for a simple rod-like molecul dihexadecylaminostyryl)-N-methylpyridinium (DiA) (see Figure 4), the f per-polarizability tensor of which is dominated by a single element. This latte nating element is β m,zzz (i) where the molecular Oz axis is taken along the ax charge transfer. Since DiA is rod-like, this molecular axis is nearly oriented along direction of DiA. Since DiA is also amphiphilic due to the two long alkyl chains pected to form micelle-like aggregates, with the Oz axis oriented towards the c the aggregate. The orientation distribution within the molecular aggregate is t expected to closely resemble that of spherical micelles. DiA dissolves in methano in water, and therefore aggregation was induced by varying the volume fraction into methanol at a fixed DiA concentration of 12.5 μM [25]. A typical vertically p resolved SHS intensity for DiA dispersed in pure methanol is given in Figure   graph allows to immediately identify a V as  In order to have a quantified characterization of the aggregate, the following ζ V parameter is defined: also found as [6]: Note in particular the sign inversion between I 4 and ζ V . As examples, it is immediate to see that the pure dipolar response of a non-centrosymmetric aggregate leads to ζ V = 0, whereas the purely retarded response of a centrosymmetric aggregate in the shape of a spherical micelle, i.e., with a radial distribution, leads to ζ V = 1. Interestingly, an azimuthal distribution of the dyes forming a centrosymmetric aggregate leads to a negative value for the retarded parameter ζ V .
As an illustration of a quantitative analysis of the aggregate structure, one may discuss the results reported in Ref. [25], comparing the vertically polarized SHS intensity obtained for DiA dispersed in pure methanol and in a mixed 5:1 v/v water-methanol solution. The depolarization ratio D = c V / a V is found to be 0.25, along with a vanishing retardation parameter ζ V for DiA dispersed in pure methanol (see Figure 5). Oppositely, in a mixed 5:1 v/v water-methanol solvent, the depolarization ratio D = c V / a V is found to be 0.28, along with a retardation parameter ζ V = 0.41 [25]. In order to have a quantified characterization of the aggregate, the following ζ V parameter is defined: also found as [6]: Note in particular the sign inversion between I 4 and ζ V . As examples, it is immediate to see that the pure dipolar response of a non-centrosymmetric aggregate leads to ζ V = 0, whereas the purely retarded response of a centrosymmetric aggregate in the shape of a spherical micelle, i.e., with a radial distribution, leads to ζ V = 1. Interestingly, an azimuthal distribution of the dyes forming a centrosymmetric aggregate leads to a negative value for the retarded parameter ζ V .
As an illustration of a quantitative analysis of the aggregate structure, one may discuss the results reported in Ref. [25], comparing the vertically polarized SHS intensity obtained for DiA dispersed in pure methanol and in a mixed 5:1 v/v water-methanol solution. The depolarization ratio D = c V /a V is found to be 0.25, along with a vanishing retardation parameter ζ V for DiA dispersed in pure methanol (see Figure 5). Oppositely, in a mixed 5:1 v/v water-methanol solvent, the depolarization ratio D = c V /a V is found to be 0.28, along with a retardation parameter ζ V = 0.41 [25].
The depolarization ratio obtained for the aggregate is different from the value of +0.25 measured for DiA dissolved in pure methanol and closely associated with a single first hyper-polarizability tensor element compound [25]. It is concluded that some distortion in the electronic structure of DiA occurs in the aggregate, leading to the appearance of contributions from other tensor elements. Furthermore, the retardation parameter, vanishing in the pure methanol solution where DiA is ideally dispersed, now increases, indicating an organization tending to the spherical micellar-like structure, as expected for a compound with two hydrophobic long alkyl chains.
In the case of aqueous solutions of Crystal Violet, the same experiments were performed, this time as a function of SDS concentration, as this anionic surfactant is known to form micelles beyond the critical micelle concentration (CMC) of 7 mM [24]. It is expected that these molecular structures favor adsorption [26]. Figure 6 shows the depolarization ratio D along with the retardation parameter ζ V as a function of the SDS concentration, below and above the SDS CMC.
The depolarization ratio obtained for the aggregate is different from the value of +0.25 measured for DiA dissolved in pure methanol and closely associated with a single first hyper-polarizability tensor element compound [25]. It is concluded that some distortion in the electronic structure of DiA occurs in the aggregate, leading to the appearance of contributions from other tensor elements. Furthermore, the retardation parameter, vanishing in the pure methanol solution where DiA is ideally dispersed, now increases, indicating an organization tending to the spherical micellar-like structure, as expected for a compound with two hydrophobic long alkyl chains.
In the case of aqueous solutions of Crystal Violet, the same experiments were performed, this time as a function of SDS concentration, as this anionic surfactant is known to form micelles beyond the critical micelle concentration (CMC) of 7 mM [24]. It is expected that these molecular structures favor adsorption [26]. Figure 6 shows the depolarization ratio D along with the retardation parameter ζ V as a function of the SDS concentration, below and above the SDS CMC. The depolarization ratio is 0.6 in the absence of SDS, as expected for the close to planar octupolar symmetry of CV [27]. It increases when SDS is present, but does not change much, indicating a change in symmetry of the whole aggregate compared to the CV monomer. Oppositely, the retardation parameter ζ V is initially zero, as expected for a compound ideally dissolved in its solvent, but then decreases as SDS is introduced into the solution. The negative value of this parameter is not in agreement with a spherical micellar-like organization, but rather with an azimuthal distribution, i.e., adsorption onto the SDS aggregates. Interestingly, this organization already starts at 3 mM, a value below the SDS CMC, suggesting the appearance of aggregates below the CMC.

Conclusions
We have developed a model to describe assemblies of nonlinear dipoles exhibiting two dominating contributions, one of electric dipole origin associated to the non-centrosymmetric organization of the compounds into the aggregates, and one asso- The depolarization ratio is 0.6 in the absence of SDS, as expected for the close to planar octupolar symmetry of CV [27]. It increases when SDS is present, but does not change much, indicating a change in symmetry of the whole aggregate compared to the CV monomer. Oppositely, the retardation parameter ζ V is initially zero, as expected for a compound ideally dissolved in its solvent, but then decreases as SDS is introduced into the solution. The negative value of this parameter is not in agreement with a spherical micellar-like organization, but rather with an azimuthal distribution, i.e., adsorption onto the SDS aggregates. Interestingly, this organization already starts at 3 mM, a value below the SDS CMC, suggesting the appearance of aggregates below the CMC.

Conclusions
We have developed a model to describe assemblies of nonlinear dipoles exhibiting two dominating contributions, one of electric dipole origin associated to the noncentrosymmetric organization of the compounds into the aggregates, and one associated with phase retardation, providing a deep insight at the nanoscale into the aggregate spatial organization. Illustration of this model with the two CV and DiA compounds forming aggregates in different conditions shows that, in particular, the sign of the retardation parameter may be crucial in disentangling the aggregate spatial organization.
On a final note, we point out that this model may be applied to other systems, like metallic nanoparticles [28].