Phase transitions in femtosecond laser ablation

In this study we simulate an interaction of femtosecond laser pulses (100 fs, 800 nm, 0.1-10 J/cm²) with metal targets of Al, Au, Cu, and Ni. For analysis of laser-induced phase transitions, melting and shock waves propagation as well as material decomposition we use an Eulerian hydrocode in conjunction with a thermodynamically complete two-temperature equation of state with stable and metastable phases. Isochoric heating, material evaporation from the free surface of the target and fast propagation of the melting and shock waves are observed. On rarefaction the liquid phase becomes metastable and its lifetime is estimated using the theory of homogeneous nucleation. Mechanical spallation of the target material at high strain rates is also possible as a result of void growth and confluence. In our simulation several ablation mechanisms are taken into account but the main issue of the material is found to originate from the metastable liquid state. It can be decomposed either into a liquid-gas mixture in the vicinity of the critical point, or into droplets at high strain rates and negative pressure. The simulation results are in agreement with available experimental findings.     

Effects of laser operating parameters on metals micromachining with ultrafast lasers

180 femtoseconds (1 kHz) and 10 picoseconds (1-50 kHz) ultrafast laser micro-structuring of the metals Ti alloy, Al and Cu have been studied for the purpose of industrial application. The effects of some key laser operating parameters were investigated. The evolution of surface morphology revealed that laser pulses overlap in a range around the spatial FWHM can help to achieve optimal residual surface roughness. While observed ablation rate (unit: μm³ per pulse) changed dramatically with repetition rate due to the combined effects of plasma absorption, residual thermal energy and phase transition, higher throughput can be achieved with higher repetition rate. This study also indicated that residual surface roughness is almost independent of repetition rate at 10 ps temporal pulse length. The ablation depth is approximately proportional to the number of overscan; however, machining accuracy deteriorates, especially for femtosecond laser processing and metals with low thermal conductivity and short electron-phonon coupling time.     

Electron-phonon and vibronic couplings in the FMO bacteriochlorophyll a antenna complex studied by difference fluorescence line narrowing

The electron-phonon and vibronic couplings governing the spectral properties have been studied in the Fenna-Matthews-Olson (FMO) bacteriochlorophyll a (BChl a)-protein complex at 4.5 K using a spectrally selective difference fluorescence line-narrowing technique. The complex is a part of the light-harvesting system of green photosynthetic bacteria. Its lowest-energy absorption band, peaking at 826 nm and responsible for the fluorescence, is believed to be due to Q_y transitions of largely isolated molecules. One of the main merits of the used method compared with the more common fluorescence line narrowing is that the zero-phonon lines (ZPL) resonant with the excitation laser can be accurately measured, allowing precise determination of the Huang-Rhys (HR) factors, the main characteristics of the linear electron-phonon and vibronic coupling strengths. Over 60 individual vibrational modes of intra- and intermolecular origin have been resolved in the energy range of 45-1600 cm⁻¹. The small HR factors for these modes, ranging between 0.001 and 0.018, add up to a value of S_vib=0.38±0.07. The effective HR factor for the phonons, S_ph, was found clearly wavelength-dependent, varying from ∼0.7 at short wavelengths to ∼0.3 at the long-wavelength tail of the absorption spectrum. Coupling between the BChl a molecules is likely responsible for this wavelength dependence.     

Modeling of fast phase transitions dynamics in metal target irradiated by pico- and femtosecond pulsed laser

We investigate laser pulse influence on aluminum target in irradiance range 10⁹ to 10¹⁶ W/cm², pulse duration between 10⁻⁸ and 10⁻¹⁵ s, Gaussian time profile with wavelength of 0.8 μm. For all computations energy density was 10 J/cm². Plasma in the evaporated material is generated at the energy density above 10 J/cm²as the modeling showed.Long and short laser pulses distinguish by the mechanisms of energy transformation. For short laser pulses there is volumetric energy absorption, together with rapid phase transitions it lead to overheating in solid and liquid states, overheated solid temperature rises up to (6-8)T_m. Under influence of the energy saved in overheated solid, duration of the phase transitions becomes nanosecond, which is several orders of magnitude longer than laser pulse.     

Temperature dependences of the electron-phonon coupling, electron heat capacity and thermal conductivity in Ni under femtosecond laser irradiation

