Study of Diffusible Behavior of Hydrogen in First Generation Advanced High Strength Steels

: This study deals with microstructural inﬂuence on the H permeation behavior of Dual-Phase (DP) and Complex Phase (CP) steels using electrochemical permeation studies. The H diffusion coefﬁcients in DP steels (DP800: 1.65 × 10 − 10 m 2 · s − 1 , DP1000: 1.58 × 10 − 10 m 2 · s − 1 ) are half of that found in CP steels (3.07 × 10 − 10 m 2 · s − 1 ).The banded microstructure along the specimen thickness and higher C content of the DP led to high H diffusivity of DP steels. The lower total H concentration along with a higher fraction of H was present in the stronger traps in CP steels suggest a better HE resistance of this steel. The H distribution in the specimens was non-uniform, with a higher H concentration speculated near the charging surface.


Introduction
Reducing vehicle weight is a crucial challenge for the automobile industry to increase fuel efficiency and reduce vehicular emission without compromising safety. Advanced High Strength Steels (AHSS), which offer superior strength, are therefore are gaining attention as the body in white material (BIW), with its use increasing from 7% in 2010 to 21% in 2017 (and projected to increase up to 40% by 2025) [1]. These materials are also under interest in other fields such as pipeline, where the corrosion performances versus the chemical composition of steels are investigated in various environments to develop predictive models and to prevent further failures [2,3]. In the latter case, bainite growth in the weld part improved corrosion resistance [4]. Metallurgy of the AHSS is a key parameter that must be controlled along the processing route since it controls the mechanical properties and corrosion resistance of engineering structures [5][6][7][8].
The AHSS are also known to be sensitive to hydrogen embrittlement (HE) and continues to pose an essential bottleneck for their utilization [9][10][11][12][13][14][15]. Hydrogen present in steels, even to an extent of less than 1 ppm, diffuses towards triaxially stressed regions, and thereby cause premature failures. Such failure can occur instantaneously or delayed, making failure prediction difficult, increasing liability, and repair costs [16][17][18]. Since H diffusion directly impacts H transport to stressed areas, studying H permeation behavior in steels is essential in understanding the overall embrittlement behavior [19][20][21][22].

Experimental
For the present study, TATA Steel (Oosterhout, Netherlands) provided two DP steels having a tensile strength of about 800 MPa (DP800) and 1000 MPa (DP1000) and a CP steel having a tensile strength of about 1000 MPa (CP1000) in the form of 30 cm × 30 cm sheet of different thicknesses. The DP steels had an average thickness of 0.98 mm, while the CP steel had 1.60 mm average thickness. The chemical compositions of these steels are shown in Table 1.

Sample Preparation for Metallography
For microstructural evaluation, two 15 mm × 20 mm coupons (one along the specimen surface and the other along the cross-section) from each steel were ground up to 4000 SiC emery paper and then polished up to 1 µm diamond paste. Such selection helped to verify the existence of banding behavior along the rolling direction. Subsequently, they were washed twice with acetone and then with ethanol in an ultrasonic bath and air-dried. The polished surfaces were etched with a 2% Nital solution and observed under a Jeol JSM-7800F PRIME (Jeol, Akishima, Japan) Scanning Electron Microscope (SEM).

Hydrogen Permeation Specimen Preparation
35 mm × 35 mm coupons were ground up to 600 grit SiC, cleaned twice in an ultrasonic bath with acetone and then with ethanol and air-dried. Pd layer was electrodeposited on one surface (to avoid Fe oxidation at the anodic surface) while other surfaces were masked with a conductive Al tape to avoid Pd deposition. The specimen was immersed in a Pd bath (1L solution: 2.5 g PdCl 2 , 120 g Na 3 PO 4 , 27.5 g NH 4 PO 4 , 1.7 g C 6 H 5 COOH) and a current density of 2.8 mA·cm −2 was applied for 180 s at room temperature, with a Pt mesh as an anode. ChimiePlus (Saint-Paul-de-Varax, France) supplied all chemicals. After Pd deposition, the tape was removed, both specimen surfaces were cleaned with alcohol and air-dried.

