Random lasing and replica symmetry breaking in GeO2

Oct 19, 2023

Scientific Reports volume 12, Article number: 19438 (2022) Cite this article

605 Accesses

2 Citations

1 Altmetric

Metrics details

We investigated the random lasing process and Replica Symmetry Breaking (RSB) phenomenon in neodymium ions (Nd3+) doped lead-germanate glass–ceramics (GCs) containing MgO. Glass samples were fabricated by conventional melt-quenching technique and the GCs were obtained by carefully devitrifying the parent glasses at 830 °C for different time intervals. The partial crystallization of the parent glasses was verified by X-ray diffraction. Photoluminescence (PL) enhancement of \(\approx\) 500% relative to the parent glasses was observed for samples with a higher crystallinity degree (annealed during 5 h). Powders with grains having average size of 2 µm were prepared by griding the GCs samples. The Random Laser (RL) was excited at 808 nm, in resonance with the Nd3+ transition 4I9/2 → {4F5/2, 2H9/2}, and emitted at 1068 nm (transition 4F3/2 → 4I11/2). The RL performance was clearly enhanced for the sample with the highest crystallinity degree whose energy fluence excitation threshold (EFEth) was 0.25 mJ/mm2. The enhanced performance is attributed to the residence-time growth of photons inside the sample and the higher quantum efficiency of Nd3+ incorporated within the microcrystals, where radiative losses are reduced. Moreover, the phenomenon of Replica Symmetry Breaking (RSB), characteristic of a photonic-phase-transition, was detected by measuring the intensity fluctuations of the RL emission. The Parisi overlap parameter was determined for all samples, for excitation below and above the EFEth. This is the first time, for the best of the authors knowledge, that RL emission and RSB are reported for a glass–ceramic system.

Laser action in disordered media, without optical cavities, has been the object of intense theoretical and experimental studies since the pioneering work of Letokhov1,2,3,4,5,6,7,8,9,10,11. In this kind of laser systems, currently referred as Random Lasers (RLs), the feedback mechanism contributing for optical amplification is not achieved by a well-engineered optical cavity as in conventional lasers. Instead, optical feedback is accomplished by light scattering due to refractive index inhomogeneities within a disordered medium12,13.

RLs were reported for several systems up to now. For example, the disordered scattering particles may be embedded within the gain medium, such as in laser-dyes liquid solutions containing high refractive index particles in suspensions14,15. The dyes may also be incorporated into solid matrices, such as polymer membranes16, biological tissues17, and glasses produced by sol–gel18, among others.

Random lasing was also extensively reported for rare-earth ions (REI) doped crystalline powders19,20,21,22,23,24,25. In these systems the particles act both as gain medium and as scatterers. In particular, random lasing can also be obtained from REI doped optical fibers, where the feedback may be obtained due to light reflections in random Bragg gratings written in fibers with nonuniform refractive indices26,27, or fibers with phase separated glass cores28.

Curiously, RL reports based on glassy particles doped with REI are very scarce. Several years ago, an upconversion RL emitting in the UV, was reported based on a fluoroindate glass-powder doped with neodymium ions29. More recently, we demonstrated RL action in neodymium (Nd3+) doped zinc-tellurite glass powders30. The RL feedback mechanism was provided by the light reflections on the glassy grains-air interfaces.

Works on RL in glass–ceramics (GCs) are rare as well31,32. Nevertheless, GCs are interesting media for photonic devices, since they can withstand high power excitations and have a high thermal threshold. Furthermore, GCs can be heavily doped with rare-earth ions to alter its emission characteristics31. Despite that, to the best of our knowledge, the present article is the first report of a RL based on Nd3+ doped glass–ceramics. Our aim was to evaluate and characterize the influence of the crystallization degree of the GCs powder on the performance of RLs. For this research, the choice of the parent lead-germanate glass was made for several reasons discussed below.

Lead-germanate glasses are strong candidates for RLs operation because they present high solubility for REI, high refractive index (~ 2.0), high transmittance window (400–5000 nm)33, and a high resistance to optical damage34. Addition of MgO in the glass composition contributed as an intermediate oxide in the glass structure, acting not only as a modifier, but also as a network former35,36.

Besides the basic characterization of the RLs, the statistical study of the intensity fluctuations was conducted in this work. This was motivated by several works37,38,39,40,41,42,43 demonstrating that the RL transition is a photonic phase-transition, analogous to the paramagnetic to spin-glass phase-transition in disordered magnetic materials44. In RLs this phase-transition is characterized by a Replica Symmetry Breaking (RSB) where intensity fluctuations are highly correlated. However, although RSB has been identified and studied in RLs based on liquids, crystalline powders and optical fibers, the observation of RSB in GCs was not reported.

Differential Scanning Calorimetry (DSC) results for the Nd3+ doped GPM powder, obtained before the annealing process to prepare GCs, are shown in Fig. 1. A typical glass-transition is observed at 605 °C. Additionally, the thermogram shows an exothermic peak with an onset temperature of about 850 °C. This feature is attributed to the growth of crystalline phases within the glass. A temperature, slightly under the onset of crystallization, at 830 °C, was used in this work to induce a controlled crystallization of the samples, with the aim to obtain GCs with different degrees of crystallization. A temperature slightly under 850 °C was used to minimize the adherence of the samples to the crucible used for the thermal annealing.

DSC of GPM (10 wt% of Nd2O3) glass sample at heating rate of 10 °C/min (Tg = glass transition temperature, Tc = crystallization temperature).

