Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DENSITY OF STATES”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Stochastically Realized Observables for Excitonic Molecular Aggregates

In this work, we show that a stochastic approach enables calculations of the optical properties of large 2-dimensional and nanotubular excitonic molecular aggregates. Previous studies of such systems relied on numerically diagonalizing the dense and disordered Frenkel Hamiltonian, which scales approximately as $\mathcal{O}$($N$ 3 ) for $N$ dye molecules. Our approach scales much more efficiently as $\mathcal{O}$($N$ log($N$)), enabling quick study of systems with a million of coupled molecules on the micrometer size scale. We calculate several important experimental observables, including the optical absorption spectrum and density of states, and develop a stochastic formalism for the participation ratio. Quantitative agreement with traditional matrix diagonalization methods is demonstrated for both small- and intermediate-size systems. The stochastic methodology enables the study of the effects of spatial-correlation in site energies on the optical signatures of large 2D aggregates. Our results demonstrate that stochastic methods present a path forward for screening structural parameters and validating experiments and theoretical predictions in large excitonic aggregates.

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

Toward Quantum Confinement in Graphitic Carbon Nitride-Based Polymeric Monolayers

raphitic carbon nitride (g-C 3 N 4 ) has garnered much attention due to its potential as an efficient metal-free photocatalyst. This study examines the evolution of properties in zero-dimensional quantum dots up to sizable clusters that mimic extended g-C 3 N 4 monolayers. We employ density functional theory to investigate systematically the structural, electronic, and optical properties of the g-C 3 N 4 -based melamine and heptazine building blocks using a “bottom-up” construction of polymeric monolayers. The results from our computations indicate that the melamine- and heptazine-based polymeric g-C 3 N 4 systems must be reduced to at least 2.74 and 4.00 nm, respectively, to observe an increase of its optical gap with a size reduction. The present study also examines the nature of the electronic transitions exhibited by g-C 3 N 4 -based monolayers through full natural transition orbital and density of state analyses. Furthermore, the most promising sites for water splitting and subsequent chemical doping studies are identified, which generally correspond to the nitrogen and carbon atoms, respectively.

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

Cation, Anion, and Radical Isomers of C4H4N: Computational Characterization and Implications for Astrophysical and Planetary Environments

Nitrogen-containing ions and molecules in the gas phase have been detected in non-Earth environments such as dark molecular clouds and more recently in the atmosphere of Saturn's moon Titan. These molecules may serve as precursors to larger heterocyclic structures that provide the foundation of complex biological molecules. On Titan, molecules of m/z 66 have been detected by the Cassini mission, and species of the empirical formula C 4 H 4 N may contribute to this signature. We have characterized seven isomers of C 4 H 4 N in anionic, neutral radical, and cationic states using density functional theory. Structures were optimized using the range-separated hybrid ωB97X-V with the cc-pVTZ and aug-cc-pVTZ basis sets. Anionic and radical C4H4N favor cyclic structures with aromatic and quasi-aromatic electron arrangements, respectively. Interestingly, ionization from the radical surface to the cation induces significant changes in structural stability, and the global minimum for positively charged isomers is CH 2 CCHCNH + , a pseudo-linear species reminiscent of cyanoallene. Select formation pathways to these structures from Titan's existing or postulated gas-phase species, reactions that are also relevant for other astrophysical environments, are discussed. By characterizing C 4 H 4 N isomers, we have identified energetically stable anionic, radical, and cationic structures that may be present in Titan's atmosphere and dark molecular clouds.

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

Anharmonic phonon dispersion in polyethylene

The second-order Green's function method for anharmonic crystals has been applied to an infinite, periodic chain of polyethylene taking into account up to quartic force constants. The frequency-independent approximation to the Dyson self-energy gives rise to numerous divergent resonances, which are fortuitous. Instead, solving the Dyson equation self-consistently with a frequency-dependent self-energy resists divergences from resonances or zero-frequency acoustic vibrations. Here, the calculated anharmonic phonon dispersion, which nonetheless displays many true resonances, and anharmonic phonon density of states furnish hitherto unknown details that explain smaller features of observed vibrational spectra.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Defect Energetics in Pseudo-Cubic Mixed Halide Lead Perovskites from First-Principles

