Search NASA⌕ Search

SEARCH · Search NASA

Results for “numerical simulations”

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.

631 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS2026

The Role of Disruptive Impacts on Ocean Generation and Longevity in Icy Moons

Among icy moons of the outer Solar System, subsurface water oceans and large collisions both seem to be common, but the impact of the latter on the presence and persistence of the former is unclear and has rarely been investigated. Here we interface a smoothed-particle hydrodynamics model to simulate collisions with a thermal-structural evolution model to simulate the evolution of moons pre-collision, post-collision and without a collision. Overall, even such large-scale collisions affect only the ocean thickness or longevity, and the presence or absence of an ocean is affected for only a part of a moon’s history. In reaccreted moons, the ocean survives the impact and becomes much thicker inside larger moons with radius near 1,000 km, whereas an ocean that would otherwise arise inside moons with radius near 500 km is absent because the collision promotes ice–rock differentiation. Our simulations have not yielded an ocean developed post-impact—whether directly via collisional or reaccretional heating or indirectly through tidal heating due to collision-induced orbital changes—in a moon that would otherwise have remained frozen. The ocean-enhancing effect is pronounced only for late disruptive impacts onto large, 1,000-km-class targets, which are unlikely in recent solar system history.

Marc Neveu↗
NASA NTRS1990

A Laboratory Breadboard System for Dual-Arm Teleoperation

The computing architecture of a novel dual-arm teleoperation system is described. The novelty of this system is that: (1) the master arm is not a replica of the slave arm; it is unspecific to any manipulator and can be used for the control of various robot arms with software modifications; and (2) the force feedback to the general purpose master arm is derived from force-torque sensor data originating from the slave hand. The computing architecture of this breadboard system is a fully synchronized pipeline with unique methods for data handling, communication and mathematical transformations. The computing system is modular, thus inherently extendable. The local control loops at both sites operate at 100 Hz rate, and the end-to-end bilateral (force-reflecting) control loop operates at 200 Hz rate, each loop without interpolation. This provides high-fidelity control. This end-to-end system elevates teleoperation to a new level of capabilities via the use of sensors, microprocessors, novel electronics, and real-time graphics displays. A description is given of a graphic simulation system connected to the dual-arm teleoperation breadboard system. High-fidelity graphic simulation of a telerobot (called Phantom Robot) is used for preview and predictive displays for planning and for real-time control under several seconds communication time delay conditions. High fidelity graphic simulation is obtained by using appropriate calibration techniques.

A K Bejczy↗
DOE OSTI2026

Enter the AHU (36th Chamber of ASHRAE): A Multi-site Field Study of ASHRAE G36

Despite being recognized as the best practice for advanced building controls, ASHRAE Guideline 36 (G36) has seen slow adoption in retrofit cases. Decisionmakers lack credible field evidence to justify the time and person-power investment. Most prior analyses have relied on software simulations, which overlook implementation challenges and fail to persuade owners to move from models to real-world deployment. This paper presents a multi-site field study of G36 performance, drawing on measured results from 17 projects across diverse building types and climate zones. The analysis disaggregates outcomes by the most widely adopted air handling unit (AHU) based G36 strategies, including trim-and-respond approaches to supply air temperature (SAT) and duct static pressure (DSP) reset, and economizer controls. Results for controls re programming implementations are encouraging, with HVAC savings ranging from 2% - 49%, with a median of 18%. The range aligns with simulation study findings showing 1% - 46% savings through these three strategies. These findings show that even without capital investment in control infrastructure upgrades, existing building owners are reaping significant benefits from updating HVAC sequences of operation to industry best-practice solutions. The results give practitioners and decisionmakers a reference point on what to expect, which strategies deliver, and how field performance may compare to simulation.

Deshpande, Reva↗
DOE OSTI2026

Laser-heated diamond anvil cell synthesis and recovery of metastable MnSb 2 and YbZn 2 for post-synthesis transport studies

