BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Achille Giacometti

Achille Giacometti appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.softcond-mat.stat-mechBiomoleculesphysics.chem-phBiological Physicscond-matMolecular Networks

Catalog footprint

What is connected

38works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

38 published item(s)

preprint2026arXiv

Information Transmission and Processing in G-Protein-Coupled-Receptor Complexes

G-protein-coupled receptors (GPCRs) are central to cellular information processing, yet the physical principles governing their switching behavior remain incompletely understood. We present a first principles theoretical framework, grounded in nonequilibrium thermodynamics, to describe GPCR switching as observed in light-controlled impedance assays. The model identifies two fundamental control parameters: (1) ATP/GTP-driven chemical flux through the receptor complex, and (2) the free-energy difference between phosphorylated and dephosphorylated switch states. Together, these parameters defin the switch configuration. The model predicts that GPCRs can occupy one of three quasi-stable configurations, each corresponding to a local maximum in information transmission. Active states support chemical flux and exist in an on or off switch configuration, whereas inactive states lack flux, introducing a distinction absent in conventional phosphorylation models. The model takes two ligand-derived inputs: fixed structural features and inducible conformations (e.g. cis or trans). It shows that phosphatase activity, modeled as an energy barrier, chiefly governs on/off occupancy, whereas the kinase sustains flux without directly determining the switch configuration. Comparison with experimental data confirms the predicted existence of multiple quasi-stable states modulated by ligand conformation. Importantly, this framework generalizes beyond GPCRs to encompass a wider class of biological switching systems driven by nonequilibrium chemical flux.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Local sequence-structure relationships in proteins

We seek to understand the interplay between amino acid sequence and local structure in proteins. Are some amino acids unique in their ability to fit harmoniously into certain local structures? What is the role of sequence in sculpting the putative native state folds from myriad possible conformations? In order to address these questions, we represent the local structure of each C-alpha atom of a protein by just two angles, theta and mu, and we analyze a set of more than 4000 protein structures from the PDB. We use a hierarchical clustering scheme to divide the 20 amino acids into six distinct groups based on their similarity to each other in fitting local structural space. We present the results of a detailed analysis of patterns of amino acid specificity in adopting local structural conformations and show that the sequence-structure correlation is not very strong compared to a random assignment of sequence to structure. Yet, our analysis may be useful to determine an effective scoring rubric for quantifying the match of an amino acid to its putative local structure.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Proteins -- a celebration of consilience

Proteins are the common constituents of all living cells. They are molecular machines that interact with each other as well as with other cell products and carry out a dizzying array of functions with distinction. These interactions follow from their native state structures and therefore understanding sequence-structure relationships is of fundamental importance. What is quite remarkable about proteins is that their understanding necessarily straddles several disciplines. The importance of geometry in defining protein native state structure, the constraints placed on protein behavior by mathematics and physics, the need for proteins to obey the laws of quantum chemistry, and the rich role of evolution and biology all come together in defining protein science. Here we review ideas from the literature and present an interdisciplinary framework that aims to marry ideas from Plato and Darwin and demonstrates an astonishing consilience between disciplines in describing proteins. We discuss the consequences of this framework on protein behavior.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

An integrated DEM-FEM approach to study breakage in packing of glass cartridges on a conveyor belt

The use of glass for pharmaceutical new applications such as high-technology drugs, requires the strictest container inertness. A common theme of paramount importance in glass container integrity preservation is the detailed mechanism driving the sudden failure due the crack propagation. Using a combination of Discrete Element Method (DEM) and Finite Element Method (FEM), a stress map for glass cartridges packed into an accumulation table and transported by a conveyor belt at a fixed velocity is obtained under realistic conditions. The DEM calculation provides a full description of the dynamics of the cartridges, as approximated by an equivalent sphere, as well as the statistics of the multiple collisions. The FEM calculation exploits this input to provide the maximum principal stress of different pairs as a function of time. Our analysis shows that, during their transportation on the conveyor belt, the cartridges are subject to several shocks of varying intensities. Under these conditions, a crack may originate inside the cartridge in the area of maximal tensile stress, and propagate outward. Estimated stresses are found in good agreement with real systems.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Effective stiffness and formation of secondary structures in a protein-like model

