This disclosure generally relates to electrocatalysts and, more particularly, to transition metal-doped platinum-based catalysts.
Proton-exchange membrane (PEM) fuel cells are desirable energy conversion devices for applications such as transportation vehicles and portable electronic devices, due to their high-energy density and low environmental impact in addition to being light-weight and affording low-temperature operation. PEM fuel cells operate based on reactions of a fuel (such as hydrogen or an alcohol) at an anode and an oxidant (molecular oxygen) at a cathode. Both cathode and anode reactions include catalysts to lower their electrochemical over-potential for high-voltage output, and so far, platinum (Pt) has been the leading choice. To fully realize the commercial viability of fuel cells, the following challenges should be addressed: the high cost of Pt, the sluggish kinetics of the oxygen reduction reaction (ORR), and the low durability of Pt-based catalysts.
It is against this background that a need arose to develop the embodiments described herein.
Alloying Pt with a secondary metal reduces the usage of scarce Pt metal while at the same time provides improved performance as compared with that of pure Pt in terms of activity. In particular, bimetallic platinum-nickel (Pt—Ni) nanostructures represent a class of electrocatalysts for ORR in fuel cells, but practical applications have been constrained by catalytic activity and durability. Although an increase in ORR activity is observed for Pt—Ni nanostructures, the activity as observed on bulk Pt3Ni(111) surface has not been matched, indicating room for further improvement. At the same time, a notable constraint of Pt—Ni nanostructures is their low durability. The Ni element in these nanostructures leaches away gradually under detrimental corrosive ORR conditions, resulting in rapid performance losses. Thus, Pt-based nanostructures with simultaneously high catalytic activity and high durability have remained a challenge.
Some embodiments of this disclosure are directed to surface-doped Pt3Ni nanostructures in the form octahedra supported on carbon, with dopants corresponding to transition metals, termed M-Pt3Ni/C, where M is vanadium (V), chromium (Cr), manganese (Mn), iron (Fe), cobalt (Co), molybdenum (Mo), tungsten (W), or rhenium (Re). In some embodiments, Mo—Pt3Ni/C exhibits particularly improved ORR performance, with a specific activity of about 10.3 mA/cm2 (or greater) and a mass activity of about 6.98 A/mgPt (or greater), which are about 81- and 73-fold enhancements compared with a Pt/C catalyst (about 0.127 mA/cm2 and about 0.096 A/mgPt). Without wishing to be bound by a particular theory, calculations indicate that Mo preferentially locates at subsurface positions near nanoparticle edges in vacuum and surface vertex/edge sites in oxidizing conditions, where Mo can enhance both the activity and the stability of the Pt3Ni catalyst. The surface doping approach can be applied to the rational design of catalysts and other materials with enhanced activity and durability, for applications such as fuel cells, batteries and chemical production.
Other aspects and embodiments of this disclosure are also contemplated. The foregoing summary and the following detailed description are not meant to restrict this disclosure to any particular embodiment but are merely meant to describe some embodiments of this disclosure.
For a better understanding of the nature and objects of some embodiments of this disclosure, reference should be made to the following detailed description taken in conjunction with the accompanying drawings.
Embodiments of this disclosure are directed to improved Pt-based electrocatalysts for ORR, exhibiting a combination of high activity and high stability. Some embodiments are directed to Pt—Ni nanostructures and, in particular, Pt3Ni-based nanostructures because Pt3Ni(111) surfaces can provide efficient catalysis for ORR. In some embodiments, challenges discussed above are addressed through a surface engineering approach based on control over dopant incorporation of various transition metals on surfaces of Pt3Ni nanostructures in the form octahedra supported on carbon, termed M-Pt3Ni/C, where M is V, Cr, Mn, Fe, Co, Mo, W, or Re. Owing to the efficient in-situ and bulky capping agent-free approach of some embodiments, the coupling of Pt3Ni nanostructures to a carbon support can remain strong. A seed-mediated growth strategy and a relatively slow continuous surface dopant infusion can provide a desirable growth condition, leading to a well-maintained particle size distribution of surface-doped Pt3Ni nanostructures. Resulting surface-doped Pt3Ni nanostructures can exhibit impressive activity in ORR, and their activity can be dopant-dependent. Of note, the resulting surface-doped Pt3Ni nanostructures can simultaneously satisfy an overall criteria of high specific activity, high mass activity, and suppressed degradation of performance. Without wishing to be bound by a particular theory, the presence of a transition metal dopant, such as an electropositive Mo in Mo—Pt3Ni/C, can stabilize an alloy composition of the catalyst by inhibiting against dissolution and diffusion through formation of strong bonds with the transition metal dopant, such as strong Mo—Pt and Mo—Ni bonds, and can shift oxygen binding energies to promote enhanced catalytic activity.
In some embodiments, a Pt-based electrocatalyst is a doped intermetallic alloy of Pt and at least one secondary metal having a chemical composition that can be represented by the formula Mz—PtxNiy, where any one or any combination of two or more of the following applies: (1) Pt represents platinum as a primary metal; (2) Ni represents nickel as a secondary metal; (3) M represents at least one metal as a dopant and with M being different from Pt and Ni, such as where M is at least one transition metal selected from Group 3, Group 4, Group 5, Group 6, Group 7, Group 8, Group 9, Group 10, Group 11, and Group 12 of the Periodic Table and with M being different from Pt and Ni; (4) x represents a molar content of Pt, y represents a molar content of Ni, z represents a molar content of M, with x>y and x>z and also, in some embodiments, y>z; (5) x has a non-zero value in a range of about 51 to about 95, such as about 60 to about 90, about 68 to about 82, about 70 to about 78, or about 72 to about 76; (6) y has a non-zero value in a range of about 5 to about 49, such as about 10 to about 40, about 18 to about 32, about 20 to about 30, about 22 to about 28, or about 23 to about 27; (7) z has a non-zero value in a range of 0 to about 8, such as about 0.1 to about 8, about 0.5 to about 5, about 0.5 to about 3, or about 0.5 to about 2.5; and (8) subject to the condition that x+y+z=100 (or 100%).
In some embodiments, a ratio of x to y is about 3, and z has a non-zero value in a range of about 0.5 to about 3 or about 0.5 to about 2.5.
In some embodiments, M is at least one transition metal selected from Group 5, Group 6, Group 7, Group 8, and Group 9 of the Periodic Table. In some embodiments, M is one or more of V, Cr, Mn, Co, Mo, W, and Re. It is contemplated that a Pt-based electrocatalyst can be doped with combinations of two or more different doping elements, such as two or more transition metals selected from Group 3, Group 4, Group 5, Group 6, Group 7, Group 8, Group 9, Group 10, Group 11, and Group 12 of the Periodic Table, or two or more transition metals selected from V, Cr, Mn, Fe, Co, Mo, W, and Re. In the case of two or more doping elements, the ranges specified above for z can correspond to a combined molar content of the two or more doping elements.
It is also contemplated that a Pt-based electrocatalyst can include another secondary metal, generally termed M′, in place of, or in combination with, Ni. Other suitable secondary metals include transition metals selected from Group 3, Group 4, Group 5, Group 6, Group 7, Group 8, Group 9, Group 10, Group 11, and Group 12 of the Periodic Table. In the case of two or more secondary metals, the ranges specified above for y can correspond to a combined molar content of the two or more secondary metals. Also contemplated are combinations of Pt with one or more secondary metals in a manner other than, or in conjunction with, alloying, such as heterostructures which include a first phase and a second phase, where the phases are joined together or next to one another, and the first phase and the second phase have different chemical compositions.
In some embodiments, a Pt-based electrocatalyst includes multiple nanostructures having the above-noted chemical composition, where any one or any combination of two or more of the following applies: (1) the nanostructures have sizes (or have an average size) in a range of up to about 100 nm, up to about 50 nm, up to about 40 nm, up to about 30 nm, up to about 20 nm, up to about 15 nm, up to about 10 nm, up to about 9 nm, up to about 8 nm, up to about 7 nm, up to about 6 nm, up to about 5 nm, or up to about 4.5 nm, and down to about 4 nm, down to about 3.5 nm, or less; (2) the nanostructures have at least one dimension or extent (or have at least one average dimension or extent) in a range of up to about 100 nm, up to about 50 nm, up to about 40 nm, up to about 30 nm, up to about 20 nm, up to about 15 nm, up to about 10 nm, up to about 9 nm, up to about 8 nm, up to about 7 nm, up to about 6 nm, up to about 5 nm, or up to about 4.5 nm, and down to about 4 nm, down to about 3.5 nm, or less; (3) the nanostructures have aspect ratios (or have an average aspect ratio) in a range of up to about 3, such as about 1 to about 3, about 1 to about 2.5, about 1 to about 2, or about 1 to about 1.5, or in a range of greater than about 3, such as about 4 or greater, about 5 or greater, or about 10 or greater; and (4) the nanostructures are largely or substantially crystalline, such as with a percentage of crystallinity (by volume or weight) of at least about 50%, at least about 60%, at least about 70%, at least about 80%, at least about 90%, at least about 95%, at least about 98%, or at least about 99% or more. Nanostructures of a Pt-based electrocatalyst can have a variety of morphologies, such as in the form of octrahedra having exposed (111) facets, although other morphologies are encompassed by this disclosure, including nanoparticles, nanorods, nanowires, or other elongated nanostructures having aspect ratios greater than about 3, as well as core-shell nanostructures, core-multi-shell nanostructures, and nanoparticle-decorated cores, among others.
In some embodiments, a Pt-based electrocatalyst includes multiple nanostructures that are surface doped with M, such that at least a majority (e.g., by weight, moles, or number) of M atoms are located within a depth of 5 atomic layers from an exterior of a nanostructure, such as within 4 atomic layers, within 3 atomic layers, or within 2 atomic layers. In some embodiments, at least a majority (e.g., by weight, moles, or number) of M atoms are in an oxidized state.
In some embodiments, a Pt-based electrocatalyst includes multiple nanostructures that are loaded on, dispersed in, affixed to, anchored to, or otherwise connected to a catalyst support, such as carbon black. In place of, or in combination with, carbon black, another catalyst support having suitable electrical conductivity can be used, such as another carbon-based support in the form of carbon fiber paper or carbon cloth, as well as metallic foams, among others. A combination of a Pt-based electrocatalyst loaded on a catalyst support can be referred to as an electrode material.
In some embodiments, an electrochemically active surface area (EASA) of a Pt-based electrocatalyst loaded on a catalyst support can be at least about 55 m2/gPt, at least about 60 m2/gPt, at least about 63 m2/gPt, at least about 65 m2/gPt, or at least about 67 m2/gPt, and up to about 70 m2/gPt, up to about 75 m2/gPt, or up to about 78 m2/gPt, or more, based on hydrogen under-potential deposition (Hupd). In some embodiments, a specific activity of the Pt-based electrocatalyst can be at least about 0.5 mA/cm2, at least about 1 mA/cm2, at least about 2 mA/cm2, at least about 3 mA/cm2, at least about 4 mA/cm2, at least about 5 mA/cm2, at least about 6 mA/cm2, at least about 7 mA/cm2, at least about 8 mA/cm2, at least about 9 mA/cm2, or at least about 10 mA/cm2, and up to about 10.3 mA/cm2, up to about 10.5 mA/cm2or more, at about 0.9 V versus a reversible hydrogen electrode (RHE) and based on Hupd. In some embodiments, a mass activity of the Pt-based electrocatalyst can be at least about 0.2 A/mgPt, at least about 0.5 A/mgPt, at least about 1 A/mgPt, at least about 1.5 A/mgPt, at least about 2 A/mgPt, at least about 2.5 A/mgPt, at least about 3 A/mgPt, at least about 3.5 A/mgPt, at least about 4 A/mgPt, at least about 4.5 A/mgPt, at least about 5 A/mgPt, at least about 5.5 A/mgPt, at least about 6 A/mgPt, or at least about 6.5 A/mgPt, and up to about 7 A/mgPt, up to about 7.5 A/mgPt, or more, at about 0.9 V versus a RHE and based on Hupd. In some embodiments, at least about 75% of an initial specific or mass activity can be retained after 4000 potential cycles, such as at least about 80%, at least about 85%, at least about 90%, or at least about 95%, and up to about 97%, up to about 98%, or more, and at least about 70% of the initial specific or mass activity can be retained after 8000 potential cycles, such as at least about 75%, at least about 80%, at least about 85%, or at least about 90%, and up to about 95%, up to about 97%, or more, when cycled between about 0.6 V and about 1.1 V versus a RHE at a scan rate of about 50 mV/s.
In some embodiments, a Pt-based electrocatalyst can be formed according to a manufacturing method including: (1) providing a dispersion of Pt—Ni nanostructures affixed to a catalyst support in a liquid medium; and (2) reacting a M-containing precursor, a Pt-containing precursor, and a Ni-containing precursor in the liquid medium to form M-doped Pt—Ni nanostructures.
In some embodiments, providing the dispersion in (1) includes reacting a Pt-containing precursor (which can be the same as or different from the Pt-containing precursor used in (2)) and a Ni-containing precursor (which can be the same as or different from the Ni-containing precursor used in (2)) in the presence of the catalyst support in the liquid medium to form the dispersion of Pt—Ni nanostructures affixed to the catalyst support. Suitable Pt-containing precursors (used in (1) and (2)) include an organometallic coordination complex of Pt with an organic anion, such as acetylacetonate, and suitable Ni-containing precursors (used in (1) and (2)) include an organometallic coordination complex of Ni with an organic anion, such as acetylacetonate. The liquid medium includes one or more solvents, such as one or more organic solvents selected from polar aprotic solvents, polar protic solvents, and non-polar solvents. In some embodiments, a solvent included in the liquid medium also can serve as a reducing agent for reduction of Pt and Ni, although the inclusion of a separate reducing agent is also contemplated. In some embodiments, a structure-directing agent, such as benzoic acid or other aromatic carboxylic acid, is also included in the liquid medium to promote a desired morphology of Pt—Ni nanostructures. Multiple metal-containing precursors including secondary metals different from Pt can be used, such as to form alloys of Pt and two or more secondary metals. Reaction can be carried out under agitation and under conditions of a temperature in a range of about 100° C. to about 300° C. or about 100° C. to about 250° C., and a time duration in a range of about 2 hours to about 24 hours or about 6 hours to about 18 hours.
In some embodiments, reacting in (2) includes adding or otherwise incorporating the M-containing precursor, the Pt-containing precursor (which can be the same as or different from the Pt-containing precursor used in (1)) and the Ni-containing precursor (which can be the same as or different from the Ni-containing precursor used in (1)) to the liquid medium. Suitable M-containing precursors include an organometallic coordination complex of M with an organic anion, such as carbonyl. Multiple different dopant-containing precursors can be used, such as to form nanostructures doped with two or more doping elements. Reaction can be carried out under agitation and under conditions of a temperature in a range of about 100° C. to about 300° C. or about 100° C. to about 250° C., and a time duration in a range of about 12 hours to about 60 hours or about 24 hours to about 60 hours.
The following example describes specific aspects of some embodiments of this disclosure to illustrate and provide a description for those of ordinary skill in the art. The example should not be construed as limiting this disclosure, as the example merely provides specific methodology useful in understanding and practicing some embodiments of this disclosure.
Because surface and near-surface features of a catalyst can have a strong influence on its catalytic performance, this example sets forth a surface engineering strategy to further enhance the performance of Pt3Ni(111) nanocatalysts. Efforts focused on Pt3Ni-based nanocatalysts because the bulk extended Pt3Ni(111) surface is one of the most efficient catalytic surfaces for oxygen reduction reaction (ORR). On the basis of the control over dopant incorporation of various transition metals onto the surface of dispersive and octahedral Pt3Ni/C (termed as M-Pt3Ni/C, where M=V, Cr, Mn, Fe, Co, Mo, W, or Re), ORR catalysts are developed which exhibit both high activity and stability. In particular, Mo—Pt3Ni/C catalyst of this example has high specific activity (about 10.3 mA/cm2), high mass activity (about 6.98 A/mgPt), and substantially improved stability for about 8000 potential cycles.
Highly dispersed Pt3Ni octahedra on carbon black was prepared by an efficient one-pot approach without using any bulky capping agents, which used platinum(II) acetylacetonate [Pt(acac)2] and nickel(II) acetylacetonate [Ni(acac)2] as metal precursors, carbon black as support, N,N-dimethylformamide (DMF) as solvent and reducing agent, and benzoic acid as a structure-directing agent (
For Pt3Ni, powder x-ray diffraction (PXRD) patterns of the colloidal products displayed peaks that could be indexed as those of fcc Pt3Ni (
To assess ORR catalytic activity, cyclic voltammetry (CV) was used to evaluate the electrochemically active surface areas (ECSAs). The catalysts were loaded (with the same Pt mass loading) onto glassy carbon electrodes. A commercial Pt/C catalyst [20 weight percent (wt. %) Pt on carbon black; Pt particle size, about 2 to about 5 nm] obtained from Alfa-Aesar was used as a baseline catalyst for comparison (
The ORR polarization curves for the different catalysts, which were normalized by the area of the glassy carbon area (about 0.196 cm2), are shown in
Because Mo—Pt3Ni/C exhibited an exceptional activity toward ORR, further examination was made of the doping effects for Pt3Ni/C modified by other transition metals. Pt3Ni/C catalysts doped with seven other transition metals—V, Cr, Mn, Fe, Co, W, or Re—were synthesized in a similar fashion with metal carbonyls (
Further evaluation was made of the electrochemical durability of the Mo—Pt3Ni/C catalyst using the accelerated durability test (ADT) between about 0.6 and about 1.1 V (versus RHE, 4000 and 800 cycles) in O2-saturated about 0.1 M HClO4 at a scan rate of about 50 mV/s. The Pt3Ni/C catalyst was used as a baseline catalyst for comparison. After 4000 and 8000 potential cycles, the Mo—Pt3Ni/C catalyst largely retained its ECSA and activity (
To investigate the cause of the enhanced durability of the Mo—Pt3Ni/C catalysts, cluster expansions of Pt—Ni—Mo NCs were used in Monte Carlo simulations to identify low-energy NC and (111) surface structures for computational analysis (details of calculations are provided in the Materials and Methods section). In vacuum, the equilibrium structures predicted by the cluster expansion have a Pt skin, with Mo atoms preferring sites in the second atomic layer along the edges connecting two different (111) facets (
The calculations indicate that doping NCs with Mo directly stabilizes both Ni and Pt atoms against dissolution and may inhibit diffusion through the formation of relatively strong Mo—Pt and Mo—Ni bonds. Calculations on a representative nanoparticle with dimensions and composition comparable with those observed experimentally (
Although the exact mechanisms by which the surface-doped Pt3Ni shows exceptional catalytic performance can be further evaluated, local changes in oxygen binding energies provide a possible explanation for at least some of the observed increase in specific activity, A Sabatier volcano of ORR catalysts predicts that ORR activity will be maximized when the oxygen binding energy is about 0.2 eV less than the binding energy on Pt(111). Calculations indicate that sites near the particle edge bind oxygenated species too strongly, such as in Pt(111), and sites near the facets of the particles bind oxygenated species too weakly, such as in Pt3Ni(111) (
Materials and Methods
Chemicals:
Platinum(II) acetylacetonate (Pt(acac)2, about 97%), nickel(II) acetylacetonate (Ni(acac)2, about 95%), cyclopentadienylvanadium(0) carbonyl (C5H5V(CO)4, about 98%), chromium(0) hexacarbonyl (Cr(CO)6, about 98%), dimanganese(0) decacarbonyl (Mn2(CO)10, about 98%), iron(0) pentacarbonyl (Fe(CO)5, >about 99.99%), dicobalt(0) octacarbonyl (Co2(CO)8, >about 99.99%), molybdenum(0) hexacarbonyl (Mo(CO)6, about 98%), tungsten(0) hexacarbonyl (W(CO)6, about 97%), dirhenium(0) decacarbonyl (Re2(CO)10, about 98%), benzoic acid (C6H5COOH, ≥about 99.5%), and N,N-dimethylformamide (DMF, ≥about 99.9%) were all purchased from Sigma-Aldrich. All chemicals were used as received without further purification. The water (about 18 MΩ/cm) used in all experiments was prepared by passing through an ultra-pure purification system (Aqua Solutions).
Preparation of Octahedral Pt3Ni/C Catalyst:
In a typical preparation of octahedral Pt3Ni/C catalyst, platinum(II) acetylacetonate (Pt(acac)2, about 8.0 mg), nickel(II) acetylacetonate (Ni(acac)2, about 4.0 mg), benzoic acid (C6H5COOH, about 61 mg) and about 10 mL commercial carbon black dispersed in DMF (about 2 mg/mL, Vulcan XC72R carbon) were added into a vial (volume: about 30 mL). After the vial had been capped, the mixture was ultrasonicated for about 5 minutes. The resulting homogeneous mixture was then heated at about 160° C. for about 12 h in an oil bath, before it was cooled to room temperature. The resulting colloidal products were collected by centrifugation and washed three times with an ethanol/acetone mixture.
Preparation of Transition Metal-Doped Pt3Ni/C Catalyst (M-Pt3Ni/C Catalyst):
The Mo—Pt3Ni/C catalyst was obtained by further growth of Mo on the preformed octahedral PtiNi/C catalyst. In a typical preparation of octahedral M-Pt3Ni/C catalyst, platinum(II) acetylacetonate (Pt(acac)2, about 2.0 mg), nickel(II) acetylacetonate (Ni(acac)2, about 1.0 mg) and molybdenum hexacarbonyl (Mo(CO)6, about 0.4 mg) were added to the suspension of unpurified Pt3Ni/C catalyst prepared above. After the vial had been capped, the mixture was ultrasonicated for about 30 minutes. The resulting mixture was then heated at about 170° C. for about 48 h in an oil bath, before it was cooled to room temperature. The resulting colloidal products were collected by centrifugation and washed three times with an ethanol/acetone mixture. The surface doping approach was robust and readily extended to other metals [V, Cr, Mn, Fe, Co, W, or Re] by replacing Mo(CO)6 with other transition metal carbonyl compounds in the above described process. In these other transition metal-doped Pt3Ni/C catalysts, highly dispersive NCs with octahedral morphology anchored on carbon black were obtained (
Characterizations:
Transmission electron microscopy (TEM) images were obtained on a FEI CM120 transmission electron microscope operated at about 120 kV, High resolution TEM images (HRTEM) and the high-angle annular dark-field scanning TEM (HAADF-STEM) and energy-dispersive X-ray spectroscopy (EDS) results were obtained on a FEI TITAN transmission electron microscope operated at about 300 kV. The STEM electron energy loss spectroscopy (EELS) tests were performed on an aberration corrected transmission electron microscope. The samples were prepared by dropping ethanol dispersion of catalysts onto carbon-coated copper TEM grids (Ted Pella, Redding, Calif.) using pipettes and dried under ambient condition. Powder x-ray diffraction (PXRD) patterns were collected on a Panalytical X'Pert Pro X-ray Powder Diffractometer with Cu—Kα radiation. X-ray photoelectron spectroscopy (XPS) tests were performed with Kratos AXIS Ultra. DLD spectrometer. The concentration of catalysts was determined by inductively coupled plasma atomic emission spectroscopy (TJA RADIAL IRIS 1000 ICP-AES).
Electrochemical Measurements:
A three-electrode cell was used to perform the electrochemical measurements. The working electrode was a glassy-carbon rotating disk electrode (RDE) (diameter: about 5 mm, area: about 0.196 cm2) from Pine Instruments. Ag/AgCl (about 3 M CF) was used as reference electrode. Pt wire was used as counter electrode. The Pt loading of octahedral M-Pt3Ni/C and octahedral Pt3Ni/C were about 0.8 μg (about 4.08 μgPt/cm2 based on the geometric electrode area of about 0.196 cm2). The lowest mass loading of the catalyst is about 0.80 μg. The electrochemical active surface area (ECSA) measurements were determined by integrating the hydrogen adsorption charge on the cyclic voltammetry (CV) scans at room temperature in nitrogen saturated about 0.1 M HClO4 solution. The potential scan rate was about 100 mV/s for the CV measurement. Oxygen reduction reaction (ORR) measurements were conducted in oxygen saturated about 0.1 M HClO4 solution which was purged with oxygen during the measurement. The scan rate for ORR measurement was about 10 mV/s. The ORR polarization curves were collected at about 1600 rpm. The accelerated durability tests (ADTs) were performed at room temperature in oxygen saturated about 0.1 M HClO4 solutions by applying cyclic potential sweeps between about 0.6 and about 1.1 V versus reversible hydrogen electrode (RHE) at a sweep rate of about 50 mV/s for 8000 cycles. For comparison, commercial Pt/C catalyst (Alfa Aesar, about 20 wt. % Pt, Pt particle size: about 2-5 nm) was used as the baseline catalyst, the same procedure as described above was used to conduct the electrochemical measurement, and the Pt loading was about 4.08 μgPt/cm2 for the commercial Pt/C catalyst.
Density Functional Theory (DFT) Calculations:
DFT calculations were performed using the Vienna Ab-initio Software Package (VASP) with the revised. Perdew-Burke-Eznerhof (RPBE) exchange-correlation functional. All DFT calculations were run with spin-polarization activated. The Mo_pv, Ni, Pt_pv_GW, O_GW, and H_GW PBE projector-augmented wave (PAW) potentials provided with VASP were used, and VASP was run with high precision. A single k-point at the center of the Brillouin zone was used for each nanoparticle. For bulk materials, a 16×16×16 k-point grid was used for a fcc unit cell, and the k-point grid was scaled appropriately for larger cells. Second-order Methfessel-Paxton smearing with a width of 0.2 eV was used to set partial occupancies. Real-space projectors were used to evaluate the non-local part of the PAW potential. Calculations were stopped when the difference for the total energy in successive ionic relaxation steps was less than 1 meV. The surface d-band centers were calculated using the site-projected densities of states for surface atoms.
Cluster Expansion:
Cluster expansions are parameterized models that can be used to rapidly and accurately predict the energies of different arrangements of atoms and vacancies on a lattice of sites. Here, a single cluster expansion is used to predict the energies of nanoparticles as a function of shape, size, and internal atomic order. A quaternary cluster expansion was generated on an fcc lattice in which each site could be occupied by molybdenum, nickel, platinum, or a vacancy. Site variable values of 0, 1, 2, and 3 respectively were assigned to these species. A discrete cosine basis was used to generate the cluster functions, where the bth basis function of the site variables is given by
To create the initial 136 structures used for the training data, a “dummy” cluster expansion was generated, composed of just nearest-neighbor pair clusters, with effective cluster interactions (ECIs) chosen in a way that assigned a value of −1 eV to atom-atom interactions (regardless of the species involved) and no energy to other interactions. These cluster expansions were used in Monte Carlo simulations at 2000 4500 K to generate random snapshots of nanoparticles. Two different sets of random nanoparticles were created. The first set of nanoparticles contained just Ni and Pt, where the numbers of Pt and Ni atoms were independently and randomly selected from a uniform distribution over all integers from 0 to 100. The second set of nanoparticles contained Mo, Ni, and Pt, where the numbers of Mo, Ni, and Pt atoms were independently and randomly selected from uniform distributions over integers from 0 to 10, 0 to 50, and 0 to 150 respectively. All nanoparticles were generated under the constraint that there had to be more than 85 total atoms in the nanoparticle, as the inclusion of smaller particles was found to lead to cluster expansions with poor predictive accuracy for multi-nanometer nanoparticles (potentially due to quantum size effects). Nanoparticles that experienced significant reconstruction upon relaxation, specified as an atom traveling more than 75% of the nearest-neighbor distance from its initial site, were excluded. These particles accounted for about 20% of the random structures generated. All nanoparticles were contained in a cubic cell with a lattice parameter of 28.8 Å. The resulting set of random nanoparticles included 74 Ni—Pt nanoparticles and 62 Mo—Ni—Pt nanoparticles. In addition to these structures, the training data was composed of the pure elements Mo, Ni, and Pt in a bulk fcc crystal, vacuum (a lattice containing just vacant sites), and various low-energy structures predicted over the course of evaluation of this example, for a total of 195 unique structures, To reduce the prediction error of the cluster expansion, the pure elements and vacuum were included twice in the training set, and the ECIs were fit to the DFT-calculated formation energies of fully relaxed nanoparticles relative to these reference states.
The cluster expansion included the empty cluster, the one-body (point) cluster, all 2-body clusters up to the 10th-nearest neighbor, all 3-body clusters up to the third-nearest neighbor, and all 4-, 5-, and 6-body clusters up to the second-nearest neighbor, for a total of 374 symmetrically distinct cluster functions. The ECIs for these cluster functions were fit to the training data using the Bayesian approach with a multivariate Gaussian prior distribution. The inverse of the covariance matrix for the prior, Λ, was diagonal, with elements given by
where nα is the number of sites in cluster function α, rα is the maximum distance between sites, and the parameters λ1, λ2, λ3, and λ4 were determined by using a conjugate gradient methodology to minimize the leave-one-out cross validation score, an estimate of prediction error. The final values for these parameters were 10−8, 1.102×10−12, 6.103, and 4.312 respectively. The resulting cluster expansion had a root mean square leave-one-out cross validation error of 0.742 meV per site, corresponding to 3.87 meV per atom.
Sample Structures:
The cluster expansions were used in Monte Carlo simulations to calculate thermodynamic averages, identify ground state structures, and identify sample structures. The structures referenced in the text of the example include a nanoparticle with composition Mo6Ni41Pt178, a 9-layer (111) slab with composition Mo2Ni7Pt27, and 4573-atom nanoparticles with compositions Ni1175Pt3398 and Mo73Ni1143Pt3357.
The nanoparticle with composition Mo6Ni41Pt178 (
The 9-layer (111) slab with composition Mo2Ni7Pt27 (
In
The 4573-atom nanoparticles were created by generating octahedra with the six vertex atoms removed. The length of the remaining edges is estimated to be about 4.1 nm, consistent with nanoparticles observed experimentally. The shapes of these nanoparticles were held fixed, and just the internal atomic order was allowed to vary. The compositions of the nanoparticles were set to match the Ni0.257Pt0.743 and Mo0.016Ni0.25Pt0.734 compositions observed experimentally. The average site occupancies of these particles at 170° C. are shown in
The most favorable site for Mo surface segregation in the presence of an adsorbed oxygen atom for the Mo6Ni41Pt178 particle was determined by evaluating Mo segregation to each of the nearest face, vertex, and edge sites for each of the symmetrically distinct Mo atoms. In each case, the adsorbed oxygen atom was placed atop the surface Mo atom, as calculations indicate that this is the most favorable site for oxygen adsorption.
Surface Segregation:
To assess the energetics of surface segregation, DFT calculations were performed on both the extended (111) slab and the Mo6Ni41Pt178 particle. For the clean slab in vacuum, Mo is more stable at a subsurface site than the lowest-energy surface site by 0.881 eV per Mo atom. For the nanoparticle in vacuum, the subsurface site is favored over the lowest-energy neighboring surface site by 1.110 eV. The situation reverses in the presence of oxygen. In the presence of adsorbed oxygen on the (111) surface (with ¼ monolayer coverage), there is a driving force of 1.559 eV per Mo atom for Mo to segregate to the surface, and the oxygen preferentially adsorbs atop the surface Mo atom. For the Mo6Ni41Pt178 nanoparticle, similar results were found: in the presence of an adsorbed oxygen atom Mo preferentially segregates to a vertex site, and the driving force for this segregation is 1.533 eV. The Mo6Ni41Pt178 nanoparticle with a single Mo atom segregated to the energetically-preferred vertex site was used to assess the stability of surface Mo-oxide species against reduction to H2O. The structures used in these calculations, composed of one, two, and three oxygen atoms adsorbed on the vertex Mo atom, are shown in
Due to the relatively strong Mo—Pt and Ni—Pt nearest-neighbor bonds, both Mo and Ni prefer to occupy similar sites with many Pt nearest neighbors. However in oxidizing conditions, the energetic driving force for Mo segregation to the surface is much stronger than the driving force for Ni segregation. For the Mo2Ni7Pt27 slab with ¼ monolayer oxygen coverage, the calculated driving force for 2nd-layer Mo to migrate to the surface is 1.559 eV per atom, as opposed to 0.284 eV per atom for Ni. This indicates that in oxidizing conditions Mo atoms may “crowd out” Ni atoms on the particle surface, reducing the number of surface Ni atoms available for dissolution.
Stability Enhancements:
The cluster expansion was used to evaluate the effects of substituting a single Mo atom into all the sites in the representative Ni1175N3398 particle (
If the Mo atom is on an edge or vertex site, the energy involved to remove a Pt (Ni) atom from a neighboring edge or vertex site increases by an average of about 362 (about 201) meV, with values ranging from about 346 (about 160) to about 444 (about 214) meV. These values are supported by DPT calculations on the Mo6Ni41Pt178 nanoparticle which predict that the presence of a Mo atom on a vertex site stabilizes the Pt atom on the neighboring vertex site by about 458 (about 444) meV with (without) an oxygen atom adsorbed atop the Mo atom. The smallest increase in the energy involved to remove Pt from an edge or vertex site (about 346 meV) was observed to occur with the Mo atom at a non-vertex edge site, with the Pt atom on a neighboring non-vertex edge site, The greatest increase in the energy involved to remove Pt from an edge or vertex site (about 444 meV) was observed to occur with the Mo atom at a vertex site and the Pt atom on a neighboring vertex site. The smallest increase in the energy involved to remove Ni from an edge or vertex site (about 160 meV) was observed to occur with the Mo atom at a vertex site, with the Ni atom on a neighboring non-vertex edge site. The greatest increase in the energy involved to remove Ni from an edge or vertex site (about 214 meV) was observed to occur with the Mo atom at a non-vertex edge site and the Ni atom on a neighboring non-vertex edge site.
The effects of second- and third-nearest-neighbor interactions were also investigated, If the Mo atom is on an edge or vertex site, the energy involved to remove a Pt (Ni) atom from a second-nearest-neighbor site on an edge or vertex increases by an average of about 3 (about −15) meV, with values ranging from 0 (about −21) to about 21 (about −6) meV. If the Mo atom is on an edge or vertex site, the energy involved to remove a Pt (Ni) atom from a third nearest-neighbor site on an edge or vertex increases by an average of about −16 (about −3) meV, with values ranging from about −42 (about −26) to about 17 (about 39) meV.
As used herein, the singular terms “a,” “an,” and “the” include plural referents unless the context clearly dictates otherwise. Thus, for example, reference to an object can include multiple Objects unless the context clearly dictates otherwise.
As used herein, the term “set” refers to a collection of one or more objects. Thus, for example, a set of objects can include a single object or multiple objects.
As used herein, the terms “substantially” and “about” are used to describe and account for small variations. When used in conjunction with an event or circumstance, the terms can refer to instances in which the event or circumstance occurs precisely as well as instances in which the event or circumstance occurs to a close approximation. For example, when used in conjunction with a numerical value, the terms can refer to a range of variation of less than or equal to ±10% of that numerical value, such as less than or equal to ±5%, less than or equal to ±4%, less than or equal to ±3%, less than or equal to ±2%, less than or equal to +1%, less than or equal to +0.5%, less than or equal to ±0.1%, or less than or equal to ±0.05%.
As used herein, the terms “connect,” “connected,” and “connection” refer to an operational coupling or linking. Connected objects can be directly coupled to one another or can be indirectly coupled to one another, such as via another set of objects.
As used herein, the term “size” refers to a characteristic dimension of an object. Thus, for example, a size of an object that is spherical can refer to a diameter of the object. In the case of an object that is non-spherical, a size of the non-spherical object can refer to a diameter of a corresponding spherical object, where the corresponding spherical object exhibits or has a particular set of derivable or measurable properties that are substantially the same as those of the non-spherical object. When referring to a set of objects as having a particular size, it is contemplated that the objects can have a distribution of sizes around the particular size. Thus, as used herein, a size of a set of objects can refer to a typical size of a distribution of sizes, such as an average size, a median size, or a peak size.
While the disclosure has been described with reference to the specific embodiments thereof, it should be understood by those skilled in the art that various changes may be made and equivalents may be substituted without departing from the true spirit and scope of the disclosure as defined by the appended claims. In addition, many modifications may be made to adapt a particular situation, material, composition of matter, method, operation or operations, to the objective, spirit and scope of the disclosure. All such modifications are intended to be within the scope of the claims appended hereto. In particular, while certain methods may have been described with reference to particular operations performed in a particular order, it will be understood that these operations may be combined, sub-divided, or re-ordered to form an equivalent method without departing from the teachings of the disclosure. Accordingly, unless specifically indicated herein, the order and grouping of the operations is not a limitation of the disclosure.
This application claims the benefit of U.S. Provisional Application Ser. No. 62/146,803, filed on Apr. 13, 2015, the content of which is incorporated herein by reference in its entirety.
This invention was made with Government support under W911NF-09-1-0433, awarded by the U.S. Army, Army Research Office. The Government has certain rights in the invention.
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/US2016/027295 | 4/13/2016 | WO | 00 |
Number | Date | Country | |
---|---|---|---|
62146803 | Apr 2015 | US |