Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic polymer”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Optimal Binding Affinity for Sieving Separation of Propylene from Propane in an Oxyfluoride Anion-Based Metal–Organic Framework

Highly efficient adsorptive separation of propylene from propane offers an ideal alternative method to replace the energy-intensive cryogenic distillation technology. Molecular sieving-type separation via high-performance adsorbents is targeted for superior selectivity, but the limit in adsorption capacity remains a great challenge. Here, we report an oxyfluoride-based ultramicroporous metal–organic framework UTSA-400, [Ni(WO 2 F 4 )(pyz) 2 ] (pyz = pyrazine), featuring one-dimensional pore channels that can accommodate the propylene molecules with optimal binding affinity while specifically excluding the propane molecules. The exposed oxide/fluoride pairs in UTSA-400 serve as strong functional sites for strengthened propylene–host interactions, accounting for a significantly enhanced propylene uptake, while the propane molecules are excluded due to the regulated host framework dynamics. The strong propylene binding enables near-saturation of propylene in the pore confinement at ambient conditions, leading to full utilization of pore space and superior packing density. Combined in situ infrared spectroscopy measurements and dispersion-corrected density functional theory calculations clearly unveil the nature of boosted host–guest binding. Direct production of polymer-grade (>99.5%) propylene with remarkable dynamic productivity is demonstrated by column breakthrough experiments. Furthermore, this work presents an example of pore engineering with atomic precision to break the trade-off in adsorptive separation through guest binding optimization.

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

Improving nanoparticle superlattice stability with deformable polymer gels

The self-assembly of colloidal nanoparticles into ordered superlattices typically uses dynamic interactions to govern particle crystallization, as these non-permanent bonds prevent the formation of kinetically trapped, disordered aggregates. However, while the use of reversible bonding is critical in the formation of highly ordered particle arrangements, dynamic interactions also inherently make the structures more prone to disassembly or disruption when subjected to different environmental stimuli. Thus, there is typically a trade-off between the ability to initially form an ordered colloidal material and the ability of that material to retain its order under different conditions. Here, we present a method for embedding colloidal nanoparticle superlattices into a polymer gel matrix. This encapsulation strategy physically prevents the nanoparticles from dissociating upon heating, drying, or the introduction of chemicals that would normally disrupt the lattice. However, the use of a gel as the embedding medium still permits further modification of the colloidal nanoparticle lattice by introducing stimuli that deform the gel network (as this deformation in turn alters the nanoparticle lattice structure in a predictable manner). Moreover, encapsulation of the lattice within a gel permits further stabilization into fully solid materials by removing the solvent from the gel or by replacing the solvent with a liquid monomer that can be photopolymerized. This embedding method therefore makes it possible to incorporate ordered colloidal arrays into a polymer matrix as either dynamic or static structures, expanding their potential for use in responsive materials.

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

Characterizing laser-heated polymer foams with simultaneous x-ray fluorescence spectroscopy and Thomson scattering at the Matter in Extreme Conditions Endstation at LCLS

Understanding the behavior of polymer foams at high energy density conditions is crucial to advance inertial fusion energy research. Here, we present a new experimental platform designed to measure the thermodynamic state of these materials at megabar pressures. At the Matter in Extreme Conditions Endstation of the Linac Coherent Light Source, we heat samples using an optical, high-intensity, femtosecond laser and dynamically probe them with ultra-short, coherent x-ray pulses of high peak brightness. We perform x-ray Thomson scattering measurements in forward and backward scattering geometries to capture both collective and non-collective electron behavior in the sample. Simultaneously, x-ray fluorescence spectroscopy is used to measure the emission from a mid-Z dopant, providing complementary information on the plasma conditions. By combining these techniques, we obtain temporally resolved temperature measurements of the transient warm dense matter states. Our initial experiments designed to benchmark the platform with carbon samples yielded data resolving the ultrafast response to laser heating with sub-picosecond resolution, measuring plasma temperatures exceeding 50 eV. These findings lay the foundation for precision studies of the dynamic evolution of laser-heated polymer foams.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Effect of surface properties and polymer chain length on polymer adsorption in solution

