Two-Way Catalysis in a Diels-Alder Reaction Limits Inhibition Induced by an External Electric Field

Oriented external electric fields (EEFs) act as catalysts that can induce selectivity in chemical reactions. The responses of the Diels-Alder (DA) reaction between butadiene and ethylene (BDE-DA) as well as cyclopentadiene and ethylene (CPDE-DA) towards EEF stimuli are investigated here using density functional theory (B3LYP) calculations. EEF is a vector that catalyzes the reaction in one direction while inhibiting it in the opposite direction. Here we report that the inhibitive direction becomes rate-enhancing after some increase in the EEF. The EEF value that brings about the maximum possible inhibition for the reaction is defined as the electrostatic resistance point (ERP). The possibility of both normal and inverse electron-demand DA reactions causes catalytic activity in both directions of the EEF starting at a unique ERP value. The C5 substituents of cyclopentadiene control the ERP values depending upon the resistance power that the functional group provides against the EEF. The endo and exo diastereomeric transition states of the DA reaction have distinct ERP values and the difference (ΔERP) provides the through-space electrostatic contribution to the stereoselectivity on a relative scale. Thus, the ERP values can be used as a gauge for the electrostatic interactions between substituent groups and external stimuli.

An Ab Initio Neural Network Potential Energy Surface for the Dimer of Formic Acid and Further Quantum Tunneling Dynamics

We construct a full-dimensional ab initio neural network potential energy surface (PES) for the isomerization system of the formic acid dimer (FAD). This is based upon ab initio calculations using the DLPNO-CCSD(T) approach with the aug-cc-pVTZ basis set, performed at over 14000 symmetry-unique geometries. An accurate fit to the obtained energies is generated using a general neural network fitting procedure combined with the fundamental invariant method, and the overall energy-weighted root-mean-square fitting error is about 6.4 cm^-1. Using this PES, we present a multidimensional quantum dynamics study on tunneling splittings with an efficient theoretical scheme developed by our group. The ground-state tunneling splitting of FAD calculated with a four-mode coupled method is in good agreement with the most recent experimental measurements. The PES can be applied for further dynamics studies. The effectiveness of the present scheme for constructing a high-dimensional PES is demonstrated, and this scheme is expected to be feasible for larger molecular systems.

Scoring molecular wires subject to an ultrafast laser pulse for molecular electronic devices

A nonionizing ultrafast laser pulse of 20-fs duration with a peak amplitude electric-field ±E = 200 × 10^-4 a.u. was simulated. It was applied to the ethene molecule to consider its effect on the electron dynamics, both during the application of the laser pulse and for up to 100 fs after the pulse was switched off. Four laser pulse frequencies ω = 0.2692, 0.2808, 0.2830, and 0.2900 a.u. were chosen to correspond to excitation energies mid-way between the (S₁ ,S₂ ), (S₂ ,S₃ ), (S₃ ,S₄ ) and (S₄ ,S₅ ) electronic states, respectively. Scalar quantum theory of atoms in molecules (QTAIM) was used to quantify the shifts of the C1C2 bond critical points (BCPs). Depending on the frequencies ω selected, the C1C2 BCP shifts were up to 5.8 times higher after the pulse was switched off compared with a static E-field with the same magnitude. Next generation QTAIM (NG-QTAIM) was used to visualize and quantify the directional chemical character. In particular, polarization effects and bond strengths, in the form of bond-rigidity vs. bond-flexibility, were found, for some laser pulse frequencies, to increase after the laser pulse was switched off. Our analysis demonstrates that NG-QTAIM, in partnership with ultrafast laser irradiation, is useful as a tool in the emerging field of ultrafast electron dynamics, which will be essential for the design, and control of molecular electronic devices.

Strong External Electric Fields Reduce Explosive Sensitivity: A Theoretical Investigation into the Reaction Selectivity in NH2NO2∙∙∙NH3

Controlling the selectivity of a detonation initiation reaction of explosive is essential to reduce sensitivity, and it seems impossible to reduce it by strengthening the external electric field. To verify this, the effects of external electric fields on the initiation reactions in NH2NO2∙∙∙NH3, a model system of the nitroamine explosive with alkaline additive, were investigated at the MP2/6-311++G(2d,p) and CCSD(T)/6-311++G(2d,p) levels. The concerted effect in the intermolecular hydrogen exchange is characterized by an index of the imaginary vibrations. Due to the weakened concerted effects by the electric field along the −x-direction opposite to the “reaction axis”, the dominant reaction changes from the intermolecular hydrogen exchange to 1,3-intramolecular hydrogen transference with the increase in the field strengths. Furthermore, the stronger the field strengths, the higher the barrier heights become, indicating the lower sensitivities. Therefore, by increasing the field strength and adjusting the orientation between the field and “reaction axis”, not only can the reaction selectivity be controlled, but the sensitivity can also be reduced, in particular under a super-strong field. Thus, a traditional concept, in which the explosive is dangerous under the super-strong external electric field, is theoretically broken. Compared to the neutral medium, a low sensitivity of the explosive with alkaline can be achieved under the stronger field. Employing atoms in molecules, reduced density gradient, and surface electrostatic potentials, the origin of the reaction selectivity and sensitivity change is revealed. This work provides a new idea for the technical improvement regarding adding the external electric field into the explosive system.

How Can Static and Oscillating Electric Fields Serve in Decomposing Alzheimer's and Other Senile Plaques?

