Spectral dependence of femtosecond laser induced circular optical properties in silica

Transmission Mueller-matrix spectroscopic ellipsometry is applied to study femtosecond laser induced nanogratings in silica glass in a wide spectral range 250 – 1800 nm. By using differential decomposition of the Mueller matrix, the circular birefringence and dichroism of femtosecond laser irradiated SiO2 are quantified for the first time in the UV and Near-IR range. A maximum value of the effective specific rotation of α ~ 860°/mm at 290 nm is found. In the near-IR range, we found a linear and circular dichroism band peaking around 1240 nm, which might be attributed to the formation of anisotropic species like the formation of oriented OH species and Si-O-Si bond. © 2018 Optical Society of America under the terms of the OSA Open Access Publishing Agreement


Introduction
Femtosecond laser direct writing (FLDW) starts with multiphoton ionization [1,2], which induces various permanent modifications in silica glass depending on the laser parameters [3], Such changes include refractive index changes [4], the formation of porous nanogratings and related linear birefringence [5][6][7], or even voids.The processing of glass with low-energy femtosecond laser pulses (sometimes combined with chemical etching) defines a flexible manufacturing platform suitable for technologies like optical data storage, optofluidics, optomechanics, marking, and photonics such as 3D optical waveguides and their combination [8,9], fiber Bragg gratings [10,11] and polarization devices [12,13].Manufacturing monolithically integrated devices is attracting great interest because it is possible to induce strong refractive index increases (Δn from 10 -3 to more than 10 -2 in a large range of transparent glasses), which are localized in the bulk due to strong non-linear effect.Controlling these modifications could potentially allow exceeding current fs laser applications opening additional new possibilities in material science.In particular, the industrial demand for sophisticated all-optical integration is increasingly large as the volume of exchanged data and information increases demanding more complex devices and therefore greater nuance and complexity of laser treatment.Presently, one crucial optical property is still missing from the component library: the integration of circular optical properties (optical rotation and circular dichroism) in glass optical materials and devices.
In a previous study done in 2008, based on measurements done with a phase shift interferometer, we revealed that femtosecond laser interactions can shear matter like a scissor giving rise to a chiral strain if the light propagation axis is considered [14,15].But shearing does not give rise to optical circular properties by itself like it is with torsion.Taylor et al.
reported highly ordered « chiral-like » nanostructures using circular polarized laser light but no circular optical properties [16].In this direction, we have shown in 2016 that it is possible to produce circular diattenuation by inducing and controlling chirality within an achiral glass namely SiO 2 by controlling the laser writing polarization with respect to the scanning direction [17].Last estimates of equivalent circular diattenuation using a circular dichroism spectrophotometer [18,19] show values up to -30°/mm in the UV wavelength range and -20°/mm at λ = 550 nm.However it was difficult to extract the real spectral shape of the circular diattenuation due to the co-existence of strong linear properties [19].
Recently Mueller-matrix spectroscopic ellipsometry was used in the visible range to investigated silica glass modified by femtosecond laser within nanogratings regime [20,21].This method is a powerful tool to examine structural and optical properties of surfaces, thin films, and multilayered materials that exhibit both linear and circular optical properties from a single measurement.We applied the differential matrix formalism [22,23] on transmission Mueller matrices measured for femtosecond laser-induced modifications in silica within the type II laser-processing regime.We showed how values of the effective structural circular birefringence can be extracted in the presence of strong linear optical properties and found a maximum value of the effective specific rotation α D ~ 143°/mm at 550 nm.
In this letter, the spectral dependence of both linear and circular optical properties are investigated according to pulse energy and in a wide spectral range extending from λ = 250 nm up to λ = 1800 nm.In particular, we reveal a large effective optical rotation in the UV range and a specific dichroism band in the near-IR range when the laser polarization is oriented at 45° off the scanning direction.