In polymer nanoparticle composites (PNCs) with attractive interactions between nanoparticles (NPs) and polymers, a bound layer of the polymer forms on the NP surface, with significant effects on the macroscopic properties of the PNCs. The adsorption and wetting behaviors of polymer solutions in the presence of a solid surface are critical to the fabrication process of PNCs. In this study, we use both classical density functional theory (cDFT) and molecular dynamics (MD) simulations to study dilute and semi-dilute solutions of short polymer chains near a solid surface. Using cDFT, we calculate the equilibrium properties of polymer solutions near a flat surface while varying the solvent quality, surface–fluid interactions, and the polymer chain lengths to investigate their effects on the polymer adsorption and wetting transitions. Using MD simulations, we simulate polymer solutions near solid surfaces with three different curvatures (a flat surface and NPs with two radii) to study the static conformation of the polymer bound layer near the surface and the dynamic chain adsorption process. Additionally, we find that the bulk polymer concentration at which the wetting transition in the poor solvent system occurs is not affected by the difference in surface–fluid interactions; however, a threshold value of surface–fluid interaction is needed to observe the wetting transition. We also find that with good solvent, increasing the chain length or the difference in the surface–polymer interaction relative to the surface–solvent interaction increases the surface coverage of polymer segments and independent chains for all surface curvatures. Finally, we demonstrate that the polymer segmental adsorption times are heavily influenced only by the surface–fluid interactions, although polymers desorb more quickly from highly curved surfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1983

Dynamic fracture toughness of glass

Experimentally determined dynamic crack propagation histories in wedge-loaded, modified tapered, and rectangular double cantilever beam specimens were used to drive a dynamic finite element code in its generation phase. The resultant dynamic fracture toughness versus crack velocity relation, during the initial crack acceleration phase of these dynamic fracture specimens, was erratic but followed the standard Gamma-shaped curves of brittle polymers and metals during subsequent crack propagation at terminal velocity and crack deceleration phases. The distinct initiation phase of dynamic crack propagation, which was not observed in dynamic fracture of brittle polymer and metal specimens, is attributed to the lower stored energy in the glass specimen.

Kobayashi, A. S.↗
DOE OSTI2024

Unraveling Morphology and Chemistry Dynamics in Fluoroethylene Carbonate Generated Silicon Anode Solid Electrolyte Interphase Across Delithiated and Lithiated States: Relative Cycling Stability Enabled by an Elastomeric Polymer Matrix

The silicon solid electrolyte interphase (SEI) faces cyclical cracking and reconstruction due to the ∼350% volume expansion. Understanding the SEI dynamic morphology and chemistry evolution from delithiated to lithiated states is thereby paramount to engineering a stable Si anode. Fluoroethylene carbonate (FEC) is a preferred additive with widely demonstrated enhancement of the Si cycling. Thus, insights into the dynamics of the FEC-SEI may provide hints toward engineering the Si interface. Herein, complementary ATR-FTIR, AFM, tip IR, and XPS probing reveal the presence of an elastomeric polycarbonate-like matrix in the FEC-generated SEI which is absent from the FEC-free SEI. Adding FEC to the baseline 1 M LiPF 6 in EC:EMC (1:1) electrolyte promotes formation of a thinner and more conformal SEI, and subdues morphology and chemistry changes between consecutive half-cycles. From AFM, morphological stabilization of the FEC-SEI occurs earlier. Furthermore, conventional SEI biproducts such as Li 2 CO 3 and LiEDC appear in reduced quantities in the FEC-SEI implying a reduced quantity of Li-consuming species. The thin polymeric FEC-SEI enables deeper (de)lithiation of silicon. In conclusion, the enhanced mechanical compliance, chemical invariance, and reduced Li inventory consumption of the FEC-SEI are logically the key features underlying the Si cycling enhancement by FEC.

25 ENERGY STORAGE↗
NASA NTRS2022

Building and Breaking Carbon Composites with REACTER

