Next Article in Journal
Numerical Study on Surface Roughness Measurement Based on Nonlinear Ultrasonics in Through-Transmission and Pulse-Echo Modes
Next Article in Special Issue
Acid Solution Processed VO2-Based Composite Films with Enhanced Thermochromic Properties for Smart Windows
Previous Article in Journal
Impact of Polyethylene Glycol Functionalization of Graphene Oxide on Anticoagulation and Haemolytic Properties of Human Blood
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Towards Scalable Large-Area Pulsed Laser Deposition

1
Federal Research Centre, Southern Scientific Centre of the Russian Academy of Sciences (SSC RAS), 41 Chekhov St., 344006 Rostov-on-Don, Russia
2
Research Laboratory of Functional Nanomaterials Technology, Southern Federal University, 2 Shevchenko St., 347922 Taganrog, Russia
3
Department of Electrical & Computer Engineering, National University of Singapore, 4 Engineering Drive 3, Singapore 117583, Singapore
4
FORS Development Center, 3, Trifonovskiy Tupik, 129272 Moscow, Russia
5
Institute of Nanotechnologies, Electronics and Equipment Engineering, Southern Federal University, 2 Shevchenko St., 347922 Taganrog, Russia
6
Research and Education Centre ‘Nanotechnologies’, Southern Federal University, 2 Shevchenko St., 347922 Taganrog, Russia
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Materials 2021, 14(17), 4854; https://doi.org/10.3390/ma14174854
Submission received: 22 July 2021 / Revised: 12 August 2021 / Accepted: 20 August 2021 / Published: 26 August 2021
(This article belongs to the Special Issue Preparation and Properties of Thin Films)

Abstract

:
One of the significant limitations of the pulsed laser deposition method in the mass-production-technologies of micro- and nanoelectronic and molecular device electronic fabrication is the issue of ensuring deposition of films with uniform thickness on substrates with large diameter (more than 100 mm) since the area of the laser spot (1–5 mm2) on the surface of the ablated target is incommensurably smaller than the substrate area. This paper reports the methodology that allows to calculate the distribution profile of the film thickness over the surface substrate with a large diameter, taking into account the construction and technological parameters of the pulsed laser deposition equipment. Experimental verification of the proposed methodology showed that the discrepancy with the experiment does not exceed 8%. The modeling of various technological parameters influence on the thickness uniformity has been carried out. Based on the modeling results, recommendations and parameters are proposed for manufacturing uniform thickness films. The results allow for increasing the film thickness uniformity with the thickness distribution < 5% accounts for ~ 31% of 300 mm diameter substrate.

1. Introduction

Currently, the pulsed laser deposition (PLD) method is widely used to form epitaxial and single-crystal complex oxides films with ferroelectric, ferromagnetic, dielectric, and superconducting properties [1,2,3,4,5]. The advantages of this method include the possibility of maintaining the stoichiometric composition of the ablated material [6], good film adhesion [7], versatility in choosing the deposited material, as well as the possibility of forming film coatings on the surface of thermosensitive materials [8]. These advantages allow to use of the PLD method in the formation of energy harvesting devices [9], optical structures [10], sensor elements [11], ferroelectric films [12], memristor structures [13,14], and medical coatings [15]. Thus, the study of the PLD process and the development of methods for solving issues limiting its implementation in mass production is highly urgent.
In the PLD method, a laser beam is focused on a target placed in a vacuum chamber. The target material is ablated under the influence of laser radiation of high-power density, as the result the ablated particles are deposited on the substrate. Although the PLD method is widely used in scientific research [16,17,18,19,20,21,22], its industrial application is limited [1,23]. One of the significant issues of PLD in mass-production-technologies of micro- and nanoelectronic device fabrication is ensuring uniform thickness films deposition on large diameter substrates (more than 100 mm), since the area of the laser spot (1–5 mm2) on the surface of the ablated target is incommensurably smaller than the substrate area [24].
In the modern micro- and nanoelectronic industry, the standard is the non-uniformity of the film thickness over the surface of the entire substrate ~5%, except for 5 mm at the edge. Since the PLD method is relatively new [25,26] and its integration with industrial technologies of micro- and nanoelectronics is still at its initial stages. The study of regularities and the search for methods to achieve the indicated values of non-uniformity on substrates of large diameter (more than 100 mm) requires additional research. In [27,28,29,30,31,32,33,34], several approaches for forming films with a uniformity of 70–95% on substrates with a diameter of 100–200 mm are described in detail. However, obtaining uniform of films by PLD on 300 mm-diameter substrates is still challenging and often requires using scanning systems that ensure the movement of the laser beam over the target surface and the movement of the substrate [1]. One of the possible ways to overcome this limitation is to use a laser beam scanning system to move the laser beam over the target surface, rotate the target and substrate, and optimize the parameters of these processes, considering the geometry of the growth chamber.
The purpose of this work is to study the processes affecting film thickness uniformity deposited by the PLD method on substrates of large diameter (100 mm and more) and developing the methodology that allows calculating the film thickness distribution profile over the substrate surface, taking into account the design and technological parameters of PLD equipment. The results presented in this work are partially based on the previous studies, conducted both in our laboratory as well as by other researchers [27,28,29,30,31,32,33,34]. We implemented a new approach to the calculation of a trajectory of a laser beam along the target surface (Video S1 and S2). These improvements made it possible to describe a trajectory of a laser beam along the target surface much more correctly and to take into account the influence of a larger number of scanning parameters. As a result, we obtained quantitatively more accurate results of the film thickness spatial distribution.

2. Materials and Methods

