Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Bounds”

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,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Unimolecular and bimolecular formic acid decomposition routes on dispersed Cu nanoparticles

The elementary steps and site requirements in formic acid (HCOOH) dehydrogenation on Cu surfaces remain of keen interest because formate species act as intermediates or spectators in methanol synthesis and water–gas shift reactions. Steady-state and transient kinetic data, isotopic effects, infrared spectra during catalytic and stoichiometric reactions, and theoretical treatments based on density functional theory (DFT) provide evidence for bimolecular reactions, in which saturated bidentate formate (*HCOO*) adlayers, present at 0.25 ML (0.25 *HCOO* per surface Cu atom), react with undissociated species (HCOOH$^\square$) bound at interstices within formate adlayers ($^\square$) to form H-bonded bimolecular HCOOH$^\square$-*HCOO* adducts. The co-existence of vicinal HCOOH$^\square$ and *HCOO* moieties is evident from antisymmetric infrared bands for *HCOO* that become stronger as a result of their H-bonding that perturbs the induced dipole moment of *HCOO* upon vibration, consistent with DFT-derived vibrational frequencies and intensities for such perturbed species. The *HCOO* moiety in this complex undergoes C-H activation via a transition state that is preferentially stabilized through H-bonding with the vicinal HCOOH$^\square$ relative to its *HCOO* precursor. DFT-derived HCOOH dehydrogenation activation barriers and those determined from the evolution of CO 2 from pre-adsorbed *HCOO* species are about 10 kJ mol -1 smaller in the presence of gaseous HCOOH reactants (because of HCOOH$^\square$-*HCOO* interactions) than those for the unimolecular decomposition of bound *HCOO* species. Such bimolecular routes are consistent with measured effects of HCOOH, H 2 , and CO pressures and of H/D isotopic substitution on dehydrogenation turnover rates and represent the predominant channel for the formation of CO 2 and H 2 during catalytic HCOOH dehydrogenation on Cu nanoparticles. A saturated *HCOO* adlayer that retains binding interstices and the presence of HCOOH(g) enable a sequence of elementary steps unavailable for *HCOO* species, thus circumventing unassisted unimolecular routes that exhibit higher activation barriers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Intersection-distribution-based remapping between arbitrary meshes for staggered multi-material arbitrary Lagrangian-Eulerian hydrodynamics

In this paper, we present a new intersection-distribution-based remapping method between arbitrary polygonal meshes for indirect staggered multi-material arbitrary Lagrangian-Eulerian hydrodynamics. All cell-centered material quantities are conservatively remapped using intersections between the Lagrangian (old, source) mesh and the rezoned (new, target) mesh. The new nodal masses are obtained by conservative distribution of all material masses in each new cell to the cell's corners and then collecting those corner masses at new nodes. This distribution is done using a local constrained optimization approach for each cell in the new mesh. In order to remap nodal momentum we first define cell-centered momentum for each cell in the old mesh, conservatively remap this to the new mesh and then conservatively distribute the new zonal momentum to each cell's bounding nodes, again using local constrained optimization. Our method also conserves total energy by applying a new nodal kinetic energy correction that relies on a process similar to that used for remapping nodal mass and momentum. Cell-centered kinetic energy is computed, conservatively remapped and then distributed to nodes. The discrepancy between this conservatively remapped and actual nodal kinetic energy is then conservatively distributed to the internal energies of the materials in the cells surrounding each node. Unlike conventional cell-based corrections of this type, this new nodal kinetic energy correction has not been observed to drive material internal energy negative in any of our testing. Unlike flux based remapping, our new intersection-distribution method can be applied to remapping between source and target meshes that are arbitrarily different, which provides superior flexibility in the rezoning strategy. Our method is accurate, essentially conservative and essentially bounds preserving.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A positivity-preserving and conservative intersection-distribution-based remapping algorithm for staggered ALE hydrodynamics on arbitrary meshes

