Thermal Population Fluctuations in Two-Dimensional Infrared

Subscriber access provided by READING UNIV

Article

Thermal Population Fluctuations in Two-Dimensional Infrared Spectroscopy Captured with Semiclassical Mechanics Prashanth Ramesh, and Roger F. Loring J. Phys. Chem. B, Just Accepted Manuscript • DOI: 10.1021/acs.jpcb.7b12122 • Publication Date (Web): 31 Jan 2018 Downloaded from http://pubs.acs.org on February 1, 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 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 B 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

Thermal Population Fluctuations in Two-Dimensional Infrared Spectroscopy Captured with Semiclassical Mechanics Prashanth Ramesh and Roger F. Loring Department of Chemistry and Chemical Biology, Baker Laboratory, Cornell University, Ithaca, New York, 14853, USA (Dated: January 29, 2018)

Abstract Time-resolved two-dimensional (2D) infrared spectra of the asymmetric stretch mode of solvated CO2 show distinct features corresponding to ground and excited state thermal populations of the bend modes. The time-dependence of these peaks arises in part from solvent-driven thermal fluctuations in populations of the lower frequency bend modes through their coupling to the higher frequency asymmetric stretch. This observation illustrates the capacity of multidimensional vibrational spectroscopy to reveal details of the interactions among vibrational modes in condensed phases. The optimized mean-trajectory (OMT) method is a trajectory-based semiclassical approach to computing the vibrational response functions of multidimensional spectroscopy from a classical Hamiltonian. We perform an OMT calculation of the 2D vibrational spectrum for two coupled anharmonic modes, with the lower frequency mode undergoing stochastic transitions in energy to mimic solvent-induced fluctuations in quantum populations. The semiclassical calculation reproduces the influence of thermal fluctuations in the low-frequency mode on the 2D spectrum of the high-frequency mode, as in measured spectra of solvated CO2 .

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

I.

INTRODUCTION

As with their analogs in magnetic resonance, two-dimensional (2D) infrared spectra reflect not only the frequencies present in a mechanical system, but the interactions among degrees of freedom characterized by those frequencies. Diagonal peaks in the 2D spectrum record frequencies, while cross-peaks can reveal dynamics induced by interactions. The dependence of the shape and amplitude of cross peaks on the delay time in a four-wave-mixing measurement can reveal the dynamics of processes such as incoherent and coherent energy transfer among vibrational modes1,2 and the vibrational analog of chemical exchange in magnetic resonance.3–14 Infrared spectroscopy has been extensively applied15–22 to carbon dioxide as a solute in room temperature ionic liquids, because of the potential application of these solvents in carbon capture. Two-dimensional spectroscopy in particular probes structure and dynamics of solvation20,21 in these systems. Brinzer, et al.15,16 and Giammanco, et al.17 have measured 2D vibrational spectra of CO2 in various ionic liquids. The chromophore in these studies is the infrared-active asymmetric stretch, which undergoes an absorption red shift upon excitation of the bend mode. The 2D spectrum of the asymmetric stretch therefore shows distinct diagonal peaks corresponding to ground-state and first-excited-state populations of the bend modes. With increasing evolution time, diagonal peaks diminish in amplitude and cross-peaks appear, reflecting thermal fluctuations in population of the bend modes, driven by the solvent. This process is well-described by a two-state kinetic model15–17 of incoherent population transfer for the bend modes. Such studies exemplify that 2D spectroscopy of a high-frequency mode can report on the detailed dynamics of a multimode system to which it is coupled, here the bend modes and the solvent. Calculating 2D infrared spectra23,24 requires propagating nuclear degrees of freedom on a single (ground-state) electronic surface. Fully classical mechanical calculations25–31 of spectroscopic response functions can deviate qualitatively from quantum mechanical results at long times, although relaxation processes in condensed phases can limit accessible time scales to a regime in which classical dynamics is approximately applicable. Numerically correct quantum mechanical treatments32 of nonlinear infrared spectroscopy are not generally feasible for many degrees of freedom. Several approximation strategies have been employed to incorporate quantum dynamics into the calculation of nonlinear vibrational spectra, includACS Paragon 2Plus 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


ing the partition of a large system into a quantum subsystem interacting with a classical or semiclassical environment,33–38 the application of a tight-binding Hamiltonian with couplings dependent on structures generated in a molecular dynamics simulation,33,39–41 and mapping protocols connecting vibrational frequencies to collective coordinates, determined from quantum studies of subsystems.20,21,42–48 Another approach is to apply semiclassical approximations49–53 to quantum dynamics that, in principle, can be applied uniformly to all degrees of freedom. The optimized mean-trajectory approximation (OMT)54–58 is a semiclassical procedure for computing nonlinear vibrational response functions such as the third-order response function relevant to 2D spectroscopy.59 The method employs classical trajectories subject to semiclassical quantization rules, and is based on a mapping between the double-sided Feynman diagrams59–61 customarily used to represent processes contributing to a quantum spectroscopic response function and semiclassical diagrams54 representing radiation-induced discontinuous transitions between classical trajectories. The OMT circumvents62 the inherent numerical challenge of applying semiclassical propagators to nonlinear response functions: evaluating multidimensional integrals over highly oscillatory integrands.63 This significant benefit carries the cost of assuming the existence of good action and angle variables, even for anharmonic systems of many degrees of freedom for which such variables may not exist. However in the application to ultrafast condensed phase spectroscopy in which time scales are set by population relaxation or dephasing rates, application of the OMT only requires that the system behaves on short time scales approximately as though good action and angle variables existed. The OMT has been demonstrated to incorporate quantum anharmonic effects,54 line broadening from both dissipation and pure dephasing,56 and vibrational coherence transfer between near-resonant anharmonic vibrations.57 The OMT has been extended to include electronic as well as vibrational transitions.64 In the present work we apply the OMT to a model inspired by solvated carbon dioxide.15,17 This model is based on a classical Hamiltonian for two nonlinearly coupled anharmonic modes: a higher frequency mode qualitatively representing the asymmetric stretch and a lower frequency mode qualitatively representing one bend mode. Hamiltonian parameters are assigned by comparison to measured spectra of solvated CO2 . The solvent is not explicitly represented; rather, its effects are incorporated through a stochastic process that adds and removes energy from the bend mode. Brinzer, et al.15 analyzed 2D spectra of solvated CO2 ACS Paragon 3Plus Environment


