Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Electrostatic Potential”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Surface Reactivity Analysis of the Crude Oil–Brine–Limestone Interface for a Comprehensive Understanding of the Low-Salinity Waterflooding Mechanism

Low-salinity waterflooding (LSWF) has proven to improve oil recovery in carbonate formations through rock wettability alteration, although the underlying mechanism remains elusive. Multivalent ionic exchange and calcite dissolution have usually been investigated using geochemical analysis in secondary coreflooding. In this work, coreflooding, in tertiary mode, coupled with a surface reactivity analysis approach was employed to investigate the interplay of wettability alteration mechanisms such as mineral dissolution, electrostatic bond attraction, and the effect of pH at in situ conditions. Improved oil recovery (IOR) in tertiary mode observed by coreflooding in Indiana limestone rocks showed an ionic strength dependence, that is, reducing brine ionic strength resulted in an increase in oil recovery. Coreflooding results showed that the seawater and low-salinity brines deprived of Mg 2+ ions resulted in the lowest IOR in tertiary mode, indicating the significance of Mg 2+ on IOR in limestone rocks. Similar results were observed through the contact angle measurement showing the limestone rock wettability state dependence on ionic strength and the effect of Mg2+ ions. Surface reactivity analysis showed an increase in solution pH, Ca 2+ and Mg 2+ ions concentration in the effluent solution from the coreflooding in tertiary mode using low salinity brines (about 40 and 20% increase in the effluent composition for Ca 2+ and Mg 2+ , respectively). These changes in solution composition were used to calculate the in situ oil–brine and rock–brine zeta potential using a validated surface complexation model, showing the changes of zeta potential as brine is injected into limestone rocks. The results show that using seawater-like brine in tertiary mode resulted in no mineral dissolution or ionic exchange. However, improved oil recovery (IOR) using such seawater-like brine was due to wettability alteration caused by reduced electrostatic bond attraction associated with Mg 2+ ions [from 2.6 × 10 –13 (mol/m 2 ) 2 for formation water salinity to 1.5 × 10 –13 (mol/m 2 ) 2 for seawater salinity]. Using low-salinity brines in tertiary mode improved oil recovery by mineral dissolution, resulting in oil desorption and an increase in solution pH. Finally, the increase in solution pH also resulted in reduced electrostatic bond attraction which lead to rock wettability alteration using low-salinity brines.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2026

Pore-Scale Transport Effects in Electrochemical CO 2 Reduction on Gold via Coupled Microkinetic-Transport Modeling

A pore-resolved modeling framework is developed to quantify how pore-scale transport affects the intrinsic microkinetics of CO 2 -to-CO on Au. A DFT-informed microkinetic model is coupled self-consistently to a Generalized-Modified Poisson–Nernst–Planck (GMPNP) transport description in a single, electrolyte-filled cylindrical pore, allowing local concentrations and potential to feed back into site-specific reaction rates. FIB-SEM is used to determine pore sizes within realistic electrode materials. Across pore diameters, d p = 10–6000 nm, the surface-averaged CO 2 reduction rate is systematically reduced relative to the ideal microkinetic baseline where mass transport is not accounted for; the effectiveness factor 𝜂 𝑠,CO 2 , which quantifies this ratio, decreases rapidly at more negative potentials and is about 1% near −1.0 V vs SHE due to reactant depletion. Spatial maps reveal pore-bulk alkalization that emerges at higher cathodic bias, with a small, near-wall pH dip due to electrostatic repulsion of hydroxide at the cathode interface. For a fixed aspect ratio L p /d p , narrower pores exhibit larger 𝜂 𝑠,CO 2 by shortening diffusion paths, whereas variations in the aspect ratio L p /d p play a secondary role. A dimensionless analysis (surface/bulk Damköhler numbers) delineates operating regimes. In conclusion, this work offers a concept for incorporating microkinetic models into homogenized porous-electrode models through effectiveness factors and pore-size distribution.

Au-catalyst↗
DOE OSTI2025

Structure of molten ytterbium aluminum garnet