In this study, we introduce new intersection-distribution-based remapping tools for indirect staggered arbitrary Lagrangian-Eulerian (ALE) simulations of multi-material shock hydrodynamics on arbitrary meshes. In addition to conserving momentum and total energy, the three-stage remapper proposed in this work preserves non-negativity of the internal energy. At the first stage, we construct slope-limited piecewise-linear reconstructions of all conserved quantities on zones of the source mesh and perform intersection-based remap to obtain bound-preserving zonal quantities on the target mesh. At the second stage, we define bound-preserving nodal quantities of the staggered ALE discretization as convex combinations of corner quantities. The nodal internal energy is corrected in a way which keeps it non-negative, while providing exact conservation of total energy. At the final stage, we distribute the non-negative nodal internal energy to corners, zones and materials using non-negative weights. Proofs of positivity preservation are provided for each stage. This work is a natural extension of our paper [14] in which a similar intersection-distribution-based remapping procedure was employed. The original version used a nodal kinetic energy fix which did not provably ensure positivity preservation for the zonal internal energy after the final distribution stage. The new algorithm cures this potential drawback by using ‘coordinated’ limiters for piecewise-linear reconstructions, remapping the internal energy to nodes and correcting it before redistribution. The effectiveness of the new nodal fix is illustrated by numerical examples.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Higher-order particle representation for particle-in-cell simulations

In this paper we present an alternative approach to the representation of simulation particles for unstructured electrostatic and electromagnetic PIC simulations. In our modified PIC algorithm we represent particles as having a smooth shape function limited by some specified finite radius, r 0 . A unique feature of our approach is the representation of this shape by surrounding simulation particles with a set of virtual particles with delta shape, with fixed offsets and weights derived from Gaussian quadrature rules and the value of r 0 . As the virtual particles are purely computational, they provide the additional benefit of increasing the arithmetic intensity of traditionally memory bound particle kernels. The modified algorithm is implemented within Sandia National Laboratories' unstructured EMPIRE-PIC code, for electrostatic and electromagnetic simulations, using periodic boundary conditions. We show results for a representative set of benchmark problems, including electron orbit, a transverse electromagnetic wave propagating through a plasma, numerical heating, and a plasma slab expansion. In this work, good error reduction across all of the chosen problems is achieved as the particles are made progressively smoother, with the optimal particle radius appearing to be problem-dependent.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Parallel transport dynamics for mixed quantum states with applications to time-dependent density functional theory

Direct simulation of the von Neumann dynamics for a general (pure or mixed) quantum state can often be expensive. One prominent example is the real-time time-dependent density functional theory (rt-TDDFT), a widely used framework for the first principle description of many-electron dynamics in chemical and materials systems. Practical rt-TDDFT calculations often avoid the direct simulation of the von Neumann equation, and solve instead a set of Schrödinger equations, of which the dynamics is equivalent to that of the von Neumann equation. However, the time step size employed by the Schrödinger dynamics is often much smaller. Here, in order to improve the time step size and the overall efficiency of the simulation, we generalize a recent work of the parallel transport (PT) dynamics for simulating pure states [An, Lin, Multiscale Model. Simul. 18, 612, 2020] to general quantum states. The PT dynamics provides the optimal gauge choice, and can employ a time step size comparable to that of the von Neumann dynamics. Going beyond the linear and near adiabatic regime in previous studies, we find that the error of the PT dynamics can be bounded by certain commutators between Hamiltonians, density matrices, and their derived quantities. Such a commutator structure is not present in the Schrödinger dynamics. We demonstrate that the parallel transport-implicit midpoint (PT-IM) method is a suitable method for simulating the PT dynamics, especially when the spectral radius of the Hamiltonian is large. The commutator structure of the error bound, and numerical results for model rt-TDDFT calculations in both linear and nonlinear regimes, confirm the advantage of the PT dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

