Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Implicit”

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 1,801 records · Page 100

Quantum Chemical Prediction of p K a Values of Cationic Ion-Exchange Groups in Polymer Electrolyte Membranes

The durability of anion-exchange membranes (AEMs) in alkaline fuel cells is largely determined by the stability of anion-exchange groups. Despite continuous research efforts, the commonly employed cations still have limited stability against hydroxide that can act as a strong base and nucleophile. This work is concerned with base-catalyzed degradation of organic cations initiated by proton abstraction to form reactive ylides or carbenes. We report on the performance of 24 methods combining density functional theory and electronic structure methods with implicit solvation calculations for predicting pK a values of organic cations in water and DMSO. The most accurate computational protocols are obtained using a combination of M06-2X/6-311++G** with the SMD solvation model for water (the mean absolute error of 0.4 pK a units) and B3LYP/aug-cc-pVTZ with the IEFPCM solvation model for DMSO (the mean absolute error of 1.4 pK a units). The aqueous pK a calculation protocol is cross-validated against the experimental C–H acidity constants outside the conventional range of 0–14 pK a values. In conclusion, this study rationalizes alkaline degradation of imidazolium cations with C2-alkyl substituents and provides a theoretical scale of C–H acidity for potential anion-exchange groups in AEMs.

25 ENERGY STORAGE↗

Computational Investigation of Aqueous Phase Effects on the Dehydrogenation and Dehydroxylation of Polyols over Pt(111)

Prediction of solvation effects on the kinetics of elementary reactions occurring at metal–water interfaces is of high importance for the rational design of catalysts for the biomass and electrocatalysis communities. A lack of knowledge of the reliability of various computational solvation schemes for processes on metal surfaces is currently a limiting factor. Using a multilevel quantum mechanical/molecular mechanical (QM/MM) description of the potential energy surface, we determined characteristic time and length scales for typical free-energy perturbation (FEP) calculations of bond cleavages in ethylene glycol, a sugar surrogate molecule, over Pt(111). Our approach is based on our explicit solvation model for metal surfaces and the repetition of FEP calculations to estimate confidence intervals. Results indicate that aqueous phase effects on the free energies of elementary processes can be determined with 95% confidence intervals from limited configuration space sampling and the fixed charge approximation used in the QM/MM-FEP methodology of smaller 0.1 eV. Next, we computed the initial O–H, C–H, and C–OH bond cleavages in ethylene glycol over Pt(111) in liquid water utilizing two different metal–water interaction potentials. Our calculations predict that aqueous phase effects are small (<0.1 eV) for the C–H bond cleavage and the activation barrier of the C–OH bond cleavage. In contrast, solvation effects are large (>0.35 eV) for the O–H bond cleavage and the reaction free energy of the C–OH bond scission. While the choice of a different Pt–water force field can lead to differences in predicted solvation effects of up to 0.2 eV, the differences are usually smaller (<0.1 eV), and the trends are always the same. In contrast, implicit solvation methods appear to currently not be able to reliably describe solvation effects originating from hydrogen bonding for metal surfaces even qualitatively.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Oxygen Vacancy Linear Clustering in a Perovskite Oxide

Oxygen vacancies have been implicitly assumed isolated ones, and understanding oxide materials possibly containing oxygen vacancies remains elusive within the scheme of the isolated vacancies, although the oxygen vacancies have been playing a decisive role in oxide materials. We report the presence of oxygen vacancy linear clusters and their orientation along a specific crystallographic direction in SrTiO 3 , a representative of a perovskite oxide. The presence of the linear clusters and associated electron localization was revealed by an electronic structure represented in the increase in the Ti 2+ valence state or corresponding Ti 3d 2 electronic configuration along with divacancy cluster model analysis and transport measurement. The orientation of the linear clusters along the [001] direction in perovskite SrTiO 3 was verified by further X-ray diffuse scattering analysis. And because SrTiO 3 is an archetypical perovskite oxide, the vacancy linear clustering with the specific aligned direction and electron localization can be extended to a wide variety of the perovskite oxides.