We use Wang-Landau and replica exchange techniques to study the effect of an increasing stiffness on the formation of secondary structures in protein-like systems. Two possible models are considered. In both models, a polymer chain is formed by tethered beads where non-consecutive backbone beads attract each other via a square-well potential representing the tendency of the chain to fold. In addition, smaller hard spheres are attached to each non-terminal backbone bead along the direction normal to the chain to mimic the steric hindrance of side chains in real proteins. The two models, however, differ in the way bending rigidity is enforced. In the first model, partial overlap between consecutive beads is allowed. This reduces the possible bending angle between consecutive bonds thus producing an effective entropic stiffness that competes with a short-range attraction, and leads to a formation of secondary structures characteristic of proteins. We discuss the low-temperature phase diagram as a function of increasing interpenetration, and find a transition from a planar, beta-like structure, to helical shape. In the second model, an energetic stiffness is explicitly introduced by imposing an infinitely large energy penalty for bending above a critical angle between consecutive bonds, and no penalty below it. The low-temperature phase of this model does not show any sign of protein-like secondary structures. At intermediate temperatures, however, where the chain is still in the coil conformation but stiffness is significant, we find the two models to predict a quite similar dependence of the persistence length as a function of the stiffness. This behaviour is rationalized in terms of a simple geometrical mapping between the two models. Finally, we discuss the effect of shrinking side chains to zero, and find the above mapping to still hold true.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

From polymers to proteins: effect of side chains and broken symmetry in the formation of secondary structures within a Wang-Landau approach

We study the equilibrium properties of a single flexible homopolymers where consecutive monomers are represented by impenetrable hard spherical beads tangential to each other, and non-consecutive monomers interact via a square-well potential. To this aim, we use both replica exchange canonical simulations and micro-canonical Wang-Landau techniques. We perform a close comparative analysis of the corresponding results, and find perfect agreement between the two methods as well as the past results. The model is then refined in two different directions. By allowing partial overlapping between consecutive beads, we break the spherical symmetry and thus provide a severe constraint on the possible conformations of the chain. This leads to a single helix or a double helix ground state, depending on the range of the interactions. Alternatively, we introduce additional spherical beads at specific positions to represent the steric hindrance of the side chains in real proteins. The ground state in this case is found to be again a globule, thus showing that contrary to some recent suggestions, side chains alone are not sufficient to insure a protein-like structure. Finally, we show that a combination of the above two effects increases the stability of the obtaining secondary structures, in agreement with previous results. The fundamental role played by the range of the square-well attraction is highlighted, and it is shown that they play a role similar to that found in simple liquids and polymers. Perspectives in terms of protein folding are finally discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Phase diagram of heteronuclear Janus dumbbells

Using Aggregation-Volume-Bias Monte Carlo simulations along with Successive Umbrella Sampling and Histogram Re-weighting, we study the phase diagram of a system of dumbbells formed by two touching spheres having variable sizes, as well as different interaction properties. The first sphere ($h$) interacts with all other spheres belonging to different dumbbells with a hard-sphere potential. The second sphere ($s$) interacts via a square-well interaction with other $s$ spheres belonging to different dumbbells and with a hard-sphere potential with all remaining $h$ spheres. We focus on the region where the $s$ sphere is larger than the $h$ sphere, as measured by a parameter $1\le α\le 2 $ controlling the relative size of the two spheres. As $α\to 2$ a simple fluid of square-well spheres is recovered, whereas $α\to 1$ corresponds to the Janus dumbbell limit, where the $h$ and $s$ spheres have equal sizes. Many phase diagrams falling into three classes are observed, depending on the value of $α$. The $1.8 \le α\le 2$ is dominated by a gas-liquid phase separation very similar to that of a pure square-well fluid with varied critical temperature and density. When $1.3 \le α\le 1.8$ we find a progressive destabilization of the gas-liquid phase diagram by the onset of self-assembled structures, that eventually lead to a metastability of the gas-liquid transition below $α=1.2$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Bridging and depletion mechanisms in colloid-colloid effective interactions: A reentrant phase diagram

