Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “polymer theory”

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 379 records · Page 21

Multi-fidelity learning for interatomic potentials: low-level forces and high-level energies are all you need

The promise of machine learning interatomic potentials (MLIPs) has led to an abundance of public quantum mechanical (QM) training datasets. The quality of an MLIP is directly limited by the accuracy of the energies and atomic forces in the training dataset. Unfortunately, most of these datasets are computed with relatively low-accuracy QM methods, e.g. density functional theory with a moderate basis set. Due to the increased computational cost of more accurate QM methods, e.g. coupled-cluster theory with a complete basis set (CBS) extrapolation, most high-accuracy datasets are much smaller and often do not contain atomic forces. The lack of high-accuracy atomic forces is quite troubling, as training with force data greatly improves the stability and quality of the MLIP compared to training to energy alone. Because most datasets are computed with a unique level of theory, traditional single-fidelity (SF) learning is not capable of leveraging the vast amounts of published QM data. In this study, we apply multi-fidelity learning (MFL) to train an MLIP to multiple QM datasets of different levels of accuracy, i.e. levels of fidelity. Specifically, we perform three test cases to demonstrate that MFL with both low-level forces and high-level energies yields an extremely accurate MLIP—far more accurate than a SF MLIP trained solely to high-level energies and almost as accurate as a SF MLIP trained directly to high-level energies and forces. Therefore, MFL greatly alleviates the need for generating large and expensive datasets containing high-accuracy atomic forces and allows for more effective training to existing high-accuracy energy-only datasets. Indeed, low-accuracy atomic forces and high-accuracy energies are all that are needed to achieve a high-accuracy MLIP with MFL.

36 MATERIALS SCIENCE↗

Real-Space Local Dynamics in 1,2,3-Triazole Using Inelastic Neutron Scattering

Enhancing proton transport in polymer electrolytes is crucial for advancing next-generation solid-state batteries, yet our understanding of proton conductivity in nonaqueous environments remains limited due to a lack of atomic-scale insights. Here, in this study, we investigated the atomic-scale dynamics of 1,2,3-triazole, a small molecule capable of dynamic hydrogen bonding, as a model system for proton hopping in nonaqueous environments. Using the real-space correlation function determined by the double Fourier transformation of inelastic neutron scattering spectra, we identified that the self-motion of protons and intermolecular dynamics occur on comparable time scales. Furthermore, we observed that the activation energy associated with the intermolecular dynamics matches the energy barrier for molecular rotations determined through Density Functional Theory calculations. These findings underscore the importance of controlling molecular dynamics at the atomic scale to control proton transport. Additionally, we demonstrated that intermolecular dynamics in systems involving protons can be studied using inelastic neutron scattering even without deuteration, thereby providing a broader avenue for studying atomic-scale dynamics in soft matter systems.

