Search NASASearch

DOE OSTI · 3400818

Analytic Nuclear Gradients Including Oriented External Electric Fields in a Molecule-Fixed Frame

Abstract

Electric-field-assisted chemistry has attracted much attention in recent years, particularly in the context of oriented external electric fields for controlling molecular structure and reactivity. Such fields have been explored in a wide range of applications, including switching materials, nanoparticles, controllable catalysts, medicines, and clinical therapies. However, the determination of fixed fields in the laboratory frame becomes ineffective for flexible molecules, as conformational changes can significantly alter the relative orientation between the applied field and molecular structure. In this work, we propose two molecular reference frames─the principal axis frame and the local reference frame─to define oriented electric fields within the molecular framework. These coordinate systems powerfully eliminate ambiguities in the relative orientation between the applied field and the molecule. Analytic nuclear gradients in the presence of external electric fields are derived and implemented, with an initial application to field-dependent geometry optimizations of cis - and trans -formanilide. Analysis of the resulting field-induced equilibrium structures reveals distinct structural responses, validating the accuracy and robustness of the proposed formalism. The analytic gradient framework enables systematic investigations of molecular properties and reactivity under arbitrarily oriented electric fields, opening new opportunities for computational modeling and rational design in electric-field-controlled chemistry.

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Lai, Duc Anh [Southern Methodist Univ., Dallas, TX (United States)], Matthews, Devin A. [Southern Methodist Univ., Dallas, TX (United States)] (ORCID:0000000327955483). 2026-07-03. Analytic Nuclear Gradients Including Oriented External Electric Fields in a Molecule-Fixed Frame. https://doi.org/10.1021/acs.jctc.6c00639

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related reports

pyEF: A Python Framework for QM and QM/MM Atom-Wise Electric Field Analysis

We introduce pyEF, a software package for computing molecular electric fields, electrostatic interaction energies, and electrostatic potentials from quantum mechanical (QM) atom-centered multipole expansions with atom-wise decomposable contributions. We demonstrate the computational efficiency and accuracy of this QM-derived electric field evaluation tool through several tests. To assess the influence of the underlying QM method and charge partitioning scheme on these electrostatic quantities, we analyze over 250 configurations of an acetone solute molecule in five solvents of variable polarity. We find that electric field calculations are highly sensitive to the choice of charge partitioning method. Even among real-space charge schemes, acetone Stark tuning rates differ by up to a factor of 2. Benchmarking computed solvent dipole moments against experimental bulk values, we conclude that the CM5, ADCH, and Hirshfeld-I charge schemes most reliably capture solvent electrostatics and therefore provide a more faithful foundation for computing electric fields. When constructed from these real-space charges, electric fields are nearly insensitive to basis set size and monotonically increase in magnitude with higher Fock exchange. We also demonstrate efficient convergence of QM electrostatics when more distant molecules are represented solely by MM point charges, reducing computational overhead. Leveraging these findings, we demonstrate the use of pyEF to deduce environmental effects on a transition metal complex from a Ga 4 L 6 12– nanocage and quantify the dominant role of organic linkers in orchestrating electrostatic preorganization.

electric fields

On the formation of edge electric fields in negative triangularity plasmas

The Negative Triangularity (NT) plasma edge is studied with full orbit simulations and compared with similar simulations using a Positive Triangularity (PT) edge geometry. It is found that the thermal edge ions in NT are less well confined than in PT. The edge ion losses in NT set-up an electric field inside the plasma in the loss region. This electric field leads to counter-current edge plasma flows and are consistent with observations in NT plasmas. It could also explain the good edge confinement in NT discharges via the E × B flow shear turbulence suppression mechanism.

electric fields