Blog

Filament stretching during micro-extrusion of silver pastes enables an improved fine-line silicon solar cell metallization | Scientific Reports

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

Scientific Reports volume  12, Article number: 12318 (2022 ) Cite this article Nylon6 Textured Yarn

Filament stretching during micro-extrusion of silver pastes enables an improved fine-line silicon solar cell metallization | Scientific Reports

The metallization of heterojunction solar cells requires a further reduction of silver consumption to lower production costs and save resources. This article presents how filament stretching of polymer-based low-temperature curing Ag pastes during micro-extrusion enables this reduction while at the same time offering a high production throughput potential. In a series of experiments the relationship between the printing velocity and the filament stretching, thus the reduction of Ag-electrode widths and Ag laydown is evaluated. Furthermore, an existing filament stretching model for the parallel dispensing process is advanced further and utilized to calculate the elongational viscosity. The stretching effect enables a reduction of the Ag-electrode width by down to Δwf = − 40%rel. depending on the nozzle diameter and paste type. The Ag laydown has been reduced from mAg,cal. = 0.84 mg per printed line to only mAg,cal. = 0.54 mg per printed Ag-electrode when 30 µm nozzle openings are used, demonstrating the promising potential of parallel dispensing technology for the metallization of silicon heterojunction solar cells.

The International Technology Roadmap for Photovoltaic (ITRPV) predicts a world market share of silicon heterojunction (SHJ) solar cells of 10% in 2024 and 17% in 2030 which corresponds to a substantial rise compared to 3% in 20191. In the last 15 years, several research groups worked towards a further reduction of the Ag-electrode width wf and Ag laydown per cell mAg to save silver, thus further minimizing cell production costs. Lorenz et al. illustrated this trend for flatbed screen-printed Ag-electrodes (in photovoltaic industry referred to as ‘fingers’) and indicated that intense industrial optimization of pastes, screens and machine technology were the main reasons for decreasing the Ag-electrode width over the years2. In 2020, Tepner et al. presented a flatbed screen-printed line electrode with a width of wf = 19 µm and an electrode height of hf = 18 µm on a passivated emitter and rear cell (PERC)3. Besides the decrease in Ag-electrode widths, the ITRPV predicts a total silver consumption of only 50 mg silver per cell in 20301,4. In order to achieve that, the parallel dispensing technology as an alternative printing process has emerged in recent years. Pospischil et al. demonstrated a dispensed line electrode with a width of wf = 17 µm on a PERC solar cell. In that study, they showed that the Ag laydown as well as the electrical cell performance were improved compared to the reference5,6,7,8,9. These publications show an impressively successful development over the last years for PERC metallization.

However, the state of research and development of low-temperature curing Ag pastes for SHJ solar cell metallization is far from those results, especially regarding obtainable process velocities and achievable Ag-electrode widths. High throughput rates and low silver consumption per SHJ cell are requirements to increase the market share of this high-efficiency solar cell concept. Erath et al. recently published applicable flooding and printing velocities of up to v = 400 mm s−1 for flatbed screen printing10. Descoeudres et al. presented a wf = 16 µm wide, screen-printed Ag-electrode by using a special, knotless screen with screen openings of wn = 12 µm11. Our latest results for SHJ metallization by parallel dispensing showed optimized line electrode widths of wf = 34 µm and an increased optical aspect ratio of ARo = 0.55 when using 25 µm nozzle openings. In that study, the Ag laydown of a 156 mm single line was mAg = 0.30 mg electrode−112. To achieve further progress in low-temperature curing Ag pastes, an understanding of the paste’s inner state during micro-extrusion is needed to solve the limitations of applicable process velocities and line electrode widths. Until now, the reduction of Ag-electrode widths is limited because the paste shows a significant spreading tendency, thus further adaptation of the formulations is potentially required.

However, one way to achieve narrow Ag-electrode shapes may be the utilization of filament stretching during micro-extrusion13,14. In this study, we present a detailed description and experimental evaluation on how the paste threads are stretched in uniaxial dimension, resulting in significant necking and therefore reducing the thread diameter before getting in contact with the substrate. Usually, to quantify this effect, extensional properties of fluids are determined for uniaxial extension by capillary breakup extensional rheometer (CaBER)15,16 or filament stretching extensional rheometer (FiSER)17,18,19. At this point, no reliable methods to determine the extensional properties of highly filled yield stress pastes are known, therefore, we analyze the impact of this phenomenon by using the dispensing process itself. We evaluate two low-temperature curing Ag pastes regarding the maximal obtainable process velocity and its correlation to the electrode geometry. For that, we implement a variation of nozzle openings (25 µm ≤ D ≤ 45 µm) and process velocities (50 mm s−1 ≤ vprocess ≤ 500 mm s−1) and measure the Ag laydown and printed Ag-electrode widths to quantify the filament stretching effect. Further, we determine the paste’s rheological properties in shear and elongation and correlate the data to the effects observed during printing.

This paragraph summarizes the state-of-the-art literature regarding the behavior of paste threads during micro-extrusion. Clasen et al. presented the complexity of dispensing fluids with a rather intricate rheology by describing the ‘map of misery’. The ‘map of misery’ includes nondimensional numbers, e.g. Ohnesorge number Oh, Elasto-Capillary number Ec and Intrinsic Deborah number Deo, which describe the relationship between various material properties. Further characterization can be accomplished through dynamic nondimensional numbers such as capillary number Ca, Weber number We and Weissenberg number Wi20. Whenever a fluid is extruded through a nozzle, the relationship between these parameters dictate the physical evolution of the fluid after the nozzle exit. For suspensions, e.g. pastes, the thread may decrease its diameter over time due to its increasing weight. This thinning of the paste thread is mainly reinforced by the surface tension and balanced by the resisting force. It can be differentiated between viscosity-controlled thinning21, inertia-controlled thinning22,23 and elasticity-controlled thinning24.

