Atomistic-continuum modeling of ultrafast laser-induced melting of silicon targets

In this work, we present an atomistic-continuum model for simulations of ultrafast laser-induced melting processes in semiconductors on the example of silicon. The kinetics of transient non-equilibrium phase transition mechanisms is addressed with MD method on the atomic level, whereas the laser light absorption, strong generated electron-phonon nonequilibrium, fast heat conduction, and photo-excited free carrier diffusion are accounted for with a continuum TTM-like model (called nTTM). First, we independently consider the applications of nTTM and MD for the description of silicon, and then construct the combined MD-nTTM model. Its development and thorough testing is followed by a comprehensive computational study of fast nonequilibrium processes induced in silicon by an ultrashort laser irradiation. The new model allowed to investigate the effect of laser-induced pressure and temperature of the lattice on the melting kinetics. Two competing melting mechanisms, heterogeneous and homogeneous, were identified in our big-scale simulations. Apart from the classical heterogeneous melting mechanism, the nucleation of the liquid phase homogeneously inside the material significantly contributes to the melting process. The simulations showed, that due to the open diamond structure of the crystal, the laser-generated internal compressive stresses reduce the crystal stability against the homogeneous melting. Consequently, the latter can take a massive character within several picoseconds upon the laser heating. Due to the large negative volume of melting of silicon, the material contracts upon the phase transition, relaxes the compressive stresses, and the subsequent melting proceeds heterogeneously until the excess of thermal energy is consumed. A series of simulations for a range of absorbed fluences allowed us to find the threshold fluence value at which homogeneous liquid nucleation starts contributing to the classical heterogeneous propagation of the solid-liquid interface. A series of simulations for a range of the material thicknesses showed that the sample width we chosen in our simulations (800 nm) corresponds to a thick sample. Additionally, in order to support the main conclusions, the results were verified for a different interatomic potential. Possible improvements of the model to account for nonthermal effects are discussed and certain restrictions on the suitable interatomic potentials are found. As a first step towards the inclusion of these effects into MD-nTTM, we performed nanometer-scale MD simulations with a new interatomic potential, designed to reproduce ab initio calculations at the laser-induced electronic temperature of 18946 K. The simulations demonstrated that, similarly to thermal melting, nonthermal phase transition occurs through nucleation. A series of simulations showed that higher (lower) initial pressure reinforces (hinders) the creation and the growth of nonthermal liquid nuclei. For the example of Si, the laser melting kinetics of semiconductors was found to be noticeably different from that of metals with a face-centered cubic crystal structure. The results of this study, therefore, have important implications for interpretation of experimental data on the kinetics of melting process of semiconductors.

Femtosecond imaging-mode laser-induced breakdown spectroscopy

Many nonlinear optical microscopy techniques based on the high-intensity nonlinear phenomena were developed recent years. A new technique based on the minimal-invasive in-situ analysis of the specific bound elements in biological samples is described in the present work. The imaging-mode Laser-Induced Breakdown Spectroscopy (LIBS) is proposed as a combination of LIBS, femtosecond laser material processing and microscopy. The Calcium distribution in the peripheral cell wall of the sunflower seedling (Helianthus Annuus L.) stem is studied as a first application of the imaging-mode LIBS. At first, several nonlinear optical microscopy techniques are overviewed. The spatial resolution of the imaging-mode LIBS microscope is discussed basing on the Point-Spread Function (PSF) concept. The primary processes of the Laser-Induced Breakdown (LIB) are overviewed. We consider ionization, breakdown, plasma formation and ablation processes. Water with defined Calcium salt concentration is used as a model of the biological object in the preliminary experiments. The transient LIB spectra are measured and analysed for both nanosecond and femtosecond laser excitation. The experiment on the local Calcium concentration measurements in the peripheral cell wall of the sunflower seedling stem employing nanosecond LIBS shows, that nanosecond laser is not a suitable excitation source for the biological applications. In case of the nanosecond laser the ablation craters have random shape and depth over 20 µm. The analysis of the femtosecond laser ablation craters shows the reproducible circle form. At 3.5 µJ laser pulse energy the diameter of the crater is 4 µm and depth 140 nm for single laser pulse, which results in 1 femtoliter analytical volume. The experimental result of the 2 dimensional and surface sectioning of the bound Calcium concentrations is presented in the work.

