Measuring electric fields and noncovalent interactions using the vibrational stark effect

Forced Interactions: Ionic Polymers at Charged Surfactant Interfaces

Characterizing electrostatic interactions at heterogeneous interfaces is critical for developing a fundamental description of the dynamic processes at charged interfaces. Water-in-oil reverse micelles (RMs) offer a high degree of tunability across composition, polarity, and temperature, making them ideal systems for studying interactions at heterogeneous liquid-liquid interfaces. In the present study, we use a combination of ultrafast two-dimensional infrared spectroscopy and molecular dynamics (MD) simulations to determine the picosecond interfacial dynamics in RMs containing binary compositions of sorbitan monostearate and anionic or cationic cosurfactants, which are used to tune the ratio of charged to nonionic surfactants at the interface. The positively charged polyethylenimine (PEI) polymer is encapsulated within the RMs, and the carbonyl stretching mode of sorbitan monostearate reports on the interfacial hydrogen-bond populations and dynamics. The results show that hydrogen-bond populations are altered through the inclusion of both negatively and positively charged cosurfactants. Charged surfactants increase interfacial water penetration into the surfactant layer, and the surface localization of polymers decreases water penetration. Local hydrogen-bond dynamics undergo a slowdown with the inclusion of charged surfactants, and the encapsulation of polymers results in similar effects, irrespective of the charge.

Theoretical Investigation of the X-ray Stark Effect in Small Molecules

We have studied the Stark effect in the soft x-ray region for various small molecules by calculating the field-dependent x-ray absorption spectra. This effect is explained in terms of the response of molecular orbitals (core and valence), the molecular dipole moment, and the molecular geometry to the applied electric field. A number of consistent trends are observed linking the computed shifts in absorption energies and intensities with specific features of the molecular electronic structure. We find that both the virtual molecular orbitals (valence and/or Rydberg) as well as the core orbitals contribute to observed trends in a complementary fashion. This initial study highlights the potential impact of x-ray Stark spectroscopy as a tool to study electronic structure and environmental perturbations at a submolecular scale.

Quantification of the Real Plasmonic Field Transverse Distribution in a Nanocavity Using the Vibrational Stark Effect

Quantifying the real plasmonic field strength experimentally has been long pursued in expanding the applications related to plasmonic enhancement. However, it is still an enormous challenge to determine the inhomogeneous plasmonic field distribution. Here, self-assembled monolayers (SAMs) of 4-mercaptobenzonitrile (MBN) are sandwiched as a gap spacer in a nanoparticle-on-mirror (NPoM) structure, effectively forming ultrahigh field enhancement to observe Stark shifts of the chemical bond. Transverse position-dependent Stark shifts of ν(C═C) and ν(C≡N) in the individual nanocavity measured by surface-enhanced Raman scattering (SERS) experiment combined with the Stark tuning rate by density functional theory (DFT) simulation accurately revealed the inhomogeneous plasmonic field transverse distribution and quantified the transverse plasmonic field strength up to ∼1.9 × 10⁹ V/m, which matches the value predicted by finite element method (FEM) simulation. This work deepens the insight into plasmon-based technologies and will coordinate high-resolution techniques such as tip-enhanced Raman spectroscopy (TESR) to reveal the real plasmonic field distribution.

Force Decomposition Analysis: A Method to Decompose Intermolecular Forces into Physically Relevant Component Contributions

Computational quantum chemistry can be more than just numerical experiments when methods are specifically adapted to investigate chemical concepts. One important example is the development of energy decomposition analysis (EDA) to reveal the physical driving forces behind intermolecular interactions. In EDA, typically the interaction energy from a good-quality density functional theory (DFT) calculation is decomposed into multiple additive components that unveil permanent and induced electrostatics, Pauli repulsion, dispersion, and charge-transfer contributions to noncovalent interactions. Herein, we formulate, implement, and investigate decomposing the forces associated with intermolecular interactions into the same components. The resulting force decomposition analysis (FDA) is potentially useful as a complement to the EDA to understand chemistry, while also providing far more information than an EDA for data analysis purposes such as training physics-based force fields. We apply the FDA based on absolutely localized molecular orbitals (ALMOs) to analyze interactions of water with sodium and chloride ions as well as in the water dimer. We also analyze the forces responsible for geometric changes in carbon dioxide upon adsorption onto (and activation by) gold and silver anions. We also investigate how the force components of an EDA-based force field for water clusters, namely MB-UCB, compare to those from force decomposition analysis.

Electroabsorption and Stark Fluorescence Spectroscopies of Thioflavin T

Thioflavin T (ThT) is a typical fluorescent marker for detecting the formation of amyloid fibrils, because its fluorescence intensity increases by more than 2 orders of magnitude upon complexation with the fibrils. Strong electrostatic fields on protein surfaces are known to be a significant factor in chemical reactions and biological functions. Therefore, ThT bound to amyloid fibrils must experience strong electric fields. This study employed electroabsorption and Stark fluorescence spectroscopies to clarify the effects of external electric fields on the photophysics of ThT. The absorption spectrum shows two bands ascribed to locally excited (LE) and charge transfer (CT) states. Coupling between the LE and CT states is enhanced in the presence of an external electric field, resulting in fluorescence quenching. The electric field strength of the amyloid fibril surface was inferred from the fluorescence quenching efficiency of ThT.

TFRegNCI: Interpretable Noncovalent Interaction Correction Multimodal Based on Transformer Encoder Fusion

The interpretability is an important issue for end-to-end learning models. Motivated by computer vision algorithms, an interpretable noncovalent interaction (NCI) correction multimodal (TFRegNCI) is proposed for NCI prediction. TFRegNCI is based on RegNet feature extraction and a transformer encoder fusion strategy. RegNet is a network design paradigm that mainly focuses on local features. Meanwhile, the Vision Transformer is also leveraged for feature extraction, because it can capture global features better than RegNet while lowering the computational cost. Using a transformer encoder as the fusion strategy rather than multilayer perceptron can enhance model performance, due to its emphasis on important features with less parameters. Therefore, the proposed TFRegNCI achieved high accurate prediction (mean absolute error of ∼0.1 kcal/mol) comparing with the coupled cluster single double (triple) (CCSD(T)) benchmark. To further improve the model efficiency, TFRegNCI applies two-dimensional (2D) inputs transformed from three-dimensional (3D) electron density cubes, which saves time (30%), while the model accuracy remains. To improve model interpretability, a visualization module, Gradient-weighted Regression Activation Mapping (Grad-RAM) has been embedded. Grad-RAM is promoted from the classification algorithm, Gradient-weighted Class Activation Mapping, to perform feature visualization for the regression task. With Grad-RAM, the visual location map for features in deep learning models can be displayed. The feature map visualizations suggest that the 2D model has the similar performance as the 3D model, because of equally effective feature extractions from electron density. Moreover, the valid feature region on the location map by the 3D model is consistent with the NCIPLOT NCI isosurface. It is confirmed that the model does extract significant features related to the NCI interaction. The interpretable analyses are carried out through molecular orbital contribution on effective features. Thereby, the proposed model is likely to be a promising tool to reveal some essential information on NCIs, with regard to the level of electronic theory.