Kunpai et al. demonstrated the filament thinning of Ag paste formulations during micro-extrusion caused by different amounts of graphite nanofibres. The widths of the dispensed structures were smaller than the applied nozzle openings because of the elongational properties of the paste13. Figure 1 depicts the stretching of paste threads in a micro-extrusion process using a parallel print head developed at Fraunhofer ISE25. Here, all paste threads exit the nozzle outlet simultaneously8. To describe the filament stretching in more detail, three distinguished velocities are defined. The paste flows from a paste reservoir through the print head and finally through the nozzle opening with an extrusion velocity vextrusion. The extrusion velocity vextrusion depends on the paste formulation and subsequently the rheological behavior as well as the process pressure p, the process temperature T and the nozzle diameter D. The process velocity vprocess defines the pace at which the substrate moves underneath the print head, mainly dictating the throughput rate of the metallization process. The specific range for process velocities that allow a stable and homogeneous printing result depends on the same influencing factors as the extrusion velocity vextrusion and additionally on the dispensing gap dgap. The extruded paste threads are free-hanging between the nozzle outlet and the substrate, thus hanging within the dispensing gap. The paste threads are then stretched, which means that the diameter of the paste threads constantly decreases from the nozzle outlet to the contact point of the paste on the substrate. The length of these paste threads can be quantified by the parameter characteristic length lc13. The wetting behavior of the strengthened paste threads onto the substrate is not considered in Kunpai’s modeling approach. However, especially the impact of textured surfaces like solar cell substrates might be significant. The extensional velocity vextension is the difference between the process velocity vprocess and the extrusion velocity vextrusion, hence the filament stretching can be enhanced by increasing the process velocity vprocess or decreasing the extrusion velocity vextrusion in relation to each other.

Overview on the parallel dispensing approach for solar cell metallization. (a) Illustration of filament stretching during micro-extrusion of low-temperature curing Ag pastes. The print head is positioned above the substrate with the dispensing gap dgap; in this gap the filament stretching occurs. The paste threads are stretched from the initial diameter d1 to a smaller diameter d2. The extensional velocity vextension is determined by the difference in process velocity vprocess and the extrusion velocity vextrusion (figure according to13). (b) Simultaneous extrusion of ten paste threads through 35 µm nozzle openings. A uniform paste distribution inside the print head ensures a homogeneous paste extrusion8,25. (c) Scanning electron microscope (SEM) image of a dispensed Ag-electrode (cross-sectional view). Here, a low-temperature curing Ag paste was applied onto SHJ solar cell by parallel dispensing. The SEM image is taken out of ref.12.

A significant filament stretching requires sufficient elongation properties of the extruded suspension, therefore a high Hencky deformation as well as a high Hencky deformation rate. The elongational viscosity for uniaxial extension is estimated by CaBER15,17,26,27 or FiSER19. The uniaxial extensional deformation corresponds to the flow of the paste thread during the dispensing process. The suspension properties in extensional flows and in shear flows can differ significantly, as shown in several publications for polymer solutions or polymer melts. Further publications show results for extensional rheological characterization of complex suspension with significant yield stress28,29,30.

A cylindrical geometry of the sample with a length L and a diameter D is assumed for determining the elongation properties. The Hencky strain ε is defined in Eq. (1), whereby L0 is the initial length and D0 the original diameter, respectively.

The change in diameter of the cylindrical sample geometry over time is defined as Hencky deformation rate \(\dot{\varepsilon }\) (Eq. (2)).

Schuemmer and Tebel define the uniaxial elongational viscosity ηe as follows in Eq. (3), whereby σzz is the axial normal stress in the fluid filament and σrr is the radial normal stress16,31.

The normal stresses σzz and σrr are defined as follows in Eqs. (4) and (5), where Γ is the surface tension and F is the axial force inside the filament.

Using the definition of the elongational viscosity (see Eq. (3)) in combination with the equations of the normal stresses (see Eqs. (4) and (5)) and the Hencky deformation rate (Eq. (2)), results in the elongational viscosity ηe in Eq. (6).

When σzz = 0 is assumed, Eq. (6) simplifies to the apparent elongational viscosity ηe,app in Eq. (7). This assumption for the axial normal stress is necessary in CaBER experiments because the axial force F in the filament cannot be determined26.

In this study, two different low-temperature curing Ag pastes A and B are analyzed regarding their rheological properties in shear and elongation as well as their printability in the dispensing process. The highly filled suspensions are developed to be applied onto transparent conducting oxide layers of SHJ solar cells and thus require a curing temperature of Tc = 200 °C to form contacts with low contact resistivity. The paste formulations consist of solvents, an epoxy-phenolic resin system, the same proportion of spherical silver nano-powders and polymer-based additives, if necessary. The solvent mixture is made up of 2-phenoxyethanol and 2-(2-hexoxyethoxy)ethanol. Above 90 wt% of the paste components are nonvolatile after the curing process.

The shear viscosity and the thixotropy of the low-temperature curing Ag pastes A and B are measured with the commercial rotational rheometer MCR 702 using the TwinDrive technology by Anton Paar GmbH, Germany. A parallel plate geometry with a diameter of D = 25 mm and a roughness of Rq = 2–4 µm is used. The roughened surface of the geometries should reduce the wall slip during the characterization of highly filled suspensions, however, the roughed surface cannot completely suppress this phenomenon32. After applying the sample onto the bottom plate geometry, the upper plate moves down to a gap distance of dtrim = 0.235 mm. Following, the excess suspension at the geometry edges is removed. The measuring gap between the upper and bottom plate is dmeasure = 0.2 mm, the measuring temperature is set to T = 25 °C. A waiting time of t = 5 min is chosen before each beginning of the measurement to achieve a homogeneous temperature distribution within the sample and to recover the inner structure of the suspension. To determine the shear viscosity, a stepwise controlled shear rate mode between \(\dot{\gamma }\)  = 10–2 s−1 to \(\dot{\gamma }\)  = 104 s−1 is applied. The measuring time per measuring point is reduced logarithmically from t = 55 s to t = 0.3 s. In total, 45 measuring points are measured. The thixotropic behavior of the highly filled suspension is estimated using the ‘three interval thixotropy test (3ITT)’33. The 3ITT method is well described in the literature34,35,36. The measuring profile is divided into three intervals whereby the frequency is kept constant at f = 1 Hz. In the first and third interval, the paste is charged at a constant low deformation with an oscillating vibration amplitude of γ = 0.1%. In this experiment, 30 measuring points are chosen with a constant measuring time of t = 10 s for the first interval and 40 measuring points with a constant measuring time of t = 15 s for each measuring point for the third interval. The third interval corresponds to the reconstruction behavior of the fluid. In the second interval, a deformation of γ = 80% is applied. The measuring time for each measuring point is set to t = 5 s for a total of 60 measuring points. In addition, the yield stress of both pastes is determined by an impeller setup. Here, the measuring gap is set to dmeasure = 0.3 mm and a shear stress between τ = 10–2 Pa to τ = 5·103 Pa is applied. The measuring time for each measuring point is set to t = 5 s. The yield stress is determined by the tangent intersection point method32,37.