Alzheimer's disease is one of the most common neurodegenerative conditions, which are ascribed to extracellular accumulation of β-amyloid peptides into plaques. This phenomenon seems to typify other related neurodegenerative diseases. The present study uses classical molecular-dynamics simulations to decipher the aggregation-disintegration behavior of β-amyloid peptide plaques in the presence of static and oscillating oriented external electric fields (OEEFs). A long-term disintegration of such plaques is highly desirable since this may improve the prospects of therapeutic treatments of Alzheimer's disease and of other neurodegenerative diseases typified by senile plaques. Our study illustrates the spontaneous aggregation of the β-amyloid, its prevention and breakdown when OEEF is applied, and the fate of the broken aggregate when the OEEF is removed. Notably, we demonstrate that the usage of an oscillating OEEF on β-amyloid aggregates appears to lead to an irreversible disintegration. Insight is provided into the root causes of the various modes of aggregation, as well as into the different fates of OEEF-induced disintegration in oscillating vs static fields. Finally, our simulation results are compared to the well-established TTFields and the Deep Brain Stimulation (DBS) therapies, which are currently used options for treatments of Alzheimer's disease and other related neurodegenerative diseases.

Study on the Anti-Aging Performance of Different Nano-Modified Natural Ester Insulating Oils Based on Molecular Dynamics

In order to investigate the anti-aging performance of nano-modified natural ester insulating oils, in this paper, two different types of nanoparticles are selected to modify insulating oils. We studied the microscopic mechanism of nano-modified models using molecular simulation techniques. Three models were established: an oil-water model without the addition of nanoparticles and two which contained nano-Fe₃O₄ and nano-Al₂O₃ particles, where the concentration of water was 1 wt.%. The research found that the diffusion of water molecules in the nano-modified model was slow, and the water molecules generated from transformer insulation aging were adsorbed around the nanoparticles, which inhibited the diffusion of water molecules, reduced the hydrolysis of ester molecules, and effectively enhanced the anti-aging performance of natural ester insulating oil. Compared with two different types of nano-modified models, the interface compatibility between nano-Fe₃O₄ and natural ester insulating oil is better, the composite model is stable, the change rate of the diffusion coefficient with temperature is small, there are more hydrogen bonds generated by nano-Fe₃O₄ and water molecules, and the anti-aging performance of the nano-Fe₃O₄-modified oil model is better.

External Electric Field Effect on the Strength of σ-Hole Interactions: A Theoretical Perspective in Like⋯Like Carbon-Containing Complexes

For the first time, σ-hole interactions within like⋯like carbon-containing complexes were investigated, in both the absence and presence of the external electric field (EEF). The effects of the directionality and strength of the utilized EEF were thoroughly unveiled in the (F-C-F₃)₂, (F-C-H₃)₂, and (H-C-F₃)₂ complexes. In the absence of the EEF, favorable interaction energies, with negative values, are denoted for the (F-C-F₃)₂ and (H-C-F₃)₂ complexes, whereas the (F-C-H₃)₂ complex exhibits unfavorable interactions. Remarkably, the strength of the applied EEF exhibits a prominent role in turning the repulsive forces within the latter complex into attractive ones. The symmetrical nature of the considered like⋯like carbon-containing complexes eradicated the effect of directionality of the EEF. The quantum theory of atoms in molecules (QTAIM), and the noncovalent interaction (NCI) index, ensured the occurrence of the attractive forces, and also outlined the substantial contributions of the three coplanar atoms to the total strength of the studied complexes. Symmetry-adapted perturbation theory (SAPT) results show the dispersion-driven nature of the interactions.

σ-Hole and LP-Hole Interactions of Pnicogen···Pnicogen Homodimers under the External Electric Field Effect: A Quantum Mechanical Study

σ-Hole and lone-pair (lp)-hole interactions within σ-hole···σ-hole, σ-hole···lp-hole, and lp-hole···lp-hole configurations were comparatively investigated on the pnicogen···pnicogen homodimers (PCl₃)₂, for the first time, under field-free conditions and the influence of the external electric field (EEF). The electrostatic potential calculations emphasized the impressive versatility of the examined PCl₃ monomers to participate in σ-hole and lp-hole pnicogen interactions. Crucially, the sizes of σ-hole and lp-hole were enlarged under the influence of the positively directed EEF and decreased in the case of reverse direction. Interestingly, the energetic quantities unveiled more favorability of the σ-hole···lp-hole configuration of the pnicogen···pnicogen homodimers, with significant negative interaction energies, than σ-hole···σ-hole and lp-hole···lp-hole configurations. Quantum theory of atoms in molecules and noncovalent interaction index analyses were adopted to elucidate the nature and origin of the considered interactions, ensuring their closed shell nature and the occurrence of attractive forces within the studied homodimers. Symmetry-adapted perturbation theory-based energy decomposition analysis alluded to the dispersion force as the main physical component beyond the occurrence of the examined interactions. The obtained findings would be considered as a fundamental underpinning for forthcoming studies pertinent to chemistry, materials science, and crystal engineering.

Electromagnetic bioeffects: a multiscale molecular simulation perspective

