Some comments on Monte Carlo and molecular dynamics methods

We highlight some links between molecular dynamics and Monte Carlo algorithms used to simulate condensed matter systems. Special attention is paid to the question of sampling the desired statistical ensemble

Folding kinetics of a polymer [corrigendum]

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

Enhancement of island size by dynamic substrate disorder in simulations of graphene growth

We demonstrate a new mechanism in the early stages of sub-monolayer epitaxial island growth, using Monte Carlo simulations motivated by experimental observations on the growth of graphene on copper foil. In our model, the substrate is “dynamically rough”, by which we mean (i) the interaction strength between Cu and C varies randomly from site to site, and (ii) these variable strengths themselves migrate from site to site. The dynamic roughness provides a simple representation of the near-molten state of the Cu substrate in the case of real graphene growth. Counterintuitively, the graphene island size increases when dynamic roughness is included, compared to a static and smooth substrate. We attribute this effect to destabilisation of small graphene islands by fluctuations in the substrate, allowing them to break up and join larger islands which are more stable against roughness. In the case of static roughness, when process (ii) is switched off, island growth is strongly inhibited and the scale-free behaviour of island size distributions, present in the smooth-static and rough-dynamic cases, is destroyed. The effects of the dynamic substrate roughness cannot be mimicked by parameter changes in the static cases

Testing the transferability of a coarse-grained model to intrinsically disordered proteins

The intermediate-resolution coarse-grained protein model PLUM [T. Bereau and M. Deserno, J. Chem. Phys., 2009, 130, 235106] is used to simulate small systems of intrinsically disordered proteins involved in biomineralisation. With minor adjustments to reduce bias toward stable secondary structure, the model generates conformational ensembles conforming to structural predictions from atomistic simulation. Without additional structural information as input, the model distinguishes regions of the chain by predicted degree of disorder, manifestation of structure, and involvement in chain dimerisation. The model is also able to distinguish dimerisation behaviour between one intrinsically disordered peptide and a closely related mutant. We contrast this against the poor ability of PLUM to model the S1 quartz-binding peptide

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

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

Unbiased molecular analysis of T cell receptor expression using template-switch achored RT-PCR

A detailed knowledge of the principles that guide clonal selection within the memory and effector T cell pools is essential to further our understanding of the factors that influence effective T cell-mediated immunity and has direct implications for the rational design of vaccines and immunotherapies. This unit provides methods for the unbiased quantification and characterization of all expressed T cell receptor (TCR) gene products within any defined T cell population. The approach is based on a template-switch anchored reverse transcription–polymerase chain reaction (RT-PCR) and is optimized for the analysis of antigen-specific T cells isolated directly ex vivo

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

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

Who participates in local government? Evidence from meeting minutes

Scholars and policymakers have highlighted institutions that enable community participation as a potential buffer against existing political inequalities. Yet these venues may bias policy discussions in favor of an unrepresentative group of individuals. To explore who participates, we compile a novel data set by coding thousands of instances of citizens speaking at planning and zoning board meetings concerning housing development. We match individuals to a voter file to investigate local political participation in housing and development policy. We find that individuals who are older, male, longtime residents, voters in local elections, and homeowners are significantly more likely to participate in these meetings. These individuals overwhelmingly (and to a much greater degree than the general public) oppose new housing construction. These participatory inequalities have important policy implications and may be contributing to rising housing costs.Accepted manuscrip

A gp41 MPER-specific llama VHH requires a hydrophobic CDR3 for neutralization but not for antigen recognition

The membrane proximal external region (MPER) of the HIV-1 glycoprotein gp41 is targeted by the broadly neutralizing antibodies 2F5 and 4E10. To date, no immunization regimen in animals or humans has produced HIV-1 neutralizing MPER-specific antibodies. We immunized llamas with gp41-MPER proteoliposomes and selected a MPER-specific single chain antibody (VHH), 2H10, whose epitope overlaps with that of mAb 2F5. Bi-2H10, a bivalent form of 2H10, which displayed an approximately 20-fold increased affinity compared to the monovalent 2H10, neutralized various sensitive and resistant HIV-1 strains, as well as SHIV strains in TZM-bl cells. X-ray and NMR analyses combined with mutagenesis and modeling revealed that 2H10 recognizes its gp41 epitope in a helical conformation. Notably, tryptophan 100 at the tip of the long CDR3 is not required for gp41 interaction but essential for neutralization. Thus bi-2H10 is an anti-MPER antibody generated by immunization that requires hydrophobic CDR3 determinants in addition to epitope recognition for neutralization similar to the mode of neutralization employed by mAbs 2F5 and 4E10

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

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