Resolving Precursor Deligation, Surface Species Evolution ... - Lei Lab

May 3, 2013 - species under practical ALD operation conditions in real time.38. It provides crucial .... base pressure at ∼3 Torr with 50 sccm ultra...
2 downloads 0 Views 2MB Size
Article pubs.acs.org/JPCC

Resolving Precursor Deligation, Surface Species Evolution, and Nanoparticle Nucleation during Palladium Atomic Layer Deposition Yu Lei,† Junling Lu,† Haiyan Zhao,‡ Bin Liu,§ Ke-Bin Low,∥ Tianpin Wu,‡ Joseph A. Libera,† Jeffrey P. Greeley,⊥ Peter J. Chupas,‡ Jeffrey T. Miller,# and Jeffrey W. Elam*,† †

Energy Systems Division, Argonne National Laboratory, Lemont, Illinois 60439, United States X-ray Science Division, Advanced Photon Source, Argonne National Laboratory, Lemont, Illinois 60439, United States § Center for Nanoscale Materials, Argonne National Laboratory, Lemont, Illinois 60439, United States ∥ Research Resources Center, University of Illinois at Chicago, Chicago, Illinois 60607, United States ⊥ School of Chemical Engineering, Purdue University, West Lafayette, Indiana 47907, United States # Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, United States ‡

S Supporting Information *

ABSTRACT: The synthesis of highly dispersed palladium nanoparticles on TiO2 surfaces from palladium hexafluoroacetylacetonate (Pd(hfac)2) was investigated using in situ infrared (IR) spectroscopy, in situ X-ray absorption spectroscopy (XAS), and in situ pair distribution function (PDF) measurements under practical atomic layer deposition conditions. Residual surface chlorine was found to directly participate in the transformation of organometallic compounds to nanoparticles. Deligation of the Pd(hfac)2, evolution of the surface species, and nucleation of the Pd nanoparticles were precisely resolved. This knowledge can help direct the future design of advanced heterogeneous catalysts from organometallic compounds.



INTRODUCTION Controlling the dispersion and structure of noble metal nanoparticles of uniform size has been critical to achieving efficient catalysts.1−8 Atomic layer deposition (ALD) enables conformal coatings of thin films and highly dispersed nanoparticles on high aspect ratio surfaces using sequential, self-limiting surface reactions.9−14 ALD allows the nanoparticle size and composition to be controlled precisely by adjusting the number and sequence of ALD cycles. This unique feature is attractive for synthesizing heterogeneous catalysts comprised of highly dispersed catalytic materials on high surface area, mesoporous supports. By combining the ALD for metal oxides with noble metal ALD, it is possible to engineer nanocatalysts with unique structures and properties. As a result, ALD has emerged as a promising method to prepare advanced, nanostructured catalysts.6,7,15−17 Inspired by its successful application in metal−organic chemical vapor deposition (MOCVD),18 palladium hexafluoroacetylacetonate (Pd(hfac)2) has been adapted to Pd ALD by alternating exposures to this precursor and formalin. This process has been shown to yield well-controlled Pd nanoparticles on metal oxides under mild synthesis conditions.19−21 Nanostructured ALD Pd catalysts prepared by this method are highly active for alkane dehydrogenation,7,22 methanol decomposition,21,23 and alcohol oxidation for fuel cells.24 The continued development of synthetic methods for preparing supported Pd nanoparticles with controlled size and structure relies on a fundamental understanding of the © XXXX American Chemical Society

mechanism for Pd ALD on metal oxide surfaces. However, the detailed evolution of the Pd surface species during the individual half-reactions, as well as the subsequent nucleation and growth of palladium nanoparticles, is still not clear. Transformation of organometallics to supported precious metal nanoparticles has been intensively studied on both planar supports under ultrahigh-vacuum conditions25−28 and high surface area supports as heterogeneous catalysts.29−31 The combination of in situ characterization techniques and density function theory (DFT) calculations provides a powerful means to determine the atomic arrangements in nanostructured materials and elucidate the surface reaction mechanisms.6,32−35 For instance, in situ Fourier transform infrared spectroscopy (FT-IR) is useful for identifying the surface species formed during each ALD half-reaction.36,37 X-ray absorption spectroscopy (XAS) allows one to probe the evolution of surface species under practical ALD operation conditions in real time.38 It provides crucial information about nearest neighbors, bond distances, and coordination numbers. In addition to XAS, total scattering and pair distribution function (PDF) provide quantitative structural information on both short- and longrange order, making it a powerful tool to probe nanoparticle nucleation and growth.39,40 Here we combine in situ FT-IR, Received: February 1, 2013 Revised: May 3, 2013

A

dx.doi.org/10.1021/jp401196f | J. Phys. Chem. C XXXX, XXX, XXX−XXX

The Journal of Physical Chemistry C

Article