Hydrogen Permeation Cell
A modified Devanathan and Stachurski cell (D.S. cell or permeation cell), consisting of two compartments, namely charging and detection, separated by a test specimen (working electrode, i.e., W.E.), as outlined in the schematic of Figure 1 was used. Each compartment was equipped with an auxiliary electrode (A.E.: graphite at the detection side and a Pt mesh at the charging side) and a saturated Ag/AgCl reference electrode (R.E.). Two independent Gamry Interface 1000 potentiostats, one to electrochemically reduce H + at the cathodic/charging side and the other to subsequently oxidize H at the detection/anodic side, were used. A fraction of the H + reduced at the cathode surface of the specimen was absorbed as H into the specimen, diffused across the specimen, and finally oxidized at the anode surface (Pd coated) of the specimen. In this study, the H + was reduced galvanostatically with varying current densities and the oxidized potentiostatically at +0.2 V w.r.t. RE. The potential represented an overpotential of +250 mV, ensuring complete ionization of H at the detection side as per H ads → H + + e − . Heated water circulated in the double jacket to maintain a constant temperature of 35 • C at each compartment.

Experimental Procedure
The electrochemical hydrogen permeation experiment was carried out according to the ISO 17081 specifications [10,46,54]. Specimens with an exposed surface area of 7 cm 2 were clamped between the two compartments. The anodic compartment was filled with

Experimental Procedure
The electrochemical hydrogen permeation experiment was carried out according to the ISO 17081 specifications [10,46,54]. Specimens with an exposed surface area of 7 cm 2 were clamped between the two compartments. The anodic compartment was filled with 0.2 M KOH solution (pH = 13; Conductivity = 45.3 mS·cm −1 ) and anodic potential (+0.2 V w.r.t. R.E.) was applied. Two Ref1000 potentiostat (GAMRY, Warminster, PA, USA) were used to apply and record the current and potentials for each cell. The results were measured using Gamry Framework data acquisition software v6.32. Once the measured current density stabilized below 0.3 µA·cm −2 , de-aerated 0.1 M K 2 SO 4 , 0.02 M H 2 SO 4 solution (pH = 2.08; Conductivity = 26.5 mS·cm −1 ) was added slowly at the charging side and the cathodic surface was charged galvanostatically. Ar gas flowed through both the electrolytes throughout the experiment for de-aerating the electrolytes. In the absence of oxygen, H oxidation is the only electrochemical reaction at the anodic surface and provides a direct measure of the H atoms permeating through the specimen.
The galvanostatic charging at the cathodic surface involved the following three steps. As an example, Figure 2a represents the permeation current density as obtained due to the charging schedule at the cathodic surface for CP1000. Similar testing was conducted for other specimens.

•
Step I (Cathodic Pre-charging): When the background current density at the detection side was below 0.3 µA·cm −2 , a cathodic current density of −10 mA·cm −2 was applied on the charging surface until an observable steady-state current permeation current density occurred (Figure 2b).

•
Step II (Subsequent Partial Transient): After step I, the charging current density was raised to −20 mA·cm −2 . Once the permeation current density reached a new steadystate, the charging current density was lowered back to −10 mA·cm −2 . The applied current density was held for 5000 s for each step. The steps were repeated to verify the permeation curves' reproducibility (as shown in Figure 2c).