Figure 2 exhibits the particle size distribution of Nd3+ doped GPM glass–ceramic powder annealed for 5 h at 830 °C (GPM5) obtained from optical microscopy measurements. The particles average size was \(\approx\) 2 µm, but it was observed a small fraction of larger particles, with dimensions up to \(\approx\) 24 µm. The other samples present similar sizes distribution because the grains were collected after sedimentation of coarser particles mixed with isopropyl alcohol as described in "Optical characterization".

Particle size distribution for the sample GPM5 (glass–ceramic powder annealed for 5 h at 830 °C).

Figure 3 shows X-ray diffractograms obtained for Nd3+ doped glass and GC samples. Notice that the GPM sample, not submitted to annealing, did not present evidence of crystallization. However, samples annealed at 830 °C for 0.5 h and 1 h show sharp crystallization peaks attributed to the growth of nanocrystalline Nd2Ge2O7 phase (PDF: 42-0208). On the other hand, GPM3 powder presents additional peaks related to the presence of MgPb3Ge5O14 (PDF: 39-1265) and GeO2 (PDF: 36-1463). Finally, in the GPM5 powder, peaks related to the growth of PbO (PDF: 3-0600) crystalline phase were also present. Scherrer equation indicated that the crystallite sizes grown in the sample had dimensions of about 50 nm.

X-ray diffraction for Nd3+ doped samples GPMx (x = 0.5, 1, 3, 5). The inset presents the crystallinity degree as a function of annealing time.

The kinetics of the isothermal crystallization of the lead-germanate glass may be described using the Johnson–Mehl–Avrami-Komogorov (JMAK) equations45. In this context, the volume fraction, \(x\left(t\right),\) of a crystal grown from the glass-phase during isothermal reactions can be described by the expression:

where n is the Avrami exponent and K is the crystallization rate constant, expressed as:

where \({E}_{eff}\) is the effective activation energy describing the overall crystallization process, \({K}_{0}\) is the isothermal JMAK parameters, T is the heating temperature, and \({k}_{b}\) is the Boltzmann constant. Both n and K reflect the nucleation and growth mechanisms of the sample.

The plot of the crystallinity degree as a function of time (inset of Fig. 3) presented a typical sigmoidal curve, as expected from the JMAK model. The general equation for the Avrami exponent is \(n=a+bc\), where a is the nucleation index (\(a=0\) for zero nucleation rate during phase change transformation, \(0<a<1\) for a decreasing nucleation rate, a = 1 for a constant nucleation rate, and \(a>1\) for an increasing nucleation rate). The parameter \(b\) is the dimensionality of the crystal grown (b = 1 for 1D, b = 2 for 2D and b = 3 for 3D crystal) and c is the growth index (c = 0.5 for diffusion-controlled growth and c = 1 for interface-controlled growth)42. In the present case, the Avrami exponent was obtained from the fitting of the experimental data to the JMAK model (Eq. 2) and \(n\) was found to be approximately 0.89. In this case, the values for b and c must be 1 and 0.5, respectively, which means that the dimensionality of the crystals were 1D and the growth process was diffusion-controlled. The parameter \(a\) is 0.39 which corresponds to a decreasing nucleation rate as a function of time46. The effective activation energy, \({E}_{eff},\) was not estimated in this work since additional crystallinity degree versus annealing time curves for other temperatures would be necessary. Nevertheless, the investigation of the details concerning the crystals growth kinetics were out of the scope of the present work.

Coherent Back Scattering (CBS) results are shown in Fig. 4a. The transport mean-free-path \(({l}_{t})\) as a function of the annealing times were determined using \({l}_{t}=0.7\frac{\lambda }{2\pi W}\) as described in "Optical characterization" and are shown in Fig. 4b. We observe an expressive reduction of the mean-free-path, as the annealing time increases. Hence, it is evident that crystallites that are grown within the host glass are also acting as scattering media. This is an interesting result from a materials characterization point of view, since our results show that CBS measurements may also be a suitable technique to access the crystallization behavior of transparent solids, such as glasses or polymers.

(a) CBS cones of scattered light intensity (at 1064 nm) as a function of detection angles for all samples prepared. (b) Mean-free-path as a function of annealing time.

Figure 5a shows the absorption cross section of the precursor GPM glass doped with 10 wt% Nd2O3, and Fig. 5b shows the diffuse reflectance spectrum of the glass and GCs powders. Absorption peaks could be observed at about 532 nm, 585 nm, 690 nm, 750 nm, 808 nm, and 880 nm, which are related to the Nd3+ transitions 4I9/2 → {4G7/2, 2G9/2, 2K13/2}, 4I9/2 → {4G5/2, 2G7/2}, 4I9/2 → 4F9/2, 4I9/2 → {4F7/2, 4S3/2}, 4I9/2 → {4F5/2, 2H9/2}, and 4I9/2 → 4F3/2, respectively. Evidencing the incorporation of neodymium in the trivalent form (Nd3+) in the glass and glass–ceramics studied.

(a) Absorption cross section of the precursor GPM glass with neodymium ions. The absorption bands are due to transitions from the ground state to the Nd3+ excited states. (b) Diffuse reflectance spectra of all prepared powders. The spectra are shifted in the vertical axis to prevent overlap.

The PL spectra and the temporal behavior of the samples excited with low pump intensity using a CW laser at 808 nm are shown in Fig. 6. Notice in Fig. 6a that the emission band at ~ 1068 nm (due to transitions 4I9/2 → {4F5/2, 2H9/2}) grows as the annealing time increases. The sample annealed for 5 h showed a PL intensity ~ 500% higher in comparison with the sample that was not subjected to the crystallization process. The temporal behavior of the PL emission at 1068 nm of the glass and GCs is shown in Fig. 6b. The decay lifetime rises from 30 µs for the GPM glass to approximately 120 µs for the GCs samples (see inset in Fig. 6b).