A hybrid Monte Carlo-deterministic second moment method with efficient variance reduction

In this work, we present a hybrid method that combines Monte Carlo with deterministic finite element methods to solve a linear Boltzmann transport equation. Our hybrid method runs orders of magnitude faster than Monte Carlo, without sacrificing accuracy, for a proxy problem from radiative transfer that contains both optically-thick and optically-thin material. We believe that this is the first demonstration of a hybrid Second Moment Method in more than one spatial dimension, the first to consider more than one material, and the first to use variance reduction. Our variance reduction approach arises from an asymptotic analysis in which we show that the magnitude of the scattering source grows without bound. We transform the problem to compute the deviation of the radiation intensity from isotropy. The magnitude of the source in the transformed problem is bounded, and the quality of the hybrid method solution is dramatically improved by a substantial reduction in the variance.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Reprint of "Selective conversion of CO into ethanol on Cu(511) surface reconstructed from Cu(pc): Operando studies by electrochemical scanning tunneling microscopy, mass spectrometry, quartz crystal nanobalance, and infrared spectroscopy"

A polycrystalline copper, surface-terminated by a well-defined (511)-oriented facet, was electrochemically generated by a series of step-wise surface reconstruction and iterations of mild oxidative-reductive processes in 0.1 M KOH. The electrochemical reduction of CO on the resultant stepped surface was investigated by four surface-sensitive operando methodologies: electrochemical scanning tunneling microscopy (STM), electrochemical quartz crystal nanobalance (EQCN), differential electrochemical mass spectrometry (DEMS), and polarization-modulation infrared spectroscopy (PMIRS). The stepped surface catalyzed the facile conversion of CO into ethanol, the exclusive alcohol product at a low overpotential of -1.06 V (SHE) or - 0.3 V (RHE). Additionally, the chemisorption of CO was found to be a necessary prelude to ethanol production; i.e. the surface coverages, rather than solution concentrations, of CO and its surface-bound intermediates primarily dictate the reaction rates (current densities). Contrary to the expected predominance of undercoordinated step-site reactivity over the coordination chemistry of vicinal surfaces, vibrational spectroscopic evidence reveals the involvement of terrace-bound CO adsorbates during the multi-atomic transformations associated with the production of ethanol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Re-evaluating probable maximum precipitation estimates: sensitivity to transposition domains and storm rotation using modern datasets

This study examines the sensitivity of Probable Maximum Precipitation (PMP) estimates to key methodological decisions embedded in the legacy approach adopted in the U.S. National Weather Service Hydrometeorological Reports No. 51 and No. 52. Although widely used for infrastructure design and risk regulation, fundamental aspects of PMP estimation—such as storm sample size, transposition domain, maximization procedures, and storm rotation—remain poorly constrained and lack formal guidance. Using the Red Rock watershed in Iowa as a case study, and leveraging the 2002–2023 NOAA Analysis of Record for Calibration (AORC) precipitation dataset, we systematically evaluate how each methodological choice, individually and in combination, influences PMP estimates. Our findings demonstrate that PMP is not a fixed physical upper bound but rather a modeling construct shaped heavily by user-defined assumptions. Notably, PMP values derived from modern gridded rainfall datasets can be substantially higher than the legacy estimate used in the original spillway design for Red Rock Dam. Decisions regarding storm sample size, domain extent, climatological window, and particularly storm rotation all contributed to higher PMP estimates. Storm rotation alone—a loosely constrained element in the current PMP practice—can amplify PMP by more than 25%. These results reveal the lack of standardized bounds in current PMP workflows and the need for systematic sensitivity and uncertainty analysis. As PMP estimation shifts toward probabilistic approaches, incorporating physically meaningful storm attributes will be key to developing more transparent, defensible methods for dam safety and climate-resilient infrastructure.

Probable maximum precipitation↗
DOE OSTI2022