The creation and exploration of new materials under extreme pressure–temperature conditions has become increasingly reliant on laser-heated diamond anvil cell (LHDAC) techniques, which provide direct access to previously unexplored regions of multinary phase diagrams. Whereas numerous high-pressure phases have been identified in situ, systematic recovery and post-synthesis physical property characterization of these materials remain significant challenges. In this work, we describe the setup and implementation of an LHDAC-based synthesis and recovery workflow and demonstrate its application to metastable MnSb 2 and YbZn 2 phases. Synchrotron x-ray diffraction and spatial mapping confirm dominant formation of the targeted phases, whereas laboratory-based refinement quantifies phase fractions despite intrinsic microstrain and minor secondary phases. High-pressure transport measurements on recovered samples reveal pressure-tunable electronic instabilities in both systems. In MnSb 2 , pressure suppresses two high-temperature magnetic ordering anomalies, observed in transport, by ∼5 GPa and, for higher pressures, induces a new low-temperature feature that increases with further pressure increase. In hexagonal high-pressure YbZn 2 , an electronic reconstruction emerges at ∼11 GPa, characterized by semiconducting-like behavior from ∼30 to 300 K and a broad low-temperature coherence crossover near 30 K. Our results establish LHDAC synthesis not only as a structural discovery tool but also as an experimental platform for investigating correlated quantum states stabilized far from equilibrium thermodynamic conditions.

Huyan, S. [Iowa State University, Ames, IA (United↗
DOE OSTI2026

Neutron detector response modeling in NOvA

Neutrons can present a significant challenge for neutrino experiments in which energy reconstruction is critical. With the ability to escape detection completely and with a weak correlation between their kinetic energy and any eventual energy deposition, it is difficult to fully account for neutrons produced in neutrino interactions. This in turn leads to significant model dependence when evaluating neutron-related systematic uncertainties. The NOvA experiment is a long-baseline neutrino oscillation experiment with a high-statistics sample of antineutrino data collected by its near detector. We report an excess relative to data of simulated neutron candidates with low energy depositions when using standard Geant4 physics lists. The simulation excess is traced to an overabundance of secondary photons produced from interactions of neutrons with kinetic energy greater than \SI{20}{\mega\eV}. Improved agreement with data is obtained by applying the data-driven neutron-on-carbon \menate model for neutrons between \SI{20}{\mega\eV} and ${\sim}$\SI{100}{\mega\eV}. With \menate, the residual oversimulation is more uniform across the calorimetric neutron energy spectrum, suggesting possible overproduction of primary neutrons by the GENIE neutrino interaction generator. These results motivate the adoption of \menate-supplemented Geant4 simulation as the nominal simulation in the production of future \nova simulation.

Abubakar, S.↗
DOE OSTI2026

Realization of two-dimensional discrete time crystals with anisotropic Heisenberg coupling

A discrete time crystal (DTC) is an out-of-equilibrium phase of matter that spontaneously breaks discrete time-translation symmetry. Previous studies have been limited to a set of models with Ising-like couplings - and mostly only in one dimension - thus precluding our understanding of the existence (or not) of DTCs in models with more realistic interactions. In this work, by combining the latest generation of IBM quantum processors with state-of-the-art tensor network methods, we demonstrate the existence of a DTC in a two-dimensional system governed by anisotropic Heisenberg interactions. We uncover a rich phase diagram encompassing spin-glass, ergodic, and time-crystalline phases, and identify the interplay of initialization, interaction anisotropy, and driving protocols in stabilizing the DTC phase. By extending the study of Floquet matter beyond simplified models, we lay the groundwork for exploring how driven systems bridge the gap between quantum coherence and emergent non-equilibrium thermodynamics.

Phase transitions and critical phenomena↗
NASA NTRS1997

Analysis of Impact Induced Damage and its Effect on Structural Integrity of Space Flight Composite Overwrapped Pressure Vessels

The objective of this research work has been to provide analytical background and support to the ongoing experimental program at NASA, White Sands Test Facility, involving testing composite overwrapped pressure vessels (COPV) for impact damage and cyclic pressurization. Preliminary theoretical basis, including the governing equations for a shallow shell subjected to internal pressure, has been established. Effects of the Griffith type cracks on the structural integrity of the cylindrical vessel were evaluated by methods of Fracture Mechanics. The results indicate that the effective mass of the pressure vessel is an important factor influencing the response to impact events. We also have found that the material properties of the target, contained in the constitutive equations of the composite attached to the Aluminum liner, dominate the impact event in the low velocity range, the material properties become less important, while the target mass distribution and the impactor mass become more significant as the velocity of the impactor increases. Therefore, at high-velocity impact it is not only the kinetic energy of the impactor but also its mass which has a significant effect on the dynamics of the event, and consequently on the induced damage. This work also suggests a methodology for an assessment of the rate of loading effects on the degradation of the material toughness associated with a high-velocity impact where the rate effects become significant. To model the rate dependence of the material response a viscoelastic-plastic constitutive equations were assumed, and on this basis predictions are made regarding the rate dependent material resistance curve. Other dynamic phenomena associated with the impact event have been treated in the framework of the Computational Mechanics using the courtesy of Prof. P. Guebelle and his graduate student at University of Illinois at Urbana-Champaign who have an access to a super-fast computer located on their campus. Finally, the guidelines for a follow-up research program are provided in the body of this report. They address three major areas: theoretical research, numerical studies, and further experimental work.

Michael P Wnuk↗
DOE OSTI2026

Machine Learning for Multipactor Susceptibility Prediction in Planar RF Gaps

Multipactor discharge is a nonlinear electron avalanche that limits the performance of high-power radio-frequency (RF) and vacuum electronic devices. Predicting multipactor susceptibility traditionally relies on Monte Carlo or particle-in-cell (PIC) simulations, which become computationally expensive for large parametric studies. In this work, we present a supervised machine-learning (ML) framework for prediction of multipactor susceptibility in a two-surface planar geometry. The models are trained using high-fidelity PIC simulation generated susceptibility data and learn the relationship between operational parameters, geometry, and material-dependent secondary electron emission properties. The proposed approach enables rapid reconstruction of susceptibility charts while preserving the physical structure of multipactor growth regions.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Prospects of Nanoscience with Nanocrystals: 2025 Edition

Abstract Nanocrystals (NCs) of various compositions have made important contributions to science and technology, with their impact recognized by the 2023 Nobel Prize in Chemistry for the discovery and synthesis of semiconductor quantum dots (QDs). Over four decades of research into NCs has led to numerous advancements in diverse fields, such as optoelectronics, catalysis, energy, medicine, and recently, quantum information and computing. The last 10 years since the predecessor perspective “Prospect of Nanoscience with Nanocrystals” was published in ACS Nano have seen NC research continuously evolve, yielding critical advances in fundamental understanding and practical applications. Mechanistic insights into NC formation have translated into precision control over NC size, shape, and composition. Emerging synthesis techniques have broadened the landscape of compounds obtainable in colloidal NC form. Sophistication in surface chemistry, jointly bolstered by theoretical models and experimental findings, has facilitated refined control over NC properties and represents a trusted gateway to enhanced NC stability and processability. The assembly of NCs into superlattices, along with two-dimensional (2D) photolithography and three-dimensional (3D) printing, has expanded their utility in creating materials with tailored properties. Applications of NCs are also flourishing, consolidating progress in fields targeted early on, such as optoelectronics and catalysis, and extending into areas ranging from quantum technology to phase-change memories. In this perspective, we review the extensive progress in research on NCs over the past decade and highlight key areas where future research may bring further breakthroughs.

Chemistry↗
NASA NTRS2026

The Canadian Galactic Emission Mapper: A New 8-10 GHz Telescope To Map the Polarization of the Northern Sky

The Canadian Galactic Emission Mapper (CGEM) radio telescope is mapping polarized Galactic foregrounds from 8-10 GHz to aid in the search for B-modes in the Cosmic Microwave Background (CMB). Here we describe the design and early on-sky characterization of the CGEM optics. CGEM employs an on-axis, hat feed optical design that is a body of revolution (BOR) and exhibits excellent intrinsic polarization purity. We describe how we have optimized the optics with a novel framework that can directly minimize simulated intensity to polarized intensity (T → P) leakage in angular power spectrum space. The optimized optics exhibit simulated T → P leakage that is orders of magnitude below anticipated CℓBB for r = 10 −3 when scaled from 8-10 GHz to the CMB observing window near 95 GHz. We go on to describe the mechanical design of the optics, including a novel secondary mirror support made from Astroquartz composite that has low loss, low dielectric constant, and preserves the BOR symmetry of the optics. We then showcase the early on-sky performance of the optics with observations of the Sun and satellites, which probe the beam to 40dB down from the peak. A beam model for the deployed CGEM optics based on electromagnetic simulations is in excellent agreement with the on-sky data.

cosmic microwave background↗
DOE OSTI1995

Nuclear Safety [Vol. 36, No. 2, July-December 1995]

Nuclear Safety is a journal that covers significant issues in the field of nuclear safety. Its primary scope is safety in the design, construction, operation, and decommissioning of nuclear power reactors worldwide and the research and analysis activities that promote this goal, but it also encompasses the safety aspects of the entire nuclear fuel cycle, including fuel fabrication, spent-fuel processing and handling, and nuclear waste disposal, the handling of fissionable materials and radioisotopes, and the environmental effects of all these activities. Table of Contents for this issue follows. THE CHORNOBYL ACCIDENT: 195 The Chornobyl Accident Revisited, Part III: Chernobyl Source Term Release Dynamics and Reconstruction of Events During the Active Phase, A. R. Sich; GENERAL SAFETY CONSIDERATIONS: 218 Second ANS Workshop on the Safety of Soviet-Designed Nuclear Power Plants, R. A. Bari; 234 Elements of a Nuclear Criticality Safety Program, C. M. Hopper; 243 Rickover, Excellence, and Criticality Safety Programs, R. E. Wilson; ACCIDENT ANALYSIS: 249 Transient Analysis of the PIUS Advanced Reactor Design with the TRAC-PF1/MOD2 Code, B. E. Boyack, J. L Steiner, S. C. Harmony, H. J. Stumpf, and J. F. Lime; 278 The Hierarchy-By-Interval Approach to Identifying Important Models that Need Improvement in Severe-Accident Simulation Codes, T. J. Heames, M. Khatib-Rahbar, J. E. Kelly, R. P. Jenks-Johnson, and Y.-S. Chen; 290 RELAP5/MOD3 Code Coupling Model, R. P. Martin; 299 Missiles Caused by Severe Pressurized-Water Reactor Accidents, R. Krieg; DESIGN FEATURES: 310 Validation of COMMIX with Westinghouse AP-600 PCCS Test Data, J. G. Sun, T. H. Chien, J. Ding, and W T. Sha; ENVIRONMENTAL EFFECTS: 321 Spent Nuclear Fuel Characterization for a Bounding Reference Assembly for the Receiving Basin for Off-Site Fuel, S. D. Kahook, R. L. Garrett, L. R. Canas, and M J. Beckum; OPERATING EXPERIENCES: 332 Reactor Shutdown Experience, Compiled by J. W. Cletcher; U.S. NUCLEAR REGULATORY COMMISSION INFORMATION AND ANALYSES: 335 Reactor Coolant System Blowdown at Wolf Creek on September 17, 1994, J. V. Kauffman and S. L. Israel; RECENT DEVELOPMENTS: 344 Reports, Standards, and Safety Guides, D. S. Oueener; 349 Proposed Rule Changes as of June 30, 1995; ANNOUNCEMENTS: 320 Symposium on Acceptability of Risk From Radiation—Application to Manned Space Flight; 320 24th DOE/NRC Nuclear Air Cleaning and Treatment Conference; 361 Radiation Biology and Radiation Protection— Modern Developments and Tendencies in Radiation Biology; 362 1997 IEEE Sixth Conference on Human Factors and Power Plants; 354 The Authors; 360 Reviewers of Nuclear Safety, Vol. 36.

05 NUCLEAR FUELS↗
DOE OSTI2026

Capturing Secondary Kinetic Instabilities in Three‐Dimensional Dayside Reconnection Using an Improved Gradient‐Based Closure

Magnetic reconnection is a highly dynamic process that excites a wide variety of kinetic waves and instabilities. Transverse current sheet instabilities such as the lower-hybrid drift and secondary drift-kink instabilities in particular have been shown by kinetic simulations to modify the reconnection and introduce significant turbulence and mixing to the reconnection layer. Past studies using the ten-moment fluid model to capture important kinetic physics such as the electron inertia and full representation of the pressure tensor proved advantageous to a two-fluid representation of reconnection, but the model struggled when using a local relaxation closure for the heat flux to replicate the current sheet instabilities and subsequent mixing seen in kinetic simulations. This work uses the Gkeyll software framework to perform simulations of asymmetric reconnection based on the 16 October 2015 MMS crossing of a diffusion region, the Burch event. An improved gradient-based heat flux closure is implemented, showing significant improvement in secondary kinetic instabilities that grow in the current sheet. These instabilities generate turbulence which leads to growth of secondary magnetic islands and flux ropes.

Bradshaw, K. [Princeton University, NJ (United Sta↗
DOE OSTI2026

String-Breaking Dynamics in Quantum Adiabatic and Diabatic Processes

Confinement prohibits isolation of color charges, e.g., quarks, in nature via a process called string breaking : the separation of two charges results in an increase in the energy of a color flux, visualized as a string, connecting those charges. Eventually, creating additional charges is energetically favored, hence breaking the string. Such a phenomenon can be probed in simpler models, including quantum spin chains, enabling enhanced understanding of string-breaking dynamics. A challenging task is to understand how string breaking occurs as time elapses, in an out-of-equilibrium setting. This work establishes the phenomenology of dynamical string breaking induced by a gradual increase of string tension over time. It, thus, goes beyond instantaneous quench processes and enables tracking the real-time evolution of strings in a more controlled setting. We focus on domain-wall confinement in a family of quantum Ising chains. Our results indicate that, for sufficiently short strings and slow evolution, string breaking can be described by the transition dynamics of a two-state quantum system akin to a Landau-Zener process. For longer strings, a more intricate spatiotemporal pattern emerges: the string breaks by forming a superposition of bubbles (domains of flipped spins of varying sizes), which involve highly excited states. We finally demonstrate that string breaking driven only by quantum fluctuations can be realized in the presence of sufficiently long-ranged interactions. This work holds immediate relevance for studying string breaking in quantum-simulation experiments.

Ising model↗
DOE OSTI2023

The brighter-fatter effect in the JWST MIRI Si:As IBC detectors: I. Observations, impact on science, and modeling

Context.The Mid-Infrared Instrument (MIRI) on board theJames WebbSpace Telescope (JWST) uses three Si:As impurity band conduction (IBC) detector arrays. The output voltage level of each MIRI detector pixel is digitally recorded by sampling up the ramp. For uniform or low-contrast illumination, the pixel ramps become nonlinear in a predictable way, but in areas of high contrast, the nonlinearity curve becomes much more complex. The origin of the effect is poorly understood and currently not calibrated out of the data. Aims.We provide observational evidence of the brighter-fatter effect (BFE) in MIRI conventional and high-contrast coronagraphic imaging, low-resolution spectroscopy, and medium-resolution spectroscopy data, and we investigate the physical mechanism that gives rise to the effect on the MIRI detector pixel raw voltage integration ramps. Methods.We used public data from the JWST/MIRI commissioning and Cycle 1 phase. We also developed a numerical electrostatic model of the MIRI detectors using a modified version of the publicPoisson_CCDcode. Results.We find that the physical mechanism behind the BFE manifesting in MIRI data is fundamentally different to that of charge-coupled devices and photodiode arrays such as the Hawaii-XRG near-infrared detectors used by the NIRISS, NIRCam, and NIRSpec instruments on board JWST. Observationally, the BFE makes the JWST MIRI data yield 10–25% larger point sources and spectral line profiles as a function of the relative level of de-biasing of neighboring detector pixels. This broadening impacts the MIRI absolute flux calibration, time-series observations of faint companions, and point spread function modeling and subtraction. We also find that the intra-pixel 2D profile of the shrinking Si:As IBC detector depletion region directly impacts the accuracy of the pixel ramp nonlinearity calibration model.

Astronomy & Astrophysics↗
DOE OSTI2026

Modeling Flow in a Microreactor Core: From High-Fidelity CFD to Subchannel Analysis

The accurate prediction of pressure drop in tightly packed, low–Reynolds number (Re) bare rod bundles is essential for the thermal-hydraulic design of the Microreactor Applications Research Validation and EvaLuation (MARVEL) reactor and other microreactor concepts. However, existing friction factor correlations, particularly the upgraded Cheng-Todreas (UCTD) correlation, have limited validation for the small pitch-to-diameter ratios (P/Ds) and transitional flow conditions characteristic of these systems. In this work, we perform high-fidelity large-eddy simulations (LESs) of both an infinite bare rod bundle and a finite bare 37-pin scalloped rod bundle across the range 1000≤ 𝑅𝑒 ≤5500. The simulations reveal strong gap vortex–driven transitional behavior and indicate that the UCTD may underpredict the friction factor by up to 28% at 𝑃/𝐷 = 1.05. Using the LES-calculated pressure drops, we formulate a new friction factor correlation that follows the Cheng-Todreas functional form but is calibrated for low-Re and tightly packed geometries representative of a MARVEL-like reactor. The correlation is implemented in the MOOSE (Multiphysics Object-Oriented Simulation Environment) subchannel module and compared against both the LES and UCTD predictions. Across all subchannel types, the proposed model reduces the streamwise velocity differences from as high as 44% (UCTD) to below 9%, and decreases the pressure gradient differences from 13% to 25% (UCTD) to 0.7% to 7% relative to the LES results. These results suggest that the new correlation has the potential to improve the pressure drop and flow field predictions for such geometries, highlighting the importance of high-fidelity simulations in supporting microreactor thermal-hydraulic model development and motivating future pressure drop experiments for compact rod bundles to further validate these findings.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
NASA NTRS1992

Monte Carlo modeling of atomic oxygen attack of polymers with protective coatings on LDEF

Characterization of the behavior of atomic oxygen interaction with materials on the Long Duration Exposure Facility (LDEF) will assist in understanding the mechanisms involved, and will lead to improved reliability in predicting in-space durability of materials based on ground laboratory testing. A computational simulation of atomic oxygen interaction with protected polymers was developed using Monte Carlo techniques. Through the use of assumed mechanistic behavior of atomic oxygen and results of both ground laboratory and LDEF data, a predictive Monte Carlo model was developed which simulates the oxidation processes that occur on polymers with applied protective coatings that have defects. The use of high atomic oxygen fluence-directed ram LDEF results has enabled mechanistic implications to be made by adjusting Monte Carlo modeling assumptions to match observed results based on scanning electron microscopy. Modeling assumptions, implications, and predictions are presented, along with comparison of observed ground laboratory and LDEF results.

Bruce A. Banks↗
NASA NTRS2026

Scientific Highlights and Status of the Mars Science Laboratory Mission, Years 9 through 12

We present the scientific highlights and status from years 9-12 of NASA’s Mars Science Laboratory mission, with its Curiosity rover. A thermochemolysis-based sample analysis revealed a variety of organic molecules preserved in clay mineral-rich lake/lake margin deposits of the Glen Torridon region. The stratigraphically higher clay-sulfate mineral transition on Aeolis Mons was found to be associated with a lithological change from lacustrine to dry aeolian depositional environments. Water intermittently returned to the surface after the increase in aridity, including during the deposition of wave ripple strata which indicate a climate capable of sustaining an ice-free lake in the early Hesperian. The overlying Mg sulfate-bearing unit has both poly- and monohydrated Mg sulfates and is also enriched in iron carbonate of a purity consistent with chemical sedimentation and therefore atmospheric sequestration. A channel-incised canyon within Gediz Vallis and distal alluvial deposits recorded numerous fluvial and debris flow events after the deposition, lithification, and erosion of the clay- and sulfate-bearing units. Bright clasts found in Gediz Vallis canyon are composed of native sulfur. The investigation of the modern atmosphere and environment focused on the influence of the crater and mountain topography and captured high-altitude noctilucent and iridescent clouds. Radiation monitoring now spans an 11-year solar cycle. Rover systems and science instruments remain capable of addressing mission scientific objectives. All ten scientific instruments continue to return high-value data despite degradation of the Chemistry and Camera instrument’s high-voltage system, loss of multispectral imaging capability, and loss of active neutron spectroscopy. The mission has mitigated the loss of redundant brakes on two arm joints, additional wheel wear, and highly reduced capability of the redundant flight computer. Efficiency improvements in rover energy utilization and operations staffing have allowed the mission team to maintain a high level of productivity despite declining rover power generation and funding.

Mars↗
DOE OSTI2026

Development of an Amine Oxide Polyzwitterion Brush Martini Model with Polarizable Water and Ions

Abstract In order to accurately simulate the fouling process of proteins onto polyzwitterion brushes, models that accurately capture the hydration properties and chain conformations of such brushes must first be established. We developed a Martini coarse-grained (CG) model for amine oxide polyzwitterion (PNOMA) brushes, a promising class of antifouling materials, in polarizable water and ions by fitting to all-atom bond and angle distributions, monomer hydration free energy, monomer–monomer distance potential of mean force (PMF), and monomer–salt radial distribution functions (RDFs). Martini 2.2P was selected for compatibility with the established polarizable water and ion models. For comparison with PNOMA, we also constructed models for conventional sulfobetaine (PSBMA) and phosphorylcholine (PMPC) polyzwitterions and the polycation PMETAC using established nonbonded bead types from the literature and refitting bond and angle potentials. We simulated each polymer brush chemistry for varying grafting density and chain length, validating brush height scaling relations against experimental data. The CG models captured the relative hydration strengths among different polyzwitterion chemistries, and brush heights extrapolated to higher molecular weights are in agreement with experimental ellipsometry data. We find that chain swelling of the superhydrophilic PNOMA brushes lies between that of the traditional polyzwitterions PSBMA/PMPC and the polycation PMETAC. For PNOMA brushes in NaCl solution, simulated brush height decreases with salt concentration due to the selectively strong interactions between amine oxide and sodium ions.

Walker, Christopher C. [Oak Ridge National Laborat↗