Rethinking Vibrational Stark Spectroscopy: Peak Shifts, Line Widths, and the Role of Non-Stark Solvent Coupling

A vibration's transition frequency is partly determined by the first-order Stark effect, which accounts for the electric field experienced by the mode. Using ultrafast infrared pump-probe and FT-IR spectroscopies, we characterized both the 0 → 1 and 1 → 2 vibrational transitions' field-dependent peak positions and line widths of the CN stretching mode of benzonitrile (BZN) and phenyl selenocyanate (PhSeCN) in ten solvents. We present a theoretical model that decomposes the observed line width into a field-dependent Stark contribution and a field-independent non-Stark solvent coupling contribution (NSC). The model demonstrates that the field-dependent peak position is independent of the line width, even when the NSC dominates the latter. Experiments show that when the Stark tuning rate is large compared to the NSC (PhSeCN), the line width has a field dependence, albeit with major NSC-induced excursions from linearity. When the Stark tuning rate is small relative to the NSC (BZN), the line width is field-independent. BZN's line widths are substantially larger for the 1 → 2 transition, indicating a 1 → 2 transition enhancement of the NSC. Additionally, we examine, theoretically and experimentally, the difference in the 0 → 1 and 1 → 2 transitions' Stark tuning rates. Second-order perturbation theory combined with density functional theory explain the difference and show that the 1 → 2 transition's Stark tuning rate is ∼10% larger. The Stark tuning rate of PhSeCN is larger than BZN's for both transitions, consistent with the theoretical calculations. This study provides new insights into vibrational line shape components and a more general understanding of the vibrational response to external electric fields.

4-Cyanotryptophan as a Sensitive Fluorescence Probe of Local Electric Field of Proteins

Electrostatic interactions are key determinants of protein structure, dynamics, and function. Since protein electrostatics are nonuniform, assessment of the internal electric fields (EFs) of proteins requires spatial resolution at the amino acid residue level. In this regard, vibrational Stark spectroscopy, in conjunction with various unnatural amino acid-based vibrational probes, has become a common method for site-specific interrogation of protein EFs. However, application of this method is often limited to proteins with relatively high solubility, due to the intrinsically low oscillator strength of vibrational transitions. Therefore, it would be useful to develop an alternative method that can overcome this limitation. To this end, we show that, using solvatochromic study and molecular dynamics simulations, the frequency of maximum emission intensity of the fluorophore of 4-cyanotryptophan (4CN-Trp), 3-methyl-1H-indole-4-carbonitrile, exhibits a linear dependence on the local EF. Since the absorption and emission spectra of 4CN-Trp are easily distinguishable from those of naturally occurring aromatic amino acids, we believe that this linear relationship provides an easier and more sensitive means to determine the local EF of proteins.

Electrostatics and Chemical Reactivity at the Air-Water Interface

It has been recently discovered that chemical reactions at aqueous interfaces can be orders of magnitude faster compared to conventional bulk phase reactions, but despite its wide-ranging implications, which extend from atmospheric to synthetic chemistry or technological applications, the phenomenon is still incompletely understood. The role of strong electric fields due to space asymmetry and the accumulation of ions at the interface has been claimed as a possible cause from some experiments, but the reorganization of the solvent around the reactive system should provide even greater additional electrostatic contributions that have not yet been analyzed. In this study, with the help of first-principles molecular dynamics simulations, we go deeper into this issue by a careful assessment of solvation electrostatics at the air-water interface. Our simulations confirm that electrostatic forces can indeed be a key factor in rate acceleration compared to bulk solution. Remarkably, the study reveals that the effect cannot simply be attributed to the magnitude of the local electric field and that the fluctuations of the full electrostatic potential resulting from unique dynamical behavior of the solvation shells at the interface must be accounted for. This finding paves the way for future applications of the phenomenon in organic synthesis, especially for charge transfer or redox reactions in thin films and microdroplets.

Frequency Changes in Terminal Alkynes Provide Strong, Sensitive, and Solvatochromic Raman Probes of Biochemical Environments

The C≡C stretching frequencies of terminal alkynes appear in the "clear" window of vibrational spectra, so they are attractive and increasingly popular as site-specific probes in complicated biological systems like proteins, cells, and tissues. In this work, we collected infrared (IR) absorption and Raman scattering spectra of model compounds, artificial amino acids, and model proteins that contain terminal alkyne groups, and we used our results to draw conclusions about the signal strength and sensitivity to the local environment of both aliphatic and aromatic terminal alkyne C≡C stretching bands. While the IR bands of alkynyl model compounds displayed surprisingly broad solvatochromism, their absorptions were weak enough that alkynes can be ruled out as effective IR probes. The same solvatochromism was observed in model compounds' Raman spectra, and comparisons to published empirical solvent scales (including a linear regression against four meta-aggregated solvent parameters) suggested that the alkyne C≡C stretching frequency mainly reports on local electronic interactions (i.e., short-range electron donor-acceptor interactions) with solvent molecules and neighboring functional groups. The strong solvatochromism observed here for alkyne stretching bands introduces an important consideration for Raman imaging studies based on these signals. Raman signals for alkynes (especially those that are π-conjugated) can be exceptionally strong and should permit alkynyl Raman signals to function as probes at very low concentrations, as compared to other widely used vibrational probe groups like azides and nitriles. We incorporated homopropargyl glycine into a transmembrane helical peptide via peptide synthesis, and we installed p-ethynylphenylalanine into the interior of the Escherichia coli fatty acid acyl carrier protein using a genetic code expansion technique. The Raman spectra from each of these test systems indicate that alkynyl C≡C bands can act as effective and unique probes of their local biomolecular environments. We provide guidance for the best possible future uses of alkynes as solvatochromic Raman probes, and while empirical explanations of the alkyne solvatochromism are offered, open questions about its physical basis are enunciated.

Inconsistent hydrogen bond-mediated vibrational coupling of amide I