The electron temperature dependences of the electron-phonon coupling factor, electron heat capacity and thermal conductivity are investigated for Ni in a range of temperatures typically realized in femtosecond laser material processing applications, from room temperature up to temperatures of the order of 10⁴ K. The analysis is based on the electronic density of states obtained through the electronic structure calculations. Thermal excitation of d band electrons is found to result in a significant decrease in the strength of the electron-phonon coupling, as well as large deviations of the electron heat capacity and the electron thermal conductivity from the commonly used linear temperature dependences on the electron temperature. Results of the simulations performed with the two-temperature model demonstrate that the temperature dependence of the thermophysical parameters accounting for the thermal excitation of d band electrons leads to higher maximum lattice and electron temperatures achieved at the surface of an irradiated Ni target and brings the threshold fluences for surface melting closer to the experimentally measured values as compared to the predictions obtained with commonly used approximations of the thermophysical parameters.     

Ablation process of silica glass induced by laser plasma soft X-ray irradiation

Silica glass can be machined by irradiation with laser plasma soft X-rays on nano- and micrometer scale. We have investigated the ablation process of silica glass induced by laser plasma soft X-ray irradiation. We observed ionic and neutral species emitted from silica surfaces after irradiation. Dominant ions and neutrals are O⁺ and Si⁺ ions and Si, O, SiO and Si₂ neutrals, respectively. The ions have kinetic energies of 13 and 25 eV, which are much higher than those of particles emitted by evaporation. The energy of laser plasma soft X-rays absorbed to silica glass at a fluence of 1.4 J/cm² is estimated to be 380 kJ/cm³, which is higher than the binding energy of SiO₂ of 76 kJ/cm³. These results suggest that the most of the bonds in silica glass are broken by absorption of laser plasma soft X-rays, that several percent of the atoms are ionized, and that neutral atoms are emitted together with repulsive ions. The process possibly enables us to fabricate nano structures.     

Substrate influence in electron-phonon coupling measurements in thin Au films

Accurate understanding and measurement of the energy transfer mechanisms during thermal nonequilibrium between electrons and the surrounding material systems is critical for a wide array of applications. With device dimensions decreasing to sizes on the order of the thermal penetration depth, the equilibration of the electrons could be effected by boundary effects in addition to electron-phonon coupling. In this study, the rate of electron equilibration in 20 nm thick Au films is measured with the Transient ThermoReflectance (TTR) technique. At very large incident laser fluences which result in very high electron temperatures, the electron-phonon coupling factors determined from TTR measurements deduced using traditional models are almost an order of magnitude greater than predicted from theory. By taking excess electron energy loss via electron-substrate transport into account with a proposed three temperature model, TTR electron-phonon coupling factor measurements are more in line with theory, indicating that in highly nonequilibrium situations, the high temperature electron system looses substantial energy to the substrate in addition to that transferred to the film lattice through coupling.     

Two-temperature relaxation and melting after absorption of femtosecond laser pulse

The theory and experiments concerned with the electron-ion thermal relaxation and melting of overheated crystal lattice constitute the subject of this paper. The physical model includes two-temperature (2T) equation of state, many-body interatomic potential, the electron-ion energy exchange, electron thermal conductivity, and optical properties of solid, liquid, and two phase solid-liquid mixture. Two-temperature hydrodynamics and molecular dynamics codes are used. An experimental setup with pump-probe technique is used to follow evolution of an irradiated target with a short time step 100 fs between the probe femtosecond laser pulses. Accuracy of measurements of reflection coefficient and phase of reflected probe light are 1% and ∼1 nm, respectively. It is found that, firstly, the electron-electron collisions make a minor contribution to a light absorption in solid Al at moderate intensities; secondly, the phase shift of a reflected probe results from heating of ion subsystem and kinetics of melting of Al crystal during  ps, where t is time delay between the pump and probe pulses measured from the maximum of the pump; thirdly, the optical response of Au to a pump shows a marked contrast to that of Al on account of excitation of d-electrons.     

Femtosecond laser ablation of aluminium

We investigate femtosecond laser ablation of aluminium using a hybrid simulation scheme. Two equations are solved simultaneously: one for the electronic system, which accounts for laser energy absorption and heat conduction, the other for the dynamics of the lattice where the ablation process takes place. For the electron-temperature a generalized heat-conduction equation is solved by applying a finite difference scheme. For the lattice properties, e.g. pressure, density or temperature, we use common molecular dynamics. Energy transfer between the subsystems is allowed by introducing an electron-phonon coupling term. This combined treatment of the electronic and atomic systems is an extension of the well known two-temperature model [Anisimov, Kapeliovich, Perel’man, Electron emission from metal surfaces exposed to ultra short laser pulses, JETP Lett. 39 (2)].     

