banner

Blog

Apr 22, 2024

Parametric study of pulsed nanosecond laser interaction with carbon

Scientific Reports volume 13, Article number: 2048 (2023) Cite this article

1309 Accesses

2 Citations

2 Altmetric

Metrics details

A laser processing technique is proposed for the processing of a 2.5 mm thick carbon nanotube (CNT) composite bipolar plate for proton exchange membrane fuel cells (PEMFCs). This study aims to understand laser interaction with the CNT composite plate experimentally using a pulsed nanosecond laser. Penetration depth, top width, spatter width, and overall physical morphologies are studied. Scanning electron microscope (SEM) and 3D Scanning Confocal Microscope were used for observation and measurements. Based on that, a parametric investigation is conducted and reported systematically. Most importantly, the pulse repetition rate presents a unique nature of interaction that resulted in a critical repetition rate distinguishing three operational regimes. The physical and chemical properties of the regimes are further analyzed by Vickers microhardness testing and energy dispersive X-ray (EDX) analyses performed on the surface and cross-section of each specimen. The results reveal that the pulse repetition rate introduces changes in mechanical properties and chemical compositions in the vicinity of the processed region. In conclusion, lower pulse repetition should be favored for less impact on mechanical properties, chemical composition, and morphological aspects.

Carbon nanotubes (CNTs) have remarkably outstanding mechanical properties (elastic modulus and tensile strength) and superior electrical and thermal conductivities, making them stiff and strong conducting materials with low weight compared to steel and other structural materials1. This stimulated a great deal of attention in the advanced composite society in their use as reinforcement materials for the advancement of composite materials2. These composites are being used for various applications in wearable systems (smart textiles), robotics, and next-generation electronics devices and energy conversion systems3,4,5. In addition to the remarkable mechanical, electrical, and thermal properties, CNTs possess high surface-to-volume ratios resulting from their small diameters, which are order of few nanometers. This has created a great opportunity for CNT composites, particularly in battery and energy conversion devices where increased effective contact surface area of electrodes per volume plays a significant role in energy conversion efficiency6,7,8. CNTs are identified as potentially useful enforcement in lithium-ion battery systems, fuel cells, and solar cells9,10,11. A CNT composite with a metal nanoparticle as an electrode doubles the performance of hydrogen fuel cells due to the increased catalytic activity of carbon nanotube-based electrodes12,13. Other studies pointed out the relevance and applications of CNTs in lithium-ion batteries14, elastic and transparent conductive films15, and flat panel displays16.

Following the increasing demand and applicability of CNTs in various technological fields, developing effective manufacturing processes is vital for processing CNT composites to the desired size, shape, and quality. Any manufacturing method selected for processing CNT composites must ensure minimum damage to the CNT structure that could result from pressure, heat, or chemical reaction with the matrix material. Conventional processes like machining and molding are associated with certain disadvantages. Since CNTs have high strength and hardness, conventional machining methods would cause intensive tool wear, reducing the tool life and increasing production cost17. The alignment of CNTs in molded composites is significantly affected by the shear flow in the molding process leading to undesired changes in its structures and properties18.

The continuous advancements in the performance of lasers during the past decades have improved their capacity in various fields, including energy, biotechnology, electronics, and mechanical engineering19. In polymeric composite cutting, lasers offer variety of advantages, including high production speeds without shortcomings associated with tool wear and vibration20,21. Lasers are specifically advantageous in the processing of difficult-to-machine materials22, like carbon fiber composite and graphite composites due to their brittleness and hardness.

Many researchers have studied the use of lasers in processing CNTs and their composites. Wu et al.23 performed the ablation and patterning of CNT film using a femtosecond laser for applications in electroluminescent and flexible electronic components. Raman spectroscopy and scanning electron microscope (SEM) were used to characterize the performance of the pattern groove. The research indicated the influence of process parameters on ablation and pattering at different pulse energies. In addition, appropriate cutting parameters that introduce minimum defects were suggested. Chen et al.24 studied the structural modification and transformation of CNT by using ND: YAG laser operating at 266 and 1067 nm wavelengths, with different energy fluences and the number of passes to study the structural modification of CNT. Accordingly, it was demonstrated that transformation on a selective area could be effectively controlled by laser fluence and the number of passes. Effects of processing parameters on laser cutting of multi-walled carbon nanotubes/polymethyl methacrylate nanocomposites are investigated by Azmats et al.25 for applications of CNTs as a reinforcement for plastics to take advantage of their good mechanical, geometrical, electronic, and electromagnetic properties. Furthermore, the research suggested that the number of CNT is an influential factor in reducing HAZ by up to 50%. Overall, these researches emphasized the need of packing of CNTs into microscopic materials: films, sheets, and ribbons that are flexible and thin structures and need to be cut or structured into nanometer and micrometer dimensions using lasers. On the other hand, CNTs are also packed into materials at the macro level to produce CNT composites for applications such as bipolar plates of PEMFCs for their effective charge transfer and thermal management capabilities26. Developing effective manufacturing methods for bipolar plate is very crucial because bipolar plate is a key component in fuel cells with multiple functional roles, contributing up to 40–45% of the total fuel cell cost27,28. The current technological advancement has introduced a variety of advanced lasers with a variety of capacities and applications that can remedy the challenges in the machinability of CNT composite bipolar plates due to their brittleness and hardness. The use of direct laser melting (DLM) technology to fabricate the flow channels/patterns of metallic bipolar plates of direct methanol fuel cell (DMFC) was identified from the work of Moon et al.29 who compared its performance to numerical control (NC) machined bipolar plate. Results show that the overall fuel cell performance increases for DLM micro-patterned bipolar plate because of rough side walls from the laser operation that causes more rapid chemical reactions. Despite these earnest efforts and achievements, no attempt has been made to overcome the challenge of flow channel fabrication on CNT composite bipolar plates using lasers.

In this study, laser processing is proposed for the manufacturing of fuel and oxidizer flow channels on a CNT composite bipolar plate made by blending graphite, epoxy, and CNTs. However, understanding the laser-material interaction is an essential step and should be a prerequisite condition before proceeding with the suggested solution. Therefore, this paper aims to clearly understand the interaction of nanosecond laser with CNT composite bipolar plate. To assess the performance of laser in the processing of CNT composites, the effects of power, scanning speed, the number of passes, repetition rate, and pulse duration on the penetration depth, top width, spatter width (if spatter forms) and physical morphology of the specimen are systematically analyzed and discussed. Furthermore, chemical characterizations and mechanical property measurements are performed for parametric settings that yield special morphological appearances.