(a) PL spectra of the powders showing the Nd3+ transitions 4F3/2 → 4I9/2, 4F3/2 → 4I11/2 and 4F3/2 → 4I13/2. (b) Temporal behavior of the 1068 nm PL signal, that corresponds to transition (4F3/2 → 4I11/2), for all samples.

The PL enhancement observed in Fig. 6a is correlated to the degree of crystallization of the samples. To verify if the photons mean-free-path was correlated to the PL results, a Monte Carlo simulation was performed. The algorithm used considered that each incident pump photon would follow a 3D random walk inside the sample, with an exponential distribution for the path lengths between successive scattering events, with the mean-free-path within the range of transport lengths determined from the CBS measurements. There were no constraints for the size of the simulation box, except for the plane of incidence of photons, located at z = 0. To simulate an incident photon on the sample, the initial polar and azimuthal angles of the displacement vector of the photon are chosen to be (θ0, φ0) = (0, 0), where θ and φ are the polar and azimuthal angles, respectively. In the model, the pump photon can either be scattered through the sample or be absorbed and give rise to a secondary photon. At each scattering event, the direction of propagation of the photons changed through θ and φ. The absorption probability was set in accordance with the absorption cross section of the sample at 808 nm. The secondary photon was allowed to randomly scatter through the sample with the same initial mean-free-path. The simulation was finished if either the secondary photon was able to reach the surface or if the secondary photon was scattered inside the sample for more than 10,000 times. The simulation was repeated for 5000 pump photons. The ratio between secondary photons that could reach the surface and would be able to be detected (Ne) to those generated within the sample (N0) is shown in Fig. 7, as a function of the transport mean-free-path. The simulation shows that the Ne/N0 ratio is independent of the mean-free-path. Thus, spontaneous emission would not benefit from the reduction of photons mean transport length, that was observed from CBS characterization. For this reason, PL enhancement may be due to the reduction of nonradiative losses, that may be attributed to the incorporation of Nd3+ ions into the structure of crystallites that were grown during the annealing treatments. This is corroborated by XRD results, since the Nd2Ge2O7 phase is present in the glass–ceramics even for samples annealed for only 0.5 h. In addition, samples annealed above 3 h, clearly show small sharp features in the PL spectra, which are overlapped to the broad emission band characteristic of the glassy phase. The emerging sharp bands are also an indication that Nd3+ ions were incorporated into NdO7 and NdO8 polyhedra47 within the structure of the Nd2Ge2O7 crystallites that were grown during the annealing treatment. The rise on the lifetime after the annealing treatment is also consistent to the fact the Nd3+ ions are migrating to specific sites within the crystalline structure of Nd2Ge2O7. In this case, for the parent glass, the average distance among Nd3+ ions is probably, much shorter than in the GCs. Thus, the higher proximity of Nd3+ ions in the glass enhances the ion-ion interactions causing a reduction in the observed lifetime48,49,50.

Monte Carlo Simulation of the ratio between spontaneous photons that can leave the sample to spontaneous photons that are generated within the host material as a function of the mean-free-path.

We observed and characterized in the present work the photonic analogue behavior of the paramagnetic to spin-glass phase in magnetic materials. The photonic-phase-transition, which corresponds to the transition from the spontaneous emission regime (below the laser threshold) to the RL glassy behavior (above the threshold) was characterized as described in "Optical characterization".

To analyze the RL modes correlations the Parisi overlap parameter \({q}_{\alpha \beta }\) of Eq. (3), discussed in "Optical characterization", was calculated from N = 200 acquired PL spectra, obtained for each energy fluence excitation (EFE) used. The probability distribution P(q), \(q= {q}_{\alpha \beta }\), was determined for all samples and for various EFEs.

Figure 8 shows P(q) curves obtained for the GPM5 sample. As can be seen in Fig. 8a \({q}_{\mathrm{max}}\approx 0\) when EFE is below the threshold for RL emission (EFEth). This indicates absence of correlations among the output intensity fluctuations. As the EFE increases, the shape of P(q) changes, as shown in Fig. 8b–f. When EFE \(\cong\) EFEth, \({q}_{\mathrm{max}}\) shifts to near − 1 and + 1, which indicates correlated and anti-correlated fluctuations, respectively. Nevertheless, a significant number of replicas are not correlated as indicated by the large magnitude of P(q = 0) in Fig. 8c. Above the RL threshold the RSB occurs, and the maximum of P(q) consolidates at \(\left|{\mathrm{q}}_{\mathrm{max}}\right|\cong 1\), as shown in Fig. 8d–f. It is worth mentioning that the RSB phenomenon was largely reported for different RLs systems2,37,38,39,40,41,42,43. Yet, it is the first time that RSB phenomenon is reported for a RL based on a glass–ceramic system, for the best of the authors knowledge.

Probability distribution of the Parisi overlap parameter for different EFEs. Below threshold: (a) EFE = 25 μJ/mm2, and (b) EFE = 125 μJ/mm2. In both cases the maximum of P(q) occurs at qmax ≈ 0. (c) On the threshold: EFE = 250 μJ/mm2; the maximum of P(q) occurs at \(\left|{\mathrm{q}}_{\mathrm{max}}\right|\cong 1\). Above the RL threshold (d) EFE = 375 μJ/mm2, (e) EFE = 750 μJ/mm2 and (f) EFE = 6250 μJ/mm2; the maximum of P(q) consolidates at \(\left|{\mathrm{q}}_{\mathrm{max}}\right|\cong 1\).

