Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
Scientific Reports volume 12, Article number: 18303 (2022 ) Cite this article China Graphite
The study employed Fick’s second law of diffusion to discover some unknown aspect of carbon diffusion in grey cast iron during carburisation process. Emphasis on the experiments and theoretical modelling were established for better accomplishments. Pulverised palm kernel and eggshell additives of 70 (wt.%) and 30 (wt.%) according to the Voige law of mixture was considered as a continuous medium without considering the atomic nature of the mixture. Furthermore, a kinetic approach was described where a physical model of the substrate immersed in the carbon mixture was established while diffusion equations were modelled to establish the mechanism of carbon diffusion during carburisation. Initial composition and concentration of diffused atom remained constant which are 2.68 and 6.67% carbon. While the carburizing time used varied from 60 min, 90 min, 120 min, 150 min, 180 min and 210 min respectively at constant carburising temperature of 900° The results revealed varying composition gradient of carbon atom ranging from 5.4%, 5.42%, 5.44%, 5.46%, 5.51%, and 5.65 compared to the initial carbon content of 2.68%. The concentration of carbon atom on the substrate surface at varying time implies that the process was non-steady state diffusion which verified Fick’s second law of diffusion. Hence, the composition achieved is a function of boundary conditions such as time position and temperature. This novel study will enhance the understanding of heat treat treatment of metals such that their applications in the industry will be numerous.
Metallic materials which have undergone heat treatment via carburisation process have surfaces characterised with improved mechanical properties1. They are basically modified for advanced engineering applications using diffusion mechanism2. Diffusion involved the squeezing of carbon atoms past its surrounding atoms in other to reach a new position. Diffusion process can be best understood from the equation parameter of Fick’s law as well as the knowledge of the activation energy required for the diffusion process3. For instance, Fick’s second law established a non-steady state diffusion of atoms as described by the differential equation \(\frac{dc}{{dt}} = \frac{{Dd^{2} }}{{dx^{2} }}\) of which the solution is a function of a particular diffusion process described by the boundary parameters in the Eq. 14.
The solution to Fick’s second law permits the evaluation of concentration of an atom diffused near the surface of the coupon material as a function of time and distance, provided that the coefficient of diffusion D remains constant and the concentration of the atom at the surface \(C_{s}\) as well as within the material \(C_{0}\) remain unchanged5. Recent study on the diffusion of palm kernel and eggshell additives to grey cast iron resulted in increase in the hardness of the material6. The tribological properties of the treated material via diffusion process was excellent which made it suitable for advanced engineering material7,8. The principle of Fick’s second law had been limitedly used in evaluating the depth of the mechanical properties that had been diffused into these materials, thus making the analysis or establishing the statistical significance of the diffused carbon atom a major problem9,10,11. One major problem in diffusion analysis is the determination of the temperature field and depth of carbon imposed on the surroundings of the substrate metal12,13. Study has shown that the knowledge of the temperature distribution could be a pointer to understanding the carbon diffusion mechanism as well as the depth of diffusion14. For grey cast iron material, knowledge of the carbon diffusion is important in analysing the structural integrity. Also, knowledge of carbon diffusion is critical to the optimisation of the coating thickness as well as the compatibility of the carburising agents15,16,17,18,19,20. However, movement of atoms is an essential factor for the diffusion process to take place in metals. Therefore, understanding the dynamics of the diffusion process remain a crucial problem in the determination of the carbon depth in solid materials21.
Furthermore, carbon content such as that from palm kernel, coconut, wood charcoal and eggshells have been the major carbon content that are usually deployed during carburisation process. The eco-friendliness and their environmental friendliness make them most promising materials for advanced heat treatment applications. However, the commercial applications of these organic materials are limited by some drawbacks including the determination of the depth of carbon diffusion22. To address these problems and improve on the applications of these organic materials, various material, methods and concepts have been designed and developed. For instance, Chen et al.23 developed a light-weight thermal energy concrete and reinforced it with carbon from palm kernel shell. This helped in improving the thermal lag and lowered the peak temperature as a composite. More so, alloying of steel with borides under high temperature resulted in increase in the hardness of the substrate24 and the thickness of the layers increased with increased temperature as reported by Hu et al.25. It has also been reported that carbon diffuses into the surface and settles at the face centred cubic region during the carburisation of austenitic steel. This therefore, caused an increase in the hardness of the layers formed and leaving chromium in its free form and allows an increase in corrosion resistance as well as improvements in tribological and mechanical properties. However, it is important to understand that in an ideal diffusion process, a diffusing atom squeezes past the surrounding atoms to get to the new position. This implies that energy will be utilized to force the atom to its new position. Thus, the energy barrier needed to move the atom to its new position is referred to as the activation energy26,27,28,29,30.
In our present research, diffusion in grey cast iron using pulverised palm kernel and eggshell additives showed an increase in the hardness of the substrate material. Therefore, the aim of this study is to deploy the Fick’s Second Law in the determination of the depth of carbon in the material that resulted in increased hardness after the carburisation process. This method had become prominent in the understanding of diffusion in solids, liquids and gases.
The experimental approach involved the use of pulverised palm kernel and eggshell additives of 70 (wt.%) of pulverised palm kernel and 30 (wt.%) of pulverised eggshells according to the Voige law of mixture. Grey cast iron substrates of dimensions (20 mm × 20 mm × 10 mm) and chemical composition (wt.%) of 2.68 C, 1.42 Si, 0.63 Mn, 0.13 S, 0.28 P was prepared using different grades of silicon carbide abrasives to obtain a polished and smooth surface for easy carbon diffusion. Prepared grey cast iron substrate were embedded into some stainless steel containers to reduce the rate of carbon absorption and were finally charged into a muffle furnace of 1200 ℃ capacity. The carburisation process was carried out at a temperature of 900 ℃ for 60, 90,120, 150, and 180 min after which it was stopped and cooling was done using water. Water was used as a quenching medium due to its natural convective heat transfer ability and the tendency of increasing the hardness property of the substrate. Initial composition and concentration of diffused atom remained constant which are 2.68 and 6.67%. This is to be able to determine the percentage of carbon atom diffused into the material surface, depth of diffusion at a particular time.
The main aim of this analyses is to determine the temperature field impacted by the surrounding medium (carburisers) on the carburised material (grey cast iron). Since, conduction process takes place for proper diffusion of carbon into the material, it is important to understand and analyse the temperature distribution around the material. To achieve this, Fourier’s law of heat conduction was employed.
Consider a grey cast iron substrates of dimensions (20 mm × 20 mm × 10 mm) immersed in a homogenous medium in which there is no bulk motion as demonstrated in Fig. 1. Let the temperature distribution be described by a Cartesian coordinates in the form, \(T_{G } \left( {x_{G} , y_{G} , z_{G} } \right)\) . Assuming an infinitesimally control volume defined by \(dx_{G} \times dy_{G } \times dz_{G}\) . Conduction heat transfer occurred at each control surface due to temperature variation, therefore the rate of heat conducted perpendicular to the coordinates \(x_{G} - , y_{G} - , z_{G}\) can be denoted by \( Q_{{x_{G} }} ,\;Q_{{y_{G} }} ,\;Q_{{z_{G} }}\) respectively.
Conduction analysis in Cartesian coordinate.
Therefore, the rate of heat conduction is given by the equations;
The thermal energy generated within the medium is given by;
Based on the thermal energy generated, there was variation in the internal thermal energy stored by the grey cast iron material which is being carburised. However, in the absence of phase transformation, latent energy is minimal and the stored energy is given by;
The conservation energy required within the medium is given by;
Since the rate of heat conduction involve energy input and output and substituting Eqs. 5 and 6 into Eq. 7 yield
Substituting from Eqs. 2–4, into Eq. 9 yields;
Assuming that the material become isotropic after the carburisation process, then the conduction heat rate can be established from Fourier’s law;
Each equation from Eq. 11 represents the heat flux resulting from diffusion across the substrate
Substituting the heat flux Eqs. 11, 12, and 13 into Eq. 9 and dividing through the control volume \(\left( {dx_{G} dy_{G } dz_{G} } \right)\) , we obtain
Equation 15 is the heat transfer equation when the substrate is represented in the Cartesian coordinate. This model establishes the fundamentals of the heat conducted through the substrate material. Thus, it is possible to obtain the temperature distributed during the carburisation process.
Assuming the temperature distribution across the substrate is ven by the equation;
This equation is expressed as a function of time and also establishes the energy conservation. Thus, \(\frac{\partial }{{\partial x_{G} }}\left( {k\frac{\partial T}{{\partial x_{G} }}} \right)\) is related to the net heat conducted into the substrate in x-coordinate direction
The same can be expressed in the y- and z-coordinates to obtain Eqs. 18 and 19
Therefore, for a constant thermal conductivity, Eq. 15 can be rewritten as;
where \(\alpha = \frac{K}{{\rho C_{p} }}\) is the thermal diffusivity.
According to Fick’s Second law, carbon atom diffused during the carburisation can be defined by the differential equation in the form [32];
And the boundary conditions for the carburisation process depends on the equation [32];
where \(c_{0}\) the carbon content for the as-received grey cast iron which is given by 2.68%, \(c_{s}\) assumed value from iron-carbon alloy system diagram between pure iron an interstitial compound, iron carbide (Fe3C), containing 6.67% carbon. \(c_{x}\) Implies the concentration of the diffused carbon at a depth denoted by x in millimetre below the material surface at time t. From the carburisation result, the value for \(c_{x} = 5.40\%\) at t = 60 min which is 3600 s, D is the diffusion coefficient and it remain constant for constant, \(c_{s}\) provided \(c_{0} \) also remain constant. D = 2 × 10−11 (m2 s−1) of carbon in FCC iron interstitial carbon diffusion [33].
Thus, based on these conditions, the solution to Fick’s second law enabled the study to determine the concentration of a diffused carbon atom as a function of carburising time and distance (depth). erf = error function which is given by 0.71 [34], t = time of carburisation which is given as 3600 s, x = is the depth (mm)?
Substituting the values into Eq. 22
Figures 2, 3, 4, 5, 6 and 7 illustrate the variation in percentage composition of carbon atoms diffused at various time at a temperature of 900°. From Fig. 2, it was observed that 5.4% of carbon diffused into the material at a depth of 0.0144 mm at 60 min. Comparing this with the initial composition of 2.68%, there was an increase of about 2.72% in carbon content on the material surface. Thus, the driving force for carbon diffusion led to microstructural changes in the substrate carbon content.
Diffusion of carbon atoms into the surface of grey cast iron at 900° for 60 min.
Diffusion of carbon atoms into the surface of grey cast iron at 900° for 90 min.
Diffusion of carbon atoms into the surface of grey cast iron at 900° for 120 min.
Diffusion of carbon atoms into the surface of grey cast iron at 900° for 150 min.
Diffusion of carbon atoms into the surface of grey cast iron at 900° for 180 min.
Diffusion of carbon atoms into the surface of grey cast iron at 900° for 210 min.
In addition, Fig. 3 illustrates the carbon atom diffusion at 90 min. About 5.40% was deposited at a depth of 0.000294655 mm. about 0.02% increase compared to the diffusion at 60 min. Further to this, with the temperature of diffusion still kept at 900° while increasing the diffusion time (holding time) to 120 and 150 min resulted in increased composition in carbon atom to about5.44 and 5.46 respectively as described by Figs. 4 and 5. The corresponding depth of deposition into the material surface were 9.13E-54 and 0.080498 as well. Thus, the result established the fact that, increased diffusion time will lead to increase in carbon deposition as well.
Similarly, Figs. 6 and 7 depicts the composition of carbon atom at 180 and 210 min respectively at the same temperature of carburisation. The composition at this time were 5.51 and 5.65 respectively at a distance of 0.043305317 and 0.0577404231. In comparison with the initial carbon composition, it was observed that an increase of about 2.83 and 2.97 was added at various time of carburisation respectively. Thus, the solution to Fick’s second law permits the study to establish the composition of diffusing carbon atom near the substrate surface as a function of time and distance respectively.
Hence, carburisation was defined in this study as the diffusion of carbon into the substrate metal. The purpose was to determine the concentration gradient of carbon at varying time since the entire process involved thermal effects. Thus, the processing parameters such as temperature and time are the key influencers of diffused carbon potentials. From the results, it is noteworthy to say that irrespective of the interaction between the carbon and other alloying elements present in the carbon, diffusion of carbon into the substrate (grey cast iron) was a perfect description of Fick’s Second Law of diffusion. This is because the composition gradient of carbon atom near the surface of the substrate vary with time due to carbon accumulation. This is the reason why it is called a non-steady state diffusion as expressed in the modelled equations.
Figure 8 showed the SEM microstructure of the as-received grey cast iron before carburisation. Similarly, Fig. 9 showed the microstructure of grey cast iron carburised at 900 °C at 60 min. From the Figure, it was observed that graphite dominates the metal surface which is traceable to the nature of carbon content used in the carburisation experiment. More so, Fig. 10 also depict the microstructure of the carburised sample at 900 °C at ninety (90) minutes. The surface was observed to have an increased precipitate of carbon, the same characteristics was observed in Figs. 11–12 for samples carburised at 120 and 150 min respectively. However, at 180 min (Fig. 13), the metal surface was characterised with graphite precipitates. Thus, the presence of graphite indicates that diffusion actually took place during the carburisation process.
SEM microstructure of As-received grey cast iron.
SEM microstructure of grey cast iron carburized at 900 °C for 60 min.
SEM microstructure of grey cast iron carburized at 900 °C for 90 min.
SEM microstructure of grey cast iron carburized at 900 °C for 120 min.
SEM microstructure of grey cast iron carburized at 900 °C for 150 min.
SEM microstructure of grey cast iron carburized at 900 °C for 180 min.
The initial micro-hardness of the as-received cast iron was 288.41HV, which increased to 355.8 HV after carburisation at 180 min. This showed that there was an increase of about 67.39HV compared to the initial hardness value. This increment is a function of the adjustment in the characteristic parameters of the carburised material. Furthermore, the fatigue and wear resistance of grey cast iron component has been improved.
From the study, it was established that atoms of the carbon used during carburisation moved in different fashion to eliminate concentration differences and eventually produce a homogeneous deposit on the grey cast iron material. These deposits are of varying composition. Thus, the understanding of atomic movements during the diffusion process has been established. More so, it was established that a diffusing atom of the carbon squeezed past the other surrounding atoms until it is deposited on the material. This process requires energy supply to force the atom to be deposited. This the reason for the temperature of carburisation.
In addition, the presence of graphite in the microstructures as observed in the carburised samples would help in the improvement of wear resistance because graphite presence serves as lubricant to the material. Thus, thermal fatigue is reduced. More so, increase in microhardness will help in reduction of surface abrasion of the grey cast iron during application.
Furthermore, Fick’s second law had been used to establish the non-steady state diffusion of carbon atoms into the substrate material (grey cast iron). However, a consequence of Fick’s second law is the possibility of achieving the constant composition profile for varying conditions, provided Dt remain constant. For several heat treatment and for most industrial applications, it will be possible to determine the temperature effect at varying time.
All data generated or analysed during this study are included in this published article.
No code was used for the computation of the data reported in this study.
Temperature distribution in the grey cast iron sample
X-coordinate in the grey cast iron sample
Y-coordinate in the grey cast iron sample
Z-coordinate in the grey cast iron sample
Infinitesimal control volume in the x, y and z-coordinate in the grey cast iron sample
Rate of heat conduction perpendicular to the coordinates
Thermal energy generated in the grey cast iron sample
Thermal energy stored by grey cast iron
Conduction heat rate in various directions
The carbon content for the as-received grey cast iron
Assumed value from iron-carbon alloy system diagram between pure iron an interstitial compound, iron carbide (Fe3C)
The concentration of the diffused carbon at a depth denoted by xn millimetre below the material surface at time t
Error function from Fick’s second law
Time of carburisation in seconds
Karunaratne, M. S. A., Yan, S., Thomson, R. C., Coghlan, L. & Higginson, R. L. Modelling carburisation in 9Cr-1Mo ferritic steel tube substrates in experimental CO2 atmospheres. Corros. Sci. 163, 108248 (2020).
Biglari, F. R. & Nikbin, K. M. A diffusion driven carburisation combined with a multiaxial continuum creep model to predict random multiple cracking in engineering alloys. Eng. Fract. Mech. 146, 89–108 (2015).
Askeland, D. R., Phulé, P. P., Wright, W. J. & Bhattacharya, D. K. The science and engineering of materials. Science 2, 558 (2003).
Dos Santos, T. J., Tavares, F. W. & Abreu, C. R. Fick diffusion coefficients via molecular dynamics: An alternative approach in the Fourier domain. J. Mol. Liq. 329, 115460 (2021).
Zhao, X., Luo, T. & Jin, H. A predictive model for self-, Maxwell-Stefan, and Fick diffusion coefficients of binary supercritical water mixtures. J. Mol. Liq. 324, 114735 (2021).
Salawu, E. Y., Ajayi, O. O., Inegbenebor, A., Akinlabi, S. & Akinlabi, E. Influence of pulverized palm kernel and eggshell additives on the hardness, coefficient of friction and microstructure of grey cast iron material for advance applications. Results Eng. 3, 100025 (2019).
Salawu, E. Y. et al. Investigation of the effects of selected bio-based carburising agents on mechanical and microstructural characteristics of gray cast iron. Heliyon 6(2), e03418 (2020).
Article PubMed PubMed Central Google Scholar
Bailey, R. & Sun, Y. Pack carburisation of commercially pure titanium with limited oxygen diffusion for improved tribological properties. Surf. Coat. Technol. 261, 28–34 (2015).
Nicolin, D. J., Rossoni, D. F. & Jorge, L. M. M. Study of uncertainty in the fitting of diffusivity of Fick’s Second Law of Diffusion with the use of Bootstrap Method. J. Food Eng. 184, 63–68 (2016).
Rauma, M. & Johanson, G. Assessment of dermal absorption by thermogravimetric analysis: Development of a diffusion model based on Fick’s second law. J. Pharm. Sci. 98(11), 4365–4375 (2009).
Article CAS PubMed Google Scholar
Diao, J. & Hess, D. W. Use of angle-resolved XPS to determine depth profiles based on Fick’s second law of diffusion: Description of method and simulation study. J. Electron Spectrosc. Relat. Phenom. 135(2–3), 87–104 (2004).
Paul, A., Laurila, T., Vuorinen, V., & Divinski, S. V. (2014). Fick’s laws of diffusion. In Thermodynamics, diffusion and the kirkendall effect in solids (pp. 115–139). Springer, Cham.
Cantor, B. Fick’s laws: Diffusion. In The equations of materials (pp. 141–161). Oxford University Press.
Rashidi, N. A. & Yusup, S. Co-valorization of delayed petroleum coke–palm kernel shell for activated carbon production. J. Hazard. Mater. 403, 123876 (2021).
Article CAS PubMed Google Scholar
Young, D. et al. Non-steady state carburisation of martensitic 9–12% Cr steels in CO2 rich gases at 550 ℃. Corros. Sci. 88, 161–169 (2014).
Liu, C., Heard, P. J., Griffiths, I., Cherns, D. & Flewitt, P. E. J. Carbide precipitation associated with carburisation of 9Cr–1Mo steel in hot CO2 gas. Materialia 7, 100415 (2019).
Zhang, J., Li, H., Kong, C. & Young, D. J. Oxidation and carburisation of Fe–6Al/Fe–6Al–3Si in dry and wet CO2 gases. Corros. Sci. 74, 256–264 (2013).
Article ADS CAS Google Scholar
He, Y. et al. Formation of hollow nanofiber rolls through controllable carbon diffusion for Li metal host. Carbon 157, 622–630 (2020).
Tseng, S. C., Lee, T. C. & Tsai, H. Y. Diffusion effect for the catalytic growth of carbon nanotubes on metal alloys substrate. Diam. Relat. Mater. 96, 112–117 (2019).
Article ADS CAS Google Scholar
Zeng, Y., Li, Q. & Bai, K. Prediction of interstitial diffusion activation energies of nitrogen, oxygen, boron and carbon in bcc, fcc, and hcp metals using machine learning. Comput. Mater. Sci. 144, 232–247 (2018).
Pongsajanukul, P. et al. Theoretical study of carbon dioxide adsorption and diffusion in MIL-127 (Fe) metal organic framework. Chem. Phys. 491, 118–125 (2017).
Lobo, L. S., Figueiredo, J. L. & Bernardo, C. A. Carbon formation and gasification on metals. Bulk diffusion mechanism: A reassessment. Catal. Today 178(1), 110–116 (2011).
Chin, C. O., Yang, X., Paul, S. C., Wong, L. S. & Kong, S. Y. Development of thermal energy storage lightweight concrete using paraffin-oil palm kernel shell-activated carbon composite. J. Clean. Prod. 261, 121227 (2020).
Buijnsters, J. G. et al. Diffusion-modified boride interlayers for chemical vapour deposition of low-residual-stress diamond films on steel substrates. Thin Solid Films 426(1–2), 85–93 (2003).
Article ADS CAS Google Scholar
Hu, X. J., Zhang, B. M., Chen, S. H., Fang, F. & Jiang, J. Q. Oxide scale growth on high carbon steel at high temperatures. J. Iron. Steel Res. Int. 20(1), 47–52 (2013).
Sun, Y. & Haruman, E. Tribocorrosion behaviour of low temperature plasma carburised 316L stainless steel in 0.5 M NaCl solution. Corros. Sci. 53(12), 4131–4140 (2011).
Sarrade, S. et al. Overview on corrosion in supercritical fluids. J. Supercrit. Fluids 120, 335–344 (2017).
Yang, P., Liu, C., Guo, Q. & Liu, Y. Variation of activation energy determined by a modified Arrhenius approach: Roles of dynamic recrystallization on the hot deformation of Ni-based superalloy. J. Mater. Sci. Technol. 72, 162–171 (2021).
Kalaivanan, R., Ganesh, N. V. & Al-Mdallal, Q. M. An investigation on Arrhenius activation energy of second grade nanofluid flow with active and passive control of nanomaterials. Case Stud. Therm. Eng. 22, 100774 (2020).
Baruah, B., Anand, R. & Behera, S. K. Master sintering curve and activation energy of sintering of ZrO2-doped Al2O3. Ceram. Int. 47(5), 7253–7257 (2021).
The authors did not receive support from any organization for the submitted work.
Department of Mechanical Engineering, Covenant University, P.M.B 1023, Ota, Ogun State, Nigeria
Enesi Y. Salawu, Oluseyi O. Ajayi, Anthony O. Inegbenebor & Joseph O. Dirisu
Department of Mechanical Engineering, Landmark University, P.M.B 1001, Omu-Aran, Kwara State, Nigeria
Department of Mechanical Engineering Science, University of Johannesburg, Johannesburg, South Africa
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
E.Y.S., O.O.A., A.O.I., A.A.A., J.O.D., and R.O.L. had the idea for the article, E.Y.S., O.O.A., A.O.I., A.A.A., J.O.D., and R.O.L. performed the literature search and data analysis. E.Y.S., O.O.A., A.O.I., A.A.A., J.O.D., and R.O.L. drafted and/or critically revised the work. E.Y.S., O.O.A., A.O.I., A.A.A., J.O.D., and R.O.L. critically read and approved the final manuscript.
Correspondence to Enesi Y. Salawu or Adeolu A. Adediran.
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/.
Salawu, E.Y., Adediran, A.A., Ajayi, O.O. et al. On the analyses of carbon atom diffused into grey cast iron during carburisation process. Sci Rep 12, 18303 (2022). https://doi.org/10.1038/s41598-022-22136-w
DOI: https://doi.org/10.1038/s41598-022-22136-w
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
Scientific Reports (Sci Rep) ISSN 2045-2322 (online)
Calcined Petroleum Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.