Rare earth aluminum garnets are important materials in optical, dielectric, and thermal barrier applications. To advance the understanding of their melt processing and glass forming ability, we report the atomic structure of molten Yb 3 Al 5 O 12 over 1770–2630 K, which spans the equilibrium and supercooled liquid regimes. The melt density at T m = 2283 K is 5.50 g cm –3 , measured via silhouette imaging of electrostatically levitated drops over 1010–2420 K. Four separate structure measurements were made with aerodynamically levitated melts using x-ray and neutron diffraction with isotope substitution of Yb ( 172 Yb, 174 Yb, or nat Yb). Empirical potential structure refinement models were developed, which are in excellent agreement with the experiments. Coordination environments for Al–O are predominantly 4- and 5-coordinate, with a mean coordination of n AlO = 4.43(8), while Yb–O environments mostly range from 5- to 8-coordinate, with n YbO = 6.26(8). The cation–oxygen polyhedra are connected primarily by corner-sharing, with edge-sharing constituting up to ~1/3 of the connectivity among polyhedra with Yb or higher-coordinated Al–O. Structurally, the –Al–O– network in molten Yb 3 Al 5 O 12 appears conducive to glass formation: n OAl = 1.85(3), there are 1.86 AlO x –AlO x connections per Al atom (e.g., a mixture of Q 3 and Q 4 units), and the modal ring size is six cations. These characterize a network that is somewhat less constrained compared to SiO 2 glass, yet Yb 3 Al 5 O 12 cannot be quenched into crystal-free glass. Here, aluminum garnet compositions with larger rare earth cations do form glass, so these characterizations help reveal the structural characteristics corresponding to the limit of glass forming ability in rare earth aluminates.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Transformer Analysis with ANSYS Maxwell and COMSOL

This is an extension of work described by Rodriguez et al . (2021). It continues analyses of a generic transformer design by Wes Greenwood. In this report, we summarize that work and add comparable results using the ANSYS Maxwell software (henceforward, “ANSYS”), and with COMSOL . We found the ANSYS and COMSOL calculations of inductance agreed well with previous results for simplified coils in air, and with a ferromagnetic core. We then describe the ANSYS and COMSOL approach and show results for a full transformer model based on magnetic field analyses. Finally, we present electrostatic analyses of E field enhancement, once again resolving individual wires. The purpose is to assess the electrostatic fields in order to locate where electric breakdown is likely to originate. We found the maximum enhancement between the secondary and either the primary or the tertiary at the end of the windings with a large potential difference.

42 ENGINEERING↗
DOE OSTI2025

A self-consistent electrostatic electrified shallow water model for charged liquid surface dynamics

The dynamics of an electrified liquid surface are investigated using a shallow water model that is self-consistently coupled with an electrostatic solver. To account for the spatiotemporal variation of a curved liquid surface, the electric field on curved surfaces is calculated by solving a two-dimensional electrostatic equation using the weighted least squares (WLSQ) interpolation method. First, the WLSQ implementation is verified with analytical theory obtained from a Laplace solution assuming a half-cylinder liquid surface profile. Then, the coupled electrified shallow water model is used to study the instability of liquid surface perturbation as a function of the potential drop between an electrode and conducting liquid, surface tension, and gravity. We present a linear dispersion theory of liquid surface instability for long-wavelength perturbations, including the effects of liquid viscosity. Furthermore, the effects of multiple sinusoidal surface perturbations on the electrified liquid surface instability are investigated.

Capillary waves↗
DOE OSTI2025

Reversible Ferroelectric Polarization Modulation of Chiral Molecular Ferroelectrics by Circularly Polarized Light

Abstract The optical modulation of ferroelectric polarization constitutes a transformative, non‐contact strategy for the precise manipulation of ferroelectric properties, heralding advancements in optically stimulated ferroelectric devices. Despite its potential, progress in this domain is constrained by material limitations and the intricate nature of the underlying mechanisms. Recent studies have achieved efficient regulation of ferroelectric polarization and thermal conductivity in chiral ferroelectric thin films through the application of left‐ and right‐handed circularly polarized light (LCP and RCP). Differential absorption of circularly polarized light (CPL) induces nonequilibrium carrier dynamics, generating distinctive interfacial electrostatic fields that enable precise control of ultrathin ferroelectric films. For (R)‐BINOL−DIPASi and (S)‐BINOL−DIPASi (C 26 H 26 O 2 Si), polarization changes surpass 23%, exhibiting opposite response under LCP and RCP excitation. In R chiral films, remnant polarization decreases from 1.05 µC cm − 2 under LCP to 0.85 µC cm − 2 under RCP, whereas in S chiral films, polarization increases from 0.85 µC cm − 2 under LCP to 0.98 µC cm − 2 under RCP. This reversible modulation facilitates reliable switching between ON and OFF states, presenting the potential of chiral ferroelectric materials for flexible, high‐speed integrated photonic sensor technologies.