36 MATERIALS SCIENCE↗

Multireference Density Functional Theory with Generalized Auxiliary Systems for Ground and Excited States

To describe static correlation, we create a new approach to density functional theory (DFT), which uses a generalized auxiliary system that is of a different symmetry, such as particle number or spin, from that of the physical system. The total energy of the physical system consists of two parts: the energy of the auxiliary system, which is determined with a chosen density functional approximation (DFA), and the excitation energy from an approximate linear response theory that restores the symmetry to that of the physical system, thus rigorously leading to a multideterminant description of the physical system. The electron density of the physical system is different from that of the auxiliary system and is uniquely determined from the functional derivative of the total energy with respect to the external potential. Our energy functional is thus an implicit functional of the physical system density, but an explicit functional of the auxiliary system density. We demonstrate that the total energy minimum and stationary states, describing the ground and excited states of the physical system, can be obtained by a self-consistent optimization with respect to the explicit variable, the generalized Kohn–Sham noninteracting density matrix. We have developed the generalized optimized effective potential method for the self-consistent optimization. Among options of the auxiliary system and the associated linear response theory, reformulated versions of the particle–particle random phase approximation (pp-RPA) and the spin-flip time-dependent density functional theory (SF-TDDFT) are selected for illustration of principle. Numerical findings show that our multireference DFT successfully describes static correlation in bond dissociation and double bond rotation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electronic Origin and Kinetic Feasibility of the Lattice Oxygen Participation During the Oxygen Evolution Reaction on Perovskites

Density functional theory is employed to investigate the electronic origin and feasibility of surface lattice oxygen (O surf ) participation during the oxygen evolution reaction (OER) on perovskites. O surf participation occurs via the nonelectrochemical pathway in which adsorbed atomic oxygen (O*) diffuses from the transition-metal site to the oxygen site, and then O surf shifts out of the surface plane to react with O* to form O surf –O* and a surface oxygen vacancy. The different thermodynamic driving forces of O surf participation on LaMO 3-δ (M = Ni, Co, and Cu) are explained by the changes in the oxidation state of the transition-metal site throughout the reaction. We show that Osurf participation on LaNiO 3 cannot be hindered by Osurf protonation in the OER potential range. By including the coverage effect and utilizing the implicit solvent model, we finally show that lattice oxygen mechanism is more feasible than the conventional mechanism for OER on LaNiO 3 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Understanding the MXene Pseudocapacitance

MXenes have attracted great attention as next-generation capacitive energy-storage materials, but the mechanisms underlying their pseudocapacitive behavior are not well understood. Here we provide a theoretical description of the surface redox process of Ti3C2Tx (T = O, OH), a prototypical MXene, in 1 M H2SO4 electrolyte, based on joint density functional theory with an implicit solvation model and the analysis of Gibbs free energy under a constant-electrode potential. From the dependence of the O/OH ratio (or the surface H coverage) and the surface charge on the applied potential, we obtain a clear picture of the capacitive energy-storage mechanism of Ti3C2Tx that shows good agreement with previous experimental findings in terms of the integral capacitance and Ti valence change. We find a voltage-dependent redox/double-layer co-charging behavior: the capacitive mechanism is dominated by the redox process, but the electric double-layer charge works against the redox process. This new insight may be useful in improving the capacitance of MXenes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Stratification in Drying Films Containing Bidisperse Mixtures of Nanoparticles

