Next Article in Journal
Compact 5G Nonuniform Transmission Line Interdigital Bandpass Filter for 5G/UWB Reconfigurable Antenna
Next Article in Special Issue
Sawing Force Prediction Model and Experimental Study on Vibration-Assisted Diamond Wire Sawing
Previous Article in Journal
Evaluation of Three-Dimensional Surface Roughness in Microgroove Based on Bidimensional Empirical Mode Decomposition
Previous Article in Special Issue
Effect of Punch Surface Microtexture on the Microextrudability of AA6063 Micro Backward Extrusion
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Simulation of Laser-Heating and Energetic Plasma Plume Expansion in Pulsed Laser Deposition of Y3Fe5O12

by
Joko Suwardy
1,*,
Muzakkiy Putra Muhammad Akhir
1,
Robi Kurniawan
2,
Beni Rio Hermanto
3,
Isa Anshori
3 and
Mohammad Hamzah Fauzi
1,4
1
Research Center for Quantum Physics, National Research and Innovation Agency (BRIN), South Tangerang 15314, Indonesia
2
Department of Physics, Faculty of Mathematics and Natural Sciences, Universitas Negeri Malang, Jl. Semarang No. 5, Malang 65145, Indonesia
3
Biomedical Engineering Department, School of Electrical Engineering and Informatics, Bandung Institute of Technology, Bandung 40132, Indonesia
4
Research Collaboration Center for Quantum Technology 2.0, Bandung 40132, Indonesia
*
Author to whom correspondence should be addressed.
Micromachines 2022, 13(11), 2012; https://doi.org/10.3390/mi13112012
Submission received: 27 September 2022 / Revised: 10 November 2022 / Accepted: 11 November 2022 / Published: 18 November 2022
(This article belongs to the Special Issue Advanced Manufacturing Technology and Systems, 2nd Edition)

Abstract

:
In the present study, numerically iterative models are employed to study two processes involved in the pulsed laser deposition of an Y3Fe5O12 target. The 1D conduction heat model is used to evaluate the temperature of the target irradiated by a nano-second pulse laser, taking into account the plasma shielding effect. Further, the gas dynamics model is employed to simulate the kinetic of plasma plume expansion. The results may be important in obtaining high-quality Y3Fe5O12 thin films.

1. Introduction

The ferrimagnetic insulator yttrium iron garnet (Y3Fe5O12), often abbreviated as YIG, has become an attractive material in the field of quantum spintronics. Having a low Gilbert damping value, YIG is widely used to investigate various quantum phenomena, such as spin waves, spin pumping, the inverse spin Hall effect, and spin caloritronic [1,2,3,4]. A high-quality YIG thin film is required to observe these quantum effects.
One of the most powerful techniques to grow complex oxide thin films, such as YIG, is pulsed laser deposition (PLD) [5]. A recent report claimed a successful deposition of YIG thin film on a gallium gadolinium garnet (GGG) substrate using room temperature PLD followed by a post-annealing treatment. This YIG has a monocrystalline structure with a very low damping parameter ( α = 6 × 10 5 ) [6]. However, the use of garnet substrates such as GGG limits its practical application. Considerable efforts have been put forward to grow YIG on commercially compatible chip substrates, e.g., silicon [7,8,9]. However, its properties still fall behind that of YIG/GGG due to the smaller grain size and cracking of the film [7,8].
Although the procedure of obtaining laser-ablated material is simple, involving a focus pulsed laser with an intensity above the threshold ablating the material target, the physics incorporated therein is overly complex. Various interdependence processes take place such as target-heating, melting, vaporization, ionization, plasma formation, plasma hydrodynamic expansion, shock wave generation, and so on [10,11,12]. In general, based on the interactions, the processes can be divided into three major zones:
  • Interaction between the laser and the target material, which results in evaporation on the target surface [13,14,15,16,17,18],
  • Interaction between the material gas and the laser, which produces plasma [19,20], and
  • Plasma plume expansion [21,22,23].
The first two interactions happen during the laser pulse, while the latter zone occurs when the laser irradiation has stopped.
In this study, we employed numerically iterative models to investigate two crucial processes in the pulsed laser deposition: laser-heating of the target and plasma plume expansion. These two processes govern the kinetic of particles deposited on the substrate. In the calculation, we employed the thermophysical properties of the YIG target. The results are applicable for experimentally grown YIG thin film on commercially compatible chip substrates, e.g., silicon substrate.

2. Theoretical Formulation

2.1. Laser-Target Interaction

