SEARCH · Engineering Papers

Results for “dynamical correlation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,621 records · Page 90

All sources NASA NTRS DOE OSTI

NASA NTRS1995

State Estimation of Main Rotor Flap and Lead-Lag Using Accelerometers and Laser Transducers on the RASCAL UH-60 Helicopter

Modern rotorcraft flight control system designs which promise to yield high vehicle response bandwidth and good gust rejection can benefit from the use of rotor-state feedbacks. The measurement of main rotor blade motions is also desirable to validate and improve rotorcraft simulation models, to identify high-order linear flight dynamics models, to provide rotor system health monitoring; during flight test, and to provide for correlation with acoustic measurements from wind tunnel and flight tests. However, few attempts have been made to instrument a flight vehicle in this manner, and no previous system has had the robustness and accuracy required for these diverse applications. A rotor blade motion measurement and estimation system has been developed by NASA and the U.S. Army for use on the Rotorcraft Aircrew Systems Concepts Airborne Laboratory (RASCAL) helicopter. RASCAL is a UH-60 Blackhawk which is being modified at Ames Research Center in a phased development program for use in flight dynamics and controls, navigation, airspace management, and rotorcraft human factors research. The aircraft will feature a full-authority, digital, fly-by-wire research flight control system; a coupled ring laser gyro, differential GPS based navigation system; a stereoscopic color wide field of view helmet, mounted display; programmable panel mounted displays; and advanced navigation sensors. The rotor blade motion system is currently installed for data acquisition only, but will be integrated with the research flight control system when it is installed later this year.

Fletcher, Jay W.↗

DOE OSTI2022

Applying Two-Dimensional Correlation Spectroscopy and Principal Component Analysis to Understand How Temperature Affects the Neptunium(V) Absorption Spectrum

The visible-near infrared (Vis-NIR) electronic absorption spectrum of neptunium(V) (NpO 2 + ) comprises numerous f-f electronic transitions with mostly undocumented temperature dependencies. The effect of temperature on the absorption spectrum of the pentavalent neptunyl dioxocation (NpO 2 + ) is an important factor to consider with spectrophotometric applications but has often been overlooked. Optical Vis-NIR absorption spectra (400–1700 nm) of Np(V) (0.017–0.89 M) in 1 M nitric acid were evaluated with varying temperatures (T = 10–80 °C). The intensity, position, and overall shape of the bands were sensitive to interactions with the solvent and coordination environment. Numerous temperature-induced isosbestic points were identified resulting from dynamic, overlapping peak shifts. Spectral variations were characterized using principal component analysis (PCA) and 2D correlation spectroscopy (COS). 2D-COS revealed that the absorption band near 1095 nm likely consists of two bands centered near 1087 and 1096 nm, which cannot be explained by current computational methods. 2D-COS analysis also provided an unambiguous assignment of unresolved peaks in the visible region for comparison with computational predictions. PCA was used to identify nonlinearity in the spectral response at elevated Np(V) concentrations ≥ 0.5 M. This unique experimental data and interpretation will foster a deeper understanding of the absorption spectra for complex actinyl ions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Sparse non-Markovian Noise Modeling of Transmon-Based Multi-Qubit Operations

The influence of noise on quantum dynamics is one of the main factors preventing current quantum processors from performing accurate quantum computations. Sufficient noise characterization and modeling can provide key insights into the effect of noise on quantum algorithms and inform the design of targeted error protection protocols. However, constructing effective noise models that are sparse in model parameters, yet predictive can be challenging. In this work, we present an approach for effective noise modeling of multi-qubit operations on transmon-based devices. Through a comprehensive characterization of seven devices offered by the IBM Quantum Platform, we show that the model can capture and predict a wide range of single- and two-qubit behaviors, including non-Markovian effects resulting from spatiotemporally correlated noise sources. The model’s predictive power is further highlighted through multi-qubit dynamical decoupling demonstrations and an implementation of the variational quantum eigensolver. As a training proxy for the hardware, we show that the model can predict expectation values within a relative error of 0.5%; this is a sevenfold improvement over default hardware noise models. Through these demonstrations, we highlight key error sources in superconducting qubits and illustrate the utility of reduced noise models for predicting hardware dynamics.

open quantum systems & decoherence↗

NASA NTRS1989

Digital Signal Processor For GPS Receivers