S K-edge XAS of Cu Ⅱ , Cu Ⅰ , and Zn Ⅱ oxidized Dithiolene complexes: Covalent contributions to structure and the Jahn-Teller effect

We report reduced dithiolene ligands are bound to high valent Mo centers in the active site of the oxotransferase family of enzymes. Related model complexes have been studied with great insight by Prof. Holm and his colleagues. This study focuses on the other limit of dithiolene chemistry: an investigation of the 2-electron oxidized dithiolene bound to low-valent late transition metal (TM) ions (Zn II , Cu I , and Cu II ). The bonding descriptions of the oxidized dithiolene [N,N-dimethyl piperazine 2,3-dithione (Me 2 Dt 0 )] complexes are probed using S K-edge X-ray absorption spectroscopy (XAS) and the results are correlated to density functional theory (DFT) calculations. These experimentally supported calculations are then extended to explain the different geometric structures of the three complexes. The Zn II (Me 2 Dt 0 ) 2 complex has only ligand-ligand repulsion so it is stabilized at the D 2d symmetry limit. The Cu I (Me 2 Dt 0 ) 2 complex has additional weak backbonding thus distorts somewhat from D 2d toward D 2h symmetry. The Cu II (Me 2 Dt 0 ) 2 complex has a strong σ donor bond that leads to both a large Jahn-Teller stabilization to D 2h and an additional covalent contribution to the geometry. The combined strong stabilization results in the square planar, D 2h structure. This study quantifies the competition between the ligand-ligand repulsion and the change in electronic structures in determining the final geometric structures of the oxidized dithiolene complexes, and provides quantitative insights into the Jahn-Teller stabilization energy and its origin.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Carbon monoxide chemistry of α-V70I Mo-nitrogenase: Evidence from EPR- and IR-monitored photolysis – or, what a difference a methyl makes

A critical step in the global nitrogen cycle is the conversion of dinitrogen into biologically accessible ammonia. In Nature this is accomplished by the nitrogenase (N 2 ase) family of enzymes. Carbon monoxide (CO) has long been known as an inhibitor of dinitrogen reduction by N 2 ase, but it can also be a substrate of the enzyme, when it is catalytically reduced to hydrocarbons. Understanding the CO interactions with N 2 ases are thus relevant to both dinitrogen fixation and Fischer-Tropsch-like chemistry. Here, in this work, the interaction of CO with the α-V70I variant of Azotobacter vinelandii MoFe N 2 ase was investigated using electron paramagnetic resonance (EPR) and infrared (IR) monitored photolysis of bound CO under cryogenic conditions. This was supplemented by further analysis of stopped-flow Fourier transform IR (SF-FT-IR) data under turnover conditions. The α-V70I variant adds a single methyl group close to the FeMo-cofactor active site, and the results show that this inhibits and slows, but does not substantially chemically change, the binding of CO to the FeMocofactor. The EPR spectra of both the hi-CO and lo-CO states closely resemble those from the wild-type enzyme. Similarly, the SF-FT-IR spectra of CO inhibited α-V70I and wild-type enzyme are strikingly similar, showing only small shifts in band energies which allow better interpretation of the published wild-type spectra. The extra carbon does, however, impact and inhibit the photochemical release and migration of CO at cryogenic temperatures, resulting in novel CO-bound species. These include a product species, termed Lo-1*, which may involve CO photochemically migrating on the FeMo-cofactor.

Carbon monoxide↗
DOE OSTI2025

Efficient screening of rare large pit anomalies on polished surfaces using a minimalist sampling scheme