Modeling the PLD process, assessing technological parameters and their influence on the film thickness uniformity were carried out in the MATLAB software (MathWorks Inc., Natick, MA, USA). The process of scanning a target with a diameter of up to 50 mm by a laser beam is considered on the example of cluster nanotechnological complex NANOFAB NTK-9 (NT-MDT, Zelenograd, Russia), comprising of a Pioneer 180 PLD module (Neocera LCC, Beltsville, MD, USA). The parameters used in the modeling are presented in Table 1. The studied parameter took the values indicated in the table, while all the others remained unchanged (highlighted in Table 1), to estimate the influence of one of the scanning system parameters on the films thickness uniformity.
The movement of the laser beam over the target surface is carried out using a scanning system consisting of a pair of mirrors and a focusing lens fixed on a movable stage. The laser beam hits the mirrors and re-reflected on the focusing lens. Then the laser beam hits the target surface through the window in the growth chamber. The movable stage is part of the module’s scanning system, moves along the x-axis, thereby scanning the laser beam along the target surface (Figure 1).
The Upper Limit and Lower Limit parameters are the lower and upper limits of the laser beam movement and allow localizing the area of action of the laser beam on the target surface by a segment (Lower Limit, Upper Limit). The point in this segment at which the irradiation process starts is determined by the Origin parameter. Values of all of the aforementioned parameters lie within [0, 2R], where R—target radius. The Max Velocity and Min Velocity parameters determine the maximum and minimum speed of movement of the laser beam on the target surface. Near the point defined by the Origin parameter, the scanning speed varies within (Min Velocity, Max Velocity). On the (Origin, Upper Limit) section, the change in scanning speed is determined by the Upper Coefficient parameter, and on the (Lower Limit, Origin) section by the Lower Coefficient parameter.
The speed of scanning system movement is described by two functions that intersect at a point with a coordinate defined by the Origin parameter. The relationship between the functions is described by the expression:
V a ( x ) = V max R u C υ ( x O r i g i n ) ,
where V max —maximum speed of the scanning system, R u is defined as the largest of two values |Lower Limit–Origin| and |Upper Limit–Origin|, and the coefficient C υ is equal to Lower Coefficient for all x < Origin and Upper Coefficient for all x > Origin.
The mutual arrangement of the centers of the target and the substrate is not coaxial; therefore, to form a continuous film, the substrate rotates with a certain angular velocity, since the scanning system ensures the movement of the laser beam only along the diameter of the target. The target also rotates with an angular velocity to ensure uniform irradiation. The motion control of the scanning system is carried out using the Pioneer 180 PLD software (Neocera LCC, Beltsville, MD, USA).
It is necessary to describe the spatial distribution of molecules in directions during target ablation to estimate the uniformity of film deposition on a substrate. Since the laser spot has a small area, and the crater on the target surface after exposure to the laser pulse has a depth of about several nanometers [35,36], the volume of the material removed in one pulse is also relatively small. Thus, the area of laser action can be considered as a Knudsen cell, consisting of an isothermal shell with an infinitely small hole d A e and infinitely thin walls (supplementary materials Figure S1) [37].
The ablation region contains N atoms that collide with the cell walls and are reflected from them without changing the velocity. Atoms moving towards the hole leave the cell at the same speed. The expression determines the distribution of the velocities of atoms in the flow of matter for a small number of atoms within a small spatial angle d ω s . The angle φ determines the direction of movement with respect to the normal of the hole d A e . It is possible to obtain the mass of the substance deposited per unit area taking into account the distribution of the velocities of atoms in directions [36]:
d M r ( φ , θ ) d A r = M e π r 2 cos φ cos θ ,
where M e —mass of vaporized matter.
From Equation (2), it follows that the propagation of matter occurs mainly in directions close to the normal to the evaporated surface (cosφ → max). Equation (2) can be simplified, since in the PLD module, the substrate is located parallel to the target, the angles φ and θ are equal, and the cosines of these angles are equal to h/l, and the distribution of the film thickness on the substrate is described by the expression [36]:
d = M e π ρ u h 2 [ 1 + ( l / h ) 2 ] 2 ,
where l—the current coordinate along the substrate radius, h—target-substrate distance [38], mass of ablated substance M e = ρ u V a ; V a = S l a s . s p o t Δ z , Δ z —target material ablation depth, S l a s . s p o t —laser spot area, and ρ u —target material density.
It is convenient to characterize the uniformity of film deposition over the diameter of the substrate using the ratio:
d d max = d ( l ) d ( 0 ) = [ 1 + ( l / h ) 2 ] 2 ,
where d max —maximum film thickness. The value of the relative thickness allows one to compare films obtained at different values of technological parameters, regardless of the absolute values of the film thickness.
We introduce a parameter Ω 95 that is defined as:
Ω 95 = S 95 S ,
where S95—the substrate area where relative thickness lies within [0.95, 1] and S—area of the entire substrate. If the relative thickness allows one to conclude the uniformity of the film, then the value of the Ω 95 parameter makes it possible to characterize the distribution of the relative film thickness from the point of view of industrial application (Figure S2). In particular, using the Ω 95 parameter, it is possible to estimate the size of the substrate for the deposition of the films with thickness uniformity is 5% by the PLD method under specific scanning modes.
The methodology considered above does not take into account the size of the laser spot on the target surface and the associated with “flip-over effect” [37]. In [39], the equation describes the thickness of a film on a substrate, taking into account the gas dynamics of three-dimensional vapor expansion during PLD and the associated evolution of the shape of the laser plume was proposed:
d = d max ( 1 + 1 p t g 2 θ x + q 2 p t g 2 θ y ) 3 / 2 ,
where p and q—coefficients are obtained from the solution of the system of gas-dynamic equations and θ x and θ y —the corresponding components of the particle scattering angle. However, solving the system of equations for finding the coefficients p and q significantly complicates the calculations.
The film growth process during PLD depends on many different factors [40,41,42,43,44,45,46]. To simplify the methodology, the following assumptions have been made: laser ablation occurs in a vacuum; the target surface is considered ideally smooth under the PLD process (there is no laser modification of the target surface); there is no desorption of the deposited material from the substrate surface; the shape of the laser plume in the process of moving towards the substrate remains unchanged.
Figure 2 shows the calculation of the films thickness algorithm based on the developed methodology (1)–(6).
In the block for input of modeling parameters, technological parameters that determine the modes of scanning the laser beam over the target surface, as well as the rotation speed of the target and substrate are set.
In the block for processing of modeling parameters, a check for the consistency of the given parameters and their transformation to a single frame of reference is made.
In the block for calculating of laser trajectory on target surface, the coordinates corresponding to each individual pulse are calculated, considering the operating modes of the scanning system and the parameters of the roll and rotation of the target. Figure 3 shows an example of the calculated trajectories of the laser beam along the target surface for various simulation parameters. Each point corresponds to one laser pulse, with the pulse repetition rate of 10 Hz, and the number of pulses is 6000.
The block for constructing the projection of the laser trajectory on the target surface onto the substrate is necessary to translate the array of coordinates of laser pulses from the reporting system associated with the center of the target to the reporting system associated with the center of the substrate, taking into account the parameters of the substrate rotation. It is necessary to calculate the distances r and l (Figure S1).
The film thickness calculation block specifies a set of points in a polar coordinate system. In the radial direction, the grid step is d r = r 2 r 1 = r 3 r 2 = = r n r n 1 , and in the corner d φ = φ 2 φ 1 = φ 3 φ 2 = = φ m φ m 1 . The resulting set of points can be represented as a matrix:
[ 0 r 1 r 2 r n 1 r n φ 1 d 1 , 1 d 1 , 2 d 1 , n 1 d 1 , n φ 2 d 2 , 1 d 2 , 2 d 2 , n 1 d 2 , n φ m 1 d m 1 , 1 d m 1 , 2 d m 1 , 3 d m 1 , n 1 d m 1 , n φ m d m , 1 d m , 2 d m , 3 d m , n 1 d m , n ]
where d m n —thickness value at point ( r n , φ m ).
Then Equation (3) for determining the film thickness at each point with coordinates ( r n , φ m ) on the substrate surface has the form
d m , n = M e π ρ u h 2 i = 1 m j = 1 n [ 1 + ( l i , j / h ) 2 ] 2 .
According to the Equation (8), the film thickness d is determined, and the obtained value is written into the corresponding cell of the matrix. After receiving an array of thicknesses, the search for the most considerable value is performed. Using the maximum film thickness d max , one can calculate the relative thickness d = d n / d max , which allows comparing the simulation results obtained at different values of technological parameters.
In the output modeling results block, the results are saved in text and graphic formats.
Based on the described methodology, we studied the influence of the scanning system modes on the film thickness distribution profile over a substrate with a diameter of 100 mm and 300 mm.
Experimental studies were carried to prove the proposed methodology. Nanocrystalline LiNbO3 films were obtained on a silicon substrate with a diameter of 100 mm using the Pioneer 180 PLD module of the NANOFAB NTK-9 cluster nanotechnological complex. An excimer KrF laser (λ = 248 nm) was used to ablate a LiNbO3 target with a purity of 99.9% (Kurt J. Lesker Company Ltd., East Sussex, UK). Energy density on the target surface is maintained at 1.5 J/cm2. The target-substrate distance (75 mm), number of pulses (36,000), pulse repetition rate (10 Hz), and laser pulses energy on the target surface (150 mJ) are kept constant. The morphology of the obtained films was studied by atomic force microscopy (AFM) using the Ntegra Probe Nano Laboratory (NT-MDT, Zelenograd, Russia). The experimental data were processed using the Image Analysis 3.5 software (NT-MDT, Zelenograd, Russia).
The study of the thickness of the LiNbO3 films was carried out by measuring the LiNbO3/Si structure height obtained by plasma-chemical etching. For this purpose, an FP-383 photoresist was applied to the surface of the LiNbO3 films by spinning at a rotation speed of 3000 rpm. After that, the photoresist film was pre-cured in the air for 10 min in order to prevent defect formation. Then, the sample was cured in an oven at 90 °C for 30 min. The sample was exposed to the UV radiation for 2 min through a photomask and then developed in a 5% aqueous solution of KOH, and hardbacked at 110 °C for 25 min on a hotplate. The obtained structure was processed in a module for plasma-chemical etching in a combined plasma of capacitive and inductive STE ICPe68 (SemiTEq St. Petersburg, Russia). Etching took place at a pressure of 2 Pa, a capacitive plasma source power of 35 W, an inductively coupled plasma source power of 400 W, a bias voltage of 75 V, a fluorinated SF6 gas flow of 15 sccm, and the etching time was 1 min. Photoresist residues were removed in dimethylformamide. In total, ten samples were used to fabricate structures in LiNbO3 films, the height of which was investigated by the Ntegra Laboratory in the semicontact AFM mode using NSG11 cantilevers.