Figure 1 illustrates the schematic of the laser heating the target material. A single pulse of nano-second laser was focused onto the YIG target surface. Upon receiving the pulse, the target absorbed some portion of laser energy and reflected the rest. The absorbed laser energy can be expressed as the following source term:
S ˙ = μ ( 1 R ) I ( t ) e μ z ,
where R is the optical reflectivity of YIG and μ = 4 π k / λ , is the absorption coefficient of YIG. For simplicity, R is assumed to be independent of temperature. Here, λ = 532 nm is the laser wavelength and k is the extinction coefficient. The temporal pulse intensity I ( t ) is expressed in a Gaussian distribution:
I ( t ) = A e x p ( 4 ln ( 2 ) ( t t c ) 2 w 2 ) w π 4 l n ( 2 ) ,
where w is the full-width half maximum (FWHM) of the laser pulse, A is the area and t c is the center of the Gaussian profile. Here, we used the property of Nd:YAG Continuum laser Surelite III, which is available in our PLD system, with a laser energy of 446 mJ, beam diameter ( d ) = 9.5 mm, and w is 5 ns. The calculated maximum intensity of 1.26 × 10 8 W/cm2 and pulse duration, t L = 15 ns is obtained.
The heat transfer of the nano-second laser-radiated YIG before melting (tm) can be expressed by a 1D heat conduction equation as follows [13]:
C s ρ s T ( t , z ) t = z   ( K s T ( t , z ) z ) + S ˙ ,
where C s , ρ s ,   K s , are the density, specific heat, and thermal conductivity of the solid phase YIG, respectively. As shown in Figure 1, z = 0 at the surface of the target. We note that in the above model we assumed that the enthalpy of fusion Δ H m = 0 and the interface laser-target interaction was static. Those assumptions are safe considering the aforementioned laser properties and that only a single pulse are used in this simulation. Moreover, the following initial and boundary conditions are used:
T   ( t , z ) t = 0 = T 0 ,
K s T ( t , z ) z | z = 0 = S ˙ μ ,               T T m
K s T ( t , z ) z | z = d = 0 ,             T T m
where Equation (4) is the initial condition, T 0 is the initial temperature, and Equations (5) and (6) are the boundary conditions at the top ( z = 0 ) and rear ( z = d ) surfaces of the YIG target, respectively.
Above the melting temperature, T > T m , the target surface will melt and change its physical properties to a liquid phase. Furthermore, if the laser intensity is sufficient, the target surface will reach the boiling point and evaporation will occur. The fluence threshold for YIG can be calculated as follows [24]:
F t h ρ s C s ( T m T 0 ) L t h ( 1 R ) ,
by taking the maximum of the laser penetration depth, L t h = 1   μ m , we obtained the fluence threshold ( F t h ) of YIG 0.54 J/cm2, which is lower than the laser fluence of our system. Thus, we can expect that evaporation will occur during the single-pulsed laser irradiation. The heat conduction on the target for t > t m until the end of the pulse t = t L can be expressed as followed:
C l ρ l T ( t , z ) t C l ρ l V r T ( t , z ) z = z   ( K s T ( t , z ) z ) + S ˙ ,
where C l , ρ l ,   K l , are the density, specific heat, and thermal conductivity of the liquid phase of YIG, respectively. Above the boiling point, the recession rate of evaporated material, V r , can be described by the Hertz–Knudsen equation [13]:
V r   P ρ s 2 π k B T S m ( 1 β ) ,
where T S is the surface target temperature, m is the average atom mass (here, we assumed m = 89 by the mass of the Y atom), β is the fraction of vapor particles that return to the surface, and k B is the Boltzmann constant. Moreover, the pressure ( P ) of ablated material as a function of surface target temperature can be expressed as [13]:
P = P b e x p { Δ H v m k B ( 1 T b 1 T ) } ,
where Δ H v is the enthalpy of vaporization at boiling point T b and atmospheric pressure P b (1 atm.). Furthermore, the boundary conditions are expressed as:
K l T ( z , t ) z | z = 0 = Δ H v ρ l V r + S ˙   μ ,
K l T ( z , t ) z | z = d = 0
Upon vaporization, the vapor material will also absorb the laser energy, causing its temperature to increase, and ionization will occur, leading to the formation of the plasma plume. The absorption of laser energy by the vapor material is dominated by inverse bremsstrahlung (IB) and the photoionization mechanism (PI) [24]. The IB process involves the absorption of laser energy by a free electron, while the PI happens when the electron is excited directly due to laser energy absorption by neutral atoms. These two mechanisms lead to various species incorporated within the plasma plume, such as electrons, ions, neutrals, clusters, and particulates.
At the end of the laser pulse, plasma with a particular height H will be formed. Due to the density and absorption of the plasma plume, now the laser energy that reaches the surface of the target will be reduced, which is known as the plasma shielding effect [20]. Therefore, the source term can be written as followed:
S ˙ = μ ( 1 R ) I ( t ) e μ z e ( α I B + α P I ) H ,
where α I B and α P I are the absorption coefficients of the IB and PI absorption mechanisms, respectively. The absorption coefficient of IB is given by [25]:
α I B = 3.69 × 10 8 ( Z I B 3 N I B 2 T v 0.5 ν 3 ) { 1 e ( h ν k b T v ) } ,
where Z I B , N I B , h, T v , and ν are average charge, ion density, Plank constant, vapor temperature, and frequency of laser, respectively. Moreover, the absorption coefficient of IB is given by [25]:
α P I = σ P I N P I ,
where N P I is the neutral atom density and σ P I is the cross-section between an excited neutral atom and a photon that is involved in the photoionization process. The typical value of σ P I is in the range of 10−21 m2 [26].

2.2. Plasma Plume Expansion