Theory of laser-induced ultrafast structural changes in solids

The present thesis is a contribution to the study of laser-solid interaction. Despite the numerous applications resulting from the recent use of laser technology, there is still a lack of satisfactory answers to theoretical questions regarding the mechanism leading to the structural changes induced by femtosecond lasers in materials. We provide here theoretical approaches for the description of the structural response of different solids (cerium, samarium sulfide, bismuth and germanium) to femtosecond laser excitation. Particular interest is given to the description of the effects of the laser pulse on the electronic systems and changes of the potential energy surface for the ions. Although the general approach of laser-excited solids remains the same, the potential energy surface which drives the structural changes is calculated with different theoretical models for each material. This is due to the difference of the electronic properties of the studied systems. We use the Falicov model combined with an hydrodynamic method to study photoinduced phase changes in cerium. The local density approximation (LDA) together with the Hubbard-type Hamiltonian (LDA+U) in the framework of density functional theory (DFT) is used to describe the structural properties of samarium sulfide. We parametrize the time-dependent potential energy surface (calculated using DFT+ LDA) of bismuth on which we perform quantum dynamical simulations to study the experimentally observed amplitude collapse and revival of coherent $A_{1g}$ phonons. On the basis of a time-dependent potential energy surface calculated from a non-orthogonal tight binding Hamiltonian, we perform molecular dynamics simulation to analyze the time evolution (coherent phonons, ultrafast nonthermal melting) of germanium under laser excitation. The thermodynamic equilibrium properties of germanium are also reported. With the obtained results we are able to give many clarifications and interpretations of experimental results and also make predictions.

Many-body theory of laser-induced ultrafast demagnetization and angular momentum transfer in ferromagnetic transition metals

An electronic theory is developed, which describes the ultrafast demagnetization in itinerant ferromagnets following the absorption of a femtosecond laser pulse. The present work intends to elucidate the microscopic physics of this ultrafast phenomenon by identifying its fundamental mechanisms. In particular, it aims to reveal the nature of the involved spin excitations and angular-momentum transfer between spin and lattice, which are still subjects of intensive debate. In the first preliminary part of the thesis the initial stage of the laser-induced demagnetization process is considered. In this stage the electronic system is highly excited by spin-conserving elementary excitations involved in the laser-pulse absorption, while the spin or magnon degrees of freedom remain very weakly excited. The role of electron-hole excitations on the stability of the magnetic order of one- and two-dimensional 3d transition metals (TMs) is investigated by using ab initio density-functional theory. The results show that the local magnetic moments are remarkably stable even at very high levels of local energy density and, therefore, indicate that these moments preserve their identity throughout the entire demagnetization process. In the second main part of the thesis a many-body theory is proposed, which takes into account these local magnetic moments and the local character of the involved spin excitations such as spin fluctuations from the very beginning. In this approach the relevant valence 3d and 4p electrons are described in terms of a multiband model Hamiltonian which includes Coulomb interactions, interatomic hybridizations, spin-orbit interactions, as well as the coupling to the time-dependent laser field on the same footing. An exact numerical time evolution is performed for small ferromagnetic TM clusters. The dynamical simulations show that after ultra-short laser pulse absorption the magnetization of these clusters decreases on a time scale of hundred femtoseconds. In particular, the results reproduce the experimentally observed laser-induced demagnetization in ferromagnets and demonstrate that this effect can be explained in terms of the following purely electronic non-adiabatic mechanism: First, on a time scale of 10–100 fs after laser excitation the spin-orbit coupling yields local angular-momentum transfer between the spins and the electron orbits, while subsequently the orbital angular momentum is very rapidly quenched in the lattice on the time scale of one femtosecond due to interatomic electron hoppings. In combination, these two processes result in a demagnetization within hundred or a few hundred femtoseconds after laser-pulse absorption.