Three innovative components combined to produce all-digital signal processor with superior characteristics: outstanding accuracy, high-dynamics tracking, versatile integration times, lower loss-of-lock signal strengths, and infrequent cycle slips. Three components are digital chip advancer, digital carrier downconverter and code correlator, and digital tracking processor. All-digital signal processor intended for use in receivers of Global Positioning System (GPS) for geodesy, geodynamics, high-dynamics tracking, and ionospheric calibration.

Thomas, J. B.↗

DOE OSTI2019

Liquid-like and rigid-body motions in molecular-dynamics simulations of a crystalline protein

To gain insight into crystalline protein dynamics, we performed molecular-dynamics (MD) simulations of a periodic 2 × 2 × 2 supercell of staphylococcal nuclease. We used the resulting MD trajectories to simulate X-ray diffraction and to study collective motions. The agreement of simulated X-ray diffraction with the data is comparable to previous MD simulation studies. We studied collective motions by analyzing statistically the covariance of alpha-carbon position displacements. The covariance decreases exponentially with the distance between atoms, which is consistent with a liquidlike motions (LLM) model, in which the protein behaves like a soft material. To gain finer insight into the collective motions, we examined the covariance behavior within a protein molecule (intraprotein) and between different protein molecules (interprotein). The interprotein atom pairs, which dominate the overall statistics, exhibit LLM behavior; however, the intraprotein pairs exhibit behavior that is consistent with a superposition of LLM and rigid-body motions (RBM). Our results indicate that LLM behavior of global dynamics is present in MD simulations of a protein crystal. They also show that RBM behavior is detectable in the simulations but that it is subsumed by the LLM behavior. Finally, the results provide clues about how correlated motions of atom pairs both within and across proteins might manifest in diffraction data. Overall, our findings increase our understanding of the connection between molecular motions and diffraction data and therefore advance efforts to extract information about functionally important motions from crystallography experiments.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

Exploring origins for correlations between flow harmonics and transverse momentum in small collision systems

High statistics data sets from experiments at the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC) with small and large collision species have enabled a wealth of new flow measurements, including the event-by-event correlation between observables. One exciting such observable $ρ(v^2_n, [p_T])$ gauges the correlation between the mean transverse momentum of particles in an event and the various flow coefficients (v n ) in the same event [P. Bozek, Phys. Rev. C 93, 044908 (2016)]. Recently it has been proposed that very low multiplicity events may be sensitive to initial-state glasma correlations [G. Giacalone, B. Schenke, and C. Shen, Phys. Rev. Lett. 125, 192301 (2020)] rather than flow-related dynamics. We find utilizing the ip-jazma framework that the color domain explanation for the glasma results are incomplete. We then explore predictions from pythia8, and the version for including nuclear collisions called pythia-angantyr, which have only nonflow correlations, and from the ampt model which has both nonflow and flow-type correlations. We find that pythia-angantyr has nonflow contributions to $ρ(v^2_n, [p_T])$ in p + O, p + Pb, and O + O collisions that are positive at low multiplicity and comparable to the glasma correlations. It is striking that in pythia8 in p + p collisions there is actually a sign change from positive to negative $ρ(v^2_n, [p_T])$ as a function of multiplicity. The ampt results match the experimental data general trends in Pb+Pb collisions at the LHC, except at low multiplicity where ampt has the opposite sign. In p + Pb collisions, ampt has the opposite sign from experimental data and we explore this within the context of parton geometry. Predictions for p + O, O + O, and Xe + Xe are also presented.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Regulation of Side‐Chain Symmetry for Delaying Triplet Formation and Suppressing Non‐Radiative Loss in Organic Solar Cells

Abstract The structural revolutions of non‐fullerene acceptors (NFAs) have driven continuous efficiency breakthroughs in organic solar cells (OSCs). Rational regulation of NFA structures toward efficient exciton dissociation and mitigated non‐radiative recombination is pivotal for OSCs. The incorporation of asymmetric side chains on NFAs can often achieve these goals by inducing a desirable aggregate state. However, it lacks the studies to directly correlate the side‐chain symmetry of NFAs with the exciton delocalization and triplet dynamics in OSCs. Herein, The influence of structural symmetry on the aggregate properties is systematically investigated and exciton/charge dynamics based on two developed biaxial‐conjugated NFAs with varied side‐chain symmetry. The symmetric NFA having diverse molecular packing orientations can form multiple charge transfer channels in its blend with polymer donor, which cannot be found in that comprising the asymmetric ones. Moreover, a slower rate and lower ratio of the spin‐triplet state are formed in the blend of symmetric NFA, resulting in a much lower non‐radiative voltage loss in corresponding OSCs. This study reveals the distinct advantages of symmetric NFAs in both aggregate properties and exciton/charge dynamics over those of asymmetric ones, paving the way for developing high‐performance OSCs using easier‐to‐prepare, low‐cost symmetric materials.