Large scale molecular dynamics simulations for bidisperse nanoparticle suspensions with an explicit solvent are used to investigate the effects of evaporation rates and volume fractions on the nanoparticle distribution during drying. Our results show that “small-on-top” stratification can occur when Pe s Φ s ≳ c with c ~1, where Pe s is the Péclet number and Φ s is the volume fraction of the smaller particles. This threshold of Pe s Φ s for “small-on-top” is larger by a factor of ~α 2 than the prediction of the model treating solvent as an implicit viscous background, where α is the size ratio between the large and small particles. Our simulations further show that when the evaporation rate of the solvent is reduced, the “small-on-top” stratification can be enhanced, which is not predicted by existing theories. This unexpected behavior is explained with thermophoresis associated with a positive gradient of solvent density caused by evaporative cooling at the liquid/vapor interface. For ultrafast evaporation the gradient is large and drives the nanoparticles toward the liquid/vapor interface. This phoretic effect is stronger for larger nanoparticles, and consequently the “small-on-top” stratification becomes more distinct when the evaporation rate is slower (but not too slow such that a uniform distribution of nanoparticles in the drying film is produced), as thermophoresis that favors larger particles on the top is mitigated. A similar effect can lead to “large-on-top” stratification for Pe s Φ s above the threshold when Pe s is large but Φ s is small. Lastly, our results reveal the importance of including the solvent explicitly when modeling evaporation-induced particle separation and organization and point to the important role of density gradients brought about by ultrafast evaporation.

36 MATERIALS SCIENCE↗

Calculations of the pH-Dependent Onset Potential for CO Electrooxidation on Au(111)

CO electrooxidation on the Au(111) surface exhibits an onset potential that depends strongly on the pH of the electrolyte. In alkaline solution, the onset potential for CO electrooxidation is found, experimentally, to be 0.5 V lower than in acidic media on the reversible hydrogen electrode scale. This phenomenon is explained here with density functional theory which is used to calculate adsorbate binding energies including the electric double layer of the Au(111)/aqueous interface. Our model consists of a charged Au(111) slab and implicit solvation for the electrolyte. The double reference method is used to determine the potential-dependent CO electrooxidation reaction energetics. A microkinetic model, based upon the calculated reaction energetics, confi rms the lower onset potential in alkaline media. Our results show that there are three factors contributing to the lower onset potential in base: (1) stronger CO adsorption, (2) attraction between adsorbed CO and OH - , and (3) the high concentration of OH - in base.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nonlocal Effect of Percolated Particle Networks on Viscoelasticity of Polymer–Filler Nanocomposites: A Mesoscale Simulation Study

With nanoparticles (NPs) as fillers, polymer nanocomposites (PNCs) usually exhibit enhanced mechanical properties. However, a direct connection between the microscopic structural relaxation and macroscopic mechanical properties of PNCs remains to be established. To investigate the micro-to-macro connection, we develop a mesoscale model, in which the NPbridging polymer chains are represented by a dynamic bonded interaction between NPs, and the bulk polymer matrix is implicitly modeled by overdamped Langevin dynamics. Extensive equilibrium simulations are performed to quantify the microscopic dynamics of model PNCs. Systematic analyses of modified Rouse dynamics, dynamic structure factor, and relaxation modulus uncover that the microscopic relaxation dynamics of PNCs are significantly decelerated across different length scales because of nonlocal effects of percolated particle networks a phenomenon that has not been adequately captured in prior simulation studies. We find that NPvolume- fraction and NP-bonding-energy barrier are the two critical variables that affect bulk viscoelasticity the most. The proposed mesoscale model is versatile and provides a powerful framework for studying structure−property relations of different PNCs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Synergic Effects of Surface Chemistry and Applied Potentials on the Kinetics of the Electrocatalytic Oxygen Evolution Reaction in IrO 2

Improving the efficacy of the oxygen evolution reaction (OER) through water oxidation is critical for advancing photoelectrochemical water splitting. Among the catalysts, IrO 2 exhibits a high OER catalytic activity and stability under acidic conditions. The OER mechanism in this system has been a topic of intense research; however, many mechanistic understandings are lacking. In this work, we carried out first-principles calculations with an implicit solvation model at a constant potential to study the OER behavior on the IrO 2 (110) surface. We find that the surface hydrogen coverage has a significant effect on OER kinetics and transition states. We then develop a microkinetics model that accounts for the continuous evolution of both OER activation energy and hydrogen coverage as a function of an applied potential. We show that this inclusion leads to significant improvement in the simulated Tafel plot compared to available experiments. In conclusion, our results point to a complex interplay between surface chemistry and the applied potential on OER kinetics.