A general class of nonadditive sticky-hard-sphere binary mixtures, where small and large spheres represent the solvent and the solute, respectively, is introduced. The solute-solute and solvent-solvent interactions are of hard-sphere type, while the solute-solvent interactions are of sticky-hard-sphere type with tunable degrees of size nonadditivity and stickiness. Two particular and complementary limits are studied using analytical and semi-analytical tools. The first case is characterized by zero nonadditivity, lending itself to a Percus-Yevick approximate solution from which the impact of stickiness on the spinodal curves and on the effective solute-solute potential is analyzed. In the opposite nonadditive case, the solvent-solvent diameter is zero and the model can then be reckoned as an extension of the well-known Asakura-Oosawa model with additional sticky solute-solvent interaction. This latter model has the property that its exact effective one-component problem involves only solute-solute pair potentials for size ratios such that a solvent particle fits inside the interstitial region of three touching solutes. In particular, we explicitly identify the three competing physical mechanisms (depletion, pulling, and bridging) giving rise to the effective interaction. Some remarks on the phase diagram of these two complementary models are also addressed through the use of the Noro-Frenkel criterion and a first-order perturbation analysis. Our findings suggest reentrance of the fluid-fluid instability as solvent density (in the first model) or adhesion (in the second model) is varied. Some perspectives in terms of the interpretation of recent experimental studies of microgels adsorbed onto large polystyrene particles are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Chiral Self-Assembly of Helical Particles

The shape of the building blocks plays a crucial role in directing self-assembly towards desired architectures. Out of the many different shapes, helix has a unique position. Helical structures are ubiquitous in nature and a helical shape is exhibited by the most important biopolymers like polynucleotides, polypeptides and polysaccharides as well as by cellular organelles like flagella. Helical particles can self-assemble into chiral superstructures, which may have a variety of applications, e.g. as photonic (meta)materials. However,a clear and definite understanding of these structures has not been entirely achieved yet. We have recently undertaken an extensive investigation on the phase behaviour of hard helical particles, using numerical simulations and classical density functional theory. Here we present a detailed study of the phase diagram of hard helices as a function of their morphology. This includes a variety of liquid-crystal phases, with different degrees of orientational and positional ordering. We show how, by tuning the helix parameters, it is possible to control the organization of the system. Starting from slender helices, whose phase behaviour is similar to that of rodlike particles, the increase in curliness leads to the onset of azimuthal correlations between the particles and the formation of phases specific to helices. These phases feature a new kind of screw order, of which there is experimental evidence in colloidal suspensions of helical flagella.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Phase diagram of the ground states of DNA condensates

Phase diagram of the ground states of DNA in a bad solvent is studied for a semi-flexible polymer model with a generalized local elastic bending potential characterized by a nonlinearity parameter $x$ and effective self-attraction promoting compaction. $x=1$ corresponds to the worm-like chain model. Surprisingly, the phase diagram as well as the transition lines between the ground states are found to be a function of $x$. The model provides a simple explanation for the results of prior experimental and computational studies and makes predictions for the specific geometries of the ground states. The results underscore the impact of the form of the microscopic bending energy at macroscopic observable scales.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Reference Interaction Site Model and Optimized Perturbation theories of colloidal dumbbells with increasing anisotropy

We investigate thermodynamic properties of anisotropic colloidal dumbbells in the frameworks provided by the Reference Interaction Site Model (RISM) theory and an Optimized Perturbation Theory (OPT), this latter based on a fourth-order high-temperature perturbative expansion of the free energy, recently generalized to molecular fluids. Our model is constituted by two identical tangent hard spheres surrounded by square-well attractions with same widths and progressively different depths. Gas-liquid coexistence curves are obtained by predicting pressures, free energies, and chemical potentials. In comparison with previous simulation results, RISM and OPT agree in reproducing the progressive reduction of the gas-liquid phase separation as the anisotropy of the interaction potential becomes more pronounced; in particular, the RISM theory provides reasonable predictions for all coexistence curves, bar the strong anisotropy regime, whereas OPT performs generally less well. Both theories predict a linear dependence of the critical temperature on the interaction strength, reproducing in this way the mean-field behavior observed in simulations; the critical density~--~that drastically drops as the anisotropy increases~--~turns to be less accurate. Our results appear as a robust benchmark for further theoretical studies, in support to the simulation approach, of self-assembly in model colloidal systems.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A unified description of solvent effects in the helix-coil transition

