DFT and SERS Study of L-Cysteine Adsorption on the Surface of Gold

techniques, X-ray photo-emission and CK-shell X-ray absorption spectra, and the molecules were found to lay flat on the surface with both the Cβ-S bo...
0 downloads 0 Views 2MB Size
Subscriber access provided by Kaohsiung Medical University

C: Surfaces, Interfaces, Porous Materials, and Catalysis

DFT and SERS Study of L-Cysteine Adsorption on the Surface of Gold Nanoparticles Guohua Yao, and Qing Huang J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.8b00949 • Publication Date (Web): 24 May 2018 Downloaded from http://pubs.acs.org on May 24, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

DFT and SERS Study of L-cysteine Adsorption on the Surface of Gold Nanoparticles Guohua Yao1 and Qing Huang1,2* 1

Key Laboratory of High Magnetic Field and Ion Beam Physical Biology, Institute of Technical Biology and Agriculture Engineering, Hefei Institutes of Physical Science, Chinese Academy of Sciences, Hefei 230031, China 2 University of Science & Technology of China, Hefei 230026, China *Corresponding author: [email protected]

Abstract Both surface enhanced Raman spectroscopy (SERS) experiments and density functional theory (DFT) calculations have been carried out to investigate adsorption of cysteine on gold nanoparticle surface. Cysteine is one of the amino acids which have a rotating dihedral angle and several possible adsorption sites, so there could be many possible adsorption conformations. At present, the adsorption conformation of cysteine on the gold surface has not been solved, although previous work had explored the conformations based on DFT method, but no explicit experimental evidence was provided. Nevertheless, the previous work had indicated that the adsorption conformations of some simple rigid molecules on the surface of metal nanoparticles could be identified through SERS experimental data. Therefore, in this work, we attempted to figure out the conformation and adsorption sites of the flexible typical amino acid cysteine on gold nanoparticle surface by combining both DFT and SERS methods. We found that the adsorption of cysteine on the surface gold nanoparticle adopted the most predominant configuration in which cysteine chemically interacted with Au cluster through both S and O sites with cysteine taking the configuration of a gauche position for the protonated amino group and an anti-position for the carboxylate group in the molecule.

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Introduction The strong adsorption of cysteine on gold surface has been widely used in many fields. Proteins or DNA strands can be attached to gold nanostructures through cysteine-gold interactions.1 These modified metallic surfaces may find significant technological applications in chemical-specific detection, sensing and catalysis.2 For example, Gold nanostars coated with cysteine were synthesized and applied in SERS-based copper ions (Cu2+) detection in aqueous media.3 Gold nanoparticle-based colorimetric competition assay was applied in the detection of cysteine based on the interaction between cysteine and gold nanoparticles.4 So, it is necessary to have a better understanding of the molecular properties of the adsorption of cysteine and other amino acids on metal surfaces, and in fact, many theoretical or experimental works have been reported.1-2 Canepa et al. have studied the vapor-phase deposition of L-cysteine (L-cys) on the Au (110) surface by means of synchrotron-based techniques, X-ray photo-emission and CK-shell X-ray absorption spectra, and the molecules were found to lay flat on the surface with both the Cβ-S bond and the carboxylic group almost parallel to the surface.5 Marek Graff et al. reported the adsorption of enantiomeric and racemic cysteine on a silver electrode–SERS sensitive to chirality of the adsorbed molecules.6 In addition, electrochemistry of L-cys on several types of electrode materials has been extensively investigated.2, 7 And the characterization of L-cys adlayers on Au (111) has also been provided using in situ (in an electrochemical environment) scanning tunneling microscopy (STM).8

Figure 1 The Newman projections (left) and the labeling of atoms for the structures of cysteine-gold cluster model (right). Until now, however, the microstructures of the interactions between the amino acids (include cysteine) and the metal nanoparticles are basically still elusive. Since cysteine has a rotable χ1(S-Cβ-Cα-Cγ) torsional angle (Figure 1), three conformers (PH, PC, PN) of cysteine may exist. Considering that the sulfur atom is attached to the surface, we consider that three conformations are possible for the relative orientations of the amino and the carboxylate groups. The Newman projections of these rotamers are presented in Figure 1. But it is not clear which conformation is present when cysteine is adsorbed onto the gold surface. It is known that there is strong chemisorption between cysteine and metal surface through thiol group. Stewart and Fredericks concluded that the ionized carboxylate group and the protonated amino group were also in close contact with the surface.9 But it has not been determined whether the amino group or carboxyl group is chemically adsorbed on the substrate. When cysteine is adsorbed on the gold nanoparticles, SERS signals can be detected, and the SERS

ACS Paragon Plus Environment

Page 2 of 22

Page 3 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