Carbon-based composites have become indispensable materials in aerospace and other high-performance applications, yet obtaining a detailed, nanoscale understanding of their morphology and failure mechanisms using only experimental methods remains a difficult challenge. REACTER is a versatile computational modeling tool for atomistic molecular dynamics simulations designed to model chemical reactions at the speed and length scales of classical force fields. In this work, several recent features of REACTER were applied to the creation and subsequent mechanical testing of two classes of carbon composites: carbon nanotube (CNT) composites and carbon fiber reinforced polymers (CFRP). A network of CNTs was grown dynamically using the new ‘create atoms’ feature of REACTER. The CNT filler was embedded into a polyarylacetylene (PAA) matrix by simulated in situ polymerization to obtain the final composite model. To generate the second class of carbon composite, fully carbonized (graphitic) carbon fiber morphologies were created by the method of Desai et al. [1], but using the advanced reaction constraints framework of REACTER. Two fiber models were created, representing a circular carbon fiber core and a flat surface, and similarly infiltrated with resin to obtain the final CFRP structure. Failure mechanisms were elucidated by simulating mechanically induced bond breaking, as characterized by third order DFT-based tight-binding simulations, via a reaction constraint on the total potential energy of the involved atoms.

Polymer↗
DOE OSTI2020

Mechanisms of Asymmetric Membrane Formation in Nonsolvent-Induced Phase Separation

We report the first simulations of nonsolvent-induced phase separation (NIPS) that predict membrane microstructures with graded asymmetric pore size distribution. In NIPS, a polymer solution film is immersed in a nonsolvent bath, enriching the film in nonsolvent, and leading to phase separation that forms a solid polymer-rich membrane matrix and polymer-poor membrane pores. We demonstrate how mass-transfer-induced spinodal decomposition, thermal fluctuations, and glass-transition dynamics implemented with mobility contrast between the polymer-rich and polymer-poor phases are essential to the formation of asymmetric membrane microstructures. Specifically, we show that the competition between the propagation of the phase-separation and glass-transition fronts determines the degree of pore-size asymmetry. In conclusion, we also explore the sensitivity of these microstructures to the initial film composition, and compare their formation in 2D and 3D.

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

Uncovering Backbone Conformation for Rigid DPP-Based Donor–Acceptor Conjugated Polymer Using Deuterium Labeling and Neutron Scattering

The conjugated polymer’s backbone conformation dictates the delocalization of electrons, ultimately affecting its optoelectronic properties. Most conjugated polymers can be viewed as semirigid rods with their backbone embedded among long alkyl side chains. Thus, it is challenging to experimentally quantify the conformation of a conjugated backbone. Here, we performed contrast variation neutron scattering on rigid conjugated donor–acceptor (D–A) diketopyrrolopyrrole (DPP) polymers with selectively deuterated side chains to measure the conjugated backbone conformation. We first synthesized DPPbased polymers with deuterated side chains, confirmed by NMR and FTIR. Using contrast variation neutron scattering, we found that the DPP-based conjugated polymers are much more rigid than poly(3-alkylthiophenes), with persistence length (L p ) at 16–18 nm versus 2–3 nm. More importantly, in contrast to the relatively flexible poly(3-alkylthiophenes) whose backbone is more flexible than the whole polymer, we found that the backbone of DPP-based polymers has the same L p value compared to the whole polymer chain. This indicates that side chain interference on backbone conformation is not present for the semirigid polymer, which is further confirmed by coarse-grained molecular dynamics (CG-MD) simulations. Our work provides a novel protocol to probe polymer’s backbone conformation and paradigm-shifting understanding of the backbone conformation of semirigid conjugated polymers.

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

Command of active and responsive elastomers by topological defects and patterns