Optimierung der laserinduzierten Plasmaspektroskopie an Metallen durch Femtosekunden-Doppelpulse für die Entwicklung eines hochauflösenden 3D-Rasterabbildungsverfahrens

Die laserinduzierte Plasmaspektroskopie (LIPS) ist eine spektrochemische Elementanalyse zur Bestimmung der atomaren Zusammensetzung einer beliebigen Probe. Für die Analyse ist keine spezielle Probenpräparation nötig und kann unter atmosphärischen Bedingungen an Proben in jedem Aggregatzustand durchgeführt werden. Femtosekunden Laserpulse bieten die Vorteile einer präzisen Ablation mit geringem thermischen Schaden sowie einer hohen Reproduzierbarkeit. Damit ist fs-LIPS ein vielversprechendes Werkzeug für die Mikroanalyse technischer Proben, insbesondere zur Untersuchung ihres Ermüdungsverhaltens. Dabei ist interessant, wie sich die initiierten Mikrorisse innerhalb der materialspezifschen Struktur ausbreiten. In der vorliegenden Arbeit sollte daher ein schnelles und einfach zu handhabendes 3D-Rasterabbildungsverfahren zur Untersuchung der Rissausbreitung in TiAl, einer neuen Legierungsklasse, entwickelt werden. Dazu wurde fs-LIPS (30 fs, 785 nm) mit einem modifizierten Mikroskopaufbau (Objektiv: 50x/NA 0.5) kombiniert, welcher eine präzise, automatisierte Probenpositionierung ermöglicht. Spektrochemische Sensitivität und räumliches Auflösungsvermögen wurden in energieabhängigen Einzel- und Multipulsexperimenten untersucht. 10 Laserpulse pro Position mit einer Pulsenergie von je 100 nJ führten in TiAl zum bestmöglichen Kompromiss aus hohem S/N-Verhältnis von 10:1 und kleinen Lochstrukturen mit inneren Durchmessern von 1.4 µm. Die für das Verfahren entscheidende laterale Auflösung, dem minimalen Lochabstand bei konstantem LIPS-Signal, beträgt mit den obigen Parametern 2 µm und ist die bislang höchste bekannte Auflösung einer auf fs-LIPS basierenden Mikro-/Mapping-Analyse im Fernfeld. Fs-LIPS Scans von Teststrukturen sowie Mikrorissen in TiAl demonstrieren eine spektrochemische Sensitivität von 3 %. Scans in Tiefenrichtung erzielen mit denselben Parametern eine axiale Auflösung von 1 µm. Um die spektrochemische Sensitivität von fs-LIPS zu erhöhen und ein besseres Verständnis für die physikalischen Prozesse während der Laserablation zu erhalten, wurde in Pump-Probe-Experimenten untersucht, in wieweit fs-Doppelpulse den laserinduzierten Abtrag sowie die Plasmaemission beeinflussen. Dazu wurden in einem Mach-Zehnder-Interferometer Pulsabstände von 100 fs bis 2 ns realisiert, Gesamtenergie und Intensitätsverhältnis beider Pulse variiert sowie der Einfluss der Materialparameter untersucht. Sowohl das LIPS-Signal als auch die Lochstrukturen zeigen eine Abhängigkeit von der Verzögerungszeit. Diese wurden in vier verschiedene Regimes eingeteilt und den physikalischen Prozessen während der Laserablation zugeordnet: Die Thermalisierung des Elektronensystems für Pulsabstände unter 1 ps, Schmelzprozesse zwischen 1 und 10 ps, der Beginn des Abtrags nach mehreren 10 ps und die Expansion der Plasmawolke nach über 100 ps. Dabei wird das LIPS-Signal effizient verstärkt und bei 800 ps maximal. Die Lochdurchmesser ändern sich als Funktion des Pulsabstands wenig im Vergleich zur Tiefe. Die gesamte Abtragsrate variiert um maximal 50 %, während sich das LIPS-Signal vervielfacht: Für Ti und TiAl typischerweise um das Dreifache, für Al um das 10-fache. Die gemessenen Transienten zeigen eine hohe Reproduzierbarkeit, jedoch kaum eine Energie- bzw. materialspezifische Abhängigkeit. Mit diesen Ergebnissen wurde eine gezielte Optimierung der DP-LIPS-Parameter an Al durchgeführt: Bei einem Pulsabstand von 800 ps und einer Gesamtenergie von 65 nJ (vierfach über der Ablationsschwelle) wurde eine 40-fache Signalerhöhung bei geringerem Rauschen erzielt. Die Lochdurchmesser vergrößerten sich dabei um 44 % auf (650±150) nm, die Lochtiefe um das Doppelte auf (100±15) nm. Damit war es möglich, die spektrochemische Sensitivität von fs-LIPS zu erhöhen und gleichzeitig die hohe räumliche Auflösung aufrecht zu erhalten.