08 HYDROGEN↗

Quantum Mechanical Screening of Single-Atom Bimetallic Alloys for the Selective Reduction of CO 2 to C 1 Hydrocarbons

Electrocatalytic reduction of CO 2 to energy-rich hydrocarbons such as alkanes, alkenes, and alcohols is a very challenging task. So far, only copper has proven to be capable of such a conversion. Here, we report density functional theory (DFT) calculations combined with the Poisson-Boltzmann implicit solvation model to show that single-atom alloys (SAAs) are promising electrocatalysts for CO 2 reduction to C 1 hydrocarbons in aqueous solution. The majority component of the SAAs studied is either gold or silver, in combination with isolated single atoms, M (M = Cu, Ni, Pd, Pt, Co, Rh, and Ir), replacing surface atoms. We envision that the SAA behaves as a one-pot tandem catalyst: first gold (or silver) reduces CO 2 to CO, and the newly formed CO is then captured by M and is further reduced to C 1 hydrocarbons such as methane or methanol. We studied 28 SAAs, and found about half of them selectively favor the CO 2 reduction reaction over the competing hydrogen evolution reaction. Most of those promising SAAs contain M = Co, Rh, or Ir. The reaction mechanism of two SAAs, Rh@Au(100) and Rh@Ag(100), is explored in detail. Both of them reduce CO 2 to methane but via different pathways. For Rh@Au(100), reduction occurs through the pathway *CO → *CHO → *CHOH → *CH + H 2 O (l) → *CH 2 + H 2 O (l) → *CH 3 + H 2 O (l) → * + H 2 O (l) + CH 4(g) ; whereas, for Rh@Ag(100), the pathway is *CO → *CHO → *CH 2 O→ *OCH 3 → *O + CH 4(g) → *OH + CH 4(g) →* + H 2 O (l) + CH 4(g) . The minimum applied voltages to drive the two electrocatalytic systems are -1.01 and -1.12 V RHE for Rh@Au(100) and Rh@Ag(100), respectively, at which the Faradaic efficiencies for CO 2 reduction to CO are 60% for gold and 90% for silver. This suggests that SAA can efficiently reduce CO 2 to methane with as small as 40% loss to the hydrogen evolution reaction for Rh@Au(100) and as small as 10% for Rh@Ag(100). Lastly, we hope these computational results can stimulate experimental efforts to explore the use of SAA to catalyze CO 2 electrochemical reduction to hydrocarbons.

36 MATERIALS SCIENCE↗

Mechanism of CO 2 Reduction at Copper Surfaces: Pathways to C 2 Products

On the basis of constraints from reported experimental observations and density functional theory simulations, in this paper we propose a mechanism for the reduction of CO 2 to C 2 products on copper electrodes. To model the effects of an applied potential bias on the reactions, calculations are carried out with a variable, fractional number of electrons on the unit cell, which is optimized so that the Fermi level matches the actual chemical potential of electrons (i.e., the applied bias); an implicit electrolyte model allows for compensation of the surface charge so that neutrality is maintained in the overall simulation cell. Our mechanism explains the presence of the seven C 2 species that have been detected in the reaction, as well as other notable experimental observations. Furthermore, our results shed light on the difference in activities toward C 2 products between the (100) and (111) facets of copper. Finally, we compare our methodologies and findings with those in other recent mechanistic studies of the copper-catalyzed CO 2 reduction reaction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Challenges in Modeling Electrochemical Reaction Energetics with Polarizable Continuum Models