with a stochastic model in which the frequency of the asymmetric stretch mode fluctuates between two values. The present model treats the effects of the solvent on the bend mode phenomenologically but fully includes the mechanics of the interacting bend and stretch modes within the model Hamiltonian. The OMT treatment of this model reproduces the phenomena observed in the measured spectra. The OMT approach is briefly summarized in Sec. II and the model Hamiltonian is defined. Calculated 2D spectra are presented and discussed in Sec. III, with conclusions given in Sec. IV. Technical details are reserved for Supporting Information.

II.

METHOD AND MODEL

In the OMT approximation, a vibrational response function is calculated from classical trajectories subject to semiclassical quantization conditions.54 Discontinuous transitions between trajectories represent radiation-matter interactions. For a single nuclear degree of freedom, the OMT method approximates the dynamics of a quantum operator proportional to |nl ihnr | with classical trajectories of action ~(nl + nr + 1)/2. Dynamics of a quantum population map onto classical trajectories with action equal to a half-odd-integral multiple of ~, while single quantum coherences are represented by trajectories with action equal to an integral multiple of ~. The three OMT diagrams representing the rephasing (echo) component of a third-order response function for a chromophore initially in its ground vibrational state are shown in Fig. 1. Each diagram is a schematic plot of chromophore action variable versus time. Horizontal solid lines represent continuous classical trajectories with time durations labeled and vertical dashed lines indicate action changes of ±~/2 resulting from the radiation-matter interaction. Solid vertices indicate phase space points at which values of momenta and coordinates are collected. Labels + and − indicate phases as defined below in Eq. (6). For the present model of an anharmonic chromophore mode c coupled to an anharmonic

ACS Paragon 4Plus Environment

Page 4 of 22

t2

- t1

t3

-

+

+

- t3 - t + + 1

t2

+ t3

-

t2

- - t1

FIG. 1:

+

Semiclassical OMT diagrams for the amplitude of the rephasing (echo) signal,

R++− (t3 , t2 , t1 ) in Eq. (1). The chromophore vibration is assumed to be thermally unexcited.

dark mode d, the sum of diagrams in Fig. 1 has the value56 R++− (t3 , t2 , t1 ) = ×γ(Jc(1)

−

i ~7

X 3

Z r

dz

(1)

r=1 (1) ~/2)γ(Jd )∆F (z (1) )δ(φ(2) (2)

(2)

Z dz

(2)

Z

dz (3)

(1)

− φt1 ) (1)

(1) ×δ(φ(3) − φt2 )∆(1) σr ∆σr0 Q− (z )Q+ (z t1 ) (3)

×Q+ (z (3) )Q− (z t3 ), ∞ X γ(J) = ~ δ(J − (n + 1/2)~),

(1) (2)

n=0

∆F (z) = F (z)|Jc →Jc −~/2 − F (z)|Jc →Jc +~/2 , exp(−βH(z)) F (z) = R 0 , dz exp(−βH(z 0 ))γ(Jd0 )γ(Jc0 ) (s+1) (s) (s) ∆± = ~2 δ Jc(s+1) − (Jc(s) ± ~/2) δ Jd − Jd ,

(3)

Q± (z) = (qc ∓ ipc /(mc ωc )) /2.

(6)

5

(4) (5)


Page 6 of 22

The index r in Eq. (1) labels the diagrams in Fig. 1. The factor r takes the value -1 for the first two diagrams in Fig. 1 and is 1 for the third diagram. This sign originates from quantum commutators within the response function. A point in the phase space of both degrees of freedom is denoted z. The factors γ(J) in Eq. (1), defined in Eq. (2), constrain initial (1)

action values Jc

(1)

and Jd (1)

of the dark mode Jd

with semiclassical quantization conditions. The initial action

is quantized to a half-odd-integral value of ~ representing an initial (1)

quantum population, while the initial chromophore action Jc

is quantized to an integral

value of ~ to represent the single-quantum coherence produced by the earliest radiationmatter interaction.54 Initial phase space points are further weighted58 by the finite difference ∆F in Eq. (3), with F (z) a modified58 classical thermal phase space distribution in Eq. (4). The finite-difference distribution arises from the effect on the thermal density operator of (s)

the first radiation-matter interaction. The factors ∆σ in Eq. (1) with σ → ±, defined in Eq. (5), depend on the identity of the diagram and constrain the chromophore action to change by ±~/2 with each radiation-matter interaction. Each vertex in a diagram in Fig. 1 labeled + or − corresponds to a phase space point at which Q± in Eq. (6) is determined. In Eq. (6), the effective mass and harmonic angular frequency of the chromophore are denoted mc and ωc . These factors give a semiclassical representation of the transition dipole, under the assumption that this operator is proportional to the chromophore coordinate, with the coefficient of proportionality suppressed. The amplitude of the absorptive 2D spectrum is defined by59 h

ˆ ++− (ω3 , −ω1 ; t2 ) Rabs (ω3 , ω1 ; t2 ) = Im R i ˆ +−+ (ω3 , ω1 ; t2 ) . +R

(7)

ˆ ++− (ω3 , ω1 ; t2 ) denotes the two-dimensional one-sided Fourier transform of R++− (t3 , t2 , t1 ) R in Eq. (1) with respect to t1 and t3 . The particular nonrephasing response function relevant to the absorptive spectrum R+−+ (t3 , t2 , t1 ) is obtained from Eq. (1) with the replacement (1)

(1)

Q− (z (1) )Q+ (z t1 ) → Q+ (z (1) )Q− (z t1 ). The time-evolution in Eq. (1) is carried out with a classical Hamiltonian for two anharmonically coupled normal modes with potential energy V (qc , qd ) = VM (qc ) + Vd (qd ) + ccd qc2 qd2 , Vd (qd ) = md ωd2 qd2 /2 + cd qd4 . ACS Paragon 6Plus Environment