3. Results and Discussion

Figure 4 shows the results of AFM studies of LiNbO3 film. Based on the obtained experimental results, the dependence of the distribution of the thickness of the LiNbO3 film over the surface of the silicon substrate is plotted (Figure 5).
To determine the LiNbO3 films thickness we used the approach described in detail in [47].
The obtained experimental results are in good agreement with the calculations based on the proposed Equations (1)–(8) for a substrate with a diameter of 100 mm. The thickness of the obtained LiNbO3 films increases from (67.2 ± 5.1) nm, reaching a maximum ((90.4 ± 7.9) nm) in the center of the substrate. The calculated thickness of the films varies from 63.6 nm to 89.1 nm, respectively. The change in the relative film thickness is less than 37% for the experiment and ~29% for the simulation results. The deviation of the simulation results from the experimental data does not exceed 8%. The parameter Ω 95 is ~15%. The region with the thickness difference is 10% of the maximum, accounts for ~31% of the substrate surface. The discrepancy between the simulation and experiment may be associated with an uneven erosion of the target surface over multiple deposition cycles and a tilt of an ablation plume, as it is described in [23]. Thus, experimental studies confirm that the proposed methodology allows to calculate the film thickness’s value obtained by the PLD method.
Figure 6 shows the relative thickness of the LiNbO3 films on the Upper Limit and Lower Limit parameters, calculated based on the developed Equations (1)–(8). Figure 6 shows the results of the study of the influence of scanning parameters on the uniformity of the thickness of the ZnO films. Insets in Figure 6, Figure 7, Figure 8 and Figure 9 provide a comparison of relative thicknesses for different cases at a higher scale. The appearance and change in the size of the “plateau” in the upper part of the dependence illustrate the physical meaning of the Ω 95 parameter.
As the Lower Limit increases with the Upper Limit unchanged, the Ω 95 parameter value increases from 9% to ~31%, while when the Upper Limit decreases with the Lower Limit unchanged, it decreases from ~8% to 4%. It can be associate with the fact that in the first case, we move the projection of the area of scanning of the target surface with a laser beam from the center of the substrate to the edge, and in the second, we bring it closer (Figure S1).
The results of studying the effect of the Origin parameter (Figure 7a) on the relative thickness of the films showed when the value of Origin decreases from 50 to 0, the value of Ω 95 the parameter increases from 6% to ~20%. This effect is linked to the distance of the projection of the target area with laser beam moves at a speed of Min Velocity from the center of the substrate to the edge (Figure S1). Thus, the Origin parameter affects the relative film thickness by moving the projections of the regions with the maximum and minimum scanning speeds of the target surface (parameters Min Velocity and Max Velocity) relative to the center of the substrate. Figure 7b shows the influence of Lower Limit, Upper Limit, and Origin on Ω 95 parameter and substrate radius.
A change in the maximum and minimum speeds of the laser beam movement over the target surface has practically no effect on the distribution of the relative thickness of the resulting films (Figure 8). Since the Min Velocity and Max Velocity parameters do not affect the projections of areas with high and low scanning speeds of the target surface relative to the center of the substrate, but only determines the scanning speed.
Unlike Max Velocity and Min Velocity, the Upper Coefficient and Lower Coefficient parameters have a more significant effect on the distribution of the relative thickness of the films (Figure 8). It is linked to the possibility of influencing the smoothness of the transition from Min Velocity to Max Velocity in the (Lower Limit, Origin) and Max Velocity in Min Velocity on the (Origin, Upper Limit) segment and, as a consequence, change the size of the areas that have the maximum and minimum scanning speeds. Therefore, with the values of Upper Coefficient and Lower Coefficient close to unity, the scanning speed practically does not change when switching from Max Velocity to Min Velocity, while at values above 60, the speed changes abruptly. Thus, increasing the Upper Coefficient value in the (Origin, Upper Limit) segment increases the area with the scanning speed is equal to the Min Velocity, which leads to the displacement of the projection of the target area with the minimum scanning speed from the center to the edge of the substrate. Figure 9c shows the Ω 95 and substrate diameter dependencies on the value of the parameters Upper Coefficient and Lower Coefficient. As the Upper Coefficient parameter increases with a constant Lower Coefficient, Ω 95 increases from ~5% to ~20%, and with an increase in the Lower Coefficient parameter, with an unchanged Upper Coefficient  Ω 95 decreases from ~24% to ~5%.
Based on the results of the performed calculations, recommendations are proposed that allow increasing the thickness uniformity of the films formed by the PLD method.
It has been established that increasing the operating mode of the scanning system with the projections of the target area with a low scanning speed will be shifted from the center to the edge of the substrate. This effect can be achieved in several ways:
changing the scanning range—Upper Limit and Lower Limit parameters;
changing the coordinate in which the scanning speed reaches its maximum—Origin parameter;
varying the smoothness of changes in scanning speeds—Upper Coefficient и Lower Coefficient parameters.
It should be noted that the first method should be used to a limited extent since the Upper Limit and Lower Limit parameters strongly affect the uniformity of the target erosion.
Using the calculation results, the operating modes of the scanning system with a high Ω 95 value when films are deposited on substrates with a diameter of 300 mm were determined: Upper Limit = 50, Lower Limit = 10, Origin = 30, Max Velocity = 5, Min Velocity = 0.1, Upper Coefficient = 60, and Lower Coefficient = 1. For these parameters, the projection onto the surface of the substrate of the laser beam trajectory along the target surface and the distribution of the relative thickness of the LiNbO3 film on the substrate were calculated (Figure 10).
The projection of the laser beam trajectory onto the surface of the substrate (Figure 10a) makes it possible to estimate the location where ablation proceeds relative to the center of the substrate. Figure 10b shows the spatial distribution of the relative thickness of the obtained film over the surface of a substrate with a diameter of 300 mm. It was found that the change in the relative thickness of the film does not exceed 50%, but Ω 95 is ~31%, while the region with the difference in thickness is 10%, of which the maximum is ~38%.