Using infrared spectroscopy and density functional theory (DFT) calculations, we scrutinized an amide (dimethylformamide) as a "model" compound to interpret the interactions of amide 1 with different phenol derivatives (para-chlorophenol (PCP) and para-cresol (CP)) as "model guest molecules". We established the involvement of amide I in vibrational coupling with symmetric and asymmetric C[double bond, length as m-dash]C modes of different phenolic derivatives and how their coupling was dependent upon different guest aromatic phenolic compounds. Interestingly, substitution of phenol perturbed the pattern of vibrational coupling with amide I. The symmetric and asymmetric C[double bond, length as m-dash]C modes of PC were coupled significantly with amide 1. For PCP, the symmetric C[double bond, length as m-dash]C mode coupled significantly, but the asymmetric mode coupled negligibly, with amide I. Here, we reveal the nature of vibrational coupling based on the structure of a guest molecule hydrogen-bonded with amide I. Our conclusions could be valuable for depiction of the unusual dynamics of coupled amide-I modes as well as the dependency of vibrational coupling on altered factors.

Blue Fluorescence of Cyano-tryptophan Predicts Local Electrostatics and Hydrogen Bonding in Biomolecules

Cyano-tryptophan is an unnatural fluorescent amino acid that emits in the visible region. Along with the structural similarity with tryptophan, the unique photophysical properties of this fluorophore make it an ideal probe for biophysical research. Herein, combining fluorescence spectroscopy, infrared spectroscopy, and molecular dynamics simulations, we show that the cyano-tryptophan's emission energy quantifies the underlying bond-specific noncovalent interactions in terms of the electric field. We further report the use of fluorophore's emission energy to predict its hydrogen bond characteristics. We demonstrate that combining experiments with molecular dynamics simulations can provide the hydrogen bonding status of the nitrile moiety. In addition, we report a method to differentiate between aqueous and nonaqueous hydrogen-bonding partners. Using a phenomenological approach, we demonstrate that the presence of the cyano-indole moiety is responsible for the distinct correlations between the fluorophore's emission and the electrostatic forces on the nitrile bond. As indole is a privileged scaffold for both native amino acids and nucleobases, cyano-indoles will have many multifaceted applications.

Brunner syndrome caused by point mutation explained by multiscale simulation of enzyme reaction

Brunner syndrome is a disorder characterized by intellectual disability and impulsive, aggressive behavior associated with deficient function of the monoamine oxidase A (MAO-A) enzyme. These symptoms (along with particularly high serotonin levels) have been reported in patients with two missense variants in MAO-A (p.R45W and p.E446K). Herein, we report molecular simulations of the rate-limiting step of MAO-A-catalyzed serotonin degradation for these variants. We found that the R45W mutation causes a 6000-fold slowdown of enzymatic function, whereas the E446K mutation causes a 450-fold reduction of serotonin degradation rate, both of which are practically equivalent to a gene knockout. In addition, we thoroughly compared the influence of enzyme electrostatics on the catalytic function of both the wild type MAO-A and the p.R45W variant relative to the wild type enzyme, revealing that the mutation represents a significant electrostatic perturbation that contributes to the barrier increase. Understanding genetic disorders is closely linked to understanding the associated chemical mechanisms, and our research represents a novel attempt to bridge the gap between clinical genetics and the underlying chemical physics.

Local Electric Fields: From Enzyme Catalysis to Synthetic Catalyst Design

This Mini-Review Article outlines recent advances in the study of local electric field (LEF) governed enzyme catalysis and the application of the LEF principle in synthetic catalyst design. We start by discussing the electrostatics principles that drive enzyme catalysis, and its experimental verifications through vibrational Stark spectroscopy. Subsequently, we describe aspects of LEFs other than catalysis, i.e., induction of mechanistic crossovers, among others. Here, we focus on the early work done using computational tools, along with some recent contributions. Following an in-depth discussion of the role of LEFs in enzyme catalysis, we then highlight some recent works on designed local electric fields (D-LEF) and their applications in organic synthesis. Subsequently, we turn to D-LEFs in synthetic enzymes and supramolecular systems (cf. the work by the Head-Gordon group). We end by discussing some of the software packages that have been developed to analyze local electric fields computationally. Overall, the present Mini-Review Article paints an insightful picture of the current state of the art using LEF in enzyme catalysis and its application for further bioengineering and synthetic organic frameworks in a broad perspective.

Predicting Frequency from the External Chemical Environment: OH Vibrations on Hydrated and Hydroxylated Surfaces

Robust correlation curves are essential to decipher structural information from IR-vibrational spectra. However, for surface-adsorbed water and hydroxides, few such correlations have been presented in the literature. In this paper, OH vibrational frequencies are correlated against 12 structural descriptors representing the quantum mechanical or geometrical environment, focusing on those external to the vibrating molecule. A nonbiased fitting procedure based on Gaussian process regression (GPR) was used alongside simple analytical functional forms. The training data consist of 217 structurally unique OH groups from 38 water/metal oxide interface systems for MgO, CaO and CeO₂, all optimized at the DFT level, and the fully anharmonic and uncoupled OH vibrational signatures were calculated. Among our results, we find the following: (i) The intermolecular R(H···O) hydrogen bond distance is particularly strong, indicating the primary cause of the frequency shift. (ii) Similarly, the electric field along the H-bond vector is also a good descriptor. (iii) Highly detailed machine learning descriptors (ACSF, SOAP) are less intuitive but were found to be more capable descriptors. (iv) Combinations of geometric and QM descriptors give the best predictions, supplying complementary information.

Modulation of Metal Carbonyl Stretching Frequencies in the Second Coordination Sphere through the Internal Stark Effect

Spectroscopic and computational examination of a homologous series of rhodium(I) pybox carbonyl complexes has revealed a correlation between the conformation of the flanking aryl-substituted oxazoline donors and the carbonyl stretching frequency. This relationship is also observed experimentally for octahedral rhodium(III) and ruthenium(II) variants and cannot be explained through the classical, Dewar-Chatt-Duncanson, interpretation of metal-carbonyl bonding. Instead, these findings are reconciled by local changes in the magnitude of the electric field that is projected along the metal-carbonyl vector: the internal Stark effect.

Resolving the Impact of Hydrogen Bonding on the Phylloquinone Cofactor through Two-Dimensional Infrared Spectroscopy