Pd ALD was performed in a stainless steel continuous-flow reactor (1 in. o.d., 10 in. length) attached to a portable ALD system at the beamline. Kapton foil was used as the X-ray window material. All gas lines were 1/4 in. stainless steel tubes and heated accordingly. A vacuum pump was used to keep a base pressure at ∼3 Torr with 50 sccm ultrahigh-purity N2 flow. The high surface area support used for the in situ XAS experiment was five ALD cycles of TiO2-coated SiO2 gel (5c TiO2/SiO2). The as-prepared TiO2-coated SiO2 gel was pressed as a pellet in a cylindrical stainless steel sample holder (8 mm i.d.) for the in situ XAS experiment. Two Pd ALD cycles were conducted for the in situ XAS study. Each cycle was comprised of four steps: (1) Pd(hfac)2 exposure at 100 °C for 60 min, (2) purge with ultrahigh-purity N2 at 100 °C for 10 min, (3) reduction with formalin vapor at increasing temperature from 100 to 225 °C at every 25 °C, (4) purge with N2 and cooldown to 100 °C. XANES and EXAFS spectra were taken in real time. The final Pd loading on the 5c TiO2/SiO2 was determined using XRF. The absolute amounts of Pd species during the Pd ALD process were back-calculated by multiplying the Pd fraction and the final Pd loading. XAS Data Analysis. The Pd surface species fraction was calculated by conducting Pd XANES linear combination fittings using Athena in the IFEFFIT (version 1.2.11) package. PdO, Pd(hfac)2, Na2PdCl4, and Pd nanoparticles (Pd0) were considered as references and measured at the beamline. Appropriate XANES references were used for different reaction regimes. Pd XANES linear combination fittings at each data point were also cross-checked by EXAFS data fittings. EXAFS data fittings were performed using WinXAS (version 3.1). Reference phase and amplitude files for Pd−Pd were obtained from the Pd foil EXAFS spectrum. Reference phase and amplitude files for Pd−O and Pd−C were both taken from PdO EXAFS spectra. Reference phase and amplitude files for Pd−Cl were obtained from Na2PdCl4. A single shell model fit of the forward and inverse k2-weighted EXAFS data was obtained between k = 2.9−11.7 Å−1 and r = 1.3−2.8 Å, respectively. Data quality was checked by performing k1- and k3-weighted EXAFS fittings with the same parameters. In Situ Pair Distribution Function Measurements. After Pd(hfac)2 was dosed to 5c TiO2/SiO2 at 100 °C in the first half cycle of Pd ALD, the sample was carefully removed from the ALD chamber and transferred to the X-ray beamline for in situ pair distribution function (PDF) measurements. This sample was loaded in a Kapton capillary as a flow cell which allows flow of reactant gas during the in situ PDF measurements.43 Temperature was precisely controlled by an Oxford Cryosystems unit. PDF data was collected: (1) in He flow at room temperature (RT); (2) elevating temperature at heating rate of 6 °C/min from RT to 227 °C in 20 sccm 3.5% H2; (3) in He at RT. Scattering data for PDF were collected at beamline 11-ID-B at the APS at ANL. High-energy X-rays (58 keV, λ = 0.2127 Å) were used in combination with a large amorphous silicon-based area detector. The 2-D images were transformed to 1-D scattering intensity data using software Fit2D. The structure function S(Q) was obtained within software PDFgetX2. Direct Fourier transform of the reduced structure function F(Q) = Q[S(Q) − 1] led to the pair distribution function G(r). Contributions from support 5c TiO2/SiO2 were measured and subtracted to get differential PDFs. The resulting PDFs showed contributions only from Pd-atom correlations.

XAS, and PDF along with DFT calculations to compose a detailed picture of Pd ALD on TiO2.



MATERIALS AND METHODS ALD Pd/TiO2/SiO2 Synthesis. The ALD samples for ex situ characterization were synthesized in a hot-walled viscous flow reactor that has been described in detail elsewhere.41 The high surface area support used in this work was Silicycle S10040 M silica gel with ∼100 m2/g surface area. For each batch, 0.5 g of SiO2 gel was used for TiO2 ALD. The SiO2 gel surface was modified using five ALD cycles of TiO2 with a timing sequence (100−200−120−200). The TiO2 ALD used alternating exposures to TiCl4 (Sigma-Aldrich, 99.9%) and deionized water at 150 °C. Five TiCl4−H2O ALD cycles yielded ∼8 wt % TiO2 over SiO2. The thickness of TiO2 is equivalent to ∼2 Å with density of 4.23 g/cm3 and a growth rate of ∼0.4 Å/cycle, consistent with the TiO2 ALD growth rate reported using TiCl4 and water as reactants on ZrO2 particles.42 Pd ALD used alternating exposures to palladium hexafluoroacetylacetonate (Pd(hfac)2, Sigma-Aldrich, 99.9%) and formalin (Sigma-Aldrich, HCHO 37 wt % in H2O with methanol added to limit the formation of paraformaldehyde) at 200 °C. Pd/TiO2/SiO2 catalysts were prepared using one Pd ALD cycle on the ALD TiO2-coated SiO2. For each preparation, 0.5 g of TiO2/SiO2 was used as support for Pd ALD (300−300−300−300). The Pd loading after one Pd ALD cycle is typically about ∼1 wt %. The 2 wt % Pd/SiO2 samples were prepared to serve as a reference for the XAS measurements. The substrate used for this synthesis was Davisil silica gel (Sigma-Aldrich, 35−60 mesh), and incipient wet impregnation (IWI) was the method of synthesis and Pd(NH3)4(NO3)2 was used as the precursor. In Situ FT-IR Experiments. Fourier transform infrared spectroscopy (FT-IR) experiments were conducted in a specially designed in situ ALD reaction cell, using a substrate comprised of zirconia nanopowder pressed into a metal grid. The FT-IR data were recorded using a Thermo Scientific Nicolet 6700 FT-IR spectrometer. A liquid nitrogen cooled IR detector was housed in a metal box purged with ultrahighpurity nitrogen to prevent water condensation. A typical FT-IR spectrum was recorded at 4 cm−1 resolution using 256 scans. The ZrO2 nanopowder was first coated with TiO2 ALD using six cycles of TiCl4−H2O at 150 °C. Next, Pd ALD was performed on the TiO2-coated ZrO2 nanopowder. Pd(hfac)2 was first deposited on the TiO2 surface at 100 °C until the surface was fully saturated. After a N2 purge, HCHO was dosed at each temperature for 10 s from 100 to 225 °C at every 25 °C step. In the final step, the sample was exposed to HCHO for 10 min at 225 °C. FT-IR spectra were recorded at each temperature. Pd ALD was also performed on chlorine-free TiO2 /ZrO2 surface as a control experiment. This TiO2 was prepared using 12 cycles of titanium isopropoxide (TTIP) and H2O at 150 °C. In Situ X-ray Absorption Spectroscopy. X-ray absorption spectroscopy (XAS) measurements at the Pd K edge (∼24.35 keV), including extended X-ray absorption fine structure spectroscopy (EXAFS) and X-ray absorption near edge structure spectroscopy (XANES), were performed at the bending-magnet beamline (10-BM) of the Materials Research Collaborative Access Team (MR-CAT) at the Advanced Photon Source (APS) at Argonne National Laboratory (ANL). XAS data were collected in the transmission mode. B

dx.doi.org/10.1021/jp401196f | J. Phys. Chem. C XXXX, XXX, XXX−XXX

The Journal of Physical Chemistry C

Article

