BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Patrick B. Warren

Patrick B. Warren contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.stat-mechcond-mat.softQuantitative MethodsBiological PhysicsBiomoleculesMolecular Networksphysics.comp-phphysics.hist-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2014arXiv

Screening properties of four mesoscale smoothed charge models, with application to dissipative particle dynamics

We extend our previous study [J. Chem. Phys. 138, 204907 (2013)] to quantify the screening properties of four mesoscale smoothed charge models used in dissipative particle dynamics. Using a combination of the hypernetted chain integral equation closure and the random phase approximation, we identify regions where the models exhibit a real-valued screening length, and the extent to which this agrees with the Debye length in the physical system. We find that the second moment of the smoothed charge distribution is a good predictor of this behaviour. We are thus able to recommend a consistent set of parameters for the models.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Malliavin weight sampling: a practical guide

Malliavin weight sampling (MWS) is a stochastic calculus technique for computing the derivatives of averaged system properties with respect to parameters in stochastic simulations, without perturbing the system's dynamics. It applies to systems in or out of equilibrium, in steady state or time-dependent situations, and has applications in the calculation of response coefficients, parameter sensitivities and Jacobian matrices for gradient-based parameter optimisation algorithms. The implementation of MWS has been described in the specific contexts of kinetic Monte Carlo and Brownian dynamics simulation algorithms. Here, we present a general theoretical framework for deriving the appropriate MWS update rule for any stochastic simulation algorithm. We also provide pedagogical information on its practical implementation.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

No-go theorem in many body dissipative particle dynamics

Many body dissipative particle dynamics (MDPD) is a particle-based simulation method in which the interaction potential is a sum of self energies depending on locally-sampled density variables. This functional form gives rise to density-dependent pairwise forces, however not all such force laws are derivable from a potential and the integrability condition for this to be the case provides a strong constraint. A strategy to assess the implications of this constraint is illustrated here by the derivation of a useful no-go theorem for multicomponent MDPD.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Phase behaviour and the random phase approximation for ultrasoft restricted primitive models

Phase separation of the ultrasoft restricted primitive model (URPM) with Gaussian charges is re-investigated in the random phase approximation (RPA)---the 'Level A' approximation discussed by Nikoubashman, Hansen and Kahl [J. Chem. Phys. 137, 094905 (2012)]. We find that the RPA predicts a region of low temperature vapour-liquid coexistence, with a critical density much lower than that observed in either simulations or more refined approximations (we also remark that the RPA critical point for a related model with Bessel charges can be solved analytically). This observation suggests that the hierarchy of approximations introduced by Nikoubashman et al. should be analogous to those introduced by Fisher and Levin for the restricted primitive model [Phys. Rev. Lett. 71, 3826 (1993)], which makes the inability of these approximations to capture the observed URPM phase behaviour even more worthy of investigation.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Revisiting Riccioli's free fall calculations

In 1651 Giovanni Riccioli reported the earliest accurate measurements the acceleration due to gravity, g, from pendulum-timed free fall experiments. The use of Huygens' pendulum formula (published 1673) allows one to deduce the pendulum length from this data, free from assumptions about the conversion to modern units, and independent of the actual value of g. When this length is compared to the reported pendulum length, a 15% systematic error is revealed. This could perhaps be attributed to the difficulty Riccioli faced in subdividing the contemporary unit of length (the Roman foot) to the requisite millimeter accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Screening properties of Gaussian electrolyte models, with application to dissipative particle dynamics

We investigate the screening properties of Gaussian charge models of electrolyte solutions by analysing the asymptotic behaviour of the pair distribution functions. We use a combination of Monte-Carlo simulations with the hyper-netted chain integral equation closure, and the random phase approximation, to establish the conditions under which a screening length is well defined and the extent to which it matches the expected Debye length. For practical applications, for example in dissipative particle dynamics, we are able to summarise our results in succinct rules-of-thumb which can be used for mesoscale modeling of electrolyte solutions. We thereby establish a solid foundation for future work, such as the systematic incorporation of specific ion effects.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Computing sensitivity coefficients in Brownian dynamics simulations by Malliavin weight sampling