Electromagnetic bioeffects remain an enigma from both the experimental and theoretical perspectives despite the ubiquitous presence of related technologies in contemporary life. Multiscale computational modelling can provide valuable insights into biochemical systems and predict how they will be perturbed by external stimuli. At a microscopic level, it can be used to determine what (sub)molecular scale reactions various stimuli might induce; at a macroscopic level, it can be used to examine how these changes affect dynamic behaviour of essential molecules within the crowded biomolecular milieu in living tissues. In this review, we summarise and evaluate recent computational studies that examined the impact of externally applied electric and electromagnetic fields on biologically relevant molecular systems. First, we briefly outline the various methodological approaches that have been employed to study static and oscillating field effects across different time and length scales. The practical value of such modelling is then illustrated through representative case-studies that showcase the diverse effects of electric and electromagnetic field on the main physiological solvent - water, and the essential biomolecules - DNA, proteins, lipids, as well as some novel biomedically relevant nanomaterials. The implications and relevance of the theoretical multiscale modelling to practical applications in therapeutic medicine are also discussed. Finally, we summarise ongoing challenges and potential opportunities for theoretical modelling to advance the current understanding of electromagnetic bioeffects for their modulation and/or beneficial exploitation in biomedicine and industry.

Geometrical, electrical, and energetic parameters of hetero-disubstituted cumulenes and polyynes in the presence and absence of the external electric field

Cumulenes and polyynes have the potential to be applied as linear, sp-hybridized, one-dimensional all-carbon nanowires in molecular electronics and optoelectronics. The delocalization and conductivity descriptors of the two π-conjugated systems, heterodisubstituted with the NO2, CN, NH2, and OH groups, were studied using the B3LYP, B3LYP/D3, CAM-B3LYP, and ωB97XD DFT functionals, combined with the aug-cc-pVTZ basis set. Three independent types of molecular descriptors, based on geometry (the HOMA index), electrical properties (trace of the polarizability tensor), and energetic (the HOMO-LUMO energy gap) were shown to be mutually correlated and provided concordant indication that communication through the cumulene chain was considerably better than through the polyyne one. The communication can be tuned by using substituents of significantly different π-electron donor-acceptor properties as well as by the external electric field directed along the carbon chain.

Theoretical investigation into the solvent effect on the thermal decomposition of RDX in tetrahydrofuran, acetone, toluene, and benzene

In order to clarify the solvent effect on the thermal decomposition of explosive, the N-NO₂ trigger-bond strengths and ring strains of RDX (cyclotrimethylenetrinitramine) in its H-bonded complexes with solvent molecules (i.e., tetrahydrofuran, acetone, toluene, and benzene), and the activation energies of the intermolecular hydrogen exchanges between the solvent molecules and C₃H₈O₂N₄ or CH₄O₂N₂, as the model molecule of RDX, were investigated by the BHandHLYP, B3LYP, MP2(full), and M06-2X methods with the 6-311 +  + G(2df,2p) basis set, accompanied by a comparison with the calculations by the integral equation formalism polarized continuum model. The solvent effects ignore the ring strain while strengthening the N-NO₂ bond, leading to a possible decreased sensitivity, as is opposite to the experimental results. However, the activation energies are in the order of C₃H₈O₂N₄/CH₄O₂N₂∙∙∙acetone < C₃H₈O₂N₄/CH₄O₂N₂∙∙∙THF < C₃H₈O₂N₄/CH₄O₂N₂∙∙∙toluene < C₃H₈O₂N₄/CH₄O₂N₂∙∙∙benzene < C₃H₈O₂N₄/CH₄O₂N₂, suggesting that the order of the critical explosion temperatures might be RDX∙∙∙acetone < RDX∙∙∙THF < RDX∙∙∙toluene < RDX∙∙∙benzene < RDX, as is roughly consistent with the experimental results. Therefore, the intermolecular hydrogen exchange with the HONO elimination is a possible mechanism of the solvent effect on the initial thermal decomposition of RDX. The solvent effect on the sensitivity is analyzed by the surface electrostatic potentials.

Hydrogen and Halogen Bonding in Homogeneous External Electric Fields: Modulating the Bond Strengths

Recent years have witnessed various fascinating phenomena arising from the interactions of noncovalent bonds with homogeneous external electric fields (EEFs). Here we performed a computational study to interpret the sensitivity of intrinsic bond strengths to EEFs in terms of steric effect and orbital interactions. The block-localized wavefunction (BLW) method, which combines the advantages of both ab initio valence bond (VB) theory and molecular orbital (MO) theory, and the subsequent energy decomposition (BLW-ED) approach were adopted. The sensitivity was monitored and analyzed using the induced energy term, which is the variation in each energy component along the EEF strength. Systems with single or multiple hydrogen (H) or halogen (X) bond(s) were also examined. It was found that the X-bond strength change to EEFs mainly stems from the covalency change, while generally the steric effect rules the response of H-bonds to EEFs. Furthermore, X-bonds are more sensitive to EEFs, with the key difference between H- and X-bonds lying in the charge transfer interaction. Since phenylboronic acid has been experimentally used as a smart linker in EEFs, switchable sensitivity was scrutinized with the example of the phenylboronic acid dimer, which exhibits two conformations with either antiparallel or parallel H-bonds, thereby, opposite or consistent responses to EEFs. Among the studied systems, the quadruple X-bonds in molecular capsules exhibit remarkable sensitivity, with its interaction energy increased by -95.2 kJ mol^-1 at the EEF strength 0.005 a.u.

Effect of External Electric Field on Tetrel Bonding Interactions in (FTF3···FH) Complexes (T = C, Si, Ge, and Sn)