Chemistry↗
DOE OSTI2025

Electron-scale Magnetic Holes Generation Driven by Whistler-to-Bernstein Mode Conversion in Fully Kinetic Plasma Turbulence

Magnetic holes (MHs) are coherent structures characterized by a strong and localized magnetic field amplitude dip, commonly observed in the heliosphere. These structures come in different sizes, from magnetohydrodynamic to kinetic scales. Subion-scale MHs are usually sustained by an electron current vortex and exhibit a strong electron temperature anisotropy, with higher temperatures perpendicular to the background magnetic field. Magnetospheric multiscale observations (MMSs) have revealed electron-scale MHs to be ubiquitous in the turbulent Earth’s magnetosheath and the solar wind, potentially playing an important role in the energy cascade and dissipation. Despite abundant observations, the origin of electron-scale MHs is still unclear and debated. In this work, we use fully kinetic simulations to investigate the role of plasma turbulence in generating electron-scale MHs. We find that the turbulence spontaneously produces electron-scale MHs via the following mechanism: first, large-scale turbulent velocity shears produce regions with high electron temperature anisotropy; these localized regions become unstable, generating oblique electron-scale whistler waves; as they propagate over the inhomogeneous turbulent background, whistler fluctuations develop an electrostatic component, turning into Bernstein-like modes; the strong electrostatic fluctuations produce current filaments that merge into an electron-scale current vortex; the resulting electron vortex locally reduces the magnetic field amplitude, finally evolving into an electron-scale MH. We show that MHs generated by this mechanism have properties consistent with MMSs and nontrivial kinetic features with a “mushroom”-shaped electron velocity distribution function. Our results have potential implications for understanding the formation and occurrence of electron-scale MHs in astrophysical turbulent and space environments, such as the Earth’s magnetosheath and the solar wind.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Impact of proton irradiation on 4H-SiC Low Gain Avalanche Detectors (LGADs)

Silicon carbide (SiC) particle detectors have the potential to provide time resolutions and robust performance in extreme environments which exceed that of silicon detectors. In this work 4H-SiC low gain avalanche detectors (LGADs) and complementary PiN diodes were irradiated with 2.5 GeV protons at fluences up to 3.33 × 10 14 p/cm 2 . The electrostatic performance of both irradiated and non-irradiated devices was evaluated using I-V and C-V measurements. Moreover, charge collection measurements using α particles were also conducted. SiC LGADs displayed a loss in rectification, high ON-state resistances >10 10 Ω-cm 2 , and the complete removal of gain when exposed to the highest proton fluence. Additionally, the reductions in capacitance and OFF-state current pointed to compensation of the gain layer as a gain reducing mechanism. The introduction of radiation-induced defects also hinders carrier acceleration, which reduces impact ionization and further gain reduction. LGADs exposed to a fluence of 1 × 10 14 p/cm 2 experienced a partial recovery in gain (originally having a value of ∼ 2) with a value of ∼ 1.65 when applying a bias of 500 V. Charge trapping was also seen in the reduction of charge collection efficiency (CCE) to ∼ 65% for the PiN diode. However, despite the reduction in device performance, the demonstration of a measurable signal after irradiation points to the potential of SiC LGAD detectors for future high energy physics applications.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Spatial Mapping of Electrostatic Fields in 2D Heterostructures

In situ electron microscopy is an effective tool for understanding the mechanisms driving novel phenomena in 2D structures. However, due to practical challenges, it is difficult to address these technologically relevant 2D heterostructures with electron microscopy. Here, we use the differential phase contrast (DPC) imaging technique to build a methodology for probing local electrostatic fields during electrical operation with nanoscale spatial resolution in such materials. We find that, by combining a traditional DPC setup with a high-pass filter, we can largely eliminate electric fluctuations emanating from short-range atomic potentials. Using a method based on this filtering algorithm, a priori electric field expectations can be directly compared with experimentally derived values to readily identify inhomogeneities and potentially problematic regions. Furthermore, we use this platform to analyze the electric field and charge density distribution across layers of hBN and MoS 2 .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Molecular simulation of liquid–vapor coexistence for NaCl: Full-charge vs scaled-charge interaction models