After the termination of the laser pulse ( t = t L ), the existing plume will expand away from the target surface due to the conversion of high thermal energy and energy stored as excitation and ionization in the plasma to kinetic energy [10]. In nano-second laser irradiation and with the laser properties under consideration, where the repetition rate is 10 Hz, we can assume the plume expansion time to be far longer than its formation time (close to pulsed duration, ~ 15 ns). Thus, the plume expansion can be considered independent of its formation.
Figure 2 displays the schematic of the plasma plume expansion after the termination of the laser pulse, showing the change in plume dimension over time. The initial plasma with dimensions of height: Z 0 and radius: X 0 and Y 0 expands as a semi-ellipsoid with a front determined by the axes X ( t ) ,   Y ( t ) and Z ( t ) . Here, the radius can be approximated by the laser beam radius, X 0 = Y 0 = 4.75 mm. The height of the initial plume is much smaller than the radial dimension ( Z 0 X 0 ,   Y 0 ) . The height of the plume in z-direction can be approximated by Z 0 υ s × t L where υ s is the sound velocity in the plume [27]. When the plasma plume is modeled as an ideal gas, the speed velocity can be calculated as [25]:
υ s = γ k B T P m ,
where γ is the specific heat ratio or adiabatic index, k B is the Boltzmann constant, and T P is the initial plasma plume temperature. For a monoatomic ideal gas, γ = 5/3; however, in laser-produced plasma, the value of γ is expected to decrease because of the high plasma temperature and high degree of excitation, as well as the ionization of the plasma species. The value of γ is estimated in the range of 1.2–1.3 [21,28].
In this study, we simulate the plume expansion under vacuum conditions. The plasma plume is modeled as an ideal gas that undergoes adiabatic expansion. The expansion of the plasma plume in Cartesian coordinate can be expressed as [27]:
Z d 2 Z d t 2 = X d 2 X d t 2 = X d 2 X d t 2 = ( 5 γ 3 ) ( E M ) ( X 0 Y 0 Z 0 X Y Z ) γ 1 ,
where γ = C p / C v is the adiabatic constant, and X, Y, and Z are the dimensions of the plume as a function of t in the axes of x, y, and z, respectively. E is the initial energy and M is the mass of the plume, where the ratio can be approximated by E M υ s 2 = γ k B T P m . The model is based on Lie group transformation theory, where the solution of the gas dynamic equation is simplified. In this model, the density and the pressure of the plume are constant on the ellipsoidal surface, e.g., x 2 X ( t ) 2 +   y 2 Y ( t ) 2 + z 2 Z ( t ) 2 = constant. Furthermore, the hydrodynamic motion of all particles in the plume is governed by self-similarity, such that the velocity is controlled by the relative position of the edge [10,23]. Therefore, the velocity distribution in x, y, and z-axes can be calculated as:
v x = x X d X ( t ) d t , v y = y Y d Y ( t ) d t , v z = z Z d Z ( t ) d t

3. Numerical Implementation

3.1. Laser Parameters and Material Properties

Table 1 and Table 2 summarize the laser parameters, and the physical properties of YIG target materials.
Several properties of liquid-YIG that cannot be found in references are set to be the same values as in the solid phase. To the best of our knowledge, the enthalpy of vaporization Δ H v has not been determined experimentally. Therefore, we set the value to be in the order of the corresponding oxides.

3.2. Simulation of the Temperature of the YIG Target