•
Step III (Partial Transient Loop): After completing step II, the charging current density was raised from −10 mA·cm −2 to −50 mA·cm −2 and after that lowered to −10 mA·cm −2 (thus creating a loop), with a step size of −10 mA·cm −2 ( Figure 2d) and holding time of 5000 s for each step.
The above charging schedule reduces the surface modification influence in the permeating current density occurring at the charging surface, which is one of the key challenges in applying ISO 17081. Additionally, as a cathodic potential is always associated with the charging surface, it minimizes the possibility of corrosion issues, particularly during fall transients [19,27,55,56]. It is reported that surface modifications (dissolution/formation of oxides/hydrides/reaction products) occurring during H charging can significantly facilitate/block hydrogen entry in the material during the test and thus influence the observed permeation behavior [57][58][59][60][61]. According to Turnbull, during the early stages of the test, there will always be a partial surface reaction, and its effect could be minimized by constant charging [54]. A long-term cathodic pre-charging, as carried out in step I, will, therefore, stabilize the specimen surface and minimize the influence of surface modification on the permeation curves [62][63][64].
To verify pre-charging's effectiveness in reducing surface modification's influence on permeation current density, the steady-state permeation current density was plotted against the square root of applied cathodic current densities for pre-charged specimens, as shown in Figure 3. The linear relationship between the two parameters validate that the permeation current density is indeed controlled by the impressed cathodic current density and not influenced by surface modifications and is attributed to the activation of cathodic surface facilitating hydrogen absorption [54,65]. Zakroczymski et al. [55] compared this relation in non-pre-charged and pre-charged specimens and observed a sigmoidal curve in case of non-pre-charged specimens, which followed a linear relationship in the case of pre-charged specimens. the permeation current density is indeed controlled by the impressed cathodic current density and not influenced by surface modifications and is attributed to the activation of cathodic surface facilitating hydrogen absorption [54,65]. Zakroczymski et al. [55] compared this relation in non-pre-charged and pre-charged specimens and observed a sigmoidal curve in case of non-pre-charged specimens, which followed a linear relationship in the case of pre-charged specimens.  Sq. root of charging current density (mA.cm -2 ) 0.5 Figure 3. The variation of steady-state permeation current density with the square root of charging current density. A linear relationship is assumed, suggesting a permeation current controlled by the impressed cathodic current density.
Step II followed step I, where the partial rise and fall transients' reproducibility were analyzed, as shown in Figure 4a,b. Both curves overlap each other, validating their reproducibility and increasing the reliability of the data obtained from the subsequent steps.
Steps I and II are the surface conditioning steps, whereas step III would be used to provide the permeation data. The subsequent partial transients, as shown in Figure 4a,b, were in good agreement with each other. They attributed the good reproducibility to (a) pre- Step II followed step I, where the partial rise and fall transients' reproducibility were analyzed, as shown in Figure 4a,b. Both curves overlap each other, validating their reproducibility and increasing the reliability of the data obtained from the subsequent steps. Steps I and II are the surface conditioning steps, whereas step III would be used to provide the permeation data. The subsequent partial transients, as shown in Figure 4a,b, were in good agreement with each other. They attributed the good reproducibility to (a) pre-charging, filling up all the irreversible traps leading to the interaction of hydrogen with only the reversible traps and lattice sites after a long pre-charging; (b) the trapping effect from the reversible traps to be similar for both transients; and (c) the surface effects for both the transients to be similar [57,59,66]. The reproducibility of the curves also indicates that there was no permanent damaged caused to the material (also confirmed by observing the specimen under an optical microscope) due to high current density of the cathodic charging, a concern which has been raised by some authors in the past [56,67]. current density. A linear relationship is assumed, suggesting a permeation current controlled by the impressed cathodic current density.
Step II followed step I, where the partial rise and fall transients' reproducibility were analyzed, as shown in Figure 4a,b. Both curves overlap each other, validating their reproducibility and increasing the reliability of the data obtained from the subsequent steps.
Steps I and II are the surface conditioning steps, whereas step III would be used to provide the permeation data. The subsequent partial transients, as shown in Figure 4a,b, were in good agreement with each other. They attributed the good reproducibility to (a) precharging, filling up all the irreversible traps leading to the interaction of hydrogen with only the reversible traps and lattice sites after a long pre-charging; (b) the trapping effect from the reversible traps to be similar for both transients; and (c) the surface effects for both the transients to be similar [57,59,66]. The reproducibility of the curves also indicates that there was no permanent damaged caused to the material (also confirmed by observing the specimen under an optical microscope) due to high current density of the cathodic charging, a concern which has been raised by some authors in the past [56,67].  Figure 4. The normalized curve for two successive (a) rise transients for CP1000 steel when the charging current density was raised from −10 mA/cm 2 to −20 mA/cm 2 and (b) fall transients for CP1000 steel when the charging current density was lowered from −20 mA/cm 2 to −10 mA/cm 2 . Figure 5 presents the scanning electron images (SEM) of the three alloys investigated in this study. DP and CP steels cross section images (right column) reveal if there exits any banding behavior along the rolling direction. Such microstructural anisotropy is expected to play a role in the diffusion mechanism of hydrogen (diffusion coefficient, trapping concentration). The phase fractions and the average grain sizes were determined using EBSD and the same are presented along with the mechanical properties in Table 2. Figure 5a,c presents the rolled-surface microstructures of DP800 and DP1000 specimens, respectively that reveal the martensitic islands and ferrite [68]. The high martensite phase fraction in the alloys makes the ferrite discontinuous, an observation similar to the works of Davies et al. [69,70]. Between the two DP steels, DP1000 contains a higher martensite phase fraction than DP800, which is responsible for the higher strength and reduced ductility of DP1000. With a higher martensitic phase fraction and similar size of the martensitic islands in DP1000, the average ferritic grain size in DP800 was significantly larger than DP1000. The higher ferritic grain size can also lower the tensile strength of an alloy. Metals 2021, 11, x FOR PEER REVIEW 8 of 22 Figure 5. SEM images of the surfaces of (a) DP800 steel (c) DP1000 steel, (e) CP1000 steel, and cross-sections of (b) DP800 (d) DP1000, and (f) CP1000 steel sheets. The orange, blue and green arrows represent ferrite, martensite, and bainite, respectively.

