Engineering Papers⌕ Search

Engineering topics

Schwartz, Benjamin J.

Publications and source records attributed to Schwartz, Benjamin J..

Evaluating the Chemical Reactivity of DFT-Simulated Liquid Water with Hydrated Electrons via the Dual Descriptor

Modeling the various properties of liquid water, particularly its reactivity, has been a longstanding challenge for simulation methods. Recently, ab initio simulations based on density functional theory (DFT) have come to the fore as tenable methods for calculating the properties and reactivity of water, with varying degrees of success for different exchange-correlation functionals. In particular, hybrid-GGA and meta-GGA functionals have been shown to reproduce many of the structural, dynamical, and energetic properties of water to a high degree of accuracy relative to their computational cost. Here, we show that the dual descriptor (DD) measure of nucleophilicity and electrophilicity, which is sometimes used to elucidate organic chemistry reaction mechanisms, can also be used to characterize the reactivity of DFTsimulated liquid water. The DD is especially apt for understanding the reactivity of excess electrons with water as its calculation explicitly involves adding and removing an excess electron from a reference system. We use the DD to explore the reactivity of water simulated using three different DFT functionals: the LDA functional (LDA), a hybrid-GGA functional (PBE0), and a hybrid meta- GGA functional (SCAN0). Using the DD, we show that the SCAN0 functional with the standard 25% Hartree–Fock exchange produces simulated liquid water with many regions that are far more reactive than either PBE0 or LDA. To understand the implications of these highly reactive regions, we then add a strong nucleophile in the form of an excess electron and find that although PBE0 and LDA predict stable hydrated electrons, the excess electron reacts nearly instantaneously with SCAN0 water via proton abstraction to form a hydrogen atom and hydroxide ion. We show that the DD provides the ability to not only predict whether or not liquid water will react with a hydrated electron but also which particular waters will be involved solely from analyzing pure water configurations generated with each functional. We rationalize this result in terms of the known trap-seeking behavior of injected hydrated electrons, which are able to find the most electronegative region in bulk water. These results highlight the utility of the dual descriptor as a fast and interpretable method for investigating condensed-phase reactivity with excess electrons.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Tuning counterion chemistry to reduce carrier localization in doped semiconducting carbon nanotube networks

Understanding and controlling the impact that electrostatic interactions have on the transport of injected charge carriers is important for the utilization of pi-conjugated semiconductors in opto-electronic applications. Here, we explore the impact of dopant chemical and electronic structure on the doping efficacy and charge carrier transport in semiconducting single-walled carbon nanotube (s-SWCNT) networks using molecular charge-transfer dopants based on functionalized icosahedral dodecaborane (DDB) clusters. Calculations indicate that localization of electron density on the DDB core reduces the coulombic interactions that contribute to hole localization in the s-SWCNTs, thereby improving charge carrier transport. The enhanced delocalization produces an increase in the electrical conductivity and thermopower at lower charge carrier densities, yielding enhanced thermoelectric transport and a thermoelectric power factor that surpasses the previous best in class for enriched s-SWCNT thin-film networks. This strategy can be applied broadly across pi-conjugated semiconductors to tune and enhance performance in a variety of energy harvesting devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Counterion Control and the Spectral Signatures of Polarons, Coupled Polarons, and Bipolarons in Doped P3HT Films

When an electron is removed from a conjugated polymer, such as poly(3-hexylthiophene-2,5-diyl) (P3HT), the remaining hole and associated change in the polymer backbone structure from aromatic to quinoidal are referred to as a polaron. Bipolarons are created by removing the unpaired electron from an already-oxidized polymer segment. In electrochemically-doped P3HT films, polarons, and bipolarons are readily observed, but in chemically-doped P3HT films, bipolarons rarely form. This is explained by studying the effects of counterion position on the formation of polarons, strongly coupled polarons, and bipolarons using both spectroscopic and X-ray diffraction experiments and time-dependent density functional theory calculations. The counterion positions control whether two polarons spin-pair to form a bipolaron or whether they strongly couple without spin-pairing are found. When two counterions lie close to the same polymer segment, bipolarons can form, with an absorption spectrum that is blueshifted from that of a single polaron. Otherwise, polarons at high concentrations do not spin-pair, but instead J-couple, leading to a redshifted absorption spectrum. The counterion location needed for bipolaron formation is accompanied by a loss of polymer crystallinity. Finally, these results explain the observed formation order of single polarons, coupled single polarons, and singlet bipolarons in electrochemically- and chemically-doped conjugated polymers.