Lawrence Livermore National Laboratory (LLNL) has made significant strides in generating clean energy through its inertial confinement fusion (ICF) experiments. These experiments rely on high-density carbon (HDC) coated shells to encapsulate the fusion fuel. The success of these experiments is heavily dependent on the surface quality of these shells, as even minor imperfections, such as deep pits, can negatively impact fusion yield. Ensuring the required smoothness involves an extensive surface-finishing process that spans approximately 20 stages, making it both time-intensive and resource-demanding. A critical challenge in this process is the need for high-resolution scans to detect rare deep pits, which can be costly and impractical if performed on every shell. This highlights the necessity of developing more efficient scanning methods to optimize time and cost without compromising accuracy. To address these challenges, we introduce a novel approach that employs the multivariate Dvoretzky–Kiefer–Wolfowitz (DKW) inequality to provide a probabilistic upper bound on the error in estimating pit distribution characteristics via a Kernel Density Estimator (KDE). This error bound enables efficient and reliable estimation of pit distribution characteristics at a specified statistical confidence level using a minimal number of surface scans. The integrated DKW-KDE approach was validated through surface-finishing experiments across two batches of HDC-coated shells, demonstrating consistent and robust performance across multiple stages of the surface-finishing experiments. The validation studies suggest that the integrated DKW-KDE approach achieves comparable accuracy in estimating the risk of deleterious large pits with six scans, thus conserving time and resources. Further evaluations show that performance remains consistent across batches and over multiple polishing stages. In conclusion, based on these findings, one can leverage the minimal-scan insights to strategically improve the bottleneck inspection process, thus enhancing the productivity and quality of shell polishing and similar challenging manufacturing processes.

Inertial confinement fusion↗
DOE OSTI2020

Calcitonin Receptor N-Glycosylation Enhances Peptide Hormone Affinity by Controlling Receptor Dynamics

The class B G protein-coupled receptor (GPCR) calcitonin receptor (CTR) is a drug target for osteoporosis and diabetes. N-glycosylation of asparagine 130 in its extracellular domain (ECD) enhances calcitonin hormone affinity with the proximal GlcNAc residue mediating this effect through an unknown mechanism. Here, we present two crystal structures of salmon calcitonin-bound, GlcNAc-bearing CTR ECD at 1.78 and 2.85 Å resolutions and analyze the mechanism of the glycan effect. The N130 GlcNAc does not contact the hormone. Surprisingly, the structures are nearly identical to a structure of hormone-bound, N-glycan-free ECD, which suggested that the GlcNAc might affect CTR dynamics not observed in the static crystallographic snapshots. Hydrogen-deuterium exchange mass spectrometry and molecular dynamics simulations revealed that glycosylation stabilized a β-sheet adjacent to the N130 GlcNAc and the N-terminal α-helix near the peptide-binding site while increasing flexibility of the peptide-binding site turret loop. These changes due to N-glycosylation increased the ligand on-rate and decreased its off-rate. The glycan effect extended to RAMP-CTR amylin receptor complexes and was also conserved in the related CGRP receptor. These results reveal that N-glycosylation can modulate GPCR function by altering receptor dynamics.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Conformational Changes of RORγ During Response Element Recognition and Coregulator Engagement

The retinoic acid receptor-related orphan receptor γ (RORγ) is a ligand-dependent transcription factor of the nuclear receptor super family that underpins metabolic activity, immune function, and cancer progression. Despite being a valuable drug target in health and disease, our understanding of the ligand-dependent activities of RORγ is far from complete. Like most nuclear receptors, RORγ must recruit coregulatory protein to enact the RORγ target gene program. To date, a majority of structural studies have been focused exclusively on the RORγ ligand-binding domain and the ligand-dependent recruitment of small peptide segments of coregulators. Herein, we examine the ligand-dependent assembly of full length RORγ:coregulator complexes on cognate DNA response elements using structural proteomics and small angle x-ray scattering. The results from our studies suggest that RORγ becomes elongated upon DNA recognition, preventing long range interdomain crosstalk. We also determined that the DNA binding domain adopts a sequence-specific conformation, and that coregulatory protein may be able to ‘sense’ the ligand- and DNA-bound status of RORγ. We propose a model where ligand-dependent coregulator recruitment may be influenced by the sequence of the DNA to which RORγ is bound. Overall, the efforts described herein will illuminate important aspects of full length RORγ and monomeric orphan nuclear receptor target gene regulation through DNA-dependent conformational changes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