4. Conclusions

The methodology development results allow calculating the film thickness distribution profile over the surface of the substrate, taking into account the design and technical parameters of the PLD equipment. The film thickness distribution profile on a substrate with a diameter of 100 mm is calculated considering the design features and operating values of the technological parameters of the Pioneer 180 PLD module of cluster nanotechnological complex NANOFAB NTK-9. The thickness of the obtained LiNbO3 films increases from (67.2 ± 5.1) nm, reaching a maximum ((90.4 ± 7.8) nm) in the center of the substrate. The calculated thickness of the films varies from 63.6 nm to 89.1 nm, respectively. The change in the relative film thickness is less than 37% for the experiment and ~ 29% for the simulation results. A good correlation of the results of theoretical calculations with experimental data is shown: the discrepancy between the model and experiment does not exceed 8%. Moreover, the developed methodology allows to solve the inverse task and scale it, overcoming relevant issues in modern micro- and nanoelectronic technology. It is possible to determine the technological parameters of the PLD module for obtaining a film with a controlled uniformity of thickness distribution using the modeling results. The analysis of the influence of the operating parameters of the scanning system (Table 1) on the uniformity of obtaining films on a substrate with a diameter of 300 mm is carried out. Recommendations for obtain films with controlled irregularity by PLD method are proposed: (Upper Limit = 50, Lower Limit = 10, Origin = 30, Max Velocity = 5, Min Velocity = 0.1, Upper Coefficient = 60, and Lower Coefficient = 1). Based on obtained results, the operating modes of the scanning system are proposed, which allows obtaining films with a relative thickness difference not exceeding 50%. In this case, the region in which the difference in thickness is less than 5% ( Ω 95 ) accounts for 31% of the substrate surface, and the region in which the difference in thickness is 10% of the maximum is ~38%.
The obtained results can be used to expand the capabilities of the large-area PLD technology, which will speed up combining the technology of multicomponent oxide films laser deposition with the silicon technology of micro- and nanoelectronics and molecular device electronic to create MEMS and new generation of energy conversion devices and sensors. Moreover, the developed methodology might be modified and scaled to calculate the thickness profile of different materials films obtained using various PLD equipment (Solmates B.V., PVD Products Inc, and Neocera LCC).

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/ma14174854/s1, Figure S1: Target material evaporation from a cell d A e by a laser beam on substrate surface element d A r , (a) and Pioneer 180 PLD software (b), Figure S2: Example of practical use of the obtained dependencies, Video S1: Scanning mode by the target with a stationary laser beam, Video S2: Scanning mode of the laser beam along the surface of a rotating target.

