2D THz-THz-Raman Photon-Echo Spectroscopy of ... - ACS Publications

Sep 6, 2017 - interactions. The structural specificity and ultrafast time resolution of 2D infrared spectroscopy has ... bandwidth by the sample cuvet...
0 downloads 39 Views 3MB Size
Subscriber access provided by - Access paid by the | UCSB Libraries

Letter

2D THz-THz-Raman Photon-Echo Spectroscopy of Molecular Vibrations in Liquid Bromoform Ian A. Finneran, Ralph Welsch, Marco A. Allodi, Thomas Francis Miller, and Geoffrey A. Blake J. Phys. Chem. Lett., Just Accepted Manuscript • DOI: 10.1021/acs.jpclett.7b02106 • Publication Date (Web): 06 Sep 2017 Downloaded from http://pubs.acs.org on September 7, 2017

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 free 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 accessible to all readers and 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.

The Journal of Physical Chemistry Letters 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 9

The Journal of Physical Chemistry Letters

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 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

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 ACS Paragon Plus Environment

Page 2 of 9

Page 3 of 9

The Journal of Physical Chemistry Letters

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 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

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 ACS Paragon Plus Environment

Page 4 of 9

Page 5 of 9

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 Letters

tion that matches the experimental spectrum (labeled X in Fig. 3(a)). These peaks are much weaker, since the Raman detection scheme is depolarized, while the ν3 vibration is polarized and ν6 is depolarized. Finally, a third set of peaks near f2 =5.3 THz (labeled IX in Fig. 3(a)) is observed in the experimental spectrum that is absent from the simulation; it is speculated that these features come from the coupling of vibrational levels with E/h>15 THz. The off-diagonal couplings recovered from the fit confirm the differences between the THz and Raman transitions (see SI). All of the multiquantum polarizability couplings fit to values of near zero, while the multi-quantum dipole couplings are relatively large. To further illustrate the effect of the couplings on the 2D spectrum, a simulation with the fit multi-quantum dipole and polarizability couplings switched is shown in Fig. 3(c)). In this case, strong peaks are predicted at f2 = ±1.9 THz, in complete disagreement with the experiment. Although the RDM simulations indicate the existence of multi-quantum THz (dipole) transitions, they are agnostic as to the atomistic origins of these couplings. The simplest cause of multi-quantum transitions is perhaps vibrational (mechanical) anharmonicity, which arises from the shape of the molecular potential energy surface as a function of the normal mode coordinates. 8,16 However, the intramolecular modes measured here may be expected to be relatively harmonic, 6 and indeed, previous combination and overtone measurements in the linear Raman spectrum of CHBr3 have shown that the vibrational anharmonicity is small for the ν6 and ν3 manifolds. 23 Additional calculations employing vibrational second-order perturbation theory (VPT2) 24 at the Hartree-Fock level also hint to this assumption (see SI). Other possible origins of multi-quantum transitions are coherence transfer (CT) and the ‘through space’ coupling of nearby monomers in the liquid. 25,26 However, such intermolecular coupling only leads to multi-quantum transitions in highly anharmonic modes, which is not the case for the ν6 and ν3 manifolds. 23,26 CT processes result from the spontaneous transfer of vibrational coherences from

many off-diagonal peaks in the first and second quadrants. We note that peaks VII and VIII were observed in our previous low-bandwidth experiment, 7 while all other features are new. Peaks in the first quadrant, (f1 =±,f2 =±), are due to non-rephasing pathways, while peaks in the second quadrant, (f1 =±,f2 =∓), arise from rephasing pathways. The frequencies of the fundamental 1-quantum transitions ν6 and ν3 have been marked with dashed gray lines in Fig. 3(a). Qualitatively, we find that the strongest peaks align with the expected 1-quantum transitions on the f2 axis, while the f1 positions are distributed between 0.5-8.5 THz, likely due to multi-quantum transitions. To interpret these results, we fit the experimental data to a previously developed reduced density matrix (RDM) simulated spectrum with vibrational eigenstates up to 15 THz (Fig. 3(b,d)). 7 The fundamental eigenstate energies and relative intensities were taken from the linear spectrum, 7 while the excited eigenstate energies and 1-quantum dipole/polarizability couplings were fixed at values consistent with a harmonic oscillator. Multi-quantum dipole and polarizability couplings were fit with a basin hopping algorithm 21 and quasi-Newton local minimizer. 22 The fitness function was calculated with an L2-norm and regularization. Alternative models for the couplings were additionally considered and found to poorly fit the experimental spectra. The full details of the RDM modeling are given in the SI. The simulated RDM spectrum, Fig. 3(b), is in qualitative agreement with the experimental spectrum. A representative set of the pathways that contribute to the experimental peaks are highlighted with Feynman diagrams in Fig. 3(e). This analysis indicates that the strongest non-rephasing and rephasing peaks I-VIII are due to multi-quantum ν6-ν3 difference band transitions with the THz pulses, and 1-quantum ν6 transitions with the Raman pulse. No ν6 overtone transitions are observed, likely due to the limited bandwidth of the THz pulses. If the THz bandwidth were to exceed 9.4 THz, these transitions should appear at f1 =4.7, 9.4 THz and f2 =4.7 THz. The simulation also shows a set of peaks near the 1-quantum ν3 transi-