A Model for the Signal Initiation Complex Between Arrestin-3 and the Src Family Kinase Fgr

Arrestins regulate a wide range of signaling events, most notably when bound to active G protein-coupled receptors (GPCRs). Among the known effectors recruited by GPCR-bound arrestins are Src family kinases, which regulate cellular growth and proliferation. Here, we focus on arrestin-3 interactions with Fgr kinase, a member of the Src family. Previous reports demonstrated that Fgr exhibits high constitutive activity, but can be further activated by both arrestin-dependent and arrestin-independent pathways. We report that arrestin-3 modulates Fgr activity with a hallmark bell-shaped concentration-dependence, consistent with a role as a signaling scaffold. Here, we further demonstrate using NMR spectroscopy that a polyproline motif within arrestin-3 interacts directly with the SH3 domain of Fgr. To provide a framework for this interaction, we determined the crystal structure of the Fgr SH3 domain at 1.9 Å resolution and developed a model for the GPCR-arrestin-3-Fgr complex that is supported by mutagenesis. This model suggests that Fgr interacts with arrestin-3 at multiple sites and is consistent with the locations of disease-associated Fgr mutations. Collectively, these studies provide a structural framework for arrestin-dependent activation of Fgr.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Monobody Inhibitor Selective to the Phosphatase Domain of SHP2 and its Use as a Probe for Quantifying SHP2 Allosteric Regulation

SHP2 is a phosphatase/adaptor protein that plays an important role in various signaling pathways. Its mutations are associated with cancers and developmental diseases. SHP2 contains a protein tyrosine phosphatase (PTP) and two SH2 domains. Selective inhibition of these domains has been challenging due to the multitude of homologous proteins in the proteome. Here, we developed a monobody, synthetic binding protein, that bound to and inhibited the SHP2 PTP domain. It was selective to SHP2 PTP over close homologs. A crystal structure of the monobody-PTP complex revealed that the monobody bound both highly conserved residues in the active site and less conserved residues in the periphery, rationalizing its high selectivity. Its epitope overlapped with the interface between the PTP and N-terminal SH2 domains that is formed in auto-inhibited SHP2. By using the monobody as a probe for the accessibility of the PTP active site, we developed a simple, nonenzymatic assay for the allosteric regulation of SHP2. The assay showed that, in the absence of an activating phospho-Tyr ligand, wild-type SHP2 and the “PTP-dead” C459E mutant were predominantly in the closed state in which the PTP active site is inaccessible, whereas the E76K and C459S mutants were in the open, active state. It also revealed that previously developed monobodies to the SH2 domains, ligands lacking a phospho-Tyr, weakly favored the open state. These results provide corroboration for a conformational equilibrium underlying allosteric regulation of SHP2, provide powerful tools for characterizing and controlling SHP2 functions, and inform drug discovery against SHP2.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Kink band and shear band localization in anisotropic perfectly plastic solids

Shear-driven strain localization has been observed in a wide variety of materials and may take the form of shear bands or kink bands. Based on observations of kink bands in plastically anisotropic metallic nanolaminates and single crystal metals, here we posit that, for the specific case of isochoric deformation, the kinematics of kink band formation are indistinguishable from those of plane strain shear band formation. The only distinction between shear bands and kink bands in these systems would then be that kink bands ‘lock up’ at a particular value of material rotation while shear bands may progress to arbitrarily high strains. In order to investigate whether strong material anisotropy is sufficient to arrest shear localization at a geometry that matches the classic kink band geometry, we model the development of a band of simple shear within an anisotropic perfectly plastic material. The resulting analytical model provides the stress state needed to maintain the kinematics of simple shear as a function of material anisotropy, deformation band orientation, and shear strain (or equivalently, material rotation). It is found that plastic anisotropy can promote either kink band or shear band formation depending on the loading orientation. When the deviatoric stress is positive parallel to the plane of anisotropy, shear localization may progress without bound and a shear band is produced. When the deviatoric stress is negative parallel to the plane of anisotropy, shear localization is arrested after a certain material rotation, resulting in a kink band. Examination of the requisite applied stress state during kink band formation provides an explanation for the experimentally-observed ‘lock up’ geometry. Solutions for the band boundary inclination angle are obtained and used to provide bounds on permissible band angles for both shear bands and kink bands.

