Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reproducibility”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Lattice simulations of inflation

The scalar field theory of cosmological inflation constitutes nowadays one of the preferred scenarios for the physics of the early universe. In this paper we aim at studying the inflationary universe making use of a numerical lattice simulation. Various lattice codes have been written in the last decades and have been extensively used for understating the reheating phase of the universe, but they have never been used to study the inflationary phase itself far from the end of inflation (i.e. about 50 e-folds before the end of inflation). In this paper we use a lattice simulation to reproduce the well-known results of some simple models of single-field inflation, particularly for the scalar field perturbation. The main model that we consider is the standard slow-roll inflation with an harmonic potential for the inflaton field. We explore the technical aspects that need to be accounted for in order to reproduce with precision the nearly scale invariant power spectrum of inflaton perturbations. We also consider the case of a step potential, and show that the simulation is able to correctly reproduce the oscillatory features in the power spectrum of this model. Even if a lattice simulation is not needed in these cases, that are well within the regime of validity of linear perturbation theory, this sets the basis to future work on using lattice simulations to study more complicated models of inflation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Investigation of fast particle redistribution induced by sawtooth instability in NSTX-U

The effects of sawtooth on fast ion transport have been examined in reproducible, 2s long sawtoothing L-mode discharges during the 2016 experimental campaign on National Spherical Torus Experiment Upgrade (NSTX-U) [J.E. Menard et al. Nucl. Fusion (2012) 083015]. Analysis of the discharges showed that standard sawtooth models (full/partial reconnection models) in the TRANSP code were not capable to fully reproduce the fast ion redistribution induced by sawtooth crashes. Some global parameters such as neutron rate can be recovered while detailed features, e.g. distribution functions, estimated using the models were different from the experimental observation. The standard sawtooth models in TRANSP do not take into account the different effect of sawtooth crashes depending on fast ion energy and orbit type and that may cause the disagreement between experiments and simulations. In this work, the newly developed kick model has been applied to replace the standard sawtooth models for the fast ion transport. TRANSP simulation results using the kick model, taking into account the characteristics of fast ion such as energy and pitch angle, can reproduce experimental neutron rates within 10\% difference. The qualitative comparison of the measurements and synthetic diagnostics of FIDA and SSNPA using the TRANSP simulation results with kick model demonstrates good agreements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Modeling of ExB effects on tungsten re-deposition and transport in the DIII-D divertor

Mixed-material DIVIMP-WallDYN modelling, now incorporating ExB drifts, is presented that simultaneously reproduces tungsten (W) erosion and deposition patterns observed during the DIII-D Metal Rings Campaign, in which a toroidally symmetric set of W-coated tiles were installed in the carbon (C) DIII-D divertor. Since most reactor plasma facing component (PFC) designs call for mixed-material environments, including ITER’s W/Be enviroment, the divertor targets will quickly evolve into reconstitued surfaces of multiple elements. This work identifies controlling physics that affects material migration patterns in the divertor, which impact PFC lifetimes and impurity leakage from the divertor to the core. These simulations indicate that radial and poloidal ExB transport dominates over parallel force balance for high-Z impurities such as W in the divertor region of DIII-D. It is demonstrated that ExB drifts are required to reproduce the experimental observation of non-local W and C co-accumulation in a band ~7-9 cm outboard of the outer-strike-point (OSP) W source, for attached Lmode conditions in the unfavorable ion grad-B drift direction. In addition, W gross erosion is localized to the region outboard of the OSP, as the formation of C co-deposits suppresses W erosion at the strike point. Time-dependent simulations with scaled ExB impurity drifts (60% of the OEDGE-calculated drift velocity) and W re-erosion quantitatively reproduce these features, including depth-resolved W/C ratios, within a factor of 2 over ~115 seconds of accumulated plasma exposure. The location of co-deposition regions is shown to be well represented by an analytic leakage model, driven largely by poloidal ExB drifts. Qualitative agreement is also found between campaign-integrated W deposition measurements and simulations for the favorable ion grad-B drift direction, the standard mode of operation for most tokamaks. Furthermore, these results imply that a longterm inward radial migration of material from the outer divertor through the private flux region may occur in future devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Modelling the Alfvén eigenmode induced fast-ion flow measured by an imaging neutral particle analyzer

An imaging neutral particle analyzer (INPA) provides energy and radially resolved measurements of the confined fast-ion population ranging from the high-field side to the edge on the midplane of the DIII-D tokamak. In recent experiments, it was used to diagnose fast-ion flow in the INPA-interrogated phase-space driven by multiple, marginally unstable Alfvén eigenmodes (AEs). Here, the key features of this measured fast-ion flow are: (I) a fast-ion flow from q min and the injection energy (81 keV) towards lower energies and plasma periphery.(II) A flow from the same location towards higher energies and the plasma core, (III) a phase-space 'hole' at the injected energy and plasma core and (IV) a pile-up at the plasma core at lower energies (~60 keV). Ad hoc energetic particle diffusivity modelling of TRANSP significantly deviates from the observation. Comparably, a reduced modelling, i.e. a combination of NOVA-K and ASCOT5 code with the measured mode structure and amplitude, generally reproduce some key features of the observed phase-space flow, but largely failed to interpret fast ion depletion near the plasma axis. At last, self-consistent, first-principle multi-phase hybrid simulations that include realistic neutral beam injection and collisions are able to reproduce most features of the time-resolved phase-space flow. During consecutive hybrid phases, an RSAE consistent with the experiment grows and saturates, redistributing the injected fast ions. The resulting synthetic INPA images are in good agreement with the measurement near the injection energy. The simulations track the fast-ion redistribution within the INPA range, confirming that the measured fast-ion flow follows streamlines defined by the intersection of phase-space surfaces of constant magnetic moment μ and constant E' = nE + ωP φ , where n and ω are the instability toroidal mode number and frequency, and E and P φ the ion energy and toroidal canonical momentum. Nonperturbative effects are required to reproduce the depletion of fast ions near the magnetic axis at the injection energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Radiation induced non-linear oscillations in ITER baseline scenario plasmas in DIII-D

Abstract This work shows how the radiation brought about by metals or metal-equivalent radiators such as Kr and Xe produces non-linear dynamics on otherwise stationary β N flattops of DIII-D ITER Baseline Scenario demonstration discharges. The Kr and Xe gases are used to reproduce the radiative loss rates of W in present machines that operate at core temperatures much lower than the expected ITER temperature. Experiments on DIII-D with injection of Kr and Xe, as well as with sources of intrinsic metals reach the range of radiated fraction values expected in the ITER core and experience slow oscillations in temperature and radiated power. In many cases of high radiated fraction, the core temperature decreases enough for the safety factor profile to rise above the 1/1 rational surface, naturally eliminating sawteeth and occasionally producing a persistent helical core. The oscillations can be reproduced by a modified Lotka–Volterra system for temperature and radiated fraction if diffusion and noise are included, which indicates that the interplay between temperature and radiation can be the main cause of the cyclic nature of the system. A new physics based model which includes equations for temperature, density and input power can also reproduce the oscillations observed in the experiments. The present results suggest that the non-linearity of the system can be increased by the inclusion of the inherently non-linear alpha heating term, which is proportional to ∼ n e 2 T i 2 , and obtains oscillations in the model when added to an otherwise more stationary system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Validation of an integrated modeling framework for investigating 3D plasma responses in tokamak plasmas

As contemporary experimental tokamaks are pushed toward reactor-relevant operation, they provide essential testbeds for demonstrating and exploring ELM control strategies for deployment in future fusion pilot plants. Accurate predictions of full plasma responses are essential to guide and optimize these demonstrations. This paper introduces and validates an integrated modeling framework over a historical range of DIII-D operational space. The integrated modeling framework uses only scalar plasma parameters and optional reference boundary to self-consistently and flexibly scan through tokamak operational space and estimate the corresponding ELM-suppression relevant plasma response over 3D coil phase space. The framework generates tightly converged equilibria that satisfy a target set of plasma parameters (I p , β N , l i ), with kinetic profiles constrained by an EPED(NN)-computed pedestal and empirical core model. The plasma response of these modeled equilibria is calculated with GPEC. Validated against 55 distinct DIII-D equilibria, the framework consistently reproduces experimentally constrained equilibrium pressure, q, and other representative profiles using the time-varying information from only evolving scalar plasma parameters; plasma response validation is performed for the n=3 perturbation. In particular, the plasma response of a DIII-D discharge scanning q 95 in search of ELM suppression windows is quantitatively reproduced. By accurately reproducing experimental equilibria and plasma responses across wide parameter variation, the framework supports potential for synthetic parameter scans in key operational, stability, and plasma response dimensions to investigate RMP ELM-suppression experiments and inform predictive RMP scenario optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

DESI mock challenge: constructing DESI galaxy catalogues based on FastPM simulations

Together with larger spectroscopic surveys such as the Dark Energy Spectroscopic Instrument (DESI), the precision of large-scale structure studies and thus the constraints on the cosmological parameters are rapidly improving. Therefore, one must build realistic simulations and robust covariance matrices. We build galaxy catalogues by applying a halo occupation distribution (HOD) model upon the FastPM simulations, such that the resulting galaxy clustering reproduces high-resolution N-body simulations. While the resolution and halo finder are different from the reference simulations, we reproduce the reference galaxy two-point clustering measurements – monopole and quadrupole – to a precision required by the DESI Year 1 emission line galaxy sample down to non-linear scales, i.e. $k\lt 0.5\, h\, \mathrm{Mpc}^{-1}$ or $s\gt 10\, \mathrm{Mpc}\, h^{-1}$. Furthermore, we compute covariance matrices based on the resulting FastPM galaxy clustering – monopole and quadrupole. We study for the first time the effect of fitting on Fourier conjugate (e.g. power spectrum) on the covariance matrix of the Fourier counterpart (e.g. correlation function). We estimate the uncertainties of the two parameters of a simple clustering model and observe a maximum variation of 20 per cent for the different covariance matrices. Nevertheless, for most studied scales the scatter is between 2 and 10 per cent. Consequently, using the current pipeline we can precisely reproduce the clustering of N-body simulations and the resulting covariance matrices provide robust uncertainty estimations against HOD fitting scenarios. We expect our methodology will be useful for the coming DESI data analyses and their extension for other studies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

A Novel Dosimetry Method for Small Animal Irradiators Using 3D-printed Mouse Phantoms and Alanine Dosimeters

Abstract Accurate dosimetry is a crucial component of small animal and preclinical irradiation studies. Various dosimetry options are available but fail to characterize complex geometries and variable energy spectra of modern x-ray irradiators accurately. These options also lack national/international standards of recognition. This paper presents a novel dosimetry system, Dosequate, which uses murine phantoms embedded with alanine dosimeters and paired with x-ray energy spectra corrections to deliver consistent and accurate dosimetry measurements. This study compares Dosequate measurements against a Precision X-RAD 320 internal ion chamber and the treatment planning system of an Xstrahl Small Animal Radiation Research Platform. Results demonstrate the accuracy and reproducibility of the Dosequate system, highlighting its potential for standardizing dosimetry in preclinical research throughout the industry. Results demonstrate the accuracy and reproducibility of the Dosequate system, highlighting its potential for standardizing dosimetry in preclinical research throughout the industry. The Dosequate dosimetry method provides a robust and standardized approach for measuring absorbed dose in small animal irradiators. Its accuracy, reproducibility, and ability to account for complex irradiation geometries make it a valuable tool for preclinical research. This system has the potential to significantly improve the intercomparability of studies across different facilities and enhance the reliability of results.

Environmental Sciences & Ecology↗
DOE OSTI2025

Towards a Deeper Fundamental Understanding of (Al,Sc)N Ferroelectric Nitrides

Density functional theory (DFT) calculations, within the virtual crystal alloy approximation, are performed, along with the development of a Landau-type model employing a symmetry-allowed analytical expression of the internal energy and having parameters determined from first principles, to investigate properties and energetics of Al1-xScxN ferroelectric nitrides in their hexagonal forms. These DFT computations and this model predict the existence of two different types of minima, namely, the fourfold-coordinated wurtzite (WZ) polar structure and a five-fold coordinated paraelectric hexagonal phase (denoted as H5), for any Sc composition up to 40%. The H5 minimum progressively becomes the lowest-energy state within hexagonal symmetry as the Sc concentration increases from 0 to 0.4. Furthermore, the model points to several key findings. Examples include the crucial role of the coupling between polarization and strains to create the WZ minimum, in addition to polar and elastic energies, and that the origin of the H5 state overcoming the WZ phase as the global minimum within hexagonal symmetry when increasing the Sc composition mostly lies in the compositional dependency of only two parameters-one linked to the polarization and another one being purely elastic in nature. Other examples are that forcing Al1-xScxN systems to have no or a weak change in lattice parameters when heating them allows us to reproduce their finite-temperature polar properties well and that a value of the axial ratio close to that of the ideal WZ structure implies a large polarization at low temperatures but not necessarily at high temperatures because of the ordered-disordered character of the temperature-induced formation of the WZ state. Such findings should allow for a better fundamental understanding of (Al,Sc)N ferroelectric nitrides, which may be used to design efficient devices having, e.g., low operating voltages.

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

Spin-orbit coupling effects on predicting defect properties with hybrid functionals: A case study in CdTe

Defect formation energies and transition levels are critical in determining doping behavior and recombination in semiconductor applications. Hybrid functionals are often used to overcome the band gap and delocalization errors of standard density-functional theory, and it is tempting to presume that the defect properties are correctly predicted once the hybrid functional mixing parameter reproduces the experimental band gap. However, pronounced spin-orbit coupling (SOC) effects can have an additional important role, which is clearly shown in this work by analyzing SOC effects originating from the Te-p orbitals in CdTe. In this work, we therefore use a hybrid functional that reproduces the experimental band gap when SOC is included, requiring a larger mixing parameter a=0.33 compared to the conventional choice of a=0.25. This hybrid functional was then used to predict defect properties, e.g., formation energy, transition level, and defect equilibrium. For defect states that do not directly involve the Te-p orbitals, such as the Cd interstitial (Cdi), we find that the effect of SOC on the defect levels can be captured by simply considering the SOC-induced band-edge shift. This is not the case for the A center (ClTe-VCd defect pair), where the localized acceptor state formed by Te-p orbitals is more directly affected by the SOC. For this defect, a mixing parameter as large as a=0.40 is required to reproduce the experimental acceptor level. Regarding the implications for photovoltaics, we suggest that the Cdi, which is the dominant compensating donor, could play an important role as a recombination center. While Cdi is usually thought of as a benign shallow donor, our predicted defect levels in the fully band-gap-corrected calculations are deep enough to raise a concern, and we propose a recombination mechanism for electron capture by Cdi.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Prompt-fission-neutron spectra in the 239 Pu( n,f ) reaction

Prompt-fission-neutron spectra from Pu 239 ( n , f ) were measured with respect to Cf 252 spontaneous fission for incident neutron energies from 0.7 to 700 MeV at the Weapons Neutron Research facility (WNR) of the Los Alamos Neutron Science Center. A newly designed high-efficiency fission chamber was coupled to the highly segmented Chi-Nu neutron liquid scintillator array to detect neutrons emitted in fission events. The double time-of-flight technique was used to deduce the incident neutron energies from the spallation target and the outgoing-neutron energies from the fission chamber. Prompt-fission-neutron spectra (PFNS) were measured with respect to Cf 252 spontaneous fission down to 200 keV and up to about 12 MeV for all the incident neutron energies with typical total uncertainties well below 2 % up to about 7-MeV outgoing-neutron energy. The general trend of PFNS is well reproduced by JEFF3.3 and ENDF evaluations, although a better agreement is found with JEFF3.3. Discrepancies were observed for the low-energy part of the spectra, especially around the opening of the second-, third- and fourth-chance fission. Additionally, neutron average kinetic energies as a function of incident neutron energy are obtained experimentally with reported total uncertainties below 0.5 % . The measured values agree with the most recent data. The trend is fairly well reproduced by the JEFF3.3 evaluation, although it fails to reproduce the experimental values within their uncertainties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Extracting the number of short-range correlated nucleon pairs from inclusive electron scattering data

The extraction of the relative abundances of short-range correlated (SRC) nucleon pairs from inclusive electron scattering is studied using the generalized contact formalism (GCF) with several nuclear interaction models. GCF calculations can reproduce the observed scaling of the cross-section ratios for nuclei relative to deuterium at high xB and large Q 2 , a 2 =(σ A /A)/(σ d /2). In the nonrelativistic instant-form formulation, the calculation is very sensitive to the model parameters and only reproduces the data using parameters that are inconsistent with ab initio many-body calculations. Using a light-cone GCF formulation significantly decreases this sensitivity and improves the agreement with ab initio calculations. The ratio of similar mass isotopes, such as 40 Ca and 48 Ca, should be sensitive to the nuclear asymmetry dependence of SRCs, but is found to also be sensitive to low-energy nuclear structure. Thus the empirical association of SRC pair abundances with the measured a2 values is only accurate to about 20%. Finally, improving this will require cross-section calculations that reproduce the data while properly accounting for both nuclear structure and relativistic effects.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Unifying disparate rate-dependent rheological regimes in non-Brownian suspensions

A typical dense non-Brownian particulate suspension exhibits shear thinning (decreasing viscosity) at low shear rate/stress followed by a Newtonian plateau (constant viscosity) at intermediate shear rate/stress values which transitions to shear thickening (increasing viscosity) beyond a critical shear rate/stress value and finally, undergoes a second shear thinning transition at extremely high shear rate/stress values. In this study, we unify and quantitatively reproduce all the disparate rate dependent regimes and the corresponding transitions for a dense non-Brownian suspension with increasing shear rate/stress. Here, we employ discrete particle dynamics simulations based on the proposed mechanism to elucidate its accuracy. We find that a competition between inter-particle interactions of hydrodynamic and non-hydrodynamic origins and the switching in the dominant stress scale with increasing the shear rate/stress lead to each of the above transitions. Inclusion of traditional hydrodynamic interactions, attractive/repulsive DLVO (Derjaguin and Landau, Verwey and Overbeek) interactions, the inter-particle contact interactions and a constant friction (or other constraint mechanism) reproduce the initial thinning as well as the shear thickening transition. However, to quantitatively capture the intermediate Newtonian plateau and the second shear thinning, an additional non-hydrodynamic interaction of non-DLVO origin and a decreasing coefficient of friction, respectively, are essential; thus, providing the first explanation for the presence of the intermediate Newtonian plateau along with reproducing the second shear thinning in a single framework. Expressions utilized for various interactions and friction are determined from experimental measurements and hence, result in an excellent quantitative agreement between the simulations and previous experiments.

42 ENGINEERING↗
DOE OSTI2021

Transition between Instability and Seeded Self-Modulation of a Relativistic Particle Bunch in Plasma

We use a relativistic ionization front to provide various initial transverse wakefield amplitudes for the self-modulation of a long proton bunch in plasma. We show experimentally that, with sufficient initial amplitude [ ≥ ( 4.1 ± 0.4) MV / m ], the phase of the modulation along the bunch is reproducible from event to event, with 3%–7% (of 2 π ) rms variations all along the bunch. The phase is not reproducible for lower initial amplitudes. We observe the transition between these two regimes. Phase reproducibility is essential for deterministic external injection of particles to be accelerated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

In situ measurement and modelling of the growth and length scale of twins in α -uranium

Twinning is an important deformation mechanism in crystalline materials. Twins nucleate, then grow, forming a discrete pattern of deformation bands with a certain thickness. The mechanism for twin growth is not completely understood. We present a simple phase field model which captures the physics of twin growth and reproduces the details observed during in situ electron backscatter diffraction experiments on α-uranium. The interaction between residual dislocations at twin interfaces and mobile dislocations in untwinned regions increases the stress needed for twin growth. This is described by a nonlocal term in the proposed constitutive equations: the nucleation stress of a twin increases proportionally to the twin phase field in a neighborhood. Competition between slip and twinning favors the nucleation of a new twin rather than the growth of a preexisting twin. The phase field model is coupled with a crystal plasticity finite element solver that includes the plastic deformation due to twinning. Simulations are able to reproduce the number of twins and their thicknesses as a function of the strain observed during in situ electron backscatter diffraction experiments. The stress-strain curve is also reproduced. This comparison allows the values of the nucleation stress and the interaction strength between residual and mobile dislocations to be found. It also gives an estimation of the density of residual dislocations that is consistent with the observed twin thickness. Finally, this model can be applied to understand microstructural effects in materials using twinning as a strengthening mechanism.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Calculation of Young’s Modulus of MoS{sub 2}-Based Single-Wall Nanotubes Using Force-Field and Hybrid Density Functional Theory

A force field is proposed that reproduces with a high accuracy a large number of properties of the bulk crystal MoS{sub 2} phases, monolayers, and nanotubes. The reproduced values are both the experimental results and the results of quantum chemical calculations. The elaborated interaction potential can be useful primarily for investigation of multiwall MoS{sub 2} nanotubes and their thermodynamic properties, especially, since the potential is able to reproduce the frequencies of the crystal phonon spectrum. In this study the proposed potential is applied to simulate the temperature dependence of a number of properties of the armchair and zigzag nanotubes. The calculations have been performed using molecular mechanics method within the framework of quasi harmonic approximation which is carried out through the estimation of the temperature dependence of the Helmholtz free energy.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Light $$Z^\prime $$ and Dirac fermion dark matter in the $$B-L$$ model

Abstract We consider a $$U(1)_{B-L}$$ U ( 1 ) B - L model with a $$Z^\prime $$ Z ′ portal Dirac fermion dark matter (DM) $$\chi $$ χ of low mass which couples very weakly to the $$B - L$$ B - L gauge boson $$Z^\prime $$ Z ′ . An arbitrary $$B-L$$ B - L charge $$Q\ne \pm 1, \pm 3$$ Q ≠ ± 1 , ± 3 of the DM $$\chi $$ χ ensures its stability. Motivated by the sensitivity reach of forthcoming “Lifetime Frontier” experiments, we focus on the $$Z^\prime $$ Z ′ mass, $$m_{Z^\prime }$$ m Z ′ , in the sub-GeV to few GeV range. To evaluate the DM relic abundance, we examine both the freeze-out and freeze-in DM scenarios. For the freeze-out scenario, we show that the observed DM abundance is reproduced near the $$Z^\prime $$ Z ′ resonance, $$m_\chi \simeq m_{Z^\prime }/2$$ m χ ≃ m Z ′ / 2 , where $$m_\chi $$ m χ is the DM mass. For the freeze-in scenario, we focus on $$m_\chi \ll m_{Z^\prime }$$ m χ ≪ m Z ′ . We show that for a fixed value of $$m_{Z^\prime }$$ m Z ′ , $$g_{BL}$$ g BL values roughly scale as 1/ Q to reproduce the observed DM abundance. For various Q values in the range between $$10^{-6}$$ 10 - 6 and $$10^2$$ 10 2 , we show that the gauge coupling values $$g_{BL}$$ g BL needed to reproduce the observed DM abundance lie in the search reach of future planned and/or proposed experiments such as FASER, Belle-II, LDMX, and SHiP. In the freeze-in case, the Q values to realize observable $$g_{BL}$$ g BL values are found to be much smaller than that in the freeze-out case.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Explicit solvent machine-learned coarse-grained model of sodium polystyrene sulfonate to capture polymer structure and dynamics

In this article, strongly charged polyelectrolytes (PEs) demonstrate complex solution behavior as a function of chain length, concentrations, and ionic strength. The viscosity behavior is important to understand and is a core quantity for many applications, but aspects remain a challenge. Molecular dynamics simulations using implicit solvent coarse-grained (CG) models successfully reproduce structure, but are often inappropriate for calculating viscosities. To address the need for CG models which reproduce viscoelastic properties of one of the most studied PEs, sodium polystyrene sulfonate (NaPSS), we report our recent efforts in using Bayesian optimization to develop CG models of NaPSS which capture both polymer structure and dynamics in aqueous solutions with explicit solvent. We demonstrate that our explicit solvent CG NaPSS model with the ML-BOP water model [Chan et al. Nat Commun 10, 379 (2019)] quantitatively reproduces NaPSS chain statistics and solution structure. The new explicit solvent CG model is benchmarked against diffusivities from atomistic simulations and experimental specific viscosities for short chains. We also show that our Bayesian-optimized CG model is transferable to larger chain lengths across a range of concentrations. Overall, this work provides a machine-learned model to probe the structural, dynamic, and rheological properties of polyelectrolytes such as NaPSS and aids in the design of novel, strongly charged polymers with tunable structural and viscoelastic properties

77 NANOSCIENCE AND NANOTECHNOLOGY↗