Hydrogen Diffusion Coefficient
H diffusion in metals is a critical parameter to understand their sensitivity to HE as it is the diffusible H that migrates to the crack tip, causing embrittlement. The hydrogen diffusion coefficient was calculated by fitting the experimentally obtained partial permeation transient data to Equations (1) and (2), assuming that during H charging, the H concentration at the charging side attains a new value just after the cathodic current density is varied [10,26,30].   Figure 5e presents the as rolled-surface microstructure of the CP1000 specimen, revealing martensite as the predominant phase followed by ferrite and bainite. The CP steel contains significantly higher quantities of bainite than both the DP steels. Further, the average size of the ferrite grains and the martensitic islands found in CP1000 is similar but is smaller than those found in DP steels. The phase fraction of the retained austenite in all the investigated steels was below 0.25%. Figure 5b,d,f present the cross-sectional images obtained along the rolling direction for all the steel sheets under study. Notably, both the DP800 and DP1000 have banded structures aligned along the rolling direction, whereas CP1000 is free from it.

Hydrogen Diffusion Coefficient
H diffusion in metals is a critical parameter to understand their sensitivity to HE as it is the diffusible H that migrates to the crack tip, causing embrittlement. The hydrogen diffusion coefficient was calculated by fitting the experimentally obtained partial permeation transient data to Equations (1) and (2), assuming that during H charging, the H concentration at the charging side attains a new value just after the cathodic current density is varied [10,26,30].
where, i n is the normalized current density (the initial and final current densities for the partial rise and partial fall transients were 0), D (m 2 ·s −1 ) is the H diffusion coefficient, L is the specimen thickness (m), and t is the time (s) for each partial transient. The H diffusion coefficient so-measured corresponds to the apparent diffusion coefficient (D app ) and not to the lattice diffusion (D L ) as traps present in the alloy impede the hydrogen movement [27]. The experimental data fitted to the above equations for the various initial and final charging test conditions for DP800 is presented in Figure 6 for illustration. The goodness of fit of the curves suggests that the hydrogen diffusion is not affected by the slow surface processes [26,27,59,64]. A similar fitting procedure was carried out for the other steels as well.
The diffusion coefficient values obtained in this work for the tested steel specimens are within the range of diffusion coefficients reported in the literature [56,57,71]. The diffusion coefficient values obtained for DP800, DP1000 and CP1000 are in the range: 0.83 to 1.68 × 10 −10 m 2 ·s −1 ; 0.77 to 1.60 × 10 −10 m 2 ·s −1 ; and 1.07 to 3.16 × 10 −10 m 2 ·s −1 respectively as presented in Figures 7 and 8. Figure 8 presents the variation of D app with the charging current density for all three steel specimens. The increase in D app with increasing charging current density or cathodic overpotential is consistent with the literature [26,56,59,64]. The above increase in D app results from the incremental filling of traps and the consequent reduction of the trapping capacity of microstructural features with a step-wise increase in cathodic overpotential/charging current density. Another factor that might contribute to the rising D app value with the applied current density is the increase in the H concentration gradient across the membrane facilitating the H diffusion across the membrane.
Similarly, the D app value decreases with the subsequent reduction in the charging current density. With the reduction in applied current density, H is released from both the lattice sites and the trap sites, maintaining equilibrium [46]. It results in the increasing influence of trapping sites impeding the H movement and consequently decreasing the D app value. With the evident decreasing influence of the traps with increase in charging current density/cathodic overpotential, the lattice H diffusion coefficient (D L ) in the investigated steels was calculated from the data obtained at the highest overpotential performed in the study [10,30,59]. The D L values obtained for DP800, DP1000, and CP1000 are 1.65 × 10 −10 m 2 ·s −1 , 1.58 × 10 −10 m 2 ·s −1 , and 3.07 × 10 −10 m 2 ·s −1 respectively. The so-obtained D L values are assumed to be a close approximation of the actual lattice H diffusion coefficient in the investigated steels. Data shows that CP1000 has the highest D L value, while both the DP steels show nearly the same D L values. As expected, the D L is significantly higher than the corresponding D app for all the observed cases. Zakroczymski et al. [27] made similar observations, where the diffusion coefficient of H in non pre-charged pure iron was an order of magnitude smaller than the lattice diffusion coefficient.

