[Paper Review] How does the scaling for the polymer chain in the dissipative particle dynamics hold?

Jaroslav Ilnytskyi, Yurij Holovatch|ArXiv.org|Dec 3, 2007

Scientific Research and Discoveries6 references3 citations

TL;DR

This study investigates scaling laws for polymer chains in dissipative particle dynamics (DPD) simulations, analyzing end-to-end distance, radius of gyration, and hydrodynamic radius across varying solvent qualities. Using linear fits, correction-to-scaling forms, and histogram fitting, the authors find the scaling exponent ν ≈ 0.59 ± 0.01, closely matching the theoretical value ν = 0.5882, confirming that DPD accurately captures universal scaling behavior for flexible polymers in good solvents.

ABSTRACT

We performed a series of simulations for a linear polymer chain in a solvent using dissipative particle dynamics to check the scaling relations for the end-to-end distance, radius of gyration and hydrodynamic radius in three dimensions. The polymer chains of up to 80 beads in explicit solvent of various quality are studied. To extract the scaling exponent ν, the data are analyzed using linear fits, correction-to-scaling forms and analytical fits to the histograms of radius of gyration distribution. For certain combinations of the polymer characteristics and solvent quality, the correction-to-scaling terms are found to be essential while for the others these are negligibly small. In each particular case the final value for the exponent νwas chosen according to the best least-squares fit. The values of νobtained in this way are found within the interval ν=0.55-0.61 but are concentrated mostly around 0.59, which is very close to the best known theoretical result ν=0.588. The existence of this interval is attributed both to the peculiarities of the method and to the moderate chain lengths being simulated. Within this shortcoming, the polymer chain in this kind of modeling is found to satisfy the scaling relations for all three radii being considered

Motivation & Objective

• To validate whether dissipative particle dynamics (DPD) accurately reproduces universal scaling laws for flexible polymer chains in good solvents.
• To assess the influence of solvent quality and chain length on scaling exponent ν for end-to-end distance, radius of gyration, and hydrodynamic radius.
• To determine whether correction-to-scaling terms are essential for accurate estimation of ν in DPD simulations.
• To evaluate the consistency and reliability of ν values across different analysis techniques and simulation parameters.
• To establish a benchmark for using DPD in simulating more complex polymer architectures such as branched or star-like polymers.

Proposed method

• Conducted DPD simulations of linear polymer chains with up to 80 beads in explicit solvent of varying quality (athermal, good, very good solvent).
• Computed three metric properties: end-to-end distance $ R_{1N} $, radius of gyration $ R_g $, and hydrodynamic radius $ R_h $, for each chain configuration.
• Applied linear fitting to $ ig angle R_{1N}^2 \big\rangle \sim N^{2\nu} $ to estimate ν across different chain length ranges.
• Used correction-to-scaling fits incorporating $ B/N^\Delta $ terms with $ \Delta = 0.478 $, based on field-theoretic renormalization group estimates.
• Fitted the histogram of radius of gyration distributions to extract ν, enhancing robustness against finite-size effects.
• Selected the best-fit ν value using least-squares minimization across multiple fitting strategies and data subsets.

Experimental results

Research questions

• RQ1 Does the dissipative particle dynamics (DPD) method reproduce the universal scaling exponent ν ≈ 0.588 for flexible polymer chains in good solvents?
• RQ2 How significant are correction-to-scaling terms in DPD simulations, and when are they necessary for accurate ν estimation?
• RQ3 How does the scaling exponent ν vary with solvent quality, and is the variation within the expected range of statistical and methodological uncertainty?
• RQ4 To what extent do different fitting methods (linear, correction-to-scaling, histogram fitting) affect the estimated value of ν?
• RQ5 Can DPD reliably model scaling behavior for polymer characteristics such as end-to-end distance, radius of gyration, and hydrodynamic radius?

Key findings

• The scaling exponent ν was found to lie in the interval 0.55–0.61, with most values concentrated around ν ≈ 0.59, closely matching the theoretical value ν = 0.5882 ± 0.0011.
• Correction-to-scaling terms were essential for accurate ν estimation in most solvent and polymer characteristic combinations, indicating finite-size effects significantly influence results.
• No significant drift in ν was observed when solvent quality was varied within the good solvent regime, though minor variations depended on the analysis method.
• The hydrodynamic radius $ R_h $ exhibited scaling behavior consistent with the other radii, supporting the validity of DPD for hydrodynamic property modeling.
• The best-fit ν value was consistently close to 0.59 across all three radii and solvent conditions, with the histogram of ν values peaking near the theoretical prediction.
• The study confirms that DPD simulations of polymer chains satisfy universal scaling laws for static properties, validating its use for complex polymer architectures.