Ex Situ Characterization. Scanning transmission electron microscopy (STEM) and atomic number contrast high angle annular dark field (HAADF) imaging were performed using a probe aberration corrected 200 kV JEOL JEM-ARM200CF STEM/TEM with a resolution of less than 0.08 nm at Research Resources Center at University of Illinois at Chicago. 0.8 Å probe with 20 and 30 μm aperture was used for low- and highmagnification STEM images, respectively. The histogram of particle size distribution was generated from the STEM images using software ImageJ.44 Care was taken to count and measure measure only those particles that appeared in focus. Multiple images were taken, and more than 400 particles were measured to obtain reasonable statistics. Pd metal loadings were determined by X-ray fluorescence spectroscopy (XRF, Oxford ED2000) and inductively coupled plasma (ICP, Varian VistaMPX instrument). The composition of the as-prepared TiO2coated SiO2 was analyzed with an X-ray photoemission spectroscopy (XPS, Kratos AXIS-165) system using a monochromatic Al Kα (1486.6 eV) X-ray source operating at 15 kV and 10 mA. Density Functional Theory Calculation. Periodic density functional theory (DFT) calculations were carried out using the Vienna Ab Initio Simulation Package (VASP).45 The electronic exchange and correlation contributions were calculated using the generalized gradient PBE functional.46 The projector augmented wave (PAW) method was used to treat the explicit valence electrons and the ironic cores,47 and a plane-wave kinetic energy cutoff of 400 eV was used. The stoichiometric rutile TiO2(110) surface was modeled using a four-layer slab, which has a dimension of 11.76 Å along the [001] direction and 13.06 Å, along the [1i0] ̅ directions, respectively. A vacuum space of ∼50 Å was used between successive slabs in the periodic cell. The lattice parameters used for this work are a = 4.62 Å and c = 2.94 Å, respectively, from the relaxation of the bulk rutile TiO2 structure, and are in good agreement with the experimental values, where a = 4.594 Å and c = 2.959, as reported by Baur.48 The top two layers were relaxed. The optimizations of the slab and the surface intermediates were performed using the conjugate gradient algorithm until the force on each atom is converged within 0.05 eV/Å. The Brillouin zone was sampled using a 2 × 2 × 1 kpoint prid based on the Monkhorst−Pack scheme.49

Figure 1. Aberration-corrected HAADF-STEM images of as-prepared ALD 1 wt % Pd/TiO2/SiO2 (a) at low and (b) high resolution with inserted FFT.

Figure 2. FT-IR spectra after Pd(hfac)2 adsorption on TiO2 (asdeposited), and following 10 s exposures to HCHO at different temperatures as indicated.

adsorption, two types of hfac ligands are seen: Pd(hfac)* species (1605 cm−1) and Ti(hfac)* species (1652 cm−1).36,37,55 This suggests that some of the hfac ligands from the Pd(hfac)2 precursor adsorbed on the TiO2 surface. The initial 10 s formalin exposure at 100 °C completely eliminated the Pd(hfac)* species and slightly increased the Ti(hfac)* species (Figure S1). As the HCHO exposure temperature was increased, the Ti(hfac)* species slowly disappeared, but it took more than 10 min of HCHO exposure at 225 °C to completely remove the Ti(hfac)* species. Meanwhile, some of the Ti(hfac)* species decomposed to Ti(COCF3)* species and Ti(CF3)* as demonstrated by the new features that appeared at 1376 and 1205 cm−1 assigned to C−C and CF3 stretching vibrations.55,56 Under identical conditions, Pd ALD was also performed on non-chlorine containing TiO2. The non-chlorine containing TiO2 was prepared by 12 ALD cycles of of titanium isopropoxide (TTIP) and H2O at 150 °C over ZrO2 surface. IR features corresponding to both Ti(hfac)* and Pd(hfac)* features were observed. This indicates that free hfac ligands deligated from Pd(hfac)2 indeed tend to bind to the TiO2 surface regardless of the presence of chlorine. However, without surface chlorine, the hfac ligands were only completely removed from Pd(hfac)* by exposure to formalin at 150 °C, 50 °C higher than the surface with chlorine. Since Pd nucleation is most likely after the deligation of Pd(hfac)*, this suggests that the presence of surface chlorine accelerates Pd agglomeration, which is consistent with a previous report.57 The role of surface chlorine in Pd ALD is further studied by in situ EXAFS and density function theory calculations. In Situ XAS. In analogy with Pd ALD on hydroxylated Al2O3,19 we might expect Ti−O−Pd(hfac)* to be the product of the first half-reaction. As we show below, however, in situ



RESULTS AND DISCUSSION TiO2 is a one of the most frequently used support materials in catalysis. Commercial TiO2 is normally synthesized by the hydrolysis of titanium tetrachloride (TiCl4) and typically contains chlorine impurities (e.g., 0.3% HCl for Degussa P25).50,51 Residual chlorine could inhibit the reduction and lead to agglomeration of the precious metal.52,53 To model the chlorine containing TiO2, five ALD cycles of TiCl4 and H2O were used to coat a SiO2 template to serve as the substrate for Pd ALD. We measured a chlorine concentration in the ALD TiO2 of ∼2 wt % using X-ray photoelectron spectroscopy, in agreement with previous reports.54 As shown in Figure 1a, one ALD Pd cycle on the TiO2-coated SiO2 yielded highly dispersed Pd nanoparticles with a diameter of 1.0 ± 0.3 nm. The atomic resolution aberration-corrected HAADF-STEM image revealed (Figure 1b) crystalline Pd nanoparticles with 0.22 nm fringes. In Situ FT-IR. Figure 2 shows FT-IR spectra recorded during the first Pd ALD cycle on the ALD TiO2 surface at 100 °C (peak assignments are given in Table S1). Upon Pd(hfac)2 C

dx.doi.org/10.1021/jp401196f | J. Phys. Chem. C XXXX, XXX, XXX−XXX

The Journal of Physical Chemistry C

Article