(8) (9)

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 isolated chromophore mode has a Morse potential VM characterized by harmonic frequency ωc and well depth D. The dark mode has the quartic potential in Eq. (9), and the two modes are coupled with the biquadratic term in Eq. (8). Parameters in this model potential are assigned values to mimic the 2D infrared spectra of the asymmetric stretch mode of CO2 in an ionic liquid solvent measured by Brinzer, et al.15,16 These parameters are chosen with a two-step procedure, starting with an empirical truncated expansion in anharmonicity15,65 for the vibrational energy levels of the coupled chromophore and dark modes, E˜nd nc = ν˜d (nd + 1/2) + ν˜c (nc + 1/2) + χ˜dd (nd + 1/2)2 +χ˜cc (nc + 1/2)2 + χ˜cd (nc + 1/2)(nd + 1/2).

(10)

This expression contains two harmonic frequencies, two anharmonicity parameters, and one anharmonic coupling strength. These parameters are evaluated by comparison to measured15 spectroscopic transitions in solvated CO2 . The transition frequencies shown in Fig. 5 of Ref. 15 allow the evaluation of all but two of the parameters in Eq. (10): ν˜c = 2370 cm−1 , χ˜cc = −12 cm−1 , χ˜cd = −12 cm−1 , ν˜d /2 + χ˜dd = 337 cm−1 . Using a measured value65–68 of χ˜dd = 1.6 cm−1 for the anharmonicity of the bend mode of gas phase CO2 , we take ν˜d = 670 cm−1 . Next, expanding the chromophore potential in Eq. (8) to fourth order in qc and applying quantum mechanical perturbation theory65 permits the expression of potential energy anharmonicity parameters from Eqs. (8)-(9) in terms of the empirical parameters in Eq. 10, D = −(~ωc )2 /(4χcc ),

(11)

ccd = χcd mc md ωc ωd /~2 ,

(12)

cd = (2χdd /3) (md ωd /~)2 ,

(13)

with χrs = hcχ˜rs . Equations (11)-(13) together with the spectroscopic data specify the potential energy in Eqs. (8) and (9), giving a two-mode model whose 2D infrared spectrum resembles that of CO2 in the region of the frequency of the asymmetric stretch mode. The form of the model potential in Eqs. (8) and (9) was chosen to minimize the number of parameters while mimicking the relevant portion of the 2D spectrum of CO2 ; other choices could also have met these goals. The quartic form of the dark mode potential was chosen because it is even in coordinate and generates a nonzero value of χdd in the perturbation ACS Paragon 7Plus Environment


(A)

𝝎3

𝝎1 (B)

(C)

FIG. 2: (A) A schematic 2D spectrum illustrates the formation of cross peaks by spectral exchange during the t2 interval. Arrows colored red or blue indicate amplitude transfer from diagonal peaks to cross peaks that are respectively red or blue shifted in the ω3 dimension. (B) Representative OMT diagram and stochastic trajectory for the action of the dark mode producing a red shift. The t2 interval is colored red in both diagrams. (C) Representative OMT diagram and stochastic trajectory for the action of the dark mode producing a blue shift. The t2 interval is colored blue in both diagrams.

expansion in Eq. (10). The biquadratic coupling between the modes in Eq. (8) was chosen because it yields a nonzero value of χcd in Eq. (10); cubic terms proportional to qc qd2 or qc2 qd would also have contributed65 to this parameter. Our focus is the effect on the chromophore mode of thermal population fluctuations in ACS Paragon 8Plus 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 dark mode. Rather than extending the potential of Eq. (8) to include explicit solvent degrees of freedom coupled to the dark mode, we introduce a stochastic process that induces transitions in Jd , the action of the dark mode. In the OMT approximation, the evolution of the population of state nd maps to classical trajectories with action Jd = (nd + 1/2)~. Quantum transitions between nd = 0 and nd = 1 are then represented by transitions in classical action of magnitude ~. In contrast to the transitions induced by radiation in which the value of the angle is held constant,54 for these incoherent thermal transitions, the angle value is randomized after the event, as described in Supporting Information. We superpose onto time-evolution of the dark mode according to the potential energy in Eq. (8) a stochastic process in which Jd → Jd + ~ with rate constant ku and Jd → Jd − ~ with rate constant kd . The statistical averaging in the response function is generalized from the explicit thermal average over two degrees of freedom in Eq. (1) to include an average over these stochastic trajectories, which are generated by a Kinetic Monte Carlo procedure69 described in Supporting Information. The values of the rate constants ku and kd must be determined externally to this calculation. We take ku = 0.015 ps−1 , with kd set by detailed balance. In CO2 at 298 K, the relative thermal populations of the ground and excited states of one bend mode are approximately 0.96 and 0.04, so that the populations of the ground and first excited state of the degenerate pair of modes are approximately 0.93 and 0.07. To mimic these populations in our two-mode model, we adopt a higher temperature of T = 383 K to double the excited-state thermal population relative to room temperature, giving kd = 0.184 ps−1 . Brinzer, et al.15,16 fit their measured 2D spectra of solvated CO2 with a model in which the chromophore frequency undergoes transitions between two values. Their best fit value16 of ku is approximately a factor of 6 smaller than the value used here. The value of these rate constants set the evolution-time dependence of the 2D spectra, as well as contributing to line shapes of spectral peaks. Our calculated spectra presented below do not include other linebroadening mechanisms present in experimental spectra and so are given nonphysical line widths by application of window functions prior to Fourier transformation. These spectra show time evolution approximately comparable to that measured in Ref. 15 for CO2 . A schematic 2D spectrum for the chromophore mode in Fig. 2(A) illustrates the formation of cross peaks from exchange processes during the t2 interval. The higher amplitude diagonal peak reflects time evolution of the chromophore during t1 and t3 with the dark mode mostly ACS Paragon 9Plus Environment


