Ram
Subbaraman
a and
Subramanian K. R. S.
Sankaranarayanan
*b
aNuclear Engineering Division, Argonne National Laboratory, Argonne, IL 6043, USA
bCenter for Nanoscale Materials, Argonne National Laboratory, Argonne, IL 60439, USA. E-mail: skrssank@anl.gov
First published on 17th July 2012
Oxidative stabilities of single component nanoclusters are studied in the presence of an oxophilic element. Atomistic simulations are used to study the early stages of heterogeneous oxidation of a representative binary alloy nanoclusters such as Ni–Al. The role of oxophilic Al atoms in enhancing the oxidation of Ni atoms in the Ni–Al alloy is highlighted. The simulated trajectories are used to better elucidate the mechanism of oxidation for two elements with distinct oxidation rates. Increased oxidation kinetics for particles with higher Al content is consistent with higher oxophilicity of Al atoms in comparison with Ni atoms. Ni rich alloys are found to be least oxidized with significant phase segregation between the Ni and Al phases. Furthermore, the oxides formed were found to be non-stoichiometric for Ni rich alloys and near stoichiometric for Al rich alloys. Phase segregation for the alloys is also found to be induced, resulting from the relative differences in the atom mobilities during the oxidation process. Thermodynamic and kinetic considerations are used to understand the higher degree of oxidation of Ni atoms in the presence of a higher composition of Al oxides. The simulations suggest a general strategy that allows for tuning the room temperature oxidative stability of single component nanoclusters by alloying with a higher oxophilic element.
Ni–Al alloys, also known as aluminides, are regarded as promising candidates for the next generation of high temperature and high-performance structural materials because of their high melting points, relatively low densities, good strength and high temperature corrosion and oxidation resistance.3 As important intermetallics, they showed great potential applications in automobile engines, aircraft, and electricity generation and energy conversion equipment.4 The corrosion properties of such materials are crucial in these applications, and as a result understanding the oxidation of these nanomaterials is critical. Furthermore, nanoscale oxides formed on Ni–Al alloys are regarded as promising candidates for the next generation of high-temperature and high-performance structural materials because of their high-melting points, relatively low densities, good strength and high-temperature corrosion and oxidation resistance.5,6 Additionally, NiO/Al2O3 mixed oxides are used in several interesting magnetic and catalytic applications.7,8 These applications require precise control of the oxide structure, composition and quality. Besides being able to control the microstructure of these oxides, controlling the stoichiometry and composition of these mixed oxides is of utmost importance.
Although the oxidation of metal nanoparticle is of significant practical value, very little is known about them.9 Tuning the oxidation characteristics and stability of the growing oxide layer, especially in low dimensional systems, remains a major challenge. Given the substantial contribution from the surface energy term (work function) to the overall Gibbs energy, it is expected that the oxidation of nanoparticles will behave significantly differently from their bulk counterparts.10 This is in line with other experimental and computational efforts.11,12 Further complications arise when studying alloy nanoparticles, resulting from the significant compositional variations during the oxidation. For example, the nanoscale oxidation of Ni–Al alloys can result in a complex structure consisting of a mixture of Al2O3, NiO, and NiAl2O4 oxides.8 The structure and composition of this oxide layer depends on the alloy composition, oxygen pressure, and oxidation temperature.13 Gaining a theoretical insight into the low temperature formation and growth mechanism of mixed oxides has seldom been addressed. Specifically, processes related to controlling the composition of binary oxides are extremely difficult due to microstructural changes brought about by phenomena such as surface segregation and micro-mixing of dopants.14 Furthermore, differences in the oxidation kinetics of the alloy constituents can also result in the formation of either a singular oxide layer (for the element with higher oxophicility),15 or result in the formation of equimolar mixed oxides.16 Phase stability of the binary oxides is another critical issue limiting the compositional control to high temperatures.17,18
Given these complexities, developing a methodology to better treat/understand the mechanism of low temperature oxidation and oxide growth at the atomic scale over a range of alloy compositions will be vital for the progress of development strategies for synthesizing complex oxides. Besides treating the mechanism, such methodologies should be able to answer the questions regarding microstructural and compositional changes that occur during the early stages of alloy oxidation, which are particularly poorly understood thus far.13,19 While low-temperature oxidation of pure metal surfaces and thin films such as Al(100) and Ni(100) as well as nanomaterials of Ni (primarily experimentally) and Al have been thoroughly investigated using experimental and computational techniques,20–23 very little effort has been dedicated toward understanding the Ni–Al alloys, most definitely on the nanoscale level. Given the great practical importance of such materials,24 obtaining insights on a molecular level into the processes occurring during the oxide growth on the complex Ni–Al alloy nanoparticle systems is the main focus of the work presented here.
Here, we demonstrate, using a variable charge molecular dynamics simulation study, the ability to tune the atomic scale oxidative stability of a metal nanocluster by alloying it with a more oxophilic element. We specifically report on the early stages of oxidation and dynamics of oxide growth on Ni–Al alloy nanoclusters, with a focus on temperature as well as size effects. Formation of the oxide scale, and its role in stabilizing/destabilizing the surface segregation/mixing of the alloy materials were investigated. The dependence of the composition of the alloy on the overall oxidation stability of the alloys as well as the resulting morphologies is evaluated using dynamical correlation functions. The valence states of the oxide constituents are characterized using the simulated trajectories. Finally, the atomistically simulated evolution of charges and atomic diffusivities for various alloy compositions is used to clearly explain the observed differences in the oxide composition and the stoichiometry of the grown oxide layer.
(1) |
In the above equation, ρj is the electron density due to atom j. The first term is the embedding energy for atom i which is embedded into the background electron density, and the second term is the core–core pair repulsion between atoms i and j separated by a distance rij. The force field parameters used to model Al and Ni were obtained from the correlations developed by Neurock et al.30 which is the same as that used for modeling the metallic interactions in the charge transfer ionic potential used for subsequent oxidation simulations.
The standard Metropolis MC method was used to generate the minimum energy configuration at 300 K. Monte Carlo (MC) simulation methods based on the Metropolis algorithm have been successfully been used in previous surface segregation calculations for various alloys to generate equilibrium alloy configurations and provide an average thermodynamic profile over a thermodynamic equilibrium ensemble. The MC simulations were carried out in a canonical ensemble where the total number of atoms of each element and the temperature are held constant. In the MC approach, starting from an initially randomly distributed atomic configuration of the alloy, the successive configurations are generated in proportion to the probabilities of a configuration occurring in the equilibrium ensemble. During each of the MC step, one of the following two configurational transformations is tried out with an equal probability:
(1) We displace a randomly selected atom from its original position in a random direction. The magnitude of the displacement is chosen to be in the range of (0, rmax]. At a given temperature, the maximum displacement rmax is selected such that the acceptance rate of new configurations is about 0.5 during the equilibrated part of the simulations.
(2) Two randomly selected atoms with different element types are exchanged and the atoms are switched in the lattice.
Event (1) accounts for the relaxation and vibration processes whereas event (2) accounts for the inter-diffusion process in the model system. The transition probability is based on the Boltzmann distribution and is given by:31
(2) |
Here, kB is the Boltzmann constant and T is the temperature and ΔE is the change in system energy when an event is executed. If P ≥1 (decrease in energy), the new configuration is always retained, while if P < 1 (increase in energy), the new configuration is accepted with probability P. The procedure is repeated for multi-million MC steps and the convergence was obtained when the minimum energy obtained during the simulation steps becomes constant (ΔE < 0.01 eV), and the state of the system corresponding to that minimum energy was retained as shown in Fig. 1. A large number of simulations (100), with random initial configurations, each having (1000 times the total number of atoms) steps were performed to ensure convergence.
Fig. 1 Schematic showing the simulation cell of alloy nanoparticle and the vacuum surrounding it. The nanoparticle of radius R is at the center of the simulation cell. A spherical reflecting wall located at 2.5R holds the gas phase oxygen atoms. |
It is important to note that the alloying of two elements with a significant size mismatch (lattice constants for Ni and Al are 3.5 and 4.05 Å, respectively) leads to the formation of a defective nanoparticle. Besides the lattice mismatch, other factors such as surface energy differences and mixing entropies of the alloying elements dictate the final crystal structure and composition. All of these can lead to the crystallographic facets being not well-defined, which justifies the use of MC simulations to generate the initial configuration for MD simulations. Also, given the relatively small sizes of the particle, MC equilibration at 300 K can also influence the “crystallinity” of the facets.
Thus, the main purpose of the MC equilibration is to identify the energetically optimized structure, both in terms of surface segregation properties as well as relative bond distances. Given the large lattice mismatch between Ni and Al, and the large surface energy differences, randomly filling the nanoparticle co-ordinates with Ni and Al atoms, can lead to a non-equilibrium structure. Our purpose here is to understand the role of surface segregation and the objective of the MC is to mimic a simulated annealing process, where the atoms are placed in the most energetically favorable positions. In the absence of such equilibration, at relatively large compositions of Al, there might be significant phase segregation and/or very defective oxide structures. Also, using the energetically optimized structures helps us avoid any ambiguities related to variations in the starting configurations for various alloys.
Here, we study oxidation of Ni–Al alloy nanoparticle having 2.5 nm diameter over a range of alloy composition (15–75% Al). For each of the starting metal alloy compositions, the initial minimum energy configuration was generated using Monte-Carlo simulations and consisted of a surface segregated structure with Al atoms segregating and preferentially locating themselves at the exposed faces (Fig. 2).13 The potential model and MD simulation details are further discussed below.
Fig. 2 Initial configurations of alloy nanoparticles after 1M MC simulation steps. Increasing Al composition leads to Al (pink) segregation on the surface in lieu of Ni (blue). Room temperature oxidation kinetics of Ni–Al nanocluster having 2.5 nm diameters as a function of (a) alloy composition; (b) pressure for 50%Ni–Al. |
The oxidation of the alloy nanoparticles is initiated by introducing O2 molecules in the vacuum with their radial positions chosen randomly (Fig. 1). The velocities of the O2 are chosen from a Maxwell-Boltzmann distribution corresponding to the required temperature. Additionally, spherical reflecting boundary conditions are imposed to the molecules that might reach the simulation box limit. The gas pressure is maintained constant during the simulation by introducing a new O2 molecule only when the previous molecule dissociates and forms bonds with the metal atoms. The equations of motion are integrated using a leapfrog scheme with time steps of 1 fs. The charge relaxation procedure used to minimize the electrostatic energy subject to the electro-neutrality principle is very time consuming. Hence, the atomic charges were updated every 10th MD step. The influence of a more frequent charge update was found to have no influence on the observed simulation results.
We have also modeled the oxidation of pure Ni nanoparticles. Ni has much slower oxidation kinetics than Al and hence the final oxide structure in the case of pure Ni comprises of discontinuous surface oxide formation as shown in the ESI.† The oxide formation and oxidation kinetics are similar to those observed for low-Al content Ni–Al alloy.
To investigate the effect of gas pressures on the binary oxide growth, we simulated oxidation of 50% Ni–Al alloy particle with 2.5 nm diameter at three different pressures (4, 7 and 10 O2). The temperature is kept constant at 300 K. Pressure is controlled in the MD box by maintaining the gas phase number for oxygen molecules constant. Gas phase molecules were characterized by a lower negative charge (typically <−0.2), i.e. non-adsorbed oxygen atoms. The kinetic curves showing the effect of gas pressure on the oxide growth on 50% Ni–Al alloy particles are shown in Fig. 2(b). As expected, the oxidation kinetics is found to increase with increasing gas pressure. This is consistent with the enhanced diffusion as well as the number density of oxygen atoms available for the oxidation process. Experimental observations of oxidation of Al surfaces have been shown to follow such trends with partial pressure of oxygen.36–39 The effect of pressure is complemented by studying the effect of particle size for a given composition at constant temperature and pressure. This is shown in Fig. 3. The oxygen uptake is defined in terms of the ratio of the oxygen/metal atoms in the growing oxide particle. With an increase in particle size, the ratio of surface area to volume goes down. Hence, the surface area with respect to the total number of metal atoms is expected to go down with size, which explains the observed decrease in the oxophilic character. Therefore, the relative oxophilicity of the particles decreases with increasing particle size, and as a result, the largest particles considered here are found to exhibit a very low degree of oxidation at the same time in comparison with the smaller particles. This can be associated with the lower values of the metal atom/oxygen ratio. This clearly indicates that the oxidation kinetics for an infinite surface will be a lot slower in comparison with nanoparticles, further confirming that the rate of oxidation is a strong function of particle size, partial pressure of oxygen as well as the relative composition of oxophilic elements. Having established the dependence of oxidation kinetics on particle size, composition and partial pressure of oxygen, we then focus on the role of temperature on the activation energy and hence the kinetics of oxidation.
Fig. 3 Size effect on the oxidation kinetics of 50% Ni–Al nanocluster. The particles were maintained under 7 O2 molecules gas pressure. |
We have also evaluated the effect of temperature on the binary oxide growth. We simulated oxidation of 50% Ni–Al alloy particle with 2.5 nm diameter at three different temperatures (300, 400 and 500 K). The oxygen gas pressure was kept constant at 2 O2. For each of the simulation runs, the temperature was kept constant using a Nose-Hoover thermostat. The kinetic curves showing the effect of temperature on oxide growth on 50% Ni–Al alloy particles are shown in Fig. 4. The oxidation kinetics at any given temperature exhibits an initial stage of fast oxide growth followed by a low growth rate. As expected, the kinetics of oxide growth becomes faster at higher temperatures. This can be attributed to the increased mobility of the atoms of both the metal and oxygen, as a result of which the overall activation energy for the process is lowered. We refrain from any detailed kinetic analysis at this juncture and so we will not discuss the variation in activation energies as a function of temperature. At the end of 300 ps of simulation time, the oxygen uptake increased from θox = 0.5 at 300 K to θox = 0.75 at 500 K.
Fig. 4 Temperature dependence of the oxidation kinetics of 50% Ni–Al nanocluster having diameters of 2.5 nm. |
It is worth noting that the Ni–Al alloys have potential high-temperature applications (T > 1000 K). It is well known that the nanoparticles undergo melting point depression when compared to their bulk counterparts. Therefore, at these high temperatures, the nanoparticles are likely to undergo a significant change in their structure during their operation. The crystal structure is likely to be more amorphous and one would expect to see higher inter-mixing between the alloying components. The oxidation kinetics and oxide growth will be expected to be different from those reported at low temperatures (Fig. 4).
Fig. 5 Bond length defined as the separation distance of the two atoms of the oxygen molecule plotted versus the distance of the center of mass of the pair from the Ni–Al surface. (a) Simulated trajectories of four typical runs of oxygen molecules approaching the Ni–Al(100) surface are shown to illustrate the dissociative chemisorption of oxygen molecules. (b) Comparison with the trajectories of the oxygen molecule approaching the Ni–Al alloy nanoparticle is also shown. |
As a result, in the case of nanoparticle oxidation, our analysis of the simulation trajectories (Fig. 5(b)) suggests that the oxygen molecule undergoes a dissociative adsorption onto the nanoparticle surface. As seen in Fig. 5(b), the oxygen molecule dissociation starts effectively when the center of mass of the molecule is located at a position that ranges between 1.8 and 3 Å above the particle. Interestingly, we find that the nanoparticle oxidation is initiated at the low coordination sites such as vertex and edges as shown in Fig. 6. Closer examination of the MD simulation trajectories indicated that nanoparticle oxidation is favored on low co-ordinated sites such as the vertex of icosahedra. We also find that prior to adsorption on these low coordinated sites, the oxygen molecule undergoes rotation which acts as a precursor before the separation distance begins to significantly exceed the bond length. This is in contrast with surfaces, where at relatively low surface coverages the adsorption of oxygen molecule can proceed via a bridge bonding or a two-site process. In contrast, at the vertex and edges of the nanoparticles, the adsorption seems to proceed via an atop adsorption step followed by a relatively fast re-orientation and dissociation on the surface. This can be again attributed to the difference in surface energy as well as the relative inter-atomic strain that exists between the atoms in the nanoparticle.
Fig. 6 Snapshots showing initial oxidation of 2.5 nm diameter Ni–Al alloy nanoparticles at 300 K and 2O2 pressure. Configurations at (a) initial (b) 1 ps (c) 5 ps (d) 10 ps (e) 20 ps (f) 30 ps (g) 40 ps. Oxidation is initiated at the low coordinated sites such as corners and edges. The Al atoms are shown in green, Ni in blue, and oxygen in red. |
Fig. 7 Radial distribution functions (Al–O) and (Ni–O) in the Ni–Al–O system (4 nm diameter) taken at 200 ps of exposure time for room temperature oxidation of Ni–Al alloys over a range of compositions and 4O2 pressure. |
Similarly, the calculated Ni–O PDF for various alloy compositions is shown in Fig. 7. The experimentally reported bond lengths of Ni–O by Caputi et al. correspond to a value of around 1.94 Å.46 Hence, formation of a nickel oxide phase should result in a sharp peak in Ni–O PDF at that distance. At higher Ni compositions (>85% Ni), the peaks in the Ni–O RDF are not well-defined and show a distribution which indicates a lower degree of formation of nickel oxide at low temperatures. The appearance of such broad and shifted peak positions suggests the existence of both non-stoichiometric as well as amorphous oxides. This could possibly arise from oxide induced surface segregation of the Ni and Al atoms that can lead to alteration of the Ni–Ni as well as Ni–O distances within the nanocluster. Interestingly, as suggested by the peak intensities in Fig. 7, the extent of formation of Ni oxide also increases when the Al composition is increased from 15 to 30% and appears to saturate for Al composition higher than 30%. For Al rich alloy particles, the oxide scale is predominantly rich in AlxOy type oxides. This is in agreement with the experimental studies of low temperature oxygen adsorption on Ni–Al surfaces by Jaeger et al. and Franchy et al.40,47
To further investigate the chemical nature of the oxide, we analyzed the ionic charge distribution of oxygen ions across the radius of the Ni–Al nanoparticle for various alloy compositions at the end of 200 ps of simulation time (Fig. 8). We find that in the case of 15% Al alloy, the nanoparticle is only partially oxidized and the core of the particle (<7 Å), which is Ni rich in composition, remains in metallic form. The average charge values of the oxygen ions in the oxidized near surface region are still significantly below the valence value of −2.0, which suggests that the oxide skin is highly non-stoichiometric and it is predominantly comprised of weakly bound oxide anions. As the Al composition of the alloy is changed to 30%, the depth of oxidation is found to increase, resulting in shrinking of the metallic core (<2.5 Å). For higher Al compositions (≥50% Al), we observe that the entire particle is oxidized and the average charges on the oxygen ions tend to approach their valence values, which is indicative of the formation of a more stoichiometric oxide.
Fig. 8 Spatial distribution of the average charges on the oxygen atoms as a function of alloy composition in the Ni–Al–O system (4 nm diameter) at room temperature and 4O2 pressure. |
The reduced oxygen stoichiometry and oxidation characteristics of Ni rich alloys are also supported by thermodynamic considerations. Comparison of the standard enthalpy of formation of the respective oxides, Hf(NiO) = −489.5 kJ mol−1 and Hf(Al2O3) = −1675.7 kJ mol−1 suggests that it is more favorable to form the Al–O bond.48 This manifests itself in the form of reduced extent of nickel oxide formation shown in the pair distribution for Ni rich alloys calculated by averaging over many such configurations.
To further evaluate the differences in the oxide stoichiometry, we calculate the average valence state of metal ions as a function of alloy composition (% Al) at the end of the oxidation process (Fig. 9). For the Ni–Al alloy compositions lower than 50%, it can be seen that neither Al nor Ni atoms assume full valence values. The Ni atoms are weakly charged and the charge transfer takes place primarily between the O atom and the nearest Al atoms. Interestingly, for higher Al in the alloy composition (≥50% Al), we observe that a significant fraction of both Al and Ni atoms oxidizes and attains a charge value close to their valence values of “3+” and “2+”, respectively. Increasing the Al composition in the alloy is expected to enhance the oxygen uptake for these nanoparticles. It appears that the increased uptake of O atoms into the matrix, along with the presence of more oxophilic element such as Al, increases the oxidizability of Ni atoms. Clearly then, we have enhanced the oxidizability of Ni atoms in the presence of Al atoms. This could be explained either with the variation in the electronic properties of Ni in the presence of Al-oxides or with the local increase in oxygen partial pressure experienced by Ni atoms due to the enhanced oxygen uptake by the Al atoms. For example, owing to the faster oxidation kinetics, the oxygen species acting as electron acceptors for electrons from the metal will be more and hence generate stronger bonding chemisorption ions leading to Ni oxidation. We will refrain from discussing these details in this manuscript and instead will focus on the intrinsic phase segregation behavior of Ni and Al atoms during the oxidation process. Simply, our simulations suggest that the fraction of Ni oxidized in the Ni–Al alloy nanoparticles is strongly correlated with the intake of oxygen atoms. Thus, kinetic considerations play an important role in dictating the overall oxide composition for Al-rich alloy particles.
Fig. 9 Average valence state of metal ions as a function of alloy composition (% Al) at the end of the oxidation process for the Ni–Al–O system (4 nm diameter) at room temperature and 4O2 pressure. |
To further evaluate the compositional variation of the growing oxide layer on the binary alloy nanoparticle, we analyze the radial atomic density distributions of Al and Ni atoms in 50% Ni–Al alloy as a function of the exposure time (Fig. 10(a)–(b)). The initial Ni–Al metal alloy nanocluster is comprised of Al rich surface and a Ni rich core. However, as seen in initial density profiles at 1 ps (Fig. 10(a), (b)) and the snapshots in Fig. 10, the segregation is not complete. There are Al atoms that are located in the cluster core and Ni occupies a small fraction of the surface sites. As seen in Fig. 10(a), the outward movement of Al atoms characterizes the early stage oxidation and oxide growth on 50% Ni–Al alloy nanoparticle. Simultaneously, as seen in Fig. 10(b), there is an inward movement of surface Ni atoms to the cluster core. The opposite movement of the constituent cations upon oxidation of the bimetallic nanoparticle is attributed to the differences in the oxidation kinetics of Al and Ni. Such a diffusion phenomenon can lead to a phase separation of the oxidized species from the unoxidized one and might lead to a segregation profile that is more pronounced than that in the starting bimetallic alloy particle. As seen from the snapshots in Fig. 6, the oxidation of Ni–Al alloys at longer times indeed leads to phase separation of aluminum oxide from metallic Ni. Even in the case of Ni–Al(100) thin films, oxidation proceeds with the formation of surface oxide, comprised primarily of aluminum oxide, with an accompanying phase segregation of Ni atoms as seen in Fig. 11.
Fig. 10 Temporal evolution of oxide composition during the nanocluster oxidation. (a) Variation of Ni composition in the nanocluster during the oxide growth process. (b) Variation of Al composition in the nanocluster during the oxide growth process. The evolution of the density profile of Ni and Al indicates phase segregation of metallic Ni from the growing aluminum oxide. The snapshots shown provide a visual evidence of oxidation induced phase segregation phenomena in Ni–Al inter-metallic alloy. Aluminum is shown in pink, Ni in blue and oxygen in red. |
Fig. 11 Snapshots showing oxidation induced phase segregation in Ni–Al(100) thin films. Oxidation results in the formation of a surface oxide comprised primarily of alumina. Ni segregation to the oxide metal interface is also seen at longer simulation times (>100 ps). Aluminum is shown in green, Ni in blue and oxygen in red. |
The phase separation of the oxidized aluminum species from the unoxidized metallic nickel can be explained based on thermodynamic and kinetic considerations. As mentioned earlier, comparison of the standard enthalpy of formation of the respective oxides, Hf(NiO) = −489.5 kJ mol−1 and Hf(Al2O3) = −1675.7 kJ mol−1 suggests that it is more favorable to form Al–O bond.48 Also, ab initio studies of oxygen adsorption on Ni–Al surfaces also indicate the preference of O atoms to have Al as a neighbor rather than Ni.41 The presence of O atom on a Ni–Ni bridge site was found to push the Ni atoms away, while the neighboring Al atoms moving towards it, following along the lines of stronger oxophilicity of Al as compared to Ni atoms. Previous ab initio as well as experimental investigations suggest much lower Ni–O bond formation for low temperature oxide grown on Ni–Al alloys and do not rule out the presence of metallic Ni on the oxide scale.35 Our MD simulations also suggest that the diffusivity of Al and O is at least an order of magnitude higher than the Ni atoms that further contributes to the observed phase segregation. It is worth noting that the phase separation of mixed oxide grown on inter-metallic alloys is a rather complicated dynamical process and is typically involved with multi-phase transformation and competition among the various phases, which would require further investigation.
Footnote |
† Electronic Supplementary Information (ESI) available. See DOI: 10.1039/c2ra20464a/ |
This journal is © The Royal Society of Chemistry 2012 |