At least three independent measurements are performed for each configuration. The diagrams depict the mean values of all measurement repetitions and their standard deviations as error bars. Each measurement is conducted by using a new sample. The measuring profiles of the shear rheological characterization are defined based on preliminary tests and literature that analyze also highly filled suspensions32,38,39. A high-speed imaging camera observes the paste during the rotational rheometry measurement for detecting potential errors and observing wall slip effects.

The second rheological experiment focuses on the determination of elongation properties of pastes A and B. The samples are extruded through a single nozzle with diameters of D = 110 µm, D = 160 µm and D = 230 µm. These micro-precision needles from Vieweg GmbH, Germany, have conical-shaped tips and are made of a nickel-silver alloy. The applied process pressure p is preset to a constant value. The time-dependent change of the diameter of the paste thread is determined by using high-speed imaging. For that, the high-speed imaging camera IDT OS7 with 2000 to 2500 frames per second is placed in front of the experimental setup and focused on the point of necking (see Fig. 2). At least four independent measurements are performed. Further, the high-speed imaging camera is used to determine the extrusion velocity vextrusion for nozzle openings of D = 25 µm to D = 45 µm by using our ‘GECKO’ R&D print head25. This print head extrudes ten paste threads simultaneously. The geometric design of our R&D print head is based on computational fluid dynamics simulations to guarantee a homogeneous paste distribution inside the print head6,8. The applied process pressure p is preset to the equal constant value used for the dispensing tests onto substrates. The extrusion velocity vextrusion of each paste thread is calculated by tracking the pixel movement per time increment. Three independent measurements are conducted.

Experimental setup of the rheological experiment to determine the time-dependent diameter of the extruded filament.

In this experiment, we have conducted a series of experiments to evaluate the process velocity ranges of two different low-temperature curing Ag pastes. Both pastes are extruded through nozzle openings of D = 45 µm, D = 40 µm, D = 35 µm, D = 30 µm and D = 25 µm. The process pressure p and the dispensing gap dgap are kept constant at dgap = 250 µm for all experiments. The so-called ‘GECKO’ R&D print head is used to perform all dispensing experiments25. Furthermore, a commercial table robot is used as described in literature40. The process velocity vprocess is varied between vprocess = 50 mm s−1 and vprocess = 500 mm s−1 in 10 mm s−1 increments for each nozzle diameter and paste combination. Each parameter combination is repeated independently three times. After applying low-temperature curing Ag pastes on the substrate, a curing process is conducted in a convection oven R0400FC from Essemtec AG, Switzerland. A curing temperature of Tc = 200 °C for a curing duration of tc = 5 min is applied for all samples.

After the curing step, the electrode shapes are visually evaluated regarding their homogeneity. Homogenous, straight lines of category II are characterized by using the 3D confocal laser scanning microscope OLS4000 from OLYMPUS with a magnification of 50x. Nine measurements are performed for each parameter combination. The microscope images are analyzed by the Fraunhofer ISE software, the so-called ‘Dash’41. In this case, the electrode widths wshading and wcore, the maximal electrode height hf,max as well as the cross-sectional area Across are determined. The shading electrode width wshading is defined as the maximum electrode width, including any paste spreading. The core electrode width wcore is this part of the shading electrode width which shows an apparent electrode height and thus significantly impacts the lateral electrode resistance. Based on these values, the optical aspect ratio ARo and the spreading coefficient ζspreading are calculated as given in Eqs. (8) and (9)6,39. These geometrical parameters that describe an electrode shape are visualized in a SEM image in ref.39.

Figure 3 shows the results of the rotational rheometer measurements. Both pastes show strong shear-thinning behavior meaning that the shear viscosity decreases with increasing shear rates. This rheological property is well established for metal pastes used in solar cell applications38,39. Paste A has a shear viscosity of ηpaste-A = (111.9 ± 5.0) Pa s at a shear rate of \(\dot{\gamma }\)  = 101 s−1 and a decreased shear viscosity of ηpaste-A = (5.6 ± 0.1) Pa s at a shear rate of \(\dot{\gamma }\)  = 103 s−1. Paste B shows a similar low shear viscosity (ηpaste-B = (119.9 ± 4.6) Pa s at \(\dot{\gamma }\)  = 101 s−1) (see Fig. 3a). For shear rates above \(\dot{\gamma }\) > 100 s−1, the data suggests that both paste formulations show significant differences in shear viscosity from each other, indicating different behavior during printing because relevant process shear rates are between \(\dot{\gamma }\) = 103 s−1 and \(\dot{\gamma }\) = 105 s−16. In literature, the phenomenon of edge fracture is well known for rotational rheometer measurements of highly filled suspensions42,43. This effect also occurred during our rotational measurements, hence the determined viscosity values for shear rates above \(\dot{\gamma }\)  > 100 s−1 might only give relative indications rather than absolute descriptions of the viscosity. Furthermore, wall slip, shear banding and sample spillage might affect the true shear viscosity values as described in ref.32.

Shear rheological characterization of low-temperature curing Ag pastes A and B determined by TwinDrive rotational rheometer from Anton Paar GmbH, Germany, at T = 25 °C. (a) The shear viscosity η is plotted as a function of the shear rate \(\dot{\gamma }\) , showing shear-thinning behavior of both pastes. (b) The shear deformation γ over the shear stress τ is measured to determine the static yield stress τf using the tangent intersection method. (c) In addition, the result of the three interval thixotropy test is depicted.