in its ground state, while the lower amplitude diagonal peak results from chromophore dynamics with primarily the dark mode excited. At t2 = 0, these diagonal peaks dominate the spectrum. All peaks in Fig. 2(A) result from processes in which all spectroscopic transitions are between the ground and first excited state of the chromophore mode. The overtone peaks anharmonically red shifted in ω3 resulting from transitions between first and second excited states of the chromophore are not shown here. With increasing t2 , there is increasing likelihood of thermal transitions between ground and excited populations of the dark mode during this interval, leading to cross peaks. The red arrow indicates transfer of peak amplitude to a cross peak resulting from thermal transitions of the dark mode from ground to first excited state during the t2 interval, producing a red shift in ω3 . Figure 2(B) shows as an example the uppermost OMT diagram from Fig. 1 with the t2 interval colored in red. Below it is a schematic of an example stochastic trajectory of the action of the dark mode, indicating a transition during the t2 interval from Jd = ~/2 to Jd = 3~/2, representing a thermal transition from the ground to the first excited state. Incoherent thermal transitions between ground and excited-state populations are represented by action changes of ±~, in contrast to the spectroscopic transitions involving one application of a raising or lowering operator which are semiclassically represented by action changes of ±~/2. This stochastic trajectory is one example of a process that transfers amplitude as indicated by the red arrow in the spectrum. The lower amplitude diagonal peak in Fig. 2(A) results from processes in which the dark mode is excited for significant portions of the t1 and t3 intervals. The blue arrow in the spectrum shows transfer of amplitude to a cross peak resulting from processes such as that shown in the blue portions of the diagrams in Fig. 2(C). The upper image shows an OMT diagram, and the lower image shows an example of a stochastic trajectory of the action of the dark mode. The dark mode begins the t2 interval in the excited state and finishes in the ground state. An incoherent transition from excited to ground state of the dark mode results in a blue shift of the chromophore frequency in the ω3 dimension. While the OMT method avoids the “sign problem” of highly oscillatory integrands arising in other semiclassical methods,63 some numerical challenges remain. Principal among these is the necessity of calculating the time-domain response function as in Eq. (1) at sufficiently large numbers of values of t1 and t3 to permit the double Fourier transform in Eq. (7).56 If time propagation is carried out naively from the instant of the first radiation-matter interaction, i.e. from left to right in Fig. 1, each choice of a new value of t1 requires ACS Paragon10 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


repeating the calculation of the response function. Alemi and Loring56 implemented a more efficient “forward-backward” approach,70 in which the state of the system is sampled after the second interaction with the radiation field at time t1 . Starting at the second vertex in any diagram in Fig. 1, time is propagated backwards for −t1 and forwards for t2 and t3 . One initial condition then provides the response function at any desired number of t1 and t3 values. The validity of this approach relies on the response function computation being an equilibrium calculation, on Liouville’s theorem, and on the approximate assumption that “good” action variables can be identified, so that applying a semiclassical quantization condition at a later time is equivalent to applying the same condition at an earlier time. The assumption of the existence of effectively good action variables underlies the justification of the rest of the OMT procedure54,55 as well. Calculations presented in Sec. III employ the forward-backward procedure of Ref. 56, modified for the present model as described in Supporting Information. Evaluating OMT diagrams as shown in Fig. 1 requires implementing canonical transformations between Cartesian coordinates and momenta and action-angle variables. Timepropagation is carried out in the former, while quantized action jumps representing interactions with radiation or incoherent thermal transitions require the latter. Previous applications of the OMT54–57,71 have demonstrated that low-order classical mechanical perturbation theory in anharmonicity72,73 suffices to carry out these transformations in this application. Transitions in chromophore action were treated to second order in the cubic anharmonicity of the Morse potential and to first order in quartic anharmonicity of the dark mode and in the biquadratic coupling between modes. Transitions in action of the dark mode were treated to zeroth order in anharmonicity, as described in Supporting Information. This perturbation theory is additionally applied to evaluate the thermal weight ∆F in Eq. (3) as also described in Supporting Information. Classical perturbation theory is only used to perform canonical transformations; all propagation is carried out with the full anharmonic Hamiltonian.

III.

RESULTS AND DISCUSSION

Two-dimensional infrared spectra from Eqs. (1) and (7) for the two-mode model of Eq. (8) are shown in Fig. 3. The stochastic transitions in energy of the dark mode take place ACS Paragon11 Plus Environment


2350

(A)

2350

(B)

(1* ,1)

(1,1) 2340

(1,1)

2340 (1* ,1* ) 2330

ω3

2330

ω3

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 12 of 22

2320

(1,2)

2310

(1,1* ) (1* ,2)

2320

(1,2)

2310

(1* ,2* )

2300

(1,2* )

2300

2300

2310

2320

2330

2340

2350

2300

ω1

2310

2320

2330

2340

2350

ω1

FIG. 3: Two-dimensional spectra for the two-mode model of Eq. (8). In (A), evolution time t2 = 0, and in (B), t2 =55 ps. Each line is labeled with peak frequencies (ω1 , ω3 ). Labels n and n∗ respectively denote frequencies of one-quantum chromophore transitions from n − 1 to n, with the dark mode in the ground or excited states.

throughout the calculation, including during the t1 and t3 intervals, and provide one source of line broadening of the chromophore spectrum. Besides this stochastic process, the model does not include direct solvent-induced dephasing of the chromophore mode that would broaden peaks associated with this mode. To smooth the spectra resulting from a discrete Fourier transform, we employ apodization, as described in the Supporting Information, effectively broadening the lines with nonphysical dephasing. Our window function in both t1 and t3 dimensions cuts off at 7.9 ps, giving line widths approximately comparable to those of the solvated CO2 system of Ref. 15, for which the pure dephasing time is determined to be 2.6 ps. The spectra in Fig. 3 were converged using 75 sets of initial angle variables for the chromophore and dark modes and 9950 stochastic Kinetic Monte Carlo trajectories for the action of the dark mode. To identify spectral features, we label the frequency of the n − 1 → n transition in the chromophore mode by n if the low-frequency dark mode mostly occupies its ground state and by n∗ if the dark mode is primarily singly excited. From the fit of CO2 data15 to Eq. (10), frequencies 1 and 2 are approximately 2340 cm−1 and 2316 cm−1 and frequencies 1∗ and 2∗ are approximately 2328 cm−1 and 2304 cm−1 , respectively. The chromophore anharmonic ACS Paragon12 Plus Environment

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