The project resulted in the development of stimuli-responsive elastomer coatings formed by photopolymerized liquid crystal molecules. The molecular orientation of the liquid crystal elastomers is coupled to rubber-like elasticity. The project developed an approach to produce complex patterns of molecular orientation in elastomers by employing photopatterning technique based on plasmonic metamasks that convert unpolarized incoming light into a transmitted light beam with spatially-varying linear polarization. The polarization-modulated light beam aligns the film of azodye molecules, which serves as the substrate for liquid crystal elastomer coating. The alignment of the azodye molecules follows the polarization pattern with a phase shift of the alignment direction by 90 degrees. The azodye substrate imposes orientation of molecules in the adjacent liquid crystal layer which is preserved after the monomers are polymerized into a liquid crystal elastomer. By using differently aligned substrates one could produce various patterns of molecular orientation in the elastomer, which dictate the mechanical properties of the elastomers and their response to external stimuli such as temperature, humidity, and ultraviolet irradiation. In particular, the predesigned molecular orientation produces deterministically defined topography response of the liquid crystal coatings, in which their free surface changes from flat to either locally elevated or locally compressed, depending on whether the pre-inscribed molecular orientation is predesigned with a splay or bend deformation. The research established that the mechanism of this deterministic relation between the dynamic topography and pre-inscribed molecular orientation is rooted in the change of the tensorial order parameter that characterizes the degree of orientational order, and in the contraction/expansion of polymer molecules in response to the order parameter changes. In particular, lowering the order causes contraction of polymer globules. The work demonstrated that the dynamic topography of elastomer coatings with three-dimensional pattern of orientational order includes up and down motion of the material (along the normal to the coating) and also shifts in lateral directions. We demonstrate that the deterministic relationship between the orientational patterns and variations of coatings’ profile is caused by forces that are mathematically equivalent to active forces in the system of “pullers” or “pushers” in active matter, demonstrating universality of the out-of-equilibrium behavior. The project resulted in theoretical models that predictively describe the dynamic response of coatings. The topography of liquid crystal elastomers is sensitive to stimuli such as temperature, ultraviolet irradiation, and humidity. The applicability of the patterned approach to produce dynamic coatings was expanded from nematic elastomers to their smectic counterparts. The project combined a battery of imaging techniques to unveil the dynamic properties of the coatings, ranging from the state-of-the-art dynamic holographic microscopy to fluorescent confocal polarizing microscopy and polscope microscopy. The research uncovered potential applications of dynamic elastomer coatings, such as control placement of colloidal particles, engineering of living tissues with alignment patterns of cells that follow the orientational order of elastomer coatings. The approach developed in the project can be used to design programmable dynamic coatings with functionalities that mimic biological tissues such as skin.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Molecular Alignment of a Meta-Aramid on Carbon Nanotubes by In Situ Interfacial Polymerization

Molecularly organized nanocomposites of polymers and carbon nanotubes (CNTs) have great promise as high-performance materials; in particular, conformal deposition of polymers can control interfacial properties for mechanical load transfer, electrical or thermal transport, or electro/chemical transduction. However, controllability of polymer-CNT interaction remains a challenge with common processing methods that combine CNTs and polymers in melt or in solution, often leading to nonuniform polymer distribution and CNT aggregation. In this work, we demonstrate CNTs within net-shape sheets can be controllably coated with a conformal coating of meta-aramid by simultaneous capillary infiltration and interfacial polymerization. We determine that π-interaction between the polymer and CNTs results in chain alignment parallel to the CNT outer wall. Subsequent nucleation and growth of the precipitated aramid forms a smooth continuous layered sheath around the CNTs. These findings motivate future investigation of mechanical properties of the resulting composites, and adaptation of the in situ polymerization method to other substrates.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

Late Stage Domain Coarsening Dynamics of Lamellar Block Copolymers

Block copolymer (BCP) phase separation dynamics can be expected to differ significantly from that of the polymer blends due to the constraint of chain connectivity. BCP phase separation dynamics has been studied theoretically, but there has been little experimental evidence to confirm the BCP domain growth scaling laws put forward theoretically. Here, we investigate the dynamics of late-stage lamellar BCP domain coarsening and find that the scaling exponent for lamellar domain growth is ≈1/6 (0.17), irrespective of the annealing temperature, a value close to the scaling exponent of 0.2 predicted by theoretical studies. Furthermore, we show that the prefactors in the domain coarsening equation show Arrhenius dependence on temperature, indicating that the BCP domain growth dynamics is Arrhenius over the temperature range investigated.

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