The thixotropic recovery of the internal structure after excessive deformation is calculated by the last elastic storage modulus G′ value of interval III and the last elastic storage modulus G′ value of interval I as indicated in Fig. 3c). Here, the relative recovery of the storage modulus is calculated after a recovery time of t = 10 min to give an indication that describes the reconstruction of the inner paste structure. Paste A shows a relative recovery of rpaste-A = (69 ± 14)%, paste B of rpaste-B = (17 ± 1)%. Even though, a recovery time of t = 10 min does not correspond to the recovery time right after the printing process, it does give a general indication on the paste’s ability to recover its inner structure. The static yield stress of paste A is τf = (7.7·102 ± 14) Pa, the static yield stress of paste B is τf = (1.2·103 ± 140) Pa (see Fig. 3b). These rheological parameters under shear flow indicate that the low-temperature curing Ag pastes A and B have similar paste properties, but as we show in this study the applicable process velocities of category II (see “Definition of process velocity range resulting in homogeneous structures” section) and the corresponding Ag-electrode shapes (see “Impact of filament stretching on solar cell metallization” section) differ significantly. In the studies of Pospischil et al.44 and Tepner et al.39, the different Ag-electrode shapes were explicable by the shear rheological properties of high-temperature curing Ag pastes. These pastes were applied by parallel dispensing or flatbed screen printing.

Figure 4 illustrates the results of time-dependent filament necking of free-hanging paste threads for both low-temperature curing Ag pastes. At the beginning of extruding a paste thread, the filament diameter stays nearly constant. After the paste thread accumulated a critical mass for which the gravitational pull overcomes the yield stress of the suspension, the filament diameter changes abruptly. The necking of the paste threads evolves until it finally ruptures. Decreasing the nozzle diameter from D = 230 µm to D = 110 µm, increases the plateau of the curve because a longer extrusion time of the paste is required to achieve the critical mass of the paste thread to overcome the yield stress of the suspension. For paste A, the initial thread diameters are significantly larger than the nozzle diameter of D = 160 µm and D = 110 µm. In contrast, paste B does not show such a strong die swell; the initial thread diameter and the nozzle diameter are similar.

Necking of free-hanging extruded paste threads of pastes A and B. The mean thread diameter Dthread is plotted over time t for nozzle diameters of D = 230 µm, D = 160 µm and D = 110 µm. The curves show an exponential evolution.

In Clasen’s study, the filament diameter as a function of the time is split into four different regimes. In this case, the elongational viscosity is determined by CaBER approach45. Based on this method, we defined two regimes for each curve in Fig. 4. The boundary condition is set as ten percent of the initial thread diameter. Therefore, the curve section including the exponential thinning is used to determine the apparent elongational viscosity ηe,app. The graph shows an exponential trend generally described by the function D(t) = D0 + A·exp(R0·t). The derivative of the function dD/dt as well as the surface tension Γ are utilized to calculate the apparent elongational viscosity ηe,app (see Eq. (7)). The surface tension is assumed to be Γ = 20 mN m−1 for both formulations6.

The apparent elongational viscosities ηe,app over the time t and the resulting exponential fits for both pastes are depicted in Fig. 5. Paste A has an apparent elongational viscosity of ηe,app|paste-A = (18 ± 6) Pa s at t = 10 ms compared to an apparent elongational viscosity of ηe,app|paste-B = (4995 ± 5848) Pa s at t = 10 ms. This high apparent elongational viscosity of paste B is especially dictated by the measuring results of nozzle diameter D = 160 µm. The other data of paste B shows an apparent elongational viscosity of one order of magnitude less. No hypothesis explains these clearly different results between the nozzle diameters. Nevertheless, a clear difference in elongational viscosity of both low-temperature curing pastes and its exponential function can be detected. A steeper gradient of graph paste B corresponds to a stronger decrease in thread diameter. Following, the resistance to stretch the paste threads is less for paste A than paste B, resulting in a more excessive filament stretching for paste A. It is pointed out that the equations of the CaBER method are taken to determine the apparent elongational viscosity based on the raw data generated by the necking experiment.

Apparent elongational viscosity of low-temperature curing Ag pastes A and B determined by filament necking at T = 25 °C. The elongational viscosity ηe,app is plotted as a function of time t.

Extruding low-temperature curing Ag pastes A and B through micrometer nozzles result in different types of Ag-electrode shapes. Figure 6 illustrates the different dispensed line shapes, which we are using to categorize the dispensing result with the focus on solar cell metallization. These microscope images show examples of the printing results for paste A with a nozzle opening of D = 45 µm. When the applied process velocity is too slow for the specified process parameters, the printed line shape has a coiled pattern (category I) (see Fig. 6, left). This means that the volume flow rate of the paste is too high for the applied process velocity. In the center of Fig. 6, the process velocity is just right in order to achieve homogeneous, straight lines (category II). However, the homogeneous line electrode widths vary significantly depending on process velocity. Consequently, there exists a velocity range from minimal process velocity vprocess,min to maximal process velocity vprocess,max for which this type of low-temperature curing Ag pastes A and B show filament stretching during the printing. The minimal process velocity vprocess,min is the lower limit for homogeneous straight line shapes, meaning that velocities below that critical threshold result in coiled lines whereas the maximal process velocity vprocess,max defines the upper limit for printing of homogeneous, straight lines. Beyond that threshold, line interruptions emerge due to the excessive differences between the extrusion and the process velocity (category III) (see Fig. 6, right). Hence, the maximum process velocity vprocess,max results in narrower line widths than the minimum process velocity vprocess,min. In the depicted example in Fig. 6, homogeneous straight line shapes could be dispensed in a process velocity range from vprocess,min = 140 mm s−1 to vprocess,max = 410 mm s−1. The shaded line electrode width wshading varies between wshading,min = 85 µm (wcore,min = 56 µm) and wshading,max = 54 µm (wcore,max = 35 µm). It indicates that the extruded paste threads are lengthened significantly during the dispensing process.