Scaled-charge models have been recently introduced for molecular simulations of electrolyte solutions and molten salts to attempt to implicitly represent polarizability. Although these models have been found to accurately predict electrolyte solution dynamic properties, they have not been tested for coexistence properties, such as the vapor pressure of the melt. In this work, we evaluate the vapor pressure of a scaled-charge sodium chloride (NaCl) force field and compare the results against experiments and a non-polarizable full-charge force field. The scaled-charge force field predicts a higher vapor pressure than found in experiments, due to its overprediction of the liquid-phase chemical potential. Reanalyzing the trajectories generated from the scaled-charge model with full charges improves the estimation of the liquid-phase chemical potential but not the vapor pressure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Influence of magnetic angle on the E × B drift in a magnetic presheath

In this work, we provide experimental measurements of the E X B flow above a non-magnetic, absorbing boundary in a magnetized plasma. Measurements are taken as a function of the angle between the background magnetic field and the boundary normal. The measurements are compared to computational predictions for oblique magnetic fields based on Chodura's model. Ion flow measurements were obtained with laser induced fluorescence and the presheath potential structure was measured with an emissive probe. The ions were found to accelerate to nearly 30% of the sound speed parallel to the boundary at oblique angles of the magnetic field (~80°) and the ion speed at the electrostatic sheath edge was found to decrease with angle. The edge of the magnetic presheath has been experimentally determined based on the appearance of the E X B drift and is found to be independent of the angle of the magnetic field. We also demonstrate that laser induced fluorescence is suitable as a non-perturbative diagnostic to measure the electric field in a magnetized sheath.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Binding Free Energy Analysis of Colicin D, E3 and E8 to Their Respective Cognate Immunity Proteins Using Computational Simulations

Colicins are antimicrobial proteins produced by bacteria for the purpose of destroying neighboring bacteria. Colicin activity is neutralized by a specific cognate immunity protein in order to protect the host. This study investigates the structural and binding mechanisms underlying the interaction of colicin-D, -E3 and -E8 to their respective immunity proteins (ImD, Im3 and Im8) using structure prediction, molecular dynamics (MD) simulations and MM-PBSA approach of free energy calculations. High-confidence colicin-immunity (Col-Im) complex structures predicted using AlphaFold2 were subjected to MD simulations of 150 ns with GROMACS and were analyzed for the binding free energy calculation using gmx_MMPBSA. Results showed that the complex of Col_E3-Im3 exhibited the most favorable binding free energy, driven by strong van der Waals and electrostatic interactions. Col_D-ImD and Col_E8-Im8 also showed the favorable binding. Electrostatics and hydrogen bonding emerged as a key factor driving binding and stability, while polar solvation acted as a destabilizing factor across all systems. These outcomes provide an understanding of the molecular mechanisms of Col-Im systems, with potential applications for developing natural antimicrobials for food safety.

Biochemistry & Molecular Biology↗
DOE OSTI2025

Multi-scale computational screening and mechanistic insights of cyclic amines as solvents for improved lignocellulosic biomass processing