Owing to the increasing popularity of lead-based hybrid perovskites for photovoltaic (PV) applications, it is crucial to understand their defect energetics and its influence on their optoelectronic properties. In this work, we simulate various point defects in pseudocubic structures of mixed iodide-bromide and bromide-chloride methylammonium lead perovskites with the general formula MAPbI 3-y Br y or MAPbBr 3-y Cl y (where y is between 0 and 3), and use first-principles based density functional theory computations to study their relative formation energies and charge transition levels. We identify vacancy defects and Pb on MA antisite defect as the lowest energy native defects in each perovskite. Here, we observe that while the low energy defects in all MAPbI 3-y Br y systems only create shallow transition levels, the Br or Cl vacancy defects in the Cl-containing pervoskites have low energy and form deep levels which become deeper for higher Cl content. We examine the structures and density of states of pure and defect-containing perovskite systems to obtain an understanding of the nature of defect levels. Further, we study extrinsic substitution by different elements at the Pb site in MAPbBr 3 , MAPbCl 3 , and the 50-50 mixed halide perovskite, MAPbBr 1.5 Cl 1.5 , and identify some transition metals that create lower energy defects than the dominant intrinsic defects and also create midgap charge transition levels.

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

Exploring Responses of Contact Kelvin Probe Force Microscopy in Triple-Cation Double-Halide Perovskites

Hybrid organic–inorganic perovskite (HOIP) solar cells have undergone unprecedentedly rapid improvement in power conversion efficiencies during the last decade, but controversy remains regarding the role of local phenomena and heterogeneity in grain by grain. These challenges in turn necessitate understanding the fundamental mechanisms operating at individual grains and exploring how these mechanisms give rise to the collective behaviors of the films. Although numerous efforts have targeted local phenomena in these materials using classical scanning probe microscopy modalities, the measurements remain difficult to interpret due to intrinsic coupling between possible ferroelectric, ion migration, bulk and surface electrochemical phenomena, and the effects of light on these behaviors. Here, we explore the electromechanical responses in (Cs 0.05 FA 0.85 MA 0.15 )Pb(I 0.85 Br 0.15 ) 3 perovskite thin film in the dark and under illumination using band excitation-piezoresponse force microscopy (BE-PFM). To identify and separate ferroic, electrostatic, and electrochemical responses, we utilize contact-Kelvin probe force microscopy (cKPFM). The cKPFM responses show weak variability with morphological features in the dark, while such variability significantly increases under illumination. Deconvolution of the cKPFM responses via principal component analysis (PCA) shows heterogeneity in grain to grain possibly due to variation in migration and segregation of halide ions in mixed perovskites at specific grains. Further, we propose an electrostatic mechanism based on the physics of semiconductors, which can directly link the observed cKPFM behaviors with the ionic and electronic density of states (DOS) and their time dynamics. Here, these findings suggest that cKPFM is a powerful probe of local ionic behaviors on the nanoscale, both in HOIPs and other materials.

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

Adsorbate Free Energies from DFT-Derived Translational Energy Landscapes

Adsorption free energies are fundamental to surface chemistry and catalysis. Standard models combine some assumed analytical form of the translational potential energy surface, often parameterized against density functional theory (DFT) calculations, with an analytical expression for the resultant translational densities of states (DOS), free energy, and entropy. Here we compare the performance of such models against numerical evaluations of the DOS and thermodynamic functions derived from solutions to the translational Schrodinger equation. We compare results for a translational potential energy surface (PES) derived from nudged eleastic band calculations with those obtained from adsorbate rastering across a series of monatomic (O, S, C, N, H) and polyatomic (NH x ) adsorbates on (100) Pt and Au facets. We find that analytical models as commonly parameterized have mixed performance for describing the translational PES and that the consequences for computed free energies are modest but potentially significant in microkinetic models. Here, numerical solutions are possible for modest to no additional computational cost over analytical models and thus should be considered when reliable free energy estimates are needed or translational potential energy surfaces are available.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