Classification of dispensed Ag-electrode shapes. Representative microscope images of different qualities of dispensed shapes of paste A (D = 45 µm) are depicted. Dispensing paste A with a process velocity of vprocess,cat-I = 90 mm s−1 results in coiled lines (category I, left). For the process velocities between vprocess,min = 140 mm s−1 (wshading = 85 µm) and vprocess,max = 410 mm s−1 (wshading = 54 µm), homogeneous, straight line shapes with different widths are observed (category II, center). At a process speed greater than vprocess,cat-III = 490 mm s−1 interrupted printing emerges (category III, right).

Figure 7 depicts the process velocity range of category II relative to the nozzle diameter and paste formulation. Both low-temperature curing Ag pastes show the tendency that the applicable range of the process velocities shrinks with a reduction of the nozzle diameter. Paste A results in homogeneous, straight line shapes when it is extruded through nozzle openings of D = 45 µm with process velocities of vprocess,min = 140 mm s−1 to vprocess,max = 410 mm s−1. Process velocities below vprocess < 130 mm s−1 at nozzle openings of D = 45 µm result in coiled lines. In contrast to that, process velocities above vprocess > 410 mm s−1 result in interrupted line shapes. When paste A is extruded through nozzle openings of D = 30 µm, the process velocity range of category II decreases down to vprocess,min = 70 mm s−1 and vprocess,max = 140 mm s−1. Furthermore, paste A cannot be extruded through nozzle openings of D = 25 µm because of its particle size distribution and agglomerate sizes. On the other hand, paste B could be extruded up to process velocity values of vprocess,max = 500 mm s−1 by using nozzle openings of D = 45 µm and D = 40 µm. Dispensing paste B through nozzle openings of D = 25 µm results in homogeneous, straight Ag-electrodes for process velocity values of vprocess,min = 170 mm s−1 and vprocess,max = 250 mm s−1. Therefore, paste B can be dispensed with faster process velocities compared to paste A, even when small nozzle diameters below D < 30 µm are used. This result is explainable by the different extrusion velocities vextrusion of both highly filled suspensions (see “Impact of filament stretching on solar cell metallization” section, Fig. 9). One possible explanation of the different paste behaviors regarding the category II process velocity range might be the difference in polymer content and the different combinations of the two polymers in the paste’s formulations.

Process velocity ranges of category II for pastes A and B are presented for nozzle openings between D = 45 µm and D = 25 µm. The process velocity range is limited by the minimum process velocity vprocess,min and the maximum process velocity vprocess,max.

Figure 8 illustrates the impact of different process velocity ranges on the electrode geometry such as the shading electrode width wshading (see Fig. 8a) and the core electrode width wcore (see Fig. 8b). The geometric variables which describe the line electrode shapes are given in ref.6,39. Here, the results of pastes A and B are shown for nozzle openings of D = 30 µm, D = 35 µm and D = 40 µm. Further results for additional nozzle diameters are summarized in Table 1. When extruding paste A through nozzle openings of D = 40 µm with the minimal process velocity of vprocess,min = 120 mm s−1, the resulting shading electrode width is wshading,paste-A = (76 ± 1) µm and the core electrode width is wcore,paste-A = (51 ± 1) µm. Increasing the process velocity up to vprocess,max = 320 mm s−1, the shading electrode width as well as the core electrode width decrease approximately linearly. The same tendency can be observed with reduced nozzle diameters. However, the gradient of the relationship becomes steeper for smaller nozzle openings. The shading electrode width of paste A can be decreased from wshading,paste-A = (60 ± 2) µm at vprocess,min = 70 mm s−1 to wshading,paste-A = (47 ± 1) µm at vprocess,max = 140 mm s−1 for nozzle openings of D = 30 µm. One important observation is made, when using nozzle diameters of D ≥ 35 µm for paste A as the extruded paste threads can result in core electrode widths smaller than the corresponding nozzle openings. However, the corresponding shading electrode widths were always larger than the applied nozzle diameters regardless of the process velocity within category II. The significant die swell tendency of paste A also occurs in the dispensing process using small nozzle diameters which means that the initial paste thread diameter is larger than the nozzle diameters resulting in wider Ag-electrodes than the used nozzle diameters.

Ag-electrode shapes depending on homogeneous process velocities of category II for low-temperature curing Ag pastes A and B. The changes of (a) shading electrode width wshading, (b) core electrode width wcore, (c) optical aspect ratio ARo, and (d) spreading coefficient ζspreading are depicted exemplarily for nozzle openings of D = 40 µm, D = 35 µm and D = 30 µm.

Paste B shows a similar tendency to paste A regarding the linear decrease of the shading electrode width and core electrode width when increasing the process velocity. Increasing the process velocity from vprocess,min = 200 mm s−1 to vprocess,max = 320 mm s−1 at D = 30 µm corresponds to an approximately linear decrease of the core electrode width from wcore,paste-B = (36 ± 1) µm to wcore,paste-B = (29 ± 1) µm. Nearly the entire process velocity range from vprocess,min to vprocess,max of paste B enables core electrode widths below the corresponding nozzle diameter.

Furthermore, Fig. 8 shows the evolution of the optical aspect ratio ARo (see Fig. 8c) and the spreading coefficient ζspreading (see Fig. 8d) relative to the process velocity range of category II. The optical aspect ratio ARo is defined as the ratio of the maximum electrode height hf,max to the shading electrode width wshading (see Eq. (8))6. The optical aspect ratio also decreases by increasing the process velocity, e.g. at the minimum process velocity vprocess,min = 120 mm s−1 an optical aspect ratio of ARo,paste-A = (0.42 ± 0.03) is achieved while the maximum process velocity vprocess,max = 320 mm s−1 results in an optical aspect ratio of ARo,paste-A = (0.32 ± 0.02) (D = 40 µm). The data suggests that when smaller nozzle openings are used, the obtainable optical aspect ratio decreases. While the decrease of the optical aspect ratio with an increasing process velocity is significant for paste A, this trend cannot be observed for paste B as the values remain almost constant.