We present a method for computing parameter sensitivities and response coefficients in Brownian dynamics simulations. The method involves tracking auxiliary variables (Malliavin weights) in addition to the usual particle positions, in an unperturbed simulation. The Malliavin weights sample the derivatives of the probability density with respect to the parameters of interest and are also interesting dynamical objects in themselves. Malliavin weight sampling is simple to implement, applies to equilibrium or nonequilibrium, steady state or time-dependent systems, and scales more efficiently than standard finite difference methods

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Self-consistent field theory for the interactions between keratin intermediate filaments

We have developed a model for the interactions between keratin intermediate filaments based on self-consistent field theory. The intermediate filaments are represented by charged surfaces, and the disordered terminal domains of the keratins are represented by charged heteropolymers grafted to these surfaces. We estimate the system is close to a charge compensation point where the heteropolymer grafting density is matched to the surface charge density. Using a protein model with amino acid resolution for the terminal domains, we find that the terminal chains can mediate a weak attraction between the keratin surfaces. The origin of the attraction is a combination of bridging and electrostatics. The attraction disappears when the system moves away from the charge compensation point, or when excess small ions and/or NMF-representing free amino acids are added. These results are in concordance with experimental observations, and support the idea that the interaction between keratin filaments, and ultimately in part the elastic properties of the keratin-containing tissue, is controlled by a combination of the physico-chemical properties of the disordered terminal domains and the composition of the medium in the inter-filament region.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Steady-state parameter sensitivity in stochastic modeling via trajectory reweighting

Parameter sensitivity analysis is a powerful tool in the building and analysis of biochemical network models. For stochastic simulations, parameter sensitivity analysis can be computationally expensive, requiring multiple simulations for perturbed values of the parameters. Here, we use trajectory reweighting to derive a method for computing sensitivity coefficients in stochastic simulations without explicitly perturbing the parameter values, avoiding the need for repeated simulations. The method allows the simultaneous computation of multiple sensitivity coefficients. Our approach recovers results originally obtained by application of the Girsanov measure transform in the general theory of stochastic processes [A. Plyasunov and A. P. Arkin, J. Comp. Phys. 221, 724 (2007)]. We build on these results to show how the method can be used to compute steady-state sensitivity coefficients from a single simulation run, and we present various efficiency improvements. For models of biochemical signaling networks the method has a particularly simple implementation. We demonstrate its application to a signaling network showing stochastic focussing and to a bistable genetic switch, and present exact results for models with linear propensity functions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The shape of a ponytail and the statistical physics of hair fiber bundles

A general continuum theory for the distribution of hairs in a bundle is developed, treating individual fibers as elastic filaments with random intrinsic curvatures. Applying this formalism to the iconic problem of the ponytail, the combined effects of bending elasticity, gravity, and orientational disorder are recast as a differential equation for the envelope of the bundle, in which the compressibility enters through an 'equation of state'. From this, we identify the balance of forces in various regions of the ponytail, extract a remarkably simple equation of state from laboratory measurements of human ponytails, and relate the pressure to the measured random curvatures of individual hairs.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Flux networks in metabolic graphs

A metabolic model can be represented as bipartite graph comprising linked reaction and metabolite nodes. Here it is shown how a network of conserved fluxes can be assigned to the edges of such a graph by combining the reaction fluxes with a conserved metabolite property such as molecular weight. A similar flux network can be constructed by combining the primal and dual solutions to the linear programming problem that typically arises in constraint-based modelling. Such constructions may help with the visualisation of flux distributions in complex metabolic networks. The analysis also explains the strong correlation observed between metabolite shadow prices (the dual linear programming variables) and conserved metabolite properties. The methods were applied to recent metabolic models for Escherichia coli, Saccharomyces cerevisiae, and Methanosarcina barkeri. Detailed results are reported for E. coli; similar results were found for the other organisms.

0 citations0 reviewsNo code linkedOpen access