X-ray Transient Absorption Studies of Exciton Self-Trapping

Here, we present optical pump/X-ray absorption probe studies to probe changes in electronic distribution and local structure upon formation of a self-trapped exciton. The experiments were carried out on the mixed-valence halide-bridged transition metal linear chain charge density wave material [Pt(en 2 )][Pt(en 2 )Cl 2 ]·(ClO 4 ) 4 . Photoinduced changes in the Pt L III XANES spectrum were measured on a picosecond time scale following generation of excitons via optical excitation of the intervalence charge transfer transition. Ab initio FEFF9 modeling of the spectra and the associated angular momentum projected density of states was used to interpret the response. We find photoinduced changes reflecting electronic redistribution and local lattice distortions that correspond to strong localization of the self-trapped state.

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

Modulating the Electron Affinity of Small Bipyridyl Molecules on Single Gold Nanoparticles for Plasmon-Driven Electron Transfer

Developing controlled platforms for plasmon-driven chemistry is of great importance in catalytic reactions at the nanoscale. We report anion radical formation for five bipyridyl complexes of varying degrees of electron affinity utilizing optically fo-cused intraband (594 nm) and interband (532 nm) pump excitation of single gold nanoparticles. The surface-enhanced Raman scattering (SERS) of anion radicals for the five non-resonant adsorbed molecules 2,2’-bipyridine (22BPY), 4,4’-bipyridine (44BPY), trans-1,2-bis(4-pyridyl)ethylene (BPE), 1,2-bis(4-pyridyl)acetylene (BPA), and 1,2-bis(4-pyridyl)ethane (BPEt) were detected using localized surface-plasmon resonance (LSPR) excitation with 785 nm. The electron affinity of the five bipyridyl complexes were determined using electrochemistry. Molecules with low electron affinity experienced high-er instances of radical anion formation under a plasmon-coupled intraband electron transfer excitation (594 nm) whereas molecules with high electron affinity showed a preference for anion radical formation under direct interband electron trans-fer excitation (532 nm). The lowest unoccupied molecular orbital (LUMO) energy levels for low electron affinity surface-bound molecules (22BPY, BPEt) are on average ~0.43 eV higher for high electron affinity surface-bound molecules (BPA, BPE, 44BPY) as calculated using time-dependent density functional theory, elucidating the importance of plasmon coupling to energy levels that facilitate charge transfer pathways. We also show the ability to ‘activate’ high vs low electron affinity single nanoparticles with the choice of pump excitation wavelength. The findings show the complex interplay between molecular electron affinity, orbital overlap with the density of states of the plasmonic metal, and excitation energetics of the pump laser wavelength. Potential applications of this work include enhanced control over molecular scale catalysis, biosensor design, and solar energy capture.

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

Ab-Initio Investigation of Finite Size Effects in Rutile Titania Nanoparticles with Semilocal and Nonlocal Density Functionals

In this work, we employ hybrid and generalized gradient approximation (GGA) level density functional theory (DFT) calculations to investigate the convergence of surface properties and electronic gap of rutile titania nanoparticles with particle size. The surface energies and electronic gaps are calculated for cuboidal particles with minimum dimension ranging from 3.7 Angstrom (24 atoms) to 10.3 Angstrom (384 atoms) using a highly-parallel real-space DFT code to enable hybrid level DFT calculations of larger nanoparticles than are typically practical. We deconvolute the geometric and electronic finite size effects in surface energy, and evaluate the influence of defects on electronic gap and density of states (DOS). The electronic finite size effects in surface energy vanish when the minimum length scale of the nanoparticles becomes greater than 10 Angstrom. We show that this length scale is consistent with a computationally efficient numerical analysis of the characteristic length scale of electronic interactions. The surface energy of nanoparticles having minimum dimension beyond this characteristic length can be approximated using slab calculations that account for the geometric defects. In contrast, the finite size effects on the electronic gap and DOS is highly dependent on the shape and size of these particles. Furthermore, the DOS for cuboidal particles and more realistic particles constructed using the Wulff algorithm reveal that defect states within the electronic gap play a key role in determining the eigen value distribution of nanoparticles and the electronic gap does not converge to the bulk limit for the particle sizes investigated.

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

