Shaul Mukamel

Using coherence to enhance function in chemical and biophysical systems

Coherence phenomena arise from interference, or the addition, of wave-like amplitudes with fixed phase differences. Although coherence has been shown to yield transformative ways for improving function, advances have been confined to pristine matter and coherence was considered fragile. However, recent evidence of coherence in chemical and biological systems suggests that the phenomena are robust and can survive in the face of disorder and noise. Here we survey the state of recent discoveries, present viewpoints that suggest that coherence can be used in complex chemical systems, and discuss the role of coherence as a design element in realizing function.

A Quantum Chemical Interpretation of Two-Dimensional Electronic Spectroscopy of Light-Harvesting Complexes

Nonlinear electronic spectroscopies represent one of the most powerful techniques to study complex multichromophoric architectures. For these systems, in fact, linear spectra are too congested to be used to disentangle the many coupled vibroelectronic processes that are activated. By using a 2D approach, instead, a clear picture can be achieved, but only when the recorded spectra are combined with a proper interpretative model. So far, this has been almost always achieved through parametrized exciton Hamiltonians that necessarily introduce biases and/or arbitrary assumptions. In this study, a first-principles approach is presented that combines accurate quantum chemical descriptions with state-of-the-art models for the environment through the use of atomistic and polarizable embeddings. Slow and fast bath dynamics, along with exciton transport between the pigments, are included. This approach is applied to the 2DES spectroscopy of the Light-Harvesting 2 (LH2) complex of purple bacteria. Simulations are extended over the entire visible-near-infrared spectral region to cover both carotenoid and bacteriochlorophyll signals. Our results provide an accurate description of excitonic properties and relaxation pathways, and give an unprecedented insight into the interpretation of the spectral signatures of the measured 2D signals.

UV-Light-Induced Vibrational Coherences: The Key to Understand Kasha Rule Violation in trans-Azobenzene

We combine sub-20 fs transient absorption spectroscopy with state-of-the-art computations to study the ultrafast photoinduced dynamics of trans-azobenzene (AB). We are able to resolve the lifetime of the ππ* state, whose decay within ca. 50 fs is correlated to the buildup of the nπ* population and to the emergence of coherences in the dynamics, to date unobserved. Nonlinear spectroscopy simulations call for the CNN in-plane bendings as the active modes in the subps photoinduced coherent dynamics out of the ππ* state. Radiative to kinetic energy transfer into these modes drives the system to a high-energy planar nπ*/ground state conical intersection, inaccessible upon direct excitation of the nπ* state, that triggers an ultrafast (0.45 ps) nonproductive decay of the nπ* state and is thus responsible for the observed Kasha rule violation in UV excited trans-AB. On the other hand, cis-AB is built only after intramolecular vibrational energy redistribution and population of the NN torsional mode.

Monitoring molecular nonadiabatic dynamics with femtosecond X-ray diffraction

Significance X-ray crystallography has long been used to determine the structure of crystals and molecular samples. More recent advancements in light sources and computational methods made it possible to routinely determine the structure of large proteins. The introduction of X-ray free-electron lasers opens up the possibility to track the dynamics of molecular structures on a femtosecond time scale and to create molecular movies of chemical reactions. The theory of time-independent diffraction is well known. However, time-resolved diffraction techniques pose not only new challenges to experiments but also to their interpretation. In this work, we present a unified theoretical framework that will aid experimental interpretations as well as predictions of types of X-ray diffraction experiments. Ultrafast time-resolved X-ray scattering, made possible by free-electron laser sources, provides a wealth of information about electronic and nuclear dynamical processes in molecules. The technique provides stroboscopic snapshots of the time-dependent electronic charge density traditionally used in structure determination and reflects the interplay of elastic and inelastic processes, nonadiabatic dynamics, and electronic populations and coherences. The various contributions to ultrafast off-resonant diffraction from populations and coherences of molecules in crystals, in the gas phase, or from single molecules are surveyed for core-resonant and off-resonant diffraction. Single-molecule ∝N scaling and two-molecule ∝N2 scaling contributions, where N is the number of active molecules, are compared. Simulations are presented for the excited-state nonadiabatic dynamics of the electron harpooning at the avoided crossing in NaF. We show how a class of multiple diffraction signals from a single molecule can reveal charge-density fluctuations through multidimensional correlation functions of the charge density.

Utilizing Microcavities To Suppress Third-Order Cascades in Fifth-Order Raman Spectra

Nonlinear optical signals in the condensed phase are often accompanied by sequences of lower-order processes, known as cascades, which share the same phase matching and power dependence on the incoming fields and are thus hard to distinguish. The suppression of cascading in order to reveal the desired nonlinear signal has been a major challenge in multidimensional Raman spectroscopy, that is, the χ(5) signal being masked by cascading signals given by a product of two χ(3) processes. Because cascading originates from the exchange of a virtual photon between molecules, it can be manipulated by performing the experiment in an optical microcavity which modifies the density of radiation field modes. Using a quantum electrodynamical (QED) treatment, we demonstrate that the χ(3) cascading contributions can be greatly suppressed. By optimizing the cavity size and the incoming pulse directions, we show that up to ∼99.5% suppression of the cascading signal is possible.

Towards Accurate Simulation of Two-Dimensional Electronic Spectroscopy