A quantum chemical study was accomplished on the σ-hole interactions of the barely explored group IV elements, for the first time, in the absence and presence of the positively and negatively directed external electric field (EEF). The analyses of molecular electrostatic potential addressed the occurrence of the σ-hole on all the inspected tetrel atoms, confirming their salient versatility to engage in σ-hole interactions. MP2 energetic findings disclosed the occurrence of favorable σ-hole interactions within the tetrel bonding complexes. The tetrel bonding interactions became stronger in the order of C < Si < Ge < Sn for F-T-F3···FH complexes with the largest interaction energy amounting to -19.43 kcal/mol for the optimized F-Sn-F3···FH complex under the influence of +0.020 au EEF. The interaction energy conspicuously evolved by boosting the magnitude of the positively directed EEF value and declining the negatively directed EEF one. The decomposition analysis for the interaction energies was also executed in terms of symmetry-adapted perturbation theory, illuminating the dominant electrostatic contribution to all the studied complexes' interactions except carbon-based interactions controlled by dispersion forces. The outcomes that emerged from the current work reported significantly how the direction and strength of the EEF affect the tetrel-bonding interactions, leading to further improvements in the forthcoming studies of supramolecular chemistry and materials science.

Kinetic Study of Transition Mutations from G-C to A-T Base Pairs in Watson-Crick DNA Base Pairs: Double Proton Transfers

According to the Löwdin model [ Rev. Mod. Phys. 1963, 35, 724-732], the Watson-Crick guanine-cytosine (G-C) base pair is tautomerized (G*-C*) with a small probability and then replication of G*-C* produces G*-thymine (T) and adenine (A)-C* base pairs. On the basis of this model and our previous work [ J. Phys. Chem. B 2020, 124, 1715-1722], we first calculated the intrinsic reaction coordinates from G*-T to G-T* using density functional theory and evaluated the probability of G*-T tautomerization to G-T* by double proton transfer (DPT) on the basis of the transition state theory. Similarly, we calculated the probability of A-C* tautomerization to A*-C by DPT. Then, according to these probabilities, we calculated the probability of transition mutations from G-C to A-T after 2 replications. The calculated probability was 1.31 × 10^-8, a value consistent with the mutation rate previously reported by Drake et al. [ Proc. Natl. Acad. Sci. U.S.A. 1991, 88, 7160-7164]. Our results suggest that DPT is one cause of the G-C → A-T transition. To investigate differences in the optical properties between G*-T and G-T* and between A-C* and A*-C, we also evaluated the infrared absorption spectra and Raman intensities for these base pairs.

Electron ratcheting in self-assembled soft matter

Ratcheted multi-step hopping electron transfer systems can plausibly produce directional charge transport over very large distances without requiring a source-drain voltage bias. We examine molecular strategies to realize ratcheted charge transport based on multi-step charge hopping, and we illustrate two ratcheting mechanisms with examples based on DNA structures. The charge transport times and currents that may be generated in these assemblies are also estimated using kinetic simulations. The first ratcheting mechanism described for nanoscale systems requires local electric fields on the 10⁹ V/m scale to realize nearly 100% population transport. The second ratcheting mechanism for even larger systems, based on electrochemical gating, is estimated to generate currents as large as 0.1 pA for DNA structures that are a few μm in length with a gate voltage of about 5 V, a magnitude comparable to currents measured in DNA wires at the nanoscale when a source-drain voltage bias of similar magnitude is applied, suggesting an approach to considerably extend the distance range over which DNA charge transport devices may operate.

The influence of external electric fields on proton transfer tautomerism in the guanine-cytosine base pair

The Watson-Crick base pair proton transfer tautomers would be widely considered as a source of spontaneous mutations in DNA replication if not for their short lifetimes and thermodynamic instability. This work investigates the effects external electric fields have on the stability of the guanine-cytosine proton transfer tautomers within a realistic strand of aqueous DNA using a combination of ensemble-based classical molecular dynamics (MD) coupled to quantum mechanics/molecular mechanics (QM/MM). Performing an ensemble of calculations accounts for the stochastic aspects of the simulations while allowing for easier identification of systematic errors. The methodology applied in this work has previously been shown to estimate base pair proton transfer rate coefficients that are in good agreement with recent experimental data. A range of electric fields in the order of 104 to 109 V m-1 is investigated based on their real-life medicinal applications which include gene therapy and cancer treatments. The MD trajectories confirm that electric fields up to 1.00 × 109 V m-1 have a negligible influence on the structure of the base pairs within DNA. The QM/MM results show that the application of large external electric fields (1.00 × 109 V m-1) parallel to the hydrogen bonds increases the thermodynamic population of the tautomers by up to one order of magnitude; moreover, the lifetimes of the tautomers remain insignificant when compared to the timescale of DNA replication.

External electric field effects on the σ-hole and lone-pair hole interactions of group V elements: a comparative investigation

σ-hole and lone-pair (lp) hole interactions of trivalent pnicogen-bearing (ZF3) compounds were comparatively scrutinized, for the first time, under field-free and external electric field (EEF) conditions. Conspicuously, the sizes of the σ-hole and lp-hole were increased by applying an EEF along the positive direction, while the sizes of both holes decreased through the reverse EEF direction. The MP2 energetic calculations of ZF3⋯FH/NCH complexes revealed that σ-holes exhibited more impressive interaction energies compared to the lp-holes. Remarkably, the strengths of σ-hole and lp-hole interactions evolved with the increment of the positive value of the considered EEF; i.e., the interaction energy increased as the utilized EEF value increased. Unexpectedly, under field-free conditions, nitrogen-bearing complexes showed superior strength for their lp-hole interactions than phosphorus-bearing complexes. However, the reverse picture was exhibited for the interaction energies of nitrogen- and phosphorus-bearing complexes interacting within lp-holes by applying the high values of a positively directed EEF. These results significantly demonstrate the crucial influence of EEF on the strength of σ-hole and lp-hole interactions, which in turn leads to an omnipresent enhancement for variable fields, including biological simulations and material science.