Four-Terminal Electrochemistry: A Back-Gate Controls the Electrochemical Potential of a 2D Working Electrode

We demonstrate that ultrathin semiconductor working electrodes integrated into metal–insulator–semiconductor (MIS) stacks are an enabling platform for understanding non-Faradaic semiconductor electrochemistry. Furthermore, 5 nm thick ZnO electrodes were deposited on 30 nm HfO 2 dielectric on a Pd “gate” electrode. Application of a bias V G between the Pd gate and the ZnO electrode causes electrons to accumulate in the ZnO layer as measured by recording the in-plane sheet conductance. By contacting the top surface of the ZnO layer with the electrolyte in a conventional three-electrode electrochemical cell, we show that the gate voltage V G modulates the electrochemical potential V ZnO of the ZnO film with respect to a reference electrode. Electrochemical potential changes ΔV ZnO up to –1 V vs Ag/Ag + are achieved for V G = +7 V. Furthermore, by measuring V ZnO vs V G , we extract the quantum capacitance CQ of the ZnO film as a function of the Fermi-level position, which provides a direct measure of the ZnO electronic density of states (DOS). Finally, we demonstrate that the gated ZnO working electrodes can disentangle the two principal components of electrochemical potential, namely, the Fermi-level shift Δδ and the double-layer charging energy eΔΦ EDL . This disentanglement hinges on a fundamental difference between back-gating and normal electrochemical control, namely, that electrochemical control requires double-layer charging, while back-gate control does not. Collectively, the results show that the backside gate electrode is an effective fourth terminal that enables measurements that are difficult to achieve in conventional three-terminal electrochemical setups.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Strain-Induced Shifts in Defective Graphite Phonon Modes Predicted by Density Functional Theory

Carbon fiber composites have gained attention as a structural material because of their high strength-to-weight ratio, and understanding the effect of defects on reactivity and mechanical properties is important for the longevity and safety of the composite. Although it is known that strain causes the underlying graphitic vibrational modes to redshift, it is not clear how strain may alter reactivity and defect-induced vibrational changes. To investigate the strain-induced phonon changes of defective carbon fiber composites, density functional theory calculations of graphite are used, including intercalated hydrogen and fluorine defects. By comparing changes in the bond lengths, formation energies, and phonon density of states for uniaxially and biaxially strained graphite, strain was found to generally make defect formation more favorable and the specific behavior changes are dependent on the strain direction and defect identity. Specifically, intercalated fluorine phonons are more sensitive to strain than hydrogen intercalation phonons, and strain applied along the zigzag direction alters the calculated properties more than strain along the armchair direction. Finally, these results highlight the importance of understanding the microstructural effect of deviations from the ideal material because small changes in strain or defect type can significantly alter the behavior of the carbon fiber composite core.

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

Fundamental Insights into Cathode Stability: Linking Compositional Tuning and Local Coordination in Complex Metal Oxides under Aqueous Transformations

Compositional tuning of complex metal oxides in Li-ion battery materials influences their performance as well as their end-of-life behavior, in particular, the tendency to release toxic metal cations in aqueous solution. We modeled ternary variants of a parent LiCoO 2 delafossite structure by varying the metal identity and relative amounts. This yielded ten model formulations of Li(A 4/6 B 1/6 C 1/6 )O 2 , where the material is enriched with the A metal and doped with B and C, with Ni, Mn, Co, Fe, Al, V, and Ti as constituent metals. To assess their stability in aqueous conditions, metal release energetics were calculated using a combination of Density Functional Theory calculations and thermodynamics. Metal release in ternary oxides is dictated by subtle variations in the coordination environment of the leaving group. To identify governing chemical features across diverse compositions with varying local coordination environments, we leverage random forest regression and descriptor importance analysis. A key result is that metal–oxygen orbital hybridization, quantified using a projected density-of-states-derived descriptor, H d/p , provides a physically grounded measure of interaction strength that governs metal release energetics. This refined perspective goes beyond conventional oxidation state considerations and offers more robust insights for materials science. Finally, we model defect surface-bound O 2 dimer formation as a proxy for reactive oxygen species (ROS) generation. The results show that Ni-rich compositions more readily stabilize spin-polarized O 2 dimers, corroborating experimental reports of an increased ROS-driven biological response. In conclusion, our results establish a compositional and electronic basis for metal release and surface oxygen reactivity that form a rationale for complex metal oxide design principles.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Fundamental Interactions of Bimetallic Cu x Pd y ( x + y = 4) Clusters Supported on the α-WC(0001) Surface and Their Performance for CO 2 Adsorption and Dissociation

