pith. sign in

arxiv: 1907.07683 · v1 · pith:6AICS7WJnew · submitted 2019-07-17 · ❄️ cond-mat.mes-hall · physics.bio-ph

Optimal transport and colossal ionic mechano-conductance in graphene crown ethers

Subin Sahu , Justin Elenewski , Christoph Rohmann , Michael Zwolak This is my paper

Pith reviewed 2026-05-24 20:06 UTC · model grok-4.3

classification ❄️ cond-mat.mes-hall physics.bio-ph
keywords graphenecrown ethersion transportmechano-conductancenanofluidicsdiffusive transportstrain tuningelectromechanical tuning
0
0 comments X

The pith

Small strains in graphene crown ether pores produce 100% changes in ion conductance.

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

The paper establishes that crown ethers embedded in graphene form pores whose ion conductance can be tuned dramatically by mechanical strain. A 1% change in pore size produces a 100% shift in conductance because the strain alters local electrostatics and geometry, moving transport into a primarily diffusive regime that approaches barrierless flow rather than a knock-on process. This setup supplies a direct experimental handle on the balance of dehydration and electrostatic effects that biological channels achieve through mutation. A reader would care because the platform converts a mechanical input into an electrical output while remaining simple enough to isolate the physical conditions for optimal transport rates and selectivity.

Core claim

Graphene crown ether pores display colossal ionic mechano-conductance in which 1% pore strain produces a 100% change in conductance. The transport remains primarily diffusive and electromechanically tunable, trending toward barrierless conditions instead of knock-on mechanisms. Mechanical modulation of the current therefore reports directly on the local electrostatic environment inside the pore.

What carries the argument

Strain-induced modification of pore geometry and local electrostatics that shifts the system toward barrierless diffusive ion transport

If this is right

  • Mechanical current modulation measurements directly reveal the local electrostatic conditions at the pore.
  • The platform enables electromechanically tunable nanofluidic devices with high transport rates.
  • Optimal transport occurs in a diffusive regime rather than a knock-on regime.
  • The observations supply a simplified model for the physical conditions underlying high-rate, selective ion flow in biological channels.

Where Pith is reading between the lines

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

  • Similar strain-tuning could be tested in other atomically thin membranes to create mechanically gated ion channels.
  • The effect suggests a route to sensors that transduce local mechanical strain into ionic current signals without external electrodes.
  • Device arrays with spatially varying strain might allow simultaneous measurement of multiple electrostatic environments in one experiment.

Load-bearing premise

Conductance changes arise only from strain-altered electrostatics and geometry, with no significant role for defects, edge effects, or shifts in ion dehydration barriers.

What would settle it

An experiment that strains the graphene sheet while holding pore size and dehydration fixed and finds no conductance change would falsify the claimed link between strain and electrostatic tuning.

Figures

Figures reproduced from arXiv: 1907.07683 by Christoph Rohmann, Justin Elenewski, Michael Zwolak, Subin Sahu.