spectra can provide a lot of structural information. The influence on the SERS spectral pattern from the chemical effect can be crucial for identifying the molecular sites involved in the chemical interaction with the substrate. In the past, the adsorption behavior of cysteine in gold colloidal solution was studied according to the ‘‘surface selection rules’’ for analyzing of the SERS spectra.10-11 The SERS signals of cysteine including some SERS bands from sulfur, NH3+ and COOgroups are enhanced; reaching the conclusion that cysteine molecule interacts with the metal surface simultaneously through sulfur, NH3+ and COO- groups.10 Compared with the normal Raman spectra, the wavenumbers of SERS bands are normally shifted when the molecule is adsorbed to the metal surface. So it is still a challenge to assign the SERS bands for determining the involved conformations. With the help of the DFT calculation of the molecule-metal cluster, the chemical interaction between the molecule and substrate can be properly described. The structure of molecule-metal cluster system and the assignment the vibrational modes can be more accurately determined if one can figure out the molecular conformation and the interaction between the molecule and the adsorbing metal. DFT calculations of molecule-cluster models have been effectively used for the simulation of SERS spectra of small molecules adsorbed on the silver and gold surfaces, such as adenine, p-aminothiophenol, benzyl Chloride and methimazole.12-16 In principle, through comparison between the theoretical computation and the experimental SERS spectra, the adsorption of molecule on the surface of the metal can be analyzed and the adsorption conformation can be determined. Actually, there have been some former DFT studies on cysteine adsorbed on metal surfaces,17-18 though with the energy analysis mainly used for estimation of the adsorbing structure. Cordeiro et al. employed the DFT method to study of the adsorption of D‑(L-) cysteine on flat and chiral stepped gold surfaces (Au (321) R, S and Au (111)).17 They concluded that the adsorption occurred with the cysteine bound to the surface through only one contact point (by its sulfur atom). However, there was no comparison between theory and experiment, so it was still difficult to decide the corresponding conformation.17 In this work, density functional theory (DFT) with the M06-2X functional was used for the optimization of the ground state geometries and simulation of the vibrational spectrum of this amino acid. The SERS spectrum with a large silver cluster as a model metallic surface was simulated previously.19 The 18-atom neutral silver cluster was constructed as a closed-shell singlet based on the experimental bulk geometry data with 12 atoms on the surface plane and 6 atoms on a parallel lower layer, representing the Ag [110] surface. But Z-matrix of the parallel was constrained, so it could not reasonably simulate the chemisorption between cysteine and metal. To the best of our knowledge, this is the first research using this molecule-metal cluster model to simulate amino acids (include cysteine) and gold substrate, in which we also compared the DFT calculation with the experimental SERS spectra in order to determine and understand the involved chemical interactions. The focus of the present work was therefore to obtain information on the adsorption behavior of the L-cysteine molecules on gold nanoparticles, and to achieve the vibrational characterization of the adsorbed molecule as well.

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Experimental Section Materials L-cysteine (noted as “cysteine” or “cys” in this work), gold acid chloride trihydrate (HAuCl4·3H2O), sodium citrate, sodium borohydride (NaBH4) were purchased from Sigma-Aldrich. High purity water was used in this work. All the chemicals were used without further purification. Fabrication of gold nanoparticles The aqueous solutions of gold colloids were prepared by two typical methods. (a) In the first method, the gold colloids were prepared by the method described by Creighton et al.20 The aqueous solution of gold colloids were prepared by reduction of 20 mL of 2x10-3 M HAuCl4·3H2O with ice-cold 60 mL of 2.0x 10-3 M NaBH4 solution. The solutions were mixed via vigorous stirring for 1 h. The gold colloids were red-brown and showed two visible absorption bands at 218 nm and 530 nm (Figure S1). (b) In the second method, the gold colloids were prepared by the method described by Lee and Meisel.21 71 μl of 1 M HAuCl was dissolved in 50 mL of water and the solution was brought to boiling. A solution of 5 mL sodium citrate (1 wt %) was added. Boiling was continued for ca. 1 h. The gold colloids were purple and had absorption maxima at approximately 525 nm (Figure S1). SERS measurements For SERS measurements, conjugation of cysteine and gold nanoparticles was formed by mixing 100 µL gold colloid solutions with 100 µL of 2 x 10-4 M cysteine solutions. The pH values of 10-4 M cysteine in these two kinds of gold colloids were about 6.0 in this work. 5 µL of the mixture was deposited into the groove of a quartz slide. All the Raman spectra with resolution ca. 3 cm-1 were recorded using XploRA Raman spectrometer (HORIBA JOBIN YVON). The SERS spectra of the mixture of cysteine and Au colloids ware measured by using a 785 nm laser. The laser power at the sample was ca. 1.2 mW while the exposure time was typically 10 s. The spectra were measured with at least 3 repeats. The SERS spectra of cysteine and gold nanoparticles reduced by (a) sodium borohydride and (b) sodium citrate are shown in Figure S2, respectively, and these two SERS spectral patterns are similar. The SERS spectra which show the main characteristic bands agree with the previous SERS experiments on cysteine in gold colloids (Figure 7c of Ref 11)11, as well as the SERS spectra of the mixture of cysteine and gold colloids (Figure 3a of Ref 10).10 The experimental SERS spectra of cysteine in different gold colloids under different preparation conditions are very similar, indicating that cysteine should be adsorbed on the surface of different gold colloids with the same adsorption conformation. In this work the spectra from gold colloids reduced by sodium borohydride are comparable with the simulated spectra. Characterization of gold nanoparticles after the adsorption of cysteine. To figure out the nanostructure and morphology of gold colloids reduced by NaBH4 after the adsorption of 10-4 M cysteine, transmission electron microscope (TEM) imaging and energy-dispersive X-ray spectroscopy (EDS) mapping were performed on a Talos F200X microscope operated at an accelerating voltage of 200 kV, as shown in Figure 2. The average particle size of the gold colloids was about 20-40 nm, and they aggregated due to the adsorption of cysteine. The spatial distribution of S and N coincides with the spatial distribution of the gold nanoparticles, indicating that cysteine has been adsorbed on the nanoparticles.

ACS Paragon Plus Environment

Page 4 of 22

Page 5 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 2 (a) TEM images of gold colloids after the adsorption of 10-4 M cysteine. (b, c, d) The area distribution of Au (yellow), S (blue) and N (green), respectively. (e) EDS spectra of gold colloids after the adsorption of 10-4 M cysteine.

Computational Details DFT calculations were carried out using Gaussian 09 software.22 The calculations in this work were performed by applying the functionals DFT-D3 correction versions of M06-2X unless noted otherwise.23-24 The M06-2X (hybrid meta-GGA exchange-correlation functional) is a high-nonlocality functional with double the amount of nonlocal exchange (2X) for nonmetals. Corminboeuf et al. had demonstrated that standard density functionals failed to accurately describe interaction energies of charge-transfer complexes not only because of the missing long-range exchange as generally assumed but also as a result of the neglect of weak interactions.25 However, these could be properly described when using well-balanced functionals such as PBE0-dDsC, M06-2X, and LC-BOP-LRD.25 The method of dispersion correction (DFT-D) as an add-on to standard Kohn–Sham density functional theory has been refined regarding higher accuracy, broader range of applicability, and less empiricism.26 The new DFT-D3 versions have been thoroughly tested not only on common organic and noncovalently bound complexes, but also on large, infinite, “heavy”, and metallic systems. DFT-D3 performs high accuracy for dispersion corrections and superior for thermochemical, intramolecular dispersion effects.26 Grimme et al. had made a thorough benchmark of density functional methods for general main group thermo-chemistry, kinetics, and noncovalent interactions, and found that M062X-D3 could be statistically the best of all 23 hybrids-level functional.27 And in our previous work,28 we had showed that the Raman spectra and the geometry parameters of the amino acid molecules (p-tyrosine, m-tyrosine) could be best described by M06-2X-D3 functionals. The basis set used for C, N, S, O, H atoms in this work was the 6-311++G(d,p), which are suitable for amino acids.28-29 The pseudopotential basis sets LANL2DZ was used for the Au atoms. The tight convergence criterion of Gaussian 09 was used in structure optimization, and the ultrafine integration grid was used in the numerical integration of the structure optimization and vibrational frequencies calculation. The geometries of cysteine-Aun complexes were fully optimized without any

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