The influence of laser-induced surface modifications on the backside etching process

Spectroscopic measurements in the UV/VIS region show reduced transmission through laser-induced backside wet etching (LIBWE) of fused silica. Absorption coefficients of up to 10⁵ cm⁻¹ were calculated from the transmission measurements for a solid surface layer of about 50 nm. The temperatures near the interface caused by laser pulse absorption, which were analytically calculated using a new thermal model considering interface and liquid volume absorption, can reach 10⁴ K at typical laser fluences. The high absorption coefficients and the extreme temperatures give evidence for an ablation-like process that is involved in the LIBWE process causing the etching of the modified near-surface region. The confinement of the ablation/etching process to the modified near-surface material region can account for the low etch rates observed in comparison to front-side ablation.     

Empirical low energy values for methane transitions in the 5852-6181 cm region by absorption spectroscopy at 81 K

The high resolution absorption spectrum of methane has been recorded at liquid nitrogen temperature by direct absorption spectroscopy between 1.62 and 1.71 μm (5852-6181 cm⁻¹) using a newly developed cryogenic cell and a series of distributed feedback (DFB) laser diodes. The minimum value of the measured line intensities is on the order of 3 × 10⁻²⁶ cm/molecule The investigated spectral range corresponds to the high energy part of the tetradecad dominated by the 2ν₃ band for which a theoretical treatment is not yet available. The positions and strengths at 81 K of 2187 transitions were obtained from the spectrum analysis. From the values of the line strength at liquid nitrogen and room temperatures, the low energy values of 845 transitions could be determined. The obtained results are discussed in relation with the previous work of Margolis and compared to the line list provided by the HITRAN database.     

Change of exciton-phonon interaction in layered ferroelastic crystals (Cs3Bi2I9)

The half-width of exciton absorption band (n=1) of Cs₃Bi₂I₉ layered ferroelastic crystals was studied carefully as function of temperature in the range from 5 to 300 K. For the first time, we have found a new physical effect: change of exciton-phonon interaction (from weak to strong) in the same sample as temperature increases. It was established that the temperature value T^*=150 K may be considered as characteristic one, below which a crystal loses the nature of layered substance. The effect is explained using a model that takes into account the reconstruction of the crystal lattice from non-layered to layered one.     

Simulation of dynamics of phase transitions in CdTe by pulsed laser irradiation

Numerical simulation of melting and solidification processes induced in CdTe by nanosecond radiation of ruby laser (λ = 694 nm, τ = 20 and 80 ns) and KrF excimer laser (λ = 248 nm, τ = 20 ns) taking into account components diffusion in melt and their evaporation from the surface has been carried out. Cd atoms evaporation has shown to essentially affect the dynamics of phase transitions in the near-surface region. Thus, in the case of the influence of ruby laser irradiation intensive surface cooling results in the formation of nonmonotone temperature profile with maximum temperature in semiconductor volume at the distance of ∼20 nm from the surface. The melt formed under the surface extends both to the surface and to the semiconductor volume as well. As a result of cadmium telluride components evaporation and diffusion in the melt the near-surface region is enriched with tellurium. The obtained melting threshold value of irradiation energy density is in a reasonable agreement with experimental data.     

Pressure dependence study on the electron-phonon coupling of thulium hexachloroelpasolite

Raman spectra of Cs₂NaTmCl₆ have been recorded using a diamond anvil cell at ambient temperature. The vibrational energy of each of the Raman-active TmCl₆⁻³ moiety modes increases linearly with pressure. The integrated band areas of the ν₁(a_1g) and ν₂(e_g) modes are independent of applied pressure. However, the band area of the ν₅(t_2g) mode shows an anomalous behaviour, which has been qualitatively interpreted as due to electron-phonon coupling of the aΓ₅ electronic state with the Γ₁+ν₅(t_2g) vibronic state. This interaction between the coupled states is strongest between ca. 10 and 13 GPa at ambient temperature. The results serve to emphasize the specificity of the occurrence of strong electron-phonon coupling for particular transitions of a given rare earth ion.     

Successive phase transitions of [(PyO)(H/D)][AuCl4] (PyO=C5H5NO)

³⁵Cl NQR as well as heat capacity measurements of [(PyO)H][AuCl₄] and its deuterated analog [(PyO)D][AuCl₄] revealed successive phase transitions at 70.5 and 62.5 K, and at 71 and 63 K, respectively. The NQR frequency varied continuously through the upper transition point while discontinuously through the lower transition point. In the intermediate-temperature phase a remarkable decrease in the signal intensity was observed. These NQR observations as well as the feature of the heat capacity anomaly in which a broad peak is succeeded by a sharp peak with decreasing temperature suggest a possibility of normal-incommensurate-commensurate phase sequence.     