Experimental details
The laser radiation was produced by a femtosecond laser system operating at λ = 1030 nm and delivering 250 fs pulses at 100 kHz repetition rate.The beam was focused to a depth 300 µm below the front face of 1 mm thick silica glass plates (Heraeus Suprasil1 containing [OH] ~ 830 ppm and a solgel silica with [OH] ~ 100 ppm) using a 0.6 NA aspheric lens (estimated beam waist w ~ 1,5 µm).The laser energy was varied from 0,15µJ up to 4µJ so we crossed the type II threshold (typ.0.3µJ in our experimental conditions) corresponding to the formation of nanogratings in SiO 2 .Based on preliminary results, we choose to investigate here a specific writing configuration in terms of polarization orientation that maximizes the imprinting of circular optical properties.When the laser was moving along x-axis (horizontal direction) and the laser linear polarization was oriented θ = +45° from the writing direction, we define it as "X+45°" configuration of writing.Then, in this X+45° configuration and using a scanning speed of v = 1 mm/s, we created several squares 3 x 3 mm 2 , made up of a set of lines with a line spacing Δy = 1 µm to have a uniform anisotropic area and avoiding any diffraction effects.Scheme and more details about the laser writing procedure can be found in Ref. [17].
Anisotropic optical properties are investigated using a phase modulated spectroscopic ellipsometer (UVISEL+, HORIBA Scientific) over λ = 250 to 1800 nm spectral range.Since this equipment gives access to the three first columns of the Mueller matrix [24], we will explain in the results section how to estimate the values of the 4 th column based on symmetry properties of Mueller matrices.In order to validate this approach, we compared the data obtained with the UVISEL+ with data obtained with a home-made Mueller ellipsometer [24,25] (Smart-SE, HORIBA Scientific), which provides a full Mueller matrix over a shorter spectral range between λ = 450 and λ = 1000 nm.All measurements were made using a collimated probe beam in normal incidence.The probe beam size of these two instruments was fixed to 1mm for all measurements.The samples were oriented in such a way their writing/scanning axis was set horizontal +/-2° for the referential of the Mueller ellipsometers.