Switch chemistry at cryogenic conditions: quantum tunnelling under electric fields

While the influence of intramolecular electric fields is a known feature in enzymes, the use of oriented external electric fields (EEF) to enhance or inhibit molecular reactivity is a promising topic still in its infancy. Herein we will explore computationally the effects that EEF can provoke in simple molecules close to the absolute zero, where quantum tunnelling (QT) is the sole mechanistic option. We studied three exemplary systems, each one with different reactivity features and known QT kinetics: π bond-shifting in pentalene, Cope rearrangement in semibullvalene, and cycloreversion of diazabicyclohexadiene. The kinetics of these cases depend both on the field strength and its direction, usually giving subtle but remarkable changes. However, for the cycloreversion, which suffers large changes on the dipole through the reaction, we also observed striking results. Between the effects caused by the EEF on the QT we observed an inversion of the Arrhenius equation, deactivation of the molecular fluxionality, and stabilization or instantaneous decomposition of the system. All these effects may well be achieved, literally, at the flick of a switch.

Adding an external electric field to reactions driven by quantum mechanical tunneling brings a whole new dimension to the idea of switch chemistry.

External electric field reduces the explosive sensitivity: a theoretical investigation into the hydrogen transference kinetics of the NH2NO2∙∙∙H2O complex

Controlling the selectivity of detonation initiation reaction to reduce the explosive sensitivity has been a Holy Grail in the field of energetic materials. The effects of the external electric fields on the homolysis of the N-NO₂ bond and initiation reaction dynamics of NH₂NO₂∙∙∙H₂O (i.e., intermolecular and 1,3-intramolecular hydrogen transfers) were investigated at the MP2/6-311++G(2d,p) and CCSD/6-311++G(2d,p)//MP2/6-311++G(2d,p) levels. The results show that the N-NO₂ bond is not the "trigger linkage." The notable transiliences of the activation energy of the intermolecular hydrogen transfer are found with the field strength of - 0.012 a.u. along the -x-direction, leading to the conversion of the main reaction between the intermolecular and 1,3-intramolecular hydrogen transference. The activation energies of two kinds of the hydrogen transferences are increased under the external electric fields along the -y-direction. In particular, due to the conversion of the main reaction, the activation energies of the overall reaction are increased significantly along the -x-direction, leading to the significant reduced explosive sensitivities. Therefore, by controlling the field strengths and orientations between the "reaction axis" and external electric field along the y- and x-directions, the selectivity of the initiation reaction could be controlled and the explosive sensitivity could be reduced. Employing AIM (atoms in molecules) and surface electrostatic potentials, the origin of the control of reaction selectivity and the reduction of sensitivity is revealed. This work is of great significance to the improvement of the technology that the external electric fields are added safely into the energetic material system to enhance the explosive performance. Graphical abstract.

Electric field induced intra-molecular self-redox: superalkali Li3N3Mg as a candidate for NLO molecular switches

A novel intra-molecular self-redox switch, Li₃N₃Mg, is constructed theoretically. Our investigation showed that a suitably oriented external electric field (OEEF) can drive a long-range excess electron transfer from Mg atoms to Li₃ rings. And subsequently, an interesting intra-molecular self-redox from Li₃²⁺N₃^3-Mg⁺ to Li₃⁺N₃^3-Mg²⁺ accompanying the large different electronic static first hyperpolarizability (β) is exhibited. The increase of the β value constitutes an order of magnitude improvement from Li₃²⁺N₃^3-Mg⁺ (34 986 a.u.) to Li₃⁺N₃^3-Mg²⁺ (101 225 a.u.), which indicates that Li₃N₃Mg is a good candidate for a self-redox NLO molecular switch.

Efficient Quantum Mechanical Calculations of Mode-Specific Tunneling Splittings upon Fundamental Excitation in the Dimer of Formic Acid

The formic acid dimer (FAD) is an important benchmark system for understanding the double hydrogen transfer process. Most recently, Zhang et al. measured a few tunneling splittings upon fundamental excitation of FAD precisely (Zhang, Y. et al. J. Chem. Phys. 2017, 146, 244306); however, relevant theoretical studies are very limited. Here, we present a multidimensional quantum dynamics study on mode-specific tunneling splittings upon fundamental excitation in FAD with an efficient theoretical scheme developed by our group in which the process-oriented basis function customization strategy is combined with the preconditioned inexact spectral transform method. Various mode-specific tunneling splittings upon fundamental excitation are systematically calculated, and interesting mode-specific excitation effects on tunneling rate are identified. In particular, the calculated tunneling splittings for the ν₂₂ and ν₂₁ states are in good agreement with experiment, and the remarkable mode-specific suppression effects upon excitation should result from that the antisymmetric vibrational modes hinder the concerted double H-transfer. The present work is helpful to acquire a better understanding of the mode-specific excitation effects on tunneling processes.

Ordered Solvents and Ionic Liquids Can Be Harnessed for Electrostatic Catalysis