Impact of Ion–Ion Correlated Motion on Salt Transport in Solvated Ion Exchange Membranes

The influence of dynamical ion-ion correlations and ion pairing on salt transport in ion exchange membranes remain poorly understood. In this study, we use the framework of Onsager transport coefficients within atomistic molecular dynamics sim- ulations to study the impact of ion-ion correlated motion on salt transport in hydrated polystyrene sulfonate membranes and compare with the results from aqueous salt so- lutions. At sufficiently high salt concentrations, cation-anion dynamical correlations exert a significant influence on both salt diffusivities and conductivities. Anion-anion distinct correlations, arising from the imbalance between the concentration of free (mo- bile) cations and anions, and the retarding effect of the fixed charge groups on cations, proves to be an additional important feature for polymer membranes. Furthermore, our results demonstrate that dynamical correlations should become an important consideration in experimental measurements of salt diffusivities and conductivities for non-dilute salt solutions in polymer membranes.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

Machine-Learned Committor Functions for Reactive Molecular Dynamics

Reactive molecular dynamics (MD) is a powerful tool for atomistic-scale modeling of a diverse range of chemical processes. However, scaling these simulations to large systems and long times scales remains a challenge because of the complexity of the potential energy function required. The authors previously developed a heuristic approach, called REACTER, that incorporates reactivity in MD simulations in a less general but much more computationally efficient manner. REACTER uses standard, fixed valence force fields as the underlying potentialenergy surface for describing all interatomic interactions but adds a procedure for enforcing user-defined reactions that occur when certain geometric constraints on relative atomic positions are satisfied. Further, these bonding changes can be accepted or rejected with a probability related tothe local thermal energy. This work seeks to generalize this approach by replacing the set of user defined geometric constraints and energetic criteria with a committor function that specifies the probability of a reaction occurring on the basis of the local atomic configuration. The committor function is a useful mathematical tool for modeling rare events but, unfortunately, is very difficult to compute for realistic systems in a general way. This work describes a method for approximating the committor function using a machine learning approach, specifically a deep neural network trained with data from reactive MD and DFT-based dynamics simulations. This network is coupled to the existing REACTER protocol, as implemented in the LAMMPS MD package, and used to make on-the-fly predictions of reaction probabilities without the more extensive user input previously required. The new method is demonstrated using the polymerization of polystyrene as a case study. Although very dependent on the quality and quantity of training data, machine-learned committor functions show promise as a method for incorporating reaction probability from higher level calculations into highly scalable MD simulations.

polymer simulations↗
DOE OSTI2021

Internal Microstructure Dictates Interactions of Polymer-grafted Nanoparticles in Solution

Understanding the effects of polymer brush architecture on particle interactions in solution is requisite to enable the development of functional materials based on self-assembled polymer-grafted nanoparticles (GNPs). Static and dynamic light scattering of polystyrene-grafted silica particle solutions in toluene reveals that the pair interaction potential, inferred from the second virial coefficient, A 2 , is strongly affected by the grafting density, σ, and degree of polymerization, N, of tethered chains. In the limit of intermediate σ (~0.3 to 0.6 nm -2 ) and high N, A 2 is positive and increases with N. This confirms the good solvent conditions and can be qualitatively rationalized on the basis of a pair interaction potential derived for grafted (brush) particles. In contrast, for high σ > 0.6 nm -2 and low N, A 2 displays an unexpected reversal to negative values, thus indicating poor solvent conditions. These findings are rationalized by means of a simple analysis based on a coarse-grained brush potential, which balances the attractive core-core interactions and the excluded volume interactions imparted by the polymer grafts. The results suggest that the steric crowding of polymer ligands in dense GNP systems may fundamentally alter the interactions between brush particles in solution and highlight the crucial role of architecture (internal microstructure) on the behavior of hybrid materials. The effect of grafting density also illustrates the opportunity to tailor the physical properties of hybrid materials by altering geometry (or architecture) rather than a variation of the chemical composition.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Water Dynamics of Superacid Aromatic Proton Exchange Membranes for Fuel Cell Applications