Two-dimensional infrared spectroscopy (2DIR) was applied to phylloquinone (PhQ), an important biological cofactor, to elucidate the impact of hydrogen bonding on the ultrafast dynamics and energetics of the carbonyl stretching modes. 2DIR measurements were performed on PhQ dissolved in hexanol, which served as the hydrogen bonding solvent, and hexane, which served as a non-hydrogen bonding control. Molecular dynamics simulations and quantum chemical calculations were performed to aid in spectral assignment and interpretation. From the position of the peaks in the 2DIR spectra, we extracted the transition frequencies for the fundamental, overtone, and combination bands of hydrogen bonded and non-hydrogen bonded carbonyl groups of PhQ in the 1635-1680 cm^-1 region. We find that hydrogen bonding to a single carbonyl group acts to decouple the two carbonyl units of PhQ. Through analysis of the time-resolved 2DIR data, we find that hydrogen bonding leads to faster vibrational relaxation as well as an increase in the inhomogeneous broadening of the carbonyl groups. Overall, this work demonstrates how hydrogen bonding to the carbonyl groups of PhQ presents in the 2DIR spectra, laying the groundwork to use PhQ as a 2DIR probe to characterize the ultrafast fluctuations in the local environment of natural photosynthetic complexes.

β-Lactamases Evolve against Antibiotics by Acquiring Large Active-Site Electric Fields

A compound bound covalently to an enzyme active site can act either as a substrate if the covalent linkage is readily broken up by the enzyme or as an inhibitor if the bond dissociates slowly. We tracked the reactivity of such bonds associated with the rise of the resistance to penicillin G (PenG) in protein evolution from penicillin-binding proteins (PBPs) to TEM β-lactamases and with the development of avibactam (Avb) to overcome the resistance. We found that the ester linkage in PBP-PenG is resistant to hydrolysis mainly due to the small electric fields present in the protein active site. Conversely, the same linkage in the descendant TEM-PenG experiences large electric fields that stabilize the more charge-separated transition state and thus lower the free energy barrier to hydrolysis. Specifically, the electric fields were improved from -59 to -140 MV/cm in an ancient evolution dating back billions of years, contributing 5 orders of magnitude rate acceleration. This trend continues today in the nullification of newly developed antibiotic drugs. The fast linkage hydrolysis acquired from evolution is counteracted by the upgrade of PenG to Avb whose linkage escapes from the hydrolysis by returning to a low-field environment. Using the framework of electrostatic catalysis, the electric field, an observable from vibrational spectroscopy, provides a unifying physical metric to understand protein evolution and to guide the design of covalent drugs.

Protein Electric Fields Enable Faster and Longer-Lasting Covalent Inhibition of β-Lactamases

The widespread design of covalent drugs has focused on crafting reactive groups of proper electrophilicity and positioning toward targeted amino-acid nucleophiles. We found that environmental electric fields projected onto a reactive chemical bond, an overlooked design element, play essential roles in the covalent inhibition of TEM-1 β-lactamase by avibactam. Using the vibrational Stark effect, the magnitudes of the electric fields that are exerted by TEM active sites onto avibactam's reactive C═O were measured and demonstrate an electrostatic gating effect that promotes bond formation yet relatively suppresses the reverse dissociation. These results suggest new principles of covalent drug design and off-target site prediction. Unlike shape and electrostatic complementary which address binding constants, electrostatic catalysis drives reaction rates, essential for covalent inhibition, and deepens our understanding of chemical reactivity, selectivity, and stability in complex systems.

The impact of the electric field of metal ions on the vibrations and internal hydrogen bond strength in alkali metal ion di- and triglycine complexes

Using infrared predissociation spectroscopy of cryogenic ions, we revisit the vibrational spectra of alkali metal ion (Li+, Na+, K+) di- and triglycine complexes. We assign their most stable conformation, which involves metal ion coordination to all C=O groups and an internal NH⋯NH2 hydrogen bond in the peptide backbone. An analysis of the spectral shifts of the OH and C=O stretching vibrations across the different metal ions and peptide chain lengths shows that these are largely caused by the electric field of the metal ion, which varies in strength as a function of the square of the distance. The metal ion–peptide interaction also remotely modulates the strength of internal hydrogen bonding in the peptide backbone via the weakening of the amide C=O bond, resulting in a decrease in internal hydrogen bond strength from Li+ > Na+ > K+.

Electrofreezing of Liquid Ammonia

Here we prove that, in addition to temperature and pressure, another important thermodynamic variable permits the exploration of the phase diagram of ammonia: the electric field. By means of (path integral) ab initio molecular dynamics simulations, we predict that, upon applying intense electric fields on ammonia, the electrofreezing phenomenon occurs, leading the liquid toward a novel ferroelectric solid phase. This study proves that electric fields can generally be exploited as the access key to otherwise-unreachable regions in phase diagrams, unveiling the existence of new condensed-phase structures. Furthermore, the reported findings have manifold practical implications, from the safe storage and transportation of ammonia to the understanding of the solid structures this compound forms in planetary contexts.

Computational Model for Electrochemical Surface-Enhanced Raman Scattering: Key Role of the Surface Charges and Synergy between Electromagnetic and Charge-Transfer Enhancement Mechanisms

We present a computational model for electrochemical surface-enhanced Raman scattering (EC-SERS). The surface excess of charge induced by the electrode potential (V_el) was introduced by applying an external electric field to a set of clusters [Ag_n]^q with (n, q) of (19, ±1) or (20, 0) on which a molecule adsorbs. Using DFT/TD-DFT calculations, these metal-molecule complexes were classified by the adsorbate partial charge, and the main V_el-dependent properties were simultaneously studied with the aid of vibronic resonance Raman computations, namely, changes on the vibrational wavenumbers, relative intensities, and enhancement factors (EFs) for all SERS mechanisms: chemical or nonresonant, and resonance Raman with bright states of the adsorbate, charge-transfer (CT) states, and plasmon-like excitations on the metal cluster. We selected two molecules to test our model, pyridine, for which V_el has a remarkable effect, and 9,10-bis((E)-2-(pyridin-4-yl)vinyl)anthracene, which is almost insensitive to the applied bias. The results nicely reproduced most of the experimental observations, while the limitations of our approach were critically evaluated. We detected that accounting explicitly for the surface charges is key for EC-SERS models and that the highest calculated EFs, up to 10⁷ to 10⁸, are obtained by interstate coupling of bright local excitations of the metal cluster and CT states. These results highlight the importance of nonadiabatic effects in SERS and the capabilities of EC-SERS as a technique with potential to study excited-state coupling by tuning the CT and plasmon-like states by manipulating V_el.

Potential Distribution across Model Membranes