Herein, we employ classical molecular dynamics simulations using the Drude oscillator-based polarizable force field, quantum chemical calculations, and ONIOM multiscale calculations to study (a) how an external field orders the solvent environment in a chemical reaction and then (b) whether in the absence of this same applied field the ordered solvent environment alone can electrostatically catalyze a chemical reaction when compared with the corresponding disordered solvent. Our results show that a 0.2 V/Å external electric field, which is below the threshold for bond breaking of solvent molecules, leads to significant ordering of bulk methanol solvent and the ionic liquid [EMIM][BF₄]. Importantly, in the absence of this same field, the ordered solvent lowers the activation energy of the hydrogen-transfer reaction of o-alkylphenyl ketones in excess of 20 kcal/mol when the solvent is methanol and by over 30 kcal/mol for [EMIM][BF₄]. Even a 0.1 V/Å external field has effects of ca. 10 and 20 kcal/mol, respectively. This work suggests a possible strategy for scaling electrostatic catalysis by applying a pulsed external field to the reaction medium to maintain solvent ordering while allowing the reaction to proceed largely in the absence of an external field.

Solvent Organization and Rate Regulation of a Menshutkin Reaction by Oriented External Electric Fields are Revealed by Combined MD and QM/MM Calculations

When and how do external electric fields (EEFs) lead to catalysis in the presence of a (polar or nonpolar) solvent? This is the question that is addressed here using a combination of molecular dynamics (MD) simulations, quantum mechanical/molecular mechanical calculations with EEF, and quantum mechanical/(local) electric field calculations. The paper focuses on a model reaction, the Menshutkin reaction between CH3I and pyridine in three solvents of varying polarity. Using MD simulations, we find that the EEF causes the solvent to undergo organization; the solvent molecules gradually align with the applied field as the field strength increases. The collective orientation of the solvent molecules modifies the electrostatic environment around the Menshutkin species and induces a global electric field pointing in the opposite direction of the applied EEF. The combination of these two entangled effects leads to partial or complete screening of the EEF, with the extent of screening being proportional to the polarity/polarizability of the solvent. Nevertheless, we find that catalysis of the Menshutkin reaction inevitably emerges once the EEF exceeds the opposing field of the organizing solvent, i.e., once polarization of the Menshutkin complex is observed to set in. Overall, our analysis provides a lucid and pictorial interpretation of the behavior of solutions in the presence of EEFs and indicates that EEF-mediated catalysis should, in principle, be feasible in bulk setups, especially for nonpolar and mildly polar solvents. By application of the charge-transfer paradigm, it is shown that the emergence of OEEF catalysis in solution can be generalized to other reactions as well.

Manipulation of Diatomic Molecules with Oriented External Electric Fields: Linear Correlations in Atomic Properties Lead to Nonlinear Molecular Responses

Oriented external electric fields (OEEFs) have been shown to have great potential in being able to provide unprecedented control of chemical reactions, catalysis, and selectivity with applications ranging from H₂ storage to molecular machines. We report a theoretical study of the atomic origins of molecular changes because of OEEFs since understanding the characteristics of OEEF-induced couplings between atomic and molecular properties is an important step toward comprehensive understanding of the effects of strong external fields on the molecular structure, stability, and reactivity. We focus on the atomic and molecular (bond) properties of a set of homo- (H₂, N₂, O₂, F₂, and Cl₂) and heterodiatomic (HF, HCl, CO, and NO) molecules under intense external electric fields in the context of quantum theory of atoms in molecules (QTAIM). It is shown that the atomic properties (atomic charges, energies, and localization indices) correlate linearly with the field strengths, but molecular properties (bond length, electron density at the bond critical point, and electron delocalization index) exhibit nonlinear responses to the imposed fields. In particular, the changes in the electron density distribution alter the shapes and locations of the zero-flux surfaces, atomic volumes, atomic electron population, and localization/delocalization indices. The topography and topology of the molecular electrostatic potential undergo dramatic changes. External fields also perturb the covalent-polar-ionic characteristic of the studied chemical bonds, hallmarking the impact of electric fields on the stability and reactivity of chemical compounds. The findings are well-rationalized within the framework of the QTAIM and form a coherent conceptual understanding of these effects in prototypical diatomic molecules.

Next-generation QTAIM for scoring molecular wires in E-fields for molecular electronic devices

The effect of a varying, directional E _x , E _y , and E _z electric field on the ethene molecule was investigated using next-generation quantum theory of atoms in molecules (QTAIM). Despite using low E-field strengths that are within the realm of experiment and do not measurably alter the molecular geometry, significant changes to the QTAIM properties were observed. Using conventional QTAIM, the shifting of the C─C and C─H bond critical points (BCPs) demonstrates polarization through an interchange in the size of the atoms involved in a bond, since a BCP is located on the boundary between a pair of bonded atoms. Next-generation QTAIM, however, demonstrates the polarization effect more directly with a change in morphology of the 3-D envelope around the BCP. Modest increases of ≈ 2% in the ellipticity ε of the BCP were uncovered when the C─C bond was aligned parallel or anti-parallel to the applied E _x -field. Significant asymmetries were found in the response of the next-generation QTAIM 3-D paths of the C─H bonds to the applied E-field. When the E-field coincided with the C─C bond, the BCP moved in response and was accompanied by the envelope constructed from 3-D next-generation paths. The response displayed a polarization effect that increased with increasing magnitude of the E _x -field parallel and anti-parallel to the C─C bond. Our analysis demonstrates that next-generation QTAIM is a useful tool for understanding the response of molecules to E-fields, for example, for the screening of molecular wires for the design of molecular electronic devices. © 2019 Wiley Periodicals, Inc.