Chemistry↗

DOE OSTI2020

A priori analysis on deep learning of subgrid-scale parameterizations for Kraichnan turbulence

In the present study, we investigate different data-driven parameterizations for large eddy simulation of two-dimensional turbulence in the a priori settings. These models utilize resolved flow field variables on the coarser grid to estimate the subgrid-scale stresses. We use data-driven closure models based on localized learning that employs a multilayer feedforward artificial neural network with point-to-point mapping and neighboring stencil data mapping, and convolutional neural network fed by data snapshots of the whole domain. The performance of these data-driven closure models is measured through a probability density function and is compared with the dynamic Smagorinksy model (DSM). The quantitative performance is evaluated using the cross-correlation coefficient between the true and predicted stresses. We analyze different frameworks in terms of the amount of training data, selection of input and output features, their characteristics in modeling with accuracy, and training and deployment computational time. We also demonstrate computational gain that can be achieved using the intelligent eddy viscosity model that learns eddy viscosity computed by the DSM instead of subgrid-scale stresses. Here, we detail the hyperparameters optimization of these models using the grid search algorithm.

42 ENGINEERING↗

DOE OSTI2021

Unraveling compacted graphite evolution during solidification of cast iron using in-situ synchrotron X-ray tomography

In spite of many years of research, the physical phenomena leading to the evolution of compacted graphite (CG) during solidification is still not fully understood. In particular, it is unknown how highly branched CG aggregates form and evolve in the semi-solid, and how local microstructural variations at micrometer length scale affect this growth process. We present here the first time-resolved synchrotron tomography combined with a bespoke high-temperature environmental cell that allows direct observation of the evolution of CG and relates this dynamic process to the local surrounding microstructures in a cast iron sample during repeated melting and solidification. Distinct processes are identified for the formation of CG involving the nucleation, growth, development of branches and interconnection of graphite particles, ultimately evolving into highly branched graphite aggregates with large sizes and low sphericities. CG is found to nucleate with a spheroidal or a plate-like shape, developing branches induced by high carbon concentration, e.g. thin melt channels. Additionally, CG grows much faster than spheroidal graphite during subsequent cooling in solid state. The direct visualization of the dynamic solidification process provides unprecedented new insights into formation mechanisms of CG and correlating factors such as local microstructural variations, and guides the development of CG iron solidification models.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Experimental evaluation of thermal-fluids performance of Helium-Cooled Flat Plate (HCFP) divertor

Various helium (He)-cooled solid-tungsten (W) divertor concepts have been proposed for long-pulse magnetic fusion energy reactors. Among these concepts, the He-cooled flat plate divertor (HCFP) modules have the largest plasma-facing surface area of ~0.2 m 2 . Simulations have shown that the most recent version of the design can withstand heat fluxes as great as 8 MW/m 2 . Earlier experimental studies of a single shortened HCFP cooling unit with a slot length of 7.6 cm used air at ambient temperature and pressures below 0.6 MPa. Here, we present initial experimental studies of a copper alloy and steel test section modeling a shortened HCFP cooling unit using He at prototypical pressure of 10 MPa, inlet temperatures T i ≤ 200 °C and steady-state incident heat fluxes q" ≤ 1.2 MW/m 2 . Results for Nusselt number Nu as a function of Reynolds number Re were obtained for Re = 1.2×10 4 –3.4×10 4 , and used to develop a Nu(Re) correlation and validate numerical models of the test section using commercial computational fluid dynamics (CFD) software. Simulations with this model are performed to evaluate the effect of the shortened slot. Furthermore, these analyses are used to estimate the thermal-fluids performance of the HCFP under prototypical conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Octahedral Distortion and Excitonic Behavior of Cs 3 Bi 2 Br 9 Halide Perovskite at Low Temperature