Figure 9 shows the PL intensity at 1068 nm (black curves) and |\({q}_{\mathrm{max}}|\) (red curves) versus the EFE for all samples. A clear transition from spontaneous emission to RL is observed from the PL intensity and |\({q}_{\mathrm{max}}|\) curves. Notice that the EFEth obtained from the PL vs EFE curves agrees with the EFEth found for |\({q}_{\mathrm{max}}|\) transition. This was also observed for other RLs already reported2,37,38,39,40,41,42,43.

Modulus of the maximum Parisi overlap parameter (red curves) and PL intensity versus EFE at 1068 nm (black curves) for Nd3+ doped GPM samples. (a) GPM; (b) GPM05; (c) GPM1; (d) GPM3; (e) GPM5; (f) EFEth as a function of heat treatment time. The dashed blue line indicates the EFEth for each sample.

Observe that EFEth is reduced from ~ 1 to ~ 0.25 mJ/mm2 when we compare the untreated glassy GPM sample with the one submitted to the crystallization treatment for 5 h. The results of EFEth as a function of annealing time in Fig. 9f, show clearly that the degree of crystallization is also correlated to the RL performance.

Figure 10a shows that the emission intensity at 1068 nm is about 130 times higher when the EFE is increased from 0.1 to 5.0 mJ/mm2, for the GPM5 sample. We did not observe, however, any significant reduction in the width of the emission band. The inhomogeneous broadening is expected for glassy active media and is present in the glass–ceramics because of the residual parent-glass that is still present in the samples, even after heat-treatments. Figure 10b shows the PL dynamics below and above the EFEth for the GPM5 sample. Below EFEth, the lifetime was in the μs range for all samples. On the other hand, for EFE > EFEth, a fast emission was observed, in the nanosecond range following the pump laser pulse, superimposed on the slower signal (in the µs range) due to the spontaneous emission by the ions that are not participating in the stimulated emission process. The temporal behavior shown by the other samples is like the one shown in Fig. 10b.

(a) Emission spectrum and (b) temporal behavior of the signal (at ~ 1068 nm) for the GPM5 sample for EFE below and above EFEth. (c) PL intensity vs EFE for all samples.

In the GPM sample, the RL feedback mechanism is attributed to the micrometric grains. In this case, the glass particles in the powder act simultaneously as scatterers and gain media. Nevertheless, the photon transport length is roughly of the same size as the particles. On the other hand, in the glass–ceramics, crystallites that are grown within the powder grains due to the annealing treatment, may also act as scatterers. Thus, whilst in glass, the photons residence time is governed by light reflections in the particles-air interfaces, in the glass–ceramics, light reflections by the crystallites-glass interfaces may also be relevant. This additional scattering mechanism may increase the photons residence time within the samples reducing the energy fluence threshold for RL observation. In addition, Nd3+ inside the crystallites grown within the glass phase may present sharper emission lines and may not suffer from large inhomogeneous broadening51. Indeed, Fig. 10c shows that the RL slope-efficiency increases for the samples with a higher crystallinity degree, even though it was not observed reduction on the emission linewidth above EFEth. Therefore, despite the inhomogeneous broadening caused by the residual glass, the crystallites act both as scatterers and as more efficient light emitters.

Nd3+ doped GPM glass and glass–ceramics were successfully fabricated, and the crystallization degree of the samples was controlled by a careful annealing treatment at 830 °C for different time intervals. X-ray diffraction analysis showed that the crystallization degree of the samples varied from 0 to 78%, depending upon the heat treatment duration. Nd2Ge2O7 and MgPb3Ge5O14 crystallites were grown in the glass–ceramics because of the annealing treatment. CBS measurements showed that crystallization had, indeed, an impact in the photons mean-free-path that changed from ~ 3.5 µm for the glassy sample to ~ 2.2 µm for the sample heat-treated for 5 h.

Powders consisting of glass and glass–ceramic microparticles were prepared for optical characterization. Approximately 500% of PL enhancement at 1068 nm was observed for samples with a higher crystallinity degree, relative to the glassy sample. We attribute the PL enhancement to the incorporation of Nd3+ ions into the structure of the crystallites grown due to the heat treatment of the powders.

The energy fluence excitation threshold (EFEth) for RL action was smaller for samples with a higher degree of crystallization. For all samples, we observed replica symmetry breaking (RSB) with the Parisi overlap parameter changing from 0 to \(\pm 1\) for excitation above the EFEth. The RL efficiency was enhanced for the samples that were submitted to higher annealing times. The enhanced RL performance for the glass–ceramics in comparison with the glass samples may be attributed to the residence-time growth of photons inside the powder and by the higher quantum efficiency of Nd3+ located in the crystallites because the radiative losses are expected to be smaller than for the ions in the host glass. The optical feedback mechanism for the RL action in the glass–ceramics is attributed to light reflections in the microparticles-air interfaces and the reflections inside the particles with the emitted photons being reflected by the boundaries of the crystallites.

This is the first time, for the best of our knowledge, that the RSB phenomenon is reported for a RL based on rare-earth ions hosted in a glass–ceramics. We claim that investigations of glass–ceramics for random lasers may be a fruitful research field due to the possibility of fine control of the optical feedback by changing the crystallization degree of the samples, which implies that the photons mean-free-path (measured by Coherent Back Scattering) and the RL action in glass–ceramics may be used as tools to investigate the kinetics of crystallization in such materials.