Samples are prepared from a 2.5 mm thick CNT composite plate that is manufactured for bipolar plates in PEMFCs, shown in Fig. 1a. The CNT composite plate material is fabricated from graphite (77% by weight), a polymer, specifically, a thermoset (epoxy) polymer (20% by weight), and CNTs as filler materials (3% by weight). The CNTs are highly conductive multi-walled CNTs with a diameter of 11–13 nm, a length of 40–50 μm, a purity of 95 wt.%, and a bulk density of 0.025 g/ml. A molding process at 50 MPa molding pressure was employed during its production. Optical properties of the material are of immense importance for its laser processibility. Figure 1b presents the UV–Vis-NIR spectrophotometer analysis of the sample for the reflectivity rate, absorptivity rate and transmittance of the 2.5 mm thick CNT composite plate sample at various wavelengths. A chemical characterization of the plate (which can be referred as unprocessed base material) is done using EDS analysis, and elemental mappings are shown in Fig. 2. The average values of elemental compositions are tabulated in Table 1.

(a) A CNT composite plate (b) reflectivity, absorptivity and transmittance of the plate in (a).

EDX mapping of unprocessed base material.

Ytterbium nanosecond pulsed Nd:YAG fiber laser (IPG-YLPM, IPG photonics, model YLP-HP IPG photonic, Southbridge, MA, USA) was used as a laser source. The laser source has an emission wavelength of 1064 nm, generating up to 20W average power, 1000 kHz repetition rate, 200 ns pulse duration, and a scanning speed of 2000 mm/s. The source has a beam quality factor (M2) of 1.5, a collimated beam diameter of 12.8 mm and a spot diameter of 30 μm at 189 mm focal distance. Figure 3a illustrates the schematics of the experimental setup. Preliminary experiments suggested that ablation of CNT composites was achieved from low laser powers to the highest [4–20 W] and from low scanning speed to moderate values [up to 300 mm/s]. For this reason, laser powers of 4–20 W with an interval of 2W and scanning speeds of 50–300 mm/s at an interval of 50 mm/s were considered. Each pulse duration (4, 20, 50, 100, and 200 ns) has its own set of working frequency ranges on the pulsed laser. For instance, with lower pulse duration values, the operation can only be performed at higher repetition rate values, and operations on lower repetition values became possible at higher pulse duration settings. Hence, at 200 ns, all frequency settings [20–1000 kHz] are active. So, typical frequency values for the onset of new pulse duration (20, 40, 60, 105, 500, and 1000 kHz) were considered for the study. The experimental parameters are presented in Table 2. The same table can be referred to understand frequency and pulse duration couplings. To understand the interaction between the laser and the CNT composite plate, laser irradiations are performed by line exposures on the specimen in a configuration given in Fig. 3b. Line irradiations are simple and proper ways to examine the ablation characteristics and understand the laser interaction with the material.

(a) Experimental setup and (b) laser irradiation path.

The response parameters are penetration depth, top width, and spatter width. The penetration depth is the actual depth created by the laser beam, while the top width is the maximum width of the material removed by the laser. Spatter width is defined as the average width of the spatter dispersed in the vicinity of the processed region. Measurement results are obtained from a high-precision 3D scanning Leica Confocal Microscope (Leica PLANAPO FOV 3.6, DMI, USA). The digital microscope has magnification from 12x to 2340x depending on the objective selected (low or high). 2340x (high magnification) can show details down to 0.4 µm. It has 10 megapixel integrated high resolution camera with integrated ring ling and coaxial LED illumination and a tilting stand (− 60° to + 60°). 3D topology map from the confocal digital microscope can capture topological variations at inner section of the groove and also it captures possibilities of ablated particles trapping and molten spatter formation inside the groove. The 3D scanner can analyze a large area and calculate the mean depth, which provides a more reliable data there by avoiding the need for repeated experiments. Accordingly, the study is conducted on a full factorial experimental design based on an experiment performed on each parametric combination. The response parameters and the measurement device with the measurement methods are schematically illustrated in Figs. 4 and 5, respectively. In addition, SEM was used to thoroughly observe and report important morphological feature resulting from the interaction with the laser. In cases when chemical and mechanical property characterizations were important, Electron-dispersive X-ray Spectroscopy (EDX) and Vickers Microhardness Tester are employed. EDX analysis is also carried out for the characterization of the spatter with respect to the base material. The apparatus used for SEM and EDS is MIRA 3-LMH High Resolution Emission Scanning Microscope with high brightness Schottky emitter electron gun (manufactured by Tescan). It has a resolution of 1.0 nm at 30 kV and 2.0 nm at 3 kV and a magnification up to 1,000,000 × and magnification up to 4 × without any image distortion in wide field optic mode. Additional specifications include, accelerating voltage of 200–30 keV and probe current of 1 pA to 100 nA. The Vickers microhardness tests are conducted with a microhardness testing machine (HM, Mitutoyo Corporation, Japan) to determine the hardness or specimen’s resistance to deformation. The CNT composite plate manufactured with a very smooth finish (measured roughness of 0.1534 μm), therefore, it does not require polishing. The Vickers hardness tester uses a diamond indenter that applies a 0.49N force (50 g-force) on the surface of the specimen for a dwell time of 10 s.

Illustrations of top width, penetration depth and spatter width.

Measuring device and measurement methods.

Laser power (average laser power) is a measure for the energy of light delivered by the beam per unit time, which is a commonly used parameter to control a laser process30. This section presents the effect of laser power on laser processing of CNT composites from 4 to 20W at fixed parameters of 200 ns, 20 kHz, 100 mm/s and 1 pass and the effect of scanning speed from 50 to 300 mm/s at fixed parameters of 20W, 200 ns, 20 kHz and 1 pass. In laser processing of CNT composite, it was observed that low to high powers could be used depending on the extent to which the material is needed to be ablated. Experimental results show that power causes both penetration depth and top width to increase. Figure 7a illustrates the effect of power on the penetration depth and top width. Furthermore, the effect of power on the physical morphology of the sample was observed. Based on that, no significant effects were observed. Physical morphologies of the specimens at low power (4W), medium power (10W) and high power (20W) are given in Fig. 6a, b and c.

