Analytic nuclear gradients including oriented external electric fields in a molecule-fixed frame

Devin A. Matthews; Duc Anh Lai

arxiv: 2604.01189 · v2 · submitted 2026-04-01 · ⚛️ physics.chem-ph

Analytic nuclear gradients including oriented external electric fields in a molecule-fixed frame

Duc Anh Lai , Devin A. Matthews This is my paper

Pith reviewed 2026-05-13 21:38 UTC · model grok-4.3

classification ⚛️ physics.chem-ph

keywords analytic nuclear gradientsoriented external electric fieldsmolecule-fixed framesprincipal axis framelocal reference framegeometry optimizationformanilideelectric field effects

0 comments

The pith

Oriented external electric fields defined in molecule-fixed frames enable analytic nuclear gradients.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper derives analytic nuclear gradients for molecules under oriented external electric fields by introducing two molecular reference frames: the principal axis frame and the local reference frame. These frames define the field orientation relative to the molecule, overcoming the problem that laboratory-frame fields lose meaning when molecules flex and change shape. The method is implemented and tested through geometry optimizations of cis- and trans-formanilide, which exhibit different structural changes depending on the field direction. A reader should care because it makes computational study of field-controlled chemistry practical for flexible molecules where orientation matters.

Core claim

Analytic nuclear gradients in the presence of external electric fields are derived and implemented using the principal axis frame and the local reference frame to define oriented fields within the molecular framework. Application to field-dependent geometry optimizations of cis- and trans-formanilide reveals distinct structural responses, validating the formalism for systematic investigations of molecular properties under arbitrarily oriented electric fields.

What carries the argument

The principal axis frame and local reference frame, which fix the electric field orientation to the molecule to remove ambiguities from conformational changes.

If this is right

Field-dependent geometry optimizations become possible without orientation ambiguities for flexible molecules.
Distinct structural responses are observed in cis- and trans-formanilide under oriented fields.
The framework supports studies of reactivity and properties under arbitrary field orientations.
New opportunities arise for modeling electric field-controlled chemistry.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

Extending this to molecular dynamics simulations could reveal how fields influence conformational transitions over time.
Applying the method to larger biomolecules might help design field-responsive drugs or catalysts.
The reference frames could be adapted for other external perturbations like magnetic fields.

Load-bearing premise

The principal axis and local reference frames fully remove orientation ambiguities for flexible molecules without adding errors to the analytic gradient formulas.

What would settle it

A calculation on a molecule with known large conformational flexibility where the computed gradient does not match numerical differentiation of the energy under the same oriented field would falsify the claim.

Figures

Figures reproduced from arXiv: 2604.01189 by Devin A. Matthews, Duc Anh Lai.

**Figure 2.** Figure 2: Potential energy surface along the C1 – C2 – N – C3 dihedral angle of cis-formanilide under c-aligned fields in the PAF. 16 [PITH_FULL_IMAGE:figures/full_fig_p016_2.png] view at source ↗

**Figure 3.** Figure 3: Total energies of trans-formanilide under OEEFs in the LRF. [PITH_FULL_IMAGE:figures/full_fig_p020_3.png] view at source ↗

**Figure 4.** Figure 4: Out-of-plane twisting of trans-formanilide phenyl ring due applied [PITH_FULL_IMAGE:figures/full_fig_p021_4.png] view at source ↗

read the original 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.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper gives analytic nuclear gradients for electric fields oriented in principal-axis and local molecular frames, which is useful for flexible molecules, but the gradient expressions may miss chain-rule contributions from the geometry-dependent frame rotations.

read the letter

The core advance is defining the external field direction inside two molecule-fixed frames instead of the lab frame. The principal axis frame comes from the inertia tensor eigenvectors and the local frame uses atom-centered vectors; both keep the field orientation tied to the molecule even when bonds rotate or the shape changes. They then derive and code the analytic nuclear gradients under this setup and run geometry optimizations on cis- and trans-formanilide, where the two isomers show clearly different field-induced distortions. That application is the concrete evidence they offer that the approach works in practice. The analytic route is the practical payoff: it should let people optimize structures under arbitrarily oriented fields without numerical differentiation at every step. The frames themselves are a straightforward way to remove the ambiguity that appears whenever a flexible molecule reorients relative to a fixed lab vector. On the soft side, the principal-axis frame is obtained by diagonalizing a geometry-dependent inertia tensor, so its eigenvectors rotate with the nuclei. The total energy derivative must therefore pick up an extra term from the derivative of that rotation matrix. The same issue can appear in the local frame if it is built from position vectors. The formanilide test case stays asymmetric throughout, so it does not expose whether those chain-rule pieces were included or omitted. The abstract states that gradients were derived and implemented but supplies no explicit formulas, no finite-difference comparisons, and no error tables, which leaves the central claim only partly supported until the full derivation is checked. This is aimed at computational chemists who already run field calculations on non-rigid systems and want consistent orientation during optimization. A reader working on field-controlled catalysis or molecular electronics would see immediate use for the frames and the gradient code. The work is coherent enough on its own terms to deserve a serious referee who can verify the derivative expressions and any numerical validation that appears in the full manuscript.