A major challenge in the modeling of electrochemical phenomena is the accurate description of the interface between an electrolyte and a charged conductor. Polarizable continuum models (PCM) have been gaining popularity because they offer a computationally inexpensive method of modeling the electrolyte. In this Perspective, we discuss challenges from using one such model which treats the ions using a linearized Poisson–Boltzmann (LPB) distribution. From a physical perspective, this model places charge unphysically close to the surface and adsorbates, and it includes excessively steep ramping of the dielectric constant from the surface to the bulk solvent. Both of these issues can be somewhat mitigated by adjusting parameters built into the model, but in doing so, the resultant capacitance deviates from experimental values. Likewise, hybrid explicit-implicit approaches to the solvent may offer a more realistic description of hydrogen bonding and solvation to reaction intermediates, but the corresponding capacitances also deviate from experimental values. These deviations highlight the need for a careful adjustment of parameters in order to reproduce not only solvation energies but also other physical properties of solid–liquid interfaces. Continuum approaches alone also necessarily do not capture local variations in the electric field from cations at the interface, which can affect the energetics of intermediates with substantial dipoles or polarizability. Finally, since the double-layer charge can be varied continuously, LPB/PCM models provide a way to determine electrochemical barriers at constant potential. However, double-layer charging and the atomic motion associated with reaction events occur on significantly different timescales. As a result, we suggest that more detailed approaches, such as the modified Poisson–Boltzmann model and/or the addition of a Stern layer, may be able to mitigate some but not all of the challenges discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Elucidating Mass Transport Regimes in Gas Diffusion Electrodes for CO 2 Electroreduction

Gas diffusion electrodes (GDEs) have shown promising performance for the electrochemical reduction of CO 2 (CO 2 R). In this study, a resolved, pore scale model of electrochemical reduction of CO 2 within a liquid-filled catalyst layer is developed. Three CO 2 mass transport regimes are identified in which the CO 2 penetration depth is controlled by CO 2 consumption in the electrolyte, CO 2 conversion along the solid-electrolyte double-phase boundaries (DPBs), and CO 2 conversion concentrated around the gas–solid–electrolyte triple-phase boundaries (TPBs). While it is possible for CO 2 R to be localized around the TPBs, in systems with submicron pore radii operating at <1 A cm –2 CO 2 R will be distributed across the DPBs within the catalyst layer. This validates the assumption of pore-scale uniformity implicit in popular, volume-averaged GDE models. Overall, the CO 2 conversion efficiency depends strongly on the governing mass transport regime, and operational-phase diagrams are constructed to guide the catalyst layer design.

10 SYNTHETIC FUELS↗

Influence of Charge Block Length on Conformation and Solution Behavior of Polyampholytes

In this paper, we investigate the effect of charge block length on polyampholyte chain conformation and phase behavior using small-angle X-ray scattering (SAXS) and implicit-solvent molecular simulations. To this end, we use solid phase peptide synthesis to precision-tailor a series of polyampholytes consisting of l-glutamic acid (E) and l-lysine (K) monomers arranged in alternating blocks from 2 to 16 monomers. We observe that the polyampholytes tend to phase separate as block size increases. With addition of NaCl, phase separated polyampholytes exhibit a salting-in effect dependent on charge block length. Fourier-transform infrared (FTIR) spectroscopy reveals the presence of intramolecular hydrogen bonds that are disrupted upon the addition of NaCl, implicating both electrostatic interactions and hydrogen bonding in the phase behavior. SAXS spectra at no-added salt conditions show minimal dependence of charge block length on the radius of gyration (R g ) for soluble polyampholytes, but local chain stiffening is found to be dependent on charge block length. With increasing NaCl, consistent with electrostatic screening, all polyampholytes expand and behave as neutral or swollen chains in good solvent conditions. Molecular simulations are qualitatively consistent with experiments. Implications for understanding intracellular condensates and material design are noted.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Thermal Stability of π-Conjugated n -Ethylene-Glycol-Terminated Quaterthiophene Oligomers: A Computational and Experimental Study