Effect of power on physical morphology; (a) 4W, (b) 10W, (c) 20W [Fixed: 20 k Hz, 200 ns, 100 mm/s and 1pass].

It was identified that most parametric combinations in laser processing of CNT composites lead to the dispersion of spatter from the processed zone to the periphery of the laser irradiated line. Observations made on the specimens before and after the cleaning of the spatter is presented in Table 3. Cleaning of the spatter is done by careful wiping using a general cleaning wet microfiber cloth to avoid skin contact and inhaling by the worker and for the safety of the specimen as well. Specific issues such as causes of spatter formation and characterization of the spatter are discussed in the upcoming sections. According to the study conducted on the formation of spatter (in "Effect of pulse duration and repetition rate" Section), it was indicated that repetition rate is the key parameter for spatter formation. More specifically, lower repetition rates produce a spatter of larger width. Since the effect of power is studied at 20 kHz, which is a pulse repetition rate that gives off abundant spatter, it was possible to observe spatter. The relationship between power and the resulting spatter width can be seen in Fig. 7b. It could also be seen that more spatter is formed due to the increase in material removal as power increases.

Effects of power; (a) Effect on top width and penetration depth, (b) Effect on spatter width [Fixed: 20 kHz, 200 ns, 100 mm/s and 1pass].

Scanning speed represents the duration of laser interaction, implying more interaction at low speeds and less at high speeds. In the laser processing of CNT composites, the extended period of interaction at low speeds [50–100 mm/s] causes deeper penetration, but the high-speed operations [150–300 mm/s] cause a shallow penetration. From Fig. 9a, the effect of scanning speed is observed to be primarily on the penetration depth, and its effect on top width is almost negligible. Scanning speed presents no significant morphological changes on the CNT composite. As shown in Fig. 8a, low scanning speed gives a longer interaction that offers high penetration depth. Increasing the scanning speed to 100, 200, and 300 mm/s, as in of Fig. 8b, c and d respectively, the penetration depth decreases proportionally due to a reduction in ablation that resulted from short interactions with increased scanning speed. In addition, the process is partially spatter-forming (shown in Fig. 9b) which gives off the spatter at low to medium scanning speed levels (50–150 mm/s) only. In contrast, higher scanning speeds (200–300 mm/s) do not result in a significant amount of spatter and in a consistent and well-established form (in a fashion presented in Table 3).

Effect of scanning speed on physical morphology; (a) 50 mm/s, (b) 100 mm/s, (C) 200 mm/s (d) 300 mm/s [Fixed: 20 kHz, 200 ns, 20W and 1pass].

Effects of scanning speed; (a) Effect on top width and penetration depth, (b) Effect on spatter width [Fixed: 20 kHz, 200 ns, 20W and 1pass].

Generally, laser power and scanning speed can be used to control the penetration depth and top width of the laser-irradiated zone in smooth operation without having any form of significant changes in the physical morphology.

The number of passes is a parameter that measures the number of laser irradiation passes. In this section, the effect of number of scan passes is investigated from 1 to 20 passes at 4 passes interval and at a fixed parameters of 20W, 20 kHz, 200 ns and 100 mm/s. In processing CNT composite, the number of passes strongly impacts the processability. As shown in Fig. 11a, a step change in the number of passes greatly affects penetration depth and top width. The formation of wider spatter zones is an additional characteristic feature of the process (Fig. 11b). Additionally, each number of pass display unique forms of a particular morphology. Three main phenomena are observed here. First, at a lower number of passes (up to 8 passes), the ablation becomes powerful to produce deep and wide grooves, as shown in (a) and (b) of Fig. 10. Here, the heat from the laser beam diffuses significantly in the lateral (horizontal) direction and causes material removal from the edges to yield wider channels near the top. Gradual increase of the number of passes to higher values (12 passes), as in the edge view of the case (c) of Fig. 10, begins to initiate crack-like narrow grooves (sharp groove) at the bottom. This effect increases with the number of passes to 16, as in case (d) of Fig. 10, where the effect is amplified with the creation of a wider grooves at the top and sharp/pointy at the bottom. This condition cannot last longer because increasing the top width with the creation of an even narrower groove leads to the deposition of particles to the inside, as shown in case (e) of Fig. 10. These trapped particles could originate from two mechanisms: (1) due to the increase in top width leading to the formation of the large number of particles ablated from the top surface that scatter into the narrow and sharp groove, and (2) due to very high number of the pass at scanning speed of 100 mm/s, that gives less time for removal of particles from the narrow groove (Fig. 11). At this point, it will become a little difficult for the laser beam to penetrate more profoundly because of the blockage of pulses by these particles. As a result, the penetration depth is negatively affected after 16 passes. For clarity, 3D topological details and line profiles at 16 passes and 20 passes are given in Fig. 12.

Effect of number of passes; (a) 1pass, (b) 8passes, (c) 12passes, (d) 16passes and (e) 20passes [Fixed: 20 kHz, 200 ns, 20W and 100 mm/s].

Effect of number of passes at fixed values of 20 kHz, 20W and 20 mm/s; (a) Effect on top width (b) Effect on spatter width [Fixed: 20 kHz, 20W, 200 ns and 100 mm/s].

Line profile views and groove morphologies at 16 passes and 20 passes.

