Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular 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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Strong stretching theory for pH-responsive polyelectrolyte brushes in large salt concentrations

Here, we develop a theory for describing the thermodynamics, configuration, and electrostatics of strongly-stretched, pH-responsive polyelectrolyte (PE) brushes in the presence of large salt concentrations. The aim of the paper, therefore, is to study the properties of a PE brush in a salt concentration regime (namely, large concentrations of several molars) that has been hitherto unexplored theoretically in the context of PE brushes but can be routinely encountered in molecular scale simulations of the problem. The brushes are modelled using our recently developed augmented Strong Stretching Theory (SST), while the effect of the presence of the large salt concentration is accounted for by including the contributions of three different types of non-Poisson–Boltzmann (non-PB) effects to the free energy description of the PE brush induced electric double layer (EDL). These non-PB effects are ionic non-mean-field ion–ion correlations, solvent polarization, and the finite size effect of the ions and water dipoles. We study the individual influences of these different effects and show that the ion–ion correlations and solvent polarization effect reduce the brush height which consequently enhances the monomer density and leads to an electrostatic potential distribution of the brush induced EDL that has a larger magnitude at near-wall locations and becomes zero at shorter distances from the wall. The finite size effect, on the other hand, increases the brush height and therefore, weakens the monomer density and leads to a smaller near-wall magnitude of the EDL potential that becomes zero at larger distances from the wall. In time, we consider the impact of all the three non-PB effects simultaneously and show that the ion–ion correlations and solvent polarization effect dominate the size effects and dictate the overall brush configuration and the EDL electrostatics. We also point out that the influence of all the three non-PB effects becomes the largest for a larger salt concentration and a smaller bulk pH. Finally, we compare our theoretical predictions with those obtained from our recently developed all-atom MD simulation model and obtain an excellent match.

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

Charge influence on point defect interactions with xenon bubbles in uranium dioxide

The interaction of xenon (Xe) bubbles with small self-defects in uranium dioxide (UO2) has been studied using molecular statics simulations. The results show that the pressure and size of bubbles have a minimal impact on the heterogeneity of the interaction between a charge neutral Xe bubble and charged point defects, while the local charge distribution around the bubble overwhelmingly determines the interaction strength and critical interaction distance. The charge effect on the interactions between Xe bubbles and defects was further confirmed by assessing the point defect interaction energies with negatively or positively charged Xe bubbles. The Xe density, or pressure, has a much smaller effect on the interaction energy between charged bubbles and small defects at least for low and medium pressure bubbles. The interactions between charged Xe bubbles and point defects basically follow the Coulomb electrostatic interaction law, which is independent of the empirical potentials employed for the Xe–UO2 system in this work.

Yang, L. (ORCID:0000000322166071)↗
DOE OSTI2025

Probing specific ion effects at air-aqueous dibutyl phosphate interfaces using vibrational sum frequency generation spectroscopy

Molecular properties at air–liquid and liquid–liquid interface hold the key to many processes involving molecular transport across phase boundaries from aerosol formation to carbon cycling and material separation using solvent extraction techniques. Using dibutyl phosphate (DBP) as a representative for partially aqueous soluble surfactants, the specific ion effect (SIE) of the Hofmeister series cations Cs + , Na + , Li + , and Mg 2+ on the partition and interaction between surfactant molecules and water molecules in the air–aqueous interface are investigated using vibrational sum frequency generation spectroscopy and surface tension measurements. In the presence of 1 mM and 1M bulk aqueous phase ionic strength salt concentrations, fundamental qualitative relationships are observed for the salting out of DBP relative to bulk aqueous phase nitrate salt concentrations and the specific cations species. At 1 mM ionic strength, the interfacial charge and hence the interfacial potential modulates the electrostatic interactions; in particular, the counter cations partially screen the negatively charged interface induced by the DBP in a direct Hofmeister order. At 1M ionic strength, the electric field at the interface or interfacial potential is effectively neutralized, and the counter cations promote the partitioning of DBP to the interface depending on their specific interaction with the DBP head group and metal ion hydration properties. The present results lay a foundation to study SIEs of heavier metals on hydrophobic-aqueous DBP interfaces.