clearly indicates a second scatterer in the first shell. A structure of Pd coordinating with two oxygen and two chloride atoms was obtained from quantitative EXAFS fitting (see Figure S2 and Table 2). Thus, Pd(hfac)* adsorbs on the TiO2 during the first half-reaction through surface chlorine to form Pd(hfac)Cl2*. In agreement with this interpretation, DFT calculations found that formation of Pd(hfac)Cl2* is more favorable than Pd(hfac)O2* on TiO2 in the presence of surface chlorine. Shown in Figure 4, gas phase Pd(hfac)2, two clean TiO2(110) surfaces, and two adsorbed Cl atoms are used as the reference state. The numbers in the parentheses represent the calculated potential energies based on the reference. Instead of the O sites on the TiO2(110) surface, Pd(hfac)* binds much stronger with surface Cl species by forming Pd(hfac)Cl2* species. The stable Pd(hfac)Cl2* species as an intermediate indicates this reaction is a surface self-limiting reaction. A recent DFT study also showed that surface species, including H, OH, and residual Cl, can mobilize to supported Pt nanoclusters and stabilize them on the Al2O3 surface.59 In accord with the FT-IR measurements, XAS showed that the Pd easily deligated at the lowest HCHO exposure temperature of 100 °C, leaving the Pd bound to nearby surface chlorine and oxygen. During the subsequent HCHO exposures at higher temperatures, EXAFS (Table 2) and DFT calculations (Figure 4) suggested that Pd(II) tended to bind to surface Cl if available. DFT calculations also show that on the TiO2(110) surface Pd is able to bind one additional Cl atom forming PdOCl3* species, which is slightly more stable than PdO2Cl2* species, consistent with the EXAFS data fitting. Figure 3b and Table 2 show that Pd−O bond cleavage started as low as 100 °C and was complete at 150 °C, as CNPd−O decreased to null. Subsequently, Pd−Cl bond cleavage began at 150 °C and completed at 200 °C. With the loss of Pd−Cl bonding, Pd atoms became mobile and formed small nanoparticles, leading to the appearance of Pd−Pd bonding at ∼2.5 Å (uncorrected distance) at 150 °C. Figure 3d shows XANES spectra recorded during the formalin exposure in the first Pd ALD cycle. A continuous narrowing of peaks accompanied by a shift of the edge to lower energy signals the gradual conversion of surface Pd species to metallic Pd. The final spectrum is similar to the XANES reference spectrum for Pd NPs shown in Figure 3c. After complete reduction at 200 °C, the EXAFS fitting revealed a Pd−Pd bond distance of 2.78 Å, indicating lattice expansion in the Pd−Pd bond distance compared to the bulk Pd value of 2.75 Å. It is known that metal nanoparticles of a few nanometers size undergo lattice contraction to minimize under coordination.60,61 The lattice expansion of Pd in our study may result from the formation of a palladium hydride or palladium carbide. As shown in Figure 5a, XANES of the ALD Pd showed broadening at the peak of the 1s → 5p transition as compared to the XANES from metallic Pd nanoparticles synthesized from a noncarbon containing precursor. Pd K edge Fourier transforms of the ALD Pd shifted 0.13 Å to higher R in Figure 5b. Consequently, both XANES and Fourier transform suggest the formation of a Pd−C phase.62,63 The formation of Pd−C might have resulted from the very long HCHO exposures used in this experiment. The Pd−C coordination number was determined to be 0.9 after the first Pd ALD cycle. In the second Pd ALD cycle (Figure S3), the asdeposited spectrum consisted of 1−2 Å bonds and Pd−Pd bonds at 2.49 Å (uncorrected distance), representing a mixture of both adsorbed Pd precursor and agglomerated Pd nano-

XAS measurements revealed that residual chlorine on the TiO2 surface completely changes the nature of the Pd(hfac)* adsorption. To facilitate interpretation of the in situ XAS measurements, XAS reference spectra were recorded for PdO, Pd(hfac)2, Na2PdCl4, and Pd NPs as shown in Figure 3. The

Figure 3. Pd K-edge (a, b) EXAFS Fourier transform and (c, d) XANES of (a, c) reference samples and (b, d) Pd surface species during the formalin exposure in the first Pd ALD cycle on an ALD TiO2 surface.

peaks at 1.53, 1.82, 2.42, and 2.95 Å in Figure 3a are attributed to the Pd−O, Pd−Cl, Pd−Pd, and Pd−(O)−Pd bonds, respectively. The corresponding data analysis is listed in Table 1. PdO is known to have four oxygen atoms as nearest Table 1. EXAFS Data Fittings of Reference Samplesa sample reference PdO reference Pd(hfac)2 reference Na2PdCl4 reference Pd NPs

scatter

CN

R, Å

DWF (×103)

E0, eV

assigned structure

Pd−O

4.0

2.05

0.0

−0.4

Pd−O

4.0

2.00

0.7

3.4

Pd(hfac)2

Pd−Cl

4.0

2.31

0.3

0.3

Na2PdCl4

Pd−Pd

9.0

2.75

1.5

2.4

Pd NPs

PdO

a

Pd K edge was measured. CN is coordination numbers with error ±10%. R is bond distance with error ±0.02 Å. DWF is Debye−Waller factor. E0 is energy shift. A two-shell model fit of the k2-weighted EXAFS data was obtained between k = 2.9−11.7 Å−1 and r = 1.3−2.8 Å.

neighbors, with a bond distance of ∼2.05 Å. Similarly, as a monometallic organometallic compound, Pd(hfac)2 is also surrounded by four oxygen atoms as nearest-neighbors with a slightly shorter bond distance of 2.00 Å. The Pd−Cl bond distance of 2.31 Å has been obtained in Na2PdCl4, close to that observed in D4h symmetry [PdCl4]2− and octahedral symmetry [PdCl6]2−.58 A H2 reduced 2 wt % Pd/SiO2 sample prepared by incipient wet impregnation was used as reference Pd NPs where the bond distance of 2.75 Å was obtained. The Pd K-edge Fourier transform for the as-deposited Pd species (Figure 3b) shows two peaks in the 1−2 Å regime. This D

dx.doi.org/10.1021/jp401196f | J. Phys. Chem. C XXXX, XXX, XXX−XXX

The Journal of Physical Chemistry C

Article

Table 2. EXAFS Data Fittings of Samples during Pd ALDa sample

scatter

CN

R, Å

DWF (×103)

E0, eV

after 1st Pd(hfac)2 deposition

Pd−Cl Pd−O Pd−Cl Pd−O Pd−Cl Pd−O Pd−Cl Pd−Pd Pd−Cl Pd−Pd Pd−C Pd−Pd Pd−C Pd−Pd Pd−C Pd−Pd

2.0 2.0 2.7 1.3 2.6 1.0 2.3 1.9 1.6 3.8 1.0 6.2 0.9 6.1 0.9 6.9

2.29 2.00 2.29 2.02 2.30 2.00 2.29 2.76 2.30 2.79 2..20 2.78 2.19 2.78 2.19 2.77

1.0 1.0 1.0 1.0 1.0 1.0 1.0 6.0 1.0 6.0 1.0 6.0 1.0 6.0 1.0 6.0