Pulse duration refers to the time duration between the instant a pulse starts and the instant the pulse ends31. This parameter is crucial for pulsed lasers because the laser does work only during this period. The amount of energy applied during this period determines the extent to which the process removes material. Depending on the pulse durations, millisecond, nanosecond, and femtosecond lasers allow flexible adaption of laser processes to treated materials. Within ultrashort pluses, the time to heat the material is less than long pulse lasers or continuous lasers. Therefore, this process will result in the highest precision and lowest damage. This way, ultrashort pulsed lasers enable effective ablation through a cold process32,33. However, excessive (long) pulse durations cause more heat to be accumulated at the workpiece, resulting in a thermal process with a significantly larger heat-affected zone31,34. Pulse duration is strongly related to pulse repetition rate. In this study, the effects of pulse duration are investigated from 4 to 200 ns at 20W, 100 mm/s, 1 pass and according to pulse duration-frequency coupling of the laser source presented in Table 2, whereas, the effects of pulse repetition rate are studied at 20W, 100 mm/s, 1pass and 200 ns. In the laser processing of CNT composite, pulse duration has imposed an effect on the penetration depth (Fig. 13a). Still, the range of pulse durations used in the experiments (4–200 ns) imposes very little change on the top width; as seen in Fig. 13b. Figure 14 presents the effect of the repetition rate at 200 ns. In Fig. 14a–c, an increase in penetration depth is observed. Additionally, spatter formation with decreasing trend in its width characterizes the operation. This phenomenon continues until 500 kHz, where it attains the maximum penetration depth in a spatter-free process (Fig. 14d). However, further increasing of repetition rate to 1000 kHz (Fig. 14e) will eventually yield a complete change in the morphology and a drop in the penetration depth. The relationship between penetration depth and repetition rate at 200 ns can be referred from Fig. 13a.

Effect of repetition rate and pulse duration; (a) on penetration depth, (b) on top width, [Fixed: 20W, 100 mm/s, 1 pass].

SEM images showing physical morphology with repetition rate at 200 ns. (a) 20 kHz, (b) 60 kHz, (c) 105 kHz, (d) 500 kHz and (e) 1000 kHz [Fixed: 20W, 100 mm/s, 1 pass and 200 ns].

A close observation of the morphology of five samples in Fig. 14 reveals typical features of the operation. Accordingly, three distinct operational regimes are identified. Regime one consists of cases (a)–(c), 20, 60, and 105 kHz, respectively, which depict operations before the repetition rate of 500 kHz, where a maximum penetration depth is achieved. Operational regime two is identified in case (d) at the 500 kHz. The third operational regime covers the range of repetition rate values above 500 kHz. These categories are well presented in Fig. 15. In regime one, the penetration depth increases with increasing repetition rate. However, the thermal behavior of the process seems to have caused no deformations and as can be confirmed from the physical morphologies of Fig. 15a, b and c. Furthermore, increasing the pulse repetition rate to 500 kHz has resulted in an additional increase in penetration depth (Fig. 13a). Still, close observation of the physical morphology of the sample from Fig. 15d shows that the edges of the groove appear to be thermally damaged. Gradually, unique and interesting characteristics come to an effect as pulse repetition rates are increased beyond 500 kHz. This gives rise to the third operational regime, where the penetration depth drops with an increase in pulse repetition rate. The thermal effect observed for the first time at 500 kHz is advanced to the third regime with more additional effects, as can be observed from Fig. 15e at 1000 kHz. Here, surface layers of the material adjacent to the processed zone are melted, and a recast of these melted layers is formed on its periphery. As discussed above, these three identified regimes have completely unique and different responses to pulse repetition rates causing different ablation characteristics, a gradual increase and decline of penetration depth and morphological features (discovery of thermally damaged and melted regions in regimes two and three). Therefore, it is crucial to conduct further analysis on these three scenarios for a better understanding of the interaction. For this reason, cases of 105, 500, and 1000 kHz are considered representative samples for detailed inspection of the effects on the material's workability.

Operational regime 1—Operation below critical pulse repetition rate

Detailed views of the morphology with repetition rate at 200 ns. (a) 20 kHz, (b) 60 kHz, (c) 105 kHz, (d) 500 kHz and (e) 1000 kHz, [Fixed: 20W, 100 mm/s, 1 pass and 200 ns].

This scenario represents the use of repetition rate below 500 kHz. In this case, penetration depth keeps on increasing with pulse repetition rate, and the process does not result in significant thermal effects. However, thermal effects and changes in the material properties can still occur. For this reason, it is essential to investigate the occurrence of a heat-affected zone (HAZ). HAZ is a region created on a material exposed to elevated temperatures. It is a non-melted zone near the exact work area (a region between the melted material and base material) where microstructural changes occur35. In the operational regime 1, no evidence of surface melting is noticed, reducing the probability of encountering HAZ. To ensure, EDX line scan mode is performed on the surface and cross-section. On the same line, a microhardness test is conducted using Vickers Microhardness Tester. The experiments will give the chemical composition and mechanical property of the specimen in the region where heat could have an impact. The EDX line scan result of the 105 kHz case shown in Fig. 16a indicates that carbon is the dominant element making other elements miscellaneous. The distribution of carbon over a horizontal line of 150 μm from the edge does not result in a significant change in its content but minor random fluctuations. Also, the difference in other miscellaneous elements is null. On the same line, results of Vickers microhardness are presented in Fig. 16c. These local values are compared to the microhardness of unprocessed base material to determine if the thermal effect is causing hardening, softening, or neither. Since the heat from the laser diffuses laterally on the surface and penetrate down the crossection, it is important to observe the surface and the crossection simultaneously. Multiple vickers hardness measurements are conducted at a fixed intervals from the edge of the groove in a configuration given at Fig. 16b. The hardness tests indicated in Fig. 16c show that hardness values decrease as measured from point 3 to point 1; however, the hardness remains unaffected outside this range (point 3–5). The base material hardness was measured at a random point in the unprocessed region, sufficiently far from point 5 (shown in on Fig. 16b). Accordingly, the Vickers microhardness associated with the base material is determined to be \(79.56 \pm 4.2\) HV. It can be noticed that a change in chemical nature and hardness values is taking place. Nearly similar results were found in its cross-section (shown in Fig. 17), where the hardness is found to be 53.1 \(\pm 1.2\) HV. The texture and surface condition of the base material on the cross-section is totally different from the top surface. Consequently, the hardness of unprocessed base material in the cross-sectional face is measured to be 52.57 ± 0.92 HV.

Chemical and hardness tests on surface—105 kHz; (a) EDX line scan (b) measurement configuration (c) Vickers hardness.

Chemical and hardness tests on cross-section—105 kHz.

