40 research outputs found

    Folding kinetics of a polymer [corrigendum]

    Full text link
    In our original article (Phys. Chem. Chem. Phys., 2012, 14, 60446053) a convergence problem resulted in an averaging error in computing the entropy from a set of Wang-Landau Monte-Carlo simulations. Here we report corrected results for the freezing temperature of the homopolymer chain as a function of the range of the non-bonded interaction. We find that the previously reported forward-flux sampling (FFS) and brute-force (BF) simulation results are in agreement with the revised Wang-Landau (WL) calculations. This confirms the utility of FFS for computing crystallisation rates in systems of this kind.Comment: 2 pages, 4 figure

    Efficient Langevin and Monte Carlo sampling algorithms: the case of field-theoretic simulations

    Get PDF
    We introduce Langevin sampling algorithms to field-theoretic simulations (FTSs) of polymers that, for the same accuracy, are ∼10× more efficient than a previously used Brownian dynamics algorithm that used predictor corrector for such simulations, over 10× more efficient than the smart Monte Carlo (SMC) algorithm, and typically over 1000× more efficient than a simple Monte Carlo (MC) algorithm. These algorithms are known as the Leimkuhler–Matthews (the BAOAB-limited) method and the BAOAB method. Furthermore, the FTS allows for an improved MC algorithm based on the Ornstein–Uhlenbeck process (OU MC), which is 2× more efficient than SMC. The system-size dependence of the efficiency for the sampling algorithms is presented, and it is shown that the aforementioned MC algorithms do not scale well with system sizes. Hence, for larger sizes, the efficiency difference between the Langevin and MC algorithms is even greater, although, for SMC and OU MC, the scaling is less unfavorable than for the simple MC

    Determination of tube theory parameters using a simple grid model as an example

    Get PDF
    Although the tube theory is successful in describing entangled polymers qualitatively, a more quantitative description requires precise and consistent definitions of its parameters. Here we investigate the simplest model of entangled polymers, namely a single Rouse chain in a cubic lattice of line obstacles, and illustrate the typical problems and uncertainties of the tube theory. In particular we show that in general one needs 3 entanglement related parameters, but only 2 combinations of them are relevant for the long-time dynamics. Conversely, the plateau modulus can not be determined from these two parameters and requires a more detailed model of entanglements with explicit entanglement forces, such as the slipsprings model. It is shown that for the grid model the Rouse time within the tube is larger than the Rouse time of the free chain, in contrast to what the standard tube theory assumes

    Instability of the Microemulsion Channel in Block Copolymer-Homopolymer Blends

    Get PDF
    Field theoretic simulations are used to predict the equilibrium phase diagram of symmetric blends of AB diblock copolymer with A- and B-type homopolymers. Experiments generally observe a channel of bicontinuous microemulsion (BμE) separating the ordered lamellar (LAM) phase from coexisting homopolymer-rich (A+B) phases. However, our simulations find that the channel is unstable with respect to macrophase separation, in particular, A+B+BμE coexistence at high T and A+B+LAM coexistence at low T. The preference for three-phase coexistence is attributed to a weak attractive interaction between diblock monolayers

    Solid–liquid interfacial free energy of ice Ih, ice Ic, and ice 0 within a mono-atomic model of water via the capillary wave method

    Get PDF
    We apply the capillarywave method, based on measurements of fluctuations in a ribbon-like interfacial geometry, to determine the solid–liquid interfacial free energy for both polytypes of ice I and the recently proposed ice 0 within a mono-atomic model of water. We discuss various choices for the molecular order parameter, which distinguishes solid from liquid, and demonstrate the influence of this choice on the interfacial stiffness. We quantify the influence of discretisation error when sampling the interfacial profile and the limits on accuracy imposed by the assumption of quasi onedimensional geometry. The interfacial free energies of the two ice I polytypes are indistinguishable to within achievable statistical error and the small ambiguity which arises from the choice of order parameter. In the case of ice 0, we find that the large surface unit cell for low index interfaces constrains the width of the interfacial ribbon such that the accuracy of results is reduced. Nevertheless, we establish that the interfacial free energy of ice 0 at its melting temperature is similar to that of ice I under the same conditions. The rationality of a core–shell model for the nucleation of ice I within ice 0 is questioned within the context of our results. © 2017 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/). [http://dx.doi.org/10.1063/1.4975776

    OPTIma: simplifying calorimetry for proton computed tomography in high proton flux environments

    Get PDF
    OBJECTIVE: Proton computed tomography (pCT) offers a potential route to reducing range uncertainties for proton therapy treatment planning, however the current trend towards high current spot scanning treatment systems leads to high proton fluxes which are challenging for existing systems. Here we demonstrate a novel approach to energy reconstruction, referred to as "de-averaging", which allows individual proton energies to be recovered using only a measurement of their integrated energy without the need for spatial information from the calorimeter. APPROACH: The method is evaluated in the context of the OPTIma (Optimising Proton Therapy through Imaging) system which uses a simple, relatively inexpensive, scintillator based calorimeter that reports only the integrated energy deposited by all protons within a cyclotron period, alongside a silicon strip based tracking system capable of reconstructing individual protons in a high flux environment. GEANT4 simulations have been performed to examine the performance of such a system at a modern commercial cyclotron facility using a σ ≈ 10mm beam for currents in the range 10-50 pA at the nozzle. MAIN RESULTS: Apart from in low-density lung tissue, a discrepancy of less than 1% on the Relative Stopping Power is found for all other considered tissues when embedded within a 150 mm spherical Perspex phantom in the 10-30 pA current range, and for some tissues even up to 50 pA. SIGNIFICANCE: By removing the need for the calorimeter system to provide spatial information, it is hoped that the de-averaging approach can facilitate clinically relevant, cost effective and less complex calorimeter systems for performing high current pCTs

    Nucleation of NaCl from aqueous solution: critical sizes, ion-attachment kinetics, and rates

    Get PDF
    Nucleation and crystal growth are important in material synthesis, climate modeling, biomineralization, and pharmaceutical formulation. Despite tremendous efforts, the mechanisms and kinetics of nucleation remain elusive to both theory and experiment. Here we investigate sodium chloride (NaCl) nucleation from supersaturated brines using seeded atomistic simulations, polymorph-specific order parameters, and elements of classical nucleation theory. We find that NaCl nucleates via the common rock salt structure. Ion desolvation - not diffusion - is identified as the limiting resistance to attachment. Two different analyses give approximately consistent attachment kinetics: diffusion along the nucleus size coordinate and reaction-diffusion analysis of approach-to-coexistence simulation data from Aragones et al. [J. Chem. Phys., 2012, 136, 244508]. Our simulations were performed at realistic supersaturations to enable the first direct comparison to experimental nucleation rates for this system. The computed and measured rates converge to a common upper limit at extremely high supersaturation. However, our rate predictions are between 15 and 30 orders of magnitude too fast. We comment on possible origins of the large discrepancies

    LUCAS: A highly accurate yet simple risk calculator that predicts survival of COVID-19 patients using rapid routine tests

    Get PDF
    Background There is an urgent need to develop a simplified risk tool that enables rapid triaging of SARS CoV-2 positive patients during hospital admission, which complements current practice. Many predictive tools developed to date are complex, rely on multiple blood results and past medical history, do not include chest X ray results and rely on Artificial Intelligence rather than simplified algorithms. Our aim was to develop a simplified risk-tool based on five parameters and CXR image data that predicts the 60-day survival of adult SARS CoV-2 positive patients at hospital admission. Methods We analysed the NCCID database of patient blood variables and CXR images from 19 hospitals across the UK contributed clinical data on SARS CoV-2 positive patients using multivariable logistic regression. The initial dataset was non-randomly split between development and internal validation dataset with 1434 and 310 SARS CoV-2 positive patients, respectively. External validation of final model conducted on 741 Accident and Emergency admissions with suspected SARS CoV-2 infection from a separate NHS Trust which was not part of the initial NCCID data set. Findings The LUCAS mortality score included five strongest predictors (lymphocyte count, urea, CRP, age, sex), which are available at any point of care with rapid turnaround of results. Our simple multivariable logistic model showed high discrimination for fatal outcome with the AUC-ROC in development cohort 0.765 (95% confidence interval (CI): 0.738 - 0.790), in internal validation cohort 0.744 (CI: 0.673 - 0.808), and in external validation cohort 0.752 (CI: 0.713 - 0.787). The discriminatory power of LUCAS mortality score was increased slightly when including the CXR image data (for normal versus abnormal): internal validation AUC-ROC 0.770 (CI: 0.695 - 0.836) and external validation AUC-ROC 0.791 (CI: 0.746 - 0.833). The discriminatory power of LUCAS and LUCAS + CXR performed in the upper quartile of pre-existing risk stratification scores with the added advantage of using only 5 predictors. Interpretation This simplified prognostic tool derived from objective parameters can be used to obtain valid predictions of mortality in patients within 60 days SARS CoV-2 RT-PCR results. This free-to-use simplified tool can be used to assist the triage of patients into low, moderate, high or very high risk of fatality and is available at https://mdscore.net/. What is already known on this topic? Clinical prediction models such as NEWS2 is currently used in practice as mortality risk assessment. In a rapid response to support COVID-19 patient assessment and resource management, published risk tools and models have been found to have a high risk of bias and therefore cannot be translated into clinical practice. What this study adds? A newly developed and validated risk tool (LUCAS) based on rapid and routine blood tests predicts the mortality of patients infected with SARS-CoV-2 virus. This prediction model has both high and robust predictive power and has been tested on an external set of patients and therefore can be used to effectively triage patients when resources are limited. In addition, LUCAS can be used with chest imaging information and NEWS2 score

    Nucleation barrier reconstruction via the seeding method in a lattice model with competing nucleation pathways

    Get PDF
    We study a three-species analogue of the Potts lattice gas model of nucleation from solution in a regime where partially disordered solute is a viable thermodynamic phase. Using a multicanonical sampling protocol, we compute phase diagrams for the system, from which we determine a parameter regime where the partially disordered phase is metastable almost everywhere in the temperature–fugacity plane. The resulting model shows non-trivial nucleation and growth behaviour, which we examine via multidimensional free energy calculations. We consider the applicability of the model in capturing the multi-stage nucleation mechanisms of polymorphic biominerals (e.g., CaCO3). We then quantitatively explore the kinetics of nucleation in our model using the increasingly popular “seeding” method. We compare the resulting free energy barrier heights to those obtained via explicit free energy calculations over a wide range of temperatures and fugacities, carefully considering the propagation of statistical error. We find that the ability of the “seeding” method to reproduce accurate free energy barriers is dependent on the degree of supersaturation, and severely limited by the use of a nucleation driving force ∆µ computed for bulk phases. We discuss possible reasons for this in terms of underlying kinetic assumptions, and those of classical nucleation theory. C 2016 Author(s

    Non-Gaussian nature of glassy dynamics by cage to cage motion

    Get PDF
    A model based on a single Brownian particle moving in a periodic effective field is used to understand the non-Gaussian dynamics in glassy systems of cage escape and subsequent recaging, often thought to be caused by a heterogeneous glass structure. The results are compared to molecular-dynamics simulations of systems with varying complexity: quasi-two-dimensional colloidlike particles, atactic polystyrene, and a dendritic glass. The model nicely describes generic features of all three topologically different systems, in particular around the maximum of the non-Gaussian parameter. This maximum is a measure for the average distance between cages
    corecore