Simulating molecular polaritons in the collective regime using few-molecule models.

The study of molecular polaritons beyond simple quantum emitter ensemble models (e.g., Tavis-Cummings) is challenging due to the large dimensionality of these systems and the complex interplay of molecular electronic and nuclear degrees of freedom. This complexity constrains existing models to either coarse-grain the rich physics and chemistry of the molecular degrees of freedom or artificially limit the description to a small number of molecules. In this work, we exploit permutational symmetries to drastically reduce the computational cost of ab initio quantum dynamics simulations for large N. Furthermore, we discover an emergent hierarchy of timescales present in these systems, that justifies the use of an effective single molecule to approximately capture the dynamics of the entire ensemble, an approximation that becomes exact as N → ∞. We also systematically derive finite N corrections to the dynamics and show that addition of k extra effective molecules is enough to account for phenomena whose rates scale as 𝒪(N-k). Based on this result, we discuss how to seamlessly modify existing single-molecule strong coupling models to describe the dynamics of the corresponding ensemble. We call this approach collective dynamics using truncated equations (CUT-E), benchmark it against well-known results of polariton relaxation rates, and apply it to describe a universal cavity-assisted energy funneling mechanism between different molecular species. Beyond being a computationally efficient tool, this formalism provides an intuitive picture for understanding the role of bright and dark states in chemical reactivity, necessary to generate robust strategies for polariton chemistry.

Physicochemical Principles of AC Electrosynthesis: Reversible Reactions.

Electrolysis integrates renewable energy into chemical manufacturing and is key to sustainable chemistry. Controlling the waveform beyond direct current (DC) addresses the long-standing obstacle of chemoselectivity, yet it also expands the parameter set to optimize, creating a demand for theoretical predictions. Here, we report the first analytical theory for predicting chemoselectivity in an alternating current (AC) electrosynthesis. The mechanism is a selective reversal of the unwanted redox reaction during periods of opposite polarity, reflected in the final reaction outcome as a time-averaged effect. In the ideal scenario of all redox reactions being reversible, square AC waveform biases the outcome toward more overoxidation/over-reduction, whereas sine AC waveform exhibits the opposite effect. However, in a more realistic scenario of some redox reactions being quasi-reversible, sine AC may behave mostly like square AC. These predictions are in numerical agreement with model experiments employing acetophenone and align qualitatively with the literature precedent. Collectively, this study provides theoretical proof for a growing trend that promotes changing waveforms to overcome limitations challenging to address by varying canonical electrochemical parameters.

Alternant Hydrocarbon Diradicals as Optically Addressable Molecular Qubits.

High-spin molecules allow for bottom-up qubit design and are promising platforms for magnetic sensing and quantum information science. Optical addressability of molecular electron spins has also been proposed in first-row transition-metal complexes via optically detected magnetic resonance (ODMR) mechanisms analogous to the diamond-nitrogen-vacancy color center. However, significantly less progress has been made on the front of metal-free molecules, which can deliver lower costs and milder environmental impacts. At present, most luminescent open-shell organic molecules are π-diradicals, but such systems often suffer from poor ground-state open-shell characters necessary to realize a stable ground-state molecular qubit. In this work, we use alternancy symmetry to selectively minimize radical-radical interactions in the ground state, generating π-systems with high diradical characters. We call them m-dimers, referencing the need to covalently link two benzylic radicals at their meta carbon atoms for the desired symmetry. Through a detailed electronic structure analysis, we find that the excited states of alternant hydrocarbon m-diradicals contain important symmetries that can be used to construct ODMR mechanisms leading to ground-state spin polarization. The molecular parameters are set in the context of a tris(2,4,6-trichlorophenyl)methyl (TTM) radical dimer covalently tethered at the meta position, demonstrating the feasibility of alternant m-diradicals as molecular color centers.

Luminescent Organic Triplet Diradicals as Optically Addressable Molecular Qubits.