constraint on the geometry and the optimized structures had no imaginary frequencies. The calculations of the harmonic vibrational wavenumbers and relative Raman activities were carried out at the same level of theory using the same basis sets. The scaling factor for the harmonic vibration frequencies of M06-2X-D3/6–311++G(d,p)/LANL2DZ was 0.983.30-31 Then the calculated activities were converted using the following relationship derived from the basic theory of Raman scattering:32-33  =  −   ⁄ 1 − −ℎ ⁄  , where  is the exciting frequency (in unit of cm-1),  is the vibrational frequency (in unit of cm-1) of the ith normal mode, h, c, and k are fundamental constants, and f is a suitably chosen common normalization factor for all peak intensities. Relative multiples of the spectral intensity ranges are marked on the right Y axis as X1, X2, X5, et al. The Raman spectra shown in this work have a Lorentzian line width of 8 cm−1. Vibrational frequency assignments were made based on the results of the Gaussview program 5.0.8 version,34 and the potential energy distribution (PED) matrix was expressed in terms of a combination of local symmetry and internal coordinates. In the calculation work, we studied the zwitterionic forms of cysteine, because the pH value of 10-4 M cysteine in the gold colloids was about 6.0 in this work. The pKa value of the carboxyl group in cysteine was 1.9, therefore, the carboxyl group was expected to be unprotonated and the amino group is at the protonated form under the conditions used to record SERS spectrum.10, 35 It was also reported that cysteine was usually in zwitterionic form in solution, such as cysteine deposited from an aqueous solution formed a ordered film on the Au (111) surface.36 The implicit solvent models SMD under default setting in Gaussian 09 package were employed for the calculations of both the energy and vibrational properties.37 Besides, the influence of explicit one and more water molecules on the Raman spectra is shown the Figure S3. Since the Raman spectra are only slightly changed, the implicit solvent models are therefore enough accurate for us to simulate and interpret the Raman spectra in this work.

Results and Discussion 1. Conformers of cysteine-gold cluster complexes From the SERS spectra of cysteine on gold colloids (Figure S2), the absence of the band at about 490 cm-1 from the stretching vibration of S-S band of cystine (cys-cys dimer)11 indicates that cysteine was not aggregated. Therefore, we considered that the cysteine monomer was adsorbed on the gold nanoparticles.11 The cysteine molecule is usually adsorbed on gold by its sulfur atom, after cleavage of the S-H bond with its deposition from aqueous solution to the metal surface.17, 38-39 The absence of the S-H stretching vibration band in the SERS spectra (Figure S2) which should appear at about 2570 cm-1 in the ordinary Raman spectrum of cysteine surmises that the thiol was ionized at the metal surface.10-11, 38 Further evidence for the gold-sulfur interaction comes from the formation of the Au-S vibrational band at 267 cm-1 (Figure 3a of Ref. 10) of SERS spectra of cysteine on the gold colloids.10 Besides, Cavalleri et al. have investigated the electronic states of cysteine adsorbed on Au (111) from aqueous solution with synchrotron-based, high resolution X-ray photoemission spectroscopy.39 The analysis of the S 2p core level region in pristine samples prepared with purified cysteine indicates that cysteine is adsorbed onto gold through thiolate.39 The 1H NMR spectrum in aqueous solution also demonstrates that the sulfide bond of cysteine interacts with colloidal gold nanoparticles surface.40

ACS Paragon Plus Environment

Page 6 of 22

Page 7 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

In our simulation of cysteine-Aun models, the initial structures were constructed by cysteine adsorbed to the Aun cluster by the combination of the possible bonds S-Au, -CγO2--Au, and -NH3+-Au. But there are no stable structures that could be optimized in which the -NH3+ group is close to the Aun cluster. This means that the cysteine may not be adsorbed to the gold substrate through -NH3+group under normal conditions. This is also proved by the polarization-dependent SERS study of cysteine on the gold films.41 Since vibrational modes associated with functional groups lying closer to the surface of the metal and having polarizability components perpendicular to the surface will display a stronger polarization effect, hence, this result suggests that cysteine was adsorbed and had a preferred orientation with respect to the surface (i.e., it had bound to the gold by its thiol group and its carboxylic acid group).41 Therefore, the cysteine-Aun models were constructed by the possible bonds S-Au, -CγO2--Au in this work. Previous research results showed that the sulfur atom is connected to Au atom, but whether O atom is connected to Au atom has not been determined. In this work four gold atoms as a cluster simulated the SERS substrate, which had been successfully applied in DFT simulation of SERS of some molecules.13, 42 Since cysteine has three isomers of different rotable χ1(S-Cβ-Cα-Cγ) torsional angle,2, 10, 43 when cysteine is adsorbed on metal surface, the conformers with different torsion angles may also exist.2, 6 Based on the denotation of references 6 and 10, there are three possible rotational conformers of cysteine that may be stabilized by interaction with the metal surface.6 As shown in Figure 1 and 3, PH/PN/PC series refer to the conformers constructed with H atom/NH3+ group/ CγO2- group in the para-position of the S atom, respectively. For instance, the PH/gauche means both CγO2- and NH3+ on the gauche position. The PN/anti(I) is constructed with an anti-position for the NH3+ and a gauche position for the CγO2-. The PC/anti (II) is constructed with a gauche position for the NH3+ and an anti-position for the CγO2-. PH/PN/PC denote the free cysteine conformers which are not interacting with gold. The structures and Raman spectra of the free cysteine conformers (PH, PN, PC) calculated at the same level of theory are shown in Figure S4. In this work, the cases for cysteine adsorbed to the gold cluster through none, one or two oxygen atoms in adding to sulfur atom were also considered. In the case of cysteine adsorbed on gold only through sulfur atom, O0 denoted the case with none oxygen atom adsorbed, and the three conformers PHO0, PNO0, PCO0 are shown in Figure 3A. In the case of sulfur atom and only one oxygen atom was adsorbed, one oxygen atom was closer to the NH3+ group, and the other oxygen was closer to the H atom. As shown in Figure 3B, ON denotes that the Au atom is connect with the oxygen close to the NH3+ group, while OH denotes that the Au atom is connected with the oxygen close to the H atom. Therefore, six different isomers are defined in Figure 3B, as PNON, PNOH, PCON, PCOH,PHON and PHOH. The charges of these 9 complexes are 0 and spin multiplicities are 2, and the optimized structures are shown in Figure 3A and 3B. In the case of two oxygen atoms are adsorbed, as shown in Figure 3C, PHO2 means that the gauche/PH both the CγO2- and NH3+ are on the gauche position and two oxygen atoms are both adsorbed to the gold clusters, also similarly for PNO2 and PCO2. The charges of these PHO2 and PNO2 complexes are +2, and the spin multiplicities are 2. The charge of the PCO2 complex is +3, and the spin multiplicity is 1. The optimized structures are shown in Figure 3C. In the PCO2 complex, due to the distance between the sulfur atom and the CγO2- group, a cluster of at least six gold atoms is required to simulate the metal substrate. The charges, spin multiplicities and some important geometric parameters of these 12 complexes as shown in Figure 3 are listed in Table S1. The 9 complexes in Figure 3A and

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