Referee Report

3 major / 3 minor

Summary. The manuscript derives and implements analytic nuclear gradients for quantum-chemical calculations under oriented external electric fields, where the field direction is defined in molecule-fixed frames: the principal axis frame (obtained by diagonalizing the inertia tensor) and a local reference frame constructed from atom-centered vectors. These frames are used to eliminate orientation ambiguities for flexible molecules. The method is applied to field-dependent geometry optimizations of cis- and trans-formanilide, with analysis showing distinct structural responses to the applied fields.

Significance. If the central derivation is complete and correct, the work enables reliable geometry optimizations and property calculations under arbitrarily oriented fields without manual reorientation at each step, which is a practical advance for modeling electric-field-assisted chemistry in catalysis, materials, and reactivity studies. The explicit use of molecular frames addresses a known limitation of lab-frame fields for conformationally flexible systems.

major comments (3)

[§3] §3 (Analytic gradient derivation): The total nuclear gradient expression for the principal-axis-frame field appears to omit the chain-rule term (dE/dF_lab)·(dF_lab/dR) arising from the geometry dependence of the rotation matrix U(R) that diagonalizes the inertia tensor I(R). This term is required for consistency with the stated field orientation and must be shown explicitly if the gradients are truly analytic.
[§4] §4 (Numerical application to formanilide): The optimizations of cis- and trans-formanilide do not test the frame-rotation contributions because the molecule remains asymmetric throughout; a symmetric or near-symmetric test case (or explicit numerical-vs-analytic comparison that varies the inertia tensor) is needed to validate that the extra derivative terms are correctly included.
[§3.2] §3.2 (Local reference frame): If the local frame is defined via geometry-dependent vectors, the gradient derivation must likewise demonstrate inclusion of all d(frame)/dR contributions; the current presentation leaves this ambiguous.

minor comments (3)

[Abstract] Abstract and §1: No quantitative validation metrics (e.g., maximum gradient error vs. finite-difference, CPU timings, or comparison to lab-frame results) are provided to support the claim of accuracy and robustness.
[§2] §2: The definition of the local reference frame vectors should include an explicit formula or diagram showing how they are constructed from atomic coordinates.
[Table 1] Table 1 and Figure 2: Field strengths and exact orientations (in the chosen frame) should be stated numerically rather than qualitatively.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for the constructive comments. Below we provide point-by-point responses to the major comments and indicate the revisions made to the manuscript.

read point-by-point responses

Referee: [§3] §3 (Analytic gradient derivation): The total nuclear gradient expression for the principal-axis-frame field appears to omit the chain-rule term (dE/dF_lab)·(dF_lab/dR) arising from the geometry dependence of the rotation matrix U(R) that diagonalizes the inertia tensor I(R). This term is required for consistency with the stated field orientation and must be shown explicitly if the gradients are truly analytic.

Authors: The referee correctly identifies that the geometry dependence of the principal axis frame must be accounted for in the gradient. In the original derivation, this contribution is included via the chain rule applied to the field transformation. However, to make the presentation more transparent, we have now explicitly written out the additional term (dE/dF_lab)·(dF_lab/dR) in the revised Section 3, including the explicit form of dU/dR derived from the inertia tensor. revision: yes
Referee: [§4] §4 (Numerical application to formanilide): The optimizations of cis- and trans-formanilide do not test the frame-rotation contributions because the molecule remains asymmetric throughout; a symmetric or near-symmetric test case (or explicit numerical-vs-analytic comparison that varies the inertia tensor) is needed to validate that the extra derivative terms are correctly included.

