Combinatorial Materials Science for Energy Applications

R. Bruce van Dover

Department of Materials Science and Engineering, Cornell University

Materials are the fundamental basis for solutions to the most pressing issues in energy generation, transport, and utilization, as well as more general issues in sustainability. In many cases, long-term solutions to these problems will depend on breakthrough innovations in materials. As emphasized by the U.S. Department of Energy panel on “New Science for a Secure and Sustainable Energy Future”:

“…existing energy approaches-even with improvements from advanced engineering and improved technology based on known concepts will not be enough to secure our energy future. Instead, meeting the challenge will require new technologies for producing, storing and using energy with performance levels far beyond what is now possible. Such technologies spring from scientific breakthroughs in new materials and chemical processes that govern the transfer of energy between light, electricity and chemical fuels.”1

Can there be an effective strategy for finding breakthrough materials, since they are, by definition, unpredictable?

One answer is found in Combinatorial Materials Science techniques, which represent a powerful approach to identifying new and unexpected materials. The high cost of single-sample synthesis/ characterization, along with the need for reduced research and development time are driving the materials community to explore high-throughput methodologies. Increasing the number of materials that are studied can increase the understanding of composition/ property relationships and the probability of a breakthrough materials discovery. To be effective, this method must balance the need for robust and authoritative measures of properties with the need for rapid characterization.

Figure 1 delineates some of the considerations that must be met for the high-throughput synthesis/characterization approach to be viable. It is not necessary that materials made by the chosen high-throughput synthesis technique be identical to those made in one-off experiments, but it is necessary that the resulting properties are similar. Materials properties depend on morphology, microstructure and other processdependent variables. Thin films are often particularly convenient to work with and sometimes have properties that differ significantly from their bulk counterparts. Fortunately, it is found, in many cases, that compositional trends observed in thin films closely approximate those of the bulk.

Is the combinatorial/high-throughput approach the right choice to find a new material that can solve an outstanding problem? This decision tree illustrates some of the key questions that should be asked. When the right conditions are met, highthroughput experiments can be impressively effective.

Figure 1.Is the combinatorial/high-throughput approach the right choice to find a new material that can solve an outstanding problem? This decision tree illustrates some of the key questions that should be asked. When the right conditions are met, high throughput experiments can be impressively effective.

Many methods for high-throughput synthesis have been developed,2,3 but the Codeposited Composition Spread (CCS) technique has proven to be an especially versatile method for forming a wide range of compositions in a single experiment. In this method, thin films are deposited by physical vapor deposition on a substrate simultaneously, from two or more sources that are spatially separated and chemically distinct, producing a film with an inherent composition gradient and intimate mixing of the constituents. With three sources, an entire ternary phase diagram may be produced in a single experiment.

Composition spreads may also be synthesized using a traveling shutter4 or shaped mask5 to create a film with a thickness gradient (wedge). A composition gradient can then be obtained by rotating the sample with respect to the shutter and depositing a new overlapping wedge of a second or third film. Atomic mixing is achieved by depositing many wedges each of submonolayer thickness. This approach has the advantage of the composition/position dependence being well-defined (ignoring resputtering effects), although surface reorganization during the short time between wedge depositions can lead to adventitious artifacts.