Shinohara, Yuya [Oak Ridge National Laboratory (OR↗

The Molecular Structures of Liquid and Glassy Nifedipine and Felodipine and Their Incorporation into PVP

Background: Amorphous drug formulations are commonly used to improve the solubility and bioavailability of poorly soluble molecular pharmaceuticals, yet less is known about their molecular conformations and local bonding interactions than their crystalline phases. Methods: High-energy X-ray diffraction structure factor measurements have been made on liquid and glassy nifedipine (NIF), felodipine (FEL), NIF 1:3 polyvinylpyrrolidone (PVP), and FEL 1:3 PVP wt.% mixtures. The corresponding X-ray pair distribution functions have been interpreted using empirical potential structure refinement using different models and density functional theory conformer calculations. Results: In both NIF and FEL, the NH···O inter-molecular hydrogen bonds between the pyridyl nitrogen and ester carbonyls are found to be considerably weaker than those observed in the crystalline polymorphs. For nifedipine, it is proposed that either inter-molecular NH…ON nitro bonds are present and/or a fraction (<20%) of conformational changes, with the aryl ring flipped, occur in the liquid state. For felodipine, the models indicate significant disorder associated with the methyl and ethyl side chains in the liquid state, with the main peak intensity at 3.0 Å arising from intra-molecular Cl-Cl atom pairs. When nifedipine molecules are incorporated into PVP, our models show they possess stronger NH···O bonds to the PVP polymer than felodipine molecules, which have stronger affinity for bonding to the polymer than to other felodipine molecules. Conclusions: The amorphous forms of both NIF and FEL show much weaker hydrogen bonding than found in their crystalline phases. Liquid NIF also exhibits configurations which are not observed in the crystal phases.

Amorphous Solid Dispersion↗

Dynamically morphing microchannels in liquid crystal elastomer coatings containing disclinations

Liquid crystal elastomers (LCEs) hold a major promise as a versatile material platform for smart soft coatings, since their orientational order can be predesigned to program a desired dynamic profile. In this work, we introduce temperature-responsive dynamic coatings based on LCEs with arrays of singular defects-disclinations that run parallel to the surface. The disclinations form in response to antagonistic patterns of the molecular orientation at the top and bottom surfaces, imposed by the plasmonic mask photoalignment. Upon heating, an initially flat LCE coating develops linear microchannels located above each disclination. The stimulus that causes a non-flat profile of LCE coatings upon heating is the activation force induced by the gradients of molecular orientation around disclinations. To describe the formation of microchannels and their thermal response, we adopt a Frank-Oseen model of disclinations in patterned director field and propose a linear elasticity theory to connect the complex spatially-varying molecular orientation to the displacements of the LCE. The thermo-responsive surface profiles predicted by the theory and by the finite element modeling are in good agreement with the experimental data; in particular, higher gradients of molecular orientation produce a stronger modulation of the coating profile. The elastic theory and the finite element simulations allow us to estimate the material parameter that characterizes the elastomer coating’s response to the thermal activation. The disclinations containing LCEs show potential as soft dynamic coatings with a predesigned responsive surface profile.

36 MATERIALS SCIENCE↗

(Invited) Continuum Mathematical Modeling of Water Electrolysis: A Tutorial

Widespread use of hydrogen energy is contingent on the development of reliable and economical sources of hydrogen. Electrolysis from renewably-derived low-carbon electricity is a potentially viable method of hydrogen generation. Prime among the electrolysis technologies are those utilizing ion-conducting polymers (ionomers) including proton-exchange-membrane water electrolyzer (PEMWE). However, these technologies need to exhibit increased efficiency, performance, and durability to become commercially viable. Like most electrochemical devices, PEMWEs involve multiple components (e.g., catalyst, ionomer, transport layers, membrane, plates) and multiple phases, with phenomena occurring across different time and length scales. Furthermore, it is difficult to experimentally probe many of the species and phenomena during operation. Thus, mathematical modeling at the continuum level has been an invaluable aid in exploring, understanding, and optimizing PEMWE cell and components. Furthermore, this is especially true in the highly coupled and complex physics and chemistries that occur with the membrane-electrode assembly (MEA). The physics in a typical volume-averaged non-isothermal model include multiphase transport in porous media, concentrated-solution theory, Ohm's Law, and Butler-Volmer kinetics, and ion, gas, and water transport in the ionomer.

Dizon, Arthur↗

Fluorine-Tuned Carbon-Based Nickel Single-Atom Catalysts for Scalable and Highly Efficient CO 2 Electrocatalytic Reduction

Electrocatalytic CO 2 reduction is garnering significant interest due to its potential applications in mitigating CO 2 and producing fuel. However, the scaling up of related catalysis is still hindered by several challenges, including the cost of the catalytic materials, low selectivity, small current densities to maintain desirable selectivity. In this study, Fluorine (F) atoms were introduced into an N-doped carbon-supported single nickel (Ni) atom catalyst via facile polymer-assisted pyrolysis. This method not only maintains the high atom utilization efficiency of Ni in a cost-effective and sustainable manner but also effectively manipulates the electronic structure of the active Ni-N 4 site through F doping. The catalyst has also been further optimized by controlling the F states, including convalent and semi-ionic states, by adjusting the fluorine sources involved. Consequently, this catalyst with unique structure exhibited comparable electrocatalytic performance for CO 2 -to-CO conversion, achieving a Faradaic efficiency (FE) of over 99% across a wide potential range and an exceptional CO evolution rate of 9.5 x 10 4 h -1 at -1.16 V vs reversible hydrogen electrode (RHE). It also delivered a practical current of 400 mA cm -2 while maintaining more than 95% CO FE. Experimental analysis combined with density functional theory (DFT) calculations have also shown that F-doping modifies the electron configuration at the central Ni-N 4 sites. In conclusion, this modification lowers the energy barrier for CO 2 activation, thereby facilitating the production of the crucial *COOH intermediate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Two-temperature warm dense hydrogen as a test of quantum protons driven by orbital-free density functional theory electronic forces

We consider a steady-state (but transient) situation in which a warm dense aggregate is a two-temperature system with equilibrium electrons at temperature Te, ions at Ti, and Te ≠ Ti. Such states are achievable by pump–probe experiments. For warm dense hydrogen in such a two-temperature situation, we investigate nuclear quantum effects (NQEs) on structure and thermodynamic properties, thereby delineating the limitations of ordinary ab initio molecular dynamics. We use path integral molecular dynamics (PIMD) simulations driven by orbital-free density functional theory (OFDFT) calculations with state-of-the-art noninteracting free-energy and exchange-correlation functionals for the explicit temperature dependence. We calibrate the OFDFT calculations against conventional (explicit orbitals) Kohn–Sham DFT. We find that when the ratio of the ionic thermal de Broglie wavelength to the mean interionic distance is larger than about 0.30, the ionic radial distribution function is meaningfully affected by the inclusion of NQEs. Moreover, NQEs induce a substantial increase in both the ionic and electronic pressures. This confirms the importance of NQEs for highly accurate equation-of-state data on highly driven hydrogen. For Te > 20 kK, increasing Te in the warm dense hydrogen has slight effects on the ionic radial distribution function and equation of state in the range of densities considered. In addition, we confirm that compared with thermostatted ring-polymer molecular dynamics, the primitive PIMD algorithm overestimates electronic pressures, a consequence of the overly localized ionic description from the primitive scheme.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Manning condensation in ion exchange membranes: A review on ion partitioning and diffusion models

Here, the rational design of ion exchange membranes (IEMs) is becoming more pertinent as their usage becomes broader and as their staple applications (i.e., electrodialysis, flow batteries, and fuel cells) improve in commercial viability. Such efforts would be catalyzed by an improved fundamental understanding of ion transport in IEMs. This review discusses recent progress in modeling ion partitioning and diffusion in IEMs in an effort to relate IEM performance metrics to fundamental membrane properties over which researchers and membrane manufacturers possess direct and sometimes precise control. Central focus is given to the Donnan-Manning model for ion partitioning and the Manning-Meares model for ion diffusion in IEMs. These two frameworks, which are derived from Manning's counter-ion condensation theory for polyelectrolyte solutions, have been widely used within the IEM literature since their recent introduction. To explore this topic, the mathematical derivation of both models is revisited, followed by a survey of experimental and computational discussions of counter-ion condensation in IEMs. Alternative models which fulfill similar roles in predicting IEM transport properties are compared. This review concludes by highlighting the uniquely favorable positions of the Donnan-Manning and Manning-Meares models and discussing their prospects as leading predictors of IEM partitioning and diffusive properties.

25 ENERGY STORAGE↗

CO 2 Capture Characteristics of Hyperbranched Poly(alkylene imine): A Molecular Dynamics Simulation Approach

This study explores the CO 2 capture characteristics of hyperbranched poly(ethylenimine) (HB-PEI) and poly- (propyleneimine) (HB-PPI) through molecular dynamics simulations using density functional theory-calibrated force fields. Key features such as density, free volume, glass transition temperature, CO 2 /H 2 O distribution, and molecular diffusion are systematically investigated to elucidate structure−function relationships under dry and hydrated conditions. HB-PEI demonstrates a slightly higher density and lower free volume compared to HB-PPI yet shows superior CO 2 capture due to the high amine concentration. Glass transition analysis indicates a higher thermal mobility in HBPEI, enhancing the CO 2 diffusivity. Pair correlation and coordination analyses confirm a stronger affinity of CO 2 with primary and secondary amines, particularly in hydrated environments where water competes with CO 2 for binding sites. Despite its more compact structure, HB-PEI outperformed HB-PPI in CO 2 and H 2 O transport, as confirmed by higher diffusion coefficients across all hydration levels. These findings highlight a critical balance among polymer architecture, amine accessibility, and hydration in designing next-generation solid amine sorbents for efficient direct air capture applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Use of Renewable Alcohols in Autocatalytic Production of Aspen Organosolv Lignins

This study aimed to investigate the intrinsic efficiency of renewable alcohols, applied under autocatalytic conditions, for removing lignin from aspen and hot-water-extracted aspen while substantially preserving the lignin structure so as to facilitate various valorization strategies. Ethylene glycol (EG), propylene glycol (PG), 1,4-butanediol (BDO), ethanol (EtOH), and tetrahydrofurfuryl alcohol (THFA) were evaluated based on their lignin solubilization ability, expressed as the relative energy difference (RED) following the principles of the Hansen solubility theory. The findings indicate that alcohols with a higher lignin solubilization potential lead to increased delignification, almost 90%, and produce a lignin with a higher content of β-O-4 bonds, up to 68% of those found in aspen milled wood lignin, thereby indicating their potential for valorization through depolymerization. However, these alcohols also produce lignin with a higher content of β-β and β-5 bonds, resulting in a higher molecular weight and polydispersity, due to readily occurring homolytic reactions. Hot-water extraction (HWE) conducted prior to alcohol treatment reduced the delignification efficiency and resulted in a lignin with a lower β-O-4 bond content. The lignins produced in these experiments exhibited a superior UV-A absorption capacity compared with synthetic benzophenone, as well as a greater radical quenching ability than synthetic butylated hydroxytoluene, indicating their potential for use in the protection of polymers against degradation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quantum Accurate Prediction of Plutonium–Plutonium Dihydride Phase Equilibrium Using a Lattice Gas Model

Plutonium-based materials are vital for use as nuclear fuels and as portable power sources for space vehicles. However, elucidating their sensitivity to hydriding corrosion represents an extreme challenge because of the toxicity of Pu as well as their anomalous magnetic properties. In this work, we develop a lattice gas model of plutonium–plutonium dihydride (Pu–PuH 2 ) phase equilibrium that retains the accuracy of density functional theory, while yielding many orders of magnitude improvement in computational efficiency. Using Monte Carlo and free-energy sampling algorithms, we compute a number of Pu–PuH 2 equilibrium properties that are difficult to probe experimentally, including equilibrium pressures and phase compositions and the PuH 2 heat of formation. In conclusion, our method will have particular impact on these types of materials studies, where there is a strong need for computationally efficient approaches to bridge time and length scale gaps between quantum calculations and experiments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Investigating Electron Conductivity Regimes in the Bacterial Cytochrome Wire OmcS

The anaerobic bacterium Geobacter sulfurreducens produces extracellular, electronically conductive cytochrome polymer wires that are conductive over micron length scales. Structure models from cryo-electron microscopy data show OmcS wires form a linear chain of hemes along the protein wire axis, which is proposed as the structural basis supporting their electronic properties. However, the mechanism by which this heme arrangement supports long-range electronic conduction remains unknown. Structure models from cryo-electron microscopy data show these wires form a linear chain of hemes along the protein wire axis, which is proposed as the structural basis supporting their electronic properties. Existing computational models using static heme redox potentials and coupling energies fail to explain experimental observations, predicting conductances 10,000 to 100,000 times lower than measured values. Here, we investigate how dynamic disorder affects site energies, interheme coupling, and long-range electronic conductivity within these cytochrome wires. We introduce an approach to extract charge carrier site information directly from Kohn–Sham density functional theory, without employing projector schemes, and show that site and coupling energies are highly sensitive to changes in interheme geometry and the surrounding electrostatic environment. Unlike models that incorporate dynamic disorder as a thermally averaged quantity, our quantum charge carrier model incorporates proxies for dynamic disorder through decoherence corrections, yielding predicted diffusion coefficient closer to what is expected from experiment and comparable with other organic-based electronic materials. Based on these simulations, we propose that the instantaneous fluctuations of the local electrostatic environment can transiently lift energy degeneracies and delocalize charge carriers. Furthermore, these studies reveal how incorporating dynamic fluctuations associated with the environment resolves the discrepancy between theory and experiment in microbial cytochrome wires and highlight design principles for bioinspired, heme-based conductive materials.

Bioinorganic chemistry↗

Electroresponse of weak polyelectrolyte brushes

Abstract End-tethered polyelectrolytes are widely used to modify substrate properties, particularly for lubrication or wetting. External stimuli, such as pH, salt concentration, or an electric field, can induce profound structural responses in weak polyelectrolyte brushes, which can be utilized to further tune substrate properties. We study the structure and electroresponsiveness of weak polyacid brushes using an inhomogeneous theory that incorporates both electrostatic and chain connectivity correlations at the Debye–Hückel level. Our calculation shows that a weak polyacid brush swells under the application of a negative applied potential, in agreement with recent experimental observation. We rationalize this behavior using a scaling argument that accounts for the effect of the surface charge. We also show that the swelling behavior has a direct influence on the differential capacitance, which can be modulated by the solvent quality, pH, and salt concentration. Graphical Abstract

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Revisiting Multi-Material Composite Structures with Homogenized Composite Properties

Composite structures inherently develop residual stresses during their curing process. Driven predominately by mismatched thermal strains between differing materials or ply orientations, but also affected by curing process phenomena like polymer shrinkage, these residual stresses can lead to failure within composite structures. There are several methods varying in complexity that can be used to model the development of residual stresses, all of which are capable of capturing sufficient detail to understand the residual stress state at the ply level. However, explicitly modeling all plies of a layup in a composite structure can be prohibitively expensive based on the number of plies, structure size, and required element size. The computational cost can be reduced through the homogenization of the composite layup without losing much fidelity of the overall response of the structure. The homogenization process reduces the many plies of a laminate to a single lamina that reduces complexity and increases the mesh size where a single element can span multiple plies. This report focuses on verification and validation efforts for a homogenization process using a suite of finite element simulations rather than an analytic solution derived from classical laminate theory. Initial verification using representative element volumes indicated there was minimal error in the homogenization process; however, this compounded to a small, but acceptable error in strip and split ring experimental composite structures. The error does under predict the residual stress state in the strip and split ring and should be accounted for when simulating composite structures with homogenized properties.

36 MATERIALS SCIENCE↗

The Fermi–Löwdin self-interaction correction for ionization energies of organic molecules

(Semi)-local density functional approximations (DFAs) suffer from self-interaction error (SIE). When the first ionization energy (IE) is computed as the negative of the highest occupied orbital (HO) eigenvalue, DFAs notoriously underestimate them compared to quasi-particle calculations. Here, the inaccuracy for the HO is attributed to SIE inherent in DFAs. We assessed the IE based on Perdew-Zunger self-interaction correction on 14 small to moderate-sized organic molecules relevant in organic electronics and polymer donor materials. Though self-interaction corrected DFAs were found to significantly improve the IE relative to the uncorrected DFAs, they overestimate. However, when the self-interaction correction is interiorly scaled using a function of the iso-orbital indicator zs , only the regions where SIE is significant get a correction. We discuss these approaches and show how these methods significantly improve the description of the HO eigenvalue for the organic molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Star Block Copolymers at Homopolymer Interfaces: Conformation and Compatibilization

Star block copolymers (s-BCPs), comprised of multiple linear diblock copolymers joined at a central point, are shown to segregate to the interface between two immiscible homopolymers that are identical to the blocks of the s-BCPs. The s-BCPs undergo a configurational transition at the interface, with different blocks of copolymers being embedded in their respective homopolymers, thereby bridging the interface and promoting adhesion. A series of 4-arm s-BCPs were synthesized with hydrogenated or deuterated polystyrene (PS/dPS) as the core block and poly(2-vinylpyridine) (P2VP) as the corona block, which was directly placed at the interface between the two homopolymers. Neutron reflectivity (NR) was used to determine the concentration profiles of the PS homopolymer, s-BCP core blocks, and P2VP total segments under equilibrium. The investigation varies the molecular weight (MW) and the total number of s-BCPs at the interface. Self-consistent-field theory (SCFT) was also employed to calculate the concentration profiles of the components at the interface, which were in excellent agreement with experimental results. The NR showed that the interfacial width between the homopolymers increased with the increasing number of s-BCPs at the interface up to a saturation limit. Beyond this limit, additional s-BCPs were released into the corona-miscible phase as unimolecular micelles. For a comparable interlayer thickness of s-BCPs at the interface, lower MW s-BCPs generated a broader interface. SCFT analysis suggested that, at the same packing density, the arms of the low MW s-BCPs align more parallel to the interface, while the arms of high MW s-BCPs adopt a more normal orientation, like their linear BCP counterparts. Furthermore, it was also observed that the core blocks, constrained by the junction points, were oriented more parallel and closer to the interface than the corona blocks. The phase behavior of the polymer blends revealed that s-BCP additives can efficiently reduce the domain size, with the low MW yielding smaller domain sizes due to the greater reduction in the interfacial energy and the high MW arresting phase separation due to their higher binding energy and a jamming of the interfacial assemblies. Asymmetric double cantilever beam (ADCB) tests demonstrated that s-BCPs promoted adhesion more efficiently than their linear BCP counterparts due to stronger binding energy per molecule, suggesting a more efficient compatibilizer for polymer upcycling. In conclusion, the results from these studies provide fundamental insights into the assembly of s-BCPs at homopolymer interfaces, the reduction of domain size, and promotion of adhesion, providing a strategy for the use of s-BCPs as stealth surfactants and universal compatibilizers.

interface↗

Modeling of surface mechanical behaviors of soft elastic solids: theory and examples

Surfaces of soft solids can have significant surface stress, extensional modulus and bending stiffness. Previous theoretical studies have usually examined cases in which both the surface stress and bending stiffness are constant, assuming small deformation. In this work we consider a general formulation in which the surface can support large deformation and carry both surface stresses and surface bending moments. We show that the large deformation theory can be reduced to the classical linear theory (Shuttleworth equation). We obtain exact solutions for problems of an inflated cylindrical shell and bending of a plate with a finite thickness. Our analysis illustrates the different manners in which surface stiffening and surface bending stabilize these structures. We discuss how the complex surface constitutive behaviors affect the stress field of the bulk. Our calculation provides insights into effects of strain-dependent surface stress and surface bending in the large deformation regime, and can be used as a model to implement surface finite elements to study large deformation of complex structures.

36 MATERIALS SCIENCE↗

Modeling of crystallization fouling on a horizontal-tube falling-film evaporator for thermal desalination

Scale formation within horizontal-tube falling-film evaporators is an important issue for thermal desalination, because of the deterioration of heat transfer performance and added maintenance costs. Predicting the fouling process is crucial for desalination plant design, operation, and maintenance. In this paper, a model is developed to predict the crystallization fouling process during seawater falling-film flow over a horizontal tube bundle. The scaling is estimated based on deposition theory, with couple heat transfer including the in-tube steam condensation, conduction through the tube wall and the scale layer, and falling-film evaporation. The spatial and temporal variations of temperature, heat transfer coefficient, and scale thickness are predicted, and the effects of variations in process parameters of steam, feed seawater, and tube properties on the scale thickness, evaporation rate, and heat transfer coefficient for falling-film evaporation are analyzed. Comparisons to existing experimental data show the scaling layer thickness falls in the general range of reported measurements, with a slight overestimation due to neglecting the delay of scaling onset under real conditions. The scale layer thickness increases dramatically as the in-tube steam pressure is increased or the seawater flow rate is decreased. In general, the scaling layer becomes thicker on the lower tubes in the tube bundle, due to the heat and mass transfer of the falling film. The effect of tube material on scaling appears mainly dependent on thermal conductivity; thus, using polymer rather than a stainless steel tube decreases scaling and evaporation rate 96% and 88%, respectively. Finally, this work has potential to guide thermal desalination plant design, material selection, and operating parameter optimization.

42 ENGINEERING↗