Optical-spin interfaces that enable the photoinitialization, coherent microwave manipulation, and optical read-out of ground state spins have been studied extensively in solid-state defects such as diamond nitrogen vacancy (NV) centers and are promising for quantum information science applications. Molecular quantum bits (qubits) offer many advantages over solid-state spin centers through synthetic control of their optical and spin properties and their scalability into well-defined multiqubit arrays. In this work, we report an optical-spin interface in an organic molecular qubit consisting of two luminescent tris(2,4,6-trichlorophenyl)methyl (TTM) radicals connected via the meta-positions of a phenyl linker. The triplet ground state of this system can be photoinitialized in its |T0⟩ state by shelving triplet populations as singlets through spin-selective excited-state intersystem crossing with 80% selectivity from |T+⟩ and |T-⟩. The fluorescence intensity in the triplet manifold is determined by the ground-state polarization, and we show successful optical read-out of the ground-state spin following microwave manipulations by fluorescence-detected magnetic resonance spectroscopy. At 85 K, the lifetime of the polarized ground state is 45 ± 3 µs, and the ground state phase memory time is Tm = 5.9 ± 0.1 µs, which increases to 26.8 ± 1.6 µs at 5 K. These results show that luminescent diradicals with triplet ground states can serve as optically addressable molecular qubits with long spin coherence times, which marks an important step toward the rational design of spin-optical interfaces in organic materials.

Linear response of molecular polaritons.

In this article, we show that the collective light-matter strong coupling regime where N molecular emitters couple to the photon mode of an optical cavity can be mapped to a quantum impurity model where the photon is the impurity that is coupled to a bath of anharmonic transitions. In the thermodynamic limit where N ≫ 1, we argue that the bath can be replaced with an effective harmonic bath, leading to a dramatic simplification of the problem into one of the coupled harmonic oscillators. We derive simple analytical expressions for linear optical spectra (transmission, reflection, and absorption) where the only molecular input required is the molecular linear susceptibility. This formalism is applied to a series of illustrative examples, showing the role of temperature, disorder, vibronic coupling, and optical saturation of the molecular ensemble, explaining that it is useful even when describing an important class of nonlinear optical experiments. For completeness, we provide Appendixes A-C that include a self-contained derivation of the relevant spectroscopic observables for arbitrary anharmonic systems (for both large and small N) within the rotating-wave approximation. While some of the presented results herein have already been reported in the literature, we provide a unified presentation of the results as well as new interpretations that connect powerful concepts in open quantum systems and linear response theory with molecular polaritonics.

Control of Photoswitching Kinetics with Strong Light-Matter Coupling in a Cavity.

Most photochemistry occurs in the regime of weak light-matter coupling, in which a molecule absorbs a photon and then performs photochemistry from its excited state. In the strong coupling regime, enhanced light-matter interactions between an optical field and multiple molecules lead to collective hybrid light-matter states called polaritons. This strong coupling leads to fundamental changes in the nature of the excited states including multi-molecule delocalized excitations, modified potential energy surfaces, and dramatically altered energy levels relative to non-coupled molecules. The effect of strong light-matter coupling on covalent photochemistry has not been well explored. Photoswitches undergo reversible intramolecular photoreactions that can be readily monitored spectroscopically. In this work, we study the effect of strong light-matter coupling on the kinetics of photoswitching within optical cavities. Reproducing prior experiments, photoswitching of spiropyran/merocyanine photoswitches is decelerated in a cavity. Fulgide photoswitches, however, show the opposite effect, with strong coupling accelerating photoswitching. While modified merocyanine switching can be explained by changes in radiative decay rates or the amount of light in the cavity, modified fulgide switching kinetics suggest direct changes to excited-state reaction kinetics.

Catalysis by Dark States in Vibropolaritonic Chemistry.

