The Structure and UV Spectroscopy of Benzene -‐ Water ( Bz -‐ W 6 ) Clusters using Time -‐ Dependent Density Functional Theory

DFT and MP2 calculations are performed to obtain optimized ground state geometries and binding energies of the cage and the prism conformers of water W6 clusters and Bz-W6 clusters using the aug-cc-pVDZ basis set. The cage conformer of Bz-W6 system is found to be more stable than prism conformer for all range of DFT functionals and MP2. Time dependent-DFT is then used to study UV spectroscopy of Bz, water W6 clusters and Bz-W6 clusters at both the MP2 and wB97XD optimized ground state geometries using the B3LYP, CAM-B3LYP and M06-2X functionals with 6-31++G(d,p) and aug-cc-pVTZ basis sets. Our results predict minor differences in the UV spectroscopy of cage and prism conformers W6 and Bz-W6 clusters that may be observable with high-resolution spectroscopy. The M06-2X and CAM-B3LYP functionals perform consistently with each other. Benzene-mediated excitations of the water W6 cluster towards longer wavelengths above 170 nm are noticed in both the cage and prism geometries of Bz-W6. Benzene is found to be influenced after interacting with the cage and prism W6 geometries, and is seen to undergo a red shift in the main π→π* electronic transition, in which the degeneracy is slightly broken. Charge transfer (CT) states and diffuse Rydberg-type states are also found to play an important role in the spectroscopy of such systems.