The composition of the glasses prepared was 40GeO2–55PbO-5MgO (wt%)—labeled as GPM; 10 wt% of Nd2O3 (corresponding to 1.87 × 1021 ions/cm3) were added to the final composition. All raw materials used were high purity (> 99.99%). Although it is expected that a large concentration of Nd2O3 in the glass samples causes luminescence concentration quenching (LCQ), we have already observed that the RL performance is enhanced for higher rare-earth ions concentrations, despite the occurrence of LCQ30. The reason is that the dynamic of the RL emission occurs in the nanosecond regime while the PL occurs is in the microsecond range25,30. The glasses were obtained by conventional melt-quenching technique. Reagents were melted at 1200 °C in a platinum crucible for 1 h, and then, quenched in water, at room temperature to prevent crystallization. The resulting GPM glasses were ground using a mortar and pestle to obtain a fine powder. Approximately 18 mg of the GPM powder were submitted to Differential Scanning Calorimetry (DSC) analysis (Labsys Evo, Setaram), to verify the most suitable temperatures for the crystallization process. DSC analysis was conducted in N2 atmosphere (100 mL/min) using an alumina crucible and the heating rate was 20 °C/min.

The GPM powder was separated in five samples labelled GPM, GPM05, GPM1, GPM3, and GPM5. Each sample was submitted to annealing treatment at 830 °C, in air, for different periods of time, as shown in Table 1. The annealing temperature was slightly below the onset of crystallization that were determined from DSC measurements. After the annealing process, each GC sample was subjected to a second pulverization process. Then, the obtained powders were submitted to X-Ray diffraction (XRD) in a Rigaku SmartLab advanced X-ray diffractometer with Cu Kα radiation (λx = 0.154059 nm; 40 kV; 30 mA) to follow the structural changes of the specimens (step size: 0.01°; time per step: 0.3 s). The average size \((D)\) of crystallites was estimated from the width of X-ray peaks according to Scherrer's equation \(D=\frac{K{\lambda }_{x}}{\beta \mathrm{cos}\theta }\)36, where \({\lambda }_{x}\) is the wavelength of X-ray radiation, \(\theta\) is the diffraction angle, \(\beta\) is the width of peak at half of its maximum intensity (FWHM), and \(K\) is a dimensionless shape factor known as the Scherrer constant (in the case of spherical particles \(K\) = 0.94). The percent crystallinity degree \((CD)\) was estimated by the ratio of the crystalline area, \({A}_{C}\), present in the diffractogram of the devitrified glass (glass–ceramics) and the total area, \({A}_{T} (\mathrm{amorphous}+\mathrm{crystalline})\), using the equation \(CD =100\frac{{A}_{C}}{{A}_{T}}\)52.

The samples preparation for the optical experiments was according to the following procedure. Initially, each annealed powder was vigorously mixed with isopropyl alcohol to obtain a homogenous suspension. Then, the suspension was left to rest for 30 s for sedimentation of coarser particles. The supernatant was collected, and after solvent evaporation, a fine powder was obtained.

Sample holders were prepared from silica glass microscope slides. Cavities of 10 × 10 mm2 and about 500 µm of depth were mechanically etched onto the surface of the glass slides, using silica microspheres sandblasting. The samples’ holders were then cleaned with isopropyl alcohol in an ultrasonic bath. Approximately 100 mg of each powder were placed in the cavity of the sample holders and gently pressed, to obtain homogeneous layers of Nd3+ doped GPM glass and GC particles. Other microscope slides were used to cover the cavities containing the samples to avoid powder leakage. Particles size distribution were obtained from optical microscopy images of each sample; the ImageJ open-source software53 was used to analyze the micrographs.

The mean-free-path of photons, for each sample, was determined by Coherent Backscattering Scattering (CBS) measurements. The setup for these experiments is described elsewhere48. The mean transport path of the photons, \({l}_{t}\), could be obtained from the equation \({l}_{t}=0.7\frac{\lambda }{2\pi W}\)54, where \(\uplambda\) is the laser wavelength (1064 nm) and W is the full width at half maximum (FWHM) of the CBS cones.

Diffuse reflectance spectra of the powders were obtained, in the VIS–NIR range, using a spectrophotometer (Ocean Optics), at room temperature. Photoluminescence (PL) experiments were performed at room temperature using two optical sources. For the PL experiments under low pump intensity a CW diode laser operating at 808 nm was used. For experiments with large intensities, it was used an Optical Parametric Oscillator (OPO), OPOTEK INC., Opollete™ HE 532 LD model, pumped by the second harmonic of a Q-switched Nd: YAG laser (7 ns, 20 Hz). The OPO wavelength was tuned to 808 nm, in resonance with the Nd3+ transition 4I9/2 → {4F5/2, 2H9/2}, to optimize the fluorescence signal due to the 4F5/2 → 4I11/2 transition, leading to a RL emission peak at 1068 nm.

