Laser Cooling and Trapping of Rare-Earth-Doped Particles

: This review focuses on optical refrigeration with the anti-Stokes ﬂuorescence of rare-earth (RE)-doped low-phonon micro- and nanocrystals. Contrary to bulk samples, where the thermal energy is contained in internal vibrational modes (phonons), the thermal energy of nanoparticles is contained in both the translational motion and internal vibrational (phonons) modes of the sample. Much theoretical and experimental research is currently devoted to the laser cooling of nanoparticles. In the majority of the related work, only the translational energy of the particles has been suppressed. In review, the latest achievements in hybrid optical refrigeration of RE-doped low-phonon micro-and nanoparticles are presented. Hybrid cooling permits the suppression of not only the translational energy of the RE-doped particles, but also their internal vibrational phonon thermal energy. Laser cooling of nanoparticles is not a simple task. Mie resonances can be used to enhance laser cooling with the anti-Stokes ﬂuorescence of nanoparticles made of low-phonon RE-doped solids. Laser-cooled nanoparticles is a promising tool for fundamental quantum-mechanical studies, nonequilibrium thermodynamics, and precision measurements of forces.


Introduction
Cooling or refrigeration is a physical process permitting the maintaining of a sample at a temperature below that of its surroundings. This technique dates back thousands of years to when people first tried to preserve food. The reliable source of coherent radiation (laser), demonstrated by Theodore Maiman in 1960, has revolutionized the cooling process. Laser cooling (optical refrigeration) and trapping have since been developed.
In late 1950 and the early 1960s, Wolfgang Paul and Hans Dehmelt introduced the first ion traps. Wolfgang Paul with colleagues developed the three-dimensional radiofrequency (RF) trap, known also as the Paul trap, from a linear quadrupole mass filter. It can be considered as a two-dimensional ion trap [1]. Hans Dehmelt and a number of groups working independently proposed a magnetic field for trapping [2][3][4][5][6][7]. Dehmelt named this trap the Penning trap after Frans Penning who, in 1936, proposed the use of a magnetic field to an electrical discharge in order to increase the lifetime of electrons using the confining effect of the field [7].
In 1968, Vladilen Letokhov proposed the use of the dipole force to trap atoms [8]. The electric dipole interaction between the atoms and a standing laser wave was used to form a periodic lattice of potential wells. In this lattice, minima coincided with the nodes or antinodes of the standing laser wave.
In 1975, Doppler cooling was proposed by two teams: Theodor Hansch and Arthur Leonard Schawlow for free atoms [9], and David Wineland and Hans Georg Dehmelt for trapped particles [10]. Doppler cooling occurs when atoms are irradiated with counterpropagating laser beams detuned below resonance. Following the Doppler mechanism, a moving atom absorbs photons into the laser wave counterpropagating its velocity rather than into the copropagating wave; thus, it encounters a force opposed to its velocity and undergoes cooling.

Laser Cooling of Bulk Rare-Earth-Doped Solids
In the majority of cases, a luminescent material interacting with light emits photons with a lower energy compared to the energy of the absorbed photons ( Figure 1a). This effect is called the Stokes shift. The law of conservation of energy applied to these redshifted emitted photons explains heating of the optical material caused by phonon generation (Figure 1a). If the energy of the emitted photons exceeds the energy of absorbed ones, the anti-Stokes shift, accompanied by phonon absorption, takes place in the system (Figure 1b).
In 1988, it was realized that, using multilevel atoms and the spatial variat polarization, a temperature lower than the Doppler cooling limit can be ach This technique is known as sub-Doppler, polarization gradient, or Sisyphus co In 1995, a fifth, low-energy form of matter (Bose-Einstein condensate) w strated by Eric Cornell and Carl Wieman in a gas of rubidium atoms, which w to 170 nK, [12] and by Wolfgang Ketterle in a gas of sodium atoms [13].
All these developed techniques cause laser cooling of the center-of-mass atoms, ions, or particles. Indeed, the energy of molecules in gases almost com long to translational degrees of freedom. Using one of the mentioned methods lational temperature of a nanocrystal can be decreased. This process is known tional cooling. In 1929, Peter Pringsheim proposed the use of anti-Stokes fluor cooling [14]. Following this approach, the internal energy included in the cry vibrations can be removed from a number of solids. Laser trapping and cooling (RE)-doped particles are considered in this review. The theory of laser cooling Stokes fluorescence is briefly discussed in Section 2. Recent achievements in th ment of optical trapping and cooling of micro-and nanoparticles are presented 3.