Author Contributions

Conceptualization and validation O.A.A., A.M. and Z.V.; methodology, A.M., D.K., V.A.S. and Z.V.; investigation, Z.V., D.K. and E.Z.; writing—Original draft preparation, Z.V., D.K. and O.A.A.; data curation and formal analysis, V.A.S., V.S.K., D.K. and R.T.; visualization R.T., D.K., E.Z. and Z.V.; lithium niobate thin films growth Z.V.; supervision, O.A.A.; project administration, writing—Review and editing, O.A.A.; plasma etching of films, V.S.K.; AFM study, R.T. and V.A.S. All authors contributed to the writing of the manuscript, which was coordinated by O.A.A. and Z.V., Z.V. and D.K. contributed equally to the work. All authors have read and agreed to the published version of the manuscript.

Funding

The reported study was funded by the Russian Foundation for Basic Research, projects No. [19-38-60052] (Development of a methodology and study of the regularities of the influence of pulsed laser deposition modes on the uniformity of material distribution during the formation of ferroelectric films), No. [19-29-03041 mk] (Development of a methodology for experimental studies of the effect of film formation modes by PLD), No. [19-37-90139] (Experimental studies on the formation of LiNbO3 films), and Russian Scientific Foundation, project No. [20-69-46076] (Development of methods for plasma etching and control of the parameters of the fabricated structures) partially supported by Grant of the President of the Russian Federation No. [МК-6252.2021.4] (Development of the computational block of the program for calculating the trajectory of the laser beam along the target surface).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing is not applicable to this article.

Acknowledgments

