The explicit bonding reaction ensemble Monte Carlo method

Kód výsledku v IS VaVaI

<a href="https://www.isvavai.cz/riv?ss=detail&h=RIV%2F00216208%3A11310%2F24%3A10485760" target="_blank" >RIV/00216208:11310/24:10485760 - isvavai.cz</a>

Výsledek na webu

<a href="https://verso.is.cuni.cz/pub/verso.fpl?fname=obd_publikace_handle&handle=Yu9s~yj_Be" target="_blank" >https://verso.is.cuni.cz/pub/verso.fpl?fname=obd_publikace_handle&handle=Yu9s~yj_Be</a>

DOI - Digital Object Identifier

<a href="http://dx.doi.org/10.1063/5.0226122" target="_blank" >10.1063/5.0226122</a>

Jazyk výsledku

angličtina

Název v původním jazyce

The explicit bonding reaction ensemble Monte Carlo method

Popis výsledku v původním jazyce

We present the explicit bonding Reaction ensemble Monte Carlo (eb-RxMC) method, designed to sample reversible bonding reactions in macromolecular systems in thermodynamic equilibrium. Our eb-RxMC method is based on the reaction ensemble method; however, its implementation differs from the latter by the representation of the reaction. In the eb-RxMC implementation, we are adding or deleting bonds between existing particles, instead of inserting or deleting particles with different chemical identities. This new implementation makes the eb-RxMC method suitable for simulating the formation of reversible linkages between macromolecules, which would not be feasible with the original implementation. To enable coupling of our eb-RxMC algorithm with molecular dynamics algorithm for the sampling of the configuration space, we biased the sampling of reactions only within a certain inclusion radius. We validated our algorithm using a set of ideally behaving systems undergoing dimerization and polycondensation reactions, for which analytical results are available. For dimerization reactions with various equilibrium constants and initial compositions, the degree of conversion measured in our simulations perfectly matched the reference values given by the analytical equations. We also showed that this agreement is not affected by the arbitrary choice of the inclusion radius or the stiffness of the harmonic bond potential. Next, we showed that our simulations can correctly match the analytical results for the distribution of the degree of polymerization and end-to-end distance of ideal chains in polycondensation reactions. Altogether, we demonstrated that our eb-RxMC simulations correctly sample both reaction and configuration spaces of these reference systems, opening the door to future simulations of more complex interacting macromolecular systems.

Název v anglickém jazyce

The explicit bonding reaction ensemble Monte Carlo method

Popis výsledku anglicky

We present the explicit bonding Reaction ensemble Monte Carlo (eb-RxMC) method, designed to sample reversible bonding reactions in macromolecular systems in thermodynamic equilibrium. Our eb-RxMC method is based on the reaction ensemble method; however, its implementation differs from the latter by the representation of the reaction. In the eb-RxMC implementation, we are adding or deleting bonds between existing particles, instead of inserting or deleting particles with different chemical identities. This new implementation makes the eb-RxMC method suitable for simulating the formation of reversible linkages between macromolecules, which would not be feasible with the original implementation. To enable coupling of our eb-RxMC algorithm with molecular dynamics algorithm for the sampling of the configuration space, we biased the sampling of reactions only within a certain inclusion radius. We validated our algorithm using a set of ideally behaving systems undergoing dimerization and polycondensation reactions, for which analytical results are available. For dimerization reactions with various equilibrium constants and initial compositions, the degree of conversion measured in our simulations perfectly matched the reference values given by the analytical equations. We also showed that this agreement is not affected by the arbitrary choice of the inclusion radius or the stiffness of the harmonic bond potential. Next, we showed that our simulations can correctly match the analytical results for the distribution of the degree of polymerization and end-to-end distance of ideal chains in polycondensation reactions. Altogether, we demonstrated that our eb-RxMC simulations correctly sample both reaction and configuration spaces of these reference systems, opening the door to future simulations of more complex interacting macromolecular systems.

Druh

J_imp - Článek v periodiku v databázi Web of Science

CEP obor

—

OECD FORD obor

10403 - Physical chemistry

Projekt

<a href="/cs/project/GC21-31978J" target="_blank" >GC21-31978J: Simulace reakčních rovnováh v polymerních systémech - vývoj simulačních metod a jejich aplikace</a><br>

Návaznosti

P - Projekt vyzkumu a vyvoje financovany z verejnych zdroju (s odkazem do CEP)

Rok uplatnění

2024

Kód důvěrnosti údajů

S - Úplné a pravdivé údaje o projektu nepodléhají ochraně podle zvláštních právních předpisů

Název periodika

Journal of Chemical Physics

ISSN

0021-9606

e-ISSN

1089-7690

Svazek periodika

161

Číslo periodika v rámci svazku

9

Stát vydavatele periodika

US - Spojené státy americké

Počet stran výsledku

12

Strana od-do

094906

Kód UT WoS článku

001304263600011

EID výsledku v databázi Scopus

2-s2.0-85203104268