The metal halide ionic octahedron, represented as [MX 6 ] $n$- (M = metal cation, X = halide anion), serves as the basic structural unit in halide perovskites and plays a crucial role in determining their optoelectronic and chemical properties. Thus, it is possible to correlate the responses of metal halide perovskites to various environmental stimuli with the dynamic behaviors of the [MX 6 ] $n$- octahedra. In this study, with the temperature-dependent single-crystal X ray diffraction (SCXRD) measurements on Cs 3 Bi 2 Br 9 2D halide perovskites, we can identify two classes of distortions through the lowering of temperature: intraoctahedral distortion, which is the off-centering of Bi 3+ cation within a [BiBr 6 ] 3– octahedron due to the Bi 3+ 6s 2 lone pair electrons, and interoctahedral distortion, which is the collective misalignments among the [BiBr 6 ] 3– building blocks. Free exciton (FE) and self-trapped exciton (STE) models are used to study the relationship between the distortion of octahedra in Cs 3 Bi 2 Br 9 and the corresponding changes in its optoelectronic properties, which transform from dominating blue emission above 100 K to red emission at 4 K. In conclusion, this work provides new insights into the excitonic behaviors of perovskites and suggests a possibility that we can design and rationalize the optical properties of halide perovskites by regulating the environmental stimuli based on the knowledge of behaviors of the individual [MX 6 ] $n$- building blocks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Emergent Nanostructure and Ion Transport in Polyzwitterion/Polyanion Blends

Solid polymer electrolytes (SPEs) hold great promise for the advancement of next-generation energy storage devices. However, the ion transport mechanism in SPEs remains poorly understood. In this study, we investigate blends of poly(1-(3-sulfonatopropyl)-2-vinylpyridinium) (P2VPPS) and poly(lithium (trifluoromethane)sulfonimide methacrylate) (P(MTFSI)Li) of varying molar ratios to develop a mechanistic understanding of ionic conductivity in a miscible polyzwitterion/polyanion system. Polyanions can act as single-ion conductors, but conductivity is often prohibitively low due to the decreased segmental mobility and ion aggregation. Here, it is hypothesized that the introduction of a polyzwitterion would competitively interact with the polyanion charge groups to realize improvements in the conductivity. Attractive interactions between the polyanions and polyzwitterions are confirmed by the blend’s increased glass transition temperature using the Gordon–Taylor equation. Notably, an ordered local nanostructure (∼24 Å) emerged in the P2VPPS/P(MTFSI)Li system at certain compositions, as characterized by small-angle X-ray and neutron scattering (SAXS/SANS). Concurrent with the emergence of this structure, broadband dielectric spectroscopy confirmed improvements in ionic conductivity. The highest conductivity is observed at a specific blend ratio P2VPPS:P(MTFSI)Li = 0.2:1 in the glassy state and 0.3:1 in the rubbery state, corresponding to the lowest effective activation energy (E*). Coarse-grained molecular dynamics simulations further emphasize the role of complexation between polyzwitterion and polyanion chains, correlating with the emergence of a new peak in SAXS and SANS for the blends. This work provides a fresh perspective on the role of local structural design in developing SPEs and offers insights into the morphological effects on ionic conductivity.

chemical structure↗

DOE OSTI2020

Medium Energy Electron Flux in Earth's Outer Radiation Belt (MERLIN): A Machine Learning Model

The radiation belts of the Earth, filled with energetic electrons, comprise complex and dynamic systems that pose a significant threat to satellite operation. While various models of electron flux both for low and relativistic energies have been developed, the behavior of medium energy (120–600 keV) electrons, especially in the MEO region, remains poorly quantified. At these energies, electrons are driven by both convective and diffusive transport, and their prediction usually requires sophisticated 4D modeling codes. In this paper, we present an alternative approach using the Light Gradient Boosting (LightGBM) machine learning algorithm. The Medium Energy electRon fLux In Earth's outer radiatioN belt (MERLIN) model takes as input the satellite position, a combination of geomagnetic indices and solar wind parameters including the time history of velocity, and does not use persistence. MERLIN is trained on >15 years of the GPS electron flux data and tested on more than 1.5 years of measurements. Tenfold cross validation yields that the model predicts the MEO radiation environment well, both in terms of dynamics and amplitudes o f flux. Evaluation on the test set shows high correlation between the predicted and observed electron flux (0.8) and low values of absolute error. The MERLIN model can have wide space weather applications, providing information for the scientific community in the form of radiation belts reconstructions, as well as industry for satellite mission design, nowcast of the MEO environment, and surface charging analysis.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Development of an ERT‐Based Framework for Bentonite Buffers Monitoring From Laboratory Tests: 1. Characterizing Thermal–Hydrological–Mechanical Processes