A computational screening workflow for the efficient deconstruction of cellulose, lignin and hemicellulose fractions of lignocellulosic biomass using cyclic amines as solvents. Lignocellulosic biomass is a promising feedstock for production of affordable fuels and chemicals from renewable resources. Effective solubilization and subsequent deconstruction of its cellulose, hemicellulose, and lignin fractions is essential for the viability of future biorefineries. This study used quantum chemistry-based equilibrium thermodynamics methods to evaluate the potential of 650 cyclic amines to solubilize cellulose, hemicellulose, and lignin. The activity coefficients of solvent - biopolymer interactions were predicted using the COSMO-RS (COnductor-like Screening MOdel for Real Solvents) method and used to identify cyclic amines that can efficiently dissolve and extract selective fractions of biopolymers during biomass pretreatment. Among the 650 cyclic amines, 1-piperazineethanmaine was predicted to be an effective solvent for extracting all three polymers and was experimentally shown to achieve the highest lignin removal (97.1%). Non-covalent interaction, reduced density gradient and quantum chemical calculations were performed to elucidate the dissolution mechanism of lignin, cellulose and hemicellulose and gain further molecular level insights into the interactions between the cyclic amines and biomass polymers that promote efficient solubilization and extraction. These analyses indicated that 1-piperazineethanmaine and 1-methylimidazole make noncovalent van der Waals, electrostatic interactions and hydrogen bonding with lignin, leading to enhanced lignin removal, while the strong intramolecular hydrogen bonding interactions in cellulose and hemicellulose result in weaker solvent-biopolymer interactions. Overall, the computational approach provided an efficient method for identifying cyclic amines tailored for optimal biomass pretreatment and resulted in the identification of a potential new class of solvents for effective biomass pretreatment.

Kumar, Nikhil↗
DOE OSTI2020

Data-Driven Many-Body Models for Molecular Fluids: CO 2 /H 2 O Mixtures as a Case Study

In this study, we extend the scope of the many-body TTM-nrg and MB-nrg potential energy functions (PEFs), originally introduced for halide ion–water and alkali-metal ion–water interactions, to the modeling of carbon dioxide (CO2) and water (H2O) mixtures as prototypical examples of molecular fluids. Both TTM-nrg and MB-nrg PEFs are derived entirely from electronic structure data obtained at the coupled cluster level of theory and are, by construction, compatible with MB-pol, a many-body PEF that has been shown to accurately reproduce the properties of water. Although both TTM-nrg and MB-nrg PEFs adopt the same functional forms for describing permanent electrostatics, polarization, and dispersion, they differ in the representation of short-range contributions, with the TTM-nrg PEFs relying on conventional Born–Mayer expressions and the MB-nrg PEFs employing multidimensional permutationally invariant polynomials. By providing a physically correct description of many-body effects at both short and long ranges, the MB-nrg PEFs are shown to quantitatively represent the global potential energy surfaces of the CO 2 –CO 2 and CO 2 –H 2 O dimers and the energetics of small clusters, as well as to correctly reproduce various properties in both gas and liquid phases. Furthermore, building upon previous studies of aqueous systems, our analysis provides further evidence for the accuracy and efficiency of the MB-nrg framework in representing molecular interactions in fluid mixtures at different temperature and pressure conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Interface property–functionality interplay suppresses bimolecular recombination facilitating above 18% efficiency organic solar cells embracing simplistic fabrication

Efficient exciton-to-charge generation from the introduction of non-fullerene acceptors (NFAs) has been an important breakthrough in organic solar cell (OSC) developments. However, low device fill factors (FFs) following significant free charge carrier recombination loss continue to undermine their marketplace potential. Previous studies have successfully uncovered the importance of donor and acceptor domains in promoting charge transport. However, the functionality of donor/acceptor (D/A) interfaces relevant to free charge recombination remains unclear despite such interfaces being present throughout the material. In this work, the disorder-induced uphill bulk-to-D/A interface transport energy landscape is unveiled to enhance the polaron recombination resistance thereby also mitigating the triplet state formation from back charge transfer. In simple words, increasing the interface disorder while keeping the purer domains highly ordered will impart greater bulk-to-interface differential energy which then serves as a recombination energy barrier. Herein, these are made possible by varying the NFA outer side chains from linear alkyls to bulkier 2D phenylalkyls which influence the donor–acceptor interaction and define the interfacial disorder. Meanwhile, the extent of interface disorder without tradeoff in geminate losses is dependent on electrostatics and nanomorphology driving efficient exciton dissociation. By understanding these interplays, remarkable FFs over 80% and power conversion efficiencies (PCEs) above 18% are revealed to remain accessible even with simple binary component device fabrication without additional components/treatments thereby maintaining the scalability. Consequently, the principles uncovered here will serve as the foundation to reach the optimum potential of binary and simple systems, a prerequisite to ultimately realize the most cost-effective OSC design strategies for practical applications.

14 SOLAR ENERGY↗
DOE OSTI2026