−2.7 1.0 −4.0 7.3 −0.3 −1.4 −1.4 0.5 −2.8 2.7 8.7 0.4 6.5 0.4 7.2 −0.1

1st ALD cycle at 100 °C 1st ALD cycle at 125 °C 1st ALD cycle at 150 °C 1st ALD cycle at 175 °C 1st ALD cycle at 200 °C 1st ALD cycle at 225 °C after 2nd ALD cycle

assigned structure Pd(hfac)Cl2* PdO1.3Cl2.7* PdO1.0Cl2.6* Pd/PdCx + PdCl2.4* Pd/PdCx + PdCl1.6* Pd/PdCx Pd/PdCx Pd/PdCx

Pd K edge was measured. CN is coordination numbers with error ±10%. R is bond distance with error ±0.02 Å. DWF is Debye−Waller factor. E0 is energy shift. A two-shell model fit of the k2-weighted EXAFS data was obtained between k = 2.9−11.7 Å−1 and r = 1.3−2.8 Å. a

Figure 4. DFT calculations of the potential energies of Pd(hfac)2 deligation during Pd ALD and optimized structures of surface intermediates on stoichiometric rutile TiO2 (110). Gray atoms represent Ti, red atoms represent O, green atoms represent Cl, light blue atoms represent F, dark blue atoms represent Pd, dark gray atoms represent C, and white atoms represent H.

Pd(hfac)Cl2* species were only completely reduced to Pd after the 200 °C HCHO treatment. In Situ PDF. In situ XAS has been successfully employed to revolve the nucleation and growth of precious metals. However, EXAFS data require intensive data fitting which is not trivial for very small targets under a complicated chemical environment. To understand particle growth during the Pd ALD, we performed in situ PDF measurements. Quantitative differential PDF (d-PDF), where the contribution from the 5c TiO2/SiO2 support is precisely subtracted from the total PDF, provides only the atom−atom correlations from the supported Pd species (e.g., Pd−Pd).39,64 Figure 7 shows no obvious changes between RT and 99 °C. Pd aggregation started at 105 °C with the appearance of peaks corresponding to the Pd−Pd shells. This aggregation temperature coincides well with the Pd(hfac)Cl2* deligation temperature seen in XAS and the loss of Pd(hfac) vibrational features observed in FT-IR. Both EXAFS data fitting and DFT calculations suggested that most Pd was

particles. After the second Pd ALD cycle, the CN for Pd−Pd increased only marginally from 6.1 to 6.9, implying a small increase in Pd particle size. This observation is consistent with the previous STEM results that the size of Pd nanoparticles increases only slightly after the first few Pd ALD cycles.20 XANES spectra of Pd(hfac)Cl2* and Pd/PdCx were used in linear combination fitting (LCF) to the XANES spectra during the second half-reaction. Figure 6 and Figure S4 illustrate the Pd surface species evolution for two cycles of Pd ALD. As expected, the Pd(hfac)Cl2* species decreased continuously with increasing HCHO exposure temperature. After the initial exposure to Pd(hfac)2 at 100 °C, the TiO2 surface was saturated with 11 μmol of Pd(hfac)Cl2. Only ∼9 μmol of Pd(hfac)Cl2 was deposited during the Pd(hfac)2 exposure in the second cycle. The lower uptake in the second cycle may be due to active sites blocked by residues from the hfac ligands as suggested by the in situ FT-IR experiments. The surface E

dx.doi.org/10.1021/jp401196f | J. Phys. Chem. C XXXX, XXX, XXX−XXX

The Journal of Physical Chemistry C

Article

Figure 7. PDFs of Pd(hfac)Cl2* annealed in 3.5% H2/He.



CONCLUSIONS In summary, we systematically investigated the initial stages of Pd ALD on an ALD TiO2 surface with residual chlorine using in situ FT-IR, in situ XAS, and in situ PDF. During the first halfreaction, Pd(hfac)2 decomposed on the TiO2 surface to form Pd(hfac)Cl2* and Ti(hfac)*. In the following HCHO exposure at 100 °C, the Pd(hfac)Cl2* rapidly lost its remaining hfac ligand and bound to the TiO2 through surface chlorine and oxygen. This transformation was accompanied by some Pd nanoparticle agglomeration. Increasing the HCHO exposure temperature caused first Pd−O bond cleavage followed Pd−Cl bond cleavage. Following the Pd−Cl bond cleavage, most of the Pd atoms rapidly agglomerated to form ∼1 nm nanoparticles. The adsorbed Pd precursors were reduced to Pd nanoparticles at 200 °C with a minor Pd−C phase. Repeated Pd ALD cycles nearly doubled the Pd loading without increasing the particle size. As a result, a detailed understanding of heterogeneous Pd/TiO2 catalyst synthesis using Pd(hfac)2 under practical ALD conditions was achieved. This knowledge can help direct the future design of advanced heterogeneous catalysts from organometallic compounds.

Figure 5. Pd K-edge XANES (a) and EXAFS Frourier transform (b) of PdCx (black) and Pd nanoparticles (red).



ASSOCIATED CONTENT

S Supporting Information *

Additional details on FT-IR spectra and assignment, EXAFS spectra, and XANES linear combination fitting. This material is available free of charge via the Internet at http://pubs.acs.org.



AUTHOR INFORMATION

Corresponding Author

*Phone 1-630-252-3520, e-mail [email protected] (J.W.E.). Notes

The authors declare no competing financial interest.



Figure 6. Absolute amount of Pd surface species determined by XANES LCF using Pd(hfac)Cl2* and PdCx as reference.

ACKNOWLEDGMENTS This material is based upon work supported as part of the Institute for Atom-efficient Chemical Transformations (IACT), an Energy Frontier Research Center funded by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences. Use of the Advanced Photon Source was supported by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, under Contract DE-AC02-06CH11357. MRCAT operations are supported by the Department of Energy and the MRCAT member institutions. We gratefully acknowledge the computing resources provided on “Fusion”, a

still bound to TiO2 through Pd−O and Pd−Cl bonds after HCHO reduction at 100 °C. However, some Pd aggregation already occurred once the hfac ligand detached from the Pd atoms. The peaks at 2.78, 4.83, and 7.37 Å are slightly larger than 2.75, 4.76, and 7.28 Å for the first, third, and sixth Pd−Pd shell of fcc bulk Pd, respectively. This is probably due to the formation of palladium hydride.65 The peak intensity further increased with increasing annealing temperature, suggesting an increase in the concentration of Pd nanoparticles. F