Fullerene-Based Switching Molecular Diodes Controlled by Oriented External Electric Fields

Employing multiscale in silico modeling, we propose switching molecular diodes on the basis of endohedral fullerenes (fullerene switching diode, FSD), encapsulated with polar molecules of general type MX (M: metal, X: nonmetal) to be used for data storage and processing. Here, we demonstrate for MX@C₇₀ systems that the relative orientation of enclosed MX with respect to a set of electrodes connected to the system can be controlled by application of oriented external electric field(s). We suggest systems with two- and four-terminal electrodes, in which the source and drain electrodes help the current to pass through the device and help the switching between the conductive states of FSD via applied voltage. The gate electrodes then assist the switching by effectively lowering the energy barrier between local minima via stabilizing the transition state of switching process if the applied voltage between the source and drain is insufficient to switch the MX inside the fullerene. Using nonequilibrium Green's function combined with density functional theory (DFT-NEGF) computations, we further show that conductivity of the studied MX@C₇₀ systems depends on the relative orientation of MX inside the cage with respect to the electrodes. Therefore, the orientation of the MX inside C₇₀ can be both enforced ("written") and retrieved ("read") by applied voltage. The studied systems thus behave like voltage-sensitive switching molecular diodes, which is reminiscent of a molecular memristor.

Double Proton Transfer in the Dimer of Formic Acid: An Efficient Quantum Mechanical Scheme

Double proton transfer plays an important role in biology and chemistry, such as with DNA base pairs, proteins and molecular clusters, and direct information about these processes can be obtained from tunneling splittings. Carboxylic acid dimers are prototypes for multiple proton transfer, of which the formic acid dimer is the simplest one. Here, we present efficient quantum dynamics calculations of ground-state and fundamental excitation tunneling splittings in the formic acid dimer and its deuterium isotopologues. These are achieved with a multidimensional scheme developed by us, in which the saddle-point normal coordinates are chosen, the basis functions are customized for the proton transfer process, and the preconditioned inexact spectral transform method is used to solve the resultant eigenvalue problem. Our computational results are in excellent agreement with the most recent experiments (Zhang et al., 2017; Li et al., 2019).

Insights into the activation process of CO2 through Dihydrogenation reaction

Based on first principle calculation, activation of CO₂ has been analyzed thoroughly by using different conceptual density functional theory based descriptors like reaction force, reaction force constant, reaction electronic flux, dual descriptor, etc. via dihydrogenation reaction of B₃N₃, H₂ and CO₂. The total reaction is a two-step reaction where initially B₃N₃H₂ is formed from the reaction between B₃N₃ and H₂ and in the second step HCOOH is form due to the reaction of CO₂ by B₃N₃H₂. It has been found that the di-hydrogen reaction for the CO₂ activation is endothermic in nature, which can be changed to exothermic reaction by applying proper external electric field. Movement of H₂ plays an important role in the CO₂ activation process. The reaction force constant, Wiberg bond index and its derivative reveal that the reaction is slightly asynchronous and concerted in nature.

A dynamic and electrostatic potential prediction of the prototropic tautomerism between imidazole 3-oxide and 1-hydroxyimidazole in external electric field

In order to obtain an optimum scheme for separating the proton-transfer tautomer, a dynamic investigation into the effect of the external electric field on the proton-transfer tautomeric conversion in imidazole 3-oxide and 1-hydroxyimidazole was carried out at the M06-2X/6-311++G** and CCSD(T)/6-311++G(2d,p) level, accompanied by the analysis of the surface electrostatic potentials. The results show that, for both the forward reaction "imidazole 3-oxide → N-hydroxyimidazole free radical → 1-hydroxyimidazole" and its reverse reaction processes, the fields parallel to the N→O or N-OH bond axis affect the barrier heights and rate constants considerably more than those parallel to the other orientations. As the field strength is increased along the orientation from the O to N atom, the chemical equilibrium moves toward the direction for the formation of 1-hydroxyimidazole, while the amount of imidazole 3-oxide is increased with the increased field strength along the opposite orientation. In the fields along the orientation consistent with the dipole moment, the electrostatic potentials and their variances "abnormally" increase for the transition states with the N→O bond in comparison with those in no field (they decrease generally), which enhances the nucleophilicity of the coordination O atom and the electrophilicity of the activated H atom. The analyses of the AIM (atoms in molecules) and NICS (nucleus-independent chemical shift) were used to explain the above anomaly. Graphical Abstract Electrostatic potentials and their variances "abnormally" increase in the external electric field, which greatly affects tautomeric conversion.

The Role of Proton Transfer on Mutations

Hydrogen bonds play a critical role in nucleobase studies as they encode genes, map protein structures, provide stability to the base pairs, and are involved in spontaneous and induced mutations. Proton transfer mechanism is a critical phenomenon that is related to the acid-base characteristics of the nucleobases in Watson-Crick base pairs. The energetic and dynamical behavior of the proton can be depicted from these characteristics and their adjustment to the water molecules or the surrounding ions. Further, new pathways open up in which protonated nucleobases are generated by proton transfer from the ionized water molecules and elimination of a hydroxyl radical in this review, the analysis will be focused on understanding the mechanism of untargeted mutations in canonical, wobble, Hoogsteen pairs, and mutagenic tautomers through the non-covalent interactions. Further, rare tautomer formation through the single proton transfer (SPT) and the double proton transfer (DPT), quantum tunneling in nucleobases, radiation-induced bystander effects, role of water in proton transfer (PT) reactions, PT in anticancer drugs-DNA interaction, displacement and oriental polarization, possible models for mutations in DNA, genome instability, and role of proton transfer using kinetic parameters for RNA will be discussed.