We analyze the problem of the helix-coil transition in explicit solvents analytically by using spin-based models incorporating two different mechanisms of solvent action: explicit solvent action through the formation of solvent-polymer hydrogen bonds that can compete with the intrinsic intra-polymer hydrogen bonded configurations (competing interactions) and implicit solvent action, where the solvent-polymer interactions tune biopolymer configurations by changing the activity of the solvent (non-competing interactions). The overall spin Hamiltonian is comprised of three terms: the background \emph{in vacuo} Hamiltonian of the "Generalized Model of Polypeptide Chain" type and two additive terms that account for the two above mechanisms of solvent action. We show that on this level the solvent degrees of freedom can be {\sl explicitly} and {\sl exactly} traced over, the ensuing effective partition function combining all the solvent effects in a unified framework. In this way we are able to address helix-coil transitions for polypeptides, proteins, and DNA, with different buffers and different external constraints. Our spin-based effective Hamiltonian is applicable for treatment of such diverse phenomena as cold denaturation, effects of osmotic pressure on the cold and warm denaturation, complicated temperature dependence of the hydrophobic effect as well as providing a conceptual base for understanding the behavior of Intrinsically Disordered Proteins and their analogues.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Cluster formation and phase separation in heteronuclear Janus dumbbells

We have recently investigated the phase behaviour of model colloidal dumbbells constituted by two identical tangent hard spheres, with the first one being surrounded by an attractive square-well interaction (Janus dumbbells, Munaò G et al 2014 Soft Matter 10 5269). Here we extend our previous analysis by introducing in the model the size asymmetry of the hard-core diameters, and study the enriched phase scenario thereby obtained. By employing standard Monte Carlo simulations we show that in such "heteronuclear Janus dumbbells" a larger hard-sphere site promotes the formation of clusters, whereas in the opposite condition a gas-liquid phase separation takes place, with a narrow interval of intermediate asymmetries wherein the two phase behaviours may compete. In addition, some peculiar geometrical arrangements, such as lamellae, are observed only around the perfectly symmetric case. A qualitative agreement is found with recent experimental results, where it is shown that the roughness of molecular surfaces in heterogeneous dimers leads to the formation of colloidal micelles.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

From rods to helices: evidence of a screw-like nematic phase

Evidence of a special chiral nematic phase is provided using numerical simulation and Onsager theory for systems of hard helical particles. This phase appears at the high density end of the nematic phase, when helices are well aligned, and is characterized by the C$_2$ symmetry axes of the helices spiraling around the nematic director with periodicity equal to the particle pitch. This coupling between translational and rotational degrees of freedom allows a more efficient packing and hence an increase of translational entropy. Suitable order parameters and correlation functions are introduced to identify this screw-like phase, whose main features are then studied as a function of radius and pitch of the helical particles. Our study highlights the physical mechanism underlying a similar ordering observed in colloidal helical flagella [E. Barry et al. \textit{Phys. Rev. Lett.} \textbf{96}, 018305 (2006)] and raises the question of whether it could be observed in other helical particle systems, such as DNA, at sufficiently high densities.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

From square-well to Janus: Improved algorithm for integral equation theory and comparison with thermodynamic perturbation theory within the Kern-Frenkel model

Building upon past work on the phase diagram of Janus fluids [Sciortino et al., Phys. Rev. Lett. \textbf{103}, 237801 (2009)], we perform a detailed study of integral equation theory of the Kern-Frenkel potential with coverage that is tuned from the isotropic square-well fluid to the Janus limit. An improved algorithm for the reference hypernetted-chain (RHNC) equation for this problem is implemented that significantly extends the range of applicability of RHNC. Results for both structure and thermodynamics are presented and compared with numerical simulations. Unlike previous attempts, this algorithm is shown to be stable down to the Janus limit, thus paving the way for analyzing the frustration mechanism characteristic of the gas-liquid transition in the Janus system. The results are also compared with Barker-Henderson thermodynamic perturbation theory on the same model. We then discuss the pros and cons of both approaches within a unified treatment. On balance, RHNC integral equation theory, even with an isotropic hard-sphere reference system, is found to be a good compromise between accuracy of the results, computational effort, and uniform quality to tackle self-assembly processes in patchy colloids of complex nature. Further improvement in RHNC however clearly requires an \emph{anisotropic} reference bridge function.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