The Surface Plasmon Resonance of Supported Noble Metal Nanoparticles: Characterization, Laser Tailoring, and SERS Application

This work deals with the optical properties of supported noble metal nanoparticles, which are dominated by the so-called Mie resonance and are strongly dependent on the particles’ morphology. For this reason, characterization and control of the dimension of these systems are desired in order to optimize their applications. Gold and silver nanoparticles have been produced on dielectric supports like quartz glass, sapphire and rutile, by the technique of vapor deposition under ultra-high vacuum conditions. During the preparation, coalescence is observed as an important mechanism of cluster growth. The particles have been studied in situ by optical transmission spectroscopy and ex situ by atomic force microscopy. It is shown that the morphology of the aggregates can be regarded as oblate spheroids. A theoretical treatment of their optical properties, based on the quasistatic approximation, and its combination with results obtained by atomic force microscopy give a detailed characterization of the nanoparticles. This method has been compared with transmission electron microscopy and the results are in excellent agreement. Tailoring of the clusters’ dimensions by irradiation with nanosecond-pulsed laser light has been investigated. Selected particles are heated within the ensemble by excitation of the Mie resonance under irradiation with a tunable laser source. Laser-induced coalescence prevents strongly tailoring of the particle size. Nevertheless, control of the particle shape is possible. Laser-tailored ensembles have been tested as substrates for surface-enhanced Raman spectroscopy (SERS), leading to an improvement of the results. Moreover, they constitute reproducible, robust and tunable SERS-substrates with a high potential for specific applications, in the present case focused on environmental protection. Thereby, these SERS-substrates are ideally suited for routine measurements.

Precision laser spectroscopy of the ground state hyperfine splitting of hydrogenlike ^209 Bi^82+

The first direct observation of a hyperfine splitting in the optical regime is reported. The wavelength of the M1 transition between the F = 4 and F = 5 hyperfine levels of the ground state of hydrogenlike ^209 Bi^82+ was measured to be \lamda_0 = 243.87(4) nm by detection of laser induced fluorescence at the heavy-ion storage ring ESR at GSI. In addition, the lifetime of the laser excited F = 5 sublevel was determined to be \tau_0 = 0.351(16) ms. The method can be applied to a number of other nuclei and should allow a novel test of QED corrections in the previously unexplored combination of strong magnetic and electric fields in highly charged ions.

Squeezed thermal phonons precurse nonthermal melting of silicon as a function of fluence