3B have the same charge, spin multiplicity with the gold cluster, therefore, the relative energies and the binding energies of these 9 complexes are also calculated and listed in Table S2. The complexes PHO0/PNO0/PCO0 which are bound to the gold cluster only by the S-Au bonds, and they are less stable than the complexes which are bound to the gold cluster by the S-Au and O-Au bonds. The PHOH and PHON complexes are most stable in the view of energetics, and the energy of complex PCOH is the third most stable in the 9 complexes, only 2.89 kcal/mol higher than PHOH. The binding energies of these complexes are about 30 to 40 kcal/mol, confirming that cysteine can be chemically adsorbed to the gold surface.

Figure 3 The optimized structures of conformers of cysteine bonding with the gold cluster through (A) sulfur atom (B1, B2) sulfur atom and one oxygen atom (C) sulfur atom and two oxygen atoms 2. Raman spectra of conformers The simulated Raman spectra of 12 complexes and experimental spectra are shown in Figures 4, 5 and 6. The frequencies of the vibrational modes of 12 complexes are also given (Table S3 in Supporting Information), and the frequencies of PH, PN and PC conformers of free cysteine can thus be compared therein. Generally, the wavenumbers of bands are shifted in varying degrees for different complexes due to the influence of the chemical interaction with gold clusters. The strong band of free cysteine conformers at about 2690 cm-1 is from the stretching vibration of S-H bond (thiol group), and the weak band at about 980 cm-1 is from the bending vibration of thiol group. In the Raman spectra of the cysteine-Aun complexes, these two bands disappear due to the formation of S-Au bonds instead. Besides, the bands in the region of 600 – 1000 cm-1 are

ACS Paragon Plus Environment

Page 8 of 22

Page 9 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

shifted considerably when cyseine interacts with gold. Because these bands of free cysteine contain more or less a part of stretching vibration of C-S bond or a part of bending vibration of C-S-H, although the bands are shifted with interaction with gold, some characteristic bands of free cysteine conformers can still be identified, while the frequency shift or difference can help to determine the rotamers. For instance, PH serial conformers (PH, PHO0, PHOH, PHON, PHO2) have a medium intensity band at about 460 - 480 cm-1, but none for PC and PN serial conformers. For the PN serial conformers (PN, PNO0, PNOH, PNON, PNO2), the bands from stretching vibration of Cβ-S bond are all located at about 750 - 770 cm-1, while the bands at located at about 680 - 710 cm-1 for PC and PH serial conformers. In the experimental spectra, the band at 673 cm-1 is rather intensive, while only the spectra of PHON, PHOH, PCOH and PHO0 can reproduce the intensity and wavenumber of this band. We found that the spectra of PNOH, PCON, PNO0, PCO0, PHO2, and PNO2 cannot reproduce intensity of this band, and the wavenumbers of the band of PNON and PCO2 are not very consistent with the experiment. There are two strong bands at about 900 cm-1 and 1050 cm-1, and some weaker bands in the regions between 700 cm-1 to 1250 cm-1. Only the results of PCOH, PCON, PNOH and PCO2 can reproduce the wavenumbers and intensities of these two strong bands. Only PCOH has three strong bands at about 1292 cm-1, 1339 cm-1 and 1391 cm-1 in the regions between 1250 cm-1 to 1500 cm-1, but the simulated wavenumbers are overestimated. Owing to the universally used frequency scaling factor, the wavenumber values of bands below 1250 cm-1 can be corrected, but the corrected wavenumbers of bands over 1250 cm-1 are usually a little overestimated when compared with the experimental values.42, 44-45 With comparison and analysis of the theoretical and experimental data, the Raman spectra of PCOH complex (Figure 5) can mostly reproduce the measured SERS spectra for gold colloids (Figure 5). We have therefore figured out the most dominant adsorption configuration for cysteine attached on the surfaces of gold nanoparticles, in which cysteine interacts with Au cluster through both S and O sites, and cysteine forms a gauche position for the NH3+ and an anti-position for the CγO2-. On the other hand, the effect of the level of theory was also checked, and the other two well evaluated functionals,27 namely, DFT-D3 correction versions of B3LYP and DFT-D correction of ωB97X-D, were employed to simulated the Raman spectra of a series of complexes (PHON, PHOH, PNON, PNOH, PCON and PCOH) and the results are shown in Figure S6 and S7, respectively. Similarly, the Raman spectra of PCOH complex (Figures S6 and S7) can mostly reproduce the measured SERS spectra for gold colloids. Their frequencies of the vibrational modes are also listed in the Table S4 in Supporting Information. The standard deviations of the same vibrational mode calculated by these three functionals are rather small (5 cm-1 on the average), indicating that the Raman spectra simulated by M06-2X-D3 in this work are credible.

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 4 The experimental SERS spectra of 10-4 M cysteine in gold colloids (blue line), and the DFT-simulated Raman spectra of the three conformers of cysteine bonding with the gold cluster through sulfur atom (black lines). Multi.: spin multiplicity.

ACS Paragon Plus Environment

Page 10 of 22

Page 11 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 5 The experimental SERS spectra of 10-4 M cysteine in gold colloids (blue line), and the DFT-simulated Raman spectra of the six conformers of cysteine bonding with the gold cluster through sulfur atom and one oxygen atom (black lines).

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 6 The experimental SERS spectra of 10-4 M cysteine in gold colloids (blue line), and the DFT-simulated Raman spectra of the three conformers of cysteine bonding with the gold cluster through sulfur atom and two oxygen atoms (black lines). In order to verify the reliability of adsorption on PCOH conformers, cysteine absorbed to six different gold clusters were simulated and the result are shown in Figure 7. The Raman spectra, the charges and the spin multiplicities of the whole molecule-cluster systems and the optimized structures are shown in the figure, and the vibrational frequencies are also listed in Table S5. The XPS measurements by Pagliai et al. showed that in the silver and gold colloids, both silver and gold surfaces presented a sizable amount of metal atoms with oxidation number +1, which allows the metal cluster take zero, one, or even more positive charges.12 From comparison of Figure 7a with Figure 7b, Figure 7c with Figure 7d, we see that the difference of the charge mainly influence the structure of the metal clusters, but do not obviously influence the chemical interaction between cysteine and metal cluster. When cysteine molecule is absorbed to significantly different clusters (Figure 7a to Figure 7f), the structures of the cysteine would have only slight differences. And the differences between the Cγ-O-Au and Cβ-S-Au angles of these structures do not exceed 10 degrees, respectively. The difference between the bond lengths of O-Au and S-Au bonds does not exceed 0.1 Å. Therefore, it’s reasonable that the simulated Raman spectra of these conformers have the similar the Raman shifts and relative intensities. The standard deviations of the same vibrational modes of the complexes with different gold clusters are less than 10 cm-1 (4 cm-1 on the average), as shown in Table S5. Therefore, the structures of