Abstract Bentonite clay is widely used in engineered barrier systems for the permanent disposal of high‐level radioactive waste due to its low permeability, high swelling capacity, and thermal stability. However, the complex thermal‐hydrological‐mechanical (THM) processes induced by heating from decaying radioactive waste and hydration from surrounding rock can lead to heterogeneous changes that are difficult to measure and predict. This study develops an Electrical Resistivity Tomography (ERT)‐based framework for monitoring THM processes, progressing from sample‐scale to bench‐scale tests, to inform field‐scale applications. Sample‐scale tests analyzed small bentonite samples under controlled variations in water content, temperature, and porosity to establish fundamental resistivity relationships. Bench‐scale tests involved larger bentonite columns subjected to heating (up to 200°C) and hydration under controlled pressure, simulating repository conditions. ERT measurements, complemented by X‐ray CT imaging, temperature monitoring, and tracing sensors, revealed coupled THM processes, such as hydration‐induced compression, swelling, and thermal gradients, leading to complex resistivity patterns. The results demonstrate the potential of ERT for capturing THM‐induced resistivity changes, though challenges remain in upscaling and quantitative analysis. This study evaluates laboratory test capabilities and proposes future improvements for understanding THM‐induced resistivity responses. A conceptual framework for ERT implementation in field‐scale monitoring is presented, synthesizing findings from both scales and exploring how ERT data can inform long‐term modeling and reduce prediction uncertainties. Overall, this ERT‐based framework offers a robust method for monitoring bentonite buffers, aiding in early issue detection and supporting the safe long‐term disposal of radioactive waste in geological repositories, while highlighting the need for future development. Plain Language Summary Bentonite clay is crucial in engineered barrier systems (EBS) for containing high‐level radioactive waste due to its ability to absorb water, swell, seal and remain stable under high temperatures. When bentonite absorbs water and heats up from radioactive decay, it experiences complex changes in its physical and mechanical properties. Understanding these changes is important for ensuring the long‐term safety and effectiveness of EBS. This study used Electrical Resistivity Tomography (ERT), a non‐invasive method that measures electrical conductivity to monitor these changes during laboratory experiments. The ERT data revealed significant variations in resistivity corresponding to changes in water content, temperature, and density, providing detailed spatial and temporal insights into the behavior of bentonite. These findings enhance our ability to predict the long‐term performance of bentonite barriers, ensuring the safe containment of radioactive waste. By improving our understanding of bentonite's behavior, this research supports the development of more reliable and effective barrier systems for radioactive waste disposal, protecting the environment and public health. Key Points ERT monitoring was employed to capture resistivity changes in bentonite during controlled heating and hydration experiments, providing insights into THM processes ERT data reveal significant resistivity changes correlated with water content, temperature, and mechanical effects, enhancing the understanding of THM dynamics in bentonite This study explores the potential of the framework for application in field‐scale EBS monitoring, emphasizing the need for integrating additional geophysical methods for comprehensive subsurface imaging

Chen, Hang↗

DOE OSTI2021

Thermoresponsive polymer assemblies via variable temperature liquid-phase transmission electron microscopy and small angle X-ray scattering

Abstract Herein, phase transitions of a class of thermally-responsive polymers, namely a homopolymer, diblock, and triblock copolymer, were studied to gain mechanistic insight into nanoscale assembly dynamics via variable temperature liquid-cell transmission electron microscopy (VT-LCTEM) correlated with variable temperature small angle X-ray scattering (VT-SAXS). We study thermoresponsive poly(diethylene glycol methyl ether methacrylate) (PDEGMA)-based block copolymers and mitigate sample damage by screening electron flux and solvent conditions during LCTEM and by evaluating polymer survival via post-mortem matrix-assisted laser desorption/ionization imaging mass spectrometry (MALDI-IMS). Our multimodal approach, utilizing VT-LCTEM with MS validation and VT-SAXS, is generalizable across polymeric systems and can be used to directly image solvated nanoscale structures and thermally-induced transitions. Our strategy of correlating VT-SAXS with VT-LCTEM provided direct insight into transient nanoscale intermediates formed during the thermally-triggered morphological transformation of a PDEGMA-based triblock. Notably, we observed the temperature-triggered formation and slow relaxation of core-shell particles with complex microphase separation in the core by both VT-SAXS and VT-LCTEM.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Dynamic crystallography reveals spontaneous anisotropy in cubic GeTe

Cubic energy materials such as thermoelectrics or hybrid perovskite materials are often understood to be highly disordered. In GeTe and related IV–VI compounds, this is thought to provide the low thermal conductivities needed for thermoelectric applications. Since conventional crystallography cannot distinguish between static disorder and atomic motions, we develop the energy-resolved variable-shutter pair distribution function technique. This collects structural snapshots with varying exposure times, on timescales relevant for atomic motions. In disagreement with previous interpretations, we find the time-averaged structure of GeTe to be crystalline at all temperatures, but with anisotropic anharmonic dynamics at higher temperatures that resemble static disorder at fast shutter speeds, with correlated ferroelectric fluctuations along the <100> c direction. We show that this anisotropy naturally emerges from a Ginzburg–Landau model that couples polarization fluctuations through long-range elastic interactions. By accessing time-dependent atomic correlations in energy materials, we resolve the long-standing disagreement between local and average structure probes and show that spontaneous anisotropy is ubiquitous in cubic IV–VI materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2026

A curated benchmark for cofolding models on kinase conformational states

Abstract Protein kinases are critical drug targets, requiring therapeutics that can modulate their active and inactive conformational states. While cofolding models can generate global folds directly from kinase sequences and ligand SMILES strings, these models have not yet been tested on their ability to recover ligand-induced-fit conformational states of the kinase proteins. Here, we introduce KinConfBench, a curated benchmark of 2225 high-quality human kinase chains to evaluate the ability of four state-of-the-art cofolding models—Boltz-2, Chai-1, Protenix, and RoseTTAFold-All-Atom—to recover both canonical and rare conformational states. We show that geometric success metrics of a ligand pose in the active site do not correlate strongly with the correct kinase conformational state, motivating a new set of dynamical benchmarks for assessing cofolding models. While all four cofolding models achieve ~60–80% prediction accuracy for kinase conformational classification, they exhibit severe mode collapse when performing multiple inferences, show negligible structural diversity in sampling induced-fit motions, and display a prevalent “apo-drift” in which most cofolding models predominantly predict the kinase to be in its ligand-free state. Our results highlight that capturing ligand-induced protein conformational diversity, not just geometric fit, is critical for next-generation structure-based drug discovery.

Sun, Kunyang↗

DOE OSTI2025

Atomistic characterization of hydration-dependent fuel cell ionomer nanostructure: validation by vibrational spectroscopy

The development of Nafion alternatives for fuel cells and electrolyzers requires a fundamental understanding of hydration-dependent ion-exchange site acid/base chemistry. We present here reactive force field (ReaxFF) molecular dynamics (MD) simulations of Nafion at varying molar water/ion-exchange-site ratios (λ), which we correlate to our experimental and density functional theory-based vibrational spectra. ReaxFF describes the formation and breaking of covalent bonds, enabling simulations of proton exchange between sulfonic acid/sulfonate groups and water/hydronium. Our MD simulations determine the λ-dependent equilibrium proportions of protonated and deprotonated sites. We find that protonated sites persist across all λ from 0 to 20, challenging the widely accepted notion that all sites are ionized above a threshold λ value (e.g., 3 or 4). Our simulations generate hundreds of realistic exchange site environments, the characterization of which are based on 6 Å radii sulfur-centered ‘inner-spheres’. These inner- and outer-sphere regions elucidate an interplay of stereoelectronic factors that influence protonation states, including the number of inner-sphere waters (Λ). Our simulations produce broad λ-dependent distributions of Λ values, representing non-uniform exchange site hydration. Moreover, we demonstrate that these distributions are specific to both protonated and deprotonated sites, with significant overlap between the two distributions for all non-zero λ. These distributions underpin the IR spectra of hydrated membranes, with each exchange site contributing an IR spectrum characteristic of its protonation state and Λ. We expect that these nanostructural characterizations of Nafion exchange sites will contribute to the development of new ionomers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