shift is -24 cm−1 , whether or not the dark mode is excited. In Fig. 3(A), the evolution time t2 = 0. This spectrum shows anharmonically shifted dominant diagonal peaks at (1, 1) and (1, 2), arising from processes with the dark mode in the ground state. The initial interaction with radiation in each case is the fundamental 0→1 transition, with the second interaction probing either the same transition or the overtone 1 → 2. The two peaks are therefore shifted along the ω3 axis by the anharmonicity of the chromophore mode. The analogous pair of diagonal peaks with the dark mode excited occurs near (1∗ , 1∗ ) and (1∗ , 2∗ ). For the absorptive 2D spectrum as defined in Eq. (7), anharmonically shifted peak amplitudes are of different sign, as indicated by red (positive amplitude) and blue (negative amplitude) in Fig. 3. Figure 3(B) shows the spectrum at t2 = 55 ps. As the model does not contain population relaxation of the chromophore mode, this spectrum resembles the asymptotic long t2 limit for this calculation. For the solvated CO2 system of Ref. 17, the vibrational lifetime of the asymmetric stretch mode has been measured to be 64 ps, longer than time scales of population dynamics of the bend mode. The diagonal peaks in Fig. 3(B) have diminished in amplitude and four cross peaks have appeared through thermal transitions in populations of the dark mode, like those indicated by blue and red arrows in Fig. 2. The diagonal peak at (1,1) has lost amplitude to a cross peak at (1,1∗ ) and the anharmonically shifted diagonal peak at (1,2) has transferred amplitude to (1,2∗ ). The diagonal peak at (1∗ ,1∗ ) in Fig. 3(A) is no longer visible in Fig. 3(B), as amplitude has transferred to the cross peak at (1∗ ,1). Similarly, the anharmonically shifted peak at (1∗ ,2∗ ) in Fig. 3(A) is no longer apparent in Fig. 3(B), because of the growth of the cross peak at (1∗ ,2). The 2D spectra in Fig. 3 at short and long times agree qualitatively with measured 2D spectra of the asymmetric stretch mode of solvated CO2 .15,17 The t2 dependence of 2D spectra as in Fig. 3 arises from incoherent excitation and deexcitation of the dark mode with respective rate constants ku and kd as specified in Sec. II. The diagonal peaks should thus decay and cross peaks should grow in with increasing t2 with the sum of these rate constants.15,17 The decay of the diagonal peak at (1∗ ,1∗ ) and the associated growth of the cross peak at (1∗ ,1) arises from the deexcitation of the excited dark mode during the t2 interval, resulting in a blue shift in ω3 . Peak volumes calculated from 2D spectra as described in Supporting Information are shown for these features in Fig. 4. Both peak volumes are scaled by the volume at t2 = 0 of the dominant fundamental peak ACS Paragon13 Plus Environment


0.08

▲ ▲ ▲ ▲ ▲ ▲ ▲ ▲ ▲ ▲ ▲ ▲ (1* ,1)

●

Peak Volume

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

▲

0.06

▲ ●

0.04

0.02

●

▲ 0.00

Page 14 of 22

0

(1* ,1* )

● ●

10

●

20

●

●

●

30

●

●

●

40

●

●

●

50

t2 (ps)

FIG. 4: Dependence of peak volumes on evolution time t2 . The diagonal peak at (1∗ ,1∗ ) reflects processes with spectroscopic transitions between ground and first excited states of the chromophore with the dark mode primarily in its excited state during t1 and t3 intervals, while the cross peak at (1∗ ,1) arises from deexcitation of the dark mode during t2 , resulting in a blue shift in ω3 . The volume of the diagonal peak (circles) decays as the volume of the cross peak (triangles) grows. Curves show fits to exponential time-dependence with decay rate equal to the sum of rate constants for excitation and deexcitation.

(1,1). Therefore the t2 = 0 value of the scaled (1∗ ,1∗ ) peak volume is approximately equal to the ratio of the equilibrium populations of excited and ground states for the dark mode. Curves show fits to exponential time-dependence with decay rate equal to ku + kd . The dependence of peak amplitudes on t2 is shown to follow the thermal population dynamics of the dark mode. The process (1∗ ,2∗ )→ (1∗ ,2) similarly arises from radiationless deexcitation of the dark mode, and the associated peak volume plots are very similar to those in Fig. 4. The processes (1,1)→ (1,1∗ ) and (1,2)→ (1,2∗ ) result from thermal excitation of the dark mode. Peak volume plots for these features follow the same exponential kinetics as in Fig. 4, but are significantly noisier for the same number of trajectories. This arises in part because the 2D spectrum for this model with initial frequency 1 is more congested than that with initial frequency 1∗ , as can be seen from Fig. 3(B), complicating the definition of peak volumes. The convergence of the calculation is also impeded by changes in line shape resulting from processes with multiple stochastic transitions in population of the dark mode during intervals t1 and t3 .

ACS Paragon14 Plus Environment

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


IV.

CONCLUSIONS

The semiclassical OMT approach has been shown to reproduce a range of phenomenology associated with 2D vibrational spectra of coupled anharmonic vibrations, including pure dephasing,56 energy relaxation,56 and coherence transfer.57 The calculations shown in Figs. 3 and 4 demonstrate that the method can capture the spectrum of a high frequency mode reporting on thermally driven population fluctuations in a low frequency mode, as observed15–17 in solvated CO2 . Both the thermal population dynamics and the radiationinduced quantum state changes are represented here by changes in values of classical action variables. The thermal dynamics are mimicked by a stochastic process. The method should be practically applicable to a completely mechanical model in which a bath of explicitly treated degrees of freedom induces classical mechanical energy fluctuations in a dark mode, thereby influencing a chromophore mode; this implementation remains for future work. The calculations presented here provide a distinct example of this classical-trajectory-based approach describing the phenomenology of multidimensional vibrational spectra of solvated anharmonic vibrations.

SUPPORTING INFORMATION

Classical-mechanical action-angle perturbation theory; Kinetic Monte Carlo implementation; apodization procedure

ACKNOWLEDGMENTS