Figure 1
Figure 1. Figure 1: Potential transport mechanisms in graphene crown ether pores. a. Using partial charges qO = −0.24 e, con￾sistent with the electrostatic potential from DFT, the ion transport mechanism is drift–diffusion. In this case, a K+ (purple sphere) finds an empty pore and translocates through it; the pore then remains empty again for several nanoseconds (see the SM). b. At larger partial charge (qO = −0.54 e or qO =… view at source ↗
Figure 2
Figure 2. Figure 2: Colossal ionic mechano–conductance. a. Schematic of graphene on (or embedded within) a polymeric matrix support (such as ∼ 50 nm thick epoxy resin) with a window where the crown ether pore is located. The membrane can be strained by stretching or bending (e.g., via a piezoelectric actuator offset from the window). Alternative experimental setups are possible, such as metallic regions that bind the graphene… view at source ↗
Figure 3
Figure 3. Figure 3: Free energy landscape. Free energy profile, ∆FK, for a K+ translocating along the z-axis of an 18- crown-6 graphene pore at different strains. The charge on the oxygen atom, qO, is presented at the top of each plot. The peaks and valleys in ∆FK(z) are due to the balance between dehydration energy penalties and electrostatic inter￾actions with the charged pore atoms. For qO = −0.54, there is an additional c… view at source ↗
Figure 4
Figure 4. Figure 4: Optimum ion transport and selectivity. Ionic current versus strain in a pore with qO = −0.24 e and an applied bias of 0.25 V. For small strain (blue line), the current increases rapidly, commensurate with a decrease in the outer barriers with increasing strain. This gives an overall flatten￾ing of the free-energy profile. Transport (and selectivity, see below) is optimal near 3 % strain (green line) when t… view at source ↗
read the original abstract

Biological ion channels balance electrostatic and dehydration effects to yield large ion selectivities alongside high transport rates. These macromolecular systems are often interrogated through point mutations of their pore domain, limiting the scope of mechanistic studies. In contrast, we demonstrate that graphene crown ether pores afford a simple platform to directly investigate optimal ion transport conditions, i.e., maximum current densities and selectivity. Crown ethers are known for selective ion adsorption. When embedded in graphene, however, transport rates lie below the drift-diffusion limit. We show that small pore strains -- 1 % -- give rise to a colossal -- 100 % -- change in conductance. This process is electromechanically tunable, with optimal transport in a primarily diffusive regime, tending toward barrierless transport, as opposed to a knock-on mechanism. Measurements of mechanical current modulation will yield direct information on the local electrostatic conditions of the pore. These observations suggest a novel setup for nanofluidic devices while giving insight into the physical foundation of evolutionarily--optimized ion transport in biological pores.

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.

Referee Report

0 major / 3 minor

Summary. The manuscript demonstrates that graphene-embedded crown ether pores provide a tunable platform for ion transport studies. It reports that 1% mechanical strain on the pore produces a 100% change in ionic conductance, operating primarily in a diffusive regime that approaches barrierless transport (distinct from knock-on mechanisms), with electromechanical tunability arising from strain-modified electrostatics and geometry. This is positioned as yielding insights into biological ion channel optimization and novel nanofluidic device concepts.

Significance. If validated, the colossal mechano-conductance sensitivity at small strains offers a direct experimental handle on local pore electrostatics via current modulation measurements, complementing mutation-based studies in biology. The diffusive-regime assignment and contrast to knock-on transport provide mechanistic clarity, while the platform's simplicity could enable systematic exploration of optimal transport conditions.

minor comments (3)
  1. [Abstract and Results] The abstract and main text state the 100% conductance change for 1% strain without accompanying error bars, convergence tests, or explicit criteria for data inclusion/exclusion from the simulations; adding these would strengthen the quantitative claim.
  2. [Discussion] Clarify the precise definition of the 'primarily diffusive regime' (e.g., via explicit comparison of drift vs. diffusion contributions or barrier heights) and how it is distinguished from residual knock-on character in the computed trajectories.
  3. [Methods] The manuscript should include a brief statement on the force-field parameters or DFT settings used to compute the strain-dependent electrostatic potential and ion-pore interactions to support reproducibility.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for their positive assessment of the manuscript, the recognition of its significance for ion transport studies, and the recommendation for minor revision. No specific major comments were provided in the report.

Circularity Check

0 steps flagged

No significant circularity in derivation chain

full rationale

The paper's central claims (colossal conductance modulation via 1% strain in a diffusive regime) are presented as outcomes of modeling strain-modified electrostatics and pore geometry, benchmarked against drift-diffusion limits. No equations or steps in the provided text reduce by construction to fitted inputs, self-definitions, or self-citation chains. The derivation remains self-contained with independent content from the simulations and comparisons to biological channels.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Only the abstract is available, so free parameters, axioms, and invented entities cannot be enumerated from the text. The central claim rests on unstated simulation or theoretical assumptions about pore electrostatics and strain response.

pith-pipeline@v0.9.0 · 5714 in / 1200 out tokens · 17806 ms · 2026-05-24T20:06:55.971269+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

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/Cost/FunctionalEquation.lean washburn_uniqueness_aczel (J uniqueness) echoes
    ?
    echoes

    ECHOES: this paper passage has the same mathematical shape or conceptual pattern as the Recognition theorem, but is not a direct formal dependency.

    The free energy of a potassium ion ... ΔF_K = ΔE_deh + E_KO + E_KC + E_KK ... I ≈ A e^{-ΔF*/kBT} ... ΔI/I ≈ e^{α a s χ / kBT} − 1 ... tending toward barrierless transport

  • IndisputableMonolith/Cost Jcost_pos_of_ne_one; Jcost_unit0 echoes
    ?
    echoes

    ECHOES: this paper passage has the same mathematical shape or conceptual pattern as the Recognition theorem, but is not a direct formal dependency.

    colossal – 100 % – change in conductance ... small pore strains – 1 % ... optimal transport in a primarily diffusive regime

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

34 extracted references · 34 canonical work pages

  1. [1]

    Perdew, J. P., K. Burke, and M. Ernzerhof, 1996. Generalized Gradient Approximation Made Simple.Phys. Rev. Lett. 77:3865–3868

    work page 1996

  2. [2]

    Semiempirical GGA-type density functional constructed with a long-range dispersion correction

    Grimme Stefan, 2006. Semiempirical GGA-type density functional constructed with a long-range dispersion correction. J. Comput. Chem. 27:1787–1799

    work page 2006

  3. [3]

    Joubert, 1999

    Kresse, G., and D. Joubert, 1999. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. Rev. B 59:1758–1775

    work page 1999

  4. [4]

    Furthm ¨uller, 1996

    Kresse, G., and J. Furthm ¨uller, 1996. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 54:11169–11186

    work page 1996

  5. [5]

    J., and J

    Monkhorst, H. J., and J. D. Pack, 1976. Special points for Brillouin-zone integrations. Phys. Rev. B 13:5188–5192

    work page 1976

  6. [6]

    Souza Filho, E

    Meunier, V ., A. Souza Filho, E. Barros, and M. Dresselhaus, 2016. Physical properties of low-dimensional s p 2-based carbon nanostructures. Rev. Mod. Phys. 88:025005

    work page 2016

  7. [7]

    Sanville, and G

    Tang, W., E. Sanville, and G. Henkelman, 2009. A grid-based Bader analysis algorithm without lattice bias.J. Phys. Condens. Matter 21:084204

    work page 2009

  8. [8]

    E., and A

    Feller, S. E., and A. D. MacKerell, 2000. An Improved Empirical Potential Energy Function for Molecular Simulations of Phospholipids. J. Phys. Chem. B 104:7510–7515

    work page 2000

  9. [9]

    C., and J

    Hariharan, P. C., and J. A. Pople, 1973. The influence of polarization functions on molecular orbital hydrogenation energies. Theoret. Chimica Acta 28:213

    work page 1973

  10. [10]

    Yang, and R

    Lee, C., W. Yang, and R. G. Parr, 1988. Development of the Colle–Salvetti correlation–energy formula into a functional of the electron density. Phys. Rev. B 37:785

    work page 1988

  11. [11]

    D., 1993

    Becke, A. D., 1993. Density–functional thermochemistry. III. The role of exact exchange. J. Chem. Phys. 98

    work page 1993

  12. [12]

    Stephens, P. J., F. J. Devlin, C. F. Chabalowski, and M. Frisch, 1994. J. Phys. Chem. 98:11623

    work page 1994

  13. [13]

    Frisch, M. J., G. W. Trucks, H. B. Schlegel, G. E. Scuseria, M. A. Robb, J. R. Cheeseman, G. Scalmani, V . Barone, B. Mennucci, G. A. Petersson, H. Nakatsuji, M. Caricato, X. Li, H. P. Hratchian, A. F. Izmaylov, J. Bloino, G. Zheng, J. L. Sonnenberg, M. Hada, M. Ehara, K. Toyota, R. Fukuda, J. Hasegawa, M. Ishida, T. Nakajima, Y . Honda, O. Kitao, H. Naka...

    work page 2009

  14. [14]

    M., and K

    Breneman, C. M., and K. B. Wiberg, 1990. Determining atom-centered monopoles from molecular electrostatic potentials. The need for high sampling density in formamide conformational analysis. J. Comput. Chem. 11:361

    work page 1990

  15. [15]

    H., 1989

    Dunning Jr, T. H., 1989. Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen. J. Chem. Phys. 90:1007–1023

    work page 1989

  16. [16]

    Kendall, R. A., T. H. Dunning Jr, and R. J. Harrison, 1992. Electron affinities of the first-row atoms revisited. Systematic basis sets and wave functions. J. Chem. Phys. 96:6796–6806

    work page 1992

  17. [17]

    Iannuzzi, F

    Hutter, J., M. Iannuzzi, F. Schiffmann, and J. VandeV ondele, 2014.Wiley Inderdiscip. Rev: Comput. Mol. Sci. 4:15

    work page 2014

  18. [18]

    Antony, S

    Grimme, S., J. Antony, S. Ehrlich, and S. Kreig, 2010. A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-d) for the 94 elements H-Pu. J. Chem. Phys. 132:154104

    work page 2010

  19. [19]

    Krack, F

    VandeV ondele, J., M. Krack, F. Mohamed, M. Parrinello, T. Chassaing, and J. Hutter, 2005. Quickstep: fast and accurate density functional calculations using a mixed Gaussian and plane waves approach,. Comp. Phys. Comm. 167:103

    work page 2005

  20. [20]

    Godbout, N., D. R. Salahub, J. Andzelm, and E. Wimmer, 1992. Can. J. Chem. 70:560

    work page 1992

  21. [21]

    Pseudopotentials for H to Kr optimized for gradient–corrected exchange–correlation functionals

    Krack, M., 2005. Pseudopotentials for H to Kr optimized for gradient–corrected exchange–correlation functionals. Theo. Chem. Account 114:145

    work page 2005

  22. [22]

    Cieplak, W

    Bayly, C., P. Cieplak, W. D. Cornell, and P. A. Kollman., 1993. A Well–Behaved Electrostatic Potential Based Method Using Charge Restraints For Determining Atom–Centered Charges: The RESP Model. J. Phys. Chem. 97:10269

    work page 1993

  23. [23]

    Molecular dynamics simulation of thermal ripples in graphene with bond-order-informed harmonic constraints

    Smolyanitsky, A., 2014. Molecular dynamics simulation of thermal ripples in graphene with bond-order-informed harmonic constraints. Nanotechnology 25:485701. 9

    work page 2014

  24. [24]

    Jorgensen, W. L., J. Chandrasekhar, J. D. Madura, R. W. Impey, and M. L. Klein, 1983. Comparison of simple potential functions for simulating liquid water. J. Chem. Phys. 79:926–935

    work page 1983

  25. [25]

    Roux, 1994

    Beglov, D., and B. Roux, 1994. Finite Representation of an Infinite Bulk System: Solvent Boundary Potential for Computer Simulations. J. Chem. Phys. 100:9050

    work page 1994

  26. [26]

    Zwolak, 2018

    Sahu, S., and M. Zwolak, 2018. Maxwell-Hall access resistance in graphene nanopores. Phys. Chem. Chem. Phys. 20:4646 – 4651

    work page 2018

  27. [27]

    Zwolak, 2018

    Sahu, S., and M. Zwolak, 2018. Golden aspect ratio for ion transport simulation in nanopores. Phys. Rev. E 98:012404

    work page 2018

  28. [28]

    Dalke, and K

    Humphrey, W., A. Dalke, and K. Schulten, 1996. VMD: Visual molecular dynamics. J. Mol. Graphics 14:33–38

    work page 1996

  29. [29]

    Phillips, J. C., R. Braun, W. Wang, J. Gumbart, E. Tajkhorshid, E. Villa, C. Chipot, R. D. Skeel, L. Kale, and K. Schulten, 2005. Scalable molecular dynamics with NAMD. J. Comput. Chem. 26:1781–1802

    work page 2005

  30. [30]

    York, and L

    Darden, T., D. York, and L. Pedersen, 1993. Particle mesh Ewald: An N log (N) method for Ewald sums in large systems. J. Chem. Phys. 98:10089–10092

    work page 1993

  31. [31]

    Martyna, G. J., D. J. Tobias, and M. L. Klein, 1994. Constant pressure molecular dynamics algorithms. J. Chem. Phys. 101:4177–4189

    work page 1994

  32. [32]

    Chipot, 2004

    H ´enin, J., and C. Chipot, 2004. Overcoming free energy barriers using unconstrained molecular dynamics simulations. J. Chem. Phys. 121:2904–2914

    work page 2004

  33. [33]

    Di Ventra, and M

    Sahu, S., M. Di Ventra, and M. Zwolak, 2017. Dehydration as a Universal Mechanism for Ion Selectivity in Graphene and Other Atomically Thin Pores. Nano Lett. 17:4719–4724

    work page 2017

  34. [34]

    H., 2002

    Nelson, P. H., 2002. A permeation theory for single-file ion channels: Corresponding occupancy states produce Michaelis–Menten behavior. J. Chem. Phys. 117:11396–11403

    work page 2002