The light beam from the OPO was focused on the sample by a 70 mm focal length lens, corresponding to an illuminated area of about 200 µm of diameter. The angle between the perpendicular direction to the sample face and the incident laser beam was 45° and the scattered light emitted from the sample was collected from its front surface and focused into a high-resolution (~ 0.01 nm) spectrometer coupled to a charge-coupled device (CCD). A long-pass filter was positioned at the entrance of the monochromator to remove the scattered light due to the incident laser beam. The output intensity was controlled by a variable neutral density filter, and a reference of the input intensity was monitored by a photodiode coupled to an oscilloscope. PL intensity versus the Energy Fluence Excitation (EFE) curves were obtained by tuning the center of the monochromator grating to the maximum emission intensity (\(\approx\) 1068 nm), and then, 200 spectra were obtained for each pumping condition to plot PL intensity vs EFE curves. The critical energy fluence for the transition between spontaneous emission to RL regime (EFEth) could be determined from the PL intensity versus EFE curves. The PL temporal evolution was determined using the same setup as described above but changing the CCD by an IR photomultiplier coupled to a 2 GHz bandwidth model oscilloscope (Tektronix MSO5204B Mixed Signal). Figure 11 illustrates the experimental setup used to characterize the RL emission.

Setup for random laser characterization.

The photonic phase-transition and the associated phenomenon of RSB were observed and characterized using the PL spectra. The probability density function, P(\({q}_{\alpha \beta }\)), where \({q}_{\alpha \beta }\) is called the overlap parameter, measures the degree of correlation among the RL spectral modes being determined for EFEs below and above the RL threshold2,37,40,41,42,43. The values of \({q}_{\alpha \beta }\) were obtained considering the intensity fluctuations by using the expression37,40,41,42,43.

where \(\alpha\), \(\beta\) = 1, 2,…, N denote the replica labels; the average intensity at the wavelength indexed by k is represented by \(\langle I\rangle \left(k\right)=\sum_{\alpha =1}^{N}{I}_{k}(k)/N\), and the intensity fluctuation is symbolized by \({\Delta }_{\alpha }\left(k\right)={I}_{\alpha }\left(k\right)-\langle I\rangle \left(k\right)\). Each output spectrum is considered a replica, i.e., a copy of the RL system under initial identical experimental conditions. Then, the distribution P(\({q}_{\alpha \beta }\)) can be determined for each value of the EFE as in references37,38,39,40,41,42,43.

The datasets used and/or analysed during the current study are available from the corresponding author on request.

Letokhov, V. S. Generation of light by a scattering medium with negative resonance absorption. Sov. J. Exp. Theor. Phys. 26, 835–840 (1968).

ADS Google Scholar

Gomes, A. S. L., Moura, A. L., de Araújo, C. B. & Raposo, E. P. Recent advances and applications of random lasers and random fiber lasers. Prog. Quantum Electron. 78, 100343 (2021).

Article Google Scholar

Kao, T. S., Hong, Y., Hong, K. & Lu, T. Perovskite random lasers: A tunable coherent light source for emerging applications. Nanotechnology 32, 282001 (2021).

Article ADS CAS Google Scholar

Chen, H. et al. Advances in random fiber lasers and their sensing application. Sensors 20, 6122 (2020).

Article ADS CAS PubMed Central Google Scholar

Cao, H., Chriki, R., Bittner, S., Friesem, A. A. & Davidson, N. Complex lasers with controllable coherence. Nat. Rev. Phys. 1, 156–168 (2019).

Article Google Scholar

Sapienza, R. Determining random lasing action. Nat. Rev. Phys. 1, 690–695 (2019).

Article Google Scholar

Du, X. et al. High-power random distributed feedback fiber laser: From science to application. Ann. Phys. 528, 649–662 (2016).

Article Google Scholar

Wang, Z., Meng, X., Kildishev, A. V., Boltasseva, A. & Shalaev, V. M. Nanolasers enabled by metallic nanoparticles: From spasers to random lasers. Laser Photon. Rev. 11, 1700212 (2017).

Article ADS Google Scholar

Luan, F. et al. Lasing in nanocomposite random media. Nano Today 10, 168–192 (2015).

Article CAS Google Scholar

Gomes, A. S. L. Nanocomposite-based random lasers: a review on basics and applications, in: Nanocomposites for Photonic and Electronic Applications (ed. Kassab, L. R. P., Ribeiro, S. J. L. & Rangel-Rojo, R.) 45–79 (Elsevier, 2020).

Moura, A. L. et al. Nonlinear effects and photonic phase-transitions in Nd3+ doped nanocrystals based random lasers. Appl. Opt. 59, D155–D162 (2020).

Article CAS PubMed Google Scholar

Lawandy, N. M., Balachandran, R. M., Gomes, A. S. L. & Sauvain, E. Laser action in strongly scattering media. Nature 368, 436–438 (1994).

Article ADS Google Scholar

Ambartsumyan, R. V., Basov, N. G., Kryukov, P. G. & Letokhov, V. S. A Laser with a nonresonant feedback. IEEE J. Quantum Electron. 2, 442–446 (1996).

Article ADS Google Scholar

Kitur, J., Zhu, G., Bahoura, M. & Noginov, M. A. Dependence of the random laser behavior on the concentrations of dye and scatterers. J. Opt. 12, 024009 (2010).

Article ADS Google Scholar

Hashem, M. R. & Hameed, M. A. High-purity narrow emission line display of titanium dioxide nanoparticles in laser dyes as random gains. Optik 260, 169058 (2022).

Article ADS CAS Google Scholar

Lü, J., Fan, T. & Chen, G. Random laser action in dye doped nanoporous polymeric film. Opt. Commun. 356, 17–20 (2015).

Article ADS Google Scholar

Lahoz, F. et al. Random lasing in brain tissues. Org. Electron. 75, 105389 (2019).

Article CAS Google Scholar

Moura, A. L. et al. Single bead near-infrared random laser based on silica-gel infiltrated with Rhodamine 640. J. Appl. Phys. 123, 133104 (2018).