The destabilization of hydrogen bonds in an external E-field for improved switch performance

The effect of an electric field on a recently proposed molecular switch based on a quinone analogue was investigated using next-generation quantum theory of atoms in molecules (QTAIM) methodology. The reversal of a homogenous external electric field was demonstrated to improve the "OFF" functioning of the switch. This was achieved by destabilization of the H atom participating in the tautomerization process along the hydrogen bond that defines the switch. The "ON" functioning of the switch, from the position of the tautomerization barrier, is also improved by the reversal of the homogenous external electric field: this result was previously inaccessible. The "ON" and "OFF" functioning of the switch was visualized in terms of the response of the most preferred directions of motion of the electronic charge density to the applied external field. All measures from QTAIM and the stress tensor provide consistent results for the factors affecting the "ON" and "OFF" switch performance. Our analysis therefore demonstrates use for future design of molecular electronic devices. © 2019 Wiley Periodicals, Inc.

Electrophilic Aromatic Substitution Reactions: Mechanistic Landscape, Electrostatic and Electric-Field Control of Reaction Rates, and Mechanistic Crossovers

This study investigates the rich mechanistic landscape of the iconic electrophilic aromatic substitution (EAS) reaction class, in the gas phase, in solvents, and under stimulation by oriented external electric fields. The study uses DFT calculations, complemented by a qualitative valence bond (VB) perspective. We construct a comprehensive and unifying framework that elucidates the many surprising mechanistic features, uncovered in recent years, of this class of reactions. For example, one of the puzzling issues which have attracted significant interest recently is the finding of a variety of concerted mechanisms that do not involve the formation of σ-complex intermediates, in apparent contradiction to the generally accepted textbook mechanism. Our VB modeling elucidates the existence of both the concerted and stepwise mechanisms and uncovers the root causes and necessary conditions for the appearance of these intermediates. Furthermore, our VB analysis offers insight into the potential applications of external electric fields as smart, green, and selective catalysts, which can control at will reaction rates, as well as mechanistic crossovers, for this class of reactions. Finally, we highlight how understanding of the electric fields effect on the EAS reaction could lead to the formulation of guiding principles for the design of improved heterogeneous catalysts. Overall, our analysis underscores the powerful synergy offered by combining molecular orbital and VB theory to tackle interesting and challenging mechanistic questions in chemistry.

Nonlinear Migration Dynamics of Excess Electrons along Linear Oligopeptides Controlled by an Applied Electric Field

Migration of an excess electron along linear oligopeptides governed by the external electric field (E_ex ) which is against the inner dipole electric field is theoretically investigated, including the effects of E_ex on the structural and electronic properties of electron migration. Two structural properties including electron-binding ability and the dipole moment of linear oligopeptides are sensitive to the E_ex values and can be largely modulated by E_ex due to the competition of E_ex and the inner electric field and electron transfer caused by E_ex . In the case of low E_ex values, two structural properties decrease slightly_, while for high E_ex values, the electron-binding ability continually increases strongly, with dipole moments firstly increasing significantly and then increasing more slowly at higher E_ex . Additionally, linear oligopeptides of different chain lengths influence the modulation extent of E_ex and the longer the chain length is, the more sensitive modulation of E_ex is. In addition, electronic properties represented by electron spin densities and singly occupied molecular orbital distributions vary with E_ex intensities, leading to an unusual electron migration behavior. As E_ex increases, an excess electron transfers from the N-terminus to the C-terminus and jumps over a neighboring dipole unit of two termini to other units, respectively, instead of transferring by means of a one-by-one dipole unit hopping mechanism. These findings not only promote a deeper understanding of the connection between E_ex and structural and electronic properties of electron transfer behavior in peptides, but also provide a new insight into the modulation of electron migration along the oligopeptides.

Oriented External Electric Fields: Tweezers and Catalysts for Reactivity in Halogen-Bond Complexes

This theoretical study establishes ways of controlling and enabling an uncommon chemical reaction, the displacement reaction, B:---(X-Y) → (B-X)⁺ + :Y^-, which is nascent from a B:---(X-Y) halogen bond (XB) by nucleophilic attack of the base, B:, on the halogen, X. In most of the 14 cases examined, these reactions possess high barriers either in the gas phase (where the X-Y bond dissociates to radicals) or in solvents such as CH₂Cl₂ and CH₃CN (which lead to endothermic processes). Thus, generally, the XB species are trapped in deep minima, and their reactions are not allowed without catalysis. However, when an oriented-external electric field (OEEF) is directed along the B---X---Y reaction axis, the field acts as electric tweezers that orient the XB along the field's axis, and intensely catalyze the process, by tens of kcal/mol, thus rendering the reaction allowed. Flipping the OEEF along the reaction axis inhibits the reaction and weakens the interaction of the XB. Furthermore, at a critical OEEF, each XB undergoes spontaneous and barrier-free reaction. As such, OEEF achieves quite tight control of the structure and reactivity of XB species. Valence bond modeling is used to elucidate the means whereby OEEFs exert their control.