The present invention relates generally to metallic glass coating materials, such as alloys of Fe, Nb, Mo and B.
The omnipresence of structural components makes mechanical wear and tear, and the resulting performance and lifetime degradation, a topic of extreme importance. Manufacturing these components consumes a sizable portion of the energy and is a major contributor to carbon emission. Wear-resistant coatings on contact-components and cutting surfaces have the potential to not only improve mechanical performance by reducing surface deformation but also prolong operating life. Wear-resistant coatings can significantly impact the economy and the environment. For example, a 2017 report by Holmberg and Erdemir suggests that better surfaces on components can easily result in over trillion USD of savings annually.
Amorphous alloys, commonly called metallic glasses (MGs), provide a potential solution to this problem by promising to be better candidates for wear-resistant is coatings than their crystalline counterparts. The lack of crystalline long-range periodicity removes traditional dislocation-based deformation pathways, allowing higher hardness, strength, and elastic strain limit. Unfortunately, finding alloy compositions with a high glass-forming likelihood is challenging. Finding one with additional property constraints such as high wear-resistance is even more difficult. Additionally, large numbers of known MGs contain scarce or expensive elements (e.g., Pd, Au, Ir), making them poor candidates for inexpensive coatings, let alone bulk structural components.
A recent computational study suggests that in multiple principal element alloys (MPEAs), entropy begins to dominate enthalpic effects as the number of principal elements increases, and disorder becomes more prevalent. Although all known MGs are MPEAs, only a small fraction of disordered MPEAs is glass-forming. The challenge is that in a search space defined by 25 inexpensive, earth-friendly elements, there are 300 binary systems, 2,300 ternary systems, and 12,650 quaternary systems and each of these higher dimensional systems contains hundreds if not thousands of possible alloy compositions. Thus, there are tens of millions of possible compositions to search for a likely glass-former. A blind search of the vast combinatorial space for MGs is intractable by traditional methods. The quasi-empirical rules for glass-formation such as higher likelihood near eutectic points, thermodynamic models, geometric packing factors, and atomic number fractions, although informative, provide insufficient predictive power to design new MGs. Moreover, there is no defined way to incorporate information from previous successful and failed experiments to improve predictions. Consequently, many MGs discovered over the last 60 years were through serendipity, which is not a reliable design strategy.
In one aspect, the invention provides a new class of metallic alloys that are four times harder and nearly two times more wear-resistant than stainless steel. These alloys can be sputtered on to a surface at room temperature. It is made from non-toxic materials and is relatively cheap to apply and has comparable electrical conductivity to graphite.
We also outline a machine-learning guided procedure to further improve the performance of these alloys. Relying on machine learning (ML) predictions of MGs alone requires a highly precise model; however, incorporating high-throughput (HiTp) experiments into the search rapidly leads to higher performing materials even from moderately accurate models. Here, we exploit this synergy between ML predictions and HiTp experimentation to discover new hard and wear-resistant MGs in the Fe—Nb—B ternary material system. Several of the new alloys exhibit hardness greater than 25 GPa, which is over three times harder than hardened stainless steel and only surpassed by diamond and diamond-like-carbon. This ability to use less than perfect ML predictions to successfully guide HiTp experiments, demonstrated here, is especially important for searching the vast Multi-Principal-Element-Alloy combinatorial space, which is still poorly understood theoretically and sparsely explored experimentally.
These techniques would be attractive for commercial applications from cutting surfaces (knives to machine tools), to contact surfaces (especially electrical connectors and brushes in electrical motors).
is Advantages and improvements over existing methods, devices or materials is that the coatings of the present invention are cheaper, room temperature processable, and electrically conducting.
In one aspect, the invention provides a metallic glass coating material comprising an alloy of Fe, B, and one of the metals Nb, Mo, Zr, or W. The ratios of Fe, B, and the metal are predetermined using machine learning predictions and high-throughput experiments.
In another aspect, the invention provides a metallic glass coating material comprising an alloy of Fe, Nb, Mo and B, of the form Fex(Nb, Mo)yBz, where x is in the range 18-28, y is in the range 35-45, and z is in the range 32-42.
For example, the metallic glass coating material may be the alloy Fe23(Nb, Mo)40B37. The alloy may be doped with Zr and/or W, where the Zr and/or W comprises at most 10% of the alloy.
Lacking reliably predictive physiochemical models for discovering new metallic glasses with desirable properties, the inventors used machine learning (ML) to guide the search for a wear-resistant metallic glass. Specifically, an iterative approach combines ML with high throughput (HiTp) experiments to accelerate the discovery of new MGs by over 100-times. By training models on the weaknesses of predecessor models (boosting), every iteration resulted in improved glass-prediction accuracy for sputtered coatings. The model for glass-forming likelihood of arbitrary composition underwent two additional iterations of HiTp synthesis, testing, and retraining for improved accuracy. The next important step was developing a model to predict wear resistance. Although there are some computational reports relating the mechanical behavior of MGs to their free volume and local stresses (which activate interdependent shear transformation zones), there are no widely accepted physiochemical theories or clear empirical rules to predict hardness and wear-resistance of MGs. The classical theory of wear suggests hardness is closely correlated with wear resistance, but it is only one of the determinants. Most wear-resistant materials are hard, but not all hard materials are is wear-resistant. Nevertheless, rejecting materials with low hardness allows shrinking of the search space for wear-resistant materials.
The inventors developed a model for hardness to guide a search towards finding wear-resistant MGs. The hardness and glass-formability models were similar in structure but trained independently. To develop the hardness model, a dataset of hardness for 491 MGs was compiled by manually extracting results from 125 publications. This dataset was used to train a Random Forest (RF) machine-learning model, using the Matminer feature set. Because of the greatly sparser training set for hardness, the feature set was pruned of invariant features to mitigate over-fitting (see methods section for details). Using both the glass-formability and MG hardness models, we extracted 4 desired characteristics for over 2000 ternary alloy systems composed of 25 earth-friendly elements. The 4 desired characteristics we picked were: the fraction of the ternaries predicted to be glass-forming, the fraction of the glass-forming region in each ternary yet to be explored, and the highest and the lowest predicted hardness for a ternary composition.
We synthesized large regions of Fe—Nb—B ternary composition space simultaneously using combinatorial co-sputtering and employed high throughput (HiTp) characterization to rapidly map both the structure (glass-formability) and mechanical properties. Besides increasing the search speed, and providing prediction error tolerance, the HiTp approach simultaneously maps both the positive (desired) and the negative (undesired) regions of properties and consequently provides a more comprehensive view of the compositional space. We divided the ternary space into 15 overlapping composition spreads. We deposited each spread as a 500-800 nm thick film on a Si-substrate using multisource magnetron co-sputtering and determined local composition using electron probe microanalysis (EPMA). HiTp x-ray diffraction (XRD) was performed at 44 discreet alloy compositions in each composition spread at a synchrotron beamline optimized for such measurements. The materials were classified into three categories by the width of the first sharp diffraction peak (FWHMFSDP) as crystal (FWHMFSDP<0.4 Å−1), glass (0.4<FWHMFSDP<0.57 Å−1), and highly amorphous alloys (FWHMFSDP>0.57 Å−1).
Since the discovery in the 1970s of metallic glasses in the vicinity of Fe80B20, Fe-based alloys have set the benchmark for commercially viable MGs for structural applications. Previous explorations of Fe-based MG's have been confined to a small region of the Fe—Nb—B ternary space near the Fe-rich Fe—B binary leg. The highest hardness previously reported for this ternary is 16 GPa (Fe56Nb8B36) and in this work, we report a hardness of 16.5 GPa for an almost identical composition (Fe57Nb7B36). Higher hardness values were reported for Co-based MGs in 2011 which were quickly surpassed by W-based MGs in 2013 (H˜ 24 GPa and H/E ˜0.07). More recently, 8 element Fe-based MG systems (SAM1651 and SAM2X5) with a hardness of 15.7 GPa and 16.3 GPa respectively, have emerged as potential candidates for naval and nuclear storage applications. Many of the alloys reported here are better than metallic glasses in commercial use today.
The ML-HiTp approach used in this study outlines a path for higher performance alloys in higher dimensions. The HiTp experiments allow us to map large swaths of the composition space simultaneously. The HiTp flood-light searches highlight trends that are often missed in a smaller and less comprehensive one-alloy-at-a-time measurement approach. One such trend becomes evident in a comparison of
Our discovery also highlights that many materials with outstanding properties often lie in composition spaces neighboring the ones that have been incrementally explored for decades. However, finding better performers in compositions “just around the corner” is not trivial, especially when the dimensionality of the search-space increases, as illustrated by the example of Nb doping of Fe—B alloys. Researchers have explored doping Fe—B MGs with Nb and other second-row TMs for decades, and in so doing have observed improved GFA but crystallization and degradation of mechanical performance sets in beyond a few at % Nb. Our results show that simply increasing Nb doping does not lead to better or harder glasses. Optimizing alloy compositions in the Fe—Nb—B ternary space requires adjusting the concentration of all three elements, and finding these multi-dimensional paths in the MPEA space without HiTp experimentation is nearly impossible.
Ultimately, the synergy between HiTp experimentation and ML enabled the discovery of these exceptional alloys within a sparsely explored combinatorial space. Employing one method without the other would have taken a prohibitive amount of time. In the absence of well-established physiochemical theories or reliable empirical rules for wear-resistant MGs, exploring the MPEA space with HiTp experimentation alone at an unsustainable rate of one ternary per day every day would take 10 years. Guidance from ML models is essential in searching the vast MPEA space. Similarly, ML models trained on sparsely populated datasets such as those currently available for MPEA are unlikely to produce highly accurate predictions. However, the wide search swaths of HiTp experimentation are error-tolerant and result in new discoveries when ML models are only moderately accurate. Even mediocre models point roughly in the right direction, and more importantly often tell where not to look. Combining this guidance with search-light sweeps of HiTp experimentation is a very efficient means to search complex unknown spaces. Moreover, as was shown in our previous work on MGs, iteration of ML predictions with HiTp experiments rapidly improves the accuracy of subsequent generations of predictions. Mediocre models become better, as the quantity and diversity of the training data increases. The iteration of ML-guided HiTp explorations will guide the search for wear resistant MGs through increasingly higher dimensional MPEA space faster and with progressively better accuracy and precision.
Given the limited size and diversity of the hardness training set, it is not surprising that ML predictions underestimate our observations considerably. However, this underestimation also reflects the tendency of Random Forest Regressors (RFR) to predict values within the bounds of its training set. They are better at interpolation than extrapolation and always underpredict record-breaking findings. The observations reported here match the RF prediction for low to moderate hardness well, but the discrepancy increases as the observed hardness increases (as shown in SM). The conservative nature of RFR predictions makes it a safer guide for searching experimentally expensive spaces, as there is a reasonable likelihood that the search will at minimum result in discoveries matching the predictions and pose a chance of exceeding predictions as was, fortunately, the case here.
In conclusion, here we report the discovery of earth-friendly MGs in the Fe—Nb—B ternary material system attractive for commercially viable wear-resistant coatings. We discovered several hard and wear-resistant compositions that exhibit hardness greater than 25 GPa, which is >3× larger than hardened stainless steel and wear resistance comparable to best commercial nitride coatings. This study also illustrates how machine-learning-guided high throughput experimentation can improve on best-known materials and quickly lead to superior ones in highly complex composition-processing spaces. The proposed approach is not limited to only finding wear-resistant MG compositions but can be readily applied to searching the complex composition spaces and synthesis pathways leading to other MG functionalities or properties such as high glass transition temperature, or exploration of crystallization kinetics. The approach can even be adapted to rapidly find new catalysts, thermoelectric, and other desired materials in the vast MPEA is design space.
Although the above examples and discussion focus on the Fe—Nb—B material system, investigations by the inventors show that certain closely related material systems have similar properties. In particular, a metallic glass coating material with desirable hardness and wear-resistance may also be composed of an alloy of Fe, B, and Mo, Zr, or W. Substituting Mo for Nb is expected to result in nearly equivalent properties. An alloy may also contain a combination of Mo and Nb with Fe and B. In particular, a metallic glass coating material comprising an alloy of Fe, Nb, Mo and B, of the form Fex(Nb, Mo)yBz, where x is in the range 18-28, y is in the range 35-45, and z is in the range 32-42, as illustrated by the region 400 in
The Fe—Nb—B system has been extensively explored for last two decades, but region 400 with far superior properties was completely overlooked, with most of the exploration focused on incremental improvement of region 200 (
Methods
Database Creation. To create a mechanical property database, we assembled ˜430 publicly available metallic glass manuscripts. Out of 430 manuscripts, the hardness training set was created based on 130 manuscripts that contain 491 hardness datapoints. The training data for the glass-forming likelihood contained 6139 observations distributed over 313 systems. These databases are available through the Materials Data Facility.
Machine Learning.
We trained two separate Random Forest models to predict glass formability and hardness in similar ways. Compositions from each dataset were given a physics-based feature set provided by the Matminer package. The model hyperparameters were optimized on this feature set using 10-fold cross-validation, and the best performing model was trained on the entire available dataset. In particular, the maximum number of features and number of estimators per split was set to 10 features, 256 estimators for the metallic glass model and 12 features, and 100 estimators for the hardness model. Invariant features were pruned from the hardness feature set before training to improve performance and prevent overfitting, while the model for glass-forming likelihood (GFL) via sputtering used a full feature set. However, the GFL model for this work went through two additional iterations of experimentation, testing, and retraining. All models were implemented in Python using scikit-learn.
Predictions were made using these models on a grid of 1326 compositions per ternary system. These compositions were equally spaced across a ternary phase diagram. Over 2000 ternary systems were predicted and ranked by both hardness and glass-forming probability.
Sample preparation. The Fe—Nb—B thin-film compositional spreads were deposited from 2-inch Fe, Nb, and B targets on unheated 2×2-inch Si-wafers using combinatorial magnetron co-sputtering instrument with <1×10−6 torr base pressure-filled with Ar (10 mTorr). The composition spread was controlled by applying different gun powers in the 20-80 W range, and the thickness of the deposited film varied in the 500-800 nm range.
Composition Measurements. The composition map was determined by wavelength dispersive spectroscopy (WDS) analysis in an electron probe microanalyzer (EPMA) JXA 8900R Microprobe, with an acceleration voltage of 15 kV. Standardization of references was carried out with pure metal references and compositions were determined to be within an experimental error of <0.3 at %. WDS was selected over energy dispersive spectroscopy (EDS) due to its higher accuracy and precision in quantifying elemental content via better energy resolution from peak/background ratio.
Structural characterization. The structural characterization was performed at beamline 1-5 at Stanford Synchrotron Radiation Lightsource (SSRL) using the two-dimensional XRD (MarCCD 165) detector. The XRD patterns were collected with 15.5-keV energy X-rays. To minimize diffraction from the silicon substrate, the samples were aligned to the incident beam at a grazing angle of 2°. The grazing incidence geometry resulted in an approximate 3-mm probe footprint on the samples. A LaB6 powder diffraction pattern was used to extract diffraction geometric parameters, for instance, direct beam position, tilting, rotation, and sample-to-detector distance. These parameters were used to transform initial 2-D raw images to Q and χ diffraction coordinate and then into 1D diffraction patterns by integrating and normalizing over the x angle by using custom python scripts.
Mechanical characterization and analysis. Experiments were performed on a Hysitron (now Bruker) TI 950 with NanoDMA III at a frequency of 200 Hz with a maximum load of 10 mN to measure hardness (H) and modulus (E) at a constant strain rate of 0.12 s−1 with a diamond, Berkovich-geometry nanoindenter tip, calibrated on fused quartz according to standard procedures. Thickness for each composition was measured by contact profilometry and used for subsequent data correction to remove the substrate compliance effect. The Li & Vlassak method was used for substrate correction of both Er and H and was implemented with some computational differences to streamline for parallel data processing. Four indentation tests were conducted on each composition and indentations were spaced at least 15 microns apart to ensure minimal overlap of indentation stress fields. We have used best practices for calibration of load frame compliance and tip area function, including using standard calibration materials, routine calibration checks throughout testing, and replacement of the tip when tip blunting causes insufficient plastic deformation in calibration standards at the indentation depths of interest. Calibration checks were performed daily in the study and the tip was replaced with a new, low radius of curvature indenter tip multiple times throughout the project, owing to tip blunting caused by repeatedly indenting materials with H>15 GPa. Elastic modulus reported here was Er, the substrate-effect-corrected reduced modulus, which is defined as:
This application claims priority from U.S. Provisional Patent Application 63/223,897 filed Jul. 20, 2021, which is incorporated herein by reference.
This invention was made with Government support under contract DE-ACO2-76SF00515 awarded by the Department of Energy. The Government has certain rights in the invention.
Number | Date | Country | |
---|---|---|---|
63223897 | Jul 2021 | US |