The work was conducted on the equipment of the Research and Education Centre “Nanotechnology”, Southern Federal University.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Blank, D.H.A.; Dekkers, M.; Rijnders, G. Pulsed laser deposition in Twente: From research tool towards industrial deposition. J. Phys. D Appl. Phys. 2013, 47, 034006. [Google Scholar] [CrossRef]
  2. Ahn, Y.; Son, J.Y. Growth and ferroelectric properties of highly c-oriented epitaxial Aurivillius BaBi4Ti4O15 thin films. Mater. Sci. Eng. C 2021, 266, 115046. [Google Scholar] [CrossRef]
  3. Li, T.; Dong, J.; Zhang, N.; Wen, Z.; Sun, Z.; Hai, Y.; Wang, K.; Liu, H.; Tamura, N.; Mi, S.; et al. Interface control of tetragonal ferroelectric phase in ultrathin Si-doped HfO2 epitaxial films. Acta Mater. 2021, 207, 116696. [Google Scholar] [CrossRef]
  4. Khan, M.A.; Braic, L.; AlSalik, Y.; Idriss, H. Growth of epitaxial strontium titanate films on germanium substrates using pulsed laser deposition. Appl. Surf. Sci. 2021, 542, 148601. [Google Scholar] [CrossRef]
  5. Najim, A.A. Synthesis and characterizations of (δ-Bi2O3)0.93(TiO2)0.07 thin films grown by PLD technique for optoelectronics. Mater. Sci. Semicond. Process. 2017, 71, 378–381. [Google Scholar] [CrossRef]
  6. Liu, G.Z.; Lei, Q.Y.; Xi, X.X. Stoichiometry of SrTiO3 films grown by pulsed laser deposition. Appl. Phys. Lett. 2012, 100, 202902. [Google Scholar] [CrossRef]
  7. Narayan, R.J. Pulsed laser deposition of functionally gradient diamondlike carbon–metal nanocomposites. Diam. Relat. Mater. 2005, 14, 1319–1330. [Google Scholar] [CrossRef]
  8. Socol, G.; Socol, M.; Stefan, N.; Axente, E.; Popescu-Pelin, G.; Craciun, D.; Duta, L.; Mihailescu, C.N.; Mihailescu, I.N.; Stanculescu, A.; et al. Pulsed laser deposition of transparent conductive oxide thin films on flexible substrates. Appl. Surf. Sci. 2012, 260, 42–46. [Google Scholar] [CrossRef]
  9. Fourmont, P.; Gerlein, L.F.; Fortier, F.X.; Cloutier, S.G.; Nechache, R. Highly efficient thermoelectric microgenerators using nearly room temperature pulsed laser deposition. ACS Appl. Mater. Interfaces 2018, 10, 10194–10201. [Google Scholar] [CrossRef]
  10. Karnati, P.; Haque, A.; Taufique, M.F.N.; Ghosh, K. A systematic study on the structural and optical properties of vertically aligned zinc oxide nanorods grown by high pressure assisted pulsed laser deposition technique. Nanomaterials 2018, 8, 62. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Vakulov, Z.E.; Zamburg, E.G.; Khakhulin, D.A.; Ageev, O.A. Thermal stability of ZnO thin films fabricated by pulsed laser deposition. Mater. Sci. Semicond. Process. 2017, 66, 21–25. [Google Scholar] [CrossRef]
  12. Ji, X.; Wang, C.; Harumoto, T.; Tian, Y.; Zhang, S.; Tu, R.; Shen, Q.; Shi, J. Deposition-temperature dependence of structure, ferroelectric property and conduction mechanism of BCZT epitaxial films. Ceram. Int. 2021, 47, 3195–3200. [Google Scholar] [CrossRef]
  13. Tominov, R.V.; Vakulov, Z.E.; Avilov, V.I.; Khakhulin, D.A.; Fedotov, A.A.; Zamburg, E.G.; Smirnov, V.A.; Ageev, O.A. Synthesis and Memristor Effect of a Forming-Free ZnO Nanocrystalline Films. Nanomaterials 2020, 10, 1007. [Google Scholar] [CrossRef] [PubMed]
  14. Tominov, R.V.; Vakulov, Z.E.; Avilov, V.I.; Khakhulin, D.A.; Polupanov, N.V.; Smirnov, V.A.; Ageev, O.A. The Effect of Growth Parameters on Electrophysical and Memristive Properties of Vanadium Oxide Thin Films. Molecules 2021, 26, 118. [Google Scholar] [CrossRef] [PubMed]
  15. Al-Kattan, A.; Grojo, D.; Drouet, C.; Mouskeftaras, A.; Delaporte, P.; Casanova, A.; Robin, J.D.; Magdinier, F.; Alloncle, P.; Constantineskus, C.; et al. Short-Pulse Lasers: A Versatile Tool in Creating Novel Nano-/Micro-Structures and Compositional Analysis for Healthcare and Wellbeing Challenges. Nanomaterials 2021, 11, 712. [Google Scholar] [CrossRef] [PubMed]
  16. Koster, G.; Blank, D.H.; Rijnders, G.A. Oxygen in Complex Oxide Thin Films Grown by Pulsed Laser Deposition: A Perspective. J. Supercond. Nov. Magn. 2020, 33, 205–212. [Google Scholar] [CrossRef] [Green Version]
  17. Tominov, R.V.; Smirnov, V.A.; Avilov, V.I.; Fedotov, A.A.; Klimin, V.S.; Chernenko, N.E. Formation of ZnO memristor structures by scratching probe nanolithography. IOP Conf. Ser. Mater. Sci. Eng. 2018, 443, 012036. [Google Scholar] [CrossRef] [Green Version]
  18. Fernández-Arias, M.; Boutinguiza, M.; Del Val, J.; Riveiro, A.; Rodríguez, D.; Arias-González, F.; Gil, J.; Pou, J. Fabrication and deposition of copper and copper oxide nanoparticles by laser ablation in open air. Nanomaterials 2020, 10, 300. [Google Scholar] [CrossRef] [Green Version]
  19. Ageev, O.A.; Gusev, E.Y.; Zamburg, E.G.; Vakulov, D.E.; Vakulov, Z.E.; Shumov, A.V.; Ivonin, M.N. Nanocrystalline ZnO films grown by PLD for NO2 and NH3 sensor. Appl. Mech. Mater. 2014, 475, 446–450. [Google Scholar]
  20. Karapetyan, G.Y.; Kaydashev, V.E.; Zhilin, D.A.; Kutepov, M.E.; Minasyan, T.A.; Kaidashev, E.M. A surface acoustic wave impedance-loaded high sensitivity sensor with wide dynamic range for ultraviolet light detection. Sens. Actuator A Phys. 2019, 296, 70–78. [Google Scholar] [CrossRef]
  21. Nguyen, M.D.; Houwman, E.P.; Yuan, H.; Wylie-van Eerd, B.J.; Dekkers, M.; Koster, G.; tel Elshof, J.E.; Rijnders, G. Controlling piezoelectric responses in Pb(Zr0.52Ti0.48)O3 films through deposition conditions and nanosheet buffer layers on glass. ACS Appl. Mater. Interfaces 2017, 9, 35947–35957. [Google Scholar] [CrossRef] [PubMed]
  22. Smirnov, V.A.; Tominov, R.V.; Avilov, V.I.; Alyabieva, N.I.; Vakulov, Z.E.; Zamburg, E.G.; Ageev, O.A. Investigation into the Memristor Effect in Nanocrystalline ZnO Films. Semiconductors 2019, 53, 72–77. [Google Scholar] [CrossRef]
  23. Greer, J.A.; Tabat, M.D. Large-area pulsed laser deposition: Techniques and applications. J. Vac. Sci. Technol. A 1995, 13, 1175–1181. [Google Scholar] [CrossRef]
  24. Greer, J.A.; Tabat, M.D.; Lu, C. Future trends for large-area pulsed laser deposition. Nucl. Instrum. Methods Phys. Res. 1997, 121, 357–362. [Google Scholar] [CrossRef]
  25. Smith, H.M.; Turner, A.F. Vacuum deposited thin films using a ruby laser. Appl. Opt. 1965, 1, 147–148. [Google Scholar] [CrossRef]
  26. Venkatesan, T.; Wu, X.; Inam, A.; Chang, C.C.; Hegde, M.S.; Dutta, B. Laser processing of high-T/sub c/superconducting thin films. IEEE J. Quantum. Electron. 1989, 25, 2388–2393. [Google Scholar] [CrossRef]
  27. Pryds, N.; Schou, J.; Linderoth, S. The spatial thickness distribution of metal films produced by large area pulsed laser deposition. Appl. Surf. Sci. 2007, 253, 8231–8234. [Google Scholar] [CrossRef]
  28. Marcu, A.; Grigoriu, C.; Jiang, W.; Yatsui, K. Pulsed laser deposition of YBCO thin films in a shadow mask configuration. Thin Solid Films 2000, 360, 166–172. [Google Scholar] [CrossRef]
  29. Lackner, J.M. Industrially-scaled large-area and high-rate tribological coating by pulsed laser deposition. Surf. Coat. Technol. 2005, 200, 1439–1444. [Google Scholar] [CrossRef]
  30. Greer, J.A. History and current status of commercial pulsed laser deposition equipment. J. Phys. D Appl. Phys. 2013, 47, 034005. [Google Scholar] [CrossRef] [Green Version]
  31. Yang, H.; Wang, W.; Liu, Z.; Li, G. Epitaxial growth of 2 inch diameter homogeneous AlN single-crystalline films by pulsed laser deposition. J. Phys. D Appl. Phys. 2013, 46, 105101. [Google Scholar] [CrossRef]
  32. Eason, R.W.; May-Smith, T.C.; Grivas, C.; Darby, M.S.B.; Shepherd, D.P.; Gazia, R. Current state-of-the-art of pulsed laser deposition of optical waveguide structures: Existing capabilities and future trends. Appl. Surf. Sci. 2009, 255, 5199–5205. [Google Scholar] [CrossRef]
  33. Christen, H.M.; Eres, G. Recent advances in pulsed-laser deposition of complex oxides. J. Phys. Condens. Matter. 2008, 20, 264005. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Ageev, O.; Konoplev, B. (Eds.) Nanotechnology in Microelectronics; Nauka Publisher: Moscow, Russia, 2019; p. 511. [Google Scholar]
  35. Inogamov, N.A.; Zhakhovskii, V.V.; Ashitkov, S.I.; Petrov, Y.V.; Agranat, M.B.; Anisimov, S.I.; Nishihara, K.; Fortov, V.E. Nanospallation induced by an ultrashort laser pulse. J. Exp. Theor. Phys. 2008, 107, 1–19. [Google Scholar] [CrossRef]
  36. Agranat, M.B.; Anisimov, S.I.; Ashitkov, S.I.; Zhakhovskii, V.V.; Inogamov, N.A.; Nishihara, K.; Petrov, Y.V.; Fortov, V.E.; Khokhlov, V.A. Dynamics of plume and crater formation after action of femtosecond laser pulse. Appl. Surf. Sci. 2007, 253, 6276–6282. [Google Scholar] [CrossRef]
  37. Eason, R. Pulsed Laser Deposition of Thin Films: Applications-Led Growth of Functional Materials, 1st ed.; Wiley: New York, NY, USA, 2007; p. 681. [Google Scholar]
  38. Ageev, O.A.; Dostanko, A.P.; Zamburg, E.G.; Konoplev, B.G.; Polyakov, V.V.; Cherednichenko, D.I. Effect of the processes in the laser ablation plume on the resistivity and morphology of nanocrystalline ZnO films. Phys. Solid State 2015, 57, 2093–2098. [Google Scholar] [CrossRef]
  39. Anisimov, S.I.; Luk’Yanchuk, B.S. Selected problems of laser ablation theory. Phys.-Uspekhi 2002, 45, 293–324. [Google Scholar] [CrossRef] [Green Version]
  40. Vakulov, Z.; Geldash, A.; Khakhulin, D.; Il’ina, M.V.; Il’in, O.I.; Klimin, V.S.; Dzhuplin, V.N.; Konoplev, B.G.; He, Z.; Ageev, O.A. Piezoelectric Energy Harvester Based on LiNbO3 Thin Films. Materials 2020, 13, 3984. [Google Scholar] [CrossRef]
  41. Lee, M.H.; Park, J.S.; Kim, M.H.; Song, T.K.; Kumar, S.; Kim, W.J.; Do, D.; Hwang, I.; Park, B.H.; Choi, K.S. Lead-free piezoelectric BiFeO3-BaTiO3 thin film with high Curie temperature. Curr. Appl. Phys. 2016, 16, 1449–1452. [Google Scholar] [CrossRef]
  42. Spreitzer, M.; Egoavil, R.; Verbeeck, J.; Blank, D.H.; Rijnders, G. Pulsed laser deposition of SrTiO3 on a H-terminated Si substrate. J. Mater. Chem. C 2013, 1, 5216–5222. [Google Scholar] [CrossRef]
  43. McGee, R.; Goswami, A.; Khorshidi, B.; McGuire, K.; Schofield, K.; Thundat, T. Effect of process parameters on phase stability and metal-insulator transition of vanadium dioxide (VO2) thin films by pulsed laser deposition. Acta Mater. 2017, 137, 12–21. [Google Scholar] [CrossRef]
  44. Socol, M.; Preda, N.; Rasoga, O.; Costas, A.; Stanculescu, A.; Breazu, C.; Gherendi, F.; Socol, G. Pulsed laser deposition of indium tin oxide thin films on nanopatterned glass substrates. Coatings 2019, 9, 19. [Google Scholar] [CrossRef] [Green Version]
  45. De Bonis, A.; Teghil, R. Ultra-Short Pulsed Laser Deposition of Oxides, Borides and Carbides of Transition Elements. Coatings 2020, 10, 501. [Google Scholar] [CrossRef]
  46. Trzciński, K.; Szkoda, M.; Gazda, M.; Karczewski, J.; Cenian, A.; Grigorian, G.M.; Sawczak, M. Pulsed Laser Deposition of Bismuth Vanadate Thin Films—The Effect of Oxygen Pressure on the Morphology, Composition, and Photoelectrochemical Performance. Materials 2020, 13, 1360. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Tominov, R.V.; Smirnov, V.A.; Chernenko, N.E.; Ageev, O.A. Study of the regimes of scratching probe nanolithography. Nanotechnologies Russ. 2017, 12, 650–657. [Google Scholar] [CrossRef]