A femtosecond-laser pulse can induce ultrafast nonthermal melting of various materials along pathways that are inaccessible under thermodynamic conditions, but it is not known whether there is any structural modification at fluences just below the melting threshold. Here, we show for silicon that in this regime the room-temperature phonons become thermally squeezed, which is a process that has not been reported before in this material. We find that the origin of this effect is the sudden femtosecond-laser-induced softening of interatomic bonds, which can also be described in terms of a modification of the potential energy surface. We further find in ab initio molecular-dynamics simulations on laser-excited potential energy surfaces that the atoms move in the same directions during the first stages of nonthermal melting and thermal phonon squeezing. Our results demonstrate how femtosecond-laser-induced coherent fluctuations precurse complete atomic disordering as a function of fluence. The common underlying bond-softening mechanism indicates that this relation between thermal squeezing and nonthermal melting is not material specific.

Electronic origin of bond softening and hardening in femtosecond-laser-excited magnesium

Many ultrafast structural phenomena in solids at high fluences are related to the hardening or softening of particular lattice vibrations at lower fluences. In this paper we relate femtosecond-laser-induced phonon frequency changes to changes in the electronic density of states, which need to be evaluated only in the electronic ground state, following phonon displacement patterns. We illustrate this relationship for a particular lattice vibration of magnesium, for which we—surprisingly—find that there is both softening and hardening as a function of the femtosecond-laser fluence. Using our theory, we explain these behaviours as arising from Van Hove singularities: We show that at low excitation densities Van Hove singularities near the Fermi level dominate the change of the phonon frequency while at higher excitations Van Hove singularities that are further away in energy also become important. We expect that our theory can as well shed light on the effects of laser excitation of other materials.

Effect of dressing on high-order harmonic generation in vibrating H2 molecules

Wir entwickeln die Starkfeldnäherung für die Erzeugung hoher Harmonischer in Wasserstoffmolekülen, wobei die Vibrationsbewegung berücksichtigt wird, sowie die laserinduzierte Kopplung zwischen den beiden untersten Born-Oppenheimer-Zuständen im Molekülion, das durch die anfängliche Ionisation des Moleküls erzeugt wird. Wir zeigen, dass die Kopplung bei längeren Laserwellenlängen (≈ 2 μm) wichtig wird und zu einer Reduzierung der Erzeugung von Harmonischen führt, sowie zu einer Änderung des Verhältnisses von Harmonischen in verschiedenen Isotopen. ----------------------------------------------------------------------- We develop the strong-field approximation for high-order harmonic generation in hydrogen molecules, including the vibrational motion and the laser-induced coupling of the lowest two Born-Oppenheimer states in the molecular ion that is created by the initial ionization of the molecule. We show that the field dressing becomes important at long laser wavelengths (≈ 2 μm), leading to an overall reduction of harmonic generation and modifying the ratio of harmonic signals from different isotopes.

Femtosecond real-time probing of reactions. 12. Vectorial dynamics of transition states

Femtosecond time-resolved techniques with KETOF (kinetic energy time-of-flight) detection in a molecular beam are developed for studies of the vectorial dynamics of transition states. Application to the dissociation reaction of IHgI is presented. For this system, the complex [I---Hg---I](++)* is unstable and, through the symmetric and asymmetric stretch motions, yields different product fragments: [I---Hg---I](++)* -> HgI(X^2/sigma^+) + I(^2P_3/2) [or I*(^2P_l/2)] (1a); [I---Hg---I](++)* -> Hg(^1S_0) + I(^2P_3/2) + I(^2P_3/2) [or I* (^2P_1/2)] (1 b). These two channels, (1a) and (1b), lead to different kinetic energy distributions in the products. It is shown that the motion of the wave packet in the transition-state region can be observed by MPI mass detection; the transient time ranges from 120 to 300 fs depending on the available energy. With polarized pulses, the vectorial properties (transition moments alignment relative to recoil direction) are studied for fragment separations on the femtosecond time scale. The results indicate the nature of the structure (symmetry properties) and the correlation to final products. For 311-nm excitation, no evidence of crossing between the I and I* potentials is found at the internuclear separations studied. (Results for 287-nm excitation are also presented.) Molecular dynamics simulations and studies by laser-induced fluorescence support these findings.