From toroidal to rod-like condensates of semiflexible polymers

The competition between toroidal and rod-like conformations as possible ground states for DNA condensation is studied as a function of the stiffness, the length of the DNA and the form of the long-range interactions between neighboring molecules, using analytical theory supported by Monte Carlo simulations. Both conformations considered are characterized by a local nematic order with hexagonal packing symmetry of neighboring DNA molecules, but differ in global configuration of the chain and the distribution of its curvature as it wraps around to form a condensate. The long-range interactions driving the DNA condensation are assumed to be of the form pertaining to the attractive depletion potential as well as the attractive counterion induced soft potential. In the stiffness-length plane we find a transition between rod-like to toroid condensate for increasing stiffness at a fixed chain length $L$. Strikingly, the transition line is found to have a $L^{1/3}$ dependence irrespective of the details of the long-range interactions between neighboring molecules. When realistic DNA parameters are used, our description reproduces rather well some of the experimental features observed in DNA condensates.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Left or right cholesterics? A matter of helix handedness and curliness

We have investigated the relationship between the morphology of helical particles and the features of the cholesteric (N$^\ast $) phase that they form. Using an Onsager-like theory, applied to systems of hard helices, we show that the cholesteric handedness and pitch depend on both the pitch and the curliness of the particles. The theory leads to the definition of pseudoscalars that correlate the helical features of the phase to the chirality of the excluded volume of the constituent particles.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Self-assembly of hard helices: a rich and unconventional polymorphism

Hard helices can be regarded as a paradigmatic elementary model for a number of natural and synthetic soft matter systems, all featuring the helix as their basic structural unit: from natural polynucleotides and polypeptides to synthetic helical polymers; from bacterial flagella to colloidal helices. Here we present an extensive investigation of the phase diagram of hard helices using a variety of methods. Isobaric Monte Carlo numerical simulations are used to trace the phase diagram: on going from the low-density isotropic to the high-density compact phases, a rich polymorphism is observed exhibiting a special chiral screw-like nematic phase and a number of chiral and/or polar smectic phases. We present a full characterization of the latter, showing that they have unconventional features, ascribable to the helical shape of the constituent particles. Equal area construction is used to locate the isotropic-to-nematic phase transition, and results are compared with those stemming from an Onsager-like theory. Density functional theory is also used to study the nematic-to-screw-nematic phase transition: within the simplifying assumption of perfectly parallel helices, we compare different levels of approximation, that is second- and third-virial expansions and Parsons-Lee correction.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Janus fluid with fixed patch orientations: theory and simulations

We study thermophysical properties of a Janus fluid with constrained orientations, using analytical techniques and numerical simulations. The Janus character is modeled by means of a Kern-Frenkel potential where each sphere has one hemisphere of square-well and the other of hard-sphere character. The orientational constraint is enforced by assuming that each hemisphere can only point either North or South with equal probability. The analytical approach hinges on a mapping of the above Janus fluid onto a binary mixture interacting via a "quasi" isotropic potential. The anisotropic nature of the original Kern-Frenkel potential is reflected by the asymmetry in the interactions occurring between the unlike components of the mixture. A rational-function approximation extending the corresponding symmetric case is obtained in the sticky limit, where the square-well becomes infinitely narrow and deep, and allows a fully analytical approach. Notwithstanding the rather drastic approximations in the analytical theory, this is shown to provide a rather precise estimate of the structural and thermodynamical properties of the original Janus fluid.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Phase diagrams of Janus fluids with up-down constrained orientations