CCS synthesis is distinct from conventional “combinatorial chemistry” as employed, for example, in drug discovery6 as well as the discrete combinatorial synthesis (DCS) approach pioneered by Xiang and Schultz.7 In the latter, mixtures of inorganic components are created by sequential deposition of discrete layers of precursors followed by moderate- or high-temperature diffusion and reaction steps. An important advantage of the DCS technique is that arbitrary compositions with a large number of constituents can be prepared as desired. A key advantage of the composition-spread approach is the opportunity to prepare materials with no subsequent processing, which means that low-temperature or metastable phases may be prepared. The CCS approach also allows properties to be determined with very fine composition resolution, often limited by the resolution of the property measurement itself. For example, it is generally straightforward to sample composition space at 1 mol% intervals, equivalent to investigating thousands of materials in a single experiment. The CCS technique has been used to create both alloys and compounds in chemical systems, including intermetallics, nitrides, oxides, and carbides. A large number of thin film deposition techniques are available for synthesis of composition spread thin films, including evaporation,8,9 sputtering,10-12 pulsed-laser deposition (PLD),13 chemical vapor deposition,14 and cold-plasma processing,15 among others. Of these, sputtering offers a unique combination of advantages:

  • In most cases sputtering rates are constant and reproducible16 when powered by modern dc, rf, or pulsed-dc regulated power supplies, allowing composition spreads with a targeted range of compositions to be synthesized.

  • When magnetron sputter guns are used for sputtering there is little interaction between the sources. That is, the deposition rate from one gun is independent of the operation of the other guns. As a result the composition profile (composition of the codeposited film as a function of position on the substrate) can be predicted quantitatively once the deposition rate for each individual source is determined. However, in some circumstances the magnetron gun may interact significantly with the growing films, leading to resputtering and deviation of the composition from that predicted by simple summation of independent deposition rates. This effect can be predicted quantitatively using known sputter yields and the conditions of the sputter system, leading to accurate adjustments and an excellent quantitative prediction of the composition.

  • The composition gradients achieved by codeposition (using two-inchdiameter sources) are typically about 1 atomic percent per mm, a convenient scale that allows many measurement techniques to resolve trends with about 1 atomic percent resolution.

  • Sputtering is convenient for deposition of metals, oxides (by reactive sputtering in an oxygen-containing ambient), nitrides (sputtering in Ar-N2), carbides (sputtering Ar-CH4), and mixed anion (e.g., oxynitride) systems. Targets of most metals and many unary and binary oxides are readily available, along with many other common materials.

Sputter codeposition also has inherent disadvantages, including:

  • It is not simple to change the composition gradient. The gradient is affected by changing the total inert gas pressure in the system, but many other aspects of sputtering change as well, confounding systematic studies.

  • Targets of the alkali metals, alkaline earth metals, and the lighter lanthanide metals are problematic, due to rapid oxidation/hydrolysis/ hydration of the metal target in air.

Determining the structure of the materials formed in a composition spread is important for understanding composition/property relations. Thin films are well suited for X-ray diffraction studies for phase identification. Using automated data acquisition, in either a conventional stationary-anode diffractometer17 or a synchrotron-based system,18 hundreds of diffraction patterns can be acquired on a single composition-spread substrate. An unsolved challenge is to develop automated techniques for the identification of unique diffraction patterns and clustering these into contiguous phase fields corresponding to regions (compositions) of the CCS film. Some progress has been made in this respect,17,19 but a fully robust algorithm has yet to be developed.

Fuel Cell Studies

Combinatorial techniques have been used in a wide range of energyrelated studies, both to identify new higher-performance materials and to elucidate the composition-dependent properties of known systems. One challenge for which the CCS technique has been found to be particularly well-suited is that of identifying superior electrocatalysts for Polymer Electrolyte Membrane (PEM) fuel cells.20 Catalysis is a complex phenomenon that is essentially impossible to predict with confidence; incremental improvements may be obtained using rational design, but the only way to identify new active compositions is by screening a large number of compositions. Mallouck and coworkers20 developed a viable technique for qualitative optical screening of catalytic activity, which has been further developed into a semiquantitative method.21 Alternative methods for identifying active materials include multiple independent electrode arrangements22 and scanning electrochemical microscopy.23

As an illustrative example, the CCS technique was recently used to gain insight into the most catalytically active material in the Pt-Ta and related systems. After exploring hundreds of chemical systems, the Pt-Ta system was identified as showing interesting activity, and it was selected for close study of the onset potential for methanol oxidation determined by optical fluorescence as a function of position on the binary composition spread. X-ray diffraction data was collected in a high-throughput automated experiment using a synchrotron beamline and used to identify phase fields, as illustrated in Figure 2. It was immediately clear that the best catalyst behavior (lowest value for the half-wave potential E1/2) is strongly correlated with the presence of the orthorhombic Pt2Ta structure, and that the activity is optimized at the approximately stoichiometric composition Pt0.71Ta0.29.24