Introduction
Water clusters are of fundamental importance in many areas of chemistry, and are also a useful model system to understand the relation of the properties of the gas phase clusters to the condensed phase, ice, and liquid water. The structures and properties of water W n clusters have been studied extensively both theoretically and experimentally because of their importance in many physical, chemical and biological fields, such as in the understanding of cloud and ice formation, bio-chemical processes, etc. [1][2][3][4][5][6][7][8][9][10][11] The interactions of aromatic molecules with water solvents, and the role of these weak interactions to determine the physical and chemical properties of such systems is an area of particular interest for many theoretical and experi-mental studies. In this context, the simplest system of this type, benzene-water complexes, considering benzene as a prototype of intermolecular interactions involving aromatic systems, has been studied extensively. 1,[12][13][14][15][16][17][18][19][20][21][22][23] The study of complexes of benzene (Bz) with water clusters is also of considerable astrophysical importance and can replicate polycyclic aromatics-ice systems to gain better understanding on the interactions of polycyclic aromatic hydrocarbons (PAHs), which may account for up to 20% of galactic carbon, 24,25 with water being the most abundant molecule in icy grain mantles. 26 Benzene has been detected in the proto-planetary nebula CRL 618 27,28 and is amongst the list of known interstellar molecules. Benzene may be thought of as a prototypical PAH compound that plays an important role as an intermediate in the formation of PAHs from acetylene, [28][29][30][31] and water clusters can be a good representation of interstellar ice surfaces. It has been observed that UV irradiation of water ice containing PAHs may play an important role in the formation of complex organic species such as alcohols, quinones and ethers. 32 Photo-desorption is an important process to account for the high gas phase abundances of water under astrophysically relevant conditions. 33,34 Benzene and its derivatives, also being environmental pollutants 35,36 and the constituents of many organic products, e.g. solvents, perfumes, etc. can be efficiently excited electronically by UV radiation present in the † Electronic supplementary information (ESI) available: UV spectra ( Fig. S1 and S2) and important electronic transitions corresponding to higher oscillator strengths for wB97XD optimized Bz-W 6 prism and cage shaped clusters (Tables S1 and S2); comparison of UV spectra obtained from TD-DFT calculations on MP2 and wB97XD optimized ground state geometries of prism and cage shaped Bz-W 6 cluster ( Fig. S3 and S4); electronic transitions obtained using B3LYP, CAM-B3LYP and M06-2X hybrid functionals with aug-cc-pVTZ basis set on MP2 optimized Bz-W 6 prism and cage shaped clusters (Tables S3 and S4). See DOI: 10.1039/c4pp00211c environment. Therefore, the computational study of the ground and excited states of the complexes of benzene with water clusters is of possible environmental as well as astrophysical relevance.
Numerous theoretical and experimental studies have been performed on the ground state properties such as binding energies and IR spectra of Bz-W n clusters, and also on noncovalent interactions such as hydrogen bonding interactions that dominate such systems. 1,6,11,[17][18][19][20][37][38][39][40][41] Zwier and their coworkers 15,[18][19][20]22,23,39 have undertaken detailed experimental studies on Bz-W n clusters using resonant ion-dip spectroscopy (RIDIRS) and resonant two-photon ionization (R2PI) techniques. It was observed that the water O-H stretch spectra depends on the size of cluster and is sensitive to the number, strength and the type of hydrogen bonds in which it participates while C-H stretch in benzene was found insensitive to the size of water cluster with n ≤ 7. 20 For n = 3-5, the water cluster was considered to bind to benzene primarily through single water molecule on or near benzene six-fold axis and the O-H spectra was consist of free O-H stretches, π H-bonded O-H stretches and the single donor O-H stretches. The significant spectral signature was observed at n ≥ 6 in O-H stretch spectra, showing new transitions that are associated with the double-donor O-H stretches, which exist, in the more compact large Bz-W n clusters. Experimental evidence predicted the change in shape of cluster from cyclic to non-cyclic at n = 6 such that Bz-W 6 cluster was found likely to exist as a cage shaped arrangement of six water molecules giving the lowest energy structure.
A detailed computational study on geometries, binding energies and infrared (IR) spectra of Bz-W n (n = 1-10) clusters has been carried out by Prakash et al. 17 using hybrid meta DFT based M05-2X method with the 6-31+G** basis set. Using Bader's theory of atoms in molecule (AIM) approach, it is found that the nonconventional H-bonding interactions, such as O-H⋯π interactions are present in all these clusters with additional contributions from C-H⋯O and lone pair (lp)⋯π interactions, leads to overall stability of these clusters. It is also noted that π H-bonded O-H stretching vibrational frequencies are red shifted in all the clusters. The inverted book conformer of water hexamer is found to have highest binding energy in compared to all other water clusters. 17 Slipchenko et al. 42 investigated the structures and bonding in water W n (n = 1-8)-Bz n (n = 1-2) complexes computationally using an effective fragment potential (EFP) method. It is predicted that benzene can act both as an H-bond donor and acceptor in the water-benzene complexes and all of the larger water-benzene complexes are dominated by H-bonding interactions i.e., O-H⋯π and C-H⋯O interactions. The interactions between such systems are considered complicated due to the existence of various possible structures that depend on the number and type of hydrogen bonds involved in these complexes as the benzene-water interactions are weaker than the water-water interactions and the very small energy differences between different isomers are involved in such complexes.
Excited state properties of such systems have received very little attention as of yet. Computational studies on low lying excited states of Bz-W n (n = 1-6) complexes have been performed by Upadhyay et al. 43 using configuration interaction method involving all the singly excited configurations (CIS). The binding energies of these clusters were also calculated and complex of water dimer with benzene (Bz-W 2 ) is found to be the most stable among various complexes, with a high stability of Bz-W 6 complex as well. It is found that the whole Bz-W n complex as well as the water cluster W n undergoes an expansion in size following excitation to lowest singlet excited states. The small blue shifts in the excitation energies are observed in the electronic absorption spectra in going from Bz to higher complexes (Bz-W n (n = 4-6)). Photon-stimulated desorption in an astrophysical context has not been studied extensively until recently. A few experiments have been performed on photo processing of laboratory models of Bz-ice complexes by UV radiation. [44][45][46] Experimental study 44 on photo-processes in model interstellar ices by modeling multilayer films of benzene & water deposited on a sapphire substrate at a temperature of around 80 K has investigated three distinct photodesorption mechanisms in such systems: (i) direct adsorbatemediated desorption of benzene; (ii) indirect adsorbatemediated desorption of water; (iii) substrate-mediated desorption of both benzene and water. The translational temperature of both desorbed species i.e., benzene and water molecule is found to be more than the ambient temperature of the complex system. Recent experimental study on photon-and electron-induced desorption from laboratory models of interstellar ice grains by mimicking the conditions found in dense interstellar clouds has been performed by Thrower et al. 46 The desorption cross sections and first order rate coefficients for the desorption processes in benzene-water (Bz-W n ) complexes have been obtained. It is observed that photon absorption by benzene can make H 2 O desorption possible at wavelengths where photon-absorption cross-section for H 2 O is negligible.
The electronically excited states play an important role in the photochemistry and electronic spectroscopy and are the subject of our investigations detailed below. The main aim of the present study is to investigate the spectroscopy and photochemistry of interstellar ice analogs i.e., Bz-W 6 cluster quantum mechanically to provide insight into elementary process involved in their processing, and to study important electronic transitions involved in such systems using time dependent DFT with range of well developed DFT functionals for response theory.
It has been observed from earlier studies that the water hexamer (H 2 O) 6 cluster is the smallest water cluster that allows non-cyclic structures and a more three dimensional structure, it can be considered as the building block of many ice forms. 1, 3 We chose a water hexamer (W 6 ) cluster since this is frequently taken as a benchmark system for many computational chemistry studies. [1][2][3][47][48][49] The hydrogen-bond arrangements of these six water molecules in the cluster determine the shape and stability of the cluster. A few low-lying isoenergetic conformers of the water hexamer cluster are ring, book, cage and prism with 6, 7, 8 and 9 number of hydrogen bonds, respectively. There have been uncertainties in literature on predicting the exact order in energetics of very close lying water hexamer conformers. Previous theoretical investigations have found the cage shaped water W 6 cluster, with four dangling hydrogen atoms (non-hydrogen bonded atoms) and eight hydrogen bonds, to be the most stable one giving the minimum energy, among the other low lying conformers i.e., prism, book, ring and chair forms. 2,5,11,49 DFT-D calculations on the water hexamer has predicted the energetic order prism < cage < book < chair. 50,51 Most recent computational studies have predicted prism, cage and book conformers as the three lowest energy conformers with ZPE correction, with prism as the lowest energy structure, and cage conformer is found to be very closer to prism in energy. 9,47 A very recent theoretical study on the water hexamer using full dimensional Diffusion Quantum Monte-Carlo simulations has predicted the presence of cage and prism conformers at low temperatures, 10 and in agreement with the experimental measurements of broadband rotational spectra of the water hexamer formed in supersonic expansion. 48 Experimental study by Pate et al. has also established the cage conformer of water hexamer as the global minimum energy structure. 48 Taking account of both theoretical 9 and experimental 48 evidences of existence of cage and prism conformers at low temperatures, that are relevant to astrophysical interstellar conditions, and are also established as lowest energy conformers of water hexamer by very recent studies, 9,47,48 we have chosen both the cage and prism structures as model systems in our present study.
The benzene (Bz) molecule is then brought closer to the water W 6 cluster and it binds to the water cluster through hydrogen bonded interactions and thus form Bz-W 6 complex system. The benzene molecule binds to the water cluster mainly through O-H⋯π hydrogen bonds where one of the free dangling hydrogen atom of the water cluster points toward the π electron cloud of the benzene ring. 17,40,42 Small contributions from C-H⋯O and lone pair (lp)⋯π hydrogen bonded interactions also stabilize the Bz-(H 2 O) 6 cluster. 17 Both cage and prism geometries of water W 6 cluster interact with Bz to give Bz-W 6 cluster and are predicted to retain their cage and prism structures in Bz-W 6 cluster too.