Membrane models assembled on electrodes are widely used tools to study potential-dependent molecular processes at or in membranes. However, the relationship between the electrode potential and the potential across the membrane is not known. Here we studied lipid bilayers immobilized on mixed self-assembled monolayers (SAM) on Au electrodes. The mixed SAM was composed of thiol derivatives of different chain lengths such that between the islands of the short one, mercaptobenzonitrile (MBN), and the tethered lipid bilayer an aqueous compartment was formed. The nitrile function of MBN, which served as a reporter group for the vibrational Stark effect (VSE), was probed by surface-enhanced infrared absorption spectroscopy to determine the local electric field as a function of the electrode potential for pure MBN, mixed SAM, and the bilayer system. In parallel, we calculated electric fields at the VSE probe by molecular dynamics (MD) simulations for different charge densities on the metal, thereby mimicking electrode potential changes. The agreement with the experiments was very good for the calculations of the pure MBN SAM and only slightly worse for the mixed SAM. The comparison with the experiments also guided the design of the bilayer system in the MD setups, which were selected to calculate the electrode potential dependence of the transmembrane potential, a quantity that is not directly accessible by the experiments. The results agree very well with estimates in previous studies and thus demonstrate that the present combined experimental-theoretical approach is a promising tool for describing potential-dependent processes at biomimetic interfaces.

Electric fields control water-gated proton transfer in cytochrome c oxidase

Aerobic life is powered by membrane-bound enzymes that catalyze the transfer of electrons to oxygen and protons across a biological membrane. Cytochrome c oxidase (CcO) functions as a terminal electron acceptor in mitochondrial and bacterial respiratory chains, driving cellular respiration and transducing the free energy from O₂ reduction into proton pumping. Here we show that CcO creates orientated electric fields around a nonpolar cavity next to the active site, establishing a molecular switch that directs the protons along distinct pathways. By combining large-scale quantum chemical density functional theory (DFT) calculations with hybrid quantum mechanics/molecular mechanics (QM/MM) simulations and atomistic molecular dynamics (MD) explorations, we find that reduction of the electron donor, heme a, leads to dissociation of an arginine (Arg438)-heme a₃ D-propionate ion-pair. This ion-pair dissociation creates a strong electric field of up to 1 V Å^-1 along a water-mediated proton array leading to a transient proton loading site (PLS) near the active site. Protonation of the PLS triggers the reduction of the active site, which in turn aligns the electric field vectors along a second, "chemical," proton pathway. We find a linear energy relationship of the proton transfer barrier with the electric field strength that explains the effectivity of the gating process. Our mechanism shows distinct similarities to principles also found in other energy-converting enzymes, suggesting that orientated electric fields generally control enzyme catalysis.

Collective Mid-Infrared Vibrations in Surface-Enhanced Raman Scattering

Surface-enhanced Raman scattering (SERS) is typically assumed to occur at individual molecules neglecting intermolecular vibrational coupling. Here, we show instead how collective vibrations from infrared (IR) coupled dipoles are seen in SERS from molecular monolayers. Mixing IR-active molecules with IR-inactive spacer molecules controls the intermolecular separation. Intermolecular coupling leads to vibrational frequency upshifts up to 8 cm^-1, tuning with the mixing fraction and IR dipole strength, in excellent agreement with microscopic models and density functional theory. These cooperative frequency shifts can be used as a ruler to measure intermolecular distance and disorder with angstrom resolution. We demonstrate this for photochemical reactions of 4-nitrothiophenol, which depletes the number of neighboring IR-active molecules and breaks the collective vibration, enabling direct tracking of the reaction. Collective molecular vibrations reshape SERS spectra and need to be considered in the analysis of vibrational spectra throughout analytical chemistry and sensing.

Unconventional interfacial water structure of highly concentrated aqueous electrolytes at negative electrode polarizations

Water-in-salt electrolytes are an appealing option for future electrochemical energy storage devices due to their safety and low toxicity. However, the physicochemical interactions occurring at the interface between the electrode and the water-in-salt electrolyte are not yet fully understood. Here, via in situ Raman spectroscopy and molecular dynamics simulations, we investigate the electrical double-layer structure occurring at the interface between a water-in-salt electrolyte and an Au(111) electrode. We demonstrate that most interfacial water molecules are bound with lithium ions and have zero, one, or two hydrogen bonds to feature three hydroxyl stretching bands. Moreover, the accumulation of lithium ions on the electrode surface at large negative polarizations reduces the interfacial field to induce an unusual “hydrogen-up” structure of interfacial water and blue shift of the hydroxyl stretching frequencies. These physicochemical behaviours are quantitatively different from aqueous electrolyte solutions with lower concentrations. This atomistic understanding of the double-layer structure provides key insights for designing future aqueous electrolytes for electrochemical energy storage devices.

Vibrational Polaritons in Disordered Molecular Ensembles

Disorder is an intrinsic attribute of any realistic molecular system. It is known to lead to localization, which hampers efficient transport. It was recently proposed that in molecular ensembles strongly coupled to photonic cavities, moderate disorder leads to delocalization and increases of the transport and chemical reaction rates. Vibrational polaritons involve molecular vibrations hybridized with an infrared cavity. When the coupling strength largely exceeds the molecular inhomogeneity, polaritons are unaffected by disorder. However, in many experiments, such a homogeneous limit does not apply. We investigated vibrational polaritons involving molecular ensembles with systematically modified disorder. Counterintuitively, moderate disorder leads to an increase in Rabi splitting and the modification of the polariton bandwidths. Experimental spectroscopic data agree with a Tavis-Cummings-like model that suggests enhanced delocalization of the reservoir states occurs via the admixture of the cavity mode. Our results provide new insights into the paradigm of disorder-induced cavity-assisted delocalization in molecular polaritons.

Solvent Organization and Electrostatics Tuned by Solute Electronic Structure: Amide versus Non-Amide Carbonyls

The ability to exploit carbonyl groups to measure electric fields in enzymes and other complex reactive environments by using the vibrational Stark effect has inspired growing interest in how these fields can be measured, tuned, and ultimately designed. Previous studies have concentrated on the role of the solvent in tuning the fields exerted on the solute. Here, we explore instead the role of the solute electronic structure in modifying the local solvent organization and electric field exerted on the solute. By measuring the infrared absorption spectra of amide-containing molecules, as prototypical peptides, and contrasting them with non-amide carbonyls in a wide range of solvents, we show that these solutes experience notable differences in their frequency shifts in polar solvents. Using vibrational Stark spectroscopy and molecular dynamics simulations, we demonstrate that while some of these differences can be rationalized by using the distinct intrinsic Stark tuning rates of the solutes, the larger frequency shifts for amides and dimethylurea primarily result from the larger solvent electric fields experienced by their carbonyl groups. These larger fields arise due to their stronger p-π conjugation, which results in larger C═O bond dipole moments that further induce substantial solvent organization. Using electronic structure calculations, we decompose the electric fields into contributions from solvent molecules that are in the first solvation shell and those from the bulk and show that both of these contributions are significant and become larger with enhanced conjugation in solutes. These results show that structural modifications of a solute can be used to tune both the solvent organization and electrostatic environment, indicating the importance of a solute-centric paradigm in modulating and designing the electrostatic environment in condensed-phase chemical processes.