ACS Paragon Plus Environment

5

The Journal of Physical Chemistry Letters

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

To gain better insight into the molecular origin(s) of these effects, we performed ab initio calculations on the dipole and polarizability surfaces of an isolated bromoform molecule. The resulting dipole coupling elements are given in the SI. In these calculations, no significant nonlinearities in either the dipole or polarizability surface were found. These results suggest that the nonlinearities are a property of the liquid and not of the single molecule. We therefore performed additional ab initio calculations of the dipole and polarizability surfaces of a bromoform molecule surrounded by seven bromoform monomers. The employed configuration was sampled from a molecular dynamics simulation. Again, no significant nonlinearities were found. Due to the large computational cost of these calculations, however, only a central bromoform and its first solvation shell were treated. Thus, no longer range effects were included. Further, only a single representative configuration of the cluster was selected from an MD simulation of the liquid, and coupled motions between the central bromoform and its first solvation shell were not considered. Since 2D TTR is only sensitive to pathways containing a multi-quantum transition, it is possible that the experiments are strongly selective to relatively rare configurations of the local liquid structure that exhibit large dipole nonlinearities; it is also possible that coupled intramolecular motions originate the 2D TTR signal in this case. To explore the effect of conformational sampling, dipole and polarizability transition matrix elements were calculated using a dipole-induced dipole (DID) 28 model averaging over 2560 cluster configurations. As before, the clusters included only a central CHBr3 and its first solvation shell, and only the motions of the central bromoform were considered. Again, no significant nonlinearities in the coupling elements were found, which might be due either to inaccuracies in the employed force field and the DID model or to the neglect of coupled intramolecular motions. In conclusion, we have measured the broadband 2D THz-THz-Raman spectrum of room temperature liquid bromoform. Many offdiagonal peaks are observed in the spectrum,

interactions with the surrounding ‘bath’, without the loss of phase memory, 25 via dynamical processes, not to their direct excitation during the light-matter interaction via modifications of the dipole and/or polarizability surfaces. The transfer energetics are limited to ∼kT, or 6.3 THz at 300 K. 27 To investigate this mechanism, CT processes were included phenomenologically in our RDM model; while they were found to yield multi-quantum transitions on both the f1 and f2 axes, it was not found that they lead to transitions along only one of the axes, as is seen experimentally. Since CT only changes the relaxation dynamics of the density matrix, it cannot prefer dipole versus polarizability coupling elements in the Hamiltonian. Therefore, CT effects cannot explain the observed multi-quantum dipole and 1-quantum polarizability transitions. Further details are given in the SI. A remaining possibility is that the TTR signatures arise from nonlinearities in the molecular dipole moment (~µ) and polarizability (˜ α) as a function of the vibrational coordinates (q), which lead to multi-quantum THz and Raman transitions, respectively. For example, dipole nonlinearity can be described with 8,16

µ ~ =~µ(q0 ) +

X  ∂~µ  i

∂qi

qi q0

1 X  ∂ 2µ ~  qi qj + ... + 2 i,j ∂qi ∂qj q0

Page 6 of 9

(2)

where the Taylor expansion summations run over the normal modes. The first term is due to the permanent dipole moment, giving rise to rotational spectra, while the second corresponds to the first dipole derivative and 1-quantum vibrational transitions. It is the higher order terms that generate multi-quantum transitions. This is a likely source of the transitions in the bromoform spectrum, assuming large nonlinearities in the dipole surface and small nonlinearities in the polarizability surface. Particular multi-quantum couplings appear to be strongly favored over others in the RDM simulations, suggesting an interesting underlying mechanism.