dx.doi.org/10.1021/jp401196f | J. Phys. Chem. C XXXX, XXX, XXX−XXX

The Journal of Physical Chemistry C

Article

Type Atomic Layer Deposition: Growth and Morphology. Langmuir 2010, 26 (21), 16486−16495. (21) Feng, H.; Elam, J. W.; Libera, J. A.; Setthapun, W.; Stair, P. C. Palladium Catalysts Synthesized by Atomic Layer Deposition for Methanol Decomposition. Chem. Mater. 2010, 22 (10), 3133−3142. (22) Lei, Y.; Liu, B.; Lu, J. L.; Lobo-Lapidus, R. J.; Wu, T. P.; Feng, H.; Xia, X. X.; Mane, A. U.; Libera, J. A.; Greeley, J. P.; et al. Synthesis of Pt-Pd Core-Shell Nanostructures by Atomic Layer Deposition: Application in Propane Oxidative Dehydrogenation to Propylene. Chem. Mater. 2012, 24 (18), 3525−3533. (23) Feng, H.; Libera, J. A.; Stair, P. C.; Miller, J. T.; Elam, J. W. Subnanometer Palladium Particles Synthesized by Atomic Layer Deposition. ACS Catal. 2011, 1 (6), 665−673. (24) Rikkinen, E.; Santasalo-Aarnio, A.; Airaksinen, S.; Borghei, M.; Viitanen, V.; Sainio, J.; Kauppinen, E. I.; Kallio, T.; Krause, A. O. I. Atomic Layer Deposition Preparation of Pd Nanoparticles on a Porous Carbon Support for Alcohol Oxidation. J. Phys. Chem. C 2011, 115 (46), 23067−23073. (25) Purnell, S. K.; Xu, X.; Goodman, D. W.; Gates, B. C. Adsorption and Reaction of [Re2(CO)10] on Ultrathin MgO Films Grown on a Mo(110) Surface: Characterization by Infrared Reflection-Absorption Spectroscopy and Temperature Programmed Desorption. J. Phys. Chem. 1994, 98 (15), 4076−4082. (26) Chen, Z. W.; Fujita, S.; Fukui, K. Scanning Tunneling Microscopy Study on the Precursor-Dependent Formation of Homogeneous Rh Clusters on Al2O3/NiAl(110). J. Phys. Chem. C 2011, 115 (29), 14270−14277. (27) Lei, Y.; Uhl, A.; Becker, C.; Wandelt, K.; Gates, B. C.; Meyer, R.; Trenary, M. Adsorption and Reaction of Rh(CO)2(acac) on Al2O3/ Ni3Al(111). Phys. Chem. Chem. Phys. 2010, 12 (6), 1264−1270. (28) Yang, F.; Trufan, E.; Adams, R. D.; Goodman, D. W. Structure of Molecular-Sized Ru3Sn3 Clusters on a SiO2 Film on Mo(112). J. Phys. Chem. C 2008, 112 (37), 14233−14235. (29) Goellner, J. F.; Gates, B. C. Synthesis and Characterization of Site-Isolated Hexarhodium Clusters on Titania Powder. J. Phys. Chem. B 2001, 105 (16), 3269−3281. (30) Laurent, P.; Veyre, L.; Thieuleux, C.; Donet, S.; Coperet, C. From Well-Defined Pt(II) Surface Species to the Controlled Growth of Silica Supported Pt Nanoparticles. Dalton Trans. 2013, 42 (1), 238−248. (31) Aydin, C.; Lu, J.; Browning, N. D.; Gates, B. C. A “Smart” Catalyst: Sinter-Resistant Supported Iridium Clusters Visualized with Electron Microscopy. Angew. Chem., Int. Ed. 2012, 51 (24), 5929− 5934. (32) Aramburo, L. R.; de Smit, E.; Arstad, B.; van Schooneveld, M. M.; Sommer, L.; Juhin, A.; Yokosawa, T.; Zandbergen, H. W.; Olsbye, U.; de Groot, F. M. F.; et al. X-ray Imaging of Zeolite Particles at the Nanoscale: Influence of Steaming on the State of Aluminum and the Methanol-To-Olefin Reaction. Angew. Chem., Int. Ed. 2012, 51 (15), 3616−3619. (33) Somorjai, G. A.; Contreras, A. M.; Montano, M.; Rioux, R. M. Clusters, Surfaces, and Catalysis. Proc. Natl. Acad. Sci. U. S. A. 2006, 103 (28), 10577−10583. (34) Sun, Y. N.; Giordano, L.; Goniakowski, J.; Lewandowski, M.; Qin, Z. H.; Noguera, C.; Shaikhutdinov, S.; Pacchioni, G.; Freund, H. J. The Interplay between Structure and CO Oxidation Catalysis on Metal-Supported Ultrathin Oxide Films. Angew. Chem., Int. Ed. 2010, 49 (26), 4418−4421. (35) Kantorovich, L. N.; Shluger, A. L.; Sushko, P. V.; Gunster, J.; Stracke, P.; Goodman, D. W.; Kempter, V. Mg Clusters on MgO Surfaces: Study of the Nucleation Mechanism with MIES and Ab Initio Calculations. Faraday Discuss. 1999, 114, 173−194. (36) Goldstein, D. N.; George, S. M. Surface Poisoning in the Nucleation and Growth of Palladium Atomic Layer Deposition with Pd(hfac)2 and Formalin. Thin Solid Films 2011, 519 (16), 5339−5347. (37) Goldstein, D. N.; George, S. M. Enhancing the Nucleation of Palladium Atomic Layer Deposition on Al2O3 Using Trimethylaluminum to Prevent Surface Poisoning by Reaction Products. Appl. Phys. Lett. 2009, 95 (14).

320-node computing cluster operated by the Laboratory Computing Resource Center at Argonne National Laboratory.



REFERENCES