Proton exchange membranes (PEMs) with high conductivity are of critical importance for the development of fuel cells, electrolyzers, and other electrochemical technologies. In this research, poly(1,1,2,2-tetrafluoro-2-phenoxyethane-1-sulfonic acid) (PTPS) with an aromatic polymer main chain and a perfluorinated superacidic polymer side chain was synthesized. The water dynamics of PTPS were characterized across various length scales using a combination of Fourier-transform infrared spectroscopy (FTIR) and nuclear magnetic resonance (NMR) and compared with Nafion, a standard perfluorinated PEM, and sulfonated poly(ether sulfone) (SPES 40), an aromatic PEM without perfluorinated superacid side chains. The T 1 and T 2 relaxation times of water in the samples probed by NMR increase from SPES 40 to PTPS to Nafion, indicating that the local motion of the water molecules becomes faster. This trend corresponds well with the relative fraction of bulk-like water determined using FTIR. At larger length scales, the diffusion coefficient of water was characterized using pulsed-field gradient NMR (PFG-NMR). At a longer diffusion time (Δ = 100 ms), PTPS has a smaller diffusion coefficient compared with both Nafion and SPES 40, due to restricted diffusion, and this effect is also evident in the proton conductivity of the hydrated membranes. From this comparison, it is apparent that the aromatic backbone and side chain type greatly influence the water dynamics in PEMs at various length scales and the water dynamics significantly impact the bulk proton conductivity. These insights will lead to new designs for aromatic PEMs and help to identify bottlenecks in current materials.

25 ENERGY STORAGE↗
NASA NTRSDate not supplied

REACTER: A Heuristic Method for Reactive Molecular Dynamics

REACTER (www.reacter.org) is a heuristic protocol that enables the simulation of complex reactions using atomistic molecular dynamics (MD) with a fixed-valence force field. Incorporating reactions into classical MD with this approach allows modeling of reactive systems over greatly-increased time scales, enabling systems to be modeled with MD that would not otherwise be feasible. One or more competing multi-step reactions or series of reactions can be invoked simultaneously. Special treatment can be applied to neighboring atoms to relax high energy configurations while the simulation progresses. The original version of REACTER, which was implemented in the open-source LAMMPS simulation package as fix bond/react, was only available for serial simulations. This work describes the expansion of fix bond/react for use in parallel simulations, as well as the addition of various new options, including deletion of reaction by-products, reversible reactions, and custom reaction constraints. These new capabilities are demonstrated through large-scale simulations (200,000+ atoms) of the polymerization of polystyrene and nylon 6,6. The morphologies of both polymers are analyzed after reaching >99% extent of polymerization. Finally, the newly-added reversible reactions feature is demonstrated by rupturing these highly-entangled systems under uniaxial strain by defining a chain scission reaction.

polymer simulations, molecular dynamics↗
NASA NTRSDate not supplied

REACTER: A Versatile Tool for Large-Scale Reactive Molecular Dynamics

Accurately describing reactive events over long length and time scales remains a grand challenge of computational materials science. REACTER is a general protocol for modeling chemical reactions using classical force fields, and is implemented in the popular molecular dynamics software LAMMPS. REACTER has a growing user base and has been used as a model-building tool for a variety of materials, including thermoplastics, thermosets, glassy materials and composites. The method has also been applied to accelerated modeling of reversible chemical reactions, such as the formation of electrochemical components for batteries. Recently, the REACTER protocol has received some major upgrades to enhance its ability to predict when reactions occur and to make it easier to use. Force field parameters can now be automatically assigned to newly created bonds, angles and other interactions. Advanced reaction constraints have also been added, including an Arrhenius constraint to enforce an effective activation energy, a root-mean-square-deviation option for complex geometrical constraints, and a constraint based on the total potential energy of the atoms involved in a reactive site. This potential energy constraint allows for the accurate reproduction of DFT-based tight-binding (DFTB3) predicted bond dissociation curves for mechanically induced bond breaking.

polymer simulations, molecular dynamics↗