The tungsten carbide α-WC(0001) surface, an active system for the activation of H 2 and important hydrogenation processes involving unsaturated hydrocarbons, can serve as a support of bimetallic clusters to produce materials with unique catalytic properties, opening routes for a wide range of technical applications. In particular, Cu x Pd y clusters are of particular interest because they combine metals with different properties. A stochastic method was applied to obtain the geometry of Cu x Pd y (x + y = 4) bare clusters, evaluating thousands of possibilities to obtain stable structures, yielding one isomer for Cu 4 , Cu 2 Pd 2 , Cu 1 Pd 3 , and Pd 4 and two isomers for Cu 3 Pd 1 . These clusters were supported on C and W terminations of the tungsten carbide (0001) surface, exploring all of the binding possibilities. The adsorption energies on the C and W terminations are in the ranges from −2.51 to −3.02 eV and from −2.26 to −3.30 eV, respectively. The strongest and weakest binding was seen for monometallic Cu 4 and Pd 4 clusters on both C and W terminations, while the Cu-Pd bimetallics have intermediate adsorption energies but lack a clear trend in terms of composition. The location of Cu x Pd y clusters over the (0001) surface induces a decrease in the work function relative to the pristine surface, while the cluster-surface Bader charge transfer and variations in the partial density of states point to changes in the electronic structure of the carbide atoms upon binding of the metallic clusters. The d-band center of the Cu x Pd y deposited on WC(0001) indicates an intermediate reactivity among Cu(111) and Pd(111) surfaces, modulating the reactivity with small numbers of Cu and Pd atoms, i.e., atom economy in catalyst design. The likelihood of existence of the most stable Cu x Pd y (x + y = 4) clusters in the temperature range of 298-400 K is 100%. The composite Cu x Pd y /α-WC(0001) (x + y = 4), is a nontrivial system since 22 isomers are needed to completely describe its structural properties. Among the isomers, seven structures are necessary to represent Cu 3 Pd 1 /α-WC(0001), five for Pd 4 /α-WC(0001), two for Cu 4 /α-WC(0001), and four for Cu 2 Pd 2 /α-WC(0001) and Cu 1 Pd 3 /α-WC(0001). The large number of cluster isomers supported on the tungsten carbide surface opens the door for several applications in the heterogeneous catalysis of the Cu x Pd y /α-WC(0001) composite, with the possibility of modulating the geometric, electronic, and chemical properties according to a desired application. Test studies for the adsorption of CO 2 indicate that the Cu x Pd y /α-WC(0001) composites are highly active for the adsorption and decomposition of the molecule, with bimetallic and admetal-carbide interactions playing a key role in the binding performance. In conclusion, this high activity indicates that these systems should be useful as catalysts for the conversion of CO 2 to oxygenates or light alkanes.

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

Electronic Structure of the Au Nanoparticle-TiO 2 Heterojunction: Influence of Nanoparticle Size, Shape, Oxygen Vacancies, and Temperature

