Nebojša
Zec
a,
Gaetano
Mangiapia
a,
Mikhail L.
Zheludkevich
bc,
Sebastian
Busch
*a and
Jean-François
Moulin
*a
aGerman Engineering Materials Science Centre (GEMS) at Heinz Maier-Leibnitz Zentrum (MLZ), Helmholtz-Zentrum Geesthacht GmbH, Lichtenbergstr. 1, 85748 Garching bei München, Germany. E-mail: sebastian.busch@hzg.de; jean-francois.moulin@hzg.de
bMagIC-Magnesium Innovation Center, Institute of Materials Research, Helmholtz-Zentrum Geesthacht, Max-Planck-Str. 1, 21502 Geesthacht, Germany
cFaculty of Engineering, University of Kiel, Kaiserstr. 2, 24143 Kiel, Germany
First published on 22nd May 2020
Deep eutectic solvents (DESs) are both green and sustainable, making them an increasingly attractive alternative to conventional solvents. One of their applications is the electrochemical deposition of metals that cannot be deposited from aqueous solution because of the limited electrochemical window of water. The electrodeposition process is influenced by the structure and dynamics of the solvent at the solid–liquid interface. Therefore,the nanoscale structure of the interface between a silicon substrate and deep eutectic solvent (choline chloride–ethylene glycol) was studied by neutron reflectometry (NR) and molecular dynamics (MD) simulations. It is not possible to model NR measurements of this system without simulating a dense DES layer at the solid–liquid interface. This study used an MD simulation trajectory to extract the density, thickness, and roughness of this DES layer. With this input, the model reproduces the reflectometry data at all measured H/D contrasts very well. The thickness of the layer does not change appreciably when applying charge or at higher temperatures. Further analysis revealed a reorganization of ions and reorientation of the choline cations in the interface layer when the electrodes are charged. These changes in ion orientation are not observed with the NR technique since they do not influence the neutron scattering length density profile due to the high number of ethylene glycol molecules at the interface. However, the agreement between measured neutron reflectometry data and model parameters obtained from MD simulations justified subnanoscale analysis of the MD trajectory and confirmed that these two complementary techniques can be successfully combined to reveal the solid/DES interface structure.
The focus of this work is to better understand the interfacial properties in DES-based systems used for the electrodeposition of metals. In particular, we are focusing on systems used for the electrodeposition of zinc (Zn). In a previous work by Starykevich et al.15 it was observed that the potentiostatic (DC) deposition of Zn on an alumina layer from an ethaline/ZnCl2 solution is kinetically hindered. It was hypothesized that the constituents of ethaline, choline cations, and ethylene glycol molecules adsorb on the alumina surface because of the interactions of the hydroxyl groups of these species with the oxygen-rich alumina surface. The hypothesized existence of this effectively passivating layer was supported by the observation that raising the deposition temperature or superimposing a rapidly alternating potential to the potentiostatic voltage made Zn deposition possible. Both increased temperature and oscillating potential were thought to be able to disrupt the interfacial solvent layer and provide the Zn ions with access to the surface. Indeed, adsorption of the ethaline components on different electrode materials has been proposed and investigated by several research groups. A compact layer on a glassy carbon (GC) electrode surface was proposed by Vieira et al.,16 investigating a ZnCl2 containing DES using polarization modulation reflection–absorption spectroscopy. Blocking of the active sites on the GC electrode by adsorbed species was also proposed by Whitehead et al.,17 based on electrochemical measurements as well as by Pereira et al.18 for a steel electrode.
In this work, we will address the very question of the liquid structure at the interface in order to understand the processes taking place during electrodeposition. Our approach is based on the complementary use of neutron reflectometry (NR) as a structural measurement technique with molecular dynamics (MD) simulations. In this two-fold approach, the simulations of the interface help to remove the ambiguity of the NR data while the simulations themselves are validated by checking the consistency of the structural profile with observed NR data.
NR is a powerful technique for probing structures along the direction perpendicular to a planar interface, i.e., for investigating composition and scattering length density profiles across interfaces. In contrast to local probes such as surface microscopy, the use of a beam to illuminate a macroscopic surface ensures that the neutron scattering method provides a statistically relevant signal that is immune to local defects. Reflectometry measurements revealing atomic ordering at surfaces can also be performed using X-rays as a probe, but using neutrons offers significant advantages. In contrast to X-rays, neutrons interact in a non-destructive fashion with the material under examination. This makes it easy to follow processes in situ over a longer period. Moreover, due to their weak interaction with matter, neutrons have a large penetration depth for most materials. Finally, the large neutron scattering length difference between hydrogen and deuterium makes the NR technique perfectly suited for the study of hydrogen-rich systems.
This is particularly exploited in soft condensed matter studies involving organic and biological samples. Using selective deuteration of a given component makes it possible to specifically adjust the scattering contrast with respect to the non-deuterated environment. It is thus possible to accurately locate each interfacial component while keeping the chemical structure at the interface essentially intact. In the work of Cooper et al.19 neutron reflectometry gave insight into the structure and composition of an IL layer formed at the oil–silica interface. They observed an 8–11 Å thick layer that remains in place even at temperatures up to 80 °C. Using the same technique, the structure of a pyrrolidinium-based ionic liquid was also investigated at a charged gold interface.20 These results indicate non-electrostatic adsorption of the cation onto the gold electrode. X-ray reflectometry (XRR) studies of different ionic liquid–solid interfaces were performed by several research groups21 confirming the layering of ions on sapphire,22,23 gold,24 and silicon25 surfaces. However, as stated in Hammond et al.,26 based on their investigation, it appears that conclusions obtained for ionic liquid systems will not necessarily be valid for DESs. Their atomic force microscopy study of a DES/platinum interface reveals that adding water to a DES has an unexpected effect on the interface, actually increasing the interfacial nanostructure by up to 40%. Therefore, more experimental, computational, and theoretical work is still needed in order to understand the interfacial properties of DESs on different solid materials.
Molecular dynamics (MD) simulations are a widely used computational technique that can give valuable information about the composition, structure, and dynamics of the layers formed at liquid–solid interfaces. Even though bulk DESs are often studied by MD simulations, literature on the DES–electrode interface is scarce. Atilhan et al.27 studied interfacial properties for different DES nanodroplets on graphene, utilizing density functional theory and classical MD. A combined X-ray reflectometry and MD simulations approach was used by Steinrück et al.28 to investigate electrolyte ordering near solid surfaces relevant to Li-ion batteries.
The main aim of this work is to reveal the structure of the layer formed at the interface between a DES and an electrode with and without applied constant electric potential using neutron reflectometry and MD simulations. As such, we have chosen to concentrate on a system where the electrode surface is well characterized and does not introduce additional complexity. Due to the availability of large single crystals whose surface can be polished, silicon was chosen as an electrode. Further development of this work will focus on other electrode materials with a more direct link to industrial applications.
Neutron reflectometry measurements were carried out at REFSANS, the horizontal time-of-flight reflectometer operated by the Helmholtz-Zentrum Geesthacht at the Heinz Maier-Leibnitz Zentrum in Garching, Germany.29 A polychromatic pulsed beam with wavelengths ranging between 3 Å and 21 Å and a wavelength resolution Δλ/λ ≈ 3% was used. A 10 m collimation system was used to define the vertical divergence of the beam. The neutron beam hit the sample from the silicon side as shown in Fig. 1. The neutron incident angles were adjusted by tilting the goniometer head on which the sample was mounted. A two-dimensional 256 × 256 pixel-based array 3He detector registered the neutrons scattered from the samples at a sample-to-detector distance of 10.38 m. Two different incident angles, namely θ1 = 0.50° and θ2 = 2.50°, allowed collecting data in the interval 0.005 Å−1 ≤ Qz ≤ 0.14 Å−1.
Beam size and measurement times were set according to the incident angle, to maximize the incident flux and increase the data statistics. Measurements of the direct beam were also performed: using the settings of the small angle, a direct beam was measured directly on the detector, whereas for the high incident angle a calibrated neutron monitor was used. Raw data were reduced and the experimental reflectivity R was obtained as a ratio of the reflected to the incident beam. The thickness of the SiO2 layer covering the wafer was determined by NR measurements of this waver vs. H2O and D2O and a simultaneous fit of both datasets to 1.0 nm. In all sample measurements, three different contrast mixtures of ethylene glycol (anhydrous, 99.8%) and choline chloride (≥99%) were used: deuterated ethylene glycol (98% atom D):hydrogenated choline chloride (D-EG:H-Cho, SLD = 3.58 × 10−6 Å−2), hydrogenated ethylene glycol:deuterated choline chloride (98% atom D) (H-EG:D-CHO SLD = 2.68 × 10−6 Å−2) and both hydrogenated compounds (H-EG:H-CHO, SLD = 0.28 × 10−6 Å−2). All chemicals were purchased from Sigma-Aldrich. For the application of the electric voltage, the two electrodes were connected to the Gamry Potentiostat: all potentials reported in this paper are measured vs. the aluminum electrode. All electrochemical experiments were performed at room temperature.
All-atom force field parameters (refined and improved OPLS parameters) for DES were taken from the work of Ferreira et al.34 while Lennard-Jones parameters for Si (σ = 0.3385 nm, ε = 2.4471 kJ mol−1) were taken from the standard OPLS force field.35,36 Knowing the bulk densities obtained from the NPT simulations, a new simulation box with fixed dimensions was constructed and confined by two silicon electrode sheets prepared with the help of VMD.37 Periodic boundary conditions (PBC) were applied in all dimensions and in order to neutralize any image charge effects, the box was expanded to 16 nm in one direction by introducing a region of vacuum beyond the silicon sheets (Fig. 2). In order to improve statistics twelve replicas of this system were simulated. Each of the replicas had the same initial configuration obtained using Packmol but the atoms were assigned different random velocities at the start of the equilibration. The first step of the simulation was an annealing process from 800–300 K for 5 ns in a stepwise fashion. After that 20 ns of NVT simulation were performed with 1 fs time steps. In the first set of simulations, the silicon electrodes were uncharged. In the second set to each silicon atom a constant charge was assigned to give wall charge of 1 e nm−2. Even though the fixed-charge method has some shortcomings including deviations in relaxation dynamics, it has been used successfully for structural studies of IL/electrode interfaces.38,39 The trajectories were analyzed using TRAVIS-1.14.0 and results averaged for all 12 replicas.40
In the application of the algorithm, the Si substrate (SLD = 2.07 × 10−6 Å−2) was modeled as an infinitely thick planar surface. This surface is always covered with a thin layer of silicon oxide (SiO2) whose actual morphology depends on the preparation method.42 The thickness and roughness of this layer were determined through separated NR measurements performed against water/heavy water mixtures. The SiO2 layer (SLD = 3.85 × 10−6 Å−2) was found to be 10 ± 2 Å thick, with a roughness of 4 ± 1 Å. The thickness and roughness of the SiO2 layer obtained from H2O/D2O measurements were fixed for the models of the sample reflectivity measurements.
NR data versus Qz for the 1:2 ChoCl/EG mixtures for different contrasts at the Si wafer are presented in Fig. 3. It can be seen that all curves are smooth, without fringes. The critical edge can only be seen for the D-EG:H-Cho contrast. Applied voltage has no influence on the reflectivity. As a first step, the Si/SiO2/DES model was used for modeling the measured reflectivity, assuming that there is no structured DES layer at the interface. Dashed lines in Fig. 3 clearly show that the values of the reflectivity as simulated for a uniform DES density lie well below the experimental values. In order to introduce a DES layer in the model, one requires its SLD, thickness, and roughness. These values can be obtained through analysis of the MD simulation trajectory. The number distribution of each atom in the mixture is extracted from the MD trajectory as function of its distance from the electrode. After binning this distribution in small distance steps, the number density was converted into the scattering length density SLD by
The difference in SLD values between DES layers and the bulk could originate from changes in composition and/or density. In Fig. 3 it can be seen that an applied voltage does not induce any change in the reflectivity patterns, i.e. there is no apparent change in SLDs. Based on the previous discussion, this means that DES layers are present at the interface even without applied voltage. Further analysis of the MD trajectory can give more information about the structure of this layer and how the applied charge and increased temperature influence this layer.
In Fig. 6 the number density distribution is shown for each component of the DES calculated for its center of mass. It can be seen that the Stern layer43 at a neutral silicon surface primarily consists of ethylene glycol molecules and Cho+ cations with alcohol groups located closer to the electrode. Excess positive charge in the Stern layer favors an anionic chloride enriched layer followed by an EG–Cho rich layer and another EG rich layer before the bulk structure is recovered. As mentioned in the experimental section the charges on Si atoms are set to obtain a charge density of 1 e nm−2 which would correspond to the order of magnitude of the experimental voltage. As expected at charged surfaces, a reorganization of choline and chloride ions occurs. However, this cation–anion reorganization will not noticeably change the SLD profile of the DES since a high number of EG molecules is still present at the interface.
Fig. 6 Number density distribution of DES components at uncharged and charged silicon surfaces calculated for their center of mass. |
In the work of Cooper et al.19 a NR study of ionic liquids confirmed the existence of the interfacial IL layer even at 80 °C. In the present work, MD simulations were performed at 25 °C and 100 °C in order to see if this would be the case with the investigated Si/DES system. No noticeable changes could be observed when the simulations were performed at 100 °C. Since the negative potential is of particular interest for the electrochemical deposition of metals, Fig. 7 shows how the orientation of choline and EG molecules is changed by introducing a negative charge of the silicon surface.
It is apparent that two peaks in the number density profile of choline center of mass correspond to two spatial orientations of choline cations. A part of choline cations lays on the electrode while the other fraction orients perpendicular to the electrode. When the Si surface is uncharged choline cations are layered with the alcohol groups closer to the surface.
A negative charge of the electrode will induce reorientation of the choline cation perpendicular to the surface so that the positively charged part of the cation stays closer to the surface. This is in agreement with the work of Hammond et al.26 where the authors proposed that application of −0.25 V promotes the orientation of Cho charged groups towards the Pt electrode and alcohol groups towards the bulk liquid. It can also be observed that neither positive nor negative potential have a noteworthy influence on EG molecules. In all cases, three distinctive peaks in the number density profiles appear around 0.5, 0.9 and 1.4 nm. It should be emphasized that these changes in molecule orientations cannot be confirmed with the NR technique since the change in the SLD is not large enough to influence the reflectivity curves.
According to the work of Starykevich et al.15 elevated temperature (up to 100 °C) allows electrodeposition of Zn from ethaline DES on an alumina surface, through disruption of the DES layer at the interface. Since the results of this work indicate the existence of the DES layer at the interface with and without applied potential even at 100 °C and no considerable change in its structure, other effects should also be taken into account. For example, higher temperature leads to a viscosity reduction and increased conductivity as well as increased void volume of the DES.44,45 Therefore further experimental and computational investigations are needed to fully understand the phenomena observed in these electrochemical experiments.15 This work shows that NR measurements and MD simulations can be used as complementary techniques to reveal the nanostructure of solid/liquid interface.
This journal is © the Owner Societies 2020 |