Ground state structures
The ground state geometries of both the cage and prism conformers of W 6 clusters and Bz-W 6 clusters were optimized using density functional theory with Truhlar's meta hybrid functional M05-2X functional and the 6-31+G(d,p) basis set. M05-2X & M06-2X are non-local functionals with double the amount of non-local (Hartree-Fock) exchange, and are found to perform better than standard hybrid functionals for systems involving non-covalent interactions, and modeling electronic excitation energies to both valence and Rydberg states. 52,53 In order to properly investigate any long range correction due to dispersion effects in such complex systems, the obtained optimized ground state geometries were further reoptimized using the following range of hybrid functionals LC-wPBE, wB97X, wB97XD, B2PLYP, B2PLYPD, in addition to second order Moller-Plesset perturbation theory (MP2), and employing the augmented correlation-consistent polarizedvalence double zeta basis set (aug-cc-pvDZ). Long range corrected functionals such as LC-wPBE, CAM-B3LYP, wB97X, wB97XD functionals account the non-Coulomb part of exchange functionals, which typically dies off too rapidly and gives inaccurate results at large distances, where wB97XD functional also includes an empirical atom-atom dispersion corrections. [54][55][56] Fig. 1(a) and (b) shows the MP2/aug-cc-pVDZ optimized geometries of cage and prism form of Bz-W 6 cluster, respectively.
The binding energy (BE) of the Bz-W 6 cluster is calculated using where E Bz-W 6 , E W 6 and E Bz denote the total energy of Bz-W 6 cluster, W 6 cluster and benzene, respectively. Table 1 shows the optimized total energies of the W 6 clusters, Bz-W 6 clusters, benzene (Bz) and absolute values of binding energies (BE's) of Bz-W 6 clusters (with and without zero point energy (ZPE) correction). It is seen that BE's of the Bz-W 6 cluster (with ZPE) ranges from 2.91-8.73 kcal mol −1 for cage conformer and 2.69-8.08 kcal mol −1 for the prism conformer. In order to get correct energetics of both the cage and prism forms of Bz-W 6 clusters, basis set superposition error (BSSE) is also calculated using the counterpoise (CP) method and corrected binding energies are also listed in Table 1. Bindings energies are reduced noticeably after applying BSSE correction and it ranges from 2.67-6.14 kcal mol −1 for the cage conformer and 2.47-5.65 kcal mol −1 for the prism conformer. Fig. 2(a) and 2(b) depict the variation of BE's of the prism and cage forms of Bz-W 6 clusters versus different DFT hybrid functionals and MP2 with ZPE and BSSE correction, respectively. It is clearly seen from these results that the cage conformer is more stable than the prism conformer, for all of the various functionals employed, as the cage form is found to have higher binding energy (BE) than the prism form, which is also consistent with previous computational studies. 17,40,49 The magnitude of BE's of Bz-W 6 cluster (with ZPE correction) for both the cage and prism forms with different computational methods vary in the following order: LC-wPBE < B2PLYP < wB97X < wB97XD < B2PLYPD < MP2. With ZPE correction, Bz-W 6 cluster is found to be most stable at the MP2 level with the highest BE, and least stable at the LC-wPBE level of calculation. However with BSSE correction, the 'wB97XD' functional gives higher binding energies than MP2, that clearly indicates that BSSE error is higher for wavefunction methods and energies are overestimated by MP2 level of calculations as shown in Fig. 3(a) and 3(b) for both prism and cage conformers, respectively. The MP2 method is used extensively in computational chemistry and is considered reliable for ground state geometry optimizations and to estimate the energies of the loosely bound hydrogen bonded and dispersion bound complexes by accounting the electron-correlation effects including dispersion. 57 Taking into account the reliability and success of MP2 method in the past, and highest BE values with ZPE correction for our Bz-W 6 clusters, we have  used MP2 optimized ground state geometries for main discussion on TD-DFT studies.

Electronic excitations in Bz-W 6 clusters
In recent years, TD-DFT linear response theory 58-61 has become the most widely used electronic structure method for calculating vertical electronic excitation energies. These TD-DFT calculations provide the Ultraviolet (UV) absorption spectra giving the information about the excitation wavelength and oscillator strength ( f ) of each excited state. The oscillator strength can be written as, 62 where ω ij is the frequency of transition from i th to j th quantum state and |〈i|μ|j〉| 2 is the associated transition dipole moment. We have performed time-dependent DFT calculations on both MP2/aug-cc-pVDZ and wB97XD/aug-cc-pVDZ optimized ground state geometries of Bz, W 6 and Bz-W 6 clusters using three different functionals i.e., B3LYP, CAM-B3LYP and M06-2X with the 6-31++G(d,p) basis set.
For Bz-W 6 prism shaped geometry, the UV spectra results obtained from TD-DFT calculations on MP2 optimized ground state Bz-W 6 prism shaped cluster, using all three functionals are presented in Fig. 4(a). Fig. 4(a) shows the diverse spectra for three different functionals where B3LYP functional generates the strongest peak at around 182 nm, with other less intense peaks within the range 174-190 nm, while CAM-B3LYP functional predicts transitions over the spectral range 157-180 nm, with strongest one centered around 180 nm and less intense peaks at around 157 nm and 175 nm. M06-2X and CAM-B3LYP functionals are consistent to generate the strongest intensity peak at 180 nm. It is also noted that B3LYP predicts few less intense peaks above 180 nm. In order to investigate the effect of Bz interaction with isolated water cluster W 6 and also to investigate how it affects excitation features in Bz-W 6 cluster, it is important to compare the UV spectra of Bz-W 6 cluster with UV spectra of isolated water cluster W 6 and benzene Bz. The performances of all three functionals i.e., B3LYP, CAM-B3LYP and M06-2X are also tested on prism-shaped W 6 & Bz-W 6 clusters and Bz molecule (see Fig. 4

(b)-4(d)).
It is also important to investigate the nature of these transitions and to understand whether Bz or W 6 cluster excitations, plays the dominant role in Bz-W 6 clusters. So, we have also analyzed the orbital transformations for those excitations that are associated with the high oscillator strengths or high intensity peaks found in the prism Bz-W 6 spectra for all three functionals, and are presented in Table 2. Fig. 4(b) clearly indicates different excitation features in Bz-W 6 prism cluster with respect to excitation features in isolated W 6 prism cluster and Bz molecule given by B3LYP functional. The excitations in Bz-W 6 cluster are shifted towards longer wavelengths (above 174 nm extended up to 193 nm) with respect to W 6 cluster where electronic excitations are dominant only in shorter wavelength region i.e., below 176 nm. The Bz excitation feature at around 181 nm is of very high intensity, with another less intense peak at 190 nm. In Bz-W 6 prism cluster, it is interesting to see that one of the strong intensity peak at around 182 nm is due to a charge transfer feature of the W 6 cluster to Bz system. Another strong peaks at around 184 nm and 185 nm is identified as the π→π* electronic transition of benzene. Few weak excitations showing partial charge transfer from benzene to W 6 moiety and Bz locally diffuse excitations are also seen at wavelengths around 179 nm and 190 nm, respectively. It is interesting to note that presence of Bz enhances the excitations in W 6 cluster towards longer wavelengths above 176 nm, not present in isolated W 6 cluster. The intensity of peaks associated with water excitations in Bz-W 6 prism cluster are also found to be higher than in isolated water W 6 cluster. Fig. 4(c) shows that the CAM-B3LYP functional predicts a few excitations below 170 nm region in Bz-W 6 cluster which  are in close proximity to excitations in W 6 cluster. Similarly peaks observed at longer wavelengths at around 175 nm and 180 nm in the Bz-W 6 cluster (completely absent in W 6 cluster) are found to be closer to Bz excitation feature at 178 nm. Using the CAM-B3LYP functional, it is found that the strongest peak at around 180 nm is the π→π* electronic transition of benzene. A second intense peak at 157 nm is predicted due to the combined influence of σ→π* transition from σ bonding orbital to π anti-bonding orbital of the benzene at 157.1 nm and the benzene charge transfer (CT) state at 157.4 nm, while third peak at around 175 nm corresponds to locally diffuse charge state associated with Bz excitation. We can also see that Bz excitations in Bz-W 6 prism cluster are red-shifted as compared to isolated Bz spectra. Thus, it is realized that Bz excitations are also influenced by presence of W 6 cluster around it.
It is seen from Fig. 4(d) that Bz-W 6 spectra are red shifted with respect to W 6 spectra, using M06-2X functional. Since, excitations observed in Bz system are above 175 nm wavelength region and therefore more likely to influence excitations at longer wavelengths in Bz-W 6 cluster rather than W 6 cluster. The M06-2X functional predicts the π→π* transition feature of benzene at around 180 nm for the strongest peak, while another strong peak at around 177 nm is due to combined effect of π→π* transition feature of benzene at 177.36 nm, and benzene charge transfer (CT) state at 176.56 nm, in which charge from Bz ring is transferred to W 6 cluster in prism shaped Bz-W 6 complex system.
It is found that most of the peak excitations in Bz-W 6 prism cluster are influenced by Bz excitations as compared to water W 6 excitations. It is also noted that intensity of the strongest peak is much higher in Bz and Bz-W 6 cluster than those in W 6 cluster, which hold for all three functionals. Among various long range corrected functionals, wB97XD functional is found to be one of the most promising DFT functional used for systems involving general non-covalent interactions. 67 It is already mentioned above that wB97XD provides highest binding energies among all DFT functionals on Bz-W 6 clusters with BSSE correction, and it is same for both cage and prism conformers (see Table 1 and Fig. 3). Therefore, TD-DFT calculations on wB97XD-optimized ground state geometries of Bz-W 6 clusters with three functionals i.e., B3LYP, CAM-B3LYP and M06-2X are also performed.
Similar UV spectral characteristics are obtained from TD-DFT calculations on the wB97XD optimized ground state Bz-W 6 prism shaped cluster, where UV spectra undergo a red shift in going from the CAM-B3LYP to the B3LYP functional (see ESI: Fig. S1 and Table S1 †). Again we clearly see that the strongest intensity peak associated with the benzene π→π* transition undergoes a small red shift in going from Bz to Bz-W 6 cluster for both CAM-B3LYP and M06-2X functionals, and the magnitude of shift is slightly larger for MP2 optimized geometries than wB97XD optimized geometries. The intensities are also found to be much stronger in Bz-W 6 cluster than those in W 6 cluster. It is again noted that presence of Bz enhances the excitations in W 6 cluster towards wavelengths above 170 nm. The strong intensity π→π* transition features of Bz has shown small blue shift of order of around 2-3 nm in wB97XD optimized ground state geometries of Bz-W 6 cluster with respect to MP2 geometries. Now, focusing on cage form of Bz-W 6 geometry, UV spectra results obtained from TD-DFT calculations on the MP2 optimized Bz-W 6 cage cluster using all three B3LYP, CAM-B3LYP and M06-2X functional are shown in Fig. 5(a). The important electronic transitions corresponding to high oscillator strengths or strong intensity peaks in UV spectra of cage shaped Bz-W 6 cluster for all three functionals, are presented in Table 3. Fig. 5(a) shows that in Bz-W 6 cage cluster, electronic excitations predicted by B3LYP functional are red shifted with respect to excitations in M06-2X and CAM-B3LYP level of calculations. Both CAM-B3LYP and M06-2X functionals are consistent to find strong intensity peaks very close to each other at wavelength around 179-180 nm. However, a few less intense peaks at shorter wavelengths below 174 nm are also observed in CAM-B3LYP level of calculations. Fig. 5(b) shows that excitations in Bz-W 6 cage cluster are shifted towards longer wavelengths (region above 180 nm) with respect to W 6 cluster using B3LYP functional. It is interesting to notice that excitations in isolated water W 6 cluster above 170 nm wavelength range are rarely observed, however the B3LYP functional generates some weaker transitions at around 180-181 nm. It is found that strongest peaks generated at around 183 nm in UV spectra corresponds to contributions from π→π* transition of benzene, and from a locally diffuse state due to Bz excitation. The another strong peak at around 182 nm is associated with π→π* transition of benzene while other peaks at around 180 nm and 181 nm correspond to Bz charge transfer excitation and to local Rydberg type state associated with W 6 cluster excitation, respectively. It is also again evident from spectra that intensity of water excitations increases in Bz-W 6 cluster as compared to the isolated water W 6 cluster. Fig. 5(c) indicates that excitations found in Bz-W 6 cage cluster do not overlap with the peaks found in isolated W 6 cage cluster for CAM-B3LYP functional. However, some weak excitation peaks at shorter wavelengths below 170 nm lie close to the W 6 cluster excitation region and some of the strong peaks at longer wavelengths (above 170 nm) lie close to Bz strongest intensity peak centered at around 178 nm. The analysis of these excitations has shown that with CAM-B3LYP, two close lying peaks with strong intensities at around 179 nm and 180 nm are associated with π→π* transition of benzene. The third peak at around 161 nm corresponds to locally diffuse Rydberg state due to W 6 cluster excitation while another peak at around 174 nm seems like Bz charge transfer state.
From Fig. 5(d) we see that peaks in Bz-W 6 cluster are red shifted (above 170 nm) with respect to W 6 cluster with W 6 excitations lies below 161 nm as given by M06-2X functional. The M06-2X functional accounts for the π→π* transition of benzene at around 180 nm and 179 nm which corresponds to strong intensities too, while another weak excitation at around 178 nm seems to be due to a benzene charge transfer excitation.
Benzene is found to dominate all the excitations in the Bz-W 6 cage cluster with M06-2X. It is again noted that once again the intensities of peaks are much stronger in Bz-W 6 cage cluster than isolated water cage W 6 cluster.
The UV spectra results obtained from TD-DFT calculations on the wB97XD optimized ground state Bz-W 6 cage shaped cluster, using all three functionals are generally consistent with the one obtained on MP2 optimized geometries. Most of the strong intensity peaks are found around the wavelength range 175-180 nm for all three hybrid functionals (see ESI: Fig. S2 and Table S2 †). UV spectra obtained by M06-2X is quite similar to those predicted by CAM-B3LYP functional. B3LYP again shows transitions towards longer wavelengths with respect to transitions given by CAM-B3LYP and M06-2X functionals.
It is interesting to see that the peaks are slightly red-shifted for TD-DFT calculations on MP2 optimized geometries as compared to wB97XD optimized ground state geometries for both prism and cage Bz-W 6 clusters (see ESI: Fig. S3 and S4 †), which is true for all three functionals. The nature of the transitions in UV spectra for both MP2 and wB97XD optimized geometries are generally quite consistent.
We again note that water excitations are observed in Bz-W 6 cage cluster at longer wavelengths, which are completely absent in isolated cage W 6 system. The intensities of water excitations are also stronger in Bz-W 6 cluster as compared to isolated W 6 cluster. Again we clearly see that the benzene π→π* transition undergoes a small red shift in Bz-W 6 cluster with respect to isolated Bz. The intensities are also found to be much stronger in Bz-W 6 cluster than those in W 6 cluster.
It is worth mentioning that the lowest valence transitions of Bz from ground to excited states 1 1 B 2U and 1 1 B 1U are strictly dipole forbidden on symmetry grounds. The low lying vertical π→π* excited states in Bz-W 6 prism clusters show small red shift in excitation wavelengths, and oscillator strengths become slightly allowed as compared to dipole forbidden 1 1 B 2U and 1 1 B 1U excited states in Bz i.e., (Bz( f = 0.0000)→ Bz-W 6 ( f = 0.0001 − 0.0006)), indicating slight allowed character in restricted lowest energy π→π* electronic transitions due to the presence of prism W 6 conformer around Bz. However, no such difference has been noticed in cage conformer of Bz-W 6 cluster.
In UV spectra of both Bz-W 6 prism and cage clusters, the strong intensity π→π* electronic transitions (bright states) with high oscillator strengths are red-shifted towards longer wavelengths with respect to corresponding dipole allowed degenerate 1 E 1U excited states of isolated Bz. The magnitude of redshifts are about 3-4 nm in Bz-W 6 prism clusters, while Bz-W 6 cage clusters give red-shifts of about 2-3 nm. The degeneracies of these states are only slightly broken by about 1.0 nm in both Bz-W 6 cage and prism clusters. It is seen that cage Bz-W 6 cluster shows higher individual transition oscillator strengths than prism Bz-W 6 cluster (see Tables 2 and 3). However due to overlapping of very close lying peaks in prism Bz-W 6 cluster, the intensity of the strongest peak is larger in prism case than in cage one for CAM-B3LYP and M06-2X functionals, as degen- eracy of strongest intensity π→π* electronic transitions is broken slightly lesser in prism cluster than in cage cluster. The oscillator strengths of these bright π→π* electronic states are lower in Bz-W 6 clusters as compared to isolated Bz, illustrating the effect of water cluster around Bz, and existence of other close lying new Bz charge transfer and diffuse states, which also increase the intensities of water excitations and shift such water excitations towards longer wavelengths as also seen in recent experimental studies. [44][45][46] It is clear from above results and discussion that the benzene interaction with W 6 cluster plays a significant role in giving new excitation features in UV spectra of Bz-W 6 clusters. In order to calibrate basis set effects on the UV spectra, we have also performed TD-DFT calculations using Dunning's correlation-consistent aug-cc-pVTZ basis set on MP2 optimized ground state geometries of Bz, W 6 and Bz-W 6 clusters. For all three employed functionals, the peak electronic excitations show the small red-shift of around 1-2 nm for both isolated Bz and Bz-W 6 clusters relative to those electronic excitations generated using 6-31++G(d,p) basis set, while small blue shift of around 1 nm is seen in W 6 clusters, and that holds for both prism and cage conformers. The nature of peak excitations associated with higher oscillator strengths in Bz-W 6 clusters (for prism and cage conformers) is generally quite consistent for both aug-cc-pVTZ and 6-31++G(d,p) calculations (see ESI: Tables S3 and S4 †).

Conclusions
TD-DFT calculations have been performed on both MP2 and wB97XD optimized geometries for Bz, W 6 and Bz-W 6 clusters. We observe minor differences in the UV spectroscopy of cage and prism conformers of W 6 and Bz-W 6 clusters. Our results and discussions also show that TD-DFT calculations are functional-dependent, with CAM-B3LYP and M06-2X generally showing good agreement. We have calculated some interesting features for water W 6 cluster excitations at longer wavelengths above 170 nm in both Bz-W 6 cage and prism geometries which are absent in isolated water W 6 cluster, indicating benzene-mediated excitations in the W 6 cluster as observed in recent experimental studies. [44][45][46] The intensities of W 6 excitations are also found to increase in Bz-W 6 cluster relative to those in isolated W 6 clusters. It is predicted that for both MP2 and wB97XD optimized geometries, benzene π→π* transition undergoes small red shift in Bz-W 6 cluster with respect to isolated benzene, for both cage and prism geometries and that the degeneracy of this transition is very slightly broken compared to the isolated gas-phase benzene. Our results also predict that charge transfer (CT) states and locally excited diffuse states play an important role in such systems due to complicated nature of benzene-water (π⋯HO) hydrogen bonding interactions prevailing in such systems, including conventional π→π* type interactions between benzene and electrostatic interactions between hydrogen bonded water W 6 cluster.