ACS Paragon Plus Environment

Page 12 of 22

Page 13 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

gold clusters may have little influence on the Raman spectra when cysteine interacts with the substrate through the same sites. Thus the metal cluster with four atoms may be enough for the simulation of the Raman spectra of cysteine on the gold surface. This agrees with Wu et al that the size of the silver cluster does not influence significantly the Raman spectra of adsorbed adenine.13 Since gold atoms on the gold colloid surfaces are not periodically regular, the spatial structure of gold atoms on the surface is different, similar to the different spatial structure of the simulated clusters. These can induce small difference of structural parameters of cysteine on the gold surfaces and the small Raman shifts of bands, which may explain that the widths of experimental SERS bands are rather wide.

Figure 7 The experimental SERS spectra of 10-4 M cysteine in gold colloids (blue line), and the DFT-simulated Raman spectra of the six PCOH conformers with different gold clusters (black lines).

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

3. Vibrational assignments A vibrational normal-mode analysis of L-cysteine has been carried out.46 According to the calculations, Raman frequencies can be used to identify specific rotamers of cysteine in solution.46 Unfortunately, the conformation-sensitive modes are the Cβ-S stretch and the S-H stretch. The oxidative adsorption of cysteine on gold breaks the S-H bond, and the Cβ-S mode is also strongly perturbed by the presence of the surface. The adsorption behaviors of cysteine in gold colloidal solution was resulted from the ‘‘surface selection rules’’ and may be not credible enough.10-11 Many previous works 10-11 generally surmised the vibrational modes of cysteine adsorbed on the gold surface without consideration the chemical interaction of metal clusters with cysteine. As seen in Figure 4, 5, 6, Figure S4 and Table S3, the wavenumbers of the bands are shifted in varying degrees for different complexes due to the influence of the chemical interaction with gold clusters. For instance, the wavenumbers of CγO2- symmetric stretching vibration of free cysteine (PH, PN, PC in Figure S4) are 1408 cm-1, 1419 cm-1 and 1412 cm-1, respectively. These bands shift only a small amount when the cysteine interacts with gold cluster not through the oxygen atom, which are 1417 cm-1 for PHO0, 1413 cm-1 for PNO0 1417 cm-1 for PCO0 while the intensities of these bands are weak. However, when the cysteine interacts with gold clusters through one or two oxygen atom, these bands are moved to 1411 cm-1 for PHON, 1420 cm-1 for PHOH, 1402 cm-1 for PNON, 1420 cm-1 for PNOH, 1370 cm-1 for PCON, 1368 cm-1 for PCOH (Figure 5), 1436 cm-1 for PHO2, 1432 cm-1 for PHO2, 1428 cm-1 for PCO2 (Figure 6), respectively. To be noted, the intensities of these bands are all rather strong. The vibrational modes of the cysteine-gold cluster complex (PCOH) are assigned in Table 1. The strongest Raman band observed at 673 cm−1 is assigned to stretching vibration of Cβ-S bond, combining with the scissoring vibration of CγO2-, same as almost all previous literatures. For the PN serial conformers (PN, PNO0, PNOH, PNON, PNO2), the simulated SERS bands of stretching vibration of Cβ-S bond are all located at about 760 cm-1. For PC serial conformers and PH serial conformers, the bands are at about 680 to 700 cm-1 region, meaning that the major χ1 rotamers may be not only in PH type,10-11 but also in PC type. Besides, in the simulated Raman spectra of free cysteine (Figure S4), this band of both PH and PC rotamers is also at about 685 cm-1. Therefore, in Ref. 10 and 11, it was not rigorous to judge that cysteine is adsorbed to the surface of gold in PH serial conformers only according to the information of this band.10-11 So we adopted additional ways, and in the simulated Raman spectra of PH serial conformers (PHO0, PHOH, PHON, PHO2 and PH), they all have a medium intensity band at about 465 cm-1 which is assigned to the out-of-plane bending of CγO2-, combining with the bending of Cα-Cβ-N. This band cannot be observed on the experimental SERS spectrum, which also indicates that cysteine is not adsorbed as the PH series conformers. The 1339 cm-1 band is mainly assigned to the CγO2- symmetric stretching vibration combining with the bending vibration of Cα-H, and the wavenumber of this band is reduced from the normal Raman band of CγO2- symmetric stretching at about 1400 cm-1, maybe caused by the chemical interaction between the oxygen atom and gold substrate.11 The band at 900 cm-1 is mainly assigned to the Cα–Cγ stretching vibration combining with the rocking the NH3+ group and scissoring of the CγO2- group.10-11, 28, 47 The enhanced of this band may also be caused by the chemical interaction between the oxygen atom and gold substrate. In the simulated Raman

ACS Paragon Plus Environment

Page 14 of 22

Page 15 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