Collective strong coupling between a disordered ensemble of N localized molecular vibrations and a resonant optical cavity mode gives rise to two polariton and N-1≫2 dark modes. Thus, experimental changes in thermally activated reaction kinetics due to polariton formation appear entropically unlikely and remain a puzzle. Here we show that the overlooked dark modes, while parked at the same energy as bare molecular vibrations, are robustly delocalized across ∼2-3 molecules, yielding enhanced channels of vibrational cooling, concomitantly catalyzing a chemical reaction. As an illustration, we theoretically show an ≈50% increase in an electron transfer rate due to enhanced product stabilization. The reported effects can arise when the homogeneous linewidths of the dark modes are smaller than their energy spacings.

Generalization of the Tavis-Cummings model for multi-level anharmonic systems: Insights on the second excitation manifold.

Confined electromagnetic modes strongly couple to collective excitations in ensembles of quantum emitters, producing light-matter hybrid states known as polaritons. Under such conditions, the discrete multilevel spectrum of molecular systems offers an appealing playground for exploring multiphoton processes. This work contrasts predictions from the Tavis-Cummings model in which the material is a collection of two-level systems, with the implications of considering additional energy levels with harmonic and anharmonic structures. We discuss the exact eigenspectrum, up to the second excitation manifold, of an arbitrary number N of oscillators collectively coupled to a single cavity mode in the rotating-wave approximation. Elaborating on our group-theoretic approach [New J. Phys. 23, 063081 (2021)], we simplify the brute-force diagonalization of N2 × N2 Hamiltonians to the eigendecomposition of, at most, 4 × 4 matrices for arbitrary N. We thoroughly discuss the eigenstates and the consequences of weak and strong anharmonicities. Furthermore, we find resonant conditions between bipolaritons and anharmonic transitions where two-photon absorption can be enhanced. Finally, we conclude that energy shifts in the polaritonic states induced by anharmonicities become negligible for large N. Thus, calculations with a single or few quantum emitters qualitatively fail to represent the nonlinear optical response of the collective strong coupling regime. Our work highlights the rich physics of multilevel anharmonic systems coupled to cavities absent in standard models of quantum optics. We also provide concise tabulated expressions for eigenfrequencies and transition amplitudes, which should serve as a reference for future spectroscopic studies of molecular polaritons.

Nonequilibrium effects of cavity leakage and vibrational dissipation in thermally activated polariton chemistry.

In vibrational strong coupling (VSC), molecular vibrations strongly interact with the modes of an optical cavity to form hybrid light-matter states known as vibrational polaritons. Experiments show that the kinetics of thermally activated chemical reactions can be modified by VSC. Transition-state theory, which assumes that internal thermalization is fast compared to reactive transitions, has been unable to explain the observed findings. Here, we carry out kinetic simulations to understand how dissipative processes, namely, those introduced by VSC to the chemical system, affect reactions where internal thermalization and reactive transitions occur on similar timescales. Using the Marcus-Levich-Jortner type of electron transfer as a model reaction, we show that such dissipation can change reactivity by accelerating internal thermalization, thereby suppressing nonequilibrium effects that occur in the reaction outside the cavity. This phenomenon is attributed mainly to cavity decay (i.e., photon leakage), but a supporting role is played by the relaxation between polaritons and dark states. When nonequilibrium effects are already suppressed in the bare reaction (the reactive species are essentially at internal thermal equilibrium throughout the reaction), we find that reactivity does not change significantly under VSC. Connections are made between our results and experimental observations.

Two-dimensional infrared spectroscopy of vibrational polaritons.

We report experimental 2D infrared (2D IR) spectra of coherent light-matter excitations--molecular vibrational polaritons. The application of advanced 2D IR spectroscopy to vibrational polaritons challenges and advances our understanding in both fields. First, the 2D IR spectra of polaritons differ drastically from free uncoupled excitations and a new interpretation is needed. Second, 2D IR uniquely resolves excitation of hybrid light-matter polaritons and unexpected dark states in a state-selective manner, revealing otherwise hidden interactions between them. Moreover, 2D IR signals highlight the impact of molecular anharmonicities which are applicable to virtually all molecular systems. A quantum-mechanical model is developed which incorporates both nuclear and electrical anharmonicities and provides the basis for interpreting this class of 2D IR spectra. This work lays the foundation for investigating phenomena of nonlinear photonics and chemistry of molecular vibrational polaritons which cannot be probed with traditional linear spectroscopy.