The carbon concentration on the cross-section and the top surface is different, giving less carbon on the surface. This condition is associated with the formation of a spatter in the regime one. Spatter formation starts in large quantities (higher spatter width) at 20 kHz, decreases in the spatter width up to 105 kHz, and eventually, spatter formation stops at a repetition rate of 500 kHz. The operation from 500 to 1000 kHz is spatter-free. These whole phenomena are indicated in Fig. 18a–e. The EDX elemental mapping and point scan quantifications (averaged) are presented in Fig. 19 and Table 4, respectively. Accordingly, the spatter consists of nearly similar elements with base material but essentially with lower carbon content. In this context, a spatter is a splash of burned and ablated particles from the base CNT composite with small carbon content due to irradiation by high laser pulse energy. Low laser pulses of 500 kHz and above are not observed to yield spatter. Therefore, the underprediction of carbon concentration in the top surface line scan of the 105 kHz case is obviously from the spatter that covers the region.

Spatter distribution with repetition rate at 200 ns. (a) 20 kHz, (b) 60 kHz, (c) 105 kHz, (d) 500 kHz and (e) 1000 kHz [fixed: 20W, 100 mm/s, 1 pass and 200 ns].

EDX mapping of the spatter.

All these results are related to the nature of pulse repetition. In the pulsed nanosecond laser used for the experiments, repetition rate values of 20–105 kHz can be considered low pulse repetition rate values. During low pulse repetition rate settings, the interval between consecutive pulses is relatively long and gives enough time for the dissipation of heat in the material causing the minimum thermal damage. Also, this allows enough time between pulses for the steady removal of ablated particles. For this reason, a significant amount of spatter is observed outside. However, when the repetition rate increases, the time interval between pulses decreases, and the effectiveness of removing ablated particles decreases. Therefore, the spatter amount keeps dropping to 105 kHz and will eventually disappear after that. Due to this, operations involving repetition rate settings below 500 kHz result in proper heat dissipation and effective conveying of ablated particles to the outside so that the changes in the chemical and mechanical (hardness) nature of the work are minimum.

Operational regime 2—Operation at a critical pulse repetition rate

Operational regime two depicts operation at a particular point of 500 kHz where the penetration depth attains its maximum value. At this point, thermal effects start to appear along the edge. The evidence of the thermally damaged region observed in Fig. 15d implies heat accumulation due to the increased number of pulses per unit time. On the other hand, the same figure shows no evidence of trapped ablated particles inside the groove as a result of the increased repetition rate. Important features of this region are obtained by inspection through EDX line scan and microhardness testing according to a configuration given in Fig. 20b. Figures 20a–c and 21 show the tests with the results.

Chemical and hardness tests on surface—500 kHz; (a) EDX line scan (b) measurement configuration (c) Vickers hardness.

Chemical and hardness tests on cross-section—500 kHz.

Carbon content is the most dominant species in the processed and unprocessed region, but carbon content is observed to be higher in the melted region as quantified and presented in Fig. 20a, thereby giving a more softened surface than the previous scenario (shown on Fig. 20c). Looking at the cross-section, the hardness value in the region was measured to be 44.6 \(\pm 1.606\) HV. Comparing this result to the hardness of the base material in the cross-section, \(52.57 \pm 0.92\) HV, it can be noted that softening is going on. Operating at regime two gives the maximum ablation depth. This implies the increase in pulse repetition rate is causing thermal effects from heat accumulation resulting from closely spaced pulses. Still, the thermal effect is not observed to hinder the increasing trend of penetration depth. This situation marks the regime as an operational regime to the onset of thermal effects. Also, this regime gives rise to the termination of spatter formation. For these reasons, 500 kHz is marked as a critical repetition rate.

Operational regime 3—Operation above the critical pulse repetition rate

Unlike regimes one and two, in this operational regime (500–1000 kHz], the penetration depth drops. Also, the regime depicts the presence of special characteristics in its morphology at 1000 kHz. Close observation of Fig. 15e reveals that the thermal effect deforms the edges of the groove. Melting of the surface layer near the processed zone has occurred with a subsequent recasting of the layer in the vicinity of the melting region. Also, larger ablated particles inside the groove are visible closely from the top. Following the same procedure, EDX line scan quantifications of the elemental compositions are carried out on the surface and the cross-section, followed by microhardness testing, as shown in Figs. 22a, b, c and 23.

Chemical and hardness tests on surface—1000 kHz; (a) EDX line scan (b) measurement configuration (c) Vickers hardness.

Chemical and hardness tests on cross-section—1000 kHz.

Like the previous case, the thermally damaged region is mainly composed of Carbon; however, the carbon content in this particular region of melting and recast is the highest. Also, the operation has deteriorated the surface hardness of the material in the same region, where it is highly softened. On the other hand, the same observations on the crossectional side yield no change in the chemical composition, but the softening effect is further penetrated inside as it was measured at 38.5 \(\pm 3.9\) HV near the processed region compared to \(52.57 \pm 0.92\) HV (the base material hardness in the cross-sectional side). As a result, melting with a subsistent formation of recast accompanies the ablation process at the third operational regime. The improper heat accumulation process significantly alters the region's chemical and mechanical nature, resulting in a softened edge with more carbon content.

Repetition rate is a measure of the number of pulses emitted per second of a regular train of pulses36. The repetition rate is an essential parameter that defines heat utilization for the laser processing of certain materials37. Operating at high repetition rate values has its advantages; for instance, in micro-hole drilling using an ultra-pulse laser, higher efficiencies are expected at higher repetition rate since more pulses are radiated at a fixed time38. However, excessive repetition rate values cause adverse effects on materials' processability, and the relative magnitudes depend on the material type37. In laser processing with high repetition rate values (> 500 kHz), the duration between individual pulses is very short in such a way that thermal energy cannot be adequately dissipated from the heated volume before the arrival of the next pulse39. Apart from pulse repetition rate, previous studies indicated that pulse energy and peak power are essential parameters that control the depth and width of a microgroove34,35. The pulse energy and peak power values for each repetition rate of 20–1000 kHz are calculated based on Eqs. (1) and (2). The results are plotted and presented in Fig. 24.

where Pavg, f, and \(\tau\) are average power, repetation rate, and pulse duration, respectively.

Pulse Energy and Peak Power versus Repetition Rate.