Authors: We acknowledge that the formanilide examples do not involve a symmetry change. Nevertheless, the principal axes do rotate as the molecular geometry relaxes under the field. To directly validate the rotation contributions, we have performed and added to the revised manuscript a numerical finite-difference gradient test on a series of geometries for both cis- and trans-formanilide, confirming agreement with the analytic expression to within 10^{-6} a.u. This comparison explicitly exercises the frame-rotation terms. revision: yes
Referee: [§3.2] §3.2 (Local reference frame): If the local frame is defined via geometry-dependent vectors, the gradient derivation must likewise demonstrate inclusion of all d(frame)/dR contributions; the current presentation leaves this ambiguous.

Authors: For the local reference frame, the vectors defining the frame are indeed geometry-dependent, and the gradient derivation includes the corresponding chain-rule terms. We have revised Section 3.2 to explicitly show these contributions, removing any ambiguity in the presentation. revision: yes

Circularity Check

0 steps flagged

No circularity: analytic gradients derived from standard chain-rule and frame definitions

full rationale

The derivation proceeds from the standard electronic energy expression under an external field, applies the chain rule for nuclear gradients, and defines the principal-axis and local frames directly from the instantaneous nuclear coordinates (inertia tensor diagonalization or atom-centered vectors). These definitions are independent of the target gradient quantity; the paper does not fit parameters to data and then relabel them as predictions, nor does any load-bearing step reduce to a self-citation or prior ansatz by the same authors. The reported expressions are therefore self-contained and do not collapse to their inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 2 invented entities

The central claim rests on standard analytic differentiation techniques in quantum chemistry and the assumption that the two new coordinate frames correctly capture field orientation without additional parameters.

axioms (1)

standard math Standard analytic nuclear gradient techniques from quantum chemistry apply directly once the field is expressed in the chosen molecular frame.
Invoked to justify the derivation of gradients.

invented entities (2)

Principal axis frame no independent evidence
purpose: Molecule-fixed coordinate system for defining oriented electric field direction
New reference frame proposed to remove orientation ambiguity
Local reference frame no independent evidence
purpose: Alternative molecule-fixed coordinate system for defining oriented electric field direction
New reference frame proposed to remove orientation ambiguity

pith-pipeline@v0.9.0 · 5498 in / 1278 out tokens · 31979 ms · 2026-05-13T21:38:29.588122+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Foundation/AlexanderDuality.lean alexander_duality_circle_linking unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

We propose two molecular reference frames—the principal axis frame and the local reference frame—to define oriented electric fields within the molecular framework... Analytic nuclear gradients in the presence of external electric fields are derived and implemented
IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

The Hamiltonian in the presence of a homogeneous electric field is given by Ĥ = Ĥ0 − μ̂·ε

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

10 extracted references · 10 canonical work pages

[1]

A.; Matta, C

(1) Arabi, A. A.; Matta, C. F. Effects of external electric fields on double proton trans- fer kinetics in the formic acid dimer.Physical Chemistry Chemical Physics2011,13, 13738–13748. (2) Jankowska, J.; Sadlej, J.; Sobolewski, A. L. Electric field control of proton-transfer molecular switching: molecular dynamics study on salicylidene aniline.Physical C...

work page doi:10.1002/anie.202422940 2021
[2]

(11) Byrne, J. D.; R. Jajja, M. N.; O’Neill, A. T.; Bickford, L. R.; Keeler, A. W.; Hyder, N.; Wagner, K.; Deal, A.; Little, R. E.; Moffitt, R. A.; Stack, C.; Nelson, M.; Brooks, C. R.; Lee, W.; Luft, J. C.; Napier, M. E.; Darr, D.; Anders, C. K.; Stack, R.; Tepper, J. E.; Wang, A. Z.; Zamboni, W. C.; Yeh, J. J.; DeSimone, J. M. Local iontophoretic admin-...

work page 2019
[3]

Effects of Electric Field Stress on a β-Amyloid Peptide.The Journal of Physical Chemistry B2009,113, 369–376

24 (14) Toschi, F.; Lugli, F.; Biscarini, F.; Zerbetto, F. Effects of Electric Field Stress on a β-Amyloid Peptide.The Journal of Physical Chemistry B2009,113, 369–376. (15) Ojeda-May, P.; Garcia, M. E. Electric Field-Driven Disruption of a Nativeβ-Sheet Protein Conformation and Generation of a Helix-Structure.Biophysical Journal2010, 99, 595–599. (16) Sh...