A class of binary mixtures of Janus fluids formed by colloidal spheres with the hydrophobic hemispheres constrained to point either up or down are studied by means of Gibbs ensemble Monte Carlo simulations and simple analytical approximations. These fluids can be experimentally realized by the application of an external static electrical field. The gas-liquid and demixing phase transitions in five specific models with different patch-patch affinities are analyzed. It is found that a gas-liquid transition is present in all the models, even if only one of the four possible patch-patch interactions is attractive. Moreover, provided the attraction between like particles is stronger than between unlike particles, the system demixes into two subsystems with different composition at sufficiently low temperatures and high densities.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The isotropic-to-nematic phase transition in hard helices: theory and simulation

We investigate the isotropic-to-nematic phase transition in systems of hard helical particles, using Onsager theory and Monte Carlo computer simulations. Motivation of this work resides in the ubiquity of the helical shape motif in many natural and synthetic polymers, as well as in the well known importance that the details of size and shape have in determining the phase behaviour and properties of (soft) condensed matter systems. We discuss the differences with the corresponding spherocylinder phase diagram and find that the helix parameters affect the phase behaviour and the existence of the nematic phase. We find that for high helicity Onsager theory significantly departs from numerical simulations even when a modifed form of the Parsons-Lee rescaling is included to account for the non-convexity of particles.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Coil-globule transition of a homopolymer chain in a square-well potential: Comparison between Monte Carlo canonical replica exchange and Wang-Landau sampling

We study the equilibrium properties of a flexible homopolymer where consecutive monomers are represented by impenetrable hard spheres that are tangent to each other, and non-consecutive monomers interact via a square-well potential. To this aim, we use both replica exchange canonical simulations and micro-canonical Wang-Landau techniques for relatively short chains, and perform a close comparative analysis of the corresponding results. These investigations are then further exploited to reproduce, at a much shorter scale and, hence, computational effort, the phase diagram previously studied with much longer chains. This opens up the possibility of improving the model and introduce specificities typical, among other examples, of protein folding.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Fluid-Fluid and Fluid-Solid transitions in the Kern-Frenkel model from Barker-Henderson thermodynamic perturbation theory

We study the Kern-Frenkel model for patchy colloids using Barker-Henderson second-order thermodynamic perturbation theory. The model describes a fluid where hard sphere particles are decorated with one patch, so that they interact via a square-well (SW) potential if they are sufficiently close one another, and if patches on each particle are properly aligned. Both the gas-liquid and fluid-solid phase coexistences are computed and contrasted against corresponding Monte-Carlo simulations results. We find that the perturbation theory describes rather accurately numerical simulations all the way from a fully covered square-well potential down to the Janus limit (half coverage). In the region where numerical data are not available (from Janus to hard-spheres), the method provides estimates of the location of the critical lines that could serve as a guideline for further efficient numerical work at these low coverages. A comparison with other techniques, such as integral equation theory, highlights the important aspect of this methodology in the present context.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

One dimensional Potts model with many-body interactions and the Generalized Model of Polypeptide Chain for the helix-coil transition

Helix-coil transition in polypeptides is an example of a spin model with a preferred spin direction, in the sense that a theoretical formulation of this problem requires to assign a preferred value of spin to the helical conformation in order to account for different symmetries of the helical {\sl vs.} the coil states. This leads to the spin Hamiltonian of the {\sl Generalized Model of Polypeptide Chain} (GMPC) variety as opposed to the Potts model variety, both with many-body interactions. We compare the explicit solution of the Potts model and the solution of the GMPC within the transfer-matrix formalism. Comparison of both secular equations reveals that the largest eigenvalue of the Potts model with $Δ$ many-body interactions is identical to the largest eigenvalue of the GMPC model with $Δ-1$ many-body interactions, indicating the equivalence of both free energies. In distinction, the second largest eigenvalues do not coincide, leading to different thermal behavior of the spatial correlation length, related to the helix-coil transition interval. Spin models with built-in spin anisotropy thus engender different physical properties in the thermodynamic limit that we explore in detail.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Osmotic pressure induced coupling between cooperativity and stability of a helix-coil transition