The ratio of the core electrode width wcore to the shading electrode width wshading is defined as the spreading coefficient ζspreading with the optimal value being ζspreading = 1 (see Eq. (9))39. The spreading coefficient of paste A is equal for all evaluated nozzle diameters and process velocities in category II (see Fig. 8d). However, the spreading coefficients of paste B show a minimal tendency depending on the process velocity which could be in the range of experimental uncertainty. Smaller spreading coefficients are reached with either increasing the process velocity or decreasing the nozzle diameter. Therefore, the data suggests that smaller nozzle diameters and higher process velocities may encourage paste spreading. This leads to the following hypothesis: for slow process velocities, the shear rheological paste properties dominate and the filament stretching does not occur to a significant extend. The minimum process velocity of paste B vprocess,min = 170 mm s−1 reaches a spreading coefficient ζspreading,paste-B = (0.70 ± 0.03) and the maximum process velocity vprocess,max = 250 mm s−1 shows a spreading coefficient ζspreading,paste-B = (0.68 ± 0.04) (D = 25 µm).

With regards to solar cell application, the applied metal grids should have narrow Ag-electrodes with no paste spreading and minimal silver consumption per cell in combination with high throughput rates, following high process velocities. Our data suggests that these goals follow a severe trade-off. The process velocity range of category II enables the reduction of the line electrode width as increasing process velocities result in decreasing core electrode widths as well as shading electrode widths. Using a smaller nozzle diameter helps to reduce the Ag-electrode width at the expense of significantly limiting the applicable process velocity while causing potential problems in process stability due to clogging of those nozzles because of the paste’s tendency to form agglomerates. In addition, smaller nozzle openings can encourage the paste spreading depending on the paste system. Nevertheless, dispensing low-temperature curing Ag pastes through small nozzle diameters results in the narrowest line electrodes.

Further, we have evaluated the Ag laydown per 156.75 mm line in relationship to the process velocity in category II (Table 1). For example, extruding paste A through nozzle openings of D = 45 µm results in a Ag laydown of mAg,paste-A = (1.90 ± 0.10) mg electrode−1, when using the minimal process velocity of vprocess,min = 140 mm s−1. Increasing the process velocity up to vprocess,max = 410 mm s−1 results in a decrease of paste laydown down to mAg,paste-A = (0.71 ± 0.03) mg electrode−1, following a delta reduction of Ag laydown by ΔmAg,paste-A = − 1.19 mg electrode−1. The Ag saving potential of paste B in combination with nozzle diameters of D = 45 µm is only ΔmAg,paste-B = − 0.66 mg electrode−1. Decreasing the nozzle diameter reduces the possible magnitude of Ag reduction as paste A shows a Ag laydown between mAg,paste-A = (0.85 ± 0.04) mg electrode−1 at vprocess,min = 70 mm s−1 and mAg,paste-A = (0.55 ± 0.03) mg electrode−1 at vprocess,min = 140 mm s−1 for 30 µm nozzle openings. Thus, filament stretching enables the control of Ag laydown per solar cell.

In conclusion, the maximum usage of filament stretching during micro-extrusion allows narrow Ag-electrodes and low Ag laydowns combined with high throughput potential. Consequently, the full usage of filament stretching should result in a gain in solar cell efficiency, as narrower Ag-electrodes decrease the shading area and thus increase the short-circuit current density. Further, the silver usage is more effective for homogeneous electrode shapes than inhomogeneous ones. We have evaluated this hypothesis in a different study, and the results are available in ref.46.

The question arises why low-temperature curing Ag pastes display such a vast process velocity range within category II allowing for such different electrode shapes as well as Ag laydown. The significant amount of polymer in the paste formulation might be one reason. For elastic materials like polymer-based suspensions, the elongational viscosity is reported as an influential property in the literature. During micro-extrusion, the paste is subject to uniaxial extensional deformation. Kunpai et al. demonstrated the different stretching of Ag pastes caused by varying amounts of graphite nanofibres. The resulting electrode widths were smaller than the applied nozzle openings because of the elongational properties of the pastes13. This elongational property of low-temperature curing Ag paste opens up the possibility to influence the Ag-electrode width and further the Ag laydown per cell without reducing the nozzle opening which can increase the process stability.

Paste A has a lower elongational viscosity than paste B (see “Behavior under uniaxial extensional flow” section) which corresponds with the results of the printing experiment. The different electrode shape parameters as well as the different corresponding Ag laydown between the minimal and the maximal process velocity is higher for paste A than paste B, therefore the filament stretching of paste A appears to be greater. The lower elongational viscosity of paste A reflects a lower resistance regarding uniaxial extensional deformation which is also evident in the different characteristic lengths lc. Figure 9b) illustrates the characteristic length in dependence on the strain rate for different nozzle diameters. When the strain rate is \(\dot{\varepsilon }\)  = 0, the characteristic length lc corresponds approximately to the corresponding dispensing gap dgap. This effect is independent of the nozzle diameter or paste formulation as the paste thread does not undergo any extensional deformation. Additionally, the measurement data allows the following assumption: vprocess,min = vextrusion. The deviations between the two experimental data are within the range of experimental uncertainty and result partly from the chosen experimental parameters; for instance the minimal process velocity vprocess,min is only determined in 10 mm s−1 increments. Increasing the strain rate results in an increased characteristic length, following the paste properties affect the characteristic length. Paste A shows a maximum characteristic length of lc = 604 µm while paste B has a maximum characteristic length of lc = 372 µm. It seems that the characteristic length is independent of the nozzle diameter for the particular paste within experimental uncertainty. Even though paste B shows a significantly higher extrusion velocity compared to paste A, paste A has a higher stretching effect (see Fig. 9a)).

(a) Extrusion velocities vextrusion for various nozzle diameters and (b) its characteristic length lc over strain rate for low-temperature curing Ag pastes A and B. The extrusion of paste threads is observed by a high-speed imaging camera to determine the extrusion velocity vextrusion for each paste depending on the nozzle diameter of 45 µm ≥ D ≥ 25 µm. Besides, the minimal process velocities vprocess,min of Fig. 6 are shown to visualize our assumption vprocess,min = vextrusion. The characteristic length lc is calculated based on the equation in13.