Figure 1. Laser beam scanning of the target surface in the Pioneer 180 PLD module.
Figure 1. Laser beam scanning of the target surface in the Pioneer 180 PLD module.
Materials 14 04854 g001
Figure 2. The main stages of films thickness calculation.
Figure 2. The main stages of films thickness calculation.
Materials 14 04854 g002
Figure 3. Trajectories of the laser beam movement along the target surface: scanning mode by the target with a stationary laser beam (a); scanning mode of the laser beam along the surface of a rotating target (b).
Figure 3. Trajectories of the laser beam movement along the target surface: scanning mode by the target with a stationary laser beam (a); scanning mode of the laser beam along the surface of a rotating target (b).
Materials 14 04854 g003
Figure 4. Experimental studies of a structure in a LiNbO3 film: AFM image (a); AFM cross-section (b) along the “white” line indicated in (a); 3D AFM reconstruction image (c); height histogram (d).
Figure 4. Experimental studies of a structure in a LiNbO3 film: AFM image (a); AFM cross-section (b) along the “white” line indicated in (a); 3D AFM reconstruction image (c); height histogram (d).
Materials 14 04854 g004
Figure 5. Comparison of the experimental results (points) and simulation (line) (a), as well as the distribution of the relative thickness of the LiNbO3 film on a substrate with a diameter of 100 mm (b).
Figure 5. Comparison of the experimental results (points) and simulation (line) (a), as well as the distribution of the relative thickness of the LiNbO3 film on a substrate with a diameter of 100 mm (b).
Materials 14 04854 g005
Figure 6. Dependencies of the relative thickness of the LiNbO3 film on the Upper Limit (a) and Lower Limit (b) parameters.
Figure 6. Dependencies of the relative thickness of the LiNbO3 film on the Upper Limit (a) and Lower Limit (b) parameters.
Materials 14 04854 g006
Figure 7. Dependencies of the relative film thickness on the Origin parameter (a) and Ω 95 and substrate radius on the Lower Limit, Upper Limit, and Origin (b).
Figure 7. Dependencies of the relative film thickness on the Origin parameter (a) and Ω 95 and substrate radius on the Lower Limit, Upper Limit, and Origin (b).
Materials 14 04854 g007
Figure 8. Dependencies of the relative thicknesses of the films on the parameters Max Velocity (a) and Min Velocity (b).
Figure 8. Dependencies of the relative thicknesses of the films on the parameters Max Velocity (a) and Min Velocity (b).
Materials 14 04854 g008
Figure 9. Dependencies of the relative film thicknesses on the parameters Upper Coefficient (a) and Lower Coefficient (b); Ω 95 and substrate radius dependencies on the parameters Upper Coefficient and Lower Coefficient (c).
Figure 9. Dependencies of the relative film thicknesses on the parameters Upper Coefficient (a) and Lower Coefficient (b); Ω 95 and substrate radius dependencies on the parameters Upper Coefficient and Lower Coefficient (c).
Materials 14 04854 g009
Figure 10. Projection of the trajectory of the laser beam along the target surface onto the substrate surface (a) and the distribution of the relative thickness (b) of the obtained film on a substrate with a diameter of 300 mm.
Figure 10. Projection of the trajectory of the laser beam along the target surface onto the substrate surface (a) and the distribution of the relative thickness (b) of the obtained film on a substrate with a diameter of 300 mm.
Materials 14 04854 g010
Table 1. Modeling parameters.
Table 1. Modeling parameters.
Technological Parameters
Number of laser pulses36,000
Laser pulse repetition rate10 Hz
Substrate diameter300 mm
Target diameter50 mm
Target-substrate distance75 mm
Target rotation speed10 °/s
Substrate rotation speed0.1 °/s
Scanning System Parameters (default parameters are highlighted)
Origin010254050
Upper Limit254050
Lower Limit01025
Maximum Velocity51020
Minimum Velocity0.115
Lower Coefficient15153060
Upper Coefficient
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Vakulov, Z.; Khakhulin, D.; Zamburg, E.; Mikhaylichenko, A.; Smirnov, V.A.; Tominov, R.; Klimin, V.S.; Ageev, O.A. Towards Scalable Large-Area Pulsed Laser Deposition. Materials 2021, 14, 4854. https://doi.org/10.3390/ma14174854

AMA Style

Vakulov Z, Khakhulin D, Zamburg E, Mikhaylichenko A, Smirnov VA, Tominov R, Klimin VS, Ageev OA. Towards Scalable Large-Area Pulsed Laser Deposition. Materials. 2021; 14(17):4854. https://doi.org/10.3390/ma14174854

Chicago/Turabian Style

Vakulov, Zakhar, Daniil Khakhulin, Evgeny Zamburg, Alexander Mikhaylichenko, Vladimir A. Smirnov, Roman Tominov, Viktor S. Klimin, and Oleg A. Ageev. 2021. "Towards Scalable Large-Area Pulsed Laser Deposition" Materials 14, no. 17: 4854. https://doi.org/10.3390/ma14174854

APA Style

Vakulov, Z., Khakhulin, D., Zamburg, E., Mikhaylichenko, A., Smirnov, V. A., Tominov, R., Klimin, V. S., & Ageev, O. A. (2021). Towards Scalable Large-Area Pulsed Laser Deposition. Materials, 14(17), 4854. https://doi.org/10.3390/ma14174854

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