Hydrogen Concentration
The sub-surface and lattice hydrogen concentrations are the other important parameters contributing to HE of steels with only a limited information is available in the published literature. Using the steady-state current density and the calculated diffusion coefficient values, both the total subsurface H concentration (C T ) and the lattice H concentration (C L ) in mol·m −3 , are determined. Under equilibrium conditions, C T and C L are related to the steady-state current density (i ss ) via the corresponding diffusion coefficients using Equations (1) and (4), respectively [19,72].
where F is the Faraday's constant.
The difference between the C T and C L corresponds to the H present in the traps (C Trapped ). The above three H concentration values are plotted against the charging current density in Figure 9, and the corresponding data is presented in Table 3. In general, all the three H concentration values increase with increase in the applied current density. The data suggest that C T and C Trapped for the steels vary in the order DP1000 > DP800 > CP1000, while a cross-over, from the above trend, occurs for DP800 at higher charging current density with respect to DP1000 for C L , indicating that DP800 has the highest value at higher charging current densities. The subsurface H concentration also has an inverse relationship with the D L . The steel with lower D L has higher sub-surface concentrations (Figures 7 and 8). It is understandable, since the H that enters the steel due to cathodic charging diffuses faster into the steel as it has a high D L value [73]. Further, it is observed that C Trapped saturates beyond 20 mA·cm −2 because of the finite trap sites. Such a distinct saturation is not noted for C T (Figure 9a). It appears that the increase in the H fugacity with applied current density enables higher solubility of hydrogen in the steel, although at high hydrogen fugacity, hydrogen recombination may be favored. Further, with an increase in the charging current density, the available free sites for H reduction are expected to decrease, due to the increased surface coverage by the adsorbed H atoms, limiting the H adsorption. As shown in Figure 9, when the charging current density is increased from −20 mA·cm −2 to −30 mA·cm −2 for DP800, the increase in total hydrogen concentration is 0.55 mol·m −3 , whereas at higher charging current densities, when it is raised from −40 mA·cm −2 to −50 mA·cm −2 , the change in hydrogen concentration is only 0.19 mol·m