Article ADS Google Scholar

Azkargorta, J. et al. Random laser properties of Nd3+ crystal powders. Opt. Express 26, 11787 (2018).

Article ADS CAS PubMed Google Scholar

Moura, A. L., Jerez, V., Maia, L. J. Q., Gomes, A. S. L. & de Araújo, C. B. Multi-wavelength emission through self-induced second-order wave-mixing processes from a Nd3+ doped crystalline powder random laser. Sci. Rep. 5, 13816/1–7 (2015).

Moura, A. L. et al. Random lasing in Nd3+ doped potassium gadolinium tungstate crystal powder. J. Appl. Phys. 117, 083102 (2015).

Article ADS Google Scholar

Moura, A. L., Maia, L. J. Q., Gomes, A. S. L. & de Araújo, C. B. Optimal performance of NdAl3(BO3)4 nanocrystals random lasers. Opt. Mater. 62, 593–596 (2016).

Article ADS CAS Google Scholar

Noginov, M. A., Fowlkes, I., Zhu, G. & Novak, J. Neodymium random lasers operating in different pumping regimes. J. Mod. Opt. 51, 2543–2553 (2004).

Article ADS CAS MATH Google Scholar

Barredo-Zuriarrain, M., Iparraguirre, I., Fernández, J., Azkargorta, J. & Balda, R. Speckle-free near-infrared imaging using a Nd3+ random laser. Laser Phys. Lett. 14, 106201 (2017).

Article ADS Google Scholar

Moura, A. L., Maia, L. J. Q., Jerez, V., Gomes, A. S. L. & de Araújo, C. B. Random laser in Nd:YBO3 nanocrystalline powders presenting luminescence concentration quenching. J. Lumin. 214, 116543 (2019).

Article CAS Google Scholar

Gagné, M. & Kashyap, R. Random fiber Bragg grating Raman fiber laser. Opt. Lett. 39, 2755–2758 (2014).

Article ADS PubMed Google Scholar

Lizárraga, N., Puente, N. P., Chaikina, E. I., Leskova, T. A. & Méndez, E. R. Single-mode Er-doped fiber random laser with distributed Bragg grating feedback. Opt. Express 17, 395–404 (2009).

Article ADS PubMed Google Scholar

Jagannathan, S. et al. Random lasing from optical fibers with phase separated glass cores. Opt. Express 28, 22049–22063 (2020).

Article ADS CAS PubMed Google Scholar

de Oliveira, M. A. S., de Araújo, C. B. & Messaddeq, Y. Upconversion ultraviolet random lasing in Nd3+ doped fluoroindate glass powder. Opt. Express 19, 5620–5626 (2011).

Article ADS PubMed Google Scholar

Câmara, J. G., da Silva, D. M., Kassab, L. R. P., de Araújo, C. B. & Gomes, A. S. L. Random laser emission from neodymium doped zinc tellurite glass-powder presenting luminescence concentration quenching. J. Lumin. 233, 117936 (2021).

Article Google Scholar

Xu, X., Zhang, W., Jin, L., Qiu, J. & Yu, S. F. Random lasing in Eu3+ doped borate glass-ceramic embedded with Ag nanoparticles under direct three-photon excitation. Nanoscale 7, 16246–16250 (2015).

Article ADS CAS PubMed Google Scholar

Li, X. et al. Phase-selective nanocrystallization of NaLnF4 in aluminosilicate glass for random laser and 940 nm LED-excitable upconverted luminescence. Laser Photon. Rev. 12, 1800030 (2018).

Article ADS Google Scholar

Bordon, C. D. S., Magalhães, E. S., da Silva, D. M., Kassab, L. R. P. & de Araújo, C. B. Influence of Al2O3 on the photoluminescence and optical gain performance of Nd3+ doped germanate and tellurite glasses. Opt. Mater. 109, 110342 (2020).

Article CAS Google Scholar

Munasinghe, H. T. et al. Lead germanate glasses and fibers: a practical alternative to tellurite for nonlinear fiber applications. Opt. Mater. Express 3, 1488–1503 (2013).

Article ADS Google Scholar

Abo-Mosallam, H. A. & Mahdy, E. A. The influence of MgO on the crystallization behaviour and properties of SrO-rich phosphosilicate glasses. Ceram. Int. 46, 12009–12014 (2020).

Article CAS Google Scholar

Nouri, M., Tavoosi, M. & Alizadeh, P. Processing and optical properties of transparent GeO2-PbO-MgO-MgF2 glass-ceramics. Ceram. Int. 42, 17524–17529 (2016).

Article CAS Google Scholar

Ghofraniha, N. et al. Experimental evidence of replica symmetry breaking in random lasers. Nat. Commun. 6, 6058 (2015).

Article ADS CAS PubMed Google Scholar

Antenucci, F., Conti, C., Crisanti, A. & Leuzzi, L. General phase diagram of multimodal ordered and disordered lasers in closed and open cavities. Phys. Rev. Lett. 114, 043901 (2015).

Article ADS CAS PubMed Google Scholar

Antenucci, F., Crisanti, A. & Leuzzi, L. The glassy random laser: replica symmetry breaking in the intensity fluctuations of emission spectra. Sci. Rep. 5, 16792 (2015).

Article ADS CAS PubMed PubMed Central Google Scholar

Tommasi, F., Ignesti, E., Lepri, S. & Cavalieri, S. Robustness of replica symmetry breaking phenomenology in random laser. Sci. Rep. 6, 37113 (2016).