spectra of PN serial conformers (PNON, PNOH, PNO0, PNO2) and PH serial conformers (PHON, PHOH, PHO0, PHO2), the bands from Cα–Cγ stretching vibration are located at about 850 cm-1. This means that the major χ1 rotamers may be not in PN and PH type. In the previous work,10-11 based on the observation that these bands involving the vibrational motions around Cα−N were enhanced, researchers deduced that the cysteine interaction with Au NPs might preferentially occur also through NH3+ terminal groups. But the bands at the 1500-1600 cm-1 region are also from the terminal groups. If the NH3+ terminal is close to the Au clusters (such as PCO2, PHO0), the simulated Raman bands at the 1500-1600 cm-1 region are very strong. But in the experimental SERS spectra, the bands at the 1500-1600 cm-1 region are very weak. This means that the cysteine interaction with Au NPs may not be through NH3+ terminal groups. The band at 1050 cm-1 is assigned to twisting vibration of CβH2 and the stretching vibration of Cα−N.11 The 1292 cm-1 band is assigned to the bending of Cα-H and wagging of CβH2. The 1391 cm-1 band is mainly assigned to scissoring of CβH2. The intensities of these bands in the simulated Raman spectra of different conformers are all not low. And the simulated wavenumbers from the same vibrational mode of different conformers are close to each other. For example, the simulated bands from twisting vibration of CβH2 of all conformers (PN, PH, PC serial conformers) are all at 1060-1070 cm-1 region. This may be due to the fact that the vibration of CβH2 group is not insensitive to χ1 rotating isomers and the CβH2 group does not directly contact with the substrate.

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Table 1. Comparison between the SERS spectra and theoretically calculated spectra of cysteine in Au colloids, the assignments and PED of vibrational modes of PCOH complex at M06-2X-D3/6-311++G(d, p) (C, N, O, S and H) level/LANL2DZ (Au) level. exp

a

1574

exp

b

1571

exp

11

--

exp

10

1599

cal

PED(%)

1584

asym str CγO2- (50), bend NH3+ (25)

1576

bend NH3+ (78), asym str CγO2- (10)

1560

bend NH3+ (72), asym str CγO2- (14)

1448

sym bend NH3+ (81)

1444

bend Cα-H(35), sym str CγO2- (15), sicss CβH2(14)

1391

1391

1392

--

1436

sicss CβH2(46)

1339

1339

1343

1337

1368

sym str CγO2- (37), bend Cα-H (34)

1292

1290

1295

1295

1313

bend Cα-H (27), wag CβH2(26)

1232

1240

--

--

1254

twist CβH2(52)

1228

wag CβH2(31), bend Cα-H (18)

1128

str Cα-Cβ(36), rock NH3+(15), str Cα-N(12)

1072

rock NH3+(40)

1125

1125

--

--

1050

1050

1048

1058

1063

twist CβH2(37), str Cα-N(24)

949

953

--

--

949

str Cα-Cβ(34), rock NH3+(26)

900

899

900

907

889

str Cα-Cγ(21), rockNH3+(14), sciss CγO2-(12)

832

838

--

832

834

wag CβH2(38), str Cα-N(23)

789

790

--

--

807

bend CγCαCβ(36), str Cα-S(21), str Cα-Cβ(11)

684

str Cβ-S(51), sciss CγO2-(25)

679

str Cβ-S(30), sciss CγO2-(23), bend-out CγO2-(19)

535

bend CαCγO(37),bend-out CγCαCβ(19), str Cα-N(15)

673

671

527

524

665

672

The wavenumbers of experimental SERS spectra of cysteine on gold nanoparticles reduced by (a) sodium borohydride and (b) sodium citrate in this work, (11) sodium citrate in Ref. 11, (10) sodium citrate in Ref. 10. Bend, bending; bend-out, out of plane bending; rock, rocking; sciss, scissoring; str, stretching; twist, twisting; wag, wagging; asym, asymmetry; sym, symmetry.

ACS Paragon Plus Environment

Page 16 of 22

Page 17 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Conclusion The paper reports the study of the adsorption behavior of L-cysteine on gold nanoparticles through the methods of experiment SERS and theoretical DFT simulation. DFT method was used to construct a series of conformations of rotational isomers and different adsorption sites so as to simulate the corresponding Raman spectra. Through comparison of spectral profiles of the experimental with the calculated Raman spectra, we found that the adsorption of cysteine on the surfaces gold nanoparticles adopted the most predominant adsorption configuration in which cysteine chemically interacts with Au cluster through both S and O sites, and cysteine forms the configuration with a gauche position for the NH3+ and an anti-position for the CγO2-. Furthermore, we also demonstrated that the method applied in this work may also be useful to study the adsorption of other amino acids on the surface of metal nanoparticles.

Supporting Information Description Figure S1. UV-visible absorption spectra of prepared gold colloids: (a) reduced by sodium borohydride; (b) reduced by sodium citrate. Figure S2. SERS spectra of cysteine adsorbed on the prepared gold colloids: (a) reduced by sodium borohydride; (b) reduced by sodium citrate. Figure S3. DFT-simulated Raman spectra of the PCOH conformers with (a) three, (b) two, (c) one and (d) none water molecules. Figure S4. Optimized structures and simulated Raman spectra of the PH, PN and PC conformers of free cysteine which are not interacting with gold at M06-2X-D3/6-311++G(d, p) (C, N, O, S and H) level. Figure S5. blue line: the experimental SERS spectra of 10-4 M cysteine in gold colloids. Black lines: the DFT-simulated Raman spectra of the six conformers of cysteine bonding with the gold cluster through sulfur atom and one oxygen atom at B3LYP-D3/6-311++G(d, p) (C, N, O, S and H) level/LANL2DZ (Au) level. The scaling factor for the vibrational frequencies is 1.0008. Figure S6. blue line: the experimental SERS spectra of 10-4 M cysteine in gold colloids. Black lines: the DFT-simulated Raman spectra of the six conformers of cysteine bonding with the gold cluster through sulfur atom and one oxygen atom at ωB97X-D/6-311++G(d, p) (C, N, O, S and H) level/LANL2DZ (Au) level. The scaling factor for the vibrational frequencies is 0.9867. Table S1. Charges, spin multiplicities and some important structural parameters of the 12 different complexes optimized at M06-2X-D3/6-311++G(d, p) (C, N, O, S and H) level/LANL2DZ (Au) level. The bond lengths (unit in angstroms, Å), the angles and dihedral angle (unit in °). Table S2. The relative energies (in kcal/mol) of different cys-Au4 complexes and the binding energy (in kcal/mol) of different cys-Au4 complexes calculated at M06-2X-D3/6-311++G(d, p) (C, N, O, S and H) level/LANL2DZ (Au) level. The Basis Set Superposition Error (BSSE) has been taken into account in the calculations of the binding energies. Table S3. The vibrational frequencies (cm-1) of the calculated Raman spectra of 3 cysteine conformers (PH, PN and PC) which are calculated at M06-2X-D3/6-311++G(d, p) (C, N, O, S and H) level and 12 cysteine-Aun complexes which are calculated at M06-2X-D3/6-311++G(d, p) (C, N, O, S and H) level/LANL2DZ (Au) level. Table S4a&S4b. The vibrational frequencies (cm-1) of the calculated Raman spectra of complexes (PHON, PHOH, PNON, PNOH, PCON and PCOH) which are calculated using the M062X-D3, B3LYP-D3,

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ωB97X-D functionals, and the basis sets are 6-311++G(d,p) for (C, N, S, O and H)/LANL2DZ for (Au). Table S5. The vibrational frequencies (cm-1) of the calculated Raman spectra of PCOH with different Aun clusters which are calculated at the M062X-D3/6-311++G(d,p) (C, N, S, O and H)/LANL2DZ(Au) level.