Most helix-coil transition theories can be characterized by a set of three parameters: energetic, describing the (free) energy cost of forming a helical state in one repeating unit; entropic, accounting for the decrease of entropy due to the helical state formation; and geometric, indicating how many repeating units are affected by the formation of one helical state. Depending on their effect on the helix-coil transition, solvents or co-solutes can be classified with respect to their action on these parameters. Solvent interactions that alter the entropic cost of helix formation by their osmotic action can affect both the stability (transition temperature) and the cooperativity (transition interval) of the helix-coil transition. A consistent inclusion of osmotic pressure effects in a description of helix-coil transition for poly(L-glutamic acid) in solution with polyethylene glycol can offer an explanation of the experimentally observed linear dependence of transition temperature on osmotic pressure as well as the concurrent changes in the cooperativity of the transition.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Competition for hydrogen bond formation in the helix-coil transition and protein folding

The problem of the helix-coil transition of biopolymers in explicit solvents, like water, with the ability for hydrogen bonding with solvent is addressed analytically using a suitably modified version of the Generalized Model of Polypeptide Chains. Besides the regular helix-coil transition, an additional coil-helix or reentrant transition is also found at lower temperatures. The reentrant transition arises due to competition between polymer-polymer and polymer-water hydrogen bonds. The balance between the two types of hydrogen bonding can be shifted to either direction through changes not only in temperature, but also by pressure, mechanical force, osmotic stress or other external influences. Both polypeptides and polynucleotides are considered within a unified formalism. Our approach provides an explanation of the experimental difficulty of observing the reentrant transition with pressure; and underscores the advantage of pulling experiments for studies of DNA. Results are discussed and compared with those reported in a number of recent publications with which a significant level of agreement is obtained.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Effect of Polydispersity and Anisotropy in Colloidal and Protein Solutions: an Integral Equation Approach

Application of integral equation theory to complex fluids is reviewed, with particular emphasis to the effects of polydispersity and anisotropy on their structural and thermodynamic properties. Both analytical and numerical solutions of integral equations are discussed within the context of a set of minimal potential models that have been widely used in the literature. While other popular theoretical tools, such as numerical simulations and density functional theory, are superior for quantitative and accurate predictions, we argue that integral equation theory still provides, as in simple fluids, an invaluable technique that is able to capture the main essential features of a complex system, at a much lower computational cost. In addition, it can provide a detailed description of the angular dependence in arbitrary frame, unlike numerical simulations where this information is frequently hampered by insufficient statistics. Applications to colloidal mixtures, globular proteins and patchy colloids are discussed, within a unified framework.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Phase diagram of the penetrable square well-model

We study a system formed by soft colloidal spheres attracting each other via a square-well potential, using extensive Monte Carlo simulations of various nature. The softness is implemented through a reduction of the infinite part of the repulsive potential to a finite one. For sufficiently low values of the penetrability parameter we find the system to be Ruelle stable with square-well like behavior. For high values of the penetrability the system is thermodynamically unstable and collapses into an isolated blob formed by a few clusters each containing many overlapping particles. For intermediate values of the penetrability the system has a rich phase diagram with a partial lack of thermodynamic consistency.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The penetrable square-well model: extensive versus non-extensive phase

The phase diagram of the penetrable square-well fluid is investigated through Monte Carlo simulations of various nature. This model was proposed as the simplest possibility of combining bounded repulsions at short scale and short-range attractions. We prove that the model is thermodynamically stable for sufficiently low values of the penetrability parameter, and in this case the system behaves similarly to the square-well model. For larger penetration, there exists an intermediate region where the system is metastable, with well defined fluid-fluid and fluid-solid transitions, at finite size, but eventually becomes unstable in the thermodynamic limit. We characterize the unstable non-extensive phase appearing at high penetrability, where the system collapses into an isolated blob of a few clusters of many ovelapping particles each.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Cluster theory of Janus particles

We apply a simple statistical mechanics cluster approximation for studying clustering in the Kern and Frenkel model of Janus fluids. The approach is motivated by recent Monte Carlo simulations work on the same model revealing that the vapor coexisting with the liquid phase contains clusters of different sizes and shapes whose equilibrium concentrations in general depend on the interaction range as well as on thermodynamic parameters. The approximation hinges on a separation between the intra- and inter-cluster contribution to thermodynamics, where only the former is explicitly computed by Monte Carlo simulations. Two levels of a simple liquid theory approximation are exploited for the description of the latter. In the first we use the ideal-gas expressions and obtain a qualitative agreement with extensive Monte Carlo bulk simulations. This can be improved to a semi-quantitative agreement, by using a hard-sphere description for the cluster-cluster correlations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Effects of patch size and number within a simple model of patchy colloids