Hydrogen Trapping
Reversible H traps also affect steels' HE tendency, as these traps act as reservoirs for hydrogen atoms and release them into the lattice over time. The experimental permeation falls transient curve is sensitive to the hydrogen in the reversible traps, with the area under the curve representing H desorbed from the lattice and reversible trap sites (C T* ). The detrapping kinetics of hydrogen from the reversible traps is calculated by subtracting a decay curve corresponding to D L (Equation (2)) from the experimentally obtained current decay curve [10,27,57,59]. Figure 10 shows the fall transients' schematic for DP1000 steel at different charging current densities and the corresponding calculated current due to oxidation of lattice H, using D L values. Assuming only one kind of trap site present in the material, and each trap site can hold only one H atom, the reversible trap density is calculated using Equation (6) [57].
where N T is the number of hydrogen trapping sites per unit volume (sites·cm −3 ), A is the difference in the area under the two curves (s), L is the thickness of the specimen (cm), i n is the steady-state current, and 6.24 × 10 12 is the number of hydrogen atoms in 1 µA·s·cm −2 .
The calculated N T summarized in Table 3 indicates that with increasing charging current density, the reversible trap sites of all the three steels decrease. Further, the calculated reversible trap density varies in the order DP800 > DP1000 > CP1000. The trap density emptied during fall transients increases with decreasing charging current densities. For example, in the DP800, lowering charging current density from −50 mA·cm −2 to −40 mA·cm −2 emptied only 12.5% of the total calculated traps, whereas lowering the charging current density from −20 mA·cm −2 to −10 mA·cm −2 , 46.2% of the total calculated traps were emptied. All the steels investigated demonstrated similar behavior. In DP steels, the calculated trap density was higher in DP800 than DP1000. It is possibly related to the energy associated with the traps. If the trap strength of DP1000 is stronger than DP800, it is understandable that fewer traps would be emptied with a reduction in charging current density (Table 3).

Trap Activation Energy
The binding energy of hydrogen with traps play a decisive role in determining the diffusion/trapping of H in the material. Following the works of Oriani [46] and later mod-

Trap Activation Energy
The binding energy of hydrogen with traps play a decisive role in determining the diffusion/trapping of H in the material. Following the works of Oriani [46] and later modified by Dong [74], the activation energy of the traps can be calculated from the obtained permeation data [75] using Equation (2): where N T is the number of hydrogen trap density (sites·cm −3 ), N L is the density of interstitial sites in the steel, E a is the activation energy for the traps, R is the gas constant and T is the absolute temperature.
Based on similar works available in the literature, the interstitial site density (N L ) is assumed 7.52 × 10 22 sites·cm −3 [20,76]. The average E a calculated for DP800 is 34.05 kJ·mol −1 , 37.61 kJ·mol −1 for DP1000 and 38.94 kJ·mol −1 for CP1000. The activation energy measured from the permeation technique is usually higher than the values obtained by TDS measurements [77]. In TDS, when the thermal energy is applied, all the hydrogen present in the material effuses out. However, in the permeation method, the decrease in the cathodic overpotential may not be sufficient to release all the hydrogen present in the traps. As trap density is inversely proportional to the activation energy, the E a calculated from permeation tests present a higher value. Though the activation energy measured from the permeation technique may not be accurate; nevertheless, it can provide a reasonable idea about the nature of traps.

