Related papers: Simulation of weak polyelectrolytes: A comparison …
We introduce a Monte-Carlo method that allows for the simulation of a polymeric phase containing a weak polyelectrolyte, which is coupled to a reservoir at a fixed pH, salt concentration and total concentration of a weak polyprotic acid.…
We present a simulation method that allows us to calculate the titration curves for systems undergoing protonation/deprotonation reactions -- such as charged colloidal suspensions with acidic/basic surface groups, polyelectrolytes,…
Variational methods are used to calculate structural and thermodynamical properties of a titrating polyelectrolyte in a discrete representation. The Coulomb interactions are emulated by harmonic repulsive forces, the force constants being…
The structural dynamics of biological macromolecules, such as proteins, DNA/RNA, or complexes thereof, are strongly influenced by protonation changes of their typically many titratable groups, which explains their sensitivity to pH changes.…
In complex systems with many degrees of freedom such as peptides and proteins there exist a huge number of local-minimum-energy states. Conventional simulations in the canonical ensemble are of little use, because they tend to get trapped…
In a recent work we proposed a coarse-grained methodology for studying the response of peptides when simulated at different values of pH; in this work we extend the methodology to analyze the pH-dependent behavior of coiled coils. This…
We study the effective interaction between differently charged polyelectrolyte-colloid complexes in electrolyte solutions via Monte Carlo simulations. These complexes are formed when short and flexible polyelectrolyte chains adsorb onto…
Redox processes are important in chemistry, with applications in biomedicine, chemical analysis, among others. As many redox experiments are also performed at a fixed value of pH, having an efficient computational method to support…
New hybrid Molecular Dynamics-Monte Carlo methods are proposed to increase the efficiency of constant-pressure simulations. Two variations of the isobaric Molecular Dynamics component of the algorithms are considered. In the first, we use…
We propose a framework to study the spatio-temporal evolution of liquid-liquid phase separation of weak polyelectrolytes in ionic solutions. Unlike strong polyelectrolytes, which carry a fixed charge, the charge state of weak…
Constant potential method molecular dynamics simulation (CPM MD) enables the accurate modelling of atomistic electrode charges when studying the electrode-electrolyte interface at the nanoscale. Here we extend the theoretical framework of…
The strongly coupled electron liquid provides a unique opportunity to study the complex interplay of strong coupling with quantum degeneracy effects and thermal excitations. To this end, we carry out extensive \textit{ab initio} path…
We use a reactive Monte Carlo simulation method and primitive model of electrolyte to study acid-base equilibrium that controls charge regulation in colloidal systems. The simulations are performed in a semi-grand canonical ensemble in…
In biomolecular systems (especially all-atom models) with many degrees of freedom such as proteins and nucleic acids, there exist an astronomically large number of local-minimum-energy states. Conventional simulations in the canonical…
The influence of the chain degree of ionization on the adsorption of weak polyelectrolytes on neutral and on oppositely and likely charged surfaces is investigated for the first time, by means of Monte Carlo simulations with the mesoscopic…
In complex systems with many degrees of freedom such as spin glass and biomolecular systems, conventional simulations in canonical ensemble suffer from the quasi-ergodicity problem. A simulation in generalized ensemble performs a random…
A better understanding of interfacial mechanisms is needed to improve the performances of electrochemical devices. Yet, simulating an electrode surface at fixed electrolyte composition remains a challenge. Here we apply a finite electric…
Conventional simulations of complex systems in the canonical ensemble suffer from the quasi-ergodicity problem. A simulation in generalized ensemble overcomes this difficulty by performing a random walk in potential energy space and other…
In a recent review~Landsgesell et al., Soft Matter {\bf 15}, 1155 (2019) stated that $\text{pH} - \text{pK}_a$ is a ``universal parameter" for titrating systems. We show that this is not the case. This broken symmetry has important…
Finding and sampling multiple reaction channels for molecular transitions remains an important challenge in physical chemistry. Here we show that the weighted ensemble (WE) path sampling method can readily sample multiple channels. In a…