Comparing the results of Fig. 24 with the results of Fig. 13a and b, it is quite evident that the effects of pulse energy and peak power are highly associated with top width. High peak power and pulse energy give larger top width and vice versa. In addition, from 20 to 105 kHz, the top width decreases to a great extent. This is from the sharp reduction of pulse energy and peak power for the same range of repetition rate values. However, the penetration depth increases while the pulse energy decreases from 1–0.04 mJ, and the penetration depth decreases while pulse energy decreases from 0.04 to 0.02 mJ. Also, a similar trend is observed with peak power. This situation is attributed to the two dominant effects of a high pulse repetition rate 40,41,42,43. The effects are associated with the shortening of the time interval between two successive pulses. This causes: (i) less time for heat diffusion in the material, which introduces a heat accumulation problem in the processed region. Accumulation of excessive heat in the material will not have a contribution to the material removal process but to the formation of morphological changes through the formation of molten spatters and melted regions near the periphery of the processed zone with the subsequent formation of the recast. (ii) less time to remove ablated particles from the processed zone. When the ablated particles are trapped inside the groove due to the shorter time between pulses for effective removal, the upcoming pulses will get shielded. This leads to the melting of these particles. Later, the resolidification of these particles will affect the penetration depth.

Generally, the characteristics observed in the interaction of the laser with the CNT composite can be mapped from the standpoint of laser peak intensity or peak pulse intensity. Laser peak intensity is defined as the maximum power the incident beam delivers per unit area44. Laser peak intensity is another crucial working parameter determining ablation rates and other physical phenomena associated with the process45. The laser peak intensity is calculated for a single pulse based on Eq. (3). Since laser peak intensity is a function of average power (\({P}_{avg}\)), pulse duration (\(\tau\)), and pulse repetition rate (\(f\)), it could give a more generalized indications on the interaction.

where Asp is the laser spot area [mm2].

Figure 25 shows the penetration depth and top width versus peak pulse intensity for the parameters studied from 50 to 200 ns. According to the results, by varying power from 4 to 20W at 20 kHz, peak intensity would be varied from 1413.7 to 7073.6.409 kW/mm2 to yield more penetration without causing any additional significant morphological changes. Figure 25 illustrates key features identified by the effect of various laser parameters. Accordingly, heat accumulation regimes discussed earlier cover a peak intensity range of 141.37–282.74 kW/mm2. Any parametric combination at a low pulse repetition rate and low number of passes is traced in the range of peak intensity range 1413.7–28,294.4 kW/mm2 and identified as a proper working regime. However, the intense peak intensities resulting from the higher number of irradiation due to increasing the number of passes (after 12 passes) results in degradation of the morphological appearance. A range of extremely high peak intensity of 845,883.2–141,472 kW/mm2 denotes these operations, suggesting that these high peak intensities introduce undesired effects to the material. Generally, peak intensities ranging from 1413.7 to 28,294.4 kW/mm2 can be favored for effective and good-quality operations.

Mapping of laser interaction according of peak intensity.

Similarly, the laser peak pulse intensities are mapped with the laser parametric combinations to see the effect of laser peak pulse intensity on the top width. It is interesting that the same pulse intensity range (1413.7–28,294.4 kW/mm2) provides the proper working regime in this case as well. In the case of the spatter (Fig. 26), it happens while working at 1413.7–141472 kW/mm2 peak intensity range on any operation that lies in operational regime one. Operational regimes two and three are already known to give spatter-free operation and cover a very narrow peak intensity range of 141.37–282.74 kW/mm2.

Mapping of laser interaction according of peak intensity and spatter.

In the laser processing of polymeric composites, intense laser beams tend to break the polymer chains and introduce process defects and failures in the material46. Incorporating CNTs in composites plays a vital role in their processability by increasing the hardness and laser absorptivity of composites47. It has been reported that laser processing of carbon-based composites suffers from thermal effects due to large differences in the thermal conductivity and vaporization temperature between the carbon and the matrix48. Besides, the latent heat capacity of carbon is 43 times that of the epoxy resin, implying that carbon will absorb most of the heat49. The absorbed thermal energy will then be conducted along carbon fibers, and the epoxy resin will evaporate by the heat even in the non-processed region. This situation marks the formation of HAZ, and material removal became complicated48. This study demonstrated that significantly high laser beam intensities could be employed for proper material removal (Fig. 25), which emphasizes the role of CNTs in absorbing part of the thermal energy input and providing additional pathways50 for effective heat conduction. This mechanism is believed to significantly reduce the thermal effects on laser processability.

This study has presented an experimental parametric study on laser interaction with CNT composite plate for use as bipolar plate in PEMFCs by considering a 2.5 mm thick plate. The basic concluding remarks outlined by the study are:

The effect from the number of passes is the highest in controlling penetration depth and top width. Additionally, care must be taken to render a high value of the number of passes as it is indicated to be a negative factor for quality operation.

Relatively smooth operations resulting from power and scanning speed indicate that these parameters are adequate for low to medium control of penetration depth and top width without significant impact on morphological aspects.

Pulse repetition rate presents a unique nature of interaction that resulted in a critical repetition rate distinguishing three different regimes: below the critical value, at the critical value, and above the critical value, each with certain distinct chemical and mechanical behaviors.

Operations that favored repetition rate values below the critical value yielded reliable results.

Typical channel dimensions of flow channels on bipolar plates for the rectangular slot were reported as 300 μm width and 300 μm height51. Overall, the ablation characteristics suggest that laser processing can be considered a potential method for manufacturing flow channels of bipolar plates using CNT composite materials.

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Esawi, A. M. K. & Farag, M. M. Carbon nanotube reinforced composites: Potential and current challenges. Mater. Des. 28(9), 2394–2401 (2007).

Article CAS Google Scholar

Gorrasi, G. et al. Effect of temperature and morphology on the electrical properties of PET/conductive nanofillers composites. Compos. Part B Eng. 135, 149–154 (2018).

Article CAS Google Scholar

Venkataraman, A., Amadi, E. V., Chen, Y. & Papadopoulos, C. Carbon nanotube assembly and integration for applications. Nanoscale Res. Lett. https://doi.org/10.1186/s11671-019-3046-3 (2019).

Article Google Scholar

Endo, M., Hayashi, T., Kim, Y. A. & Muramatsu, H. Development and application of carbon nanotubes. Jpn. J. Appl. Phys. Part 1 Regul. Pap. Short Notes Rev. Pap. (2006).