Polaritonic normal modes in transition state theory.

A series of experiments demonstrates that strong light-matter coupling between vibrational excitations in isotropic solutions of molecules and resonant infrared optical microcavity modes leads to modified thermally activated kinetics. However, Galego et al. [Phys. Rev. X 9, 021057 (2019)] recently demonstrated that, within transition state theory, effects of strong light-matter coupling with reactive modes are mostly electrostatic and essentially independent of light-matter resonance or even of the formation of vibrational polaritons. To analyze this puzzling theoretical result in further detail, we revisit it under a new light, invoking a normal mode analysis of the transition state and reactant configurations for an ensemble of an arbitrary number of molecules in a cavity, obtaining simple analytical expressions that produce similar conclusions as Feist. While these effects become relevant in optical microcavities if the molecular dipoles are anisotropically aligned, or in cavities with extreme confinement of the photon modes, they become negligible for isotropic solutions in microcavities. It is concluded that further studies are necessary to track the origin of the experimentally observed kinetics.

Computational method for highly constrained molecular dynamics of rigid bodies: Coarse-grained simulation of auxetic two-dimensional protein crystals.

The increasing number of protein-based metamaterials demands reliable and efficient theoretical and computational methods to study the physicochemical properties they may display. In this regard, we develop a simulation strategy based on Molecular Dynamics (MD) that addresses the geometric degrees of freedom of an auxetic two-dimensional protein crystal. This model consists of a network of impenetrable rigid squares linked through massless rigid rods. Our MD methodology extends the well-known protocols SHAKE and RATTLE to include highly non-linear holonomic and non-holonomic constraints, with an emphasis on collision detection and response between anisotropic rigid bodies. The presented method enables the simulation of long-time dynamics with reasonably large time steps. The data extracted from the simulations allow the characterization of the dynamical correlations featured by the protein subunits, which show a persistent motional interdependence across the array. On the other hand, non-holonomic constraints (collisions between subunits) increase the number of inhomogeneous deformations of the network, thus driving it away from an isotropic response. Our work provides the first long-timescale simulation of the dynamics of protein crystals and offers insights into promising mechanical properties afforded by these materials.

State-Selective Polariton to Dark State Relaxation Dynamics.

The modification of vibrational dynamics is essential for controlling chemical reactions and IR photonic applications. The hybridization between cavity modes and molecular vibrational modes provides a new way to control molecular dynamics. In this work, we study the dynamics of molecular vibrational polaritons in various solvent environments. We find the dynamics of the polariton system is strongly influenced by the nature of the solvents. While the relaxation from upper polariton (UP) to dark modes is always fast (<5 ps) regardless of the medium, lower polariton (LP) in low polarity solvents shows much slower transfer (10-30 ps) into dark modes, despite the fact that the LP lifetime remains within 5 ps. This result suggests that in the latter media, the energy pumped into the LP is first transferred into intermediate states which only subsequently decay into dark modes. In contrast, in solvent environments that strongly interact with the solute, the LP population relaxes into the dense dark state manifold within a much faster time scale. We propose the intermediate state to be the high-lying excited states of dark modes, which are effectively populated by LP via, e.g., ladder-climbing. Such population in the high-lying states can be retained for tens of picoseconds, which could be pertinent to recently observed cavity-modified chemistry.

Topologically protected excitons in porphyrin thin films.

The control of exciton transport in organic materials is of fundamental importance for the development of efficient light-harvesting systems. This transport is easily deteriorated by traps in the disordered energy landscape. Here, we propose and analyse a system that supports topological Frenkel exciton edge states. Backscattering of these chiral Frenkel excitons is prohibited by symmetry, ensuring that the transport properties of such a system are robust against disorder. To implement our idea, we propose a two-dimensional periodic array of tilted porphyrins interacting with a homogeneous magnetic field. This field serves to break time-reversal symmetry and results in lattice fluxes that mimic the Aharonov-Bohm phase acquired by electrons. Our proposal is the first blueprint for realizing topological phases of matter in molecular aggregates and suggests a paradigm for engineering novel excitonic materials.