36 MATERIALS SCIENCE↗

Nonequilibrium Solvent Effects during Photodissociation in Liquids: Dynamical Energy Surfaces, Caging and Chemical Identity

In the gas phase, potential energy surfaces can be used to provide insight into the details of photochemical reaction dynamics. In solution, however, it is unclear what potential energy surfaces, if any, can be used to describe even simple chemical reactions such as the photodissociation of a diatomic solute. In this paper, we use mixed quantum/classical (MQC) molecular dynamics (MD) to study the photodissociation of $Na$ $^+_2$ in both liquid Ar and liquid tetrahydrofuran (THF). We examine both the gas-phase potential surfaces and potentials of mean force (PMF), which assume that the solvent remains at equilibrium with the solute throughout the photodissociation process and show that neither resemble a nonequilibrium dynamical energy surface that is generated by taking the time integral of work. For the photodissociation of $Na$ $^+_2$ in liquid Ar, the dynamical energy surface shows clear signatures of solvent caging, and the degree of caging is directly related to the mass of the solvent atoms. For $Na$ $^+_2$ in liquid THF, local specific interactions between the solute and solvent lead to changes in chemical identity that create a kinetic trap that effectively prevents the molecule from dissociating. Here, the results show that nonequilibrium effects play an important role even in simple solution-phase reactions, requiring the use of dynamical energy surface to understand such chemical events.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Role of the Solvent in the Condensed-Phase Dynamics and Identity of Chemical Bonds: The Case of the Sodium Dimer Cation in THF

When a solute molecule is placed in solution, is it acceptable to presume that its electronic structure is essentially the same as in the gas phase? In this paper, we address this question from a simulation perspective for the case of the sodium dimer cation (Na 2 + ) molecule in both liquid Ar and liquid tetrahydrofuran (THF). In previous work, we showed that when local specific interactions between a solute and solvent are energetically on the order of a hydrogen bond, the solvent can become part of the chemical identity of the solute. Here, using mixed quantum/classical molecular dynamics simulations, we see that for the Na 2 + molecule, solute--solvent interactions lead to two stable, chemically-distinct coordination states (Na(THF) 4 -Na(THF) 5 + and Na(THF) 5 -Na(THF) 5 + ) that are not only stable themselves as gas-phase molecules but that also have a completely new electronic structure with important implications for the excited-state photodissociation of this molecule in the condensed phase. Furthermore, we show through a set of comparative classical simulations that treating the solute's bonding electron explicitly quantum mechanically is necessary to understand both the ground-state dynamics and chemical identity of this simple diatomic molecule; even use of the quantum-derived potential of mean force is insufficient to describe the behavior of the molecule classically. Lastly, we calculate the results of a proposed transient hole-burning experiment that could be used to spectroscopically disentangle the presence of the different coordination states.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Tunable Dopants with Intrinsic Counterion Separation Reveal the Effects of Electron Affinity on Dopant Intercalation and Free Carrier Production in Sequentially Doped Conjugated Polymer Films

Carrier mobility in doped conjugated polymers is limited by Coulomb interactions with dopant counterions. This complicates studying the effect of the dopant’s oxidation potential on carrier generation because different dopants have different Coulomb interactions with polarons on the polymer backbone. In this work, dodecaborane (DDB)-based dopants are used, which electrostatically shield counterions from carriers and have tunable redox potentials at constant size and shape. DDB dopants produce mobile carriers due to spatial separation of the counterion, and those with greater energetic offsets produce more carriers. Neutron reflectometry indicates that dopant infiltration into conjugated polymer films is redox-potential-driven. Remarkably, X-ray scattering shows that despite their large 2-nm size, DDBs intercalate into the crystalline polymer lamellae like small molecules, indicating that this is the preferred location for dopants of any size. These findings elucidate why doping conjugated polymers usually produces integer, rather than partial charge transfer: dopant counterions effectively intercalate into the lamellae, far from the polarons on the polymer backbone. Finally, it is shown that the IR spectrum provides a simple way to determine polaron mobility. Overall, higher oxidation potentials lead to higher doping efficiencies, with values reaching 100% for driving forces sufficient to dope poorly crystalline regions of the film.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