Another link between the paste properties and the ratio of the cross-sectional area of the dispensed line Across to the cross-sectional area of the corresponding nozzle Anozzle is presented in Fig. 10. The combination of different process velocities and extrusion velocities can result in the same extensional velocity, resulting in a comparable ratio of the two cross-sectional areas, especially for extensional velocities above vextension = 102 mm s−1 this relationship seems to appear to be independent of the suspension and the nozzle diameter. When the ratio of the cross-sectional areas becomes greater than one, the paste threads swell at the nozzle outlet or the wetting behavior of the paste onto the substrate becomes significant superimposing the filament stretching effect. When the ratio of the cross-sectional areas remains below one, the filament stretching effect dominates.

Correlation between elongation properties of the low-temperature curing Ag pastes A and B and the resulting electrode shapes. The ratio of cross-sectional area of the dispensed Ag-electrode Across to the cross-sectional area of the nozzle Anozzle are plotted over the extensional velocity vextension.

In this study, we have demonstrated the filament stretching of low-temperature curing Ag pastes during micro-extrusion, enabling a significant decrease of line electrode width and Ag laydown by increasing the corresponding process velocities. In a series of rheological experiments, we have shown that two pastes with similar shear rheological properties differ significantly in elongational behavior. Further, the elongational viscosity for uniaxial extensional deformation for both pastes differs and reflects the degree of filament stretching and the potential characteristic length of the threads and its impact on the printing process.

The stretching effect allows to reduce the Ag-electrode width by down to Δwf = − 40%rel. depending on the nozzle diameter and paste type. Further, a reduction of Ag laydown from mAg = 0.84 mg per Ag-electrode to mAg = 0.54 mg per line electrode for 30 µm nozzle openings has been demonstrated. These results show a promising way to further decrease the metallization costs while improving production throughput for SHJ solar cells and potentially enhancing solar cell performances. In order to achieve this, we suggest the development of new paste formulations which show a strong filament stretching effect and at the same time a low die swell tendency by adjustment of the polymer matrix.

The datasets generated during and/or analysed during the current study are not publicly available due to an agreement between project partners in the project ALTURA with 03EE1006C number but are available from the corresponding author on reasonable request.

VDMA Photovoltaics Equipment. International Technology Roadmap for Photovoltaics (ITRPV) 2019 Results. 11th ed (2020).

Lorenz, A. et al. Screen printed thick film metallization of silicon solar cells—recent developments and future perspectives. In Proc. 35th European Photovoltaic Solar Energy Conference and Exhibition, Brussels, Belgium, 819–824 (2018).

Tepner, S. et al. Screen pattern simulation for an improved front-side Ag-electrode metallization of Si-solar cells. Prog. Photovolt Res. Appl. 28, 1054–1062 (2020).

Baliozian, P. et al. The international technology roadmap for photovoltaics and the significance of its decade-long projections. In Proc. 37th European Photovoltaic Solar Energy Conference and Exhibition, Lisbon, Portugal, 420–426 (2020).

Pospischil, M. et al. Applications of parallel dispensing in PV metallization. AIP Conf. Proc. 2156, 20005 (2019).

Pospischil, M. Dissertation (Albert Ludwig University, 2016).

Pospischil, M. et al. high speed dispensing—a high-throughput metallization technology for >21% PERC type solar cells. In Proc. 32nd European Photovoltaic Solar Energy Conference and Exhibition, Munich, Germany https://doi.org/10.4229/EUPVSEC20162016-2CO.2.2 (2016).

Pospischil, M. et al. Development of a high-throughput fine line metallization process using CFD-simulation. In 39th Photovoltaic Specialists Conference (PVSC), 2250–2253 (2013).

Pospischil, M. et al. Optimizing fine line dispensed contact grids. Energy Procedia 55, 693–701 (2014).

Erath, D. et al. Fast screen printing and curing process for silicon heterojunction solar cells. AIP Conf. Proc. 2367, 20006 (2021).

Descoeudres, A. et al. Low-temperature processes for passivation and metallization of high-efficiency crystalline silicon solar cells. Sol. Energy 175, 54–59 (2018).

Article  ADS  CAS  Google Scholar 

Gensowski, K. et al. Dispensing of low-temperature silver pastes. AIP Conf. Proc. 2367, 020007 (2021).

Kunpai, C., Kang, M. G., Song, H. & Shin, D.-Y. Fine front side metallisation by stretching the dispensed silver paste filament with graphite nanofibres. Sol. Energy Mater. Sol. Cells 169, 167–176 (2017).

Udofia, E. N. & Zhou, W. Microextrusion based 3D printing—a review. In Proceedings of the 29th Annual International Solid Freeform Fabrication Symposium—An Additive Manufacturing Conference, Austin, TX, USA, 13–18 (2018).

Bazilevsky, A. V., Entov, V. M. & Rozhkov, A. N. Liquid filament microrheometer and some of its applications. In Third European Rheology Conference and Golden Jubilee Meeting of the British Society of Rheology (1990).

Sachsenheimer, D. Dissertation.Karlsruhe Institute of Technology, Karlsruhe (2014).

Rodd, L. E., Scott, T. P., Cooper-White, J. J. & McKinley, G. H. Capillary break-up rheometry of low-viscosity elastic fluids. Appl. Rheol. 15, 12–27 (2005).

Arnolds, O. Dissertation.Karlsruhe Institute of Technology, Karlsruhe, (2011).

Matta, J. E. & Typtus, R. P. Liquid stretching using a falling cylinder. J. Non-Newton. Fluid Mech. 35, 215–229 (1990).

Clasen, C., Phillips, P. M., Palangetic, L. & Vermant, A. J. Dispensing of rheologically complex fluids: The map of misery. AIChE J. 58, 3242–3255 (2012).

Papageorgiou, D. T. On the breakup of viscous liquid threads. Phys. Fluids 7, 1529–1544 (1995).

Article  ADS  MathSciNet  CAS  Google Scholar 

Day, R. F., Hinch, E. J. & Lister, J. R. Self-similar capillary pinchoff of an inviscid fluid. Phys. Rev. Lett. 80(4), 704–707 (1998).

Article  ADS  CAS  Google Scholar 

Tirtaatmadja, V., McKinley, G. H. & Cooper-White, J. J. Drop formation and breakup of low viscosity elastic fluids: Effects of molecular weight and concentration. Phys. Fluids 18, 43101 (2006).