The catalytic activity of Pt-Ta compositions is measured using a fluorescence test: a lower half-wave potential, E1/2, implies greater activity for methanol oxidation. This plot shows that strong activity is associated with the presence of the orthorhombic Pt2Ta phase, is best in the single phase region (18 -35 at% Ta), and is optimized at the composition Pt0.71Ta.29, close to the stoichiometric value.24

Figure 2.The catalytic activity of Pt-Ta compositions is measured using a fluorescence test: a lower half-wave potential, E1/2, implies greater activity for methanol oxidation. This plot shows that strong activity is associated with the presence of the orthorhombic Pt2Ta phase, is best in the single phase region (18 -35 at% Ta), and is optimized at the composition Pt0.71Ta.29, close to the stoichiometric value.24

Figure 2 shows that the fine compositional resolution offered by the CCS technique permits two important conclusions about the internal consistency of the data. First, the close agreement between values at adjacent compositions indicates that random variations in the measurement are small compared to the overall trend. Second, the smooth trend with composition in the Pt2Ta phase field allows the optimum composition to be identified with confidence. Internal consistency does not prove that the data is accurate in composition or in E1/2 value; the absolute accuracy must be validated by comparison with detailed one-off studies. The key advantage of the combinatorial method is that it allows the researcher to efficiently identify particular compositions for further study based on data rather than on conservative or radical speculation.

Predicting catalytic properties is not reliable—neither from first principles nor from accumulated experience—so catalyst development has always relied on an empirical approach. High-throughput studies have proven especially useful in facilitating rapid optimization of catalyst composition using factorial designs as well as more exotic approaches, such as genetic algorithms.25 Solution-based synthesis is often employed for optimization because it most closely approximates the processing of realistic formulations. The high-throughput approach has been applied successfully to primary screening and optimization for a wide range of catalyst functions, including polymerization catalysts, enantioselective catalysts, oxidation catalysts, reduction catalysts, dehydrogenation catalysts, and many others.

Transparent Conductivity Oxides

Another example of the usefulness of the combinatorial materials science approach is provided in the study of transparent conducting oxides-materials that serve a wide range of energy-related optoelectronic functions, from low-emissivity window coatings to frontside current collectors in photovoltaics.26,27 Established n-type transparent oxides, such as In1-xSnxOy, SnO2:F, and ZnO:Al, offer adequate performance, but with concomitant drawbacks, such as low thermal stability or high cost. More complex materials based on Ga2O3 (215066), SnO2 (518174), ZnO (204951), CdO (202894), and In2O3 (203424) and their multinary mixtures offer the prospect of improved performance and have been the subject of high-throughput as well as conventional studies.

ZnO is an inexpensive wide-bandgap semiconductor with excellent optical transmission. Highly conductive n-type ZnO has been achieved through doping with Al, In, or Ga. The electrical properties of ZnO are highly dependent on native point defects such as oxygen vacancies, zinc interstitials, and hydrogen, all of which act as electron donors. In order to obtain high conductivity ZnO, these point defects must be deliberately induced by doping. The highest conductivity is achieved after annealing in reducing conditions, while exposure to oxidizing conditions (typically air at moderate temperatures) degrades the conductivity significantly. Extensive studies on the substitution of Al3+ on the Zn2+ site have shown that doping at the level of a few atomic percent leads to a moderately high conductivity that is fairly stable. Doping with In3+ has a similar effect. Since Al3+ has an ionic radius 32% smaller than that of Zn2+, and In3+ has an ionic radius 9% larger, it is not surprising to find that the introduction of either impurity degrades the electron mobility of ZnO. A composition spread study of the conductivity, mobility, and crystal structure of the (Zn, Al, In, O) system allowed the effect of these impurities to be clarified. Figure 3 shows that codoping with both Al and In degraded mobility less than doping with either element alone.28 The lattice strain inferred from x-ray diffraction shows a similar trend; the average lattice constant matches that of undoped ZnO for materials with the highest mobility. The conductivity is also maximum for this condition. Use of a composition spread sample allowed the experiment to be executed without the confounding associated with typical run-to-run variations that accompany one-off studies, thereby allowing robust conclusions to be drawn.