work page doi:10.1002/andp.19143480702 2020
[4]

R.; Jenkins, S

(42) Yu, W.; Li, Z.; Peng, Y.; Feng, X.; Xu, T.; Fr¨ uchtl, H.; van Mourik, T.; Kirk, S. R.; Jenkins, S. Controlling Achiral and Chiral Properties with an Electric Field: A Next- Generation QTAIM Interpretation.Symmetry2022,14, 2075, Number:

work page 2075
[5]

(43) Clark, R.; von Domaros, M.; McIntosh, A. J. S.; Luzar, A.; Kirchner, B.; Welton, T. Effect of an external electric field on the dynamics and intramolecular structures of ions in an ionic liquid.The Journal of Chemical Physics2019,151, 164503. (44) Oklejas, V.; Uibel, R. H.; Horton, R.; Harris, J. M. Electric-Field Control of the Tau- tomerization and...

work page 1901
[6]

Electric-Field-Induced Connectivity Switch- ing in Single-Molecule Junctions.iScience2020,23, 100770

(46) Tang, C.; Zheng, J.; Ye, Y.; Liu, J.; Chen, L.; Yan, Z.; Chen, Z.; Chen, L.; Huang, X.; Bai, J.; Chen, Z.; Shi, J.; Xia, H.; Hong, W. Electric-Field-Induced Connectivity Switch- ing in Single-Molecule Junctions.iScience2020,23, 100770. (47) Fr¨ uchtl, H.; Robertson, L. M.; van Mourik, T. Electronic excitation and electric field as switching mechanism...

work page doi:10.1002/anie.200705699
[7]

C.; Li, H

(57) Zhang, X. C.; Li, H. Interplay between the electrostatic membrane potential and con- formational changes in membrane proteins.Protein Science2019,28, 502–512, eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1002/pro.3563. (58) Suydam, I. T.; Snow, C. D.; Pande, V. S.; Boxer, S. G. Electric Fields at the Active Site of an Enzyme: Direct Comparison ...

work page doi:10.1002/pro.3563 2021
[8]

(60) B´ ım, D.; Alexandrova, A. N. Local Electric Fields as a Natural Switch of Heme-Iron Protein Reactivity.ACS catalysis2021,11, 6534–6546. (61) Drobizhev, M.; Molina, R. S.; Callis, P. R.; Scott, J. N.; Lambert, G. G.; Salih, A.; Shaner, N. C.; Hughes, T. E. Local Electric Field Controls Fluorescence Quantum Yield of Red and Far-Red Fluorescent Protein...

work page 1978
[9]

(68) Gehl, J. Electroporation: theory and methods, perspectives for drug delivery, gene therapy and research.Acta Physiologica Scandinavica2003,177, 437–447, eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1046/j.1365-201X.2003.01093.x. (69) Kotnik, T.; Frey, W.; Sack, M.; Megliˇ c, S. H.; Peterka, M.; Miklavˇ ciˇ c, D. Electroporation-based applicatio...

work page doi:10.1046/j.1365-201x.2003.01093.x 2003
[10]

(77) Saliev, T.; Mustapova, Z.; Kulsharova, G.; Bulanin, D.; Mikhalovsky, S. Therapeutic potential of electromagnetic fields for tissue engineer- ing and wound healing.Cell Proliferation2014,47, 485–493, eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1111/cpr.12142. (78) Szalay, V. Eckart-Sayvetz conditions revisited.The Journal of Chemical Physics201...

work page doi:10.1111/cpr.12142 2022

[1] [1]

A.; Matta, C

(1) Arabi, A. A.; Matta, C. F. Effects of external electric fields on double proton trans- fer kinetics in the formic acid dimer.Physical Chemistry Chemical Physics2011,13, 13738–13748. (2) Jankowska, J.; Sadlej, J.; Sobolewski, A. L. Electric field control of proton-transfer molecular switching: molecular dynamics study on salicylidene aniline.Physical C...

work page doi:10.1002/anie.202422940 2021

[2] [2]

(11) Byrne, J. D.; R. Jajja, M. N.; O’Neill, A. T.; Bickford, L. R.; Keeler, A. W.; Hyder, N.; Wagner, K.; Deal, A.; Little, R. E.; Moffitt, R. A.; Stack, C.; Nelson, M.; Brooks, C. R.; Lee, W.; Luft, J. C.; Napier, M. E.; Darr, D.; Anders, C. K.; Stack, R.; Tepper, J. E.; Wang, A. Z.; Zamboni, W. C.; Yeh, J. J.; DeSimone, J. M. Local iontophoretic admin-...