ACS Paragon Plus Environment

6

Page 7 of 9

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 Letters

which we have assigned to coupling between the ν3 and ν6 modes. The most prominent features correspond to multi-quantum dipole and 1-quantum polarizability transitions, which are not predicted by electronic structure theory calculations of isolated bromoform molecules or from a cluster model of a representative solvation structure. A more extensive theoretical investigation of the ensemble of solvent configurations and coupled motions will be needed to reveal the molecular origin of the multi-quantum transitions from first principles.

(3) Shim, S.-H.; Gupta, R.; Ling, Y. L.; Strasfeld, D. B.; Raleigh, D. P.; Zanni, M. T. Two-dimensional IR Spectroscopy and Isotope Labeling Defines the Pathway of Amyloid Formation with Residue-Specific Resolution. Proc. Natl. Acad. Sci. U.S.A. 2009, 106, 6614–6619. (4) Woerner, M.; Kuehn, W.; Bowlan, P.; Reimann, K.; Elsaesser, T. Ultrafast TwoDimensional Terahertz Spectroscopy of Elementary Excitations in Solids. New J. Phys. 2013, 15, 025039.

Acknowledgement The authors acknowledge Dr. Mostafa Shalaby (Paul Scherrer Institut, Switzerland) for suggesting the use of a diamond window and polarizer beam combiner, and the National Science Foundation (Grants CHE-1214123 and CHE-1057112) for financial support. We also thank Matt Welborn, Ioan Magdau, Brett Savoie, and Feizhi Ding for helpful discussions. R.W. acknowledges financial support from the Deutsche Forschungsgemeinschaft under grant number WE 5762/1-1. M.A.A. acknowledges current support from a Yen Postdoctoral Fellowship from the Institute for Biophysical Dynamics at the University of Chicago.

(5) Lu, J.; Zhang, Y.; Hwang, H. Y.; Ofori-Okai, B. K.; Fleischer, S.; Nelson, K. A. Nonlinear Two-Dimensional Terahertz Photon Echo and Rotational Spectroscopy in the Gas Phase. Proc. Natl. Acad. Sci. U.S.A. 2016, 113, 11800– 11805. (6) Savolainen, J.; Ahmed, S.; Hamm, P. Two-Dimensional Raman-Terahertz Spectroscopy of Water. Proc. Natl. Acad. Sci. U.S.A. 2013, 110, 20402–20407. (7) Finneran, I. A.; Welsch, R.; Allodi, M. A.; Miller, T. F.; Blake, G. A. Coherent TwoDimensional Terahertz-Terahertz-Raman Spectroscopy. Proc. Natl. Acad. Sci. U.S.A. 2016, 113, 6857–6861.

Supporting Information Available: Experimental Setup, Reduced Density Matrix (RDM) Modeling, Coherence Transfer, Anharmonic Vibrational Calculations, Calculation of Dipole and Polarizability Non-Linearities. This material is available free of charge via the Internet at http://pubs.acs.org/.

(8) Tokmakoff, A.; Lang, M.; Larsen, D.; Fleming, G.; Chernyak, V.; Mukamel, S. Two-Dimensional Raman Spectroscopy of Vibrational Interactions in Liquids. Phys. Rev. Lett. 1997, 79, 2702–2705.

References

(9) Tanimura, Y.; Mukamel, S. TwoDimensional Femtosecond Vibrational Spectroscopy of Liquids. J. Chem. Phys. 1993, 99, 9496–9511.

(1) Fecko, C.; Eaves, J.; Loparo, J.; Tokmakoff, A.; Geissler, P. Ultrafast Hydrogen-Bond Dynamics in the Infrared Spectroscopy of Water. Science 2003, 301, 1698–1702.

(10) Hamm, P.; Shalit, A. Perspective: Echoes in 2D-Raman-THz Spectroscopy. J. Chem. Phys. 2017, 146, 130901.

(2) Ramasesha, K.; De Marco, L.; Mandal, A.; Tokmakoff, A. Water Vibrations have Strongly Mixed Intra- and Intermolecular Character. Nat. Chem. 2013, 5, 935–940.

(11) Ivanecky, J. E.; Wright, J. C. An Investigation of the Origins and Efficiencies of Higher-Order Nonlinear Spectroscopic