This material is based upon work supported by the National Science Foundation under CHE-1361484.

REFERENCES 1

Woutersen, S.; Mu, Y.; Stock, G.; Hamm, P. Subpicosecond Conformational Dynamics of Small Peptides Probed by Two-Dimensional Vibrational Spectroscopy . Proc. Natl. Acad. Sci. US 2001, 98, 11254–11258. ACS Paragon15 Plus Environment


2

Wong, D. B.; Giammanco, C. H.; Fenn, E. E.; Fayer, M. D. Dynamics of Isolated Water Molecules in a Sea of Ions in a Room Temperature Ionic Liquid. J. Phys. Chem. B 2013, 117, 623–635.

3

Kwak, K.; Zheng, J.; Cang, H.; Fayer, M. D. Ultrafast Two-Dimensional Infrared Vibrational Echo Chemical Exchange Experiments and Theory. J. Phys. Chem. B 2006, 110, 19998–20013.

4

Kwac, K.; Lee, C.; Jung, Y.; Han, J.; Kwak, K.; Zheng, J. R.; Fayer, M. D.; Cho, M. PhenolBenzene Complexation Dynamics: Quantum Chemistry Calculation, Molecular Dynamics Simulations, and Two Dimensional IR Spectroscopy . J. Chem. Phys. 2006, 125, 244508.

5

Jansen, T. L. C.; Knoester, J. Calculation of Two-Dimensional Infrared Spectra of Ultrafast Chemical Exchange with Numerical Langevin Simulations . J. Chem. Phys. 2007, 127, 234502.

6

Fayer, M. D. Dynamics of Liquids, Molecules, and Proteins Measured with Ultrafast 2D IR Vibrational Echo Chemical Exchange Spectroscopy . Annu. Rev. Phys. Chem. 2009, 60, 21–38.

7

Bagchi, S.; Charnley, A. K.; Smith, A. B.; Hochstrasser, R. M. Equilibrium Exchange Processes of the Aqueous Tryptophan Dipeptide . J. Phys. Chem. B 2009, 24, 8412–8417.

8

Bagchi, S.; Thorpe, D. G.; Thorpe, I. F.; Voth, G. A.; Fayer, M. D. Conformational Switching between Protein Substates Studied with 2D IR Vibrational Echo Spectroscopy and Molecular Dynamics Simulations. J. Phys. Chem. B 2010, 114, 17187–17193.

9

Ji, M. B.; Gaffney, K. J.

Orientational Relaxation Dynamics in Aqueous Ionic Solution:

Polarization-Selective Two-Dimensional Infrared Study of Angular Jump-Exchange Dynamics in Aqueous 6M NaClO4 . J. Chem. Phys. 2011, 134, 044516. 10

Lee, K. K.; Park, K. H.; Joo, C.; Kwon, H. J.; Ha, J. H.; Park, S.; Cho, M. Ultrafast Internal Rotational Dynamics of the Azido Group in (4S)-Azidoproline: Chemical Exchange 2DIR Spectroscopic Investigations. Chem. Phys. 2012, 396, 23–29.

11

Thielges, M. C.; Fayer, M. D. Protein Dynamics Studied with Ultrafast Two-Dimensional Infrared Vibrational Echo Spectroscopy . Acc. Chem. Res. 2012, 45, 1866–1874.

12

Choi, J.-H.; Kwak, K.-W.; Cho, M. Computational Infrared and Two-Dimensional Infrared Photon Echo Spectroscopy of Both Wild-Type and Double Mutant Myoglobin-CO Proteins . J. Phys. Chem. B 2013, 117, 15462–15478.

13

Borek, J. A.; Perakis, F.; Hamm, P. Testing for Memory-Free Spectroscopic Coordinates by 3D IR Exchange Spectroscopy . Proc. Natl. Acad. Sci. U. S. A. 2014, 111, 10462–10467.

14

Feng, C.-J.; Tokmakoff, A. The Dynamics of Peptide-Water Interactions in Dialanine: an UlACS Paragon16 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


trafast Amide I 2D IR and Computational Spectroscopy Study . J. Chem. Phys. 2017, 147, 085101. 15

Brinzer, T.; Berquist, E. J.; Ren, Z.; Dutta, S.; Johnson, C. A.; Krisher, C. S.; Lambrecht, D. S.; Garrett-Roe, S. Ultrafast Vibrational Spectroscopy (2D-IR) of CO2 in Ionic Liquids: Carbon Capture from Carbon Dioxide’s Point of View. J. Chem. Phys. 2015, 142, 212425.

16

Brinzer, T.; Berquist, E. J.; Ren, Z.; Dutta, S.; Johnson, C. A.; Krisher, C. S.; Lambrecht, D. S.; Garrett-Roe, S. Erratum: “Ultrafast Vibrational Spectroscopy (2D-IR) of CO2 in Ionic Liquids: Carbon Capture from Carbon Dioxide’s Point of View”[J. Chem. Phys. 142, 212425 (2015)]. J. Chem. Phys. 2017, 147, 049901.

17

Giammanco, C. H.; Kramer, P. L.; Yamada, S. A.; Nishida, J.; Tamimi, A.; Fayer, M. D. Coupling of Carbon Dioxide Stretch and Bend Vibrations Reveals Thermal Population Dynamics in an Ionic Liquid. J. Phys. Chem. B 2016, 120, 549–556.

18

Giammanco, C. H.; Kramer, P. L.; Yamada, S. A.; Nishida, J.; Tamimi, A.; Fayer, M. D. Carbon Dioxide in an Ionic Liquid: Structural and Rotational Dynamics. J. Chem. Phys. 2016, 144, 104506.

19

Giammanco, C. H.;

Yamada, S. A.;

Kramer, P. L.;

Tamimi, A.;

Fayer, M. D.

Structural and Rotational Dynamics of Carbon Dioxide in 1-alkyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide Ionic Liquids: the Effect of Chain Length. J. Phys. Chem. B 2016, 120, 6698–6711. 20

Daly, Jr., C. A.; Berquist, E. J.; Brinzer, T.; Garrett-Roe, S.; Lambrecht, D. S.; Corcelli, S. A. Modeling Carbon Dioxide Vibrational Frequencies in Ionic Liquids: II. Spectroscopic Map. J. Phys. Chem. B 2016, 120, 12633–12642.