Edge‐Driven Fringe‐Field Effects, Reduced Screening, and Bandgap Widening in Graphene Nanoribbons Enable Single‑Molecule Sensitivity

Graphene nanoribbons (GNRs) offer promising platforms for single‐molecule sensing due to their quasi‐1D channels and discrete electronic states, providing superior sensitivity toward molecular perturbations. While prior studies emphasize smoother edges as essential for optimal performance, the potential benefits of controlled edge roughness remain largely unexplored. Additionally, most investigations focus on GNR arrays, leaving critical edge‐ and width‐dependent factors, including fringe fields, bandgap widening, interactions between adsorbing molecules and GNR atoms, density of states (DOS) suppression, and electrostatic screening lengths, and their collective impact on sensitivity, poorly understood. Here, in this work, we fabricated field‐effect transistors using individual GNRs (widths: 200–20 nm) and characterized their response to molecular adsorption with perfluorooctanoic acid as the model analyte. Narrower ribbons displayed significantly enhanced sensitivity, yielding a coverage‐normalized response of 116 ± 10 mV per molecule in 20 nm‐wide GNRs (from calibrated ensemble Dirac‐point shifts). Experimental and theoretical analyses reveal that this heightened sensitivity arises from stronger fringe fields, width‐dependent quantum confinement effects, reduced DOS, and increased edge roughness that facilitates molecular anchoring, enhanced orbital overlap, and improved charge transfer efficiency. Our findings challenge the conventional assumption that smoother edges inherently enhance sensor performance, demonstrating that controlled edge disorder substantially boosts molecular sensitivity in GNR sensors.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Dynamics and Structure of Molten Neodymium Titanate

Binary rare earth titanates can be vitrified using levitation melt quenching, and their glasses are promising for optical applications. To better understand their thermophysical properties and dynamics, we have used molecular dynamics (MD) and empirical potential structure refinement (EPSR) to simulate molten neodymium titanate, 83TiO 2 –17Nd 2 O 3 . The models are compared with high-energy synchrotron X-ray diffraction of the melt structure over 1060–1960 K and thermophysical property measurements of electrostatically levitated droplets over 1500–2200 K. As temperature increases, the network periodicity increases and the correlation length decreases, both in monotonic trends that suggest continuous structural evolution through the melt and supercooled liquid regimes. At 1900 K, the mean Ti–O and Nd–O coordination numbers are 5.8(5) and 6.7(5), compared to the glass, 5.72(6) and 7.70(26). The equilibrium melting point is measured as 1860(20) K, over 100 K higher than previous estimates. Viscosity data from MD are in good agreement with experimental measurements. The dynamic propensity of Nd is similar to that of the network-forming Ti, suggesting that Nd does not participate as a conventional modifier. The viscosity of this fragile liquid is likely controlled by the breaking of cation–oxygen bonds and rearrangements of Ti–O polyhedra.

molecular dynamics↗
DOE OSTI2023

Anodic Dissolution Rates Accelerate with Decreasing MoS 2 Nanoflake Thickness

Electrochemical gating of 2D transition metal dichalcogenide (TMD) electrodes is an emerging frontier in the field of semiconductor electrochemistry. In this approach, an applied bias modifies the charge carrier concentration of the 2D TMD, causing band edge shifts and drastic changes in charge transfer rates. However, leveraging this effect for (photo)electrocatalysis is practically limited by the stability of the TMD material under gating conditions. Gerischer showed anodic dissolution of bulk TMD electrodes can occur in the dark and hypothesized that the reaction proceeds via an electron tunneling mechanism from surface states to the TMD conduction band [H. Gerischer, D. Ross, and M. Lubke, Z. Physickalische Chem., 139, 1 (1984)]. Here we investigate this possibility in single MoS 2 nanoflakes using in situ optical microscopy and explore whether Gerischer’s electron tunneling mechanism can explain anodic dissolution rates of thin 2D semiconductors. Here, spatially resolved measurements show anodic dissolution initiates at perimeter edge sites and accelerates exponentially with decreasing layer thickness, consistent with Gerischer’s tunneling mechanism. Interestingly, single layer MoS 2 is impervious to anodic dissolution at applied potentials >200 mV more positive than those required to drive dissolution in bulk and multilayer-thick nanoflakes.

2D semiconductor↗