ACS Paragon Plus Environment

7

The Journal of Physical Chemistry Letters

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

Page 8 of 9

(21) Wales, D. J.; Doye, J. P. Global Optimization by Basin-Hopping and the Lowest Energy Structures of Lennard-Jones Clusters Containing up to 110 Atoms. J. Phys. Chem. A 1997, 101, 5111–5116.

Processes. Chem. Phys. Lett. 1993, 206, 437–444. (12) Blank, D. A.; Kaufman, L. J.; Fleming, G. R. Fifth-Order Two-Dimensional Raman Spectra of CS2 are Dominated by Third-Order Cascades. J. Chem. Phys. 1999, 111, 3105–3114.

(22) Perez, R. E.; Jansen, P. W.; Martins, J. R. R. A. pyOpt: A Python-Based ObjectOriented Framework for Nonlinear Constrained Optimization. Struct. Multidiscipl. Optim. 2012, 45, 101–118.

(13) Hamm, P.; Savolainen, J. TwoDimensional-Raman-Terahertz Spectroscopy of Water: Theory. J. Chem. Phys. 2012, 136, 094516.

(23) Fernandez-Liencres, M. P.; Navarro, A.; Lopez, J. J.; Fernandez, M.; Szalay, V.; de los Arcos, T.; Garcia-Ramos, J. V.; Escribano, R. M. The Force Field of Bromoform: A Theoretical and Experimental Investigation. J. Phys. Chem. A 1996, 100, 16058–16065.

(14) Shalit, A.; Ahmed, S.; Savolainen, J.; Hamm, P. Terahertz Echoes Reveal the Inhomogeneity of Aqueous Salt Solutions. Nat. Chem. 2017, 9, 273–278. (15) Allodi, M. A.; Finneran, I. A.; Blake, G. A. Nonlinear Terahertz Coherent Excitation of Vibrational Modes of Liquids. J. Chem. Phys. 2015, 143, 234204.

(24) Bloino, J.; Barone, V. A Second-Order Perturbation Theory Route to Vibrational Averages and Transition Properties of Molecules: General Formulation and Application to Infrared and Vibrational Circular Dichroism Spectroscopies. J. Chem. Phys. 2012, 136, 124108.

(16) Park, K.; Cho, M. Time-and FrequencyResolved Coherent Two-Dimensional IR Spectroscopy: Its Complementary Relationship with the Coherent Two-Dimensional Raman Scattering Spectroscopy. J. Chem. Phys. 1998, 109, 10559–10569.

(25) Khalil, M.; Demird¨oven, N.; Tokmakoff, A. Vibrational Coherence Transfer Characterized with Fourier-Transform 2D IR Spectroscopy. J. Chem. Phys. 2004, 121, 362–373.

(17) Dexheimer, S. L. Terahertz Spectroscopy: Principles and Applications; CRC press: Boca Raton, U.S.A., 2007.

(26) Hamm, P.; Zanni, M. Concepts and Methods of 2D Infrared Spectroscopy; Cambridge University Press: New York, U.S.A., 2011.

(18) Hoffmann, M. C.; Brandt, N. C.; Hwang, H. Y.; Yeh, K.-L.; Nelson, K. A. Terahertz Kerr Effect. Appl. Phys. Lett. 2009, 95, 231105.

(27) Baiz, C. R.; Kubarych, K. J.; Geva, E. Molecular Theory and Simulation of Coherence Transfer in Metal Carbonyls and Its Signature on Multidimensional Infrared Spectra. J. Phys. Chem. B 2011, 115, 5322–5339.

(19) Sajadi, M.; Wolf, M.; Kampfrath, T. Terahertz-Field-Induced Optical Birefringence in Common Window and Substrate Materials. Opt. Express 2015, 23, 28985– 28992.

(28) Van Duijnen, P. T.; Swart, M. Molecular and Atomic Polarizabilities: Thole’s Model Revisited. J. Phys. Chem. A 1998, 102, 2399–2407.

(20) Shalaby, M.; Vicario, C.; Hauri, C. P. Extreme Nonlinear Terahertz Electro-Optics in Diamond for Ultrafast Pulse Switching. APL Photonics 2017, 2, 036106.

ACS Paragon Plus Environment

8

ournal Pageof 9 of Physical 9 Chemistry L

1 2 3 4 5 ACS 6 Paragon Plus Environmen 7 8