Sub-Nanometer Mapping of the Interfacial Electric Field Profile Using a Vibrational Stark Shift Ruler

The characterization of electrical double layers is important since the interfacial electric field and electrolyte environment directly affect the reaction mechanisms and catalytic rates of electrochemical processes. In this work, we introduce a spectroscopic method based on a Stark shift ruler that enables mapping the electric field strength across the electric double layer of electrode/electrolyte interfaces. We use the tungsten-pentacarbonyl(1,4-phenelenediisocyanide) complex attached to the gold surface as a molecular ruler. The carbonyl (CO) and isocyanide (NC) groups of the self-assembled monolayer (SAM) provide multiple vibrational reporters situated at different distances from the electrode. Measurements of Stark shifts under operando electrochemical conditions and direct comparisons to density functional theory (DFT) simulations reveal distance-dependent electric field strength from the electrode surface. This electric field profile can be described by the Gouy-Chapman-Stern model with Stern layer thickness of ∼4.5 Å, indicating substantial solvent and electrolyte penetration within the SAM. Significant electro-induction effect is observed on the W center that is ∼1.2 nm away from the surface despite rapid decay of the electric field (∼90%) within 1 nm. The applied methodology and reported findings should be particularly valuable for the characterization of a wide range of microenvironments surrounding molecular electrocatalysts at electrode interfaces and the positioning of electrocatalysts at specific distances from the electrode surface for optimal functionality.

War and Peace between Electrostatic and van der Walls Forces Regulate Translational and Rotational Diffusion

Diffusion of a Brownian particle or a molecule in a liquid solvent is caused, in the Stokes-Einstein picture, by unbalanced fluctuations of osmotic forces on different sides of the particle. When the particle carries a charge or a higher multipolar moment, this picture is amended by fluctuations of electrostatic forces producing dielectric friction. Dielectric friction slows both the translational and rotational diffusion down. While this picture is well established and is physically sound, standard theories grossly overestimate the magnitude of dielectric friction for small dipolar solutes and larger colloidal particles, such as proteins. Motivated by recent simulation studies, this perspective discusses the interplay between osmotic (van der Waals) and electrostatic forces in promoting molecular and colloidal diffusion. Much can be learned about microscopic friction mechanisms from statistical and dynamical correlations between osmotic and electrostatic forces.

A Computational Quantum-Based Perspective on the Molecular Origins of Life’s Building Blocks

The search for the chemical origins of life represents a long-standing and continuously debated enigma. Despite its exceptional complexity, in the last decades the field has experienced a revival, also owing to the exponential growth of the computing power allowing for efficiently simulating the behavior of matter—including its quantum nature—under disparate conditions found, e.g., on the primordial Earth and on Earth-like planetary systems (i.e., exoplanets). In this minireview, we focus on some advanced computational methods capable of efficiently solving the Schro¨dinger equation at different levels of approximation (i.e., density functional theory)—such as ab initio molecular dynamics—and which are capable to realistically simulate the behavior of matter under the action of energy sources available in prebiotic contexts. In addition, recently developed metadynamics methods coupled with first-principles simulations are here reviewed and exploited to answer to old enigmas and to propose novel scenarios in the exponentially growing research field embedding the study of the chemical origins of life.

Comparing interfacial cation hydration at catalytic active sites and spectator sites on gold electrodes: understanding structure sensitive CO2 reduction kinetics

Hydrated cations present in the electrochemical double layer (EDL) are known to play a crucial role in electrocatalytic CO2 reduction (CO2R), and numerous studies have attempted to explain how the cation effect contributes to the complex CO2R mechanism. CO2R is a structure sensitive reaction, indicating that a small fraction of total surface sites may account for the majority of catalytic turnover. Despite intense interest in specific cation effects, probing site-specific, cation-dependent solvation structures remains a significant challenge. In this work, CO adsorbed on Au is used as a vibrational Stark reporter to indirectly probe solvation structure using vibrational sum frequency generation (VSFG) spectroscopy. Two modes corresponding to atop adsorption of CO are observed with unique frequency shifts and potential-dependent intensity profiles, corresponding to direct adsorption of CO to inactive surface sites, and in situ generated CO produced at catalytic active sites. Analysis of the cation-dependent Stark tuning slopes for each of these species provides estimates of the hydrated cation radius upon adsorption to active and inactive sites on the Au electrode. While cations are found to retain their bulk hydration shell upon adsorption at inactive sites, catalytic active sites are characterized by a single layer of water between the Au surface and the electrolyte cation. We propose that the drastic increase in catalytic performance at active sites stems from this unique solvation structure at the Au/electrolyte interface. Building on this evidence of a site-specific EDL structure will be critical to understand the connection between cation-dependent interfacial solvation and CO2R performance.

Oriented External Electric Fields Regurating the Reaction Mechanism of CH4 Oxidation Catalyzed by Fe(IV)-Oxo-Corrolazine: Insight from Density Functional Calculations

Methane is the simplest alkane and can be used as an alternative energy source for oil and coal, but the greenhouse effect caused by its leakage into the air is not negligible, and its conversion into liquid methanol not only facilitates transportation, but also contributes to carbon neutrality. In order to find an efficient method for converting methane to methanol, CH₄ oxidation catalyzed by Fe(IV)-Oxo-corrolazine (Fe(IV)-Oxo-Cz) and its reaction mechanism regulation by oriented external electric fields (OEEFs) are systematically studied by density functional calculations. The calculations show that Fe(IV)-Oxo-Cz can abstract one H atom from CH₄ to form the intermediate with OH group connecting on the corrolazine ring, with the energy barrier of 25.44 kcal mol⁻¹. And then the product methanol is formed through the following rebound reaction. Moreover, the energy barrier can be reduced to 20.72 kcal mol⁻¹ through a two-state reaction pathway. Furthermore, the effect of OEEFs on the reaction is investigated. We found that OEEFs can effectively regulate the reaction by adjusting the stability of the reactant and the transition state through the interaction of electric field-molecular dipole moment. When the electric field is negative, the energy barrier of the reaction decreases with the increase of electric intensity. Moreover, the OEEF aligned along the intrinsic Fe‒O reaction axis can effectively regulate the ability of forming the OH on the corrolazine ring by adjusting the charges of O and H atoms. When the electric field intensity is −0.010 a.u., the OH can be directly rebounded to the CH₃· before it is connecting on the corrolazine ring, thus forming the product directly from the transition state without passing through the intermediate with only an energy barrier of 17.34 kcal mol⁻¹, which greatly improves the selectivity of the reaction.