Cesano, F., Uddin, M. J., Damin, A. & Scarano, D. Multifunctional conductive paths obtained by laser processing of non-conductive carbon nanotube/polypropylene composites. Nanomaterials 11, 604 (2021).

Article CAS Google Scholar

Ola, O., Niu, Q., Chen, Y., Xia, Y. & Zhu, Y. Carbon nanotube reinforced nanocomposites for energy conversion and storage. J. Power Sources 443, 227277 (2019).

Article CAS Google Scholar

Hu, S. & Hu, C. Carbon nanotube-based electrochemical sensors: Principles and applications in biomedical systems. J. Sensors https://doi.org/10.1155/2009/187615 (2009).

Article Google Scholar

Abdalla, S., Al-Marzouki, F., Al-Ghamdi, A. A. & Abdel-Daiem, A. Different technical applications of carbon nanotubes. Nanoscale Res. Lett. 10, 1–10 (2015).

Article CAS Google Scholar

Tan, C. W. et al. Carbon nanotubes applications: Solar and fuel cells, hydrogen storage, lithium batteries, supercapacitors, nanocomposites, gas, pathogens, dyes, heavy metals and pesticides. In Environmental Chemistry for a Sustainable World (Springer, 2012).

Google Scholar

Dai, L., Chang, D. W., Baek, J. B. & Lu, W. Carbon nanomaterials for advanced energy conversion and storage. Small 8, 1130–1166 (2012).

Article ADS CAS Google Scholar

Luo, C., Xie, H., Wang, Q., Luo, G. & Liu, C. A review of the application and performance of carbon nanotubes in fuel cells. J. Nanomater. https://doi.org/10.1155/2015/560392 (2015).

Article Google Scholar

Britto, P. J., Santhanam, K. S. V., Rubio, A., Alonso, J. A. & Ajayan, P. M. Improved charge transfer at carbon nanotube electrodes. Adv. Mater. 11, 154–157 (1999).

3.0.CO;2-B" data-track-action="article reference" href="https://doi.org/10.1002%2F%28SICI%291521-4095%28199902%2911%3A2%3C154%3A%3AAID-ADMA154%3E3.0.CO%3B2-B" aria-label="Article reference 12" data-doi="10.1002/(SICI)1521-4095(199902)11:23.0.CO;2-B">Article CAS Google Scholar

Shojaeifar, M., Askari, M. B., Hashemi, S. R., Di Bartolomeo, A. MnO 2-NiO-MWCNTs nanocomposite as a catalyst for methanol and ethanol electrooxidation. J. Phys. D. Appl. Phys. 55 (2022).

Li, L., Yang, H., Zhou, D. & Zhou, Y. Progress in application of CNTs in lithium-ion batteries. J. Nanomater. 17, 493–516 (2014).

Google Scholar

Zhou, Y. & Azumi, R. Carbon nanotube based transparent conductive films: progress, challenges, and perspectives. Sci. Technol. Adv. Mater. 17, 493–516 (2016).

Article CAS Google Scholar

Liang, X., Xia, J., Dong, G., Tian, B. & Peng, I. Carbon nanotube thin film transistors for flat panel display application. In Topics in Current Chemistry (Springer, 2016).

Google Scholar

Ma, J. et al. Non-traditional processing of carbon nanotubes: A review. Alex. Eng. J. 61, 597–617 (2022).

Article Google Scholar

Mahmoodi, M., Arjmand, M., Sundararaj, U. & Park, S. The electrical conductivity and electromagnetic interference shielding of injection molded multi-walled carbon nanotube/polystyrene composites. Carbon 50, 1455–1464 (2012).

Article CAS Google Scholar

Chang-Jian, S.-K. & Ho, J.-R. Laser patterning of carbon-nanotubes thin films and their applications. In Carbon Nanotubes Applications on Electron Devices (IntechOpen, 2011).

Google Scholar

Davim, J. P., Barricas, N., Conceição, M. & Oliveira, C. Some experimental studies on CO2 laser cutting quality of polymeric materials. J. Mater. Process. Technol. 198, 99–104 (2008).

Article CAS Google Scholar

Choudhury, I. A. & Shirley, S. Laser cutting of polymeric materials: An experimental investigation. Opt. Laser Technol. 42, 503–508 (2010).

Article ADS CAS Google Scholar

Yuan, X. Z., Wang, H., Zhang, J. & Wilkinson, D. P. Bipolar plates for PEM fuel cells - From materials to processing. J. New Mater. Electrochem. Syst. 8, 257 (2005).

CAS Google Scholar

Wu, X., Yin, H. & Li, Q. Ablation and patterning of carbon nanotube film by femtosecond laser irradiation. Appl. Sci. 9, 3045 (2019).

Article CAS Google Scholar

Kichambare, P. D. et al. Laser irradiation of carbon nanotubes. Mater. Chem. Phys. 72, 218–222 (2001).

Article CAS Google Scholar

Karimzad Ghavidel, A., Azdast, T., Shabgard, M. R., Navidfar, A. & Mamaghani Shishavan, S. Effect of carbon nanotubes on laser cutting of multi-walled carbon nanotubes/poly methyl methacrylate nanocomposites. Opt. Laser Technol. 67, 119–124 (2015).

Article ADS CAS Google Scholar

Gautam, R. K., Banerjee, S. & Kar, K. K. Bipolar plate materials for proton exchange membrane fuel cell application. Recent Patents Mater. Sci. 8, 15–45 (2015).

Article CAS Google Scholar

Wilberforce, T. et al. Numerical modelling and CFD simulation of a polymer electrolyte membrane (PEM) fuel cell flow channel using an open pore cellular foam material. Sci. Total Environ. 678, 728–740 (2019).

Article ADS CAS Google Scholar

Kakati, B. K., Sathiyamoorthy, D. & Verma, A. Electrochemical and mechanical behavior of carbon composite bipolar plate for fuel cell. Int. J. Hydrog. Energy 35, 4185–4194 (2010).

Article CAS Google Scholar

Jang, J. H., Joo, B. D., Mun, S. M., Van Tyne, C. J. & Moon, Y. H. Application of a direct laser melting process for fabricating a micropattern in bipolar plates used in direct methanol fuel cells. Proc. Inst. Mech. Eng. Part B J. Eng. Manuf. 225, 1784–1791 (2011).