(1) Haruta, M. Size- and Support-Dependency in the Catalysis of Gold. Catal. Today 1997, 36 (1), 153−166. (2) Argo, A. M.; Odzak, J. F.; Gates, B. C. Role of Cluster Size in Catalysis Spectroscopic Investigation of Gamma-Al2O3-supported Ir4 and Ir6 during Ethene Hydrogenation. J. Am. Chem. Soc. 2003, 125 (23), 7107−7115. (3) Fu, Q.; Saltsburg, H.; Flytzani-Stephanopoulos, M. Active Nonmetallic Au and Pt Species on Ceria-Based Water-Gas Shift Catalysts. Science 2003, 301 (5635), 935−938. (4) Lee, S. S.; Fan, C. Y.; Wu, T. P.; Anderson, S. L. CO Oxidation on Aun/TiO2 Catalysts Produced by Size-Selected Cluster Deposition. J. Am. Chem. Soc. 2004, 126, 5682−5683. (5) Guzman, J.; Carrettin, S.; Fierro-Gonzalez, J. C.; Hao, Y. L.; Gates, B. C.; Corma, A. CO Oxidation Catalyzed by Supported Gold: Cooperation between Gold and Nanocrystalline Rare-Earth Supports Forms Reactive Surface Superoxide and Peroxide Species. Angew. Chem., Int. Ed. 2005, 44 (30), 4778−4781. (6) Lei, Y.; Mehmood, F.; Lee, S.; Greeley, J.; Lee, B.; Seifert, S.; Winans, R. E.; Elam, J. W.; Meyer, R. J.; Redfern, P. C.; et al. Increased Silver Activity for Direct Propylene Epoxidation via Subnanometer Size Effects. Science 2010, 328 (5975), 224−228. (7) Lu, J. L.; Fu, B. S.; Kung, M. C.; Xiao, G. M.; Elam, J. W.; Kung, H. H.; Stair, P. C. Coking- and Sintering-Resistant Palladium Catalysts Achieved Through Atomic Layer Deposition. Science 2012, 335 (6073), 1205−1208. (8) Zhu, Y.; Qian, H. F.; Drake, B. A.; Jin, R. C. Atomically Precise Au25(SR)18 Nanoparticles as Catalysts for the Selective Hydrogenation of alpha,beta-Unsaturated Ketones and Aldehydes. Angew. Chem., Int. Ed. 2010, 49 (7), 1295−1298. (9) Puurunen, R. L. Surface Chemistry of Atomic Layer Deposition: A Case Study for the Trimethylaluminum/Water Process. J. Appl. Phys. 2005, 97 (12), 121301−121352. (10) Leskela, M.; Ritala, M. Atomic Layer Deposition Chemistry: Recent Developments and Future Challenges. Angew. Chem., Int. Ed. 2003, 42 (45), 5548−5554. (11) Marichy, C.; Bechelany, M.; Pinna, N. Atomic Layer Deposition of Nanostructured Materials for Energy and Environmental Applications. Adv. Mater. 2012, 24 (8), 1017−1032. (12) Elam, J. W.; Dasgupta, N. P.; Prinz, F. B. ALD for Clean Energy Conversion, Utilization, and Storage. MRS Bull. 2011, 36 (11), 899− 906. (13) Lim, B. S.; Rahtu, A.; Gordon, R. G. Atomic Layer Deposition of Transition Metals. Nat. Mater. 2003, 2 (11), 749−754. (14) Zaera, F. The Surface Chemistry of Atomic Layer Depositions of Solid Thin Films. J. Phys. Chem. Lett. 2012, 3 (10), 1301−1309. (15) King, J. S.; Wittstock, A.; Biener, J.; Kucheyev, S. O.; Wang, Y. M.; Baumann, T. F.; Giri, S. K.; Hamza, A. V.; Baeumer, M.; Bent, S. F. Ultralow Loading Pt Nanocatalysts Prepared by Atomic Layer Deposition on Carbon Aerogels. Nano Lett. 2008, 8 (8), 2405−2409. (16) Lu, J. L.; Lei, Y.; Elam, J. W. In Noble Metals; Su, Y.-H., Ed.; INTECH: Janeza Trdine, 2012; pp 159−178. (17) Liang, X. H.; Li, J. H.; Yu, M.; McMurray, C. N.; Falconer, J. L.; Weimer, A. W. Stabilization of Supported Metal Nanoparticles Using an Ultrathin Porous Shell. ACS Catal. 2011, 1 (10), 1162−1165. (18) Lin, W. B.; Warren, T. H.; Nuzzo, R. G.; Girolami, G. S. SurfaceSelective Deposition of Palladium and Silver Films from MetalOrganic Precursors: A Novel Metal-Organic Chemical Vapor Deposition Redox Transmetalation Process. J. Am. Chem. Soc. 1993, 115 (24), 11644−11645. (19) Elam, J. W.; Zinovev, A.; Han, C. Y.; Wang, H. H.; Welp, U.; Hryn, J. N.; Pellin, M. J. Atomic Layer Deposition of Palladium Films on Al2O3 Surfaces. Thin Solid Films 2006, 515 (4), 1664−1673. (20) Lu, J. L.; Stair, P. C. Nano/Subnanometer Pd Nanoparticles on Oxide Supports Synthesized by AB-Type and Low-Temperature ABCG

dx.doi.org/10.1021/jp401196f | J. Phys. Chem. C XXXX, XXX, XXX−XXX

The Journal of Physical Chemistry C

Article

Hexachloropalladate(IV), K2PdCl4 AND K2PdCl6 at 120 K. Acta Crystallogr., Sect. B: Struct. Sci. 1988, 44, 580−585. (59) Mager-Maury, C.; Chizallet, C.; Sautet, P.; Raybaud, P. Platinum Nanoclusters Stabilized on gamma-Alumina by Chlorine Used As a Capping Surface Ligand: A Density Functional Theory Study. ACS Catal. 2012, 2 (7), 1346−1357. (60) Miller, J. T.; Kropf, A. J.; Zha, Y.; Regalbuto, J. R.; Delannoy, L.; Louis, C.; Bus, E.; van Bokhoven, J. A. The Effect of Gold Particle Size on Au-Au Bond Length and Reactivity toward Oxygen in Supported Catalysts. J. Catal. 2006, 240 (2), 222−234. (61) Lei, Y.; Jelic, J.; Nitsche, L. C.; Meyer, R.; Miller, J. T. Effect of Particle Size and Adsorbates on the L3, L2 and L1 X-ray Absorption Near Edge Structure of Supported Pt Nanoparticles. Top. Catal. 2011, 54, 334−348. (62) McCaulley, J. A. Temperature Dependence of the Pd K-edge Extended X-ray Absorption Fine Structure of PdCx (x∼0.13). Phys. Rev. B 1993, 47 (9), 4873−4879. (63) Tew, M. W.; Nachtegaal, M.; Janousch, M.; Hutchwelker, T.; Van Bokhoven, J. A. The Irreversible Formation of Palladium Carbide during Hydrogenation of 1-pentyne over Silica-Supported Palladium Nanoparticles: In-Situ Pd K and L3 edge XAS. Phys. Chem. Chem. Phys. 2012, 14, 5761−5768. (64) Newton, M. A.; Chapman, K. W.; Thompsett, D.; Chupas, P. J. Chasing Changing Nanoparticles with Time-Resolved Pair Distribution Function Methods. J. Am. Chem. Soc. 2012, 134 (11), 5036−5039. (65) Bauer, M.; Schoch, R.; Shao, L. D.; Zhang, B. S.; Knop-Gericke, A.; Willinger, M.; Schlogl, R.; Teschner, D. Structure-Activity Studies on Highly Active Palladium Hydrogenation Catalysts by X-ray Absorption Spectroscopy. J. Phys. Chem. C 2012, 116 (42), 22375− 22385.