21

Berquist, E. J.; Daly, Jr., C. A.; Brinzer, T.; Bullard, K. K.; Campbell, Z. M.; Corcelli, S. A.; Garrett-Roe, S.; Lambrecht, D. S. Modeling Carbon Dioxide Vibrational Frequencies in Ionic Liquids: I. Ab Initio Calculations. J. Phys. Chem. B 2017, 121, 208–220.

22

Shin, J. Y.; Yamada, S. A.; Fayer, M. D. Carbon Dioxide in a Supported Ionic Liquid Membrane: Structural and Rotational Dynamics Measured with 2D IR and Pump-Probe Experiments. J. Am. Chem. Soc. 2017, 139, 11222–11232.

23

Cho, M. Coherent Two-Dimensional Optical Spectroscopy. Chem. Rev. 2008, 108, 1331–1418.

24

Zhuang, W.; Hayashi, T.; Mukamel, S. Coherent Multidimensional Vibrational Spectroscopy of Biomolecules: Concepts, Simulations, and Challenges. Chem. Rev. 2009, 48, 3750–3781. ACS Paragon17 Plus Environment


25

Mukamel, S.; Khidekel, V.; Chernyak, V. Classical Chaos and Fluctuation-Dissipation Relations for Nonlinear Response. Phys. Rev. E 1996, 53, R1–R4.

26

Saito, S.; Ohmine, I. Off-Resonant Fifth-Order Response Function for Two-Dimensional Raman Spectroscopy of Liquids CS2 and H2 O. Phys. Rev. Lett. 2002, 88, 207401.

27

Kryvohuz, M.; Cao, J. Classical Divergence of Nonlinear Response Functions. Phys. Rev. Lett. 2006, 96, 030403.

28

Malinin, S. V.; Chernyak, V. Y. Collective Oscillations in the Classical Nonlinear Response of a Chaotic System. Phys. Rev. E 2008, 77, 025201(R).

29

Jeon, J.; Cho, M. Direct Quantum Mechanical/Molecular Mechanical Simulations of TwoDimensional Vibrational Responses: N-Methylacetamide in Water. New. J. Phys. 2010, 12, 065001.

30

Sakurai, A.; Tanimura, Y. Does ~ Play a Role in Multidimensional Spectroscopy? Reduced Hierarchy Equations of Motion Approach to Molecular Vibrations. J. Phys. Chem. A 2011, 115, 4009–4022.

31

Ito, H.; Hasegawa, T.; Tanimura, Y. Calculating Two-Dimensional THz-Raman-THz and Raman-THz-THz Signals for Various Molecular Liquids. J. Chem. Phys. 2014, 141, 124503.

32

Ishizaki, A.; Tanimura, Y. Modeling Vibrational Dephasing and Energy Relaxation of Intramolecular Anharmonic Modes for Multidimensional Infrared Spectroscopies. J. Chem. Phys. 2006, 125, 084501.

33

Jeon, J.; Yang, S.; Choi, J.-H.; Cho, M. Computational Vibrational Spectroscopy of Peptides and Proteins in One and Two Dimensions. Acc. Chem. Res. 2009, 42, 1280–1289.

34

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.

35

Jansen, T. L. C.; Knoester, J. Waiting Time Dynamics in Two-Dimensional Infrared Spectroscopy. Acc. Chem. Res. 2009, 42, 1405–1411.

36

Shi, Q.; Geva, E. A Comparison Between Different Semiclassical Approximations for Optical Response Functions in Nonpolar Liquid Solutions. J. Chem. Phys. 2005, 122, 064506.

37

McRobbie, P. L.; Geva, E. A Benchmark Study of Different Methods for Calculating One- and Two-Dimensional Optical Spectra. J. Phys. Chem. A 2009, 113, 10425–10434.

38

Hanna, G.; Geva, E. Multidimensional Spectra via the Mixed Quantum-Classical Liouville ACS Paragon18 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


Method: Signatures of Nonequilibrium Dynamics. J. Phys. Chem. B 2009, 113, 9278–9288. 39

Hahn, S.; Ham, S.; Cho, M. Simulation Studies of Amide I IR Absorption and Two-Dimensional IR Spectra of Beta Hairpins in Liquid Water. J. Phys. Chem. B 2005, 109, 11789–11801.

40

Hayashi, T.; Mukamel, S. Vibrational-Exciton Couplings for the Amide I, II, III, and A Modes of Peptides. J. Phys. Chem. B 2007, 111, 11032–11046.

41

Woys, A. M.; Almeida, A. M.; Wang, L.; Chiu, C. C.; McGovern, M.; de Pablo, J. J.; Skinner, J. L.; Gellman, S. H.; Zanni, M. T. Parallel Beta-Sheet Vibrational Couplings Revealed by 2D IR Spectroscopy of an Isotopically Labeled Macrocycle: Quantitative Benchmark for the Interpretation of Amyloid and Protein Infrared Spectra. J. Am. Chem. Soc. 2012, 134, 19118–19128.

42

Schmidt, J. R.; Corcelli, S. A.; Skinner, J. L. Ultrafast Vibrational Spectroscopy of Water and Aqueous N-Methylacetamide: Comparison of Different Electronic Structure/Molecular Dynamics Approaches. J. Chem. Phys. 2004, 121, 8887–8896.

43

Corcelli, S. A.; Skinner, J. L. Infrared and Raman Line Shapes of Dilute HOD in Liquid H2 O and D2 O from 10 to 90◦ C. J. Phys. Chem. A 2005, 109, 6154–6165.

44

Li, F.; Skinner, J. L. Infrared and Raman Line Shapes for Ice Ih. I. Dilute HOD in H2 O and D2 O. J. Chem. Phys. 2010, 132, 204505.

45

Kwac, K.; Cho, M. H. Molecular Dynamics Simulation Study of N-Methylacetamide in Water. I. Amide I Mode Frequency Fluctuation. J. Chem. Phys. 2003, 119, 2247–2255.

46