Practical witness for electronic coherences.

The origin of the coherences in two-dimensional spectroscopy of photosynthetic complexes remains disputed. Recently, it has been shown that in the ultrashort-pulse limit, oscillations in a frequency-integrated pump-probe signal correspond exclusively to electronic coherences, and thus such experiments can be used to form a test for electronic vs. vibrational oscillations in such systems. Here, we demonstrate a method for practically implementing such a test, whereby pump-probe signals are taken at several different pulse durations and used to extrapolate to the ultrashort-pulse limit. We present analytic and numerical results determining requirements for pulse durations and the optimal choice of pulse central frequency, which can be determined from an absorption spectrum. Our results suggest that for numerous systems, the required experiment could be implemented by many ultrafast spectroscopy laboratories using pulses of tens of femtoseconds in duration. Such experiments could resolve the standing debate over the nature of coherences in photosynthetic complexes.

Quantum state and process tomography of energy transfer systems via ultrafast spectroscopy.

The description of excited state dynamics in energy transfer systems constitutes a theoretical and experimental challenge in modern chemical physics. A spectroscopic protocol that systematically characterizes both coherent and dissipative processes of the probed chromophores is desired. Here, we show that a set of two-color photon-echo experiments performs quantum state tomography (QST) of the one-exciton manifold of a dimer by reconstructing its density matrix in real time. This possibility in turn allows for a complete description of excited state dynamics via quantum process tomography (QPT). Simulations of a noisy QPT experiment for an inhomogeneously broadened ensemble of model excitonic dimers show that the protocol distills rich information about dissipative excitonic dynamics, which appears nontrivially hidden in the signal monitored in single realizations of four-wave mixing experiments.

Chiral edge waves in a dance-based human topological insulator.

Topological insulators are insulators in the bulk but feature chiral energy propagation along the boundary. This property is topological in nature and therefore robust to disorder. Originally discovered in electronic materials, topologically protected boundary transport has since been observed in many other physical systems. Thus, it is natural to ask whether this phenomenon finds relevance in a broader context. We choreograph a dance in which a group of humans, arranged on a square grid, behave as a topological insulator. The dance features unidirectional flow of movement through dancers on the lattice edge. This effect persists when people are removed from the dance floor. Our work extends the applicability of wave physics to dance.

Remarks on time-dependent [current]-density functional theory for open quantum systems.

Time-dependent [current]-density functional theory for open quantum systems (OQS) has emerged as a formalism that can incorporate dissipative effects in the dynamics of many-body quantum systems. Here, we review and clarify some formal aspects of these theories that have been recently questioned in the literature. In particular, we provide theoretical support for the following conclusions: (1) contrary to what we and others had stated before, within the master equation framework, there is in fact a one-to-one mapping between vector potentials and current densities for fixed initial state, particle-particle interaction, and memory kernel; (2) regardless of the first conclusion, all of our recently suggested Kohn-Sham (KS) schemes to reproduce the current and particle densities of the original OQS, and in particular, the use of a KS closed driven system, remains formally valid; (3) the Lindblad master equation maintains the positivity of the density matrix regardless of the time-dependence of the Hamiltonian or the dissipation operators; (4) within the stochastic Schrödinger equation picture, a one-to-one mapping from stochastic vector potential to stochastic current density for individual trajectories has not been proven so far, except in the case where the vector potential is the same for every member of the ensemble, in which case, it reduces to the Lindblad master equation picture; (5) master equations may violate certain desired properties of the density matrix, such as positivity, but they remain as one of the most useful constructs to study OQS when the environment is not easily incorporated explicitly in the calculation. The conclusions support our previous work as formally rigorous, offer new insights into it, and provide a common ground to discuss related theories.