TUPÃ: Electric field analyses for molecular simulations

We introduce TUPÃ, a Python-based algorithm to calculate and analyze electric fields in molecular simulations. To demonstrate the features in TUPÃ, we present three test cases in which the orientation and magnitude of the electric field exerted by biomolecules help explain biological phenomena or observed kinetics. As part of TUPÃ, we also provide a PyMOL plugin to help researchers visualize how electric fields are organized within the simulation system. The code is freely available and can be obtained at https://mdpoleto.github.io/tupa/.

Vibrational Stark Effects: Ionic Influence on Local Fields

Molecules containing vibrational Stark shift reporters provide a useful tool for measuring DC electric fields in situ. To quantify this effect theoretically, density functional theory (DFT) calculations are usually utilized in a uniform electric field. However, using a combined theoretical and experimental study, we demonstrate here that uniform field DFT cannot simultaneously model the behavior of the three strongest vibrational modes in molecules forming a monolayer on an electrode. We show, by directly modeling ionic movement, that the measured Stark shifts are explained by partial electrical double-layer penetration into the molecular layer. This effect is sensitive to the local environment, and the Stark shifts can be fully suppressed experimentally by introducing a mixed molecular layer that prevents ionic double-layer penetration.

A two-directional vibrational probe reveals different electric field orientations in solution and an enzyme active site

The catalytic power of an electric field depends on its magnitude and orientation with respect to the reactive chemical species. Understanding and designing new catalysts for electrostatic catalysis thus requires methods to measure the electric field orientation and magnitude at the molecular scale. We demonstrate that electric field orientations can be extracted using a two-directional vibrational probe by exploiting the vibrational Stark effect of both the C=O and C-D stretches of a deuterated aldehyde. Combining spectroscopy with molecular dynamics and electronic structure partitioning methods, we demonstrate that, despite distinct polarities, solvents act similarly in their preference for electrostatically stabilizing large bond dipoles at the expense of destabilizing small ones. In contrast, we find that for an active-site aldehyde inhibitor of liver alcohol dehydrogenase, the electric field orientation deviates markedly from that found in solvents, which provides direct evidence for the fundamental difference between the electrostatic environment of solvents and that of a preorganized enzyme active site.

Alkyne-Tagged Raman Probes for Local Environmental Sensing by Hydrogen-Deuterium Exchange

Alkyne-tagged Raman probes have shown high promise for noninvasive and sensitive visualization of small biomolecules to understand their functional roles in live cells. However, the potential for alkynes to sense cellular environments that goes beyond imaging remains to be further explored. Here, we report a general strategy for Raman imaging-based local environment sensing by hydrogen-deuterium exchange (HDX) of terminal alkynes (termed alkyne-HDX). We first demonstrate, in multiple Raman probes, that deuterations of the alkynyl hydrogens lead to remarkable shifts of alkyne Raman peaks for about 130 cm^-1, providing resolvable signals suited for imaging-based analysis with high specificity. Both our analytical derivation and experimental characterizations subsequently establish that HDX kinetics are linearly proportional to both alkyne pK_as and environmental pDs. After validating the quantitative nature of this strategy, we apply alkyne-HDX to sensing local chemical and cellular environments. We establish that alkyne-HDX exhibits high sensitivity to various DNA structures and demonstrates the capacity to detect DNA structural changes in situ from UV-induced damage. We further show that this strategy is also applicable to resolve subtle pD variations in live cells. Altogether, our work lays the foundation for utilizing alkyne-HDX strategy to quantitatively sense the local environments for a broad spectrum of applications in complex biological systems.

Nitrile Infrared Intensities Characterize Electric Fields and Hydrogen Bonding in Protic, Aprotic, and Protein Environments

Nitriles are widely used vibrational probes; however, the interpretation of their IR frequencies is complicated by hydrogen bonding (H-bonding) in protic environments. We report a new vibrational Stark effect (VSE) that correlates the electric field projected on the -C≡N bond to the transition dipole moment and, by extension, the nitrile peak area or integrated intensity. This linear VSE applies to both H-bonding and non-H-bonding interactions. It can therefore be generally applied to determine electric fields in all environments. Additionally, it allows for semiempirical extraction of the H-bonding contribution to the blueshift of the nitrile frequency. Nitriles were incorporated at H-bonding and non-H-bonding protein sites using amber suppression, and each nitrile variant was structurally characterized at high resolution. We exploited the combined information available from variations in frequency and integrated intensity and demonstrate that nitriles are a generally useful probe for electric fields.

Electrostatic vs. inductive effects in phosphine ligand donor properties and reactivity

Enhanced rates and selectivity in enzymes are enabled in part by precisely tuned electric fields within active sites. Analogously, the use of charged groups to leverage electrostatics in molecular systems is a promising strategy to tune reactivity. However, separation of the through space and through bond effects of charged functional groups is a long standing challenge that limits the rational application of electric fields in molecular systems. To address this challenge we developed a method using the phosphorus selenium coupling value (J _P-Se) of anionic phosphine selenides to quantify the electrostatic contribution of the borate moiety to donor strength. In this analysis we report the synthesis of a novel anionic phosphine, PPh₂CH₂BF₃K, the corresponding tetraphenyl phosphonium and tetraethyl ammonium selenides [PPh₄][SePPh₂CH₂BF₃] and [TEA][SePPh₂CH₂BF₃], and the Rh carbonyl complex [PPh₄][Rh(acac)(CO)(PPh₂(CH₂BF₃))]. Solvent-dependent changes in J _P-Se were fit using Coulomb's law and support up to an 80% electrostatic contribution to the increase in donor strength of [PPh₄][SePPh₂CH₂BF₃] relative to SePPh₂Et, while controls with [TEA][SePPh₂CH₂BF₃] exclude convoluting ion pairing effects. Calculations using explicit solvation or point charges effectively replicate the experimental data. This J _P-Se method was extended to [PPh₄][SePPh₂(2-BF₃Ph)] and likewise estimates up to a 70% electrostatic contribution to the increase in donor strength relative to SePPh₃. The use of PPh₂CH₂BF₃K also accelerates C-F oxidative addition reactivity with Ni(COD)₂ by an order of magnitude in comparison to the comparatively donating neutral phosphines PEt₃ and PCy₃. This enhanced reactivity prompted the investigation of catalytic fluoroarene C-F borylation, with improved yields observed for less fluorinated arenes. These results demonstrate that covalently bound charged functionalities can exert a significant electrostatic influence under common solution phase reaction conditions and experimentally validate theoretical predictions regarding electrostatic effects in reactivity.