The electronic structure of the gold nanoparticle-titanium dioxide (AuNP-TiO 2 ) heterojunction plays a critical role for charge transfer and recombination dynamics that underpin its photocatalytic function. However, building a representative model to capture the key physics remains a significant challenge. Here, we investigate the influence of the AuNP size and shape, as well as oxygen vacancy (V O ) defects at the anatase-phase TiO 2 (101) surface and the temperature of the heterojunction, on its interfacial electronic properties. Using density functional theory (DFT), we compare the closed-shell Au 20 and open-shell Au 19 clusters interfaced with pristine and V O defect TiO 2 surfaces. We find that the presence of a V O defect transforms pure TiO 2 from a p-type to an n-type semiconductor, reversing the interfacial band bending from downward to upward. For the heterosystem, density of states (DOS) analysis shows that V O minimally affects the open-shell Au 19 -TiO 2 system, but it significantly alters the closed-shell Au 20 -TiO 2 system, converting it from p-type to n-type at the Γ-point at 0 K. Furthermore, ab initio molecular dynamics (AIMD) simulations at 300 K reveal significant thermal fluctuations in AuNP positions relative to the TiO 2 surface. These fluctuations result in dynamic variations in the gap between the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) in all systems studied: Au 19 NP-Pristine TiO 2 , Au 20 NP-Pristine TiO 2 , Au 19 NP-V O defect TiO 2 , and Au 20 NP-V O defect TiO 2 . Complementary Bader charge analysis performed at both 0 K and finite-temperature AIMD snapshots supports the emergence of an upward band bending and the formation of a Schottky barrier at the V O containing heterojunctions. Notably, we find that at finite temperature, an Au atom can dynamically passivate the V O , leading to a pronounced widening of the HOMO-LUMO gap in the AuNP-V O defect TiO 2 heterojunctions. Our computational findings underscore the pivotal role of the V O defect and thermal effects in modulating interfacial band alignment, electronic states, and HOMO-LUMO gaps, providing insights for designing Au-TiO 2 heterojunctions with tailored electronic properties.

Mora Perez, Carlos [Lawrence Berkeley National Lab↗
DOE OSTI2020

A Metastable CaSH 3 Phase Composed of HS Honeycomb Sheets that is Superconducting Under Pressure

Evolutionary searches have predicted a number of ternary Ca–S–H phases that could be synthesized at pressures of 100–300 GPa. P6 3 /mmc CaSH 2 , Pnma CaSH 2 , Cmc2 1 CaSH 6 , and I$\bar{4}$CaSH 20 were composed of a Ca–S lattice along with H 2 molecules coordinated in a “side-on” fashion to Ca. The H–H bond lengths in these semiconducting phases were elongated because of H 2 σ → Ca d donation, and Ca d → H 2 σ* back-donation, via a Kubas-like mechanism. P$\bar{6}$m2 CaSH 3 , consisting of two-dimensional HS and CaH 2 sheets, was metastable and metallic above 128 GPa. The presence of van Hove singularities increased its density of states at the Fermi level and concomitantly the superconducting critical temperature, which was estimated to be as high as ~100 K at 128 GPa. Furthermore, this work will inspire the search for superconductivity in materials based upon honeycomb HX (X = S, Se, Te) and MH 2 (M = Mg, Ca, Sr, Ba) layers under pressure.

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

Spatially Resolved Investigation of Mixed Valence and Insulator-to-Metal Transition in an Organic Salt

Using scanning tunneling microscopy/spectroscopy (STM/STS), we investigate the evolution of electronic structures across the boundaries of 7,7,8,8-tetracyanoquinodimethane (TCNQ) and K-TCNQ assemblies on a weakly interacting substrate. Despite the semiconducting/insulating nature of TCNQ (TCNQ 0 ) and K-TCNQ (TCNQ -1 ), a continuum metallic-like density of states extending deep (~1.5 nm) into the TCNQ assembly is observed near domain boundary. We attribute the formation of these states to the abrupt change of molecular valence, which perturbs the electrostatics of the junction and creates local electric fields as evidenced by the band bending near domain boundary. To the best of our knowledge, this study provides the first microscopic understanding of the crucial physics occurring near domain boundaries of mixed valence in K-TCNQ, or broadly speaking charge transfer complexes, which highlights these boundaries as potential “weak” points to initiate the electric field induced insulator-to-metal transition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