ZnO with Al or In forms a transparent conductor. Compared to Zn2+, the Al3+ ion is much smaller while In3+ is much larger. Co-doping yields an “average” dopant size that is a better fit in the ZnO crystal, reducing scattering. A Zn-Al-In-oxide composition spread allows the effect of varying levels of both dopants to be measured. The binary spreads Zn-Al-oxide and Zn-In-oxide are prepared in separate experiments under nominally identical conditions. For a given overall doping level, the highest mobility is observed in co-doped compositions.

Figure 3.ZnO with Al or In forms a transparent conductor. Compared to Zn2+, the Al3+ ion is much smaller while In3+ is much larger. Co-doping yields an “average” dopant size that is a better fit in the ZnO crystal, reducing scattering. A Zn-Al-In-oxide composition spread allows the effect of varying levels of both dopants to be measured. The binary spreads Zn-Al-oxide and Zn-In-oxide are prepared in separate experiments under nominally identical conditions. For a given overall doping level, the highest mobility is observed in co-doped compositions.28

Other Energy-related Studies

Combinatorial materials science techniques have great potential for many other studies of energy-related materials. Thermoelectrics offer the potential to revolutionize technologies such as high energy refrigeration or energy scavenging from low-grade heat sources, a prospect that has driven the search for a breakthrough in thermoelectric materials.29 Many of the ingredients needed for a successful combinatorial search are in place. There is a well-accepted Figure of Merit (FOM) for thermoelectric materials that can, in principle, be measured in a thin film geometry. Recent concepts regarding the factors that might lead to high FOM materials could provide guidance regarding materials systems worthy of investigation.30 Perhaps the most challenging aspect of a search for new thermoelectrics is the sensitivity of the thermoelectric FOM to doping, which implies that a suitable dopant introduced at the optimum concentration must be identified for any prospective candidate. This both increases the number of materials combinations dramatically and introduces the possibility of strong sensitivity to process conditions.

Another example is that of piezoelectric materials, which have a variety of commercial applications and are proving useful in harvesting lowgrade energy. Only a small number of piezoelectric materials are commonly used: others are known but barely characterized. Undoubtedly, many more have yet to be synthesized or characterized. Highlytextured crystalline thin films are needed for reliable characterization of piezoelectric properties, so microstructure control may be a key synthesis issue. Screening approaches should be straightforward, using optical techniques or atomic force microscope (AFM) cantilever probing. New phases, perhaps with a wider processing window or larger piezoelectric Figure of Merit, could be quite exciting, both from a fundamental point of view as well as from the perspective of commercial importance.

Outstanding Issues

While development of synthesis and characterization techniques is basic to the concept of high-throughput experimentation, these steps lead to data and information, not knowledge or insight. In fact, the massive quantity of data, coupled with the modest level of accuracy or precision associated with the speed/quality tradeoff, quickly leads to a new challenge-how to make use of the data. Tracking only the latest topperforming material in a high-throughput investigation discounts the insights and new directions that could be gleaned from lesserperforming materials.

Efficient methods are urgently needed for transforming an overwhelmingly data-rich environment into actionable insights regarding structure-processing-property associations and relationships in materials. Machine learning and statistical approaches, such as cluster analysis and multivariate regression, are important components for this task, but they do not incorporate the relationships that are inherent in the physics and chemistry of materials. Thus, perhaps the most essential difference between conventional research and the combinatorial materials science approach is in the role of informatics—the processing, management and retrieval of information.


