Using Constrained Density Functional Theory to Track Proton Transfers and to Sample Their Associated Free Energy Surface
Not Available
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
SEARCH · Search NASA
Results for “surface charge dynamics”
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.
Substrate-Directed Dimensional and Phase Control of Peptide Assemblies on Two-Dimensional van der Waals Materials
Understanding and controlling biomolecular self-assembly on van der Waals (vdW) materials has the potential to advance hybrid bioelectronic devices by enabling precise tuning of the interface and modulation of the resulting electronic properties of the biomolecule-vdW heterostructure. However, how surface properties of vdW materials direct biomolecule assembly remains poorly understood. To fill this knowledge gap, we investigated the assembly of a peptide known to assemble into two-dimensional (2D) crystalline films on MoS 2 on three representative vdW surfaces: WS 2 , MoS 2 , and highly oriented pyrolytic graphite (HOPG). Using in situ atomic force microscopy (AFM), we find that assembly is substrate-dependent, resulting in multilayers on WS 2 , monolayers on MoS 2 , and multiple coexisting phases on HOPG. WS 2 exhibits a higher negative charge, strong long-range electrostatic interactions, and extensive hydration layering that may promote multilayer stacking. In contrast, MoS 2 has stronger short-range interactions with the peptides but much weaker long-range interactions and hydration structure, which may favor monolayer formation. Molecular dynamics simulations predict a corresponding switch from monolayer to multilayer aggregates of the adsorbed monomers, reflected in their relative mobilities. On hydrophobic HOPG, the peptides bind most strongly and remain as monomers with high surface mobility. The peptide dimers comprising the basic unit of the crystals are more compact on HOPG, which has a smaller lattice constant than WS 2 or MoS 2 , suggesting strain contributes to stabilizing multiple phases. Our results provide mechanistic insights into how surface charge and hydration structure, and the lattice structure of the substrates governs peptide assembly on vdW materials, offering a framework to rationally control the 2D peptide-vdW heterostructures.
Molecular dynamics simulations↗
Lanthanide binding peptide surfactants at air–aqueous interfaces for interfacial separation of rare earth elements
Rare earth elements (REEs) are critical materials to modern technologies. They are obtained by selective separation from mining feedstocks consisting of mixtures of their trivalent cation. We are developing an all-aqueous, bioinspired, interfacial separation using peptides as amphiphilic molecular extractants. Lanthanide binding tags (LBTs) are amphiphilic peptide sequences based on the EF-hand metal binding loops of calcium-binding proteins which complex selectively REEs. We study LBTs optimized for coordination to Tb 3+ using luminescence spectroscopy, surface tensiometry, X-ray reflectivity, and X-ray fluorescence near total reflection, and find that these LBTs capture Tb 3+ in bulk and adsorb the complex to the interface. Molecular dynamics show that the binding pocket remains intact upon adsorption. We find that, if the net negative charge on the peptide results in a negatively charged complex, excess cations are recruited to the interface by nonselective Coulombic interactions that compromise selective REE capture. If, however, the net negative charge on the peptide is −3, resulting in a neutral complex, a 1:1 surface ratio of cation to peptide is achieved. Surface adsorption of the neutral peptide complexes from an equimolar mixture of Tb 3+ and La 3+ demonstrates a switchable platform dictated by bulk and interfacial effects. The adsorption layer becomes enriched in the favored Tb 3+ when the bulk peptide is saturated, but selective to La 3+ for undersaturation due to a higher surface activity of the La 3+ complex.
Ortuno Macias, Luis E. (ORCID:0000000284342192)↗
Extended Lagrangian Born–Oppenheimer molecular dynamics for orbital-free density-functional theory and polarizable charge equilibration models
We report extended Lagrangian Born–Oppenheimer molecular dynamics (XL-BOMD) is formulated for orbital-free Hohenberg–Kohn density-functional theory and for charge equilibration and polarizable force-field models that can be derived from the same orbital-free framework. The purpose is to introduce the most recent features of orbital-based XL-BOMD to molecular dynamics simulations based on charge equilibration and polarizable force-field models. These features include a metric tensor generalization of the extended harmonic potential, preconditioners, and the ability to use only a single Coulomb summation to determine the fully equilibrated charges and the interatomic forces in each time step for the shadow Born–Oppenheimer potential energy surface. The orbital-free formulation has a charge-dependent, short-range energy term that is separate from long-range Coulomb interactions. This enables local parameterizations of the short-range energy term, while the long-range electrostatic interactions can be treated separately. The theory is illustrated for molecular dynamics simulations of an atomistic system described by a charge equilibration model with periodic boundary conditions. The system of linear equations that determines the equilibrated charges and the forces is diagonal, and only a single Ewald summation is needed in each time step. The simulations exhibit the same features in accuracy, convergence, and stability as are expected from orbital-based XL-BOMD.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Giant quantum oscillations in thermal transport in low-density metals via electron absorption of phonons
Oscillations of conductance observed in strong magnetic fields are a striking manifestation of the quantum dynamics of charge carriers in solids. The large charge carrier density in typical metals sets the scale of oscillations in both electrical and thermal conductivity, which characterize the Fermi surface. In semimetals, thermal transport at low-charge carrier density is expected to be phonon dominated, yet several experiments observe giant quantum oscillations in thermal transport. This raises the question of whether there is an overarching mechanism leading to sizable oscillations that survives in phonon-dominated semimetals. In this work, we show that such a mechanism exists. It relies on the peculiar phase-space allowed for phonon scattering by electrons when only a few Landau levels are filled. Our measurements on the Dirac semimetal ZrTe 5 support this counterintuitive mechanism through observation of pronounced thermal quantum oscillations, since they occur in similar magnitude and phase in directions parallel and transverse to the magnetic field. Our phase-space argument applies to all low-density semimetals, topological or not, including graphene and bismuth. Our work illustrates that phonon absorption can be leveraged to reveal degrees of freedom through their imprint on longitudinal thermal transport.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Orientation-dependent etching of silicon by fluorine molecules: A quantum chemistry computational study
Anisotropic etching is a widely used process in semiconductor manufacturing, in particular, for micro- and nanoscale texturing of silicon surfaces for black silicon production. The typical process of plasma-assisted etching uses energetic ions to remove materials in the vertical direction, creating anisotropic etch profiles. Plasmaless anisotropic etching, considered here, is a less common process that does not use ions and plasma. The anisotropy is caused by the unequal etching rates of different crystal planes; the etching process, thus, proceeds in a preferred direction. In this paper, we have performed quantum chemistry modeling of gas-surface reactions involved in the etching of silicon surfaces by molecular fluorine. The results confirm that orientation-dependent etch rates are the reason for anisotropy. The modeling of F 2 dissociative chemisorption on F-terminated silicon surfaces shows that Si–Si bond breaking is slow for the Si(111) surface, while it is fast for Si(100) and Si(110) surfaces. Both Si(100) and Si(110) surfaces incorporate a larger number of fluorine atoms resulting in Si–Si bonds having a larger amount of positive charge, which lowers the reaction barrier of F 2 dissociative chemisorption, yielding a higher etch rate for Si(100) and Si(110) surfaces compared to Si(111) surfaces. Molecular dynamics modeling of the same reactions has shown that the chosen reactive bond order potential does not accurately reproduce the lower reaction barriers for F 2 dissociative chemisorption on Si(100) and Si(100) surfaces. Furthermore, reparameterization is necessary to model the anisotropic etching process that occurs at lower temperatures.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Cryogenic spatial–temporal imaging of surface photocarrier dynamics in MAPbI 3 films at the single grain level
We use cryogenic spatial–temporal photoluminescence (PL) imaging measurements down to 10 K and with short-wavelength, 405 nm laser excitation to study surface charge generation, trapping, and recombination at single bright and dark grains as well as their boundaries in model methylammonium lead iodide (MAPbI 3 ) polycrystalline thin films. These salient conditions are shown to be critical for identifying both the detrimental and cooperative roles of grain microstructures where the dark grains serve as the PL quenching center, while the grain boundaries are largely benign and may promote electron–hole separation.
36 MATERIALS SCIENCE↗
Solvent-Dependent Dynamics of Cellulose Nanocrystals in Process-Relevant Flow Fields
Flow-assisted alignment of anisotropic nanoparticles is a promising route for the bottom-up assembly of advanced materials with tunable properties. While aligning processes could be optimized by controlling factors such as solvent viscosity, flow deformation, and the structure of the particles themselves, it is necessary to understand the relationship between these factors and their effect on the final orientation. In this study, we investigated the flow of surface-charged cellulose nanocrystals (CNCs) with the shape of a rigid rod dispersed in water and propylene glycol (PG) in an isotropic tactoid state. In situ scanning small-angle X-ray scattering (SAXS) and rheo-optical flow-stop experiments were used to quantify the dynamics, orientation, and structure of the assigned system at the nanometer scale. The effects of both shear and extensional flow fields were revealed in a single experiment by using a flow-focusing channel geometry, which was used as a model flow for nanomaterial assembly. Due to the higher solvent viscosity, CNCs in PG showed much slower Brownian dynamics than CNCs in water and thus could be aligned at lower deformation rates. Moreover, CNCs in PG also formed a characteristic tactoid structure but with less ordering than CNCs in water owing to weaker electrostatic interactions. The results indicate that CNCs in water stay assembled in the mesoscale structure at moderate deformation rates but are broken up at higher flow rates, enhancing rotary diffusion and leading to lower overall alignment. Albeit being a study of cellulose nanoparticles, the fundamental interplay between imposed flow fields, Brownian motion, and electrostatic interactions likely apply to many other anisotropic colloidal systems.
36 MATERIALS SCIENCE↗
Molecular Dynamics Simulation and Cryo-Electron Microscopy Investigation of AOT Surfactant Structure at the Hydrated Mica Surface
Structural properties of the anionic surfactant dioctyl sodium sulfosuccinate (AOT or Aerosol-OT) adsorbed on the mica surface were investigated by molecular dynamics simulation, including the effect of surface loading in the presence of monovalent and divalent cations. The simulations confirmed recent neutron reflectivity experiments that revealed the binding of anionic surfactant to the negatively charged surface via adsorbed cations. At low loading, cylindrical micelles formed on the surface, with sulfate head groups bound to the surface by water molecules or adsorbed cations. Cation bridging was observed in the presence of weakly hydrating monovalent cations, while sulfate groups interacted with strongly hydrating divalent cations through water bridges. The adsorbed micelle structure was confirmed experimentally with cryogenic electronic microscopy, which revealed micelles approximately 2 nm in diameter at the basal surface. At higher AOT loading, the simulations reveal adsorbed bilayers with similar surface binding mechanisms. Adsorbed micelles were slightly thicker (2.2–3.0 nm) than the corresponding bilayers (2.0–2.4 nm). Upon heating the low loading systems from 300 K to 350 K, the adsorbed micelles transformed to a more planar configuration resembling bilayers. The driving force for this transition is an increase in the number of sulfate head groups interacting directly with adsorbed cations.
58 GEOSCIENCES↗
Tailored PEDOT:PSS hole transport layer for higher performance in perovskite solar cells: Enhancement of electrical and optical properties with improved morphology
Precise control over the charge carrier dynamics throughout the device can result in outstanding performance of perovskite solar cells (PSCs). Poly(3,4-ethylenedioxythiophene):poly(styrene sulfonate) (PEDOT:PSS) is the most actively studied hole transport material in p-i-n structured PSCs. However, charge transport in the PEDOT:PSS is limited and inefficient because of its low conductivity with the presence of the weak ionic conductor PSS. In addition, morphology of the underlying PEDOT:PSS layer in PSCs plays a crucial role in determining the optoelectronic quality of the active perovskite absorber layer. This work is focused on realization of a non-wetting conductive surface of hole transport layer suitable for the growth of larger perovskite crystalline domains. This is accomplished by employing a facile solvent-engineered (ethylene glycol and methanol) approach resulting in removal of the predominant PSS in PEDOT:PSS. The consequence of acquiring larger perovskite crystalline domains was observed in the charge carrier dynamics studies, with the achievement of higher charge carrier lifetime, lower charge transport time and lower transfer impedance in the solvent-engineered PEDOT:PSS-based PSCs. Use of this solvent-engineered treatment for the fabrication of MAPbI3 PSCs greatly increased the device stability witnessing a power conversion efficiency of 18.18%, which corresponds to ~37% improvement compared to the untreated PEDOT:PSS based devices.
14 SOLAR ENERGY↗
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↗
Operando Depth-Resolved Measurement of Solvation Entropy, Interfacial Transport, and Charge-Transfer Kinetics in Lithium-Ion Batteries
Understanding and improving the performance and longevity of lithium-ion batteries critically depends on insight into the dynamic processes occurring at buried electrode-electrolyte interfaces. However, direct, depth-resolved, and operando diagnosis of these interfaces remains a longstanding challenge due to their inaccessibility beneath bulk materials, the limitations of conventional surface- and bulk-sensitive characterization tools, and the difficulty of maintaining realistic cell environments during measurement. These challenges have made it nearly impossible to uniquely resolve important interfacial properties such as charge transfer resistance, SEI (solid electrolyte interphase) resistance, and solvation entropy at the individual electrode interfaces within a working cell, information that is essential for mechanistic insight and accelerated battery design. Here, in this study, we report the development of multiharmonic electro-thermal spectroscopy (METS), an operando technique that enables depth-resolved measurement of solvation entropy, interfacial transport resistance, charge-transfer resistance, and SEI resistance at individual electrode-electrolyte interfaces within practical lithium-ion batteries. By leveraging frequency-dependent, thermal-wave sensing and interface-specific modeling, METS uniquely attributes interfacial properties to specific electrodes, as validated by comparison with traditional electrochemical impedance spectroscopy (EIS). The ability to spatially and temporally resolve interfacial processes in real time provides new diagnostic capabilities that are crucial for mechanistic studies of battery degradation and for the rapid development of next-generation energy storage systems.
Chalise, Divya [University of California, Berkeley↗
Theoretical model of dielectrophoretic water vapor condensation during multiphase flow of electrically charged water droplets in air
Water droplet growth in air, influenced by vapor diffusion and electrical charging, has significant implications for indoor air quality, cloud formation, and aerosol dynamics. Here, this study presents a comprehensive theoretical model that describes the growth dynamics of electrically charged water droplets under external electric fields, incorporating diffusion, electrostatic, and dielectrophoretic forces to enhance understanding and practical applications in air dehumidification and environmental control. A theoretical model was established to analyze the growth dynamics of the charged water droplets. The approach extended the thermodynamic framework of the Modified Kelvin-Thomson equation by integrating electrostatic and dielectrophoretic effects, enabling a comprehensive evaluation of vapor concentration near the surface of charged droplets. A sensitivity analysis was performed to assess the roles of droplet size, electric field intensity, surface charge, and vapor concentration on condensation rates. The study progressed beyond the classical diffusion-limited growth of neutral droplets, addressing the behavior of charged droplets in electrohydrodynamic environments. The model highlighted the significance of dielectrophoretic drift flux for dipolar vapor molecules, especially in micron-sized droplets ranging from 10 to 100 µm. Findings demonstrated that, for smaller droplets subjected to elevated electrical charge, average growth rates varied from 0.3 to 1.6 µm/s if the air relative humidity increased from 50 to 100 %. Analytical results are corroborated by numerical simulations, which exhibited strong concordance, thereby validating the proposed theoretical approach. This research introduced an innovative and scalable analytical methodology for electric field-enhanced droplet growth, offering valuable insights into the advancement of energy-efficient air dehumidification, water desalination, and atmospheric water harvesting technologies.
42 ENGINEERING↗
First Principles Nonadiabatic Excited-State Molecular Dynamics in NWChem
Computational simulation of non-adiabatic molecular dynamics is an indispensable tool for understanding complex photoinduced processes such as internal conversion, energy transfer, charge separation, and spatial localization of excitons, to name a few. Here, we report an implementation of the fewest-switches surface hopping algorithm in the NWChem computational chemistry program. Here the surface hopping method is combined with linear-response time-dependent density functional theory calculations of adiabatic excited state potential energy surfaces. To treat quantum transitions between arbitrary electronic Born{Oppenheimer states, we have implemented both numerical and analytical differentiation schemes for derivative non-adiabatic couplings. A numerical approach for the time-derivative non-adiabatic couplings together with an analytical method for calculating non-adiabatic coupling vectors is an efficient combination for surface hopping approaches. Additionally, electronic decoherence schemes and a state reassigned unavoided crossings algorithm are also implemented to improve the accuracy of the simulated dynamics and to handle trivial unavoided crossings. We apply our code to study the ultrafast decay of photoexcited benzene, including a detailed analysis of the potential energy surface, population decay time scales, and vibrational coordinates coupled to the excitation dynamics. The development presented in this work is a baseline for future implementations of more sophisticated frameworks for simulating of non-adiabatic molecular dynamics in NWChem.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Femtosecond dynamics and coherence of ionic retro-Diels–Alder reactions
Ultrafast tunnel ionization enables femtosecond time-resolved dynamic measurements of the retro-Diels–Alder reactions of positively charged cyclohexene, norbornene, and dicyclopentadiene. Unlike the reaction times of 500–600 ps that are observed following UV excitation of neutral species, on the ionic potential energy surfaces, these reactions occur on a single picosecond timescale and, in some cases, exhibit vibrational coherence. In the case of norbornene, a 270 cm –1 vibrational mode is found to modulate the retro-Diels–Alder reaction.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Role of generic scale invariance in a Mott transition from a U(1) spin-liquid insulator to a Landau Fermi-liquid metal
Highlights: • We look at the role played by generic scale invariance (GSI) in the phase transition between a U(1) spin liquid and a Fermi liquid. • We identify dangerously irrelevant operators in the renormalization group transformation, involved with the GSI. • We highlight novel mean-field dynamics of this Mott quantum criticality, resulting from the GSI. We investigate the role of generic scale invariance in a Mott transition from a U(1) spin-liquid insulator to a Landau Fermi-liquid metal, where there exist massless degrees of freedom in addition to quantum critical fluctuations. Here, the Mott quantum criticality is described by critical charge fluctuations, and additional gapless excitations are U(1) gauge-field fluctuations coupled to a spinon Fermi surface in the spin-liquid state, which turn out to play a central role in the Mott transition. An interesting feature of this problem is that the scaling dimension of effective leading local interactions between critical charge fluctuations differs from that of the coupling constant between U(1) gauge fields and matter-field fluctuations in the presence of a Fermi surface. As a result, there appear dangerously irrelevant operators, which can cause conceptual difficulty in the implementation of renormalization group (RG) transformations. Indeed, we find that the curvature term along the angular direction of the spinon Fermi surface is dangerously irrelevant at this spin-liquid Mott quantum criticality, responsible for divergence of the self-energy correction term in U(1) gauge-field fluctuations. Performing the RG analysis in the one-loop level based on the dimensional regularization method, we reveal that such extremely overdamped dynamics of U(1) gauge-field fluctuations, which originates from the emergent one-dimensional dynamics of spinons, does not cause any renormalization effects to the effective dynamics of both critical charge fluctuations and spinon excitations. However, it turns out that the coupling between U(1) gauge-field fluctuations and both matter-field excitations still persists at this Mott transition, which results in novel mean-field dynamics to explain the nature of the spin-liquid Mott quantum criticality. We discuss physical implications of effective one-dimensional spin dynamics and extremely overdamped gauge dynamics at the Mott quantum criticality.
36 MATERIALS SCIENCE↗
Dynamic Ligand Surface Chemistry of Excited PbS Quantum Dots
The ligand shell around colloidal quantum dots mediates the electron and energy transfer processes that underpin their use in optoelectronic and photocatalytic applications. In this work, we show that the surface chemistry of carboxylate anchoring groups of oleate ligands passivating PbS quantum dots undergoes significant changes when the quantum dots are excited to their excitonic states. Additionally, we directly probe the changes of surface chemistry using time-resolved mid-infrared spectroscopy that records the evolution of the vibrational frequencies of carboxylate groups following excitation of the electronic states. The data reveal a reduction of the Pb-O coordination of carboxylate anchoring groups to lead atoms at the quantum dot surfaces. The dynamic surface chemistry of the ligands may increase their surface mobility in the excited state and enhance the ability of molecular species to penetrate the ligand shell to undergo energy and charge transfer processes that depend sensitively on distance.
36 MATERIALS SCIENCE↗
Spatiotemporal Imaging of Thickness-Induced Band-Bending Junctions
van der Waals materials exhibit naturally passivated surfaces and an ability to form versatile heterostructures to enable an examination of carrier transport mechanisms not seen in traditional materials. Here, we report a new type of homojunction termed a "band-bending junction" whose potential landscape depends solely on the difference in thickness between the two sides of the junction. Using MoS 2 on Au as a prototypical example, we find that surface potential differences can arise from the degree of vertical band bending in thin and thick regions. Furthermore, by using scanning ultrafast electron microscopy, we examine the spatiotemporal dynamics of charge carriers generated at this junction and find that lateral carrier separation is enabled by differences in the band bending in the vertical direction, which we verify with simulations. As a result, band-bending junctions may therefore enable new optoelectronic devices that rely solely on band bending arising from thickness variations to separate charge carriers.
77 NANOSCIENCE AND NANOTECHNOLOGY↗