Article CAS Google Scholar

Roth, B., Søndergaard, R. R., & Krebs, F. C. Handbook of Flexible Organic Electronics. (Woodhead Publishing, 2015).

Hashmi, M. Saleem J. Comprehensive Materials Processing (Newnes, 2014).

Weber, R. et al. Processing constraints resulting from heat accumulation during pulsed and repetitive laser materials processing. Opt. Express 2, 3966 (2017).

Article ADS Google Scholar

Phillips, K. C., Gandhi, H. H., Mazur, E. & Sundaram, S. K. Ultrafast laser processing of materials: A review. Adv. Opt. Photonics 7, 684–712 (2015).

Article ADS CAS Google Scholar

Hamad, A. H. Effects of different laser pulse regimes (Nanosecond, Picosecond and Femtosecond) on the ablation of materials for production of nanoparticles in liquid solution. In High Energy and Short Pulse Lasers (InTech, 2016).

Google Scholar

Agrawal, R. & Wang, C. Laser beam machining. In Encyclopedia of Nanotechnology (ed. Bhushan, B.) 1739–1753 (Springer Netherlands, 2016).

Chapter Google Scholar

Weik, M. H. Pulse repetition rate. In Computer Science and Communications Dictionary (2000).

Weber, R. et al. Heat accumulation during pulsed laser materials processing. Opt. Express 22, 11312–11324 (2014).

Article ADS CAS Google Scholar

Gruner, A., Naumann, L., Schille, J. & Loeschner, U. High pulse repetition frequency micro hole drilling of silicon using ultrashort pulse laser radiation. J. Laser Micro Nanoeng. 14, 190–197 (2019).

CAS Google Scholar

Vainos, N. A. Laser Growth and Processing of Photonic Devices (Elsevier, 2012).

Book Google Scholar

Zheng, Q. et al. Mechanism and morphology control of underwater femtosecond laser microgrooving of silicon carbide ceramics. Opt. Express 27, 26264–26280 (2019).

Article ADS CAS Google Scholar

Chang, F.-Y., Hsu, C.-F. & Lu, W.-H. Nanosecond-fiber laser cutting and finishing process for manufacturing polycrystalline diamond-cutting tool blanks. Appl. Sci. 11, 5871 (2021).

Article CAS Google Scholar

Zhao, W., Shen, X., Liu, H., Wang, L. & Jiang, H. Effect of high repetition rate on dimension and morphology of micro-hole drilled in metals by picosecond ultra-short pulse laser. Opt. Lasers Eng. 124, 105811 (2020).

Article Google Scholar

Mustafa, H., Matthews, D. T. A. & Römer, G. R. B. E. The role of pulse repetition rate on picosecond pulsed laser processing of Zn and Zn-coated steel. Opt. Laser Technol. 131, 106408 (2020).

Article CAS Google Scholar

Smeenk, C. et al. Precise in-situ measurement of laser pulse intensity using strong field ionization. Opt. Express 19(10), 9336–9344 (2011).

Article ADS CAS Google Scholar

Nath, A. K. Laser drilling of metallic and nonmetallic substrates. In Comprehensive Materials Processing (2014).

Kannatey-Asibu, E. Jr. Principles of Laser Materials Processing (John Wiley & Sons, 2009).

Book Google Scholar

Wan, L. et al. Directed energy deposition of CNTs/AlSi10Mg nanocomposites: Powder preparation, temperature field, forming, and properties. Opt. Laser Technol. 139, 106984 (2021).

Article CAS Google Scholar

Tao, N. et al. Temperature-dependent material removal during pulsed laser processing of CFRP composites. Opt. Laser Technol. 144, 107445 (2021).

Article CAS Google Scholar

Miyamoto, I., Cvecek, K. & Schmidt, M. Evaluation of nonlinear absorptivity and absorption region in fusion welding of glass using ultrashort laser pulses. Phys. Procedia 12, 302–310 (2011).

Article Google Scholar

Grundler, M., Derieth, T., Beckhaus, P., Heinzel, A. & Cell, F. CarbonNanoTubes (CNT ) in bipolar plates for PEM fuel cell applications. In Proceedings of the 18th World Hydrogen Energy Conference (2010).

Bhattacharyya, B. Electrochemical Micromachining for Nanofabrication, MEMS and Nanotechnology (Elsevier, 2015).

Google Scholar

Download references

The research described herein was sponsored by the National Research Foundation of Korea (NRF) grant (No. 2021R1C1C1008671) funded by the Ministry of Science and ICT(MSIT, Korea), sponsored by the Regional Innovation Strategy (RIS) through the National Research Foundation of Korea (NRF) (2021RIS-004) funded by the Ministry of Education(MOE), and sponsored by the Enabling Technology Development Project funded by the Cheonan Institute of Science and Technology Platform (CISTEP, Cheonan, Korea), This work is also supported by the Korea Institute for Advancement of Technology (KIAT) (P0018009) funded by the Ministry of Trade, Industry and Energy(MOTIE, Korea).

Department of Future Convergence Engineering, Cheonan College of Engineering, Kongju National University, Cheonan, 31080, South Korea

Dawit Musse & Dongkyoung Lee

Department of Mechanical and Automotive Engineering, Cheonan College of Engineering, Kongju National University, Cheonan, 31080, South Korea

Dongkyoung Lee

Center for Advanced Powder Materials and Parts of Powder (CAMP2), Cheonan College of Engineering, Kongju National University, Cheonan, 31080, South Korea

Dongkyoung Lee

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

D.M.: Conceptualization, methodology, experimental preparation, formal analysis, investigation, writing original draft, writing review, and editing. D.K.L.: Supervision, administration, funding organization, paper review.

Correspondence to Dongkyoung Lee.

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/.

Reprints and Permissions

Musse, D., Lee, D. Parametric study of pulsed nanosecond laser interaction with carbon-nanotube composite bipolar plate for PEMFCs. Sci Rep 13, 2048 (2023). https://doi.org/10.1038/s41598-023-28700-2

Download citation

Received: 04 November 2022

Accepted: 23 January 2023

Published: 04 February 2023

DOI: https://doi.org/10.1038/s41598-023-28700-2

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

Scientific Reports (2023)

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.

SHARE