The high-throughput approach to inorganic materials discovery improves the likelihood of discovering new materials with useful properties because it dramatically lowers the cost-in terms of financial resources, human effort, and time-of examining unexplored regions of composition space, including regions that might be avoided as unlikely candidates for an expensinve one-off study. Experience has shown that high-throughput screening can be a useful tool for solving real-world problems in materials discovery if three broad criteria are met: In general, the approach is suitable for a well-defined problem for which samples can be prepared using parallelized synthesis and evaluated using a suitable high-throughput screen. High-throughput techniques are also valuable for investigations of known materials systems, where the goal is to elucidate the composition dependence of materials properties. The codeposited composition spread technique has proven particularly effective for exploring energy-related materials in a wide range of investigations.



Hemminger J, Crabtree G, Kastner M. 2008. New science for a secure and sustainable energy future,. Office of Basic Energy Sciences, Department of Energy
Xiang X, Takeuchi I. Combinatorial Materials Synthesis.
Beach JW. 2001. High Throughput Synthesis. Principles and Practice Edited by Irving Sucholeiki. Marcel Dekker, New York. 2001. xxi + 366 pp. 15.5 × 23.5 cm. ISBN 0-8247-0256-5. $175.00.. J. Med. Chem.. 44(20):3336-3337.
Xiang X. 1999. Bull. Am. Phys. Soc. 44103.
Dahn JR, Trussler S, Hatchard TD, Bonakdarpour A, Mueller-Neuhaus JR, Hewitt KC, Fleischauer M. 2002. Economical Sputtering System To Produce Large-Size Composition-Spread Libraries Having Linear and Orthogonal Stoichiometry Variations. Chem. Mater.. 14(8):3519-3523.
Schmitz DF. 1995. Review. 64(2):317-318.
Xiang X-, Sun X, Briceno G, Lou Y, Wang K, Chang H, Wallace-Freedman WG, Chen S, Schultz PG. 1995. A Combinatorial Approach to Materials Discovery. Science. 268(5218):1738-1740.
Kennedy K, Stefansky T, Davy G, Zackay VF, Parker ER. 1965. Rapid Method for Determining Ternary?Alloy Phase Diagrams. Journal of Applied Physics. 36(12):3808-3810.
Hammond RH, Ralls KM, Meyer CH, Snowden DP, Kelly GM, Pereue JH. 1972. Superconducting Properties of Nb3 (Al, Ge) Prepared by Vacuum Vapor Deposition of the Elements. Journal of Applied Physics. 43(5):2407-2413.
Hanak JJ. 1970. The ?multiple-sample concept? in materials research: Synthesis, compositional analysis and testing of entire multicomponent systems. J Mater Sci. 5(11):964-971.
van Dover RB, Hong M, Gyorgy EM, Dillon JF, Albiston SD. 1985. Intrinsic anisotropy of Tb?Fe films prepared by magnetron Co sputtering. Journal of Applied Physics. 57(8):3897-3899.
van Dover RB, Schneemeyer LF, Fleming RM. 1998. Discovery of a useful thin-film dielectric using a composition-spread approach. Nature. 392(6672):162-164.
Lippmaa M, Koida T, Minami H, Jin Z, Kawasaki M, Koinuma H. 2002. Design of compact pulsed laser deposition chambers for the growth of combinatorial oxide thin film libraries. Applied Surface Science. 189(3-4):205-209.
Smith RC, Hoilien N, Roberts J, Campbell SA, Gladfelter WL. 2002. Combinatorial Chemical Vapor Deposition of Metal Dioxides Using Anhydrous Metal Nitrates. Chem. Mater.. 14(2):474-476.
Terajima T, Koinuma H. 2004. Development of a combinatorial atmospheric pressure cold plasma processor. Applied Surface Science. 223(1-3):259-263.
Gregoire JM, Dale D, Kazimirov A, DiSalvo FJ, Bruce van Dover R. 2010. Cosputtered composition-spread reproducibility established by high-throughput x-ray fluorescence. Journal of Vacuum Science & Technology A: Vacuum, Surfaces, and Films. 28(5):1279-1280.
Long CJ, Hattrick-Simpers J, Murakami M, Srivastava RC, Takeuchi I, Karen VL, Li X. 2007. Rapid structural mapping of ternary metallic alloy systems using the combinatorial approach and cluster analysis. Review of Scientific Instruments. 78(7):072217.
Gregoire JM, Dale D, Kazimirov A, DiSalvo FJ, van Dover RB. 2009. High energy x-ray diffraction/x-ray fluorescence spectroscopy for high-throughput analysis of composition spread thin films. Review of Scientific Instruments. 80(12):123905.
Long CJ, Bunker D, Li X, Karen VL, Takeuchi I. 2009. Rapid identification of structural phases in combinatorial thin-film libraries using x-ray diffraction and non-negative matrix factorization. Review of Scientific Instruments. 80(10):103902.
Reddington E. 1998. Combinatorial Electrochemistry: A Highly Parallel, Optical Screening Method for Discovery of Better Electrocatalysts. 280(5370):1735-1737.