ACS Paragon Plus Environment

Page 18 of 22

Page 19 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Acknowledgements This work was partly supported by the National natural Science Foundation of China (11635013, 11475217, 11775272) and the Anhui Provincial Natural Science Foundation (1508085QB44). The computation work was supported by USTC-SCC Supercomputer Centers.

References (1) Chang, J. B.; Kim, Y. H.; Thompson, E.; No, Y. H.; Kim, N. H.; Arrieta, J.; Manfrinato, V. R.; Keating, A. E.; Berggren, K. K., The orientations of large aspect-ratio coiled-coil proteins attached to gold nanostructures. Small 2016, 12 (11), 1498-1505. (2) Brolo, A. G.; Germain, P.; Hager, G., Investigation of the adsorption of L-cysteine on a polycrystalline silver electrode by surface-enhanced Raman scattering (SERS) and surface-enhanced second harmonic generation (SESHG). J. Phys. Chem. B 2002, 106 (23), 5982-5987. (3) Ndokoye, P.; Ke, J.; Liu, J.; Zhao, Q. D.; Li, X. Y., L-cysteine-modified gold nanostars for SERS-based copper Ions detection in aqueous media. Langmuir 2014, 30 (44), 13491-13497. (4) You, J.; Hu, H. Z.; Zhou, J. P.; Zhang, L. N.; Zhang, Y. P.; Kondo, T., Novel cellulose polyampholyte-gold nanoparticle-based colorimetric competition assay for the detection of cysteine and mercury(II). Langmuir 2013, 29 (16), 5085-5092. (5) Cossaro, A.; Terreni, S.; Cavalleri, O.; Prato, M.; Cvetko, D.; Morgante, A.; Floreano, L.; Canepa, M., Electronic and geometric characterization of the L-cysteine paired-row phase on Au(110). Langmuir 2006, 22 (26), 11193-11198. (6) Graff, M.; Bukowska, J., Adsorption of enantiomeric and racemic cysteine on a silver electrode SERS sensitivity to chirality of adsorbed molecules. J. Phys. Chem. B 2005, 109 (19), 9567-9574. (7) Zhang, L. Y.; Yuan, R.; Chai, Y. Q.; Li, X. L., Investigation of the electrochemical and electrocatalytic behavior of positively charged gold nanoparticle and L-cysteine film on an Au electrode. Anal. Chim. Acta 2007, 596 (1), 99-105. (8) Zhang, J. D.; Chi, Q. J.; Nielsen, J. U.; Friis, E. P.; Andersen, J. E. T.; Ulstrup, J., Two-dimensional cysteine and cystine cluster networks on Au(111) disclosed by voltammetry and in situ scanning tunneling microscopy. Langmuir 2000, 16 (18), 7229-7237. (9) Stewart, S.; Fredericks, P. M., Surface-enhanced Raman spectroscopy of amino acids adsorbed on an electrochemically prepared silver surface. Spectrochim Acta A 1999, 55 (7-8), 1641-1660. (10) Jing, C. Y.; Fang, Y., Experimental (SERS) and theoretical (DFT) studies on the adsorption behaviors of L-cysteine on gold/silver nanoparticles. Chem. Phys. 2007, 332 (1), 27-32. (11) López-Tobar, E.; Hernández, B.; Ghomi, M.; Sanchez-Cortes, S., Stability of the disulfide bond in cystine adsorbed on silver and gold nanoparticles as evidenced by SERS data. The Journal of Physical Chemistry C 2013, 117 (3), 1531-1537. (12) Pagliai, M.; Caporali, S.; Muniz-Miranda, M.; Pratesi, G.; Schettino, V., SERS, XPS, and DFT study of adenine adsorption on silver and gold Surfaces. The Journal of Physical Chemistry Letters 2012, 3 (2), 242-245. (13) Huang, R.; Yang, H. T.; Cui, L.; Wu, D. Y.; Ren, B.; Tian, Z. Q., Structural and charge sensitivity of surface-enhanced Raman spectroscopy of adenine on silver surface: a quantum chemical study. J. Phys. Chem. C 2013, 117 (45), 23730-23737. (14) Wu, D. Y.; Zhao, L. B.; Liu, X. M.; Huang, R.; Huang, Y. F.; Ren, B.; Tian, Z. Q., Photon-driven charge transfer and photocatalysis of p-aminothiophenol in metal nanogaps: a DFT study of SERS. Chem. Commun. 2011, 47 (9), 2520-2522.

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(15) Wang, A.; Huang, Y. F.; Sur, U. K.; Wu, D. Y.; Ren, B.; Rondinini, S.; Amatore, C.; Tian, Z. Q., In situ identification of intermediates of benzyl chloride reduction at a silver electrode by SERS coupled with DFT calculations. J. Am. Chem. Soc. 2010, 132 (28), 9534-9536. (16) Muniz-Miranda, M.; Muniz-Miranda, F.; Pedone, A., Raman and DFT study of methimazole chemisorbed on gold colloidal nanoparticles. Phys. Chem. Chem. Phys. 2016, 18 (8), 5974-80. (17) Fajin, J. L. C.; Gomes, J. R. B.; Cordeiro, M. N. D. S., DFT study of the adsorption of D(L-)cysteine on flat and chiral stepped gold surfaces. Langmuir 2013, 29 (28), 8856-8864. (18) Di Felice, R.; Selloni, A., Adsorption modes of cysteine on Au(111): thiolate, amino-thiolate, disulfide. J. Chem. Phys. 2004, 120 (10), 4906-4914. (19) Fleming, G. D.; Finnerty, J. J.; Campos-Vallette, M.; Celis, F.; Aliaga, A. E.; Fredes, C.; Koch, R., Experimental and theoretical Raman and surface-enhanced Raman scattering study of cysteine. J. Raman Spectrosc. 2009, 40 (6), 632-638. (20) Creighton, J. A.; Blatchford, C. G.; Albrecht, M. G., Plasma resonance enhancement of Raman-scattering by pyridine adsorbed on silver or gold sol particles of size comparable to the excitation wavelength. J Chem Soc Farad T 2 1979, 75, 790-798. (21) Lee, P. C.; Meisel, D., Adsorption and surface-enhanced Raman of dyes on silver and gold sols. J. Phys. Chem. 1982, 86 (17), 3391-3395. (22) Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. A., et al. Gaussian 09, D. 01, Gaussian, Inc.: Wallingford CT, 2009. (23) Zhao, Y.; Truhlar, D. G., The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and transition elements: two new functionals and systematic testing of four M06-class functionals and 12 other functionals. Theor. Chem. Acc. 2008, 120 (1-3), 215-241. (24) Zhao, Y.; Truhlar, D. G., A prototype for graphene material simulation: structures and interaction potentials of coronene dimers. J. Phys. Chem. C 2008, 112 (11), 4061-4067. (25) Steinmann, S. N.; Piemontesi, C.; Delacht, A.; Corminboeuf, C., Why are the interaction energies of charge-transfer complexes challenging for DFT? J. Chem. Theory Comput. 2012, 8 (5), 1629-1640. (26) Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H., A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. J. Chem. Phys. 2010, 132 (15), 154104. (27) Goerigk, L.; Grimme, S., A thorough benchmark of density functional methods for general main group thermochemistry, kinetics, and noncovalent interactions. Phys. Chem. Chem. Phys. 2011, 13 (14), 6670-6688. (28) Yao, G. H.; Zhang, J. J.; Huang, Q., Conformational and vibrational analyses of meta-tyrosine: an experimental and theoretical study. Spectrochim Acta A 2015, 151, 111-123. (29) Bachrach, S. M., Microsolvation of glycine: a DFT study. J. Phys. Chem. A 2008, 112 (16), 3722-3730. (30) Alecu, I. M.; Zheng, J. J.; Zhao, Y.; Truhlar, D. G., Computational thermochemistry: scale factor databases and scale factors for vibrational frequencies obtained from electronic model chemistries. J. Chem. Theory Comput. 2010, 6 (9), 2872-2887. (31) Kesharwani, M. K.; Brauer, B.; Martin, J. M. L., Frequency and zero-point vibrational energy scale factors for double-hybrid density functionals (and other selected methods): can anharmonic force fields be avoided? J. Phys. Chem. A 2015, 119 (9), 1701-1714.