This work represents a joint computational and experimental study on a series of $\textit{n}$-ethylene glycol (PEO$\textit{n}$)-terminated quaterthiophene (4T) oligomers for 1 < $\textit{n}$ < 10 to elucidate their self-assembly behavior into a smectic-like lamellar phase. This study builds on an earlier study for $\textit{n}$ = 4 that showed that our model predictions were consistent with experimental data on the melting behavior and structure of the lamellar phase, with the latter consisting of crystal-like 4T domains and liquid-like PEO4 domains. The present study aims to understand how the length of the terminal PEO$\textit{n}$ chains modulates the disordering temperature of the lamellar phase and hence the relative stability of the ordered structure. In this work, a simplified bilayer model, where the 4T domains are not explicitly described, is put forward to efficiently estimate the disordering effect of the PEO domains with increasing n; this method is first validated by correctly predicting that layers of alkyl (PE)-capped 4T oligomers (for 1 < $\textit{n}$ < 10) stay ordered at room temperature. Both 4T-domain implicit and explicit model simulations reveal that the order-disorder temperature decreases with the length of the PEO capping chains, as the associated increase in conformational entropy drives a tendency toward disorder that overtakes the cohesive energy, keeping the ordered packing of the 4T domains.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Assessing the Effect of Explicit Polarizability on Models of Carbon Dioxide Solvation in Ionic Liquids

Ionic liquids are an important possible carbon capture material because of their anomalously high sorption selectivity for carbon dioxide over other gases common in air. Many research groups have investigated the molecular origins of this property and provided important insights, including using 1D and 2D-IR spectroscopy. Molecular dynamics simulations have been indispensable to the interpretation of these experiments. In prior molecular dynamics simulation work, charge-scaled force fields have typically been used to provide a mean-field treatment of effects vital to ionic liquid systems such as charge transfer and polarization. Here, we compare models of carbon dioxide solvated in ionic liquids with explicit polarization to models of the same with implicit polarizability through charge-scaling. We calculate structural, dynamical, and spectroscopic properties, and make comparisons to the same items measured in experiment. In this study, we focus on two ionic liquids: 1-butyl-3- methylimidazolium (BMIM + ) paired with bis(trifluoromethane sulfonyl imide) (Tf 2 N − ) and 1-butyl-3-methylimidazolium (BMIM+) paired with hexafluorophosphate (PF 6 − ). We find that many structural, dynamical, and spectroscopic properties are changed when polarization is modeled explicitly. We also find that explicit polarizability softens local ion cages around the carbon dioxide and that the long-time diffusion of the carbon dioxide is gated by the reorganization of the ionic liquid molecules. Comparisons to experiment show modest improvement of many observables compared with experiment for the explicitly polarizable model over the charge-scaled model. Overall, our results show that charge-scaled force fields are likely sufficient to compute spectroscopic properties of carbon dioxide in ionic liquids and suggest some interpretive rules for understanding their structural and dynamical properties. Those using charge-scaled force fields should generally assume that the ion cages around solutes such as carbon dioxide are too stiff and cation-rich in their models and adjust their interpretations and predictions accordingly.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Role of Reducible Oxide-Metal Cluster Charge Transfer in Catalytic Processes. New Insights on The Catalytic Mechanism of CO Oxidation on Au/TiO 2 from Ab Initio Molecular Dynamics

To probe metal particle/reducible oxide interactions Density Functional Theory based Ab Initio Molecular Dynamics studies were performed on a prototypical metal cluster (Au 2 0) supported on reducible oxides (rutile TiO 2 (110)) to implicitly account for finite temperature effects and the role of excess surface charge in the metal oxide. It was found that the charge state of the Au particle is negative in a reducing chemical environment whereas in the presence of oxidizing species co-adsorbed to the oxide surface the cluster obtained a net positive charge. In the context of the well known CO oxidation reaction, charge transfer facilitates the plasticization of Au 2 0 which allows for a strong adsorbate induced surface reconstruction upon addition of CO leading to the formation of mobile Au-CO species on the surface. The charging/discharging of the cluster during the catalytic cycle of CO oxidation enhances and controls the amount of O 2 adsorbed at oxide surface/cluster interface and strongly influences the energetics of all redox steps in catalytic conversions. A detailed comparison of the current findings with previous studies is presented and generalities about the role of surface-adsorbate charge transfer for metal cluster/reducible oxide interactions are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Refine your search to explore more results.