Wide bandwidth phase-locked diode laser with an intra-cavity electro-optic modulator

Two extended cavity laser diodes are phase-locked, thanks to an intra-cavity electro-optical modulator. The phase-locked loop bandwidth is on the order of 10 MHz, which is about twice larger than when the feedback correction is applied on the laser current. The phase noise reaches −120 dBrad²/Hz at 10 kHz. This new scheme reduces the residual laser phase noise, which constitutes one of the dominant contributions in the sensitivity limit of atom interferometers using two-photon transitions.     

Generation of nanospikes via laser ablation of metals in liquid environment and their activity in surface-enhanced Raman scattering of organic molecules

The formation of dense arrays of nanospikes occurs under laser ablation of bulk targets (Ag, Au, Ta, Ti) immersed in liquids such as water or ethanol. The average height of spikes is 50 nm and their density on the target amounts to 10¹⁰ cm⁻². The effect is observed with sufficiently short laser pulses. In particular, either a 350 ps or a 90 ps Nd:YAG lasers are used in their fundamental harmonics. The nanospikes are characterized by UV-Visible reflection spectrometry and atomic force microscopy. The oscillations of electrons within nanospikes result in a permanent coloration of the surface and a modification of the optical reflection spectra of the metal. Scanning the laser beam along the metal surface allows its nanostructuring over extended areas (∼1 cm²). The nanostructured Ag surface shows enhanced Raman scattering of acridine molecules at a concentration of 10⁻⁵ M/l, whereas the initial Ag targets do not show any signal within the accuracy of measurements.     

VUV spectroscopy of pure and Tm-doped LicaAlF6 crystals

Nominally pure and Tm³⁺-doped LiCaAlF₆ crystals were grown by the Czochralski technique in a reducing atmosphere. The optical properties of transparent single crystals were studied using absorption and time-resolved luminescence spectroscopy in the VUV spectral range (330-100 nm). The strongest VUV emission peaking at 60 800 cm⁻¹ with a decay time of 5.6 μs (7 μs) at 300 K (7.4 K) was assigned to the spin-forbidden 4f¹¹5d-4f¹² transition of Tm³⁺. The fine structure observed in the VUV emission and corresponding excitation spectra indicate intermediate strength of electron-phonon coupling in this system. The efficient excitation of f-f emissions above 72 000 cm⁻¹, higher than the onset of f-d absorption at 63 000 cm⁻¹, is mainly caused by the F⁻ to Tm³⁺ charge transfer absorption. The nature of various host-related excitation processes in the energy transfer to the Tm³⁺ ions is discussed.     

Magnetic-field-driven quantum critical behavior in graphite and bismuth

We study magnetotransport properties of graphite and rhombohedral bismuth samples and found that in both materials applied magnetic field induces the metal-insulator- (MIT) and reentrant insulator-metal-type (IMT) transformations. The corresponding transition boundaries plotted on the magnetic field-temperature (B − T) plane nearly coincide for these semimetals and can be best described by power laws T ∼ (B − B_c)^κ, where B_c is a critical field at T = 0 and κ = 0.45 ± 0.05. We show that insulator-metal-insulator (I-M-I) transformations take place in the Landau level quantization regime and illustrate how the IMT in quasi-3D graphite transforms into a cascade of I-M-I transitions, related to the quantum Hall effect in quasi-2D graphite samples. We discuss the possible coupling of superconducting and excitonic correlations with the observed phenomena, as well as signatures of quantum phase transitions associated with the M-I and I-M transformations.     

Pressure-induced structural phase transition in ReO3

We have investigated the pressure-induced structural phase transition in ReO₃ by neutron diffraction on a single crystal. We collected neutron diffraction intensities from the ambient and high pressure phases at P=7 kbar and refined the crystal structures. We have determined the stability of the high pressure phase as a function temperature down to T=2 K and have constructed the (P-T) phase diagram. The critical pressure is P_c=5.2 kbar at T=300 K and decreases almost linearly with decreasing temperature to become P_c=2.5 kbar at T=50 K. The phase transition is driven by the softening of the M₃ phonon mode. The high pressure phase is formed by the rigid rotation of almost undistorted ReO₆ octahedra and the Re-O-Re angle deviates from 180°. We do not see any evidence for the existence of the tetragonal (P4/mbm) intermediate pressure phase reported earlier.