Article ADS CAS PubMed PubMed Central Google Scholar

Gomes, A. S. L. et al. Observation of Lévy distribution and replica symmetry breaking in random lasers from a single set of measurements. Sci. Rep. 6, 27987 (2016).

Article ADS CAS PubMed PubMed Central Google Scholar

de Araújo, C. B., Gomes, A. S. L. & Raposo, E. P. Levy statistics and the glassy behavior of light in random fiber lasers. Appl. Sci. 7, 644 (2017).

Article Google Scholar

González, I. R. R. et al. Influence of fifth-order nonlinearities on the statistical fluctuations in emission intensities in a photonic open-cavity complex system. Phys. Rev. A 102, 063515 (2020).

Article ADS Google Scholar

Mézard, M., Parisi, G. & Virasoro, M. A. Spin glass theory and beyond (World Scientific, 1987).

Choi, H. W., Kim, Y. H., Rim, Y. H. & Yang, Y. S. Crystallization kinetics of lithium niobate glass: Determination of the Johnson–Mehl–Avrami–Kolmogorov parameters. Phys. Chem. Chem. Phys. 15, 9940–9946 (2013).

Article CAS PubMed Google Scholar

Zeng, X. et al. Dynamics of supramolecular crystal growth at the liquid-solid interface studied via scanning tunneling microscope and the Avrami equation. J. Phys. Chem. C 125, 10451–10457 (2021).

Article CAS Google Scholar

Vetter, G. & Queyroux, F. Détermination structurale de Nd2Ge2O7. J. Solid State Chem. 45, 293–302 (1982).

Article ADS CAS Google Scholar

García-Solé, J., Bausá, L. E. & Jaque, D. An introduction to the optical spectroscopy of inorganic solids (Wiley, 2005).

Santos, C. C. et al. Spectroscopic properties of Er3+-doped lead phosphate glasses for photonic application. J. Phys. D 43, 025102 (2010).

Article ADS Google Scholar

Jaque, D., Caldiño, U., Romero, J. J. & Solé, J. G. Influence of neodymium concentration on the cw laser properties of Nd doped Ca3Ga2Ge3O12 laser garnet crystal. J. Appl. Phys. 86, 6627–6633 (1999).

Article ADS CAS Google Scholar

Snitzer, E. Glass lasers. Appl. Opt. 5, 1487–1499 (1966).

Article ADS CAS PubMed Google Scholar

Daguano, J. K. M. F. et al. Effect of partial crystallization on the mechanical properties and cytotoxicity of bioactive glass from the 3CaO.P2O5–SiO2–MgO system. J. Mech. Behav. Biomed. Mater. 14, 78–88 (2012).

Schneider, C. A., Rasband, W. S. & Eliceiri, K. W. NIH Image to ImageJ: 25 years of image analysis. Nat. Methods 9, 671–675 (2012).

Article CAS PubMed PubMed Central Google Scholar

Noginov, M. A. Neodymium Random Lasers: Experimental Studies of Stimulated Emission in: Solid-State Random Lasers (ed. William T. Rhodes) 10–47 (Springer, 2005).

Download references

This work was supported by the Brazilian agencies Conselho Nacional de Desenvolvimento Científico e Tecnológico—CNPq (Grant: 431162/2018-2 and the National Institute of Photonics (INCT) program—Grant: 465.763/2014) and the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) for a doctoral scholarship for J. G. Câmara. We also acknowledge the Laboratório de Difratometria de Raios X of the Departamento de Física-UFPE for the X-ray diffraction analysis.

Manoel L. Silva-Neto

Present address: Department of Chemistry, University of Toronto, Toronto, ON, M5S 3H6, Canada

Departamento de Engenharia Elétrica, Escola Politécnica, Universidade de São Paulo, São Paulo, SP, 05508-970, Brazil

Josivanir G. Câmara

Faculdade de Tecnologia de São Paulo, Pça Cel. Fernando Prestes, 30, São Paulo, SP, 01124-060, Brazil

Davinson M. da Silva & Luciana R. P. Kassab

Departamento de Física, Universidade Federal de Pernambuco, Recife, PE, 50670-901, Brazil

Manoel L. Silva-Neto, Guillermo Palacios & Cid B. de Araújo

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

Conceptualization, D.M.d.S. and C.B.d.A.; methodology, J.G.C., D.M.d.S., M.L.S.-N. and G.P.; validation, C.B.d.A.; formal analysis, J.G.C., D.M.d.S., G.P. and C.B.d.A.; investigation, J.G.C., M.L.S.-N.; resources, D.M.d.S., L.R.P.K. and C.B.d.A.; data curation, J.G.C.; writing—original draft preparation, J.G.C. and D.M.d.S.; writing—review and editing, C.B.d.A., D.M.d.S., L.R.P.K.; supervision, C.B.d.A., D.M.d.S.; project administration, D.M.d.S., L.R.P.K. and C.B.d.A.; funding acquisition, D.M.d.S. and L.R.P.K. All authors have read and agreed to the published version of the manuscript.

Correspondence to Davinson M. da Silva.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Câmara, J.G., da Silva, D.M., Kassab, L.R.P. et al. Random lasing and replica symmetry breaking in GeO2-PbO-MgO glass–ceramics doped with neodymium. Sci Rep 12, 19438 (2022). https://doi.org/10.1038/s41598-022-23893-4

Download citation

Received: 05 July 2022

Accepted: 07 November 2022

Published: 14 November 2022

DOI: https://doi.org/10.1038/s41598-022-23893-4

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.