Results
We have seen in our preliminary experiments that our samples show asymmetric transmission when measuring ellipticity using a circular dichroism spectrometer [19].In 2018 we investigated this observation by measuring both M and M rev , i.e., the Mueller matrix in the reverse direction (from the backface) [21].The differences were that m 13 , m 31 , m 24 , and m 42 change their respective signs whereas all other elements remain invariant.This symmetry implies that M fulfills the principle of reciprocity.The difference between the elements of predicted M rev from a measured Mueller matrix, M, and the elements of the measured M rev is typically < 0.01 which corresponds to instruments accuracy.This indicates that both linear and circular anisotropic optical properties are quite similarly distributed in depth.
Fig. 1 shows the normalized (m 11 was normalized to 1) spectral Mueller matrices M obtained with the two instruments.The Mueller matrices represented in Fig. 1 correspond to a Suprasil1 sample written at pulse E = 1.5 µJ/pulse with a X+45° writing configuration.The sample exhibits low linear and circular dichroism because the elements m 12 , m 21 , m 14 and m 41 are very small.In contrast, the elements m 24 and m 42 , being close to +1 and -1 respectively highlight a significantly high linear birefringence [26].The Mueller matrix of this X+45° sample, can be firstly approximated by M = R(−45°).M ΨΔ .R(45°), where R(α) is the rotation matrix and M ΨΔ is the Mueller matrix for a perfect linear dichroic retarder.This corresponds to the slow axis direction coinciding with the chosen writing polarization orientation; that is, rotated α = 45° from the direction of writing.The element m 33 should ideally be equal to one.However, the measured values appear to be slightly smaller (0.95-0.99) in the blue indicating a slight depolarization in this range, and down to 0.85 in the UV range at λ = 250 nm.In summary, our samples for λ ∈ 450-1800nm are characterized by having high linear retardation, low dichroism, and a small depolarization.Note there is an interesting dichroic band around λ = 1300 nm.In the UV range, there is an increase of the linear dichroism and the depolarization value is around (0.85-0.95).
So as a first approximation if we consider that our sample is indeed a linear dichroic retarder at 45° (i.e.circular properties are very small compared to the linear ones), we can use symmetry properties of Mueller matrices to evaluate the values of the 4 th column i.e. m 14 = m 41 , m 24 = -m 42 , m 34 = m 43 , m 44 = m 22 .Note that both matrices coincide quite well in the common spectral range so this first approximation to deduce the values of the 4 th column of the Mueller matrix from phase-modulated ellipsometry could provide some reliable information.Note that the differences between the two ellipsometers (e.g. up to 0.1 in m 32 and m 34 ) are related to a slightly different azimuthal orientation (+/-2°) of the samples during measurements.This leads to different "projections" of the linear properties in the xy plan but this has no impact on the intrinsic polarimetric properties that we extract after decomposition.
As a second order approximation we can take into account that elements m 14 and m 41 is formed by a superposition of effects due to both circular and linear parts, with the linear contribution intervening with opposite signs in m 14 and m 41 respectively as it has been extensively study in refs [27][28][29].Therefore since linear properties are non-zero, one expects that before decomposition m 14 be different than m 41 .So the dash blue line in Fig. 1 corresponds to m 14 taking into account the linear, which agree quite well with the full Mueller ellipsometer (Smart-SE).This enables further data treatment in order to determine the complete polarimetric properties of our fs-irradiated samples.According to [23] the polarimetric optical response of a medium of length, l, with a complex refractive index , can be defined as a superposition of the following optical phase shift related to basic polarimetric properties: linear birefringence , 45° linear birefringence , linear dichroism , 45° linear dichroism , circular birefringence and circular dichroism . By using the differential decomposition described in [5][6], it is possible to extract all the above mentioned polarimetric properties from the Mueller matrix of a sample providing that the laser track is considered as homogeneous in the direction of light propagation.A differential decomposition of M to find its differential matrix, m, implies determining its logarithm L = ln(M).Since dM⁄dz = m(z).M(z), if the differential matrix m, does not depend on z, then the matrices m and L are related by L = m.l.For a non-depolarizing medium, the differential matrix L is Minkowski antisymmetric and contains the six elementary optical properties fully characterizing the medium namely LD, LB, LD', LB', CD and CB.If the medium is depolarizing as it is the case here in the UV range (up to 15% depolarization at λ = 250 nm) [19], the depolarization can be interpreted as resulting from the fluctuations of the optical properties and , which influence the elementary polarization properties of the irradiated silica glass [22,23].Result of the differential decomposition of the experimental Mueller matrix is shown in Fig. 2. Based on preliminary investigations [17,19] we know that fs irradiation modifies inhomogeneously the physical structure of the material along the pulse path across the sample.Consequently, we think that the measured birefringence parameters do not correspond to intrinsic material properties but to effective or equivalent optical parameters induced by fs laser.For sake of comparison we add the CB and CD spectra obtained from the decomposition of the Mueller matrix using the two different approximations describe above: first one (full line) considers circular properties are very small when compared to the linear ones whereas the second one (dash line) introduce a linear term correction to calculate the elements of 4 th column as described in [27][28][29].
We can observe three spectral regions in Fig. 2: (I) below λ = 450 nm, (II) Δλ = (450 -1100) nm and above λ = 1400nm, and (III) between λ = 1100 and λ = 1400nm.In region I, the absolute values of the off-diagonal elements in M are quite high which have values above 0.3 rad.For example, at λ = 250 nm, we reveal 0.5 rad for LD', -3 rad for LB' and significant circular properties as well 0.3 rad for CD and -0.75 rad for CB.Note that the reliable presence of CD in the UV range has previously been confirmed using a spectro-polarimeter in [18].In region II, all depolarization effects remain smaller than 5% and monotonically decrease with λ below 1% for λ > 800 nm.According to the data shown in Fig. 2, it is seen that the two xy linear effects (LB and LD) are small, while the crossed linear properties (LB' and LD') are quite important.This is because the polarization of the sample (X+45°) written by laser was oriented ±45° to the reference axis.If the polarization had been oriented parallel or perpendicular to the x-direction (X0° or X90°), the pair (LB and LD) would have been maximal and the pair (LD' and LB') would have been minimized.For the circular effects we notice a monotonous increase of CD towards shorter wavelengths.
In region III Δλ = (1100 -1400) nm, the sample exhibits a strong linear dichroism at 45° (LD') band.In particular we notice that m 13 ≈ m 31 with LD' values up to 0.25.For the chosen orientation we also observe that m 14 ≈ m 41 and is nonzero with CD values up to 0.10, which reveals a slight circular dichroic band.The corresponding Kramers-Kronig consistent feature in LB' and CB also creates a related feature ~ λ = 1240 nm.The spectral dependence of both linear and circular anisotropic optical properties are shown respectively in Figs. 3 and 4 as a function of the pulse energy.In fig. 3 we chose to draw the "Total LB" and the "Total LD" that correspond to the following equations and .Note that these two properties are thus independent of the samples azimuthal orientation during the ellipsometry measurements.
Below 1100 nm, the amplitudes of the anisotropic linear properties are increasing monotonously with increasing pulse energy, E, from 0.15 µJ up to 4 µJ.Above 0.5 µJ, we can observe a linear dichroism band within Δλ = (1100 -1400) nm as revealed in Fig 3 (Left).The amplitude of this LD band increases with the energy of the fs pulse and the sample shows a maximum LD of 0.27 rad at Δλ = 1240 nm for E = 2 µJ.The spectral position of the band also shifts from 1380 nm at E = 0.5 µJ down to 1240 nm at 2 µJ.This agrees with the corresponding Kramers-Kronig-consistent feature in the LB spectra that creates an increasing spectral feature around Δλ ~ 1240 nm up to 2 µJ as shown in Fig. 3(Right).However above 3 µJ, we observe a decrease of the LD amplitude and the band shifts back to longer wavelengths.Now comparing the two OH samples written with E = 1,5 µJ, we can see that the « background LB » (away from the dichroic band e.g.around 1000nm) did not almost change irrespectively of the OH species concentration in agreement with [30].In contrast the amplitude of LD band(s) around 1100-1400nm is much lower in the "low OH" sample (100ppm, full brown line).It is worth noting that the dichroic band appears to be composed of 2 or even 3 unknown contributions that still remain to be identified.
The corresponding spectra related to the circular optical properties are shown in Fig 4. Similar to the linear properties, we observed a CD band in the 1100 -1400 nm spectral range, which evolves in the same way (increasing amplitude and blue shift) as the LD band according to the pulse energy up to 2 µJ.It is worth noting that the CD amplitude in the near IR is much lower (by a factor 5) in the "low OH" sample.In addition the CB spectra are consistent with Kramers-Kronig and a related feature that can be seen in Fig. 4(Right) around λ ~ 1250 nm, especially up to E = 2-3 µJ.Then from 3 µJ, the dependence is more complex since the LD amplitude is decreasing whereas CD band in Fig. 4(Left) is strongly increasing.This indicates that CD and LD amplitudes are not simply connected to each other.Following the same trend, CD and CB are no more Kramers-Kronig consistent for E = 3 µJ.This suggests that the origin of the circular properties is likely not related to chiral molecular species but rather to a chiral arrangement of anisotropic linear layers along the probe beam pathway.In the later case, a non-parallel non-orthogonal assembly of two LB (or LD) contributions would create some effective CB (respectively CD) whose amplitude is proportional to the LB (or LD) and the misalignment angle between the two contributions [31].

Discussion
Apart from the spectral feature around λ ~ 1250 nm, the linear optical properties are mainly related to the formation of nanogratings resulting in a strong linear form birefringence.However, there is also a slight contribution of stress-induced birefringence, which has been investigated in the literature [32,33].We reported earlier that the spectral dependence of the optical pathway difference expressed as is quite flat in the investigated spectral range, which results in an increase of at short wavelengths.An increasing apparent LD at short wavelength accompanies the linear birefringence, which is related to the nanogratings.The positive LD' (and nearly zero LD) observed in Fig. 2 implies that a higher losses were measured for polarization oriented perpendicular to the nanolayers in agreement with [34].It is known that a layered medium made of alternating layers of two different isotropic materials with complex refractive indices, exhibits a linear dichroism [35].However the linear dichroism observed in the UV-Vis range should be rather called linear diattenuation since it can be mostly attributed to polarization dependent scattering [34] due to the intrinsic nanoporous nature of the nanogratings.About the observed circular optical properties in the UV-Vis range, we do not know yet if the observed circular birefringence originates from internal linear birefringence that could be related to a non-parallel non-orthogonal assembly of two (or more) linear contributions, or from a fs pulse induced molecular optical activity.However, we may, from the observed structural circular birefringence, estimate the value of an effective specific rotation; i.e. circular birefringence per unit length to facilitate comparison with naturally optically active materials.From Fig. 2 we find that CB at λ = 290 nm has an extreme value of -0.75 rad produced by the birefringent laser track that exhibit a length of l = 50 µm (in the z-propagation direction).
Therefore the specific rotation is -860°⁄mm.
Following the same view in the Near-IR range, we can deduce from Fig. 4 a maximum CB of -0.22 rad at 1240 nm for 4 µJ/pulse corresponding to an equivalent specific rotation estimated to -126°⁄mm.This could be compared with the specific rotation of natural materials, like quartz (less than 5°/mm at 1240 nm and 49°/mm at 400nm) or the structural rotary power of liquid crystals (cholesteric monocrystalline layer can reach thousands of °/mm) [36].Within this specific wavelength range, we can observe an intense dichroism band both linear and circular.
These bands might be related to vibrational absorption bands ascribed to the combination of 2.ν stretching Si-OH + ν stretching Si-O-Si.Another possibility would be the combination 2.ν stretching Si-OH + ν bending Si-OH [37].Note there are some variants Si-OH silanol groups as described in [38], which could explain the occurrence of 2-3 bands as observed in the "low OH" sample.However there is no CD band at 1400 nm corresponding to the first harmonic of SiOH stretching, 2.νOH, which indicates that CD origin might be related to Si-O-Si bonds arrangement.This was tentatively interpreted through the breaking of symmetry arising from a volume torque due to the combined action of the stress field and a DC electric field (defined by the pulse front tilt, the focusing conditions and the laser polarization).This volumetric torque may contribute to twist local matter and create chiral atomic arrangements [17].However, it is unlikely that this could create a high enough population of Si-O-Si bonds with a chiral arrangement (such as tetrahedra in α-quartz) to achieve the observed high CD.Another possibility would be that CD band originates from several internal linear dichroism and birefringence contributions with non-parallel neutral axes.Following this view, the amplitude of the observed CD band is not only related to the LD amplitude of each contributions but also to the misalignment between their neutral axes.LB neutral axis being mainly related to the nanogratings contribution.

Conclusion
We have investigated the polarimetric properties of femtosecond laser induced nanogratings in a wide spectral range extending from UV up to the NIR.The equivalent circular optical properties using Mueller-matrix spectroscopic ellipsometry in transmission mode were extracted.A high effective specific rotation of α ~ 860°/mm at λ = 290 nm is found.We also reveal the formation of a linear and circular dichroism band peaking around λ ~ 1240 nm.

Fig. 1 .
Fig. 1.Normalized Mueller matrix of Suprasil1 sample written according to the X+45° configuration.Red lines are for the liquid crystals ellipsometer whereas the one obtained with photoelastic modulator ellipsometer correspond to the black lines.The dash blue line corresponds to the second order approximation taking into account the linear contribution.