ACS Paragon Plus Environment

Page 20 of 22

Page 21 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(32) Muniz-Miranda, M.; Gellini, C.; Pagliai, M.; Innocenti, M.; Salvi, P. R.; Schettino, V., SERS and computational studies on microRNA chains adsorbed on silver surfaces. J. Phys. Chem. C 2010, 114 (32), 13730-13735. (33) Krishnakumar, V.; Keresztury, G.; Sundius, T.; Seshadri, S., Density functional theory study of vibrational spectra and assignment of fundamental vibrational modes of 1-methyl-4-piperidone. Spectrochim Acta A 2007, 68 (3), 845-850. (34) Dennington, R.; Keith, T.; Millam, J. GaussView 5, Semichem Inc.: Shawnee Mission, KS, 2009. (35) Mocanu, A.; Cernica, I.; Tomoaia, G.; Bobos, L. D.; Horovitz, O.; Tomoaia-Cotisel, M., Self-assembly characteristics of gold nanoparticles in the presence of cysteine. Colloid Surface A 2009, 338 (1-3), 93-101. (36) Dodero, G.; De Michieli, L.; Cavalleri, O.; Rolandi, R.; Oliveri, L.; Dacca, A.; Parodi, R., L-cysteine chemisorption on gold: an XPS and STM study. Colloid Surface A 2000, 175 (1-2), 121-128. (37) Marenich, A. V.; Cramer, C. J.; Truhlar, D. G., Universal solvation model based on solute electron density and on a continuum model of the solvent defined by the bulk dielectric constant and atomic surface tensions. J. Phys. Chem. B 2009, 113 (18), 6378-6396. (38) Aryal, S.; Remant, B. K. C.; Dharmaraj, N.; Bhattarai, N.; Kim, C. H.; Kim, H. Y., Spectroscopic identification of S-Au interaction in cysteine capped gold nanoparticles. Spectrochim Acta A 2006, 63 (1), 160-163. (39) Cavalleri, O.; Gonella, G.; Terreni, S.; Vignolo, M.; Floreano, L.; Morgante, A.; Canepa, M.; Rolandi, R., High resolution X-ray photoelectron spectroscopy of L-cysteine self-assembled films. Phys. Chem. Chem. Phys. 2004, 6 (15), 4042-4046. (40) Aryal, S.; Bahadur, K. C. R.; Bhattarai, N.; Kim, C. K.; Kim, H. Y., Study of electrolyte induced aggregation of gold nanoparticles capped by amino acids. J. Colloid Interface Sci. 2006, 299 (1), 189-195. (41) Anema, J. R.; Brolo, A. G., Use of polarization-dependent SERS from scratched gold films to monitor the electrochemically-driven desorption and readsorption of cysteine. J. Electroanal. Chem. 2010, 649 (1-2), 159-163. (42) Yao, G. H.; Zhai, Z. M.; Zhong, J.; Huang, Q., DFT and SERS study of N-15 full-labeled adenine adsorption on silver and gold surfaces. J. Phys. Chem. C 2017, 121 (18), 9869-9878. (43) Pawlukojc, A.; Leciejewicz, J.; Ramirez-Cuesta, A. J.; Nowicka-Scheibe, J., L-cysteine: neutron spectroscopy, Raman, IR and ab initio study. Spectrochim. Acta. A Mol. Biomol. Spectrosc. 2005, 61 (11-12), 2474-81. (44) Giese, B.; McNaughton, D., Surface-enhanced Raman spectroscopic and density functional theory study of adenine adsorption to silver surfaces. J. Phys. Chem. B 2002, 106 (1), 101-112. (45) Laury, M. L.; Carlson, M. J.; Wilson, A. K., Vibrational frequency scale factors for density functional theory and the polarization consistent basis sets. J. Comput. Chem. 2012, 33 (30), 2380-2387. (46) Li, H. M.; Wurrey, C. J.; Thomas, G. J., Structural studies of viruses by laser Raman-spectroscopy .36. Cysteine conformation and sulfhydryl Interactions in proteins and viruses .2. normal coordinate analysis of the cysteine side-chain In model compounds. J. Am. Chem. Soc. 1992, 114 (19), 7463-7469. (47) Lee, H. I.; Suh, S. W.; Kim, M. S., Raman-spectroscopy of L-tryptophan-containing peptides adsorbed on a silver surface. J. Raman Spectrosc. 1988, 19 (7), 491-495.

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

TOC GRAPHIC

ACS Paragon Plus Environment

Page 22 of 22