(38) Setthapun, W.; Williams, W. D.; Kim, S. M.; Feng, H.; Elam, J. W.; Rabuffetti, F. A.; Poeppelmeier, K. R.; Stair, P. C.; Stach, E. A.; Ribeiro, F. H.; et al. Genesis and Evolution of Surface Species during Pt Atomic Layer Deposition on Oxide Supports Characterized by in Situ XAFS Analysis and Water-Gas Shift Reaction. J. Phys. Chem. C 2010, 114 (21), 9758−9771. (39) Chupas, P. J.; Chapman, K. W.; Jennings, G.; Lee, P. L.; Grey, C. P. Watching Nanoparticles Grow: The Mechanism and Kinetics for the Formation of TiO2-Supported Platinum Nanoparticles. J. Am. Chem. Soc. 2007, 129 (45), 13822−+. (40) Zhao, H. Y.; Nenoff, T. M.; Jennings, G.; Chupas, P. J.; Chapman, K. W. Determining Quantitative Kinetics and the Structural Mechanism for Particle Growth in Porous Templates. J. Phys. Chem. Lett. 2011, 2 (21), 2742−2746. (41) Elam, J. W.; Groner, M. D.; George, S. M. Viscous Flow Reactor with Quartz Crystal Microbalance for Thin Film Growth by Atomic Layer Deposition. Rev. Sci. Instrum. 2002, 73 (8), 2981−2987. (42) Ferguson, J. D.; Yoder, A. R.; Weimer, A. W.; George, S. M. TiO2 Atomic Layer Deposition on ZrO2 Particles Using Alternating Exposures of TiCl4 and H2O. Appl. Surf. Sci. 2004, 226 (4), 393−404. (43) Chupas, P. J.; Ciraolo, M. F.; Hanson, J. C.; Grey, C. P. In Situ X-ray Diffraction and Solid-State NMR Study of the Fluorination of Gamma-Al2O3 with HCF2Cl. J. Am. Chem. Soc. 2001, 123 (8), 1694− 1702. (44) Figueroa, S. J. A.; Stewart, S. J.; Rueda, T.; Hernando, A.; de la Presa, P. Thermal Evolution of Pt-Rich FePt/Fe3O4 Heterodimers Studied Using X-ray Absorption Near-Edge Spectroscopy. J. Phys. Chem. C 2011, 115 (13), 5500−5508. (45) Kresse, G.; Furthmuller, J. Efficiency of Ab-Initio Total Energy Calculations for Metals and Semiconductors Using a Plane-Wave Basis Set. Comput. Mater. Sci. 1996, 6 (1), 15−50. (46) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77 (18), 3865− 3868. (47) Blochl, P. E. Projector Augmented Wave Method. Phys. Rev. B 1994, 50 (24), 17953−17979. (48) Grant, F. A. Properties of Rutile (Titanium Dioxide). Rev. Mod. Phys. 1959, 31 (3), 646−674. (49) Monkhorst, H. J.; Pack, J. D. Special Points for Brillouin Zone Integrations. Phys. Rev. B 1976, 13 (12), 5188−5192. (50) Lawhorne, E. R. Titanium Dioxide with High Purity and Uniform Particle Size and Method Therefore. US4944936, 1990. (51) Ettinger, M. Highly Dispersed Metallic Oxides Producted by the AEROSIL® Process; Degussa’s Technical Bulletin Pigments No. 56, Frankfurt, Germany, 1995. (52) Karhu, H.; Kalantar, A.; Vayrynen, I. J.; Salmi, T.; Murzin, D. Y. XPS Analysis of Chlorine Residues in Supported Pt and Pd Catalysts with Low Metal Loading. Appl. Catal., A 2003, 247 (2), 283−294. (53) Bus, E.; Prins, R.; van Bokhoven, J. A. Time-Resolved in Situ XAS Study of the Preparation of Supported Gold Clusters. Phys. Chem. Chem. Phys. 2007, 9 (25), 3312−3320. (54) Lu, J. L.; Kosuda, K. M.; Van Duyne, R. P.; Stair, P. C. Surface Acidity and Properties of TiO2/SiO2 Catalysts Prepared by Atomic Layer Deposition: UV-visible Diffuse Reflectance, DRIFTS, and Visible Raman Spectroscopy Studies. J. Phys. Chem. C 2009, 113 (28), 12412−12418. (55) Lin, W. B.; Wiegand, B. C.; Nuzzo, R. G.; Girolami, G. S. Mechanistic Studies of Palladium Thin Film Growth from Palladium(II) Beta-diketonates 0.1. Spectroscopic Studies of the Reactions of Bis(hexafluoroacetylacetonato)palladium(II) on Copper Surfaces. J. Am. Chem. Soc. 1996, 118 (25), 5977−5987. (56) Morris, M. L.; Moshier, R. W.; Sievers, R. E. Infrared Spectra of Metal Chelate Compounds of Hexafluoroacetylacetone. Inorg. Chem. 1963, 2 (2), 411−412. (57) Bartholomew, C. H. Mechanisms of Catalyst Deactivation. Appl. Catal., A 2001, 212 (1−2), 17−60. (58) Takazawa, H.; Ohba, S.; Saito, Y. Electron-Density Distribution in Crystals of Dipotassium Tetrachloropalladate(II) and Dipotassium H

dx.doi.org/10.1021/jp401196f | J. Phys. Chem. C XXXX, XXX, XXX−XXX