Laser Cooling of Bulk Rare-Earth-Doped Solids
In the majority of cases, a luminescent material interacting with light emi with a lower energy compared to the energy of the absorbed photons (Figur effect is called the Stokes shift. The law of conservation of energy applied to shifted emitted photons explains heating of the optical material caused by pho ation (Figure 1a). If the energy of the emitted photons exceeds the energy o ones, the anti-Stokes shift, accompanied by phonon absorption, takes place in (Figure 1b). As mentioned above, in 1929, Pringsheim proposed the use of anti-Stok cence for optical refrigeration [14]. In 1995, laser cooling with anti-Stokes fluore reached for RE-doped solids [15]. In this first proof-of-principle experiment, an (Yb 3+ )-doped fluorozirconate ZrF4-BaF2-LaF3-AlF3-NaF-PbF2 (ZBLANP) sa cooled down to only 0.3 K below the environmental temperature. As mentioned above, in 1929, Pringsheim proposed the use of anti-Stokes fluorescence for optical refrigeration [14]. In 1995, laser cooling with anti-Stokes fluorescence was reached for RE-doped solids [15]. In this first proof-of-principle experiment, an ytterbium (Yb 3+ )-doped fluorozirconate ZrF 4 -BaF 2 -LaF 3 -AlF 3 -NaF-PbF 2 (ZBLANP) sample was cooled down to only 0.3 K below the environmental temperature.
Let us consider in detail a laser cooling cycle in a Yb 3+ -doped sample. As one can see in Figure 2, it includes photon absorption at the frequency ν P , electron thermalization in the excited 2 F5 /2 manifold accompanied by phonon absorption, and spontaneous photon emission from the excited 2 F5 /2 manifold to the ground 2 F7 /2 manifold. Let us consider in detail a laser cooling cycle in a Yb 3+ -doped sample. As one ca in Figure 2, it includes photon absorption at the frequency , electron thermalizati the excited 2 F5 ∕ 2 manifold accompanied by phonon absorption, and spontaneous ph emission from the excited 2 F5 ∕ 2 manifold to the ground 2 F7 ∕ 2 manifold. All photon absorption-emission cycles where the energy of the radiated photo ceeds the energy of the absorbed one (anti-Stokes fluorescence) result in cooling, a energy difference has to be compensated by phonon absorption (Figure 2). All the ph absorption-emission cycles where the energy of the absorbed photons exceeds the en of the radiated ones (Stokes fluorescence) causes heat generation in the sample, a energy difference in these cycles has to be compensated by phonon generation (Figu Spontaneous emission is characterized by the mean fluorescence frequency, where Φ( ) is the fluorescence flux density, which can be obtained experimentally = ⁄ . Here, is the mean fluorescence wavelength and is the speed of lig vacuum. If one wants to cool a sample, the pump frequency, , is bound to be sm than the mean fluorescence frequency, , and all radiated anti-Stokes photons must the sample (not to be reabsorbed by Yb 3+ ions). Indeed, these anti-Stokes reabsorbed tons serve as the source of heat generation in the sample, which is undesirable for co The efficiency of the cooling cycle is the difference between the energy of the anti-S emitted photon and the energy of the pump photon normalized by the energy of the p photon that is The electrons excited to the 2 F5 ∕ 2 manifold can also decay nonradiatively to the gr 2 F7 ∕ 2 manifold, generating phonons. The "competition" between the radiative and n diative decays of the 2 F5 ∕ 2 electrons can be quantified by the quantum efficiency, are the radiative and nonradiative decay rate spectively, with and being the radiative and nonradiative lifetimes of the ex manifold, respectively. If the quantum efficiency of the anti-Stokes generation proc very high ( ≥ 99%), laser cooling caused by the removal of phonon energy with Stoke photons from the material takes place in the system, as the probability of nonr tive decay with phonon generation from the excited 2 F5 ∕ 2 to the ground 2 F7 ∕ 2 manif very low [16]. Only such samples are suitable for laser cooling.
To simulate cooling in a RE-doped sample, one can use a two-level model deve in [17,18]: All photon absorption-emission cycles where the energy of the radiated photon exceeds the energy of the absorbed one (anti-Stokes fluorescence) result in cooling, as this energy difference has to be compensated by phonon absorption (Figure 2). All the photon absorption-emission cycles where the energy of the absorbed photons exceeds the energy of the radiated ones (Stokes fluorescence) causes heat generation in the sample, as the energy difference in these cycles has to be compensated by phonon generation (Figure 2). Spontaneous emission is characterized by the mean fluorescence frequency, where Φ(ν) is the fluorescence flux density, which can be obtained experimentally, and ν F = c/λ F . Here, λ F is the mean fluorescence wavelength and c is the speed of light in vacuum. If one wants to cool a sample, the pump frequency, ν P , is bound to be smaller than the mean fluorescence frequency, ν F , and all radiated anti-Stokes photons must leave the sample (not to be reabsorbed by Yb 3+ ions). Indeed, these anti-Stokes reabsorbed photons serve as the source of heat generation in the sample, which is undesirable for cooling. The efficiency of the cooling cycle is the difference between the energy of the anti-Stokes emitted photon and the energy of the pump photon normalized by the energy of the pump photon that is The electrons excited to the 2 F5 /2 manifold can also decay nonradiatively to the ground 2 F7 /2 manifold, generating phonons. The "competition" between the radiative and nonradiative decays of the 2 F5 /2 electrons can be quantified by the quantum efficiency, where W r = 1/τ r and W nr = 1/τ rnr are the radiative and nonradiative decay rates, respectively, with τ r and τ nr being the radiative and nonradiative lifetimes of the excited manifold, respectively. If the quantum efficiency of the anti-Stokes generation process is very high (η q ≥ 99%), laser cooling caused by the removal of phonon energy with anti-Stoke photons from the material takes place in the system, as the probability of nonradiative decay with phonon generation from the excited 2 F5 /2 to the ground 2 F7 /2 manifold is very low [16]. Only such samples are suitable for laser cooling.
To simulate cooling in a RE-doped sample, one can use a two-level model developed in [17,18]: where N T is the RE ions density in the sample, and N 0 and N 1 are the population density in the ground and excited manifolds, respectively. I p is the pump intensity at the frequency ν p . σ a,e (ν p ) is the absorption (a) and emission (e) cross-sections at the pump frequency ν p . The cooling power density generated in the sample, P cool , is the difference between the absorbed power density and radiated power density.
The cooling power generated in a sample of volume V when the system is in the steady-state regime (dN 1 /dt = 0) is [19] where I S = hν P / η q τ r σ a (ν P ) is the saturation intensity and ν F = 2η q − 1 ν F ≈ η q ν F is the effective mean fluorescence frequency. If the sample is placed on a low-contact and low-thermal conductivity support in a vacuum chamber, only radiative heat load may take place in the system. The final equilibrium temperature of the sample, T S , can be estimated by the Stefan-Boltzmann law, where S is the total surface area of the sample, T r is room temperature, and ε is the sample emissivity.
The resonant absorption coefficient, α r (ν, T), can be estimated as a product of the absorption cross-section σ a (ν, T) and the thermal population density, N 0 (E 0i , T), of the initial state with the energy E 0ii that is α r (ν, T) = σ a (ν, T) · N 0 (E 0i , T). Following the Boltzmann distribution, the thermal population density of the initial state is N 0 (E 0i , T) ∼ exp(−E 0i /κ B T). The absorption cross-section is proportional to the Voigt spectral line σ a (ν, T) ∼ g ν (ν, T), which is a result of the convolution of two broadening mechanisms: one that causes a Gaussian broadening and the other that produces a Lorentzian profile. Consequently, the resonant absorption coefficient can be presented as The inhomogeneous broadening is dominant for glass hosts as a result of structural disorder in glasses. It almost does not change with temperature, that is, g ν (ν, T) is almost constant and In glasses, the peak absorption coefficient α r decreases exponentially with temperature. In crystal hosts, homogeneous broadening is dominant. The homogeneous linewidth scales as ∼ T 2 and is associated with the two-phonon Raman scattering in the temperature range of 40-300 K. The peak absorption coefficient is proportional to Appl. Sci. 2022, 12, 3777

of 15
As one can see in relations (9) and (10), when the temperature decreases, the peak absorption coefficient α r in crystals is higher than in glasses. Crystal hosts are the best choice for laser cooling.

Optical Refrigeration of Rare-Earth-Doped Particles
In the majority of cases, bulk, cm− and mm− sized RE-doped low-phonon samples were laser-cooled. Optical refrigeration of RE-doped samples with different shapes and sizes was investigated in [27].
Three different approaches can be used to calculate the optical force applied to a particle. These approaches depend on the size of the particles compared to the wavelength of the laser beam. If the radius r of the particle is much larger than the light wavelength λ (r > 20λ), diffractive effects can be neglected. In this case, geometrical optics (ray optics) can be used for the optical force calculation. If the particle radius is comparable to the light wavelength ( r ∼ λ), the wave's scattering should be considered. If the particle radius is much smaller than the light wavelength (r λ), dipole approximation can be used.

Laser Cooling of Rare-Earth-Doped Mie Resonant Samples
In 1908, Gustav Mie developed the theory of the plane-wave scattering by spheres of arbitrary size (Mie theory) [28]. Following Mie theory, a sphere can support eigenmodes (the Mie resonant modes) if the radius of the sphere and the pump wavelength are properly arranged. Different eigenmodes have different electric field intensity distributions inside the sphere. As noted above, if a sphere radius is much smaller than the wavelength of the incident light (r λ), the responses of the sphere to the incident light can be represented by an electric dipole and a magnetic dipole (the Rayleigh or dipole approximation). According to Mie theory, the pump plane wave can be expanded in spherical harmonics as where E 0 is the amplitude of the incident field, and → M o1n and → N e1n are the vector spherical harmonics, where r = 2πn Yb r/λ P and r is the radial distance. The index n = 1 stands for dipole, n = 2 for quadrupole, etc., j n are the spherical Bessel functions, and π n and τ n are the angle-dependent functions defined as and where P 1 n is the associated Legendre function of the first kind of degree n and the first order. The field inside the sphere is where the coefficients c n and d n are functions of the normalized radius of the sphere, x = 2πR/λ P , with R being the radius of the sphere, and where h (1) n is the spherical Hankel functions, m = n Yb The c n coefficients correspond to the magnetic multipole modes and the d n coefficients correspond to the electric multipole modes. If, for some values of R, the denominator of coefficient c n or the denominator of d n vanishes, the sphere resonates and scatters, exhibiting a peak. The radial distribution of the absolute-square of the electric field averaged over the polar angle, θ, and azimuthal angle, ϕ, of spherical coordinates looks like where m n and n n are functions of the radial distance, r, and are given as where r = 2πn Yb r/λ P . The cooling power generated in the sphere can be calculated as the difference between the pump power absorbed by the Yb 3+ ions in the sphere and the power radiated with anti-Stokes fluorescence at the mean fluorescence wavelength. The field intensity given by relation (19) can then be integrated over the nanosphere and then inserted into relation (6) as which corresponds to the cooling power generated in the Mie-resonant nanosphere. Laser cooling of RE-doped low-phonon microparticles enhanced with Mie resonant modes was considered in [29]. By way of example, an Yb 3+ :YAG nanosphere pumped with the plane electromagnetic wave propagating in vacuum at the pump wavelength λ P = 1030 nm located in the long wavelength tail of the Yb 3+ absorption spectrum was considered ( Figure 3). The mean fluorescence wavelength of the spontaneous radiation λ F < λ P . A bulk Yb 3+ :YAG sample was laser-cooled in [26].
A sample with an Yb 3+ ion concentration of~10% was investigated in [29]. This concentration is smaller than the critical concentration of Yb 3+ :YAG (17%). Such samples are free from cooperative effects such as resonant-radiative transfer (reabsorption) and resonant nonradiative transfer, which deteriorate laser cooling [19]. Indeed, these cooperative effects result in the excited energy migration from the ion to ion inside the sample. This migrating energy can excite impurities in the sample and decay nonradiatively, causing heat generation in the sample. The pump power considered in [29] was I P = 0.5 · 10 −4 W/µm 2 . The cooling power density generated in the nanosphere is equal to where V = 4 3 πR 3 is the volume of the sphere. R is the radius of the sphere. The cooling power density, ρ cool , for nanospheres with different radii in the range between 150 nm and 650 nm was estimated using relations (22) and (23). The results are plotted in Figure 4. It changes considerably on this range and has several peaks. The radii of the nanospheres corresponding to these peaks are very close to the resonant radii of the Mie resonant modes including the magnetic dipole (MD) modes, the magnetic quadrupole (MQ) modes, the electric dipole (ED) modes, and the electric quadrupole (EQ) modes.
Appl. Sci. 2022, 11, x FOR PEER REVIEW Figure 3. System under investigation. Here, and are the pump and mean fluoresce lengths, respectively [29].
A sample with an Yb 3+ ion concentration of ~10% was investigated in [29]. T centration is smaller than the critical concentration of Yb 3+ :YAG (17%). Such sam free from cooperative effects such as resonant-radiative transfer (reabsorption) a nant nonradiative transfer, which deteriorate laser cooling [19]. Indeed, these coo effects result in the excited energy migration from the ion to ion inside the sam migrating energy can excite impurities in the sample and decay nonradiatively heat generation in the sample. The pump power considered in [29] was 10 −4 W μm 2 ⁄ . The cooling power density generated in the nanosphere is equal to 3 is the volume of the sphere. is the radius of the sphere. The power density, , for nanospheres with different radii in the range between and 650 nm was estimated using relations (22) and (23). The results are plotted 4. It changes considerably on this range and has several peaks. The radii of t spheres corresponding to these peaks are very close to the resonant radii of the M nant modes including the magnetic dipole (MD) modes, the magnetic quadrupo modes, the electric dipole (ED) modes, and the electric quadrupole (EQ) modes.  A sample with an Yb 3+ ion concentration of ~10% was investigated in [29]. T centration is smaller than the critical concentration of Yb 3+ :YAG (17%). Such sam free from cooperative effects such as resonant-radiative transfer (reabsorption) a nant nonradiative transfer, which deteriorate laser cooling [19]. Indeed, these coo effects result in the excited energy migration from the ion to ion inside the samp migrating energy can excite impurities in the sample and decay nonradiatively, heat generation in the sample. The pump power considered in [29] was 10 −4 W μm 2 ⁄ . The cooling power density generated in the nanosphere is equal to 3 is the volume of the sphere. is the radius of the sphere. The power density, , for nanospheres with different radii in the range between and 650 nm was estimated using relations (22) and (23). The results are plotted in 4. It changes considerably on this range and has several peaks. The radii of th spheres corresponding to these peaks are very close to the resonant radii of the M nant modes including the magnetic dipole (MD) modes, the magnetic quadrupo modes, the electric dipole (ED) modes, and the electric quadrupole (EQ) modes. With the knowledge of the cooling power density, , it is possible to estim temperature of the nanospheres placed on a low-contact and low-thermal-cond support in a vacuum chamber, using Equation (7). These results are presented in With the knowledge of the cooling power density, ρ cool , it is possible to estimate the temperature of the nanospheres placed on a low-contact and low-thermal-conductivity support in a vacuum chamber, using Equation (7). These results are presented in Figure 5. As one might expect, the radii of the nanospheres for which the temperature reaches minima almost perfectly coincide with those associated with the highest cooling power density ( Figure 4).
Appl. Sci. 2022, 12, 3777 8 of 15 5. As one might expect, the radii of the nanospheres for which the temperature minima almost perfectly coincide with those associated with the highest coolin density (Figure 4). In [30], a nanoscale semiconductor optomechanical resonator (CdSNR) w cooled with a Yb 3+ :YLF microcrystal attached to it. The Yb 3+ :YLF microcrystal wa thermally synthesized. The cooling process with a microscopic cooler provided lo ing. This approach permits a thermal steady temperature state for nanoscale devi achieved. A CdSNR was placed at the end of a silicon substrate. A 10% Yb 3+ :YLF was placed at the free end of the CdSNR cantilever. The 10% Yb 3+ :YLF crystal wa thermally grown. The silicon substrate was located inside a cryostat chamber, wh pumped to ~10 −4 torr ( Figure 6). The temperature of the Yb 3+ :YLF placed at the end of the cantilever was m with differential luminescence thermometry (DLT) [31,32]. Emission from differe tal-field levels of Yb 3+ ions was analyzed using a Boltzmann distribution. At an power 0 = 40.1 mW and spot radius 0 = 1.15 μm, a temperature drop of 23 low room temperature was observed.
A cantilever of length , width , and thickness was considered for co simulations. The Yb 3+ :YLF crystal of = 6 μm, = 7.5 μm, and = 6 approximated as a cuboid. The volume and aspect ratio were almost the same a crystal used in experiments (Figure 7). The thermal conductivity of the YLF (~6 Wm −1 K −1 ) was relatively large. It resulted in an almost uniform temperatur In [30], a nanoscale semiconductor optomechanical resonator (CdSNR) was lasercooled with a Yb 3+ :YLF microcrystal attached to it. The Yb 3+ :YLF microcrystal was hydrothermally synthesized. The cooling process with a microscopic cooler provided local cooling. This approach permits a thermal steady temperature state for nanoscale devices to be achieved. A CdSNR was placed at the end of a silicon substrate. A 10% Yb 3+ :YLF crystal was placed at the free end of the CdSNR cantilever. The 10% Yb 3+ :YLF crystal was hydrothermally grown. The silicon substrate was located inside a cryostat chamber, which was pumped to ∼ 10 −4 torr ( Figure 6).
Appl. Sci. 2022, 11, x FOR PEER REVIEW 8 of 15 5. As one might expect, the radii of the nanospheres for which the temperature reaches minima almost perfectly coincide with those associated with the highest cooling power density ( Figure 4). In [30], a nanoscale semiconductor optomechanical resonator (CdSNR) was lasercooled with a Yb 3+ :YLF microcrystal attached to it. The Yb 3+ :YLF microcrystal was hydrothermally synthesized. The cooling process with a microscopic cooler provided local cooling. This approach permits a thermal steady temperature state for nanoscale devices to be achieved. A CdSNR was placed at the end of a silicon substrate. A 10% Yb 3+ :YLF crystal was placed at the free end of the CdSNR cantilever. The 10% Yb 3+ :YLF crystal was hydrothermally grown. The silicon substrate was located inside a cryostat chamber, which was pumped to ~10 −4 torr ( Figure 6). The temperature of the Yb 3+ :YLF placed at the end of the cantilever was measured with differential luminescence thermometry (DLT) [31,32]. Emission from different crystal-field levels of Yb 3+ ions was analyzed using a Boltzmann distribution. At an incident power 0 = 40.1 mW and spot radius 0 = 1.15 μm, a temperature drop of 23.6 K below room temperature was observed.
A cantilever of length , width , and thickness was considered for computer simulations. The Yb 3+ :YLF crystal of = 6 μm, = 7.5 μm, and = 6 μm was approximated as a cuboid. The volume and aspect ratio were almost the same as for the crystal used in experiments (Figure 7). The thermal conductivity of the YLF crystal (~6 Wm −1 K −1 ) was relatively large. It resulted in an almost uniform temperature distribution in the crystal. The temperature of the Yb 3+ :YLF placed at the end of the cantilever was measured with differential luminescence thermometry (DLT) [31,32]. Emission from different crystalfield levels of Yb 3+ ions was analyzed using a Boltzmann distribution. At an incident power P 0 = 40.1 mW and spot radius w 0 = 1.15 µm, a temperature drop of 23.6 K below room temperature was observed.
A cantilever of length L, width W, and thickness H was considered for computer simulations. The Yb 3+ :YLF crystal of H YLF = 6 µm, L YLF = 7.5 µm, and W YLF = 6 µm was approximated as a cuboid. The volume and aspect ratio were almost the same as for the crystal used in experiments (Figure 7). The thermal conductivity of the YLF crystal ∼ 6 Wm −1 K −1 was relatively large. It resulted in an almost uniform temperature distribution in the crystal. Appl. Sci. 2022, 11, x FOR PEER REVIEW 9 of 15 Figure 7. The geometry of the system [30].
As the cantilever was placed in the vacuum chamber, conduction and convection were absent. The temperature at the silicon/CdS interface (x = 0) was equal to the cryostat temperature. The authors of the paper assumed that all of the heat generated in the Yb 3+ :YLF crystal was transferred to or from the CdSNR across the interface at = . They calculated a cooling power for the sample with = 150 nm , = 2.5 μm , and = 53 μm and estimated that the coldest point in the cantilever was between 26 K and 58 K below room temperature. The calculated cooling power was smaller than the experimentally obtained cooling power. The discrepancy can be explained in the following way: the size of the YLF microcrystals corresponds to the Mie scattering regime. Internal optical fields may be enhanced considerably by cavity resonances.
In 2021, 10% Yb 3+ -doped LiLuF4 (LLF) microcrystals were grown with a hydrothermal process [33]. They were transferred to electron-transparent Si3N4 windows irradiated by a near-infrared laser for solid-state laser refrigeration of windows in vacuum (10 −3 torr) (Figure 8). Analysis of the Yb 3+ luminescence spectrum at an irradiance of 0.854 MWcm −2 revealed laser cooling by more than 20.4 K below room temperature. Using computer simulations, bulk cooling of the Si3N4 window by ~15 K was demonstrated. The silicon substrate was at room temperature. In 2022, 10%Yb 3+ :LiLuF4 (Yb 3+ :LLF) crystals synthesized using a safe and scalable hydrothermal method were laser-cooled more than 15 K below room temperature in air and 5 K in deionized water. They were pumped with a 1020 nm diode laser [34].

Laser Trapping and Cooling Rare-Earth-Doped Particles
Consider RE-doped particles under both translational and internal laser cooling. Currently available methods for nanoparticle cooling allow translational temperatures of a few hundred microkelvin to be achieved [35]. Consider new mechanisms of optical cooling that allow us to reduce the translational temperature of a sample concurrently with its internal cooling, which tends to decrease the internal temperature of the sample.
Here, 10% Yb 3+ -doped β-NaYF4 nanowires in a fluid medium (D2O) were lasertrapped and cooled with a near-infrared single-beam at the wavelength 1020 nm (inset in Figure 9) [36]. β-NaYF4 nanowires were prepared using a low-cost, scalable hydrothermal synthesis [37]. A nanowire can be considered as a cylinder (height ~ 2000 nm) with an hexagonal cross-section (edge ~ 255 nm). As one can see in Figure 9a, undoped (0% Yb 3+ ) β-NaYF4 nanowires under identical conditions (an irradiance of 73 MWcm −2 at the As the cantilever was placed in the vacuum chamber, conduction and convection were absent. The temperature at the silicon/CdS interface (x = 0) was equal to the cryostat temperature. The authors of the paper assumed that all of the heat generated in the Yb 3+ :YLF crystal was transferred to or from the CdSNR across the interface at x = L c . They calculated a cooling power for the sample with H = 150 nm, W = 2.5 µm, and L c = 53 µm and estimated that the coldest point in the cantilever was between 26 K and 58 K below room temperature. The calculated cooling power was smaller than the experimentally obtained cooling power. The discrepancy can be explained in the following way: the size of the YLF microcrystals corresponds to the Mie scattering regime. Internal optical fields may be enhanced considerably by cavity resonances.
As the cantilever was placed in the vacuum chamber, conduction and convection were absent. The temperature at the silicon/CdS interface (x = 0) was equal to the cryostat temperature. The authors of the paper assumed that all of the heat generated in the Yb 3+ :YLF crystal was transferred to or from the CdSNR across the interface at = . They calculated a cooling power for the sample with = 150 nm , = 2.5 μm , and = 53 μm and estimated that the coldest point in the cantilever was between 26 K and 58 K below room temperature. The calculated cooling power was smaller than the experimentally obtained cooling power. The discrepancy can be explained in the following way: the size of the YLF microcrystals corresponds to the Mie scattering regime. Internal optical fields may be enhanced considerably by cavity resonances.
In 2021, 10% Yb 3+ -doped LiLuF4 (LLF) microcrystals were grown with a hydrothermal process [33]. They were transferred to electron-transparent Si3N4 windows irradiated by a near-infrared laser for solid-state laser refrigeration of windows in vacuum (10 −3 torr) (Figure 8). Analysis of the Yb 3+ luminescence spectrum at an irradiance of 0.854 MWcm −2 revealed laser cooling by more than 20.4 K below room temperature. Using computer simulations, bulk cooling of the Si3N4 window by ~15 K was demonstrated. The silicon substrate was at room temperature. In 2022, 10%Yb 3+ :LiLuF4 (Yb 3+ :LLF) crystals synthesized using a safe and scalable hydrothermal method were laser-cooled more than 15 K below room temperature in air and 5 K in deionized water. They were pumped with a 1020 nm diode laser [34].

Laser Trapping and Cooling Rare-Earth-Doped Particles
Consider RE-doped particles under both translational and internal laser cooling. Currently available methods for nanoparticle cooling allow translational temperatures of a few hundred microkelvin to be achieved [35]. Consider new mechanisms of optical cooling that allow us to reduce the translational temperature of a sample concurrently with its internal cooling, which tends to decrease the internal temperature of the sample.
Here, 10% Yb 3+ -doped β-NaYF4 nanowires in a fluid medium (D2O) were lasertrapped and cooled with a near-infrared single-beam at the wavelength 1020 nm (inset in Figure 9) [36]. β-NaYF4 nanowires were prepared using a low-cost, scalable hydrothermal synthesis [37]. A nanowire can be considered as a cylinder (height ~ 2000 nm) with an hexagonal cross-section (edge ~ 255 nm). As one can see in Figure 9a, undoped (0% Yb 3+ ) β-NaYF4 nanowires under identical conditions (an irradiance of 73 MWcm −2 at the Analysis of the Yb 3+ luminescence spectrum at an irradiance of 0.854 MWcm −2 revealed laser cooling by more than 20.4 K below room temperature. Using computer simulations, bulk cooling of the Si 3 N 4 window by ∼ 15 K was demonstrated. The silicon substrate was at room temperature. In 2022, 10%Yb 3+ :LiLuF 4 (Yb 3+ :LLF) crystals synthesized using a safe and scalable hydrothermal method were laser-cooled more than 15 K below room temperature in air and 5 K in deionized water. They were pumped with a 1020 nm diode laser [34].

Laser Trapping and Cooling Rare-Earth-Doped Particles
Consider RE-doped particles under both translational and internal laser cooling. Currently available methods for nanoparticle cooling allow translational temperatures of a few hundred microkelvin to be achieved [35]. Consider new mechanisms of optical cooling that allow us to reduce the translational temperature of a sample concurrently with its internal cooling, which tends to decrease the internal temperature of the sample.
In [38], a similar approach based on single-beam laser trapping at the wavelen 1020 nm was used for optical refrigeration of Yb 3+ :YLF nanocrystals (Figure 10a) p pared using the same method as in [36]: low-cost, scalable hydrothermal synthesis [3 Using cold Brownian motion analysis, the temperature of an Yb 3+ :YLF nanocrystal i fluid medium (D2O) was estimated for different irradiances (Figure 10b). As one can see in Figure 10b, trapping at the wavelength 1064 nm gives rise to he ing of a sample. The wavelength = 1064 nm is energetically insufficient to be sorbed by Yb 3+ ions and subsequently cannot initiate upconversion laser cooling. At trapping wavelength of 1020 nm, 10%Yb 3+ :YLF crystals were observed to cool from 19 at a 5.9 MWcm −2 trapping irradiance to 4 ∘ C at a 25.5 MWcm −2 trapping irradiance (F ure 10b).
In the cases being considered, laser refrigeration of RE-doped particles has be demonstrated in optical tweezers in liquid. It was analyzed by measuring cold Brown motion. Consider optical refrigeration of RE-doped levitated particles optically trapp in a low-pressure environment. The ability to arrange levitated particles in vacuum i promising instrument for understanding macroscopic quantum mechanics [32].
The changes in fluorescence of the nanocrystal as a function of laser trap wavelen were used to measure the temperature of the particle. Indeed, fluorescence spectra directly related to the relative populations of sublevels of the excited 2 F5 ∕ 2 manifold via Boltzmann distribution. The changes in the damping of the motion of the nanocrystal w the temperature were also used for temperature measurements. When the optically refrigerated particle is levitated in gas, cooling of the gas mo In [38], a similar approach based on single-beam laser trapping at the wavelength 1020 nm was used for optical refrigeration of Yb 3+ :YLF nanocrystals (Figure 10a) prepared using the same method as in [36]: low-cost, scalable hydrothermal synthesis [37]. Using cold Brownian motion analysis, the temperature of an Yb 3+ :YLF nanocrystal in a fluid medium (D 2 O) was estimated for different irradiances (Figure 10b). l. Sci. 2022, 11, x FOR PEER REVIEW 10 o pump wavelength 1020 nm) undergo heating by 6 ∘ C above the ambient temperature an irradiance of 73 MWcm −2 , nanowires were observed to refrigerate by > 9 ∘ C be ambient temperature (Figure 9b).
In [38], a similar approach based on single-beam laser trapping at the wavelen 1020 nm was used for optical refrigeration of Yb 3+ :YLF nanocrystals (Figure 10a) p pared using the same method as in [36]: low-cost, scalable hydrothermal synthesis [ Using cold Brownian motion analysis, the temperature of an Yb 3+ :YLF nanocrystal i fluid medium (D2O) was estimated for different irradiances (Figure 10b). As one can see in Figure 10b, trapping at the wavelength 1064 nm gives rise to h ing of a sample. The wavelength = 1064 nm is energetically insufficient to be sorbed by Yb 3+ ions and subsequently cannot initiate upconversion laser cooling. At trapping wavelength of 1020 nm, 10%Yb 3+ :YLF crystals were observed to cool from 1 at a 5.9 MWcm −2 trapping irradiance to 4 ∘ C at a 25.5 MWcm −2 trapping irradiance (F ure 10b).
In the cases being considered, laser refrigeration of RE-doped particles has b demonstrated in optical tweezers in liquid. It was analyzed by measuring cold Brown motion. Consider optical refrigeration of RE-doped levitated particles optically trapp in a low-pressure environment. The ability to arrange levitated particles in vacuum promising instrument for understanding macroscopic quantum mechanics [32].
The changes in fluorescence of the nanocrystal as a function of laser trap wavelen were used to measure the temperature of the particle. Indeed, fluorescence spectra directly related to the relative populations of sublevels of the excited 2 F5 ∕ 2 manifold via Boltzmann distribution. The changes in the damping of the motion of the nanocrystal w the temperature were also used for temperature measurements. When the optically refrigerated particle is levitated in gas, cooling of the gas m cules surrounding the particle takes place in the system. A modification in the local As one can see in Figure 10b, trapping at the wavelength 1064 nm gives rise to heating of a sample. The wavelength λ p = 1064 nm is energetically insufficient to be absorbed by Yb 3+ ions and subsequently cannot initiate upconversion laser cooling. At the trapping wavelength of 1020 nm, 10%Yb 3+ :YLF crystals were observed to cool from 19 • C at a 5.9 MWcm −2 trapping irradiance to 4 • C at a 25.5 MWcm −2 trapping irradiance (Figure 10b).
In the cases being considered, laser refrigeration of RE-doped particles has been demonstrated in optical tweezers in liquid. It was analyzed by measuring cold Brownian motion. Consider optical refrigeration of RE-doped levitated particles optically trapped in a low-pressure environment. The ability to arrange levitated particles in vacuum is a promising instrument for understanding macroscopic quantum mechanics [32].
The changes in fluorescence of the nanocrystal as a function of laser trap wavelength were used to measure the temperature of the particle. Indeed, fluorescence spectra are directly related to the relative populations of sublevels of the excited 2 F5 /2 manifold via the Boltzmann distribution. The changes in the damping of the motion of the nanocrystal with the temperature were also used for temperature measurements. When the optically refrigerated particle is levitated in gas, cooling of the gas molecules surrounding the particle takes place in the system. A modification in the local gas viscosity also takes place in the system. The damping rate of the translational motion of a particle undergoes changes as a result of hot or cold Brownian motion of the particle. The effect was used to measure particle temperature. Figure 11 shows the temperatures of 13 particles trapped using laser light at 1031 nm and 1064 nm at a pressure of 600 mbar The trapping power was 200 mW. Some changes between the temperatures of particles trapped by 1031 nm light was observed. The average internal temperature was ∼ 167 K. The lowest temperature reached in this experiment was ∼ 130 K. As one can see in Figure 11, in the beginning of the experiment, particle 2 was trapped using a 1031 nm laser beam. The temperature reached was 175 K. Later, the laser wavelength was changed to 1064 nm. For this wavelength, the temperature reached was 255 K (Figure 11). For this particle, the temperature was 175 K when trapped with 1031 nm light and 255 K when trapped using 1064 nm light ( Figure 11). Using (6) and (7), one can show that lower temperatures can be achieved with a laser at 1020 nm. Lower pressure and improved crystal purity can also help to reach lower temperatures.
Appl. Sci. 2022, 11, x FOR PEER REVIEW 11 of 15 effect was used to measure particle temperature. Figure 11 shows the temperatures of 13 particles trapped using laser light at 1031 nm and 1064 nm at a pressure of 600 mbar The trapping power was 200 mW. Some changes between the temperatures of particles trapped by 1031 nm light was observed. The average internal temperature was ~ 167 K.
The lowest temperature reached in this experiment was ~ 130 K. As one can see in Figure  11, in the beginning of the experiment, particle 2 was trapped using a 1031 nm laser beam. The temperature reached was 175 K. Later, the laser wavelength was changed to 1064 nm. For this wavelength, the temperature reached was 255 K (Figure 11). For this particle, the temperature was 175 K when trapped with 1031 nm light and 255 K when trapped using 1064 nm light ( Figure 11). Using (6) and (7), one can show that lower temperatures can be achieved with a laser at 1020 nm. Lower pressure and improved crystal purity can also help to reach lower temperatures. Anti-Stokes cooling of Yb 3+ :NaYF nanocrystals was reported in [39]. A sublevel splitting of Yb 3+ ions doped in a NaYF host is smaller than that of Yb 3+ ions doped in other crystal hosts. Yb 3+ :NaYF can provide better cooling efficiency. Yb 3+ :NaYF particles appeared as spheres (insert in Figure 12). The fabrication process of Yb 3+ :NaYF nanoparticles is safer than the fabrication process of Yb 3+ :YLF nanoparticles. It is predicted that Yb 3+ :NaYF can be cooled to ~ 5 K colder than Yb 3+ :YLF crystals [38]. In 2021, 10% Yb 3+ :NaLF4 optically levitated nanoparticles were laser-cooled by 42 K to a temperature of 252 K [39]. The ytterbium ions were excited with a laser at 1020 nm.  Anti-Stokes cooling of Yb 3+ :NaYF nanocrystals was reported in [39]. A sublevel splitting of Yb 3+ ions doped in a NaYF host is smaller than that of Yb 3+ ions doped in other crystal hosts. Yb 3+ :NaYF can provide better cooling efficiency. Yb 3+ :NaYF particles appeared as spheres (insert in Figure 12). The fabrication process of Yb 3+ :NaYF nanoparticles is safer than the fabrication process of Yb 3+ :YLF nanoparticles. It is predicted that Yb 3+ :NaYF can be cooled to ∼ 5 K colder than Yb 3+ :YLF crystals [38]. In 2021, 10% Yb 3+ :NaLF 4 optically levitated nanoparticles were laser-cooled by 42 K to a temperature of 252 K [39]. The ytterbium ions were excited with a laser at 1020 nm.
Appl. Sci. 2022, 11, x FOR PEER REVIEW 11 of 15 effect was used to measure particle temperature. Figure 11 shows the temperatures of 13 particles trapped using laser light at 1031 nm and 1064 nm at a pressure of 600 mbar The trapping power was 200 mW. Some changes between the temperatures of particles trapped by 1031 nm light was observed. The average internal temperature was ~ 167 K.
The lowest temperature reached in this experiment was ~ 130 K. As one can see in Figure  11, in the beginning of the experiment, particle 2 was trapped using a 1031 nm laser beam. The temperature reached was 175 K. Later, the laser wavelength was changed to 1064 nm. For this wavelength, the temperature reached was 255 K (Figure 11). For this particle, the temperature was 175 K when trapped with 1031 nm light and 255 K when trapped using 1064 nm light ( Figure 11). Using (6) and (7), one can show that lower temperatures can be achieved with a laser at 1020 nm. Lower pressure and improved crystal purity can also help to reach lower temperatures. Anti-Stokes cooling of Yb 3+ :NaYF nanocrystals was reported in [39]. A sublevel splitting of Yb 3+ ions doped in a NaYF host is smaller than that of Yb 3+ ions doped in other crystal hosts. Yb 3+ :NaYF can provide better cooling efficiency. Yb 3+ :NaYF particles appeared as spheres (insert in Figure 12). The fabrication process of Yb 3+ :NaYF nanoparticles is safer than the fabrication process of Yb 3+ :YLF nanoparticles. It is predicted that Yb 3+ :NaYF can be cooled to ~ 5 K colder than Yb 3+ :YLF crystals [38]. In 2021, 10% Yb 3+ :NaLF4 optically levitated nanoparticles were laser-cooled by 42 K to a temperature of 252 K [39]. The ytterbium ions were excited with a laser at 1020 nm. A 1064 nm laser was used to create an optical dipole trap. A 1064 nm laser induced a dipole in the Yb 3+ :NaYF nanocrystal. The induced dipole moved to the point of highest intensity of the laser field. The beam of a 1020 nm laser was collinear to the beam of the 1064 nm laser. The 1020 nm laser was responsible for anti-Stokes cooling in the system. A 1064 nm laser was used to create an optical dipole trap. A 1064 nm laser induced a dipole in the Yb 3+ :NaYF nanocrystal. The induced dipole moved to the point of highest intensity of the laser field. The beam of a 1020 nm laser was collinear to the beam of the 1064 nm laser. The 1020 nm laser was responsible for anti-Stokes cooling in the system.
As one can see in Figure 12, anti-Stokes photoluminescence (PL) of the trapped crystal was collected with the lens used to focus the trapping and cooling laser beams. A 1000 nm dichroic beam splitter and a 1000 nm short-pass filter permitted the separation of the PL from the laser light. The filter cut-off PL at wavelengths higher than 1000 nm ( Figure 12). The collected PL was analyzed with an OSA. These spectra were used to estimate the temperature of the sample. The lowest temperatures reached by seven nanocrystals are presented in Figure 13.
Appl. Sci. 2022, 11, x FOR PEER REVIEW 12 of 15 As one can see in Figure 12, anti-Stokes photoluminescence (PL) of the trapped crystal was collected with the lens used to focus the trapping and cooling laser beams. A 1000 nm dichroic beam splitter and a 1000 nm short-pass filter permitted the separation of the PL from the laser light. The filter cut-off PL at wavelengths higher than 1000 nm ( Figure 12). The collected PL was analyzed with an OSA. These spectra were used to estimate the temperature of the sample. The lowest temperatures reached by seven nanocrystals are presented in Figure 13.  The scheme in Figure 14 consists of a Paul trap for nanocrystal spatial localization [40,44] and a 3D standing light wave for the optical cooling of nanocrystals [45]. The standing light wave can be replaced by a single light wave. Indeed, as is evident from the foregoing, RE-doped nanocrystals can be optically cooled with a single optical wave. In the theoretical analysis of optical cooling presented in Figure 14, 5% Yb 3+ :CsF2 nanocrystals were considered. It was assumed that the nanocrystal of mass = 6.4 × 10 −19 g was placed in the center of the Paul trap ( 0 = 0 = 0 = 0) at initial velocities of 0 = 0 = 0 = 10 cm ⋅ s −1 . The nanocrystal dynamics consist of the motions of two kinds: the secular motion and the small oscillations about the secular path ( Figure 14). The secular motion takes place along the smooth trajectory, which converge to the center of the Paul trap. The interaction between the charged nanocrystals and radio-frequency field results in small oscillations about the secular path. The numerical simulations presented in [42] prove that the Yb 3+ :CsF2 nanocrystal mechanical action can reach a value ~10ℏ, which is close to the quantum limit of ℏ.

Conclusions
Optical refrigeration based on anti-Stokes fluorescence has dramatically improved over last several years. Early in the development of laser cooling of RE-doped low-phonon solids, -and −sized samples were the focus of attention. These bulk samples are interesting for the development of an all-optical solid-state cryocooler, which is free from mechanical vibrations, moving parts, and liquids. The first all-solid-state optical cryocooler started its operation in 2018 [46]. Laser cooling with anti-Stokes fluorescence in bulk samples can also be used for heat mitigation in high-power lasers [47].
At the moment, optical refrigeration of solids has moved to the area of small and -sized RE-doped low-phonon particles. A hybrid cooling mechanism resulting in a decrease in translational temperature and in internal temperature of the sample has been considered in this paper. The theoretical background of optical refrigeration based on anti-Stokes fluorescence has been presented and discussed. It is shown that Mie resonances can enhance laser cooling in a number of properly designed samples. It has been shown that hybrid trapping can be used for the localization and internal cooling of RE-doped nanocrystals. This approach can be applied to the fundamental studies in biophysics, chemistry, and medicine.
It has been shown that hybrid laser cooling can be realized using a single-beam optical trapping. This approach can be useful for physiological laser refrigeration and for fu- The scheme in Figure 14 consists of a Paul trap for nanocrystal spatial localization [40,44] and a 3D standing light wave for the optical cooling of nanocrystals [45]. The standing light wave can be replaced by a single light wave. Indeed, as is evident from the foregoing, RE-doped nanocrystals can be optically cooled with a single optical wave. In the theoretical analysis of optical cooling presented in Figure 14, 5% Yb 3+ :CsF 2 nanocrystals were considered. It was assumed that the nanocrystal of mass M = 6.4 × 10 −19 g was placed in the center of the Paul trap (x 0 = y 0 = z 0 = 0) at initial velocities of v 0x = v 0y = v 0z = 10 cm · s −1 . The nanocrystal dynamics consist of the motions of two kinds: the secular motion and the small oscillations about the secular path ( Figure 14). The secular motion takes place along the smooth trajectory, which converge to the center of the Paul trap. The interaction between the charged nanocrystals and radio-frequency field results in small oscillations about the secular path. The numerical simulations presented in [42] prove that the Yb 3+ :CsF 2 nanocrystal mechanical action can reach a value ∼ 10 , which is close to the quantum limit of .

Conclusions
Optical refrigeration based on anti-Stokes fluorescence has dramatically improved over last several years. Early in the development of laser cooling of RE-doped low-phonon solids, mmand cm− sized samples were the focus of attention. These bulk samples are interesting for the development of an all-optical solid-state cryocooler, which is free from mechanical vibrations, moving parts, and liquids. The first all-solid-state optical cryocooler started its operation in 2018 [46]. Laser cooling with anti-Stokes fluorescence in bulk samples can also be used for heat mitigation in high-power lasers [47].
At the moment, optical refrigeration of solids has moved to the area of small µmand nm-sized RE-doped low-phonon particles. A hybrid cooling mechanism resulting in a decrease in translational temperature and in internal temperature of the sample has been considered in this paper. The theoretical background of optical refrigeration based on anti-Stokes fluorescence has been presented and discussed. It is shown that Mie resonances can enhance laser cooling in a number of properly designed samples. It has been shown that hybrid trapping can be used for the localization and internal cooling of RE-doped nanocrystals. This approach can be applied to the fundamental studies in biophysics, chemistry, and medicine.
It has been shown that hybrid laser cooling can be realized using a single-beam optical trapping. This approach can be useful for physiological laser refrigeration and for future levitated hybrid opto-mechanics experiments, where the control of translational and internal temperatures is important. For example, using this approach and Yb 3+ :YLF nanocrystals, one can build a levitated nanocryostat. The birefringence of the particles can be used to control their space orientation. Laser cooling of RE-doped nanocrystals can also be used to mitigate heat generated in spherical microlasers [48].
Laser hybrid cooling of micro-and nanoparticles is in its infancy. I believe that further miniaturization of the RE-doped samples that undergo hybrid laser cooling is in sight. In very small samples (several nm), quantization of phonon modes is very important. As the size of the nanocrystal decreases, its phonon cut-off frequency increases. This effect can deteriorate thermalization in very small particles [27]. Engineering of the host material can be a solution to this problem.