Entov, V. M. & Hinch, E. J. Effect of a spectrum of relaxation times on the capillary thinning of a filament of elastic liquid. J. Nonnewton. Fluid Mech. 72, 31–53 (1997).

Pospischil, M. et al. Process development for a high-throughput fine line metallization approach based on dispensing technology. Energy Procedia 43, 111–116 (2013).

Sachsenheimer, D., Hochstein, B., Buggisch, H. & Willenbacher, N. Determination of axial forces during the capillary breakup of liquid filaments—the tilted CaBER method. Rheol. Acta 51, 909–923 (2012).

Arnolds, O., Buggisch, H., Sachsenheimer, D. & Willenbacher, N. Capillary breakup extensional rheometry (CaBER) on semi-dilute and concentrated polyethyleneoxide (PEO) solutions. Rheol. Acta 49, 1207–1217 (2010).

McDonnell, A. G. et al. Extensional viscosity of copper nanowire suspensions in an aqueous polymer solution. Soft Matter 11, 8076–8082 (2015).

Article  ADS  CAS  Google Scholar 

Niedzwiedz, K., Arnolds, O., Willenbacher, N. & Brummer, R. How to characterize yield stress fluids with capillary breakup extensional rheometry (CaBER)?. Appl. Rheol. 19, 41969-1-41969–10 (2009).

Niedzwiedz, K., Buggisch, H. & Willenbacher, N. Extensional rheology of concentrated emulsions as probed by capillary breakup elongational rheometry (CaBER). Rheol. Acta 49, 1103–1116 (2010).

Schümmer, P. & Tebel, K. H. A new elongational rheometer for polymer solutions. J. Non-Newton. Fluid Mech. 12, 331–334 (1983).

Yüce, C. & Willenbacher, N. Challenges in rheological characterization of highly concentrated suspensions—a case study for screen-printing silver pastes. J. Vis. Exp. https://doi.org/10.3791/55377 (2017).

Article  PubMed  PubMed Central  Google Scholar 

Mezger, T. The Rheology Handbook. For Users of Rotational and Oscillatory Rheometers 5th edn. (Vincentz Network, 2020).

Barnes, HA Thixotropy—A review.J. Nonnewton.Fluid Mech.70, 1–33 (1997).

Mewis , J. & Wagner , NJ Thixotropy .Adv.Coll.Interface.Sci.147–148, 214–227 (2009).

Larson, R. G. & Wei, Y. A review of thixotropy and its rheological modeling. J. Rheol. 63, 477–501 (2019).

Article  ADS  CAS  Google Scholar 

Brummer, R. Rheology Essentials of Cosmetic and Food Emulsions. 18 tables (Springer, 2006).

Thibert, S. et al. Influence of silver paste rheology and screen parameters on the front side metallization of silicon solar cell. Mater. Sci. Semicond. Process. 27, 790–799 (2014).

Tepner, S. et al. The link between Ag-paste rheology and screen-printed solar cell metallization. Adv. Mater. Technol. 5, 2000654 (2020).

Gensowski, K. et al. Conductive highly filled suspensions for an electrochemical dispensing approach to pattern full-area thin metal layers by physical vapour deposition. Sci. Rep. 10, 1–10 (2020).

Wenzel, T. et al. Progress with screen printed metallization of silicon solar cells—towards 20 µm line width and 20 mg silver laydown for PERC front side contacts. SSRN 244, 111804 (2022).

Mattes, K. M., Vogt, R. & Friedrich, C. Analysis of the edge fracture process in oscillation for polystyrene melts. Rheol. Acta 47, 929–942 (2008).

Hemingway, E. J., Kusumaatmaja, H. & Fielding, S. M. Edge fracture in complex fluids. Phys. Rev. Lett. 119, 28006 (2017).

Pospischil, M. et al. Paste rheology correlating with dispensed finger geometry. IEEE J. Photovoltaics 4, 498–503 (2014).

Clasen, C. Capillary breakup extensional rheometry of semi-dilute polymer solutions. Korea-Austral. Rheol. J. 22, 331–338 (2010).

Gensowski, K. et al. Filament stretching during parallel dispensing—a way to reduce Ag consumption in SHJ metallization. Solar Energy Mater. Solar Cells. https://doi.org/10.1016/j.solmat.2022.111871 (2022).

This work was supported by the German Federal Ministry for Economic Affairs and Energy within the research project “Altura” under the contract number 03EE1006C. The authors would like to all co-workers at Fraunhofer Institute for Solar Energy Systems ISE. The authors are responsible for the content.

Open Access funding enabled and organized by Projekt DEAL.

Fraunhofer Institute for Solar Energy Systems, Heidenhofstraße 2, 79110, Freiburg im Breisgau, Germany

Katharina Gensowski, Maximilian Much, Elisabeth Bujnoch, Stefan Spahn, Sebastian Tepner & Florian Clement

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

K.G. conducted the entire literature research, wrote the main manuscript text and prepared all figures for publication. K.G. performed all printing experiments, E.B. supported these experiments. M.M. did most experimental work regarding the rheological characterization of printing pastes, supervised by K.G. and ST.. S.S. supported these rheological measurements. M.M. and S.S. assisted in the data analysis procedure. E.B. did most of the characterization of printed structures, supervised by K.G.. S.T. supervised all research activities, especially the data analysis and interpretation of results. K.G. and S.T. developed the idea and concept of the study. F.C. enabled the financial resources of this work. All authors—K.G., M.M., E.B., S.S., S.T. and F.C.—reviewed the manuscript and provided feedback on the manuscript before submission.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Gensowski, K., Much, M., Bujnoch, E. et al. Filament stretching during micro-extrusion of silver pastes enables an improved fine-line silicon solar cell metallization. Sci Rep 12, 12318 (2022). https://doi.org/10.1038/s41598-022-16249-5

DOI: https://doi.org/10.1038/s41598-022-16249-5

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

Scientific Reports (Sci Rep) ISSN 2045-2322 (online)

Filament stretching during micro-extrusion of silver pastes enables an improved fine-line silicon solar cell metallization | Scientific Reports

Variegated Yarn Knitting Patterns Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.