# Fileset

[Improving efficiency of autonomous material search via transfer learning from nontarget properties (2).pdf](https://mdr.nims.go.jp/filesets/115eb31f-0de2-4346-ba62-6df7b51f786f/download)

## Creator

Jaekyun Hwang, [Yuma Iwasaki](https://orcid.org/0000-0002-7117-277X)

## Rights

[Creative Commons BY Attribution 4.0 International](https://creativecommons.org/licenses/by/4.0/)

## Other metadata

[Improving efficiency of autonomous material search via transfer learning from nontarget properties](https://mdr.nims.go.jp/datasets/1bdec3f6-50ca-4eee-acfe-12d1ca6e53b9)

## Fulltext

Improving efficiency of autonomous material search via transfer learning from nontarget propertiesFull Terms & Conditions of access and use can be found athttps://www.tandfonline.com/action/journalInformation?journalCode=tstm20Science and Technology of Advanced Materials: MethodsISSN: (Print) (Online) Journal homepage: https://www.tandfonline.com/loi/tstm20Improving efficiency of autonomous materialsearch via transfer learning from nontargetpropertiesJaekyun Hwang & Yuma IwasakiTo cite this article: Jaekyun Hwang & Yuma Iwasaki (2023) Improving efficiency of autonomousmaterial search via transfer learning from nontarget properties, Science and Technology ofAdvanced Materials: Methods, 3:1, 2254202, DOI: 10.1080/27660400.2023.2254202To link to this article:  https://doi.org/10.1080/27660400.2023.2254202© 2023 The Author(s). Published by NationalInstitute for Materials Science in partnershipwith Taylor & Francis GroupPublished online: 28 Sep 2023.Submit your article to this journal Article views: 157View related articles View Crossmark datahttps://www.tandfonline.com/action/journalInformation?journalCode=tstm20https://www.tandfonline.com/loi/tstm20https://www.tandfonline.com/action/showCitFormats?doi=10.1080/27660400.2023.2254202https://doi.org/10.1080/27660400.2023.2254202https://www.tandfonline.com/action/authorSubmission?journalCode=tstm20&show=instructionshttps://www.tandfonline.com/action/authorSubmission?journalCode=tstm20&show=instructionshttps://www.tandfonline.com/doi/mlt/10.1080/27660400.2023.2254202https://www.tandfonline.com/doi/mlt/10.1080/27660400.2023.2254202http://crossmark.crossref.org/dialog/?doi=10.1080/27660400.2023.2254202&domain=pdf&date_stamp=28 Sep 2023http://crossmark.crossref.org/dialog/?doi=10.1080/27660400.2023.2254202&domain=pdf&date_stamp=28 Sep 2023Improving efficiency of autonomous material search via transfer learning from nontarget propertiesJaekyun Hwang and Yuma IwasakiCenter for Basic Research on Materials, National Institute for Materials Science, Tsukuba, JapanABSTRACTRecently, autonomous material search methods combining machine learning and experiments/simulations have become indispensable for exploring the extremely vast material exploration space. However, conventional autonomous material search methods focus solely on target material properties and their descriptors, leaving room for improvement in search efficiency. More efficient autonomous material search can be realized by utilizing information on nontarget properties that are dormant in databases. Here, we propose a novel method for autonomous material search using transfer learning and an ensemble neural network. This method can perform an autonomous material search to optimize the target properties while transferring information on nontarget properties. To demonstrate the usefulness of this method, we applied it to search for ternary magnetic alloys with a high Curie temperature while transferring information on magnetic moment and spin polarization. Results indicate that the proposed method improves the efficiency of autonomous material search, particularly in the early stages of the search process.ARTICLE HISTORY Received 24 April 2023  Revised 7 August 2023  Accepted 28 August 2023 KEYWORDS Machine learning; ensemble neural network; transfer learning; autonomous material search1. IntroductionRecently, the increasing complexity of materials has led to a massive expansion in the material search space. Traditional material search methods involve a laborious process of material synthesis, measurements, and analysis performed by humans. However, with the increasing size of the search space, conventional methods have proven to be inadequate. To address this challenge, autonomous material search methods based on machine learning have been developed. These methods leverage the power of machine learning to identify promising materials without requiring human intervention. One example of this is robotic autonomous material search, which combines material experiments, robotics, and machine learning to perform autonomous material synthesis and measurements [1–9]. Another example is in silico autonomous material search, which combines ab initio calculations and machine learning to efficiently search for materials [10–15]. Using machine learning to identify promising candidate materials and process conditions, these methods can considerably accelerate material discovery.Typically, conventional autonomous material search methods focus on optimizing a target property using CONTACT Yuma Iwasaki IWASAKI.Yuma@nims.go.jp Center for Basic Research on Materials, National Institute for Materials Science, 1-1 Namiki, Tsukuba-shi, Ibaraki-ken 305-0044, JapanSCIENCE AND TECHNOLOGY OF ADVANCED MATERIALS: METHODS 2023, VOL. 3, NO. 1, 2254202 https://doi.org/10.1080/27660400.2023.2254202© 2023 The Author(s). Published by National Institute for Materials Science in partnership with Taylor & Francis Group  This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The terms on which this article has been published allow the posting of the Accepted Manuscript in a repository by the author(s) or with their consent.http://orcid.org/0000-0002-1792-0044http://orcid.org/0000-0002-7117-277Xhttp://www.tandfonline.comhttps://crossmark.crossref.org/dialog/?doi=10.1080/27660400.2023.2254202&domain=pdf&date_stamp=2023-09-27Bayesian optimization with data solely based on the target property. Figure 1(a) illustrates this approach, where the autonomous material search loop uses only the data of the target property and its material descriptors (such as composition and crystal structure) to maximize the target property. In the machine learning part of Figure 1(a), a Gaussian regression model is commonly employed: YTP ¼ f XDð Þ;where YTP and XD represent the target property (objective variable) and material descriptors (explanatory variables), respectively. This model predicts the target property of new materials and quantifies their uncertainties, enabling the determination of material descriptor values for the next material synthesis or material simulation. This autonomous loop of material exploration facilitates the continuous exploration of the material space without requiring human intervention.However, this conventional approach overlooks the potential benefits of incorporating information on nontarget properties that exist in material databases, limiting the search efficiency. Figure 1(b) shows the process, whereby a machine learning model is constructed by leveraging both target and nontarget property data. Adding nontarget property data can improve the accuracy of the machine learning model and subsequently enable an efficient autonomous material search. However, in the context of autonomous material search, it is impractical to include nontarget property data as explanatory variables in the machine learning model as in YTP ¼ f XD;XNTPð Þ:Such a model would restrict the search space to those materials with information on the nontarget property XNTP. Thus, another autonomous material search approach that can effectively harness the potential of nontarget properties in existing material databases is necessary.Herein, we introduce a new approach for autonomous material search, illustrated in Figure 1(b). Our proposed method streamlines the process of optimizing target properties while simultaneously considering nontarget properties, affording an efficient search for materials with desired characteristics. To show the effectiveness of our approach, we apply it to data on the magnetic moment M, spin polarization S, and Curie temperature Tc of ternary magnetic B2 structure materials.2. Proposed methodTypically, Bayesian optimization with the Gaussian process regression model is the standard approach for autonomous material search. However, incorporating nontarget property information into the Gaussian process regression model can be challenging. Therefore, instead of the Gaussian process regression, we utilized an ensemble neural network (ENN) approach [16]. An ENN employs multiple independently trained neural networks (NNs), producing results superior to those obtained by employing a single NN. For instance, in numerical prediction tasks, taking the average of predicted values from NNs can enhance the prediction accuracy. Additionally, an ENN can approximate the distribution of predicted values, enabling us to quantify uncertainty akin to the Gaussian process. Consequently, the average of predicted values from an ENN serves as the prediction output for autonomous material search, with the variance of predicted values reflecting the degree of uncertainty.Moreover, the ENN approach offers transfer learning, a technique that enhances model performance by pretraining with diverse data beforehand and retraining the model with the target data [17]. In autonomous material search, transfer learning enables us to leverage nontarget property information. Herein, we transferred the knowledge obtained by pretraining an ENN model with nontarget property data to the ENN model utilized in the autonomous material search.We employed an ENN comprising 10 NNs, each comprising three hidden layers with 100, 100, and 10 neurons, respectively, from the input to the output Figure 1. Comparison of conventional and proposed new autonomous material search methods. (a) The conventional autonomous material search method used machine learning and material experiments/simulations. To optimize the target properties, conventional autonomous material search is performed using only the target properties and their material descriptors. (b) The proposed autonomous material search method. Target properties and nontarget properties are utilized in the autonomous material search.Sci. Technol. Adv. Mater. Meth. 3 (2023) 2                                                                                                                           J. HWANG AND Y. IWASAKIside. For transfer learning, we froze two input-side layers and fine-tuned one output-side layer. During training, we applied early stopping using a validation set comprising 20% randomly selected training data. In other words, parameters were updated through backpropagation using 80% of the data and the remaining 20% determined the optimal parameter selection for the final training results. We utilized the upper confidence bound as the acquisition function, where the predicted values and uncertainty are ensemble average and variance, respectively [18]. We set the exploration coefficient to 3.0. Further details regarding ENN settings and performance are provided in Supplementary Information S1 and S2.3. DemonstrationTo demonstrate the efficiency of ENN-based autonomous search, we applied our methodology to a dataset of ternary magnetic alloys with a B2 structure. The B2 crystal structure is shown in Figure 2(a), comprising two different sites denoted as X and Y in red and blue, respectively. The X and Y sites are populated by elements from the red and blue regions of the periodic table, respectively, as illustrated in Figure 2(b). X sites are populated by one or two elements selected from {Li, Mg, Ti, V, Mn, Fe, Co, Ni, Cu, Ru, Rh, Pd, Ag, Cd, Ir, Pt, and Au}, while Y sites are populated by one or two elements selected from {Li, Be, B, Mg, Al, Si, Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, Ga, Ge, As, Y, Zr, Nb, Mo, Ru, Ag, In, Sn, Sb, Hf, Ta, W, Pt, and Pb}. One of these sites is disordered and contains two elements, whereas the other site contains one element. Consequently, the ternary alloys are represented as X11� xX2xY1 and X1Y11� xY2x , where x denotes the composition in 0.2 increments. For instance, when two elements (e.g. Fe and Rh) are selected and disordered for the X site, only one element (e.g. Al) is selected for the Y site, the resulting ternary alloys are represented as Fe1-xRhxAl. This study only considers ternary alloys containing either Fe, Co, or Ni elements and thus comprises a dataset of approximately 17,000 magnetic ternary alloys.The magnetic moment M, spin polarization S, and Curie temperature Tc were obtained through density functional theory (DFT) calculations based on Green’s function. For this purpose, we used the Korringa – Kohn–Rostoker coherent-potential approximation (KKR-CPA) method implemented in the AkaiKKR software [19–24]. Additional information regarding the dataset and DFT calculations can be found in Supplemental Information S3.For material descriptors, we used molar ratios to represent the combination of elements and their amounts. To provide more detailed material description, we adopted the elemental fingerprints proposed by Hwang et al. [25]. The resulting 123 material descriptors and their corresponding details are listed in Supplemental Information S4.To simulate a scenario where there are abundant M and S data but no Tc data in the database, we conducted an autonomous material search. Our aim was to maximize the Tc as the target property, while leveraging information on the nontarget properties M and S.4. Results and discussionWe conducted several autonomous material searches to optimize the Tc. Figure 3(a) shows the results of one such search, which aimed to maximize the Tc using an ENN model without transfer learning (model N). The blue dots indicate the Tc of the material found for each search iteration, whereas the thick black line represents the maximum Tc achieved thus far. As the number of search iterations increased, the training data and accuracy of the ENN model improved, allowing the discovery of materials with a higher Tc.Figure 2. Material data used in this demonstration. (a) The B2 crystal structure, which comprises two different sites denoted as X and Y in red and blue, respectively. (b) A periodic table representing the elements of X and Y sites. One of these sites is disordered and contains two elements, whereas the other site remains fixed.Sci. Technol. Adv. Mater. Meth. 3 (2023) 3                                                                                                                           J. HWANG AND Y. IWASAKIFigure 3(b) shows the maximum Tc achieved till date in four different autonomous material search cases: ENN with transfer learning from M (model M), S (model S), and both M and S (model SM), and no transfer learning (model N). We conducted 10 searches each using models M, S, SM, and N, and the results are represented as thin blue, green, red, and black lines, respectively. The thick lines indicate the average of these 10 results. Models M and SM were able to discover materials with a higher Tc more rapidly than model N.In material development, it is often advantageous to propose multiple high-performance materials rather than just a top one. Figure 3(c) shows the number of autonomous material searches versus the count of the top 20 high Tc materials discovered. On average, model N identified 4.2 high Tc (Top-20) candidate materials in 30 observations, whereas models M and SM found an average of 8.8 and 8.9 high Tc candidates in 30 observations. Our findings demonstrate that transfer learning, where nontarget property information is used, can result in efficient autonomous material searches. Additionally, this approach is especially useful when only a limited amount of target property data is available.However, as shown in Figure 3(b,c), the discovery speed of model S was observed to be slower than that of model N, indicating that transfer learning is not always an efficient means for autonomous material search. To determine the optimal use of transfer learning in this context, we evaluated the predictive performance of the ENN model after transfer learning. Figure 4 shows the mean absolute error (MAE) as a function of observed data for the four transfer learning types, represented by the blue (M), green (S), red (SM), and black (N) lines, respectively.Figure 3. (a) The results of autonomous material search to maximize Tc using an ENN model without transfer learning (model N). The blue dots indicate the Tc of the material found for each search iteration, whereas the thick black line shows the maximum Tc achieved till date. (b) The maximum Tc achieved till date in four different autonomous material search cases (models M, S, MS, and N). The thick lines show the average of 10 results. (c) The count of the top 20 high Tc materials discovered by the four autonomous material searches.Figure 4. The mean absolute error (MAE) as a function of observed data for the four transfer learning types: models M, S, SM, and N. The predictive performance of model S is inferior to that of model N because of negative transfer.Sci. Technol. Adv. Mater. Meth. 3 (2023) 4                                                                                                                           J. HWANG AND Y. IWASAKIWe found that when the observed data are scarce, the predictive performances of models M and SM surpass that of model N, leading to an efficient autonomous material search, as shown in Figure 3(b,c). However, the predictive performance of model S is inferior to that of model N, and leveraging the source domain data/ knowledge in transfer learning afforded a negative transfer, which decreases learning performance in the target domain [26]. As a result of low predictive performance owing to negative transfer, autonomous material search utilizing model S exhibited low efficiency.To mitigate the potential for negative transfer during the autonomous materials search process, it is crucial to implement a verification step. Principally, there exist two strategies that can be employed for verification.The first strategy relies on insights derived from materials science. This strategy essentially involves understanding the correlation between the target property and various nontarget properties from a materials science perspective, thereby enabling the selection of an appropriate transfer learning model. For instance, considering Curie temperature’s minimal correlation with spin polarization – which reflects the electron states only near the Fermi energy – it becomes clear that model S can be excluded from the candidate list of transfer learning models for the autonomous materials search. However, it is important to note that such determinations are not always feasible from a materials science viewpoint.In situations where the first strategy proves ineffective, we rely on the second strategy, grounded in data science. This involves the generation of multiple transfer learning models using data of various nontarget properties. Their predictive performances are then evaluated using the target property data, with the highest- performing model being selected for use in the autonomous materials search. In this case study, for instance, models N, S, M, and SM were constructed at each stage of the autonomous materials search. The model demonstrating the highest predictive performance with respect to the target property data (Tc) is selected in each iteration of the autonomous search. For instance, at the stage of the search process where the observed data was 30 as illustrated in Figure 4, model M —with the highest predictive performance – is selected. The dynamic selection of transfer learning models at each stage of the search process is anticipated to enhance the overall efficiency of the materials search.5. ConclusionsWe proposed a novel approach for autonomous material search. This method optimizes target properties while leveraging transfer learning to transfer information on nontarget properties. Our results demonstrate the effectiveness of this method in the search for materials with a high Tc while utilizing information on the M and S of ternary magnetic alloys. The use of transfer learning enabled an efficient autonomous search, particularly in the early stage of the search process. However, we observed instances where negative transfer negatively impacted the search efficiency, highlighting the need for the dynamic selection of models based on their performance. Overall, our approach is widely applicable to diverse properties and material systems and can serve as a foundation for accelerating material development.AcknowledgementsWe thank K.S., Y.S, L.G., C.W., D.S., T.K., H.N., Y.G., and E.G. in the National Institute for Materials Science, M.K in Tokyo University of Science, Y.I in Tsukuba University for valuable discussions. This study was supported by JST- CREST Grant No. JPMJCR21O1. This paper is partially based on results obtained from a project, JPNP20003, commissioned by the New Energy and Industrial Technology Development Organization (NEDO).Disclosure statementNo potential conflict of interest was reported by the author(s).FundingThis study was supported by JST-CREST [Grant No. JPMJCR21O1]. This paper is partially based on results obtained from a project, JPNP20003, commissioned by the New Energy and Industrial Technology Development Organization (NEDO).ORCIDJaekyun Hwang http://orcid.org/0000-0002-1792-0044Yuma Iwasaki http://orcid.org/0000-0002-7117-277XData availability statementThe data that support the results reported herein and other findings of this study are available from the corresponding author upon reasonable request.References[1] Coley CW, Thomas DA, Lummiss JAM, et al. A robotic platform for flow synthesis of organic compounds informed by AI planning. Science. 2019;365 (6453). doi: 10.1126/science.aax1566[2] Burger B, Maffettone PM, Gusev VV, et al. A mobile robotic chemist. Nature. 2020;583:237–241. doi: 10. 1038/s41586-020-2442-2[3] Nikolaev P, Hooper D, Webber F, et al. Autonomy in materials research: a case study in carbon nanotube Sci. Technol. Adv. Mater. Meth. 3 (2023) 5                                                                                                                           J. HWANG AND Y. IWASAKIhttps://doi.org/10.1126/science.aax1566https://doi.org/10.1038/s41586-020-2442-2https://doi.org/10.1038/s41586-020-2442-2growth. npj Comput Mater. 2016;2(16031). doi: 10. 1038/npjcompumats.2016.31[4] Shimizu R, Kobayashi S, Watanabe Y, et al. Autonomous materials synthesis by machine learning and robotics. APL Mater. 2020;8(111110). doi: 10. 1063/5.0020370[5] Li Z, Najeeb MA, Alves L, et al. Robot-accelerated perovskite investigation and discovery. Chem Mater. 2020;32(13):5650–5663. doi: 10.1021/acs.chemmater. 0c01153[6] Roch AM, Hase F, Kreisbeck C, et al. ChemOS: orchestrating autonomous experimentation. Sci Rob. 2018;3(19). doi: 10.1126/scirobotics.aat5559[7] Attia PM, Grover A, Jin N, et al. Closed-loop optimization of fast-charging protocols for batteries with machine learning. Nature. 2020;578(7795):397–402. doi: 10.1038/s41586-020-1994-5[8] Granda JM, Donina L, Dragone V, et al. Controlling an organic synthesis robot with machine learning to search for new reactivity. Nature. 2018;559 (7714):377–381. doi: 10.1038/s41586-018-0307-8[9] Szymanski NJ, Zeng Y, Huo H, et al. Toward autonomous design and synthesis of novel inorganic materials. Mater Horiz. 2021;8(8):2169–2198. doi:  10.1039/D1MH00495F[10] Iwasaki Y, Sawada R, Saitoh E, et al. Machine learning autonomous identification of magnetic alloys beyond the Slater-Pauling limit. Commun Mater. 2021;2(31). doi: 10.1038/s43246-021-00135-0[11] Sawada R, Iwasaki Y, Ishida M. Boosting material modeling using game tree search. Phys Rev Mater. 2018;2(103802). doi: 10.1103/PhysRevMaterials.2. 103802[12] Seko A, Togo A, Hayashi H, et al. Prediction of low-thermal-conductivity compounds with first-principles anharmonic lattice-dynamics calculations and bayesian optimization. Phys Rev Lett. 2015;115(205901). doi: 10.1103/PhysRevLett.115. 205901[13] Jalem R, Kanamori K, Takeuchi I, et al. Bayesian- driven first-principles calculations for accelerating exploration of fast ion conductors for rechargeable battery application. Sci Rep. 2018;8(5845). doi: 10. 1038/s41598-018-23852-y[14] Iwasaki Y, Jaekyun H, Sakuraba Y, et al. Efficient autonomous material search method combining ab initio calculations, autoencoder, and multi-objective Bayesian optimization. Sci Technol Adv Mater. 2022;2(1):365–371. doi: 10.1080/27660400.2022. 2123263[15] Furuya D, Miyashita T, Miura Y, et al. Autonomous synthesis system integrating theoretical, informatics, and experimental approaches for large-magnetic-anisotropy materials. Sci Technol Adv Mater. 2022;2 (1):280–293. doi: 10.1080/27660400.2022.2094698[16] Lakshminarayanan B, Pritzel A, Blundell C. Simple and scalable predictive uncertainty estimation using deep ensembles. Adv Neural Inf Process Syst. 2017;30.[17] Yang Y, Lv H, Chen N. A survey on ensemble learning under the era of deep learning. Artif Intell Rev. 2022;56 (6):5545–5589. doi: 10.1007/s10462-022-10283-5[18] Auer P. Using confidence bounds for exploitation-exploration trade-off. J Mach Learn Res. 2002;3:397–422.[19] Akai H. Electronic structure Ni-Pd alloys calculated by the self-consistent KKR-CPA method. J Phys Soc Jpn. 1982;51:468–474. doi: 10.1143/JPSJ.51.468.[20] Khan SN, Staunton JB, Stocks GM. Statistical physics of multicomponent alloys using KKR-CPA. Phys Rev B. 2016;93(5):054206. doi: 10.1103/PhysRevB.93. 054206[21] Yang L, Liu B, Luo H, et al. Investigation of the site preference in Mn2RuSn using KKR-CPA-LDA calculation. J Magn Mater. 2015;382(15):247–251. doi: 10.1016/j.jmmm.2015.01.081[22] Akai H. Fast Korringa-Kohn-Rostoker coherent potential approximation and its application to FCC Ni-Fe systems. J Phys Condens Matter. 1989;1:8045–8063. doi: 10.1088/0953-8984/1/43/006.[23] AkaiKKR (machikaneyama). Ab-initio electronic-structure calculation code. http://kkr. issp.u-tokyo.ac.jp [24] Akai H, Dederichs PH. Local moment disorder in ferromagnetic alloys. Phys Rev B. 1993;47 (14):8739–8747. doi: 10.1103/PhysRevB.47.8739[25] Hwang J, Tanaka Y, Ishino S, et al. Prediction of viscosity behavior in oxide glass materials using cation fingerprints with artificial neural networks. Sci Technol Adv Mater. 2020;21(1):492–504. doi: 10. 1080/14686996.2020.1786856[26] Rosenstein MT, Marx Z, Kaelbling LP, et al. To transfer or not to transfer. NIPS 2005 Workshop Trans Learn. 2005;898(3):1–4.Sci. Technol. Adv. Mater. Meth. 3 (2023) 6                                                                                                                           J. HWANG AND Y. IWASAKIhttps://doi.org/10.1038/npjcompumats.2016.31https://doi.org/10.1038/npjcompumats.2016.31https://doi.org/10.1063/5.0020370https://doi.org/10.1063/5.0020370https://doi.org/10.1021/acs.chemmater.0c01153https://doi.org/10.1021/acs.chemmater.0c01153https://doi.org/10.1126/scirobotics.aat5559https://doi.org/10.1038/s41586-020-1994-5https://doi.org/10.1038/s41586-018-0307-8https://doi.org/10.1039/D1MH00495Fhttps://doi.org/10.1039/D1MH00495Fhttps://doi.org/10.1038/s43246-021-00135-0https://doi.org/10.1103/PhysRevMaterials.2.103802https://doi.org/10.1103/PhysRevMaterials.2.103802https://doi.org/10.1103/PhysRevLett.115.205901https://doi.org/10.1103/PhysRevLett.115.205901https://doi.org/10.1038/s41598-018-23852-yhttps://doi.org/10.1038/s41598-018-23852-yhttps://doi.org/10.1080/27660400.2022.2123263https://doi.org/10.1080/27660400.2022.2123263https://doi.org/10.1080/27660400.2022.2094698https://doi.org/10.1007/s10462-022-10283-5https://doi.org/10.1143/JPSJ.51.468https://doi.org/10.1103/PhysRevB.93.054206https://doi.org/10.1103/PhysRevB.93.054206https://doi.org/10.1016/j.jmmm.2015.01.081https://doi.org/10.1088/0953-8984/1/43/006http://kkr.issp.u-tokyo.ac.jphttp://kkr.issp.u-tokyo.ac.jphttps://doi.org/10.1103/PhysRevB.47.8739https://doi.org/10.1080/14686996.2020.1786856https://doi.org/10.1080/14686996.2020.1786856 Abstract 1. Introduction 2. Proposed method 3. Demonstration 4. Results and discussion 5. Conclusions Acknowledgements Disclosure statement Funding ORCID Data availability statement References