Protein Electrostatics Investigated through Paramagnetic NMR for Nonpolar Groups

Experimental validation of theoretical models for protein electrostatics remains rare. Recently, we have developed a paramagnetic NMR-based method for de novo determination of effective near-surface electrostatic potentials, which allows for straightforward examination of electrostatic models for biomolecules. In the current work, we expand this method and demonstrate that effective near-surface electrostatic potentials can readily be determined from ¹H paramagnetic relaxation enhancement (PRE) data for protein C_αH and CH₃ groups. The experimental data were compared with those predicted from the Poisson-Boltzmann theory. The impact of structural dynamics on the effective near-surface electrostatic potentials was also assessed. The agreement between the experimental and theoretical data was particularly good for methyl ¹H nuclei. Compared to the conventional pK_a-based validation, our paramagnetic NMR-based approach can provide a far larger number of experimental data that can directly be used to examine the validity of theoretical electrostatic models for proteins.

Origin of thiocyanate spectral shifts in water and organic solvents

Vibrational spectroscopy is a useful technique for probing chemical environments. The development of models that can reproduce the spectra of nitriles and azides is valuable because these probes are uniquely suited for investigating complex systems. Empirical vibrational spectroscopic maps are commonly employed to obtain the instantaneous vibrational frequencies during molecular dynamics simulations but often fail to adequately describe the behavior of these probes, especially in its transferability to a diverse range of environments. In this paper, we demonstrate several reasons for the difficulty in constructing a general-purpose vibrational map for methyl thiocyanate (MeSCN), a model for cyanylated biological probes. In particular, we found that electrostatics alone are not a sufficient metric to categorize the environments of different solvents, and the dominant features in intermolecular interactions in the energy landscape vary from solvent to solvent. Consequently, common vibrational mapping schemes do not cover all essential interaction terms adequately, especially in the treatment of van der Waals interactions. Quantum vibrational perturbation (QVP) theory, along with a combined quantum mechanical and molecular mechanical potential for solute-solvent interactions, is an alternative and efficient modeling technique, which is compared in this paper, to yield spectroscopic results in good agreement with experimental FTIR. QVP has been used to analyze the computational data, revealing the shortcomings of the vibrational maps for MeSCN in different solvents. The results indicate that insights from QVP analysis can be used to enhance the transferability of vibrational maps in future studies.

LModeA-nano: A PyMOL Plugin for Calculating Bond Strength in Solids, Surfaces, and Molecules via Local Vibrational Mode Analysis

The analysis of chemical bonding in crystal structures and surfaces is an important research topic in theoretical chemistry. In this work, we present a PyMOL plugin, named LModeA-nano, as implementation of the local vibrational mode theory for periodic systems (Tao et al. J. Chem. Theory Comput. 2019, 15, 1761) assessing bond strength in terms of local stretching force constants in extended systems of one, two, and three dimensions. LModeA-nano can also analyze chemical bonds in isolated molecular systems thus enabling a head-to-head comparison of bond strength across systems with different dimensions in periodicity (0-3D). The new code is interfaced to the output generated by various solid-state modeling packages including VASP, CP2K, Quantum ESPRESSO, CASTEP, and CRYSTAL. LModeA-nano is cross-platform, open-source and freely available on GitHub: https://github.com/smutao/LModeA-nano.

The Solvation-Induced Onsager Reaction Field Rather than the Double-Layer Field Controls CO2 Reduction on Gold

The selectivity and activity of the carbon dioxide reduction (CO₂R) reaction are sensitive functions of the electrolyte cation. By measuring the vibrational Stark shift of in situ-generated CO on Au in the presence of alkali cations, we quantify the total electric field present at catalytic active sites and deconvolute this field into contributions from (1) the electrochemical Stern layer and (2) the Onsager (or solvation-induced) reaction field. Contrary to recent theoretical reports, the CO₂R kinetics does not depend on the Stern field but instead is closely correlated with the strength of the Onsager reaction field. These results show that in the presence of adsorbed (bent) CO₂, the Onsager field greatly exceeds the Stern field and is primarily responsible for CO₂ activation. Additional measurements of the cation-dependent water spectra using vibrational sum frequency generation spectroscopy show that interfacial solvation strongly influences the CO₂R activity. These combined results confirm that the cation-dependent interfacial water structure and its associated electric field must be explicitly considered for accurate understanding of CO₂R reaction kinetics.

Vibrational Stark fields in carboxylic acid dimers

Carboxylic acids form exceptionally stable dimers and have been used to model proton and double proton transfer processes. The stabilization energies of the carboxylic acid dimers are very weakly dependent on the nature of substitution. However, the electric field experienced by the OH group of a particular carboxylic acid is dependent more on the nature of the substitution on the dimer partner. In general, the electric field was higher when the partner was substituted with an electron-donating group and lower with an electron-withdrawing substituent on the partner. The Stark tuning rate (Δ) of the O-H stretching vibrations calculated at the MP2/aug-cc-pVDZ level was found to be weakly dependent on the nature of substitution on the carboxylic acid. The average Stark tuning rate of O-H stretching vibrations of a particular carboxylic acid when paired with other acids was 5.7 cm^-1 (MV cm^-1)^-1, while the corresponding average Stark tuning rate of the partner acids due to a particular carboxylic acid was 21.9 cm^-1 (MV cm^-1)^-1. The difference in the Stark tuning rate is attributed to the primary and secondary effects of substitution on the carboxylic acid. The average Stark tuning rate for the anharmonic O-D frequency shifts is about 40-50% higher than the corresponding harmonic O-D frequency shifts calculated at the B3LYP/aug-cc-pVDZ level, much greater than the typical scaling factors used, indicating the strong anharmonicity of O-H/O-D oscillators in carboxylic acid dimers. Finally, the linear correlation observed between pK_a and the electric field was used to estimate the pK_a of fluoroformic acid to be around 0.9.