Louie, Christina [Washington State University, Pul↗
DOE OSTI2023

Polarity‐Driven Atomic Displacements at the 2D Mg 2 TiO 4 ‐MgO (001) Oxide Interface for Hosting Potential Interlayer Excitons

Abstract Interlayer excitons in solid‐state systems have emerged as candidates for realizing novel platforms ranging from excitonic transistors and optical qubits to exciton condensates. Interlayer excitons have been discovered in 2D transition metal dichalcogenides, with large exciton binding energies and the ability to form various van der Waals heterostructures. Here, an oxide system consisting of a single unit cell of Mg 2 TiO 4 on MgO (001) is proposed as a platform for hosting interlayer excitons. Using a combination of density functional theory (DFT) calculations, molecular beam epitaxy growth, and in situ crystal truncation rod measurements, it is shown that the Mg 2 TiO 4 ‐MgO interface can be precisely controlled to yield an internal electric field suitable for hosting interlayer excitons. The atoms in the polar Mg 2 TiO 4 layers are observed to be displaced to reduce polarity at the interface with the non‐polar MgO (001) surface. Such polarity‐driven atomic displacements strongly affect electrostatics of the film and the interface, resulting in localization of filled and empty band‐edge states in different layers of the Mg 2 TiO 4 film. The DFT calculations suggest that the electronic structure is favorable for localization of photoexcited electrons in the bottom layer and holes in the top layer, which may bind to form interlayer exciton states.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Evidence of Au(111) topological states in a kagome analogue lattice and their robustness beyond ultra-low temperatures and defect-free conditions

The experimental realization of kagome lattices that exhibit the predicted coexistence of topological states with high electron kinetics and non-dispersing quantum states remains challenging. Additionally, the robustness of these states against structural perturbations has rarely been explored. Here, we report on the formation of an analogue kagome structure via the electrostatic self-assembly of 4,7-dibromobenzo[c]-1,2,5-thiadiazole (2Br-BTD) molecules on Au(111). Local spectroscopic measurements, supported by theoretical calculations, reveal that the weak molecular coupling reshapes the topological-induced Shockley surface state of Au(111) by imposing a (7 × 7) periodicity, resulting in new band crossings. The molecular overlayer favours the opening of electron gaps at these positions manifested as sharp peaks in dI/dV spectra and electron localization in either the hexagonal or triangular sublattices of the kagome structure. To explore the robustness of these topological states, we monitored their stability under varied conditions, including different temperatures, the unaltered herringbone reconstruction of the Au(111) surface and local structural relaxation of the molecular assembly. These results demonstrate the degree of topological protection of these states, which holds potential for fundamental and applied research.

Austin, Dave [University of Central Florida, Orlan↗
NASA NTRS1999

Molecular Modeling and Experimental Investigations of Nonlinear Optical Compounds Monosubstituted Derivatives of Dicyanovinylbenzene

A search for potential nonlinear optical compounds was performed using the Cambridge Structure Database and molecular modeling. We investigated a series of monosubstituted derivatives of dicyanovinylbenzene, since the nonlinear optical (NLO) properties of such derivatives (o-methoxy-dicyanovinylbenzene, DIVA) were studied earlier. The molecular geometry of these compounds was investigated with x-ray analysis and discussed along with the results of molecular mechanics and ab initio quantum chemical calculations. The influence of crystal packing on the planarity of the molecules of this series has been revealed. Two new compounds from the series studied, ortho-F and para-Cl-dicyanovinylbenzene, according to powder measurements, were found to be NLO compounds in the crystal state about 10 times more active than urea. The peculiarities of crystal structure formation in the framework of balance between van der Waals and electrostatic interactions have been discussed. The crystal shape of DIVA and two new NLO compounds have been calculated on the basis of the known crystal structure.

Timofeeva, Tatiana V.↗
NASA NTRS1991

Ab initio study of the electrostatic multipole nature of torsional potentials in CH3SSCH3, CH3SSH, and HOOH

The origin of torsional potentials in H3CSSCH3, H3CSSH, and HOOH and the anisotropy of the local charge distribution has been analyzed in terms of atomic multipoles calculated from the ab initio LCAO-MO-SCF wave function in the 6-31G* basis set. The results indicate that for longer -S-S-bonds the major contribution to these torsional barriers are electrostatic interactions of the atomic multipoles located on two atoms forming the rotated bond. This finding demonstrates the important role of electrostatic 1-2 interatomic interactions, usually neglected in conformational studies. It also opens the possibility to derive directly from accurate ab initio wave functions a simple nonempirical torsional potential involving atomic multipoles of two bonded atoms defining the torsional angle. For shorter -O-O- bonds, use of more precise models and inclusion of 1-3 interactions seems to be necessary.

NASA Program Exobiology↗
DOE OSTI2019

Surface Reconstruction Limited Conductivity in Block-Copolymer Li Battery Electrolytes

Solid polymer electrolytes for lithium batteries promise improvements in safety and energy density if their conductivity can be increased. Nanostructured block-copolymer electrolytes specifically have the potential to provide both good ionic conductivity and good mechanical properties. This study shows that the previously neglected nanoscale composition of the polymer electrolyte close to the electrode surface has an important effect on impedance measurements, despite its negligible extent compared to the bulk electrolyte. Using standard stainless steel blocking electrodes, the impedance of lithium salt-doped poly(isoprene-b-styrene-b-ethylene oxide) (ISO) exhibits a marked decrease upon thermal processing of the electrolyte. In contrast, covering the electrode surface with a low molecular weight poly(ethylene oxide) (PEO) brush results in higher and more reproducible conductivity values, which are insensitive to the thermal history of the device. A qualitative model of this effect is based on the hypothesis that ISO surface reconstruction at the different electrode surfaces leads to a change in the electrostatic double layer, affecting electrochemical impedance spectroscopy measurements. As a main result, PEO-brush modification of electrode surfaces is beneficial for the robust electrolyte performance of PEO-containing block-copolymers and may be crucial for their accurate characterization and use in Li-ion batteries.

block-copolymer electrolytes↗
DOE OSTI2026

Electrolyte Organization Leads to Potential-Dependence in Thermochemical Catalysis of Nonpolar Reactions

Electrochemical polarization is now known to play a key role in thermochemical catalysis at solid–liquid interfaces. However, existing frameworks cannot account for why even nonpolar, nonfaradaic reactions are sensitive to interfacial polarization. In order to uncover the molecular basis of this phenomenon, we herein study the potential-dependent reaction kinetics of ethylene and trans-2-butene hydrogenation at Pt–liquid interfaces. Measurements were performed in aqueous and ortho-difluorobenzene (o-DFB) solutions, spontaneously polarizing the Pt–liquid interfaces by, respectively, varying the pH or dissolving distinct metallocene redox buffers into solution. Here, we find that at comparable mechanistic regimes, the rates of both ethylene and trans-2-butene hydrogenation are maximized near the same electrochemical potential, E. Moreover, the potential-dependence, defined as $\frac{∂ln 𝑟}{∂𝐸}$, of trans-2-butene hydrogenation is approximately 2.2× greater than that of ethylene hydrogenation across the full potential range studied. These observations are all consistent with a model in which polarization of the Pt surface away from the local potential of zero free charge (E PZFC ) induces electrostatic organization of the polar solvent and charged ions near the interface, which impedes olefin adsorption and surface reaction because these surface reactions induce electrolyte displacement. Accordingly, interfacial polarization alters the free energy landscape and thus the rate of nonpolar heterogeneous catalysis by controlling the degree of electrostatic organization of polar and charged spectators at the interface, which do not in general need to be specifically chemisorbed onto the surface but could simply be close enough to the surface to be perturbed by the olefin adsorption. These results point toward electrochemical design handles, namely, the electrolyte, catalyst potential, and local E PZFC of the catalyst, with which to tune interfacial catalysis of thermochemical organic transformations.

adsorption↗
DOE OSTI2019

Insight into the Mechanisms Driving the Self-Assembly of Functional Interfaces: Moving from Lipids to Charged Amphiphilic Oligomers

Polymer-stabilized liquid/liquid interfaces are an important and growing class of bioinspired materials that combine the structural and functional capabilities of advanced synthetic materials with naturally evolved biophysical systems. These platforms have the potential to serve as selective membranes for chemical separations and molecular sequencers and to even mimic neuromorphic computing elements. Despite the diversity in function, basic insight into the assembly of well-defined amphiphilic polymers to form functional structures remains elusive, which hinders the continued development of these technologies. In this work, we provide new mechanistic insight into the assembly of an amphiphilic polymer-stabilized oil/aqueous interface, in which the headgroups consist of positively charged methylimidazolium ionic liquids, and the tails are short, monodisperse oligodimethylsiloxanes covalently attached to the headgroups. We demonstrate using vibrational sum frequency generation spectroscopy and pendant drop tensiometery that the composition of the bulk aqueous phase, particularly the ionic strength, dictates the kinetics and structures of the amphiphiles in the organic phase as they decorate the interface. These results show that H-bonding and electrostatic interactions taking place in the aqueous phase bias the grafted oligomer conformations that are adopted in the neighboring oil phase. The kinetics of self-assembly were ionic strength dependent and found to be surprisingly slow, being composed of distinct regimes where molecules adsorb and reorient on relatively fast time scales, but where conformational sampling and frustrated packing takes place over longer time scales. These results set the stage for understanding related chemical phenomena of bioinspired materials in diverse technological and fundamental scientific fields and provide a solid physical foundation on which to design new functional interfaces.

Chowdhury, Azhad↗
DOE OSTI2022

The influence of surface impurities on photoelectric currents driven by intense soft x rays

In an x-ray driven cavity experiment, an intense flux of soft x rays on the emitting surface produces significant emission of photoelectrons having several kiloelectronvolts of kinetic energy. At the same time, rapid heating of the emitting surface occurs, resulting in the release of adsorbed surface impurities and subsequent formation of an impurity plasma. This numerical study explores a simple model for the photoelectric currents and the impurity plasma. In this work, attention is given to the effect of varying the composition of the impurity plasma. The presence of protons or hydrogen molecular ions leads to a substantially enhanced cavity current, while heavier plasma ions are seen to have a limited effect on the cavity current due to their lower mobility. Additionally, it is demonstrated that an additional peak in the current waveform can appear due to the impurity plasma. A correlation between the impurity plasma composition and the timing of this peak is elucidated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

From Intermolecular Interaction Energies and Observable Shifts to Component Contributions and Back Again: A Tale of Variational Energy Decomposition Analysis

Quantum chemistry in the form of density functional theory (DFT) calculations is a powerful numerical experiment for predicting intermolecular interaction energies. However, no chemical insight is gained in this way beyond predictions of observables. Energy decomposition analysis (EDA) can quantitatively bridge this gap by providing values for the chemical drivers of the interactions, such as permanent electrostatics, Pauli repulsion, dispersion, and charge transfer. These energetic contributions are identified by performing DFT calculations with constraints that disable components of the interaction. This review describes the second-generation version of the absolutely localized molecular orbital EDA (ALMO-EDA-II). The effects of different physical contributions on changes in observables such as structure and vibrational frequencies upon complex formation are characterized via the adiabatic EDA. Example applications include red- versus blue-shifting hydrogen bonds; the bonding and frequency shifts of CO, N 2 , and BF bound to a [Ru(II)(NH 3 ) 5 ] 2 + moiety; and the nature of the strongly bound complexes between pyridine and the benzene and naphthalene radical cations. Additionally, the use of ALMO-EDA-II to benchmark and guide the development of advanced force fields for molecular simulation is illustrated with the recent, very promising, MB-UCB potential.

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

Electrostatic Gating of Ionic Conductance through Heterogeneous van der Waals Nanopores

Nanofluidic ionic transistors typically require gate voltages above 1 V and operate only at submillimolar ionic strengths, limiting their biocompatible applications. We demonstrate ionic transistors consisting of single sub-10 nm nanopores drilled in van der Waals (vdW) heterostructures with internal gate electrodes made of few-layer graphene. These devices deliver up to 10-fold current modulation at gate voltages as low as 0.3 V in 10 mM KCl, and ∼2-fold modulation at near-physiological 100 mM KCl. Baseline conductance with no gate shows surface-charge-dominated transport below 100 mM KCl, consistent with negatively charged hBN walls and ∼5 nm opening of the pores. The surface charge and the electrochemical asymmetry introduced by the three-electrode configuration govern the device’s behavior: negative gate voltage (V G ) enriches ionic concentrations and enhances current, whereas positive V G induces a local depletion zone that suppresses transport. The current modulation by V G is dependent on the polarity of the transmembrane potential and leads to ion current rectification. Molecular dynamics simulations of a nanopore in a hBN–graphene–hBN stack reveal confinement and surface charge-dependent suppression of the relative permittivity of interfacial water. Continuum modeling with radially varying interfacial water permittivity reproduces the asymmetric I–V characteristics and explains how the embedded gate sculpts local potential and ion concentrations. By enabling sub-0.5 V control of ionic transport at up to 100 mM salt concentrations, these devices address a key need in nanofluidics to create low-power ionic circuits and biosensing.

Materials science↗
DOE OSTI2023

Quantum mechanics based non-bonded force field functions for use in molecular dynamics simulations of materials and systems: The nitrogen and oxygen columns

Accurate Force Fields (FFs) are essential for Molecular Dynamics (MD) simulations of the dynamics of realistic materials in terms of atomic-level interactions. The FF parameters of short-range valence interactions can be derived through Quantum Mechanical (QM) calculations on model systems practical for QM (<300 atoms). Similarly, the dynamic electrostatic interactions can be described with methods such as QEq or PQEq that allow charges and polarization to adjust dynamically. However, accurately extracting long-range van der Waals (vdW) interactions from QM calculations poses challenges due to the absence of a definitive method to distinguish between the different energetic components of electrostatics, polarization, vdW, hydrogen bonding, and valence interactions. To do this we use the Perdew–Burke–Ernzerhof flavor of Density Functional Theory, including empirical D3 vdW corrections, to predict the Equation of State for each element (keeping any covalent bonds fixed), from which we obtain the two-body vdW nonbond potential. Here, we extend these calculations to include non-bonded parameters for the N and O columns of the periodic table so that we now describe columns 15 (N), 16 (O), 17 (F), and 18 (Ne) of the periodic table. For these 20 elements, we find that the two-body vdW potentials can all be mapped to a single universal two-body curve, with just three scaling parameters: Re, De, and L. We refer to this as the Universal NonBond (UNB) potential. We expect this to be useful for new MD simulations and a helpful starting point to obtain UNB parameters for the remainder of the periodic table.

Chemistry↗
DOE OSTI2023

Revisiting the pseudo-supercritical path method: An improved formulation for the alchemical calculation of solid–liquid coexistence

Alchemical free energy calculations via molecular dynamics have been applied to obtain thermodynamic properties related to solid–liquid equilibrium conditions, such as melting points. In recent years, the pseudo-supercritical path (PSCP) method has proved to be an important approach to melting point prediction due to its flexibility and applicability. In the present work, we propose improvements to the PSCP alchemical cycle to make it more compact and efficient through a concerted evaluation of different potential energies. The multistate Bennett acceptance ratio (MBAR) estimator was applied at all stages of the new cycle to provide greater accuracy and uniformity, which is essential concerning uncertainty calculations. In particular, for the multistate expansion stage from solid to liquid, we employed the MBAR estimator with a reduced energy function that allows affine transformations of coordinates. Free energy and mean derivative profiles were calculated at different cycle stages for argon, triazole, propenal, and the ionic liquid 1-ethyl-3-methyl-imidazolium hexafluorophosphate. Comparisons showed a better performance of the proposed method than the original PSCP cycle for systems with higher complexity, especially the ionic liquid. Furthermore, a detailed study of the expansion stage revealed that remapping the centers of mass of the molecules or ions is preferable to remapping the coordinates of each atom, yielding better overlap between adjacent states and improving the accuracy of the methodology.

25 ENERGY STORAGE↗
DOE OSTI2021

Polarization Effects in Water-Mediated Selective Cation Transport across a Narrow Transmembrane Channel

We report despite the progress in modeling complex molecular systems of ever-increasing complexity, a quantitatively accurate computational treatment of ion permeation through narrow membrane channels remains challenging. An important factor to reach this goal is induced electronic polarization, which is likely to impact the permeation rate of small ions through narrow molecular pores. In this work, we extended the recently developed polarizable force field based on the classical Drude oscillators to assess the role of induced polarization effects on the energetics of sodium and potassium ion transport across the gramicidin A (gA) ion channel. The inclusion of induced polarization lowers barriers present in 1D potential of mean force (PMF) for cation permeation by ~50% compared to those obtained with the additive force field. Conductance properties calculated with 1D PMFs from Drude simulations are in better agreement with experimental results. Polarization of single-file water molecules and protein atoms forming the narrow pore has a direct impact on the free-energy barriers and cation-specific solid-state NMR chemical shifts. Sensitivity analysis indicates that small changes to water–channel interactions can alter the free energy barrier for ion permeation. These results, illustrating polarization effects present in the complex electrostatic environment of the gA channel, have broad implications for revising proposed mechanisms of ion permeation and selectivity in a variety of ion channels.

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

Molecular Dynamics Simulation Study of the Protonation State Dependence of Glutamic Acid Transport through a Cyclic Peptide Nanotube

The effect of the protonation state of glutamic acid on its translocation through cyclic peptide nanotubes (CPNs) was assessed by using molecular dynamics (MD) simulations. Anionic (GLU–), neutral zwitterionic (GLU0), and cationic (GLU+) forms of glutamic acid were selected as three different protonation states for an analysis of energetics and diffusivity for acid transport across a cyclic decapeptide nanotube. Based on the solubility-diffusion model, permeability coefficients for the three protonation states of the acid were calculated and compared with experimental results for CPN-mediated glutamate transport through CPNs. Potential of mean force (PMF) calculations reveal that, due to the cation-selective nature of the lumen of CPNs, GLU–, so-called glutamate, shows significantly high free energy barriers, while GLU+ displays deep energy wells and GLU0 has mild free energy barriers and wells inside the CPN. The considerable energy barriers for GLU– inside CPNs are mainly attributed to unfavorable interactions with DMPC bilayers and CPNs and are reduced by favorable interactions with channel water molecules through attractive electrostatic interactions and hydrogen bonding. Unlike the distinct PMF curves, position-dependent diffusion coefficient profiles exhibit comparable frictional behaviors regardless of the charge status of three protonation states due to similar confined environments imposed by the lumen of the CPN. The calculated permeability coefficients for the three protonation states clearly demonstrate that glutamic acid has a strong protonation state dependence for its transport through CPNs, as determined by the energetics rather than the diffusivity of the protonation state. In addition, the permeability coefficients also imply that GLU– is unlikely to pass through a CPN due to the high energy barriers inside the CPN, which is in disagreement with experimental measurements, where a considerable amount of glutamate permeating through the CPN was detected. To resolve the discrepancy between this work and the experimental observations, several possibilities are proposed, including a large concentration gradient of glutamate between the inside and outside of lipid vesicles and bilayers in the experiments, the glutamate activity difference between our MD simulations and experiments, an overestimation of energy barriers due to the artifacts imposed in MD simulations, and/or finally a transformation of the protonation state from GLU– to GLU0 to reduce the energy barriers. Altogether, our study demonstrates that the protonation state of glutamic acid has a strong effect on the transport of the acid and suggests a possible protonation state change for glutamate permeating through CPNs.

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

The importance of localized modes spectral contribution to thermal conductivity in amorphous polymers

Abstract Polymers are a unique class of materials from the perspective of normal mode analysis. Polymers consist of individual chains with repeating units and strong intra-chain covalent bonds, and amorphous arrangements among chains with weak inter-chain van der Waals and for some polymers also electrostatic interactions. Intuitively, this strong heterogeneity in bond strength can give rise to special features in the constituent phonons, but such effects have not been studied deeply before. Here, we use lattice dynamics and molecular dynamics to perform modal analysis of the thermal conductivity in amorphous polymers. We find an abnormally large population of localized modes in amorphous polymers, which is fundamentally different from amorphous inorganic materials. Contrary to the common picture of thermal transport, localized modes in amorphous polymers are found to be the dominant contributors to thermal conductivity. We find that a significant portion of the localization happens within individual chains, but heat is dominantly conducted when localized modes involve two chains. These results suggest localized modes generally play a key role in thermal transport for different polymers. The results provide an alternative perspective on why polymer thermal conductivity is generally quite low and gives insight into how to potentially change it.

36 MATERIALS SCIENCE↗