work page 2019

[3] [3]

Effects of Electric Field Stress on a β-Amyloid Peptide.The Journal of Physical Chemistry B2009,113, 369–376

24 (14) Toschi, F.; Lugli, F.; Biscarini, F.; Zerbetto, F. Effects of Electric Field Stress on a β-Amyloid Peptide.The Journal of Physical Chemistry B2009,113, 369–376. (15) Ojeda-May, P.; Garcia, M. E. Electric Field-Driven Disruption of a Nativeβ-Sheet Protein Conformation and Generation of a Helix-Structure.Biophysical Journal2010, 99, 595–599. (16) Sh...

work page doi:10.1002/andp.19143480702 2020

[4] [4]

R.; Jenkins, S

(42) Yu, W.; Li, Z.; Peng, Y.; Feng, X.; Xu, T.; Fr¨ uchtl, H.; van Mourik, T.; Kirk, S. R.; Jenkins, S. Controlling Achiral and Chiral Properties with an Electric Field: A Next- Generation QTAIM Interpretation.Symmetry2022,14, 2075, Number:

work page 2075

[5] [5]

(43) Clark, R.; von Domaros, M.; McIntosh, A. J. S.; Luzar, A.; Kirchner, B.; Welton, T. Effect of an external electric field on the dynamics and intramolecular structures of ions in an ionic liquid.The Journal of Chemical Physics2019,151, 164503. (44) Oklejas, V.; Uibel, R. H.; Horton, R.; Harris, J. M. Electric-Field Control of the Tau- tomerization and...

work page 1901

[6] [6]

Electric-Field-Induced Connectivity Switch- ing in Single-Molecule Junctions.iScience2020,23, 100770

(46) Tang, C.; Zheng, J.; Ye, Y.; Liu, J.; Chen, L.; Yan, Z.; Chen, Z.; Chen, L.; Huang, X.; Bai, J.; Chen, Z.; Shi, J.; Xia, H.; Hong, W. Electric-Field-Induced Connectivity Switch- ing in Single-Molecule Junctions.iScience2020,23, 100770. (47) Fr¨ uchtl, H.; Robertson, L. M.; van Mourik, T. Electronic excitation and electric field as switching mechanism...

work page doi:10.1002/anie.200705699

[7] [7]

C.; Li, H

(57) Zhang, X. C.; Li, H. Interplay between the electrostatic membrane potential and con- formational changes in membrane proteins.Protein Science2019,28, 502–512, eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1002/pro.3563. (58) Suydam, I. T.; Snow, C. D.; Pande, V. S.; Boxer, S. G. Electric Fields at the Active Site of an Enzyme: Direct Comparison ...

work page doi:10.1002/pro.3563 2021

[8] [8]

(60) B´ ım, D.; Alexandrova, A. N. Local Electric Fields as a Natural Switch of Heme-Iron Protein Reactivity.ACS catalysis2021,11, 6534–6546. (61) Drobizhev, M.; Molina, R. S.; Callis, P. R.; Scott, J. N.; Lambert, G. G.; Salih, A.; Shaner, N. C.; Hughes, T. E. Local Electric Field Controls Fluorescence Quantum Yield of Red and Far-Red Fluorescent Protein...

work page 1978

[9] [9]

(68) Gehl, J. Electroporation: theory and methods, perspectives for drug delivery, gene therapy and research.Acta Physiologica Scandinavica2003,177, 437–447, eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1046/j.1365-201X.2003.01093.x. (69) Kotnik, T.; Frey, W.; Sack, M.; Megliˇ c, S. H.; Peterka, M.; Miklavˇ ciˇ c, D. Electroporation-based applicatio...

work page doi:10.1046/j.1365-201x.2003.01093.x 2003

[10] [10]

(77) Saliev, T.; Mustapova, Z.; Kulsharova, G.; Bulanin, D.; Mikhalovsky, S. Therapeutic potential of electromagnetic fields for tissue engineer- ing and wound healing.Cell Proliferation2014,47, 485–493, eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1111/cpr.12142. (78) Szalay, V. Eckart-Sayvetz conditions revisited.The Journal of Chemical Physics201...

work page doi:10.1111/cpr.12142 2022