An official website of the United States government
Ion receptors are molecular hosts that bind ionic guests, often with great selectivity. The interplay of solvation and ion binding in anion host-guest complexes in solution governs the binding efficiency and selectivity of such ion receptors. To gain molecular-level insight into the intrinsic binding properties of octamethyl calix[4]pyrrole (omC4P) host molecules with halide guest ions, we performed cryogenic ion vibrational spectroscopy (CIVS) of omC4P in complexes with fluoride, chloride, and bromide ions. We interpret the spectra using density functional theory, describing the infrared spectra of these complexes with both harmonic and anharmonic second-order vibrational perturbation theory (VPT2) calculations. The NH stretching modes of the pyrrole moieties serve as sensitive probes of the ion binding properties, as their frequencies encode the ion-receptor interactions. While scaled harmonic spectra reproduce the experimental NH stretching modes of the chloride and bromide complexes in broad strokes, the high proton affinity of fluoride introduces strong anharmonic effects. As a result, the spectrum of F−·omC4P is not even qualitatively captured by harmonic calculations, but it is recovered very well by VPT2 calculations. In addition, the VPT2 calculations recover the intricate coupling of the NH stretching modes with overtones and combination bands of CH stretching and NH bending modes and with low-frequency vibrations of the omC4P macrocycle, which are apparent for all halide ion complexes investigated here. A comparison of the CIVS spectra with infrared spectra of solutions of the same ion-receptor complexes in d3-acetonitrile and d6-acetone shows how ion solvation changes the ion-receptor interactions for the different halide ions.
more »
« less
This content will become publicly available on November 1, 2025
Simulating anharmonic vibrational polaritons beyond the long wavelength approximation
In this work, we investigate anharmonic vibrational polaritons formed due to strong light–matter interactions in an optical cavity between radiation modes and anharmonic vibrations beyond the long-wavelength limit. We introduce a conceptually simple description of light–matter interactions, where spatially localized cavity radiation modes couple to localized vibrations. Within this theoretical framework, we employ self-consistent phonon theory and vibrational dynamical mean-field theory to efficiently simulate momentum-resolved vibrational-polariton spectra, including effects of anharmonicity. Numerical simulations in model systems demonstrate the accuracy and applicability of our approach.
more »
« less
- Award ID(s):
- 2245592
- PAR ID:
- 10581698
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 162
- Issue:
- 1
- ISSN:
- 0021-9606
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Shea, Joan-Emma (Ed.)We aim to automatize the identification of collective variables to simplify and speed up enhanced sampling simulations of conformational dynamics in biomolecules. We focus on anharmonic low-frequency vibrations that exhibit fluctuations on timescales faster than conformational transitions but describe a path of least resistance towards structural change. A key challenge is that harmonic approximations are ill-suited to characterize these vibrations, which are observed at far-infrared frequencies and are easily excited by thermal collisions at room temperature. Here, we approached this problem with a frequency-selective anharmonic (FRESEAN) mode analysis that does not rely on harmonic approximations and successfully isolates anharmonic low-frequency vibrations from short molecular dynamics simulation trajectories. We applied FRESEAN mode analysis to simulations of alanine dipeptide, a common test system for enhanced sampling simulation protocols, and compare the performance of isolated low-frequency vibrations to conventional user-defined collective variables (here backbone dihedral angles) in enhanced sampling simulations. The comparison shows that enhanced sampling along anharmonic low-frequency vibrations not only reproduces known conformational dynamics but can even further improve sampling of slow transitions compared to user-defined collective variables. Notably, free energy surfaces spanned by low-frequency anharmonic vibrational modes exhibit lower barriers associated with conformational transitions relative to representations in backbone dihedral space. We thus conclude that anharmonic low-frequency vibrations provide a promising path for highly effective and fully automated enhanced sampling simulations of conformational dynamics in biomolecules.more » « less
-
The vibrational dynamics of diborane have been extensively studied both theoretically and experimentally ever since the bridge structure of diborane was established in the 1950s. Numerous infrared and several Raman spectroscopic studies have followed in the ensuing years at ever increasing levels of spectral resolution. In parallel, ab initio computations of the underlying potential energy surface have progressed as well as the methods to calculate the anharmonic vibration dynamics beyond the double harmonic approximation. Nevertheless, even 70 years after the bridge structure of diborane was established, there are still significant discrepancies between experiment and theory for the fundamental vibrational frequencies of diborane. In this work we use para-hydrogen (pH2) matrix isolation infrared spectroscopy to characterize six fundamental vibrations of B2H6 and B2D6 and compare them with results from configuration-selective vibrational configuration interaction theory. The calculated frequencies and intensities are in very good agreement with the pH2 matrix isolation spectra, even several combination bands are well reproduced. We believe that the reason discrepancies have existed for so long is related to the large amount of anharmonicity that is associated with the bridge BH stretching modes. However, the calculated frequencies and intensities reported here for the vibrational modes of all three boron isotopologues of B2H6 and B2D6 are within ± 2.00 cm− 1 and ± 1.44 cm− 1, respectively, of the experimental frequencies and therefore a refined vibrational assignment of diborane has been achieved.more » « less
-
Plasmon-phonon coupling between metamaterials and molecular vibrations provides a new path for studying mid-infrared light-matter interactions and molecular detection. So far, the coupling between the plasmonic resonances of metamaterials and the phonon vibrational modes of molecules has been realized under linearly polarized light. Here, mid-infrared chiral plasmonic metasurfaces with high circular dichroism (CD) in absorption over 0.65 in the frequency range of 50 to 60 THz are demonstrated to strongly interact with the phonon vibrational resonance of polymethyl methacrylate (PMMA) molecules at 52 THz, under both left-handed and right-handed circularly polarized (LCP and RCP) light. The mode splitting features in the absorption spectra of the coupled metasurface-PMMA systems under both circular polarizations are studied in PMMA layers with different thicknesses. The relation between the mode splitting gap and the PMMA thickness is also revealed. The demonstrated results can be applied in areas of chiral molecular sensing, thermal emission, and thermal energy harvesting.more » « less
-
Simulations of anharmonic vibrational motion rely on computationally expedient representations of the governing potential energy surface. The n-mode representation (n-MR)—effectively a many-body expansion in the space of molecular vibrations—is a general and efficient approach that is often used for this purpose in vibrational self-consistent field (VSCF) calculations and correlated analogues thereof. In the present analysis, a lack of convergence in many VSCF calculations is shown to originate from negative and unbound potentials at truncated orders of the n-MR expansion. For cases of strong anharmonic coupling between modes, the n-MR can both dip below the true global minimum of the potential surface and lead to effective single-mode potentials in VSCF that do not correspond to bound vibrational problems, even for bound total potentials. The present analysis serves mainly as a pathology report of this issue. Furthermore, this insight into the origin of VSCF non-convergence provides a simple, albeit ad hoc, route to correct the problem by “painting in” the full representation of groups of modes that exhibit these negative potentials at little additional computational cost. Somewhat surprisingly, this approach also reasonably approximates the results of the next-higher n-MR order and identifies groups of modes with particularly strong coupling. The method is shown to identify and correct problematic triples of modes—and restore SCF convergence—in two-mode representations of challenging test systems, including the water dimer and trimer, as well as protonated tropine.more » « less