Kwac, K.; Cho, M. H. Molecular Dynamics Simulation Study of N-Methylacetamide in Water. II. Two-Dimensional Infrared Pump Probe Spectra. J. Chem. Phys. 2003, 119, 2256–2263.

47

Hayashi, T.; Jansen, T. L. C.; Zhuang, W.; Mukamel, S. Collective Solvent Coordinates for the Infrared Spectrum of HOD in D2 O Based on an Ab Initio Electrostatic Map. J. Phys. Chem. A 2005, 109, 64–82.

48

Skinner, J. L.; Pieniazek, P. A.; Gruenbaum, S. M. Vibrational Spectroscopy of Water at Interfaces. Acc. Chem. Res. 2012, 45, 93–100.

49

Thoss, M.; Miller, W. H. Generalized Forward-Backward Initial Value Representation for the Calculation of Correlation Functions in Complex Systems. J. Chem. Phys. 2001, 114, 9220– 9235.

50

Wu, J.; Cao, J. Linear and Nonlinear Response Functions of the Morse Oscillator: Classical Divergence and the Uncertainty Principle. J. Chem. Phys. 2001, 115, 5381–5391. ACS Paragon19 Plus Environment


51

Makri, N. Forward-Backward Semiclassical and Quantum Trajectory Methods for Time Correlation Functions. Phys. Chem. Chem. Phys. 2011, 13, 14442–14452.

52

Grossmann, F. A Semiclassical Hybrid Approach to Linear Response Functions for Infrared Spectroscopy. Phys. Scr. 2016, 91, 044004.

53

Buchholz, M.; Grossmann, F.; Ceotto, M. Application of the Mixed Time-Averaging Semiclassical Initial Value Representation Method to Complex Molecular Spectra. J. Chem. Phys. 2017, 147, 164110.

54

Gerace, M.; Loring, R. F. An Optimized Semiclassical Approximation for Vibrational Response Functions. J. Chem. Phys. 2013, 138, 124104.

55

Gerace, M.; Loring, R. F. Two-Dimensional Spectroscopy of Coupled Vibrations with the Optimized Mean-Trajectory Approximation. J. Phys. Chem. B 2013, 117, 15452–15461.

56

Alemi, M.; Loring, R. F. Two-Dimensional Vibrational Spectroscopy of a Dissipative System with the Optimized Mean-Trajectory Approximation. J. Phys. Chem. B 2015, 119, 8950–8959.

57

Alemi, M.; Loring, R. F. Vibrational Coherence and Energy Transfer in Two-Dimensional Spectra with the Optimized Mean-Trajectory Approximation. J. Chem. Phys. 2015, 142, 212417.

58

Moberg, D. R.; Alemi, M.; Loring, R. F. Thermal weights for semiclassical response functions. J. Chem. Phys. 2015, 143, 084101.

59

Hamm, P.; Zanni, M. Concepts and Methods of 2D Infrared Spectroscopy; Cambridge University Press: New York, 2011.

60

Yee, T. K.; Gustafson, T. K. Diagrammatic Analysis of the Density Operator for Nonlinear Optical Calculations: Pulsed and CW Responses. Phys. Rev. A 1978, 18, 1597–1617.

61

Mukamel, S. Principles of Nonlinear Optical Spectroscopy; Oxford University Press: New York, 1995.

62

Gruenbaum, S. M.; Loring, R. F. Interference and Quantization in Semiclassical Response Functions. J. Chem. Phys. 2008, 128, 124106.

63

Noid, W. G.; Ezra, G. S.; Loring, R. F. Semiclassical Calculation of the Vibrational Echo. J. Chem. Phys. 2004, 120, 1491–1499.

64

Loring, R. F. Mean-Trajectory Approximation for Electronic and Vibrational-Electronic Nonlinear Spectroscopy. J. Chem. Phys. 2017, 146, 144106.

65

Allen, W. D.; Yamaguchi, Y.; Csaszar, A. G.; Clabo, Jr., D. A.; Remington, R. B.; Schaefer III, H. F. A Systematic Study of Molecular Vibrational Anharmonicity and Vibration-Rotation ACS Paragon20 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


Interaction by Self-Consistent-Field Higher-Derivative Methods. Linear Polyatomic Molecules. Chem. Phys. 1990, 145, 427–466. 66

Teffo, J.-L.; Sulakshina, O. N.; Perevalov, V. I. Effective Hamiltonian for Rovibrational Energies and Line Intensities of Carbon Dioxide. J. Mol. Spectry. 1992, 156, 48–64.

67

Martin, J. M. L.; Taylor, P. R.; Lee, T. J. Accurate Ab Initio Quartic Force Fields for the N2 O and CO2 Molecules. Chem. Phys. Lett. 1993, 205, 535–542.

68

Dressler, S.; Thiel, W. Anharmonic Force Fields from Density Functional Theory. Chem. Phys. Lett. 1993, 205, 535–542.

69

Bortz, A.; Kalos, M.; Lebowitz, J. A New Algorithm for Monte Carlo Simulation of Ising Spin Systems. J. Comput. Phys. 1975, 17, 10–18.

70

Hasegawa, T.; Tanimura, Y. Nonequilibrium Molecular Dynamics Simulations with a BackwardForward Trajectories Sampling for Multidimensional Infrared Spectroscopy of Molecular Vibrational Modes. J. Chem. Phys. 2008, 128, 064511.

71

Gruenbaum, S. M.; Loring, R. F. Semiclassical Mean-Trajectory Approximation for Nonlinear Spectroscopic Response Functions. J. Chem. Phys. 2008, 129, 124510.

72

Schatz, G. C.; Mulloney, T. Classical Perturbation Theory of Good Action Angle Variables. Application to Semiclassical Eigenvalues and to Collisional Energy Transfer in Polyatomic Molecules. J. Phys. Chem. 1979, 83, 989–999.

73

Gruenbaum, S. M. A Practical Method for Understanding and Evaluating Semiclassical Vibrational Response Functions. Ph.D. thesis, Cornell University, 2010.



TOC Graphic


Page 22 of 22

Thermal Population Fluctuations in Two-Dimensional Infrared

Recommend Documents