We introduce the basic concepts of two-dimensional electronic spectroscopy (2DES) and a general theoretical framework adopted to calculate, from first principles, the nonlinear response of multi-chromophoric systems in realistic environments. Specifically, we focus on UV-active chromophores representing the building blocks of biological systems, from proteins to nucleic acids, describing our progress in developing computational tools and protocols for accurate simulation of their 2DUV spectra. The roadmap for accurate 2DUV spectroscopy simulations is illustrated starting with benchmarking of the excited-state manifold of the chromophoric units in a vacuum, which can be used for building exciton Hamiltonians for large-scale applications or as a reference for first-principles simulations with reduced computational cost, enabling treatment of minimal (still realistic) multi-chromophoric model systems. By adopting a static approximation that neglects dynamic processes such as spectral diffusion and population transfer, we show how 2DUV is able to characterize the ground-state conformational space of dinucleosides and small peptides comprising dimeric chromophoric units (in their native environment) by tracking inter-chromophoric electronic couplings. Recovering the excited-state coherent vibrational dynamics and population transfers, we observe a remarkable agreement between the predicted 2DUV spectra of the pyrene molecule and the experimental results. These results further led to theoretical studies of the excited-state dynamics in a solvated dinucleoside system, showing that spectroscopic fingerprints of long-lived excited-state minima along the complex photoinduced decay pathways of DNA/RNA model systems can be simulated at a reasonable computational cost. Our results exemplify the impact of accurate simulation of 2DES spectra in revealing complex physicochemical properties of fundamental biological systems and should trigger further theoretical developments as well as new experiments.

Attosecond X-ray Diffraction Triggered by Core or Valence Ionization of a Dipeptide

With the advancement of intense ultrafast X-ray sources, it is now possible to create a molecular movie by following the electronic dynamics in real time and real space through time-resolved X-ray diffraction. Here we employ real-time time-dependent density functional theory (RT-TDDFT) to simulate the electronic dynamics after an impulse core or valence ionization in the glycine-phenylalanine (GF) dipeptide. The time-evolving dipole moment, the charge density, and the time-resolved X-ray diffraction signals are calculated. The charge oscillation is calculated for 7 fs for valence ionization and 500 as for core ionization. The charge oscillation time scale is comparable to that found in a phenylalanine monomer (4 fs) [ Science 2014 , 346 , 336 ] and is slightly longer because of the elongated glycine chain. Following valence ionization, the charge migration across the GF is mediated by the delocalized lone-pair orbitals of oxygen and nitrogen of the electron-rich amide group. The temporal Fourier transform of the dipole moment provides detailed information on the charge migration dynamics and the molecular orbitals involved. Heterodyne-detected attosecond X-ray diffraction signals provide the magnitude and phase of the scattering amplitude in momentum space and can thus be inverted to yield the charge density in real space.

Linear and nonlinear frequency- and time-domain spectroscopy with multiple frequency combs

Two techniques that employ equally spaced trains of optical pulses to map an optical high frequency into a low frequency modulation of the signal that can be detected in real time are compared. The development of phase-stable optical frequency combs has opened up new avenues to metrology and spectroscopy. The ability to generate a series of frequency spikes with precisely controlled separation permits a fast, highly accurate sampling of the material response. Recently, pairs of frequency combs with slightly different repetition rates have been utilized to down-convert material susceptibilities from the optical to microwave regime where they can be recorded in real time. We show how this one-dimensional dual comb technique can be extended to multiple dimensions by using several combs. We demonstrate how nonlinear susceptibilities can be quickly acquired using this technique. In a second class of techniques, sequences of ultrafast mode locked laser pulses are used to recover pathways of interactions contributing to nonlinear susceptibilities by using a photo-acoustic modulation varying along the sequences. We show that these techniques can be viewed as a time-domain analog of the multiple frequency comb scheme.

The highly excited-state manifold of guanine: calibration for nonlinear electronic spectroscopy simulations

A computational protocol based on the complete and restricted active space self-consistent field (CASSCF/RASSCF) methods and their second-order perturbation theory extensions (CASPT2/RASPT2) is employed to benchmark the highly excited-state manifold of the DNA/RNA canonical purine nucleobase guanine in vacuo. Several RASPT2 schemes are tested, displaying a steady convergence of electronic transition energies and dipole moments upon active space enlargement toward the reference values. The outcome allows calibrating and optimizing computational efforts by considering cheaper and more approximate RAS schemes that could enable the characterization of the excited-state manifolds of multi-chromophoric systems, such as DNA/RNA nucleobase dimers or multimers. Simulations of two-dimensional electronic spectra show similar trends to those observed on the other purine nucleobase adenine, deviating from this and other pyrimidine nucleobases in featuring its main excited-state absorption signal, embodied by sizable double HOMO to LUMO excitation contributions, in the UV probing window.

Reply to Stirnemann et al.: Frame retardation is the key reason behind the general slowdown of water reorientation dynamics in concentrated electrolytes

Stirnemann et al. (1) argue that the general deceleration of water reorientation dynamics in the concentrated electrolytes is not entirely due to the frame retardation, since suppressing frame reorientation entirely leads to a maximum slowdown of 1 τ r e o r j u m p + 1 τ r e o r f r a m e 1 τ r e o r j u m p ≈ 1.6 , smaller than the deceleration in many concentrated electrolytes. This same argument, however, can be used to reason that the general deceleration does not entirely come from the retarded jump either, since the decelerations in many concentrated electrolytes are larger than 1 τ r e o r j u m p + 1 τ r e o r f r a m e 1 τ r e o … nn[↵][1]1To whom correspondence may be addressed. Email: smukamel{at}uci.edu or wzhuang{at}fjirsm.ac.cn.nn [1]: #xref-corresp-1-1