We report on a computer simulation and integral equation study of a simple model of patchy spheres, each of whose surfaces is decorated with two opposite attractive caps, as a function of the fraction $χ$ of covered attractive surface. The simple model explored --- the two-patch Kern-Frenkel model --- interpolates between a square-well and a hard-sphere potential on changing the coverage $χ$. We show that integral equation theory provides quantitative predictions in the entire explored region of temperatures and densities from the square-well limit $χ= 1.0$ down to $χ\approx 0.6$. For smaller $χ$, good numerical convergence of the equations is achieved only at temperatures larger than the gas-liquid critical point, where however integral equation theory provides a complete description of the angular dependence. These results are contrasted with those for the one-patch case. We investigate the remaining region of coverage via numerical simulation and show how the gas-liquid critical point moves to smaller densities and temperatures on decreasing $χ$. Below $χ\approx 0.3$, crystallization prevents the possibility of observing the evolution of the line of critical points, providing the angular analog of the disappearance of the liquid as an equilibrium phase on decreasing the range for spherical potentials. Finally, we show that the stable ordered phase evolves on decreasing $χ$ from a three-dimensional crystal of interconnected planes to a two-dimensional independent-planes structure to a one-dimensional fluid of chains when the one-bond-per-patch limit is eventually reached.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Liquid-vapor coexistence in square-well fluids: an RHNC study

We investigate the ability of the reference hypernetted-chain integral equation to describe the phase diagram of square-well fluids with four different ranges of attraction. Comparison of our results with simulation data shows that the theory is able to reproduce with fairly good accuracy a significant part of the coexistence curve, provided an extrapolation procedure is used to circumvent the well-known pathologies of the pseudo-spinodal line, which are more severe at reduced width of the attractive well. The method provides a useful approach for a quick assessment of the location of the liquid-vapor coexistence curve in this kind of fluid and serves as a check for the more complex problem of anisotropic "patchy" square-well molecules.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Local orientational ordering in fluids of spherical molecules with dipolar-like anisotropic adhesion

We discuss some interesting physical features stemming from our previous analytical study of a simple model of a fluid with dipolar-like interactions of very short range in addition to the usual isotropic Baxter potential for adhesive spheres. While the isotropic part is found to rule the global structural and thermodynamical equilibrium properties of the fluid, the weaker anisotropic part gives rise to an interesting short-range local ordering of nearly spherical condensation clusters, containing short portions of chains having nose-to-tail parallel alignment which runs antiparallel to adjacent similar chains.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Phase diagram and structural properties of a simple model for one-patch particles

We study the thermodynamic and structural properties of a simple, one-patch fluid model using the reference hypernetted-chain (RHNC) integral equation and specialized Monte Carlo simulations. In this model, the interacting particles are hard spheres, each of which carries a single identical, arbitrarily-oriented, attractive circular patch on its surface; two spheres attract via a simple square-well potential only if the two patches on the spheres face each other within a specific angular range dictated by the size of the patch. For a ratio of attractive to repulsive surface of 0.8, we construct the RHNC fluid-fluid separation curve and compare with that obtained by Gibbs ensemble and grand canonical Monte Carlo simulations. We find that RHNC provides a quick and highly reliable estimate for the position of the fluid-fluid critical line. In addition, it gives a detailed (though approximate) description of all structural properties and their dependence on patch size.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

Iterated Function System and Diffusion in the Presence of Disorder and Traps

The escape probability $ξ_{x}$ from a site $x$ of a one-dimensional disordered lattice with trapping is treated as a discrete dynamical evolution by random iterations over nonlinear maps parametrized by the right and left jump probabilities. The invariant measure of the dynamics is found to be a multifractal. However the measure becomes uniform over the support when the disorder becomes weak for any non-zero trapping probability. Implications of our findings in terms of diffusion are discussed.

0 citations0 reviewsNo code linkedOpen access