Discussion
The present study suggests that the lattice H diffusion coefficient in CP steel was higher compared to the DP steels. For CP steel, the bainite promotes hydrogen diffusion [44,78] and in the case of percolate phase, the hydrogen effective diffusivity could be enhanced, as discussed by Fielding et al. [79]. Therefore, H diffusion path (percolate phases, metallugical defects) influences the observed diffusion coefficient values of these steels. As shown in Figure 5b,d,f, DP steels have a banded structure across its cross-section, which is not present in CP steel. The banded structure indicates a longer mean diffusible path for hydrogen in DP steels than in CP steel. In DP steels, the grain boundaries are the preferred diffusion pathway as the matrix is predominantly ferrite [80]. Koyama et al. [81], using the Ag decoration technique, reported that grain boundaries were the preferred diffusion pathways for hydrogen in ferrite. Haq et al. [56] reported that effective hydrogen diffusivity in specimens with banded structure perpendicular to the direction of hydrogen diffusion tends to be an order of magnitude lower than specimens with bands parallel to the direction of hydrogen diffusion. Olden et al. [82] reported that the greater the diffusion path and the higher the tortuosity of hydrogen, the lower the diffusion coefficient of hydrogen is. Zakromyzski discussed that the flux of hydrogen is directed along the layers of fibres of inclusions [27].
Among the three steels investigated, the C T , C L , and the C Trapped ( Figure 9) were significantly lower for CP1000 compared to both DP1000 and DP800. Looking at the carbon content in all the steels, as presented in Table 1, the C content in CP steel is relatively less than DP steels. Carbon concentration in steels is known to influence its HE susceptibility and, consequently, the hydrogen concentration in steel and its diffusion coefficient. Noh et al. [83] observed that an increase of C in austenitic stainless steels increases the susceptibility of hydrogen embrittlement [40]. For Hadam et Zakroczymski [64], the lattice solubility of hydrogen in Armco iron is 70 times lower than in high carbon steel, which is in line with the study of Noh et al., when the influence of other elements variations is neglected in the first approach.
In the absence of the knowledge of the charging schedule, an inference could be drawn that the hard phases, particularly martensite in CP1000 steel, is softer as compared to the DP1000. This is due to the higher phase fraction of hard phases in CP steel with lower C concentration as compared to DP1000. Softer martensite contains lesser dislocation density and less stressed ferrite regions around the ferrite/martensite interface [37]. As dislocation density presented in the martensite is considered a reversible trap site, CP1000 has a lesser number of reversible trap sites compared to DP steels, even though CP1000 has a higher grain boundary area due to its refined microstructure. It is consistent with the trap density observed in Table 3, explaining the small quantity of trapped H in CP steels.
However, though the trap density is smaller in CP1000, under similar charging conditions, a higher fraction of H is trapped in CP steels as compared to the DP steels. From Figure 9, if we calculate the fraction of the total absorbed hydrogen present in traps, it is evident that a higher fraction of H is trapped in CP1000, as shown in Figure 11. In CP steels, there is a significant portion of bainite, which is composed of cementite precipitates. These cementite particles are known to be strong hydrogen traps and are the major trapping sites in the bainitic phase [37]. The size and morphology of the carbide particles also play an essential role in determining the permeability of hydrogen. Ramunni et al. [84], observed finer carbides as being more efficient in hydrogen trap sites. The higher trap activation energy of CP steel also supports the present observation (Table 3). Chun et al. [85] reported that the Ea associated with dislocations and grain boundaries to be around 35 kJ·mol −1 . Luppo et al. [86], in their study of low carbon martensitic steel, found that Ea for lath interfaces to be 46 kJ·mol −1 . Wong et al. [87] found the energy required to release hydrogen from dislocations to be around 21 kJ·mol −1 . Based on these previous results, the hydrogen released from the reversible traps may be present in the grain boundaries, dislocations, and martensite lath interfaces in DP and CP steels along with cementite particles present in CP steels. The low hydrogen uptake, the lower trap density, and the presence of stronger traps in CP steels are expected to render this steel more resistant to HE as compared to the DP steels [16,37,63,76,[88][89][90]. Figure 9 suggests that CL is higher for DP1000 at low charging current density, but at higher charging current density, CL is higher for DP800 (Figure 9b). Both DP steels have ferrite and martensite, with higher martensite fraction in DP1000. Martensite, due to its high dislocation density, acts as a hydrogen trap, and resultantly DP1000 has higher trap density than DP800. Thus, the traps in DP800 are quickly get saturated compared to the traps in DP1000, leading to enrichment of hydrogen in DP800 lattice. This observation is particularly interesting in terms of embrittlement phenomenon as Koyama et al. [91] observed that the crack initiation in DP steels occurred through martensite cracking and propagated through the ferrite/martensite interface. The martensite or the martensite/ferrite interface thus can act as strong or weak H trapping sites depending on the activation energy. Chun et al. [85] reported that the E a associated with dislocations and grain boundaries to be around 35 kJ·mol −1 . Luppo et al. [86], in their study of low carbon martensitic steel, found that E a for lath interfaces to be 46 kJ·mol −1 . Wong et al. [87] found the energy required to release hydrogen from dislocations to be around 21 kJ·mol −1 . Based on these previous results, the hydrogen released from the reversible traps may be present in the grain boundaries, dislocations, and martensite lath interfaces in DP and CP steels along with cementite particles present in CP steels. The low hydrogen uptake, the lower trap density, and the presence of stronger traps in CP steels are expected to render this steel more resistant to HE as compared to the DP steels [16,37,63,76,[88][89][90]. Figure 9 suggests that C L is higher for DP1000 at low charging current density, but at higher charging current density, C L is higher for DP800 (Figure 9b). Both DP steels have ferrite and martensite, with higher martensite fraction in DP1000. Martensite, due to its high dislocation density, acts as a hydrogen trap, and resultantly DP1000 has higher trap density than DP800. Thus, the traps in DP800 are quickly get saturated compared to the traps in DP1000, leading to enrichment of hydrogen in DP800 lattice. This observation is particularly interesting in terms of embrittlement phenomenon as Koyama et al. [91] observed that the crack initiation in DP steels occurred through martensite cracking and propagated through the ferrite/martensite interface. The martensite or the martensite/ferrite interface thus can act as strong or weak H trapping sites depending on the activation energy.
Another interesting observation is made by comparing H concentration change at the charging side and discharging side during each fall transient. The change in the hydrogen concentration at the discharge side (C T* ) is calculated by measuring the area of the experimental curve. It is observed that for all the steels, the concentration of hydrogen coming out of the discharge side is almost a third of the change in the concentration of the charging side. The ratio between the change in the hydrogen concentration at the discharging side and at the charging side is 0.36 for DP800, 0.33 for DP1000, and 0.30 for CP1000 steels (Table 3). This helps to get an insight into the possible hydrogen distribution within the metal. Hadam et al. [64] suggested that H distribution within the material is uneven along the material cross-section, with most of the H trapped at the sub-surface on the charging side of the specimen. Zakroczymski [64] showed that most of the hydrogen desorbed through the charging side of the specimen. After the permeation current reached the steady state, Zakroczymski interrupted the cathodic charging, and an anodic potential equal to the discharging side was applied to the charging side, which allowed the measurement of desorbed hydrogen from both surfaces. The hydrogen detected from the discharging side accounted for only 12.4% of the total hydrogen released. In the present test conditions, there is the possible release of hydrogen from both sides (charging and discharging). As a result, the total hydrogen desorbed from the discharge side is around one-third of the total change in the hydrogen concentration at the charging side. It must be remembered that in this case, there is no anodic potential applied at the charging side after steady-state, but the only driving force of desorption of hydrogen from the charging side is the lowering of cathodic current density. This interpretation is, however, just speculative and needs further validation as the H concentration profile is extremely difficult to be verified experimentally [92]. However, if true, such observations can provide better insight to interpret the trapping effect of different microstructural structures in the material.

Conclusions
Permeation characteristics of DP and CP steels were investigated in the study using the electrochemical hydrogen permeation technique to understand the microstructural influence. The influence of charging surface modification during permeation experiments was successfully mitigated by pre-charging the specimens before the permeation experiment. The salient findings of the work are summarized below: • D L in CP1000 steel is twice as compared to that of the DP800 and DP1000 steels. • Higher D L in CP1000 steel is due to the microstructural feature that lowered trap density and enabled easier H diffusion and decreased H concentration at the subsurface. The relatively lower carbon concentration in CP steels is a possible reason for the low H sub-surface concentration.

•
Higher fraction of the trapped hydrogen in CP steel is due to the presence of cementite particles. • CP steel had the least trap density compared to DP steels but had a higher trap strength. • C L is higher for DP800 due to easier saturation of traps in DP800 leading of H enrichment in the lattice. • H concentration effusing from the discharge side is almost a third from that of the charging side, suggesting a non-uniform distribution of H across the cross-section, and higher H concentration near the charging surface. Funding: This research received no external funding.

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Data Availability Statement:
No new data were created or analyzed in this study. Data sharing is not applicable to this article.