Gregoire JM, Kostylev M, Tague ME, Mutolo PF, van Dover RB, DiSalvo FJ, Abrun?a HD. 2009. High-Throughput Evaluation of Dealloyed Pt?Zn Composition-Spread Thin Film for Methanol-Oxidation Catalysis. J. Electrochem. Soc.. 156(1):B160.
Strasser P, Fan Q, Devenney M, Weinberg WH, Liu P, Nørskov JK. 2003. High Throughput Experimental and Theoretical Predictive Screening of Materials ? A Comparative Study of Search Strategies for New Fuel Cell Anode Catalysts. J. Phys. Chem. B. 107(40):11013-11021.
Jayaraman S, Hillier AC. 2003. Screening the Reactivity of PtxRuyand PtxRuyMozCatalysts toward the Hydrogen Oxidation Reaction with the Scanning Electrochemical Microscope. J. Phys. Chem. B. 107(22):5221-5230.
Gregoire JM, Tague ME, Cahen S, Khan S, Abrun?a HD, DiSalvo FJ, van Dover RB. 2010. Improved Fuel Cell Oxidation Catalysis in Pt1?xTax?. Chem. Mater.. 22(3):1080-1087.
Potyrailo RA, Maier WF. Combinatorial and High-Throughput Discovery and Optimization of Catalysts and Materials.
Ginley DS, Bright C. 2000. Transparent Conducting Oxides. MRS Bull.. 25(8):15-18.
Perkins J, del Cueto J, Alleman J, Warmsingh C, Keyes B, Gedvilas L, Parilla P, To B, Readey D, Ginley D. 2002. Combinatorial studies of Zn-Al-O and Zn-Sn-O transparent conducting oxide thin films. Thin Solid Films. 411(1):152-160.
Kirby S, van Dover R. 2009. Improved conductivity of ZnO through codoping with In and Al. Thin Solid Films. 517(6):1958-1960.
DiSalvo FJ. 1999. Thermoelectric Cooling and Power Generation. 285(5428):703-706.
Mahan GD. 1989. Figure of merit for thermoelectrics. Journal of Applied Physics. 65(4):1578-1583.


This work was supported by the U.S. Department of Energy (DOE), Office of Science, Office of Basic Energy Sciences (under Award Number DE-FG02-07ER46440). The author would like to gratefully acknowledge formative discussions with Lynn Schneemeyer, Héctor Abruña, Francis DiSalvo and John Gregoire.

Social Media

LinkedIn icon
Twitter icon
Facebook Icon
Instagram Icon


Research. Development. Production.

We are a leading supplier to the global Life Science industry with solutions and services for research, biotechnology development and production, and pharmaceutical drug therapy development and production.

© 2021 Merck KGaA, Darmstadt, Germany and/or its affiliates. All Rights Reserved.

Reproduction of any materials from the site is strictly forbidden without permission.