To simulate the effect of laser pulse radiation on the temperature of the YIG target, we solve the 1D heat conduction partial differential equation (Equation (3)) by using the finite different numerical (FDM) method. Two source terms are assigned to each equation, for temperature below melting point ( T T m ) and above melting point ( T > T m ) , expressed in Equation (1) and Equation (13), respectively.
Two different domains are employed, where n is for time domain (t) and j is for the depth from the surface of the target domain (z). Based on the known initial condition at t = 0, where the initial temperature T = 300 K is independent of depth, we can use a center approach for the z-domain; meanwhile, the forward approach is applicable for the t-domain. The number of partitions for the time and depth domains are assigned to M and N, respectively. The number of partitions determines the calculation results, where M >> N led to a convergence value.
Numerical expression for heat conduction equation below melting point ( t t m ) :
T n + 1 , j = T n , j + Δ t   C s ρ s { K s Δ z 2   ( T n , j + 1 2 T n , j + T n , j 1 ) + ( 1 R ) μ I n exp ( μ z ) } ,
with boundary conditions:
T n + 1 , 1 = T n + 1 , 2 + Δ z K s ( 1 R ) I n ,
T n + 1 , N = T n + 1 , N 1 ,
and initial temperature T 1 , j = 300 K. Furthermore, for temperatures above the melting point ( t M < t < t L ) , the numerical equation can be express as:
T n + 1 , j = T n , j + Δ t   C l ρ l { K s Δ z 2   ( T n , j + 1 2 T n , j + T n , j 1 ) + ( C l ρ l V r ( T n , 1 ) + ( 1 R ) μ I n exp ( μ z ) exp ( ( α I B n + α P I n ) H ) }
V r   P b e x p { Δ H v ( T b ) m k B ( 1 T b 1 T n , 1 ) } ρ ( 2 π k B T n , 1 m ) 1 2 ( 1 β )
with boundary condition:
T n + 1 , 1 = T n + 1 , 2 + Δ z K s ( 1 R ) I n exp ( μ z ) exp ( ( α I B n + α P I n ) H Δ H v ρ l V r ,
T n + 1 , N = T n + 1 , N 1 ,
and the initial temperatures for each depth are set to the melting temperature, t = t m , T 1 , j = T m , j ,. The shielding effect of the plasma plume decreases the laser incident intensity which reaches the surface target. The rate of decrease is determined by the absorption coefficient IB ( α I B ) and the PI ( α P I ) mechanism, which can be expressed as:
α I B ( n ) = 3.69 × 10 8 ( Z I B 3 N I B 2 T v n , 1 0.5 ν 3 ) { 1 e ( h ν k b T v n , 1 ) } ,
α P I ( n ) = N P I σ P I , x ,
where for such local thermal equilibrium (LTE) of the plume, the average charge Z I B 2 [20], υ = 5.64 × 10 14 , and T v is assumed to be the same with the boiling temperature Tb, while σ P I , x is set to 10−21 m2. In this calculation, the density of the ion ( N I B ) and the neutral atom density ( N P I ) are set as a function of time, where the maximum values are reached at the end of the laser pulse ( t L ) , as followed:
N I B ( n ) = N P I ( n ) = ( n 1 ) 10 ( N m a x ( M 1 ) 10 ) ,
where N m a x is set to 1021/cm3. As shown in Equation (13), the decrease in laser intensity is also a function of plume height (H), which is expressed as:
H ( n ) = ( n 1 ) ( H m a x M 1 ) ,
where H m a x is the maximum height of the plume. Just as for N I B and N P I , the height of the plume is set to reach the maximum at the end of the laser pulse.

3.3. Simulation of Plasma Plume Expansion

The geometry of plasma plume expansion after the termination of the laser pulse is depicted in Figure 2. We calculate dimension and edge velocity for the equation plume expansion (Equation (17)) by using the ordinary differential equation, ode45 function in MATLAB®. The initial plume dimension X 0 = Y 0 is estimated by a laser beam radius of 4.75 mm, while Z 0 = H m a x is estimated from sound velocity (Equation (16)), which is dependent of the initial plasma temperature (TP). Since the initial plume has much larger thermal energy than the kinetic energy, the initial condition for the edge velocity can be set as:
d X ( 0 ) d t = d Y ( 0 ) d t = d Z ( 0 ) d t = 0 ,
and the calculation is conducted for an expansion time of 10 µm with time step Δ t = 0.1 ns.
As one can see from Equation (17), the values of γ ,   E , and M are critical to obtaining the correct calculation result. However, to the best of our knowledge, those values have not been experimentally confirmed yet. Therefore, we input the estimation values for γ ,   E , and M in our calculation. The ratio of E M ν s 2 , which depends on the temperature of the initial plasma. We evaluate the expansion by varying the TP. The range for the initial plasma temperature is set to 7000 K–20,000 K, referring to the previous study of nano-second laser ablation of a YBa2Cu3O7 target [21].

4. Result and Discussion

The theoretical models consisting of laser-target heating (Equations (3) and (8)) and plume expansion (Equation (17)) were numerically simulated. The thermophysical properties of YIG used in the calculation are summarized in Table 2. Further, Table 3 summarizes the parameters of the plasma plume used in the calculation:
We choose one value of γ = 1.2 and varied the plasma temperatures, resulting in the three different values for height and the ratio of the energy and mass of the initial plasma plume ( t = t L ).
Figure 3 shows the surface temperature of the YIG target during a single pulse shot, where the left y-axis denotes the surface temperature (T), while the right y-axis denotes the intensity of laser (I(t)). For this calculation, we employed γ = 1.2 and T P = 7000 K, resulting a maximum height of the plume H m a x = 13.29 µm. The temperature rises because of the absorption of the laser energy by the target and reaches the melting point ( T m = 1828 K) at t = 3.8 ns. At this region (t = 0–3.8 ns), the calculation corresponds to the first 1D heat conduction equation (Equation (3)). After the target surface liquidizes, the surface temperature rises further and reaches the boiling point ( T m = 3611 K), when t = 4.8 ns. It can be seen that the temperature of target surface increases further, surpassing the boiling point and reaching a maximum of about T = 5250 K before decreasing exponentially. It is interesting to see that the surface temperature is saturated even before the laser reach its maximum intensity. This result implies that the shielding effect of the plasma plume reduces the intensity of laser to reach the target surface. In the calculation, we included the term for plasma absorption, which reduced the source heat term used in Equation (8) by rate of e ( α I B + α P I ) H . Furthermore, we also considered the liquid-vapor phase change by taking into account the recession rate of evaporation material v r (Equation (9)), which also contributed to the temperature decrease at the surface after evaporation occurred.
The models with the plasma shielding effect in nano-second laser solid target ablation had been previously developed [20,26]. Zhang et al. [20] evaluate the accuracy of the model with the ablation depth in the target experimentally and compare their model to the model with no plasma shielding effect of Singh et al. [21]. The model with the plasma shielding effect gives better accuracy (~3% error) compared to that of Singh et al. (~10% error) [12,20,21]. Therefore, in the above YIG target temperature model, we incorporated the plasma shielding effect, based on IB and the PI mechanism for moderate laser fluence. However, the accuracy of our calculation still needs confirmation from the experiment, which yet remains for further study.
The distribution of temperature inside the YIG target during the laser pulse duration is shown in Figure 4. Temperature distribution inside the target material during laser pulse irradiation (t ≤ 15 ns) is calculated for a maximum depth z = 10 µm. The range for depth, z, is set to 10 µm. The curve inside the graph represents the isothermal boundary line. As expected, the maximum temperature is found at the surface of the target. The melting depth reached ~ 1 µm from the surface, as shown by the isothermal line for T > 1828 K. Even though the surface temperature drops significantly after t = 7 ns, a significant region remains molten for the rest of the laser pulse. Furthermore, the evaporation window, when the temperature target surpassed the melting point, is obtained as deep as ~ 0.4 µm in the 4.8–8 ns radiation period, at which the laser radiation reaches the maximum ( t = 7.7   ns ). This result implies that the supply of evaporated material available to form a plasma is limited to that period. Finally, during laser radiation, the conduction of heat can reach a depth of approximately 7 µm, indicated by the isothermal line of 310 K.
The calculation of the plume expansion (represented by Equation (17)) for three different initial plume temperatures, TP = 20,000, 15,000, and 7000 K, is shown in Figure 5. Different initial temperatures of the plasma plume resulted in different initial plume heights ( Z 0 ), as summarized in Table 1, where a higher plasma temperature led to a higher initial plume height. The edge velocity in the z-direction, or the so-called center of mass velocity ( V z ), is shown in Figure 5a. It can be seen that the edge velocity increases rapidly at an early stage of expansion (t = 7 ns) before reaching an asymptotic value. The highest velocity is obtained when TP = 20,000 K, with a maximum of approximately 7 × 10 5 cm/s. Further, Figure 5b shows the evolution of the plume height. Similar dependency on the temperature of the plasma is also observed, where at the end of calculation time of 10 µm, a maximum plume height of 7 cm is obtained for TP = 20,000 K.
Figure 6 shows the velocity distributions inside the plume for each plasma temperature when the plume height is 4 cm. The red dashed line denotes the plume dimensions. The velocity shown here is the resultant velocity V R = s q r t   ( V x 2 + V y 2 + V z 2 ) . Self-similar gas dynamics is used to calculate the velocity distribution from the edge velocity, as express in Equation (18). It can be seen that, at a height of 4 cm, the radius of plume only reached approximately 1 cm, which implies that the edge velocity in the x, y-direction is much smaller than that of z-direction. In fact, for TP = 20,000 K, the ratio V x / V z is 0.2 when the plume height is 4 cm. This edge velocity difference in the radial and perpendicular direction leads to a semi-ellipsoid shape of expansion of the plume, elongated perpendicularly to the sample surface. Furthermore, the dependence of the plasma temperature clearly shows where a higher plasma temperature led to a higher velocity distribution inside the plume, where for TP = 20,000, 15,000, and 7000 K, the maximum resultant velocities, VR, of approximately 7, 6.1, and 4.2 × 105 cm/s are achieved, respectively.
It was previously demonstrated experimentally, in SrTiO3 grown by PLD on a Si substrate, that a high energy of plasma can induce visible burning and destroy the crystallinity of the substrate [33]. Even though, after post-annealing, the Si substrate crystallinity showed improvement, the existence of secondary phases was observed on the grown SrTiO3 film. In the study considered here, the plasma plume expansion in vacuum was simulated, as shown in Figure 6, wherein the velocity of the ablated particles was calculated. However, the kinetic energy will be dependent on the ablated species and its respective density, which should be evaluated further by direct measurement, for instance, by using laser-induced breakdown spectroscopy (LIBS) [11].

5. Summary

In summary, we have studied the effect of nano-laser radiation on the temperature of the YIG target. Within the laser parameters under consideration, a single pulse of laser increases the YIG target’s temperature to above the boiling temperature and generates a plasma. Furthermore, the expansion of the plasma plume in vacuum was also simulated. It shows that the kinetic of the expanding plasma plume strongly depends on the temperature of the initial plasma plume. These results can be used as an input to set the growth parameters of pulsed laser deposition of YIG thin films, such as the laser fluence and the substrate-target distances.

Author Contributions

J.S. is the main contributor. Conceptualization, J.S.; data curation, J.S.; formal analysis, J.S. and M.H.F.; funding acquisition, J.S., M.H.F., B.R.H., and I.A.; investigation, J.S.; methodology, J.S.; project administration, J.S, M.H.F., B.R.H., and I.A.; software, J.S., B.R.H., and I.A.; supervision, J.S. and M.H.F.; validation, J.S. and M.P.M.A.; visualization, J.S.; writing—original draft, J.S.; writing—review and editing, J.S., R.K., M.P.M.A., and M.H.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partly funded by “Rumah Program Sains Fundamental Molekuler 2022”, proposal code No. 132, Research Organization for Nanotechnology and Materials, National Research and Innovation Agency (BRIN), Indonesia. The authors thank the School of Electrical Engineering and Informatics, Bandung Institute of Technology, for funding the article processing charge.

Data Availability Statement

Data available upon request from the corresponding author.

Acknowledgments

J.S. would like to thank Donny Dwiputra of the Theoretical Physics Laboratory of the Bandung Institute of Technology and Hairil Anwar of the Bandung Institute of Technology for the discussions of code development for the programming.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kajiwara, Y.; Harii, K.; Takahashi, S.; Ohe, J.; Uchida, K.; Mizuguchi, M.; Umezawa, H.; Kawai, H.; Ando, K.; Takanashi, K.; et al. Transmission of electrical signals by spin-wave interconversion in a magnetic insulator. Nature 2010, 464, 262–266. [Google Scholar] [CrossRef] [PubMed]
  2. Uchida, K.-I.; Xiao, J.; Adachi, H.; Von Der Ohe, J.; Takahashi, S.; Ieda, J.; Ota, T.; Kajiwara, Y.; Umezawa, H.; Kawai, H.; et al. Spin Seebeck insulator. Nat. Mater. 2010, 9, 894–897. [Google Scholar] [CrossRef]
  3. Kelly, O.D.; Anane, A.; Bernard, R.; Ben Youssef, J.; Hahn, C.; Molpeceres, A.H.; Carrétéro, C.; Jacquet, E.; Deranlot, C.; Bortolotti, P.; et al. Inverse spin Hall effect in nanometer-thick yttrium iron garnet/Pt system. Appl. Phys. Lett. 2013, 103, 082408. [Google Scholar] [CrossRef] [Green Version]
  4. Jungfleisch, M.B.; Lauer, V.; Neb, R.; Chumak, A.; Hillebrands, B. Improvement of the yttrium iron garnet/platinum interface for spin pumping-based applications. Appl. Phys. Lett. 2013, 103, 022411. [Google Scholar] [CrossRef]
  5. Chrisey, D.B.; Hubler, G.K. Pulsed Laser Deposition of Thin Films; John and Wiley and Sons: Hoboken, NJ, USA, 1994. [Google Scholar]
  6. Hauser, C.; Richter, T.; Homonnay, N.; Eisenschmidt, C.; Qaid, M.; Deniz, H.; Hesse, D.; Sawicki, M.; Ebbinghaus, S.G.; Schmidt, G. Yttrium Iron Garnet Thin Films with Very Low Damping Obtained by Recrystallization of Amorphous Material. Sci. Rep. 2016, 6, 20827. [Google Scholar] [CrossRef] [Green Version]
  7. Capku, Z.; Yildiz, F. Spin wave modes observation in YIG thin films with perpendicular magnetic anisotropy. J. Magn. Magn. Mater. 2021, 538, 168290. [Google Scholar] [CrossRef]
  8. Balinskiy, M.; Ojha, S.; Chiang, H.; Ranjbar, M.; Ross, C.A.; Khitun, A. Spin wave excitation in sub-micrometer thick Y3Fe5O12 films fabricated by pulsed laser deposition on garnet and silicon substrates: A comparative study. J. Appl. Phys. 2017, 122, 123904. [Google Scholar] [CrossRef]
  9. Kumar, N.; Misra, D.; Venkataramani, N.; Prasad, S.; Krishnan, R. Magnetic properties of pulsed laser ablated YIG thin films on different substrates. J. Magn. Magn. Mater. 2004, 272–276 (Suppl. S1), E899–E900. [Google Scholar] [CrossRef]
  10. Schou, J.; Amoruso, S.; Lunney, J.G. Plume Dynamics. In Laser Ablation and Its Applications; Springer: Boston, MA, USA, 2007; Volume 129, pp. 67–95. [Google Scholar] [CrossRef]
  11. Harilal, S.S.; Phillips, M.C.; Froula, D.H.; Anoop, K.K.; Issac, R.C.; Beg, F.N. Optical diagnostics of laser-produced plasmas. Rev. Mod. Phys. 2022, 94, 035002. [Google Scholar] [CrossRef]
  12. Marla, D.; Bhandarkar, U.V.; Joshi, S.S. Critical assessment of the issues in the modeling of ablation and plasma expansion processes in the pulsed laser deposition of metals. J. Appl. Phys. 2011, 109, 021101. [Google Scholar] [CrossRef]
  13. Peterlongo, M.A.; Kelly, R. Laser-pulse sputtering of aluminum: Vaporization, boiling, superheating, and gas-dynamic effects. Phys. Rev. E 1994, 50, 4716. [Google Scholar] [CrossRef] [PubMed]
  14. Wang, B.; Huang, Y.; Jiao, J.; Wang, H.; Wang, J.; Zhang, W.; Sheng, L. Numerical Simulation on Pulsed Laser Ablation of the Single-Crystal Superalloy Considering Material Moving Front and Effect of Comprehensive Heat Dissipation. Micromachines 2021, 12, 225. [Google Scholar] [CrossRef]
  15. Li, X.; Xu, W.; Wang, Y.; Zhang, X.; Hui, Z.; Zhang, H.; Wageh, S.; Al-Hartomy, O.A.; Al-Sehemi, A.G. Optical-intensity modulators with PbTe thermoelectric nanopowders for ultrafast photonics. Appl. Mater. Today 2022, 28, 101546. [Google Scholar] [CrossRef]
  16. Feng, J.; Li, X.; Zhu, G.; Wang, Q.J. Emerging high-performance Sns/CdS nanoflowers heterojunction for ultrafast photonics. ACS Appl. Mater. Interfaces 2020, 12, 43098–43105. [Google Scholar] [CrossRef] [PubMed]
  17. Zhang, Y.; Li, X.; Qyyum, A.; Feng, T.; Guo, P.; Jiang, J.; Zheng, H. PbS Nanoparticles for Ultrashort Pulse Generation in Optical Communication Region. Part. Part. Syst. Charact. 2018, 35, 1800341. [Google Scholar] [CrossRef]
  18. Liu, J.; Li, X.; Guo, Y.; Qyyum, A.; Shi, Z.; Feng, T.; Zhang, Y.; Jiang, C.; Liu, X. SnSe2 nanosheets for subpicosecond harmonic mode-locked pulse generation. Small 2019, 15, 1902811. [Google Scholar] [CrossRef]
  19. Lunney, J.G.; Jordan, R. Pulsed laser ablation of metals. Appl. Surf. Sci. 1998, 127–129, 941–946. [Google Scholar] [CrossRef]
  20. Zhang, D.; Liu, D.; Li, Z.; Hou, S.; Yu, B.; Guan, L.; Tan, X.; Li, L. A new model of pulsed laser ablation and plasma shielding. Phys. B Condens. Matter 2005, 362, 82–87. [Google Scholar] [CrossRef]
  21. Singh, R.K.; Narayan, J. Pulsed-laser evaporation technique for deposition of thin films: Physics and theoretical model. Phys. Rev. B 1990, 41, 8843–8859. [Google Scholar] [CrossRef]
  22. Anisimov, S.I.; Luk’Yanchuk, B.S. Selected problems of laser ablation theory. Physics-Uspekhi 2002, 45, 293–324. [Google Scholar] [CrossRef] [Green Version]
  23. Stapleton, M.W.; McKiernan, A.P.; Mosnier, J.-P. Expansion dynamics and equilibrium conditions in a laser ablation plume of lithium: Modeling and experiment. J. Appl. Phys. 2005, 97, 064904. [Google Scholar] [CrossRef]
  24. Amoruso, S.; Bruzzese, R.; Spinelli, N.; Velotta, R. Characterization of laser-ablation plasmas. J. Phys. B At. Mol. Opt. Phys. 1999, 32, R131–R172. [Google Scholar] [CrossRef]
  25. Neamtu, J.; Mihailescu, I.N.; Ristoscu, C.; Hermann, J. Theoretical modelling of phenomena in the pulsed-laser deposition process: Application to Ti targets ablation in low-pressure N2. J. Appl. Phys. 1999, 86, 6096–6106. [Google Scholar] [CrossRef]
  26. Chen, Z.; Bogaerts, A. Laser ablation of Cu and plume expansion into 1 atm ambient gas. J. Appl. Phys. 2005, 97, 063305. [Google Scholar] [CrossRef] [Green Version]
  27. Anisimov, S.I.; Bäuerle, D.; Luk’Yanchuk, B.S. Gas dynamics and film profiles in pulsed-laser deposition of materials. Phys. Rev. B 1993, 48, 12076–12081. [Google Scholar] [CrossRef]
  28. Zel’Dovich, Y.B.; Raizer, Y.P.; Hayes, W.D.; Probstein, R.F.; Gill, S.P. Physics of Shock Waves and High-Temperature Hydrodynamic Phenomena; Academic Press: New York, NY, USA, 1996. [Google Scholar] [CrossRef]
  29. Haynes, W.M.; Lide, D.R.; Bruno, T.J. CRC Handbook of Chemistry and Physics, 97th ed.; CRC Press: Boca Raton, FL, USA, 2017. [Google Scholar]
  30. Grant, P.M. Reflectivity of YIG And YGG: Observation of charge transfer and crystal field transitions. Appl. Phys. Lett. 1967, 11, 166–168. [Google Scholar] [CrossRef]
  31. Metalfron. Available online: https://www.metalfron.com/product/yttrium-iron-garnet-target-3/ (accessed on 1 September 2022).
  32. Nur-E-Alam, M.; Vasiliev, M.; Alameh, K. Optical constants of rare-earth substituted amorphous oxide-mix-based layers deposited to enable synthesis of magneto-optic garnets. Opt. Mater. 2019, 98, 109309. [Google Scholar] [CrossRef]
  33. Suwardy, J.; Darma, Y. Thermal annealing effects of SrTiO3 film on Si(100). In AIP Conference Proceedings; American Institute of Physics: College Park, MD, USA, 2013; Volume 1554, pp. 112–115. [Google Scholar] [CrossRef]
Figure 1. Schematic of laser irradiation of YIG target. The target surface is denoted by z = 0. The up and down arrows represent the laser energy that is reflected and absorbed by the target.
Figure 1. Schematic of laser irradiation of YIG target. The target surface is denoted by z = 0. The up and down arrows represent the laser energy that is reflected and absorbed by the target.
Micromachines 13 02012 g001
Figure 2. Schematic of the plasma plume expansion after termination of the laser pulse ( t = t L = 0 ). The initial plasma plume, with dimensions of height ( Z 0 ) and radius ( X 0 , Y 0 ) at ( t = 0 ), expands (X(t), Y(t), Z(t)) at t > 0) as a semi-ellipsoid perpendicular to the target surface (z-axis). The time when the edge of the plasma plume reached the substrate is denoted as t S .
Figure 2. Schematic of the plasma plume expansion after termination of the laser pulse ( t = t L = 0 ). The initial plasma plume, with dimensions of height ( Z 0 ) and radius ( X 0 , Y 0 ) at ( t = 0 ), expands (X(t), Y(t), Z(t)) at t > 0) as a semi-ellipsoid perpendicular to the target surface (z-axis). The time when the edge of the plasma plume reached the substrate is denoted as t S .
Micromachines 13 02012 g002
Figure 3. Effect of laser irradiance on the surface temperature of the YIG target. The left y-axis denotes the surface temperature ( T ), while the right y-axis denotes the intensity of laser, I ( t ) . T was calculated by numerically solving 1D heat conduction equations as expressed in the Equation (3) for temperatures under the melting point ( T     T m ) and Equation (8) for T above the melting point ( T > T m ). γ = 1.2 and TP = 7000 K, resulting in maximum height of plume H m a x = 13.29 µm, are used in the calculation. The numbers of partitions are M = 12,000 and N = 300.
Figure 3. Effect of laser irradiance on the surface temperature of the YIG target. The left y-axis denotes the surface temperature ( T ), while the right y-axis denotes the intensity of laser, I ( t ) . T was calculated by numerically solving 1D heat conduction equations as expressed in the Equation (3) for temperatures under the melting point ( T     T m ) and Equation (8) for T above the melting point ( T > T m ). γ = 1.2 and TP = 7000 K, resulting in maximum height of plume H m a x = 13.29 µm, are used in the calculation. The numbers of partitions are M = 12,000 and N = 300.
Micromachines 13 02012 g003
Figure 4. Temperature distribution inside the YIG target during laser pulse irradiation (t ≤ 15 ns) calculated for maximum depth z = 10 µm. Numerical parameters used as the same as the result shown in Figure 3.
Figure 4. Temperature distribution inside the YIG target during laser pulse irradiation (t ≤ 15 ns) calculated for maximum depth z = 10 µm. Numerical parameters used as the same as the result shown in Figure 3.
Micromachines 13 02012 g004
Figure 5. Plume expansion after the termination of laser for initial plasma plume temperature, T P = 20,000, 15,000 and 7000 K. (a) Edge velocity or center mass velocity, V z and (b) The height of plume edge in the z-direction, Z ( t ) .
Figure 5. Plume expansion after the termination of laser for initial plasma plume temperature, T P = 20,000, 15,000 and 7000 K. (a) Edge velocity or center mass velocity, V z and (b) The height of plume edge in the z-direction, Z ( t ) .
Micromachines 13 02012 g005
Figure 6. Velocity distribution inside the expanding plasma plume in different initial plasma temperatures, (a) TP = 20,000 K, (b) 15,000 K, and (c) 7000 K. Red-dashed line represent the semi-ellipsoid plasma plume dimension at height of 4 cm.
Figure 6. Velocity distribution inside the expanding plasma plume in different initial plasma temperatures, (a) TP = 20,000 K, (b) 15,000 K, and (c) 7000 K. Red-dashed line represent the semi-ellipsoid plasma plume dimension at height of 4 cm.
Micromachines 13 02012 g006
Table 1. Laser parameters used for simulation. Referring to parameters of Nd:YAG Continuum Surelite III model.
Table 1. Laser parameters used for simulation. Referring to parameters of Nd:YAG Continuum Surelite III model.
ParametersSymbolValue
Wavelength λ 532 nm
Laser Energy E 446 mJ
Repetition Rate-10 Hz
Beam Diameter D 9.5 mm
Pulsed Width (FWHM)w5 ns
Laser Fluence ϕ 0.63 J/cm2
Peak Intensity I m a x 1.26 × 108 W/cm2
Table 2. Thermal and optical properties of Yttrium-Iron Garnet (YIG) used for simulation.
Table 2. Thermal and optical properties of Yttrium-Iron Garnet (YIG) used for simulation.
PropertiesSymbolValueRefs.
Solid phase thermal conductivity K s 7.4 J/smK[29]
Liquid phase thermal conductivity K l 7.4 J/smK-
Solid phase density ρ s 5170 kg/m3[29]
Liquid phase density ρ l 5170 kg/m3-
Enthalpy of vaporization Δ H v 8 × 106 J/kg-
ReflectivityR0.13[30]
Solid phase specific heat C s 590 J/kgK[29]
Liquid phase specific heat C l 590 J/kgK-
Thermal diffusivity α 2.4 × 10−6 m2/scalculated
Melting point T m 1828 K[29]
Boiling point T b 3611 K[31]
Absorption coefficient (@532 nm) μ 2.5 × 106 m−1[32]
Table 3. Plasma plume parameters used in the calculation.
Table 3. Plasma plume parameters used in the calculation.
γTP (K)Z0 = Hmax (µm)E/M (J/kg)
1.2700013.29 7.9 × 10 5
1.215,00019.45 1.7 × 10 6
1.220,00022.46 2.2 × 10 6
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Suwardy, J.; Akhir, M.P.M.; Kurniawan, R.; Hermanto, B.R.; Anshori, I.; Fauzi, M.H. Simulation of Laser-Heating and Energetic Plasma Plume Expansion in Pulsed Laser Deposition of Y3Fe5O12. Micromachines 2022, 13, 2012. https://doi.org/10.3390/mi13112012

AMA Style

Suwardy J, Akhir MPM, Kurniawan R, Hermanto BR, Anshori I, Fauzi MH. Simulation of Laser-Heating and Energetic Plasma Plume Expansion in Pulsed Laser Deposition of Y3Fe5O12. Micromachines. 2022; 13(11):2012. https://doi.org/10.3390/mi13112012

Chicago/Turabian Style

Suwardy, Joko, Muzakkiy Putra Muhammad Akhir, Robi Kurniawan, Beni Rio Hermanto, Isa Anshori, and Mohammad Hamzah Fauzi. 2022. "Simulation of Laser-Heating and Energetic Plasma Plume Expansion in Pulsed Laser Deposition of Y3Fe5O12" Micromachines 13, no. 11: 2012. https://doi.org/10.3390/mi13112012

APA Style

Suwardy, J., Akhir, M. P. M., Kurniawan, R., Hermanto, B. R., Anshori, I., & Fauzi, M. H. (2022). Simulation of Laser-Heating and Energetic Plasma Plume Expansion in Pulsed Laser Deposition of Y3Fe5O12. Micromachines, 13(11), 2012. https://doi.org/10.3390/mi13112012

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop