Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “kinetic parameters”

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

Kinetics of Carbon Condensation in Detonation of High Explosives: First-Order Phase Transition Theory Perspective

We report the kinetics of carbon condensation, or carbon clustering, in detonation of carbon-rich high explosives is modeled by solving a system of rate equations for concentrations of carbon particles. Unlike previous efforts, the rate equations account not only for the aggregation of particles but also for their fragmentation in a thermodynamically consistent manner. Numerical simulations are performed, yielding the distribution of particle concentrations as a function of time. In addition to that, analytical expressions are obtained for all the distinct steps and regimes of the condensation kinetics, which facilitates the analysis of the numerical results and allows one to study the sensitivity of the kinetic behavior to the variation of system parameters. The latter is important because the numerical values of many parameters are not reliably known at present. The theory of the kinetics of first-order phase transitions is found adequate to describe the general kinetic trends of carbon condensation, as described by the rate equations. Such physical phenomena and processes as the coagulation, nucleation, growth, and Ostwald ripening are observed, and their dependence on various system parameters is studied and reported. It is believed that the present work will become useful when analyzing the present and future results for the kinetics of carbon condensation, obtained from experiments or atomistic simulations.

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

Ultrathin Plasmonic Tungsten Oxide Quantum Wells with Controllable Free Carrier Densities

Localized surface plasmon resonances (LSPR) of nanostructures can be tuned by controlling their morphology, local dielectric environment, and free carrier concentration. We report the colloidal synthesis of an ~3 tungsten–oxygen (W-O) layer thick (~1 nm), two-dimensional (2D) WO 3-x nanoplatelets (NPLs) (x ≈ 0.55–1.03), which display tunable near-infrared LSPR properties and additionally high free electron density (N e ) that arises predominantly from the large shape factor of 2D NPLs. Importantly, the W to O composition ratios inferred from their LSPR measurements show much higher percentage of oxygen vacancies than those determined by X-ray diffraction analysis, suggesting that the aspect ratio of ultrathin WO 3-x NPLs is the key to producing an unprecedentedly large Ne, although synthesis temperature is also an independent factor. We find that NPL formation is kinetically controlled, whereas thermodynamic parameter manipulation leads to N e values as high as 4.13 × 10 22 cm –3 , which is close to that of plasmonic noble metals, and thus our oxide-based nanostructures can be considered as quasi-metallic. The unique structural properties of 2D nanomaterials along with the high N e of WO 3-x NPLs provide an attractive alternative to plasmonic noble metal nanostructures for various plasmon-driven energy conversions and design of photochromic nanodevices.

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

A DFT-based kinetic Monte Carlo simulation of multiphase oxide-metal thin film growth

Functional thin films of nanoscale metal pillars in oxide or nitride matrices known as vertically aligned nanocomposite (VAN) have gained much interest owing to their unique strain-coupled and highly anisotropic properties. So far, the deposition of these films has been explored mostly experimentally. In this work, a density functional theory (DFT)-based kinetic Monte Carlo simulation model using Bortz–Kalos–Lebowitz algorithm was developed to understand the growth of VAN films deposited by pulsed laser technique on mismatching substrates. The model has been parameterized and applied to understand the kinetics of growth thin films consisting of Au pillars in CeO2 matrix deposited on SrTiO3 substrates. The effects of pulsed laser deposition (PLD) conditions including the pulse frequency, deposition flux, and substrate temperature were explored. The simulations indicate that the Au pillar size and shape exhibit significant dependence on the PLD conditions. Namely, increasing the temperature increases the average pillar size and lowers the pillar density, and vice versa. In addition, the simulations revealed that increasing the deposition rate results in lowering the average pillar size and increasing the density. Particularly, the DFT results suggest that Au pillar size can be tuned during the initial growth of the first monolayer due to the significantly low activation barrier. Our analysis showed that the relationship between the average pillar size and pillar density is influenced by the kinetics. Furthermore, autocorrelation analysis showed that pillars self-organize in quasi-ordered patterns at certain windows of the deposition conditions, which is attributed to the complex nature of the chemical interactions in the system, the kinetics, and the deposition parameters.

Physics↗
DOE OSTI2026

Identified charged hadron production in Au+Au collisions at $\sqrt{s_{NN}}$ = 54.4 GeV with the STAR detector

Here, we present results on the production of 𝜋±, 𝐾±, 𝑝, and $\bar{𝑝}$ in Au+Au collisions at $\sqrt{𝑠_{N⁢N}}$ = 54.4~GeV using the STAR detector at RHIC, at mid-rapidity (|𝑦|< 0.1). Invariant yields of these particles as a function of transverse momentum are shown. We determine bulk properties such as integrated particle yields (𝑑 ⁢𝑁 /𝑑 ⁢𝑦 ), mean transverse momentum (⟨𝑝 𝑇 ⟩), particle ratios, which provide insight into the particle production mechanisms. Additionally, the kinetic freeze-out parameters (𝑇 kin and ⟨𝛽 𝑇 ⟩), which provide information about the dynamics of the system at the time of freeze-out, are obtained. The Bjorken energy density (𝜖 𝐵⁢𝐽 ), which gives an estimate of the energy density in the central rapidity region of the collision zone at the formation time 𝜏, is calculated and presented as a function of multiplicity for various energies. The results are compared with those from the models such as A Multi-Phase Transport (AMPT) and Heavy Ion Jet INteraction Generator (HIJING) for further insights.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Pit Stability Predictions of Additively Manufactured SS316 Surfaces Using Finite Element Analysis

Stainless steels are susceptible to localized forms of corrosion attack, such as pitting. The size and lifetime of a nucleated pit can vary, depending on a critical potential or current density criterion, which determines if the pit repassivates or continues growing. This work uses finite element method (FEM) modeling to compare the critical pit radii predicted by thermodynamic and kinetic repassivation criteria. Experimental electrochemical boundary conditions are used to capture the active pit kinetics. Geometric and environmental parameters, such as the pit shape and size (analogous to additively manufactured lack-of-fusion pores), solution concentration, and water layer thickness were considered to assess their impact on the pit repassivation criterion. The critical pit radius (the transition point from stable growth to repassivation) predicted for a hemispherical pit was larger when using the repassivation potential (E rp ) criteria, as opposed to the current density criteria (pit stability product). Including both the pit stability product and E rp into its calculations, the analytical maximum pit model predicted a critical radius two times more conservative than the FEA approach, under the conditions studied herein. The complex pits representing lack-of-fusion pores were shown to have minimal impact on the critical radius in atmospheric conditions.

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

Convective and Turbulent Motions in Nonprecipitating Cu. Part III: Characteristics of Turbulence Motions

Velocity field in a nonprecipitating Cu under BOMEX conditions, simulated by SAM with 10-m resolution and spectral bin microphysics is separated into the convective part and the turbulent part, using a wavelet filtering. In Part II of the study properties of convective motions of this Cu were investigated. Here in Part III of the study, the parameters of cloud turbulence are calculated in the cloud updraft zone at different stages of cloud development. The main points of this study are (i) application of a fine-scale LES model of a single convective cloud allowed a direct estimation of turbulence parameters using the resolved flow in the cloud and (ii) the separation of the resolved flow into the turbulence flow and the nonturbulence flow allowed us to estimate different turbulent parameters with sufficient statistical accuracy. We calculated height and time dependences of the main turbulent parameters such as turbulence kinetic energy (TKE), spectra of TKE, dissipation rate, and the turbulent coefficient. It was found that the main source of turbulence in the cloud is buoyancy whose contribution is described by the buoyancy production term (BPT). The shear production term (SPT) increases with height and reaches its maximum near cloud top, and so does BPT. In agreement with the behavior of BPT and SPT, turbulence in the lower cloud part (below the inversion level) is weak and hardly affects the processes of mixing and entrainment. The fact that BPT is larger than SPT determines many properties of cloud turbulence. For instance, the turbulence is nonisotropic, so the vertical component of TKE is substantially larger than the horizontal components. Another consequence of the fact that BPT is larger than STP manifests itself in the finding that the turbulence spectrum largely obeys the -11/5 Bolgiano–Obukhov scaling. The classical Kolmogorov -5/3 scaling dominates for the low part of a cloud largely at the dissolving stage of cloud evolution. Using the spectra obtained we evaluated an “effective” dissipation rate which increases with height from nearly zero at cloud base up to 20 cm 2 s -3 near cloud top. The coefficient of turbulent diffusion was found to increase with height and ranged from 5 m 2 s -1 near cloud base to 25 m 2 s -1 near cloud top. In conclusion, the possible role of turbulence in the process of lateral entrainment and mixing is discussed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Role of copper in nitrous oxide accumulation in wetland soils from Illinois and South Carolina and stream sediments from Tennessee

Denitrification is microbially-mediated through enzymes containing metal cofactors. Laboratory studies of pure cultures have highlighted that the availability of copper (Cu), required for the multicopper enzyme nitrous oxide reductase, can limit nitrous oxide (N2O) reduction. However, in natural aquatic systems, such as wetlands and hyporheic zones in stream beds, the role of Cu in controlling denitrification remains incompletely understood. In this study, we collected soils and sediments from three natural environments -- riparian wetlands, marsh wetlands, and a stream -- to investigate their nitrogen species transformation activity at background Cu levels and different supplemented Cu loadings. All of the systems contained solid-phase associated Cu below or around geological levels (40–280 nmol g-1) and exhibited low dissolved Cu (3–50 nM), which made them appropriate sites for evaluating the effect of limited Cu availability on denitrification.The dataset contains the variation in the nutrient (nitrate (NO3-), nitrite (NO2-), and ammonium (NH4+)) and N2O concentrations during the incubation experiments to evaluate the effect of Cu. We have also reported the variation in metal concentrations (Cu, iron(Fe) and manganese (Mn)) and dissolved organic carbon (DOC) concentrations during the incubations. The parameters obtained using a kinetic model to quantitatively report the effect of Cu on nitrogen species conversion in natural aquatic systems are also added in the dataset. The dataset also contains the labile concentration of Cu estimated using a speciation model in MINTEQ. All data are provided in text-based CSV format with header sections indicating the data contained in each file and the corresponding units.Our study suggests that high concentrations of N2O accumulated in all microcosms lacking Cu amendment except for one stream sediment sample. With Cu added to provide dissolved concentrations at trace levels (10–300 nM), the reduction rate of N2O to N2 in the wetland soils and stream sediments was enhanced. A kinetic model could account for the trends in nitrogen species by combining the reactions for microbial reduction of NO3- to NO2-/N2O/N2 and abiotic reduction of NO2 to nitrogen (N2_. The model revealed that the rate of N2O to N2 conversion increased significantly in the presence of Cu. For riparian wetland soils and stream sediments, the kinetic model also suggested that overall denitrification is driven by abiotic reduction of NO2- in the presence of inorganic electron donors. This study demonstrated that natural aquatic systems containing Cu at concentrations less than or equal to crustal abundances may display incomplete reduction of N2O to N2 that would cause N2O accumulation and release to the atmosphere.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Determining the Crosslinking and Degradation Reaction Kinetics in Photovoltaic Encapsulants Using Accelerated Aging: Preprint

Degradation of photovoltaic (PV) module encapsulant mechanical characteristics that lead to embrittlement and delamination remains a cause of failure in solar installations. A multiscale reliability model based on detailed molecular degradation reaction kinetics was previously published, connecting the encapsulant mechanical properties (elastic modulus, yield strength, and adhesion energy) to the degraded molecular structure and interfacial bond density to adjacent solar cell and glass substrates. The model, developed primarily for poly(ethylene-co-vinyl acetate) (EVA) encapsulants, remains to be experimentally validated. Determining the degradation and crosslinking kinetics of alternative encapsulants, such as polyolefin elastomer (POE) and EVA/POE/EVA composites (EPE), can further generalize the model. Importantly, the activation energy for crosslinking of fully cured PV encapsulant products from temperature or UV is presently unknown or unavailable. In this work, we subject EVA, POE, and EPE encapsulants to a set of accelerated aging conditions, varying the temperature, UV intensity, and relative humidity. We use DSC, FTIR, and Soxhlet extraction (gel content) to characterize the encapsulants' changing molecular structure. This allows for determination of the photochemical degradation and crosslinking kinetics of the encapsulants. Preliminary results show an increase in gel content (crosslinking) and a decrease in crystallinity of EVA, POE, and EPE encapsulants under hot-aerobic (90 deg C, 22% RH) and hot-anaerobic (90oC, sealed in N2 air) accelerated aging, even in the absence of UV and crosslinking initiating agents. For a least squares regression with an assumed first-order crosslinking kinetics, the crosslinking rate parameters were computed for the three encapsulants under hot-aerobic and hot-anaerobic aging conditions. FTIR analysis showed insignificant encapsulant degradation for hot-aerobic and hot-anaerobic aging, demonstrating the critical role of UV and moisture in accelerating degradation. Aging conditions with UV exposure and elevated humidity are presently in progress.

adhesion↗
DOE OSTI2026

Design and Commissioning of a Deuterium-Tritium Gas Delivery System for Muon Catalyzed Fusion in a Diamond Anvil Cell

We report the design, commissioning, and operation of deuterium-deuterium (DD) and deuterium-tritium (DT) gas delivery systems developed to load a diamond anvil cell (DAC) beam target for muon-catalyzed fusion (muCF). The DAC approach enables DT fuel to be compressed to GPa pressures at more than twice the liquid density and heated from cryogenic temperatures through 500 K, opening access to a substantially expanded parameter range for muCF kinetics and yield measurements. In this approach, DT is cryo-condensed to a liquid in a minichamber and then compressed in the DAC using a helium-driven pneumatic membrane, achieving high pressures in a millimeter-scale DT sample volume. A DD gas delivery system was designed and used to validate the experimental apparatus, measure the gas quantities needed for filling, develop operational experience, and collect kinetics and yield data with DD targets. The DT gas delivery system adds tritium-specific capabilities for inventory minimization, secondary containment, and activity monitoring. The DT system integrates depleted uranium storage beds and a liquid helium cryogenic condenser used for pressure building and cryopumping. High-purity delivery is provided by a rapid-response palladium permeator. The system is housed in a helium-atmosphere glovebox held at negative pressure with continuous cleanup. We present the process and instrumentation design, a failure modes and effects analysis (FMEA), and data from the experiment's in situ Raman spectrometer, which provides direct confirmation of target loading and composition through the optically clear diamond anvils. The 2024 and 2025 DT campaigns achieved repeatable target fills and operation with no measurable tritium releases to the stack, demonstrating safe, high-purity DT loading at novel density-temperature conditions for muCF studies.

Koukina, Elena [Acceleron Fusion]↗
DOE OSTI2024

In situ thermal conductivity measurement revealing kinetics of thermochemical reactions

Utilizing thermochemical reactions for thermal energy storage and solar fuel production has been an emerging research topic. Thermal transport properties of the materials are an important parameter that can determine the kinetics and efficiency of thermochemical reactions. With the increasing number of new thermochemical materials (TCMs); however, there is a lack of reliable techniques to monitor the thermal transport property of the materials and their changes as a function of reactions in real time. In this work, we report the in situ monitoring of thermochemical reactions using modulated photothermal radiometry (MPR). The thermal conductivities of two TCMs, namely, calcium hydroxide (Ca(OH) 2 ) and Ba 0.15 Sr 0.85 FeO 3–δ (BSF1585), were measured as a function of temperature and time using the MPR technique. The measured thermal conductivities were correlated to the reaction. The work has two significant contributions to the research communities. First, it provides a non-invasive diagnostic tool for monitoring the thermal transport properties of TCMs that can potentially be a high-throughput measurement technique conducive to optimizing TCMs, reactors, and related thermal systems. Second, for TCMs that show observable changes in thermal transport properties, a correlation between the measured thermal conductivity and the conversion fraction of the reaction can be established for monitoring the reaction kinetics based on thermal characterization.

14 SOLAR ENERGY↗
DOE OSTI2021

Pushing and Pulling: A Dual pH Trigger Controlled by Varying the Alkyl Tail Length in Heme Coordinating Peptide Amphiphiles

Some organisms in nature that undergo anaerobic respiration utilize 1D nanoscale arrays of densely packed cytochromes containing the molecule heme. The assemblies can be mimicked with 1D nanoscale fibrils composed of peptide amphiphiles designed to coordinate heme in dense arrays. To create such materials and assemblies it is critical to understand the assembly process and what controls the various aspects of hierarchical assembly. MD simulations suggest that shorter alkyl chains on the peptide lead to more dynamic structures than the peptides with longer chains that yield kinetically trapped states. The hydration parameters manifest themselves experimentally through the observation of a dual pH trigger which controls the peptide assembly rate, the heme binding affinity, and heme organization kinetics. Great strides in understanding the relative complexity of the self-assembly process in relation to incorporating a functional moiety like heme opens up many possibilities in developing abiotic assemblies for bioelectronic devices and assemblies

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

Dark Photons from Perturbative Decay of a Misaligned Higgs Field

We reconsider the production of dark photons $A'$ as dark matter, from the perturbative decay of a dark Higgs field $h$, that is stochastically misaligned from the minimum of its potential during inflation. This is a simple and predictive framework for generating the $A'$ relic abundance. It is constrained by structure formation, since the $A'$ are initially boosted, and inflationary isocurvature fluctuations, which require small quartic couplings $λh^4$. We identify $A'$ masses between 100 eV and 1 GeV and gauge couplings $g\sim 10^{-15}-10^{-10}$ that are consistent in this scenario, and which become more tightly constrained if a generic level of kinetic mixing is present. The favored parameter region could be tested through future CMB or Lyman-$α$ observations, and, in the presence of kinetic mixing, by direct detection experiments or diffuse soft gamma-ray searches.

Cline, James M. [McGill U., Montreal (main)]↗
DOE OSTI2025

A neural master equation framework for multiscale modeling of molecular processes: application to atomic-scale plasma processes

Plasma-surface interactions (PSI) play a crucial role in microelectronics fabrication; however, their multiscale nature and array of complex, often unknown interactions make computational modeling of PSIs extremely difficult. To this end, we propose a general neural master equation (NME) framework that uses master equations to describe the dynamics of a molecular process, wherein neural networks learned from atomistic simulations represent unknown transitions between different system states. By leveraging the physics-based structure of master equations and data-driven state transitions, the NME framework promotes generalizability and physics interpretability, and can bridge disparate length and time scales. The framework is demonstrated for multiscale modeling of Si atomic layer etching and reactive ion etching, where the learned NME-based surface kinetic models exhibit good predictive and extrapolative capabilities for predicting experimentally relevant observables as a function of process parameters. The NME-based surface kinetic models obey physical constraints, which are violated in models based on neural ordinary differential equations. The proposed NME framework for multiscale modeling of molecular processes can pave the way for the discovery of new chemistries and materials in atomic-scale plasma processes.

Chemical engineering↗
DOE OSTI2022

Pathways Toward High-energy Li-sulfur Batteries, Identified via Multi-reaction Chemical Modeling

Here we present a 1D model of a Li-Sulfur battery with physically derived geometric parameters and thermodynamically consistent electrochemical kinetics. The approach enables straightforward comparison of proposed Li-S mechanisms and provides insights into the influence of polysulfide intermediates on battery discharge. Comparing predictions from multiple mechanisms demonstrates the need for both lithiated and non-lithiated polysulfide species, and highlights the challenge of developing parameter estimates for complex electrochemical mechanisms. The model is also used to explore cathode design strategies. Discharge performance and polysulfide concentrations for electrolyte/sulfur ratios in the range 2–4 μ L mg −1 identifies trade-offs that limit battery energy and power density, and highlights the risk of polysulfide precipitation. New cathode and electrolyte approaches must limit polysulfide concentrations in the electrolyte, both to unlock better rate capabilities in Li-S technology and to prevent capacity fade due to polysulfide precipitation.

25 ENERGY STORAGE↗
DOE OSTI2023

Proton transfer kinetics of transition metal hydride complexes and implications for fuel-forming reactions

Proton transfer reactions involving transition metal hydride complexes are prevalent in a number of catalytic fuel-forming reactions, where the proton transfer kinetics to or from the metal center can have significant impacts on the efficiency, selectivity, and stability associated with the catalytic cycle. Here, this review correlates the often slow proton transfer rate constants of transition metal hydride complexes to their electronic and structural descriptors and provides perspective on how to exploit these parameters to control proton transfer kinetics to and from the metal center. A toolbox of techniques for experimental determination of proton transfer rate constants is discussed, and case studies where proton transfer rate constant determination informs fuel-forming reactions are highlighted. Opportunities for extending proton transfer kinetic measurements to additional systems are presented, and the importance of synergizing the thermodynamics and kinetics of proton transfer involving transition metal hydride complexes is emphasized.

08 HYDROGEN↗
DOE OSTI2020

Size Transformation of the Au 22 (SG) 18 Nanocluster and Its Surface-Sensitive Kinetics

For many applications of well-defined gold nanoclusters, it is desirable to understand their structural evolution behavior under working conditions with molecular precision. Here we report the first systematic investigation of the size transformation products of the Au 22 (SG) 18 nanocluster under representative working conditions and highlight the surface effect on the transformation kinetics. Under thermal and aerobic conditions, the consecutive and pH-dependent transformation from Au 22 to both well-defined clusters and small Au(I)SR species was identified by ESI-MS and UV-vis spectroscopy. By introducing a perturbation onto the Au 22 surface, significant changes in the activation parameters were determined from the kinetic study of the Au 22 transformation. This indicates the sensitivity of the nanocluster transformation pathway to the cluster surface. Finally, the systematic study of cluster transformation and the sensitivity of cluster transformation to the surface revealed herein has significant implications for future attempts to design gold nanoparticles with adaptation to the working environment and the regeneration of active nanoparticles.

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

Validation study of RWM stability in DIII-D high- β N plasmas

The n = 1 (n is the toroidal mode number) resistive wall mode (RWM) stability is numerically investigated for two DIII-D high-β N discharges 176440 and 172461, utilizing the MARS-F (Liu et al 2000 Phys. Plasmas 7 3681) and MARS-K (Liu et al 2008 Phys. Plasmas 15 112503) codes. Systematic validation efforts are attempted, for the first time, for discharges with very slow or vanishing toroidal flow for a large fraction of the plasma volume. While gaining physics insights in accessing stable operation regime at β N exceeding the Troyon no-wall limit in these slow-rotation experiments, the predictive capability of fluid and non-perturbative magnetohydrodynamic-kinetic hybrid models for the RWM is further confirmed. The MARS-F fluid model, with a strong but numerically tunable viscosity mimicking ion Landau damping of parallel sound waves, finds complete stabilization of the n = 1 RWM in the considered DIII-D plasmas under the experimental flow conditions. Similarly, either full stabilization (for discharge 176440) or marginal stability (for discharge 172461) of the mode is computed by the MARS-K hybrid model, which is first-principle based without free model parameters. In particular, all drift kinetic resonances, including those of thermal and energetic particles, are found to synergistically act to marginally stabilize the RWM in discharge 172461. These MARS-F/K modeling results explain the experimentally observed stable operational regime in DIII-D, as far as the RWM stability is concerned. Extensive numerical sensitivity studies, with respect to the plasma toroidal flow speed as well as the radial location of the resistive wall, are also carried out to further support the validation study.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Thermodynamic stability versus kinetic accessibility: Pareto fronts for programmable self-assembly

A challenge in designing self-assembling building blocks is to ensure the target state is both thermodynamically stable and kinetically accessible. These two objectives are known to be typically in competition, but it is not known how to simultaneously optimize them. We consider this problem through the lens of multi-objective optimization theory: we develop a genetic algorithm to compute the Pareto fronts characterizing the tradeoff between equilibrium probability and folding rate, for a model system of small polymers of colloids with tunable short-ranged interaction energies. We use a coarse-grained model for the particles' dynamics that allows us to efficiently search over parameters, for systems small enough to be enumerated. For most target states there is a tradeoff when the number of types of particles is small, with medium-weak bonds favouring fast folding, and strong bonds favouring high equilibrium probability. Additionally, the tradeoff disappears when the number of particle types reaches a value m *, that is usually much less than the total number of particles. This general approach of computing Pareto fronts allows one to identify the minimum number of design parameters to avoid a thermodynamic–kinetic tradeoff. However, we argue, by contrasting our coarse-grained model's predictions with those of Brownian dynamics simulations, that particles with short-ranged isotropic interactions should generically have a tradeoff, and avoiding it in larger systems will require orientation-dependent interactions.

36 MATERIALS SCIENCE↗