36 MATERIALS SCIENCE↗
DOE OSTI2023

NMR investigation of proton transport in polybenzimidazole/polyphosphoric acid membranes prepared via novel synthesis route

Here, in this work, we present a molecular-level view of the structural changes in polybenzimidazole (PBI) membranes doped with phosphoric acid (PA) generated by a novel membrane fabrication technique. The modified PA doped membranes displayed unprecedented ionic conductivity at elevated temperatures, comparable to the starting PBI gel membranes prepared by the commonly used PPA process, while also exhibiting enhanced mechanical properties. To elucidate the cause of these effects, we used multi-nuclear ( 1 H, 13 C, 31 P) 1D PFG, MAS and CPMAS, and 2D HETCOR magnetic resonance (NMR) techniques to characterize the structure and dynamics of the modified film and the original gel PBI membranes. 1 H diffusivity measurements show significantly enhanced proton diffusivity, both in magnitude and in lower activation energy, which were consistent with its high ionic conductivity despite the lower PA content compared to the original gel membrane. CPMAS experiments further substantiated that the location of these distinct phosphate environments as being close to the polymer network. Finally, the phosphate groups in the modified membrane were revealed to be strongly bounded to each other and to the PBI polymer backbone in contrast to the combination of weakly and strongly bounded groups of the original gel membrane.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Pattern formation in a coupled membrane-bulk reaction-diffusion model for intracellular polarization and oscillations

Reaction-diffusion systems have been widely used to study spatio-temporal phenomena in cell biology, such as cell polarization. Coupled bulk-surface models naturally include compartmentalization of cytosolic and membrane-bound polarity molecules. Here we study the distribution of the polarity protein Cdc42 in a mass-conserved membrane-bulk model, and explore the effects of diffusion and spatial dimensionality on spatio-temporal pattern formation. We first analyze a one-dimensional (1-D) model for Cdc42 oscillations in fission yeast, consisting of two diffusion equations in the bulk domain coupled to nonlinear ODEs for binding kinetics at each end of the cell. In 1-D, our analysis reveals the existence of symmetric and asymmetric steady states, as well as anti-phase relaxation oscillations typical of slow-fast systems. We then extend our analysis to a two-dimensional (2-D) model with circular bulk geometry, for which species can either diffuse inside the cell or become bound to the membrane and undergo a nonlinear reaction-diffusion process. We also consider a nonlocal system of PDEs approximating the dynamics of the 2-D membrane-bulk model in the limit of fast bulk diffusion. In all three model variants we find that mass conservation selects perturbations of spatial modes that simply redistribute mass. In 1-D, only anti-phase oscillations between the two ends of the cell can occur, and in-phase oscillations are excluded. In higher dimensions, no radially symmetric oscillations are observed. Instead, the only instabilities are symmetry-breaking, either corresponding to stationary Turing instabilities, leading to the formation of stationary patterns, or to oscillatory Turing instabilities, leading to traveling and standing waves. Codimension-two Bogdanov—Takens bifurcations occur when the two distinct instabilities coincide, causing traveling waves to slow down and to eventually become stationary patterns. Our work clarifies the effect of geometry and dimensionality on behaviors observed in ma.ss-conserved cell polarity models.

97 MATHEMATICS AND COMPUTING↗