Self-Consistent Field Theory of Multiply-Branched Block Copolymer Melts

We present a numerical algorithm to evaluate the self-consistent field theory for melts composed of block copolymers with multiply-branched architecture. We present results for the case of branched copolymers with doubly-functional groups for multiple branching generations. We discuss the stability of the cubic phase of spherical micelles, the A15 phase, as a consequence of tendency of the AB interfaces to conform to the polyhedral environment of the Voronoi cell of the micelle lattice.Comment: 12 pages, 10 includes figure

A format for phylogenetic placements

We have developed a unified format for phylogenetic placements, that is, mappings of environmental sequence data (e.g. short reads) into a phylogenetic tree. We are motivated to do so by the growing number of tools for computing and post-processing phylogenetic placements, and the lack of an established standard for storing them. The format is lightweight, versatile, extensible, and is based on the JSON format which can be parsed by most modern programming languages. Our format is already implemented in several tools for computing and post-processing parsimony- and likelihood-based phylogenetic placements, and has worked well in practice. We believe that establishing a standard format for analyzing read placements at this early stage will lead to a more efficient development of powerful and portable post-analysis tools for the growing applications of phylogenetic placement.Comment: Documents version 3 of the forma

Strong-Segregation Theory of Bicontinuous Phases in Block Copolymers

We compute phase diagrams for $A_nB_m$ starblock copolymers in the strong-segregation regime as a function of volume fraction $\phi$, including bicontinuous phases related to minimal surfaces (G, D, and P surfaces) as candidate structures. We present the details of a general method to compute free energies in the strong segregation limit, and demonstrate that the gyroid G phase is the most nearly stable among the bicontinuous phases considered. We explore some effects of conformational asymmetry on the topology of the phase diagram.Comment: 14 pages, latex, 21 figures, to appear in Macromolecule

Using genotype abundance to improve phylogenetic inference

Modern biological techniques enable very dense genetic sampling of unfolding evolutionary histories, and thus frequently sample some genotypes multiple times. This motivates strategies to incorporate genotype abundance information in phylogenetic inference. In this paper, we synthesize a stochastic process model with standard sequence-based phylogenetic optimality, and show that tree estimation is substantially improved by doing so. Our method is validated with extensive simulations and an experimental single-cell lineage tracing study of germinal center B cell receptor affinity maturation

Interfaces in Diblocks: A Study of Miktoarm Star Copolymers

We study AB$_n$ miktoarm star block copolymers in the strong segregation limit, focussing on the role that the AB interface plays in determining the phase behavior. We develop an extension of the kinked-path approach which allows us to explore the energetic dependence on interfacial shape. We consider a one-parameter family of interfaces to study the columnar to lamellar transition in asymmetric stars. We compare with recent experimental results. We discuss the stability of the A15 lattice of sphere-like micelles in the context of interfacial energy minimization. We corroborate our theory by implementing a numerically exact self-consistent field theory to probe the phase diagram and the shape of the AB interface.Comment: 12 pages, 11 included figure

The non-centrosymmetric lamellar phase in blends of ABC triblock and ac diblock copolymers

The phase behaviour of blends of ABC triblock and ac diblock copolymers is examined using self-consistent field theory. Several equilibrium lamellar structures are observed, depending on the volume fraction of the diblocks, phi_2, the monomer interactions, and the degrees of polymerization of the copolymers. For segregations just above the order-disorder transition the triblocks and diblocks mix together to form centrosymmetric lamellae. As the segregation is increased the triblocks and diblocks spatially separate either by macrophase-separating, or by forming a non-centrosymmetric (NCS) phase of alternating layers of triblock and diblock (...ABCcaABCca...). The NCS phase is stable over a narrow region near phi_2=0.4. This region is widest near the critical point on the phase coexistence curve and narrows to terminate at a triple point at higher segregation. Above the triple point there is two-phase coexistence between almost pure triblock and diblock phases. The theoretical phase diagram is consistent with experiments.Comment: 9 pages, 8 figures, submitted to Macromolecule

Reconciling taxonomy and phylogenetic inference: formalism and algorithms for describing discord and inferring taxonomic roots

Although taxonomy is often used informally to evaluate the results of phylogenetic inference and find the root of phylogenetic trees, algorithmic methods to do so are lacking. In this paper we formalize these procedures and develop algorithms to solve the relevant problems. In particular, we introduce a new algorithm that solves a "subcoloring" problem for expressing the difference between the taxonomy and phylogeny at a given rank. This algorithm improves upon the current best algorithm in terms of asymptotic complexity for the parameter regime of interest; we also describe a branch-and-bound algorithm that saves orders of magnitude in computation on real data sets. We also develop a formalism and an algorithm for rooting phylogenetic trees according to a taxonomy. All of these algorithms are implemented in freely-available software.Comment: Version submitted to Algorithms for Molecular Biology. A number of fixes from previous versio

New Method for Phase transitions in diblock copolymers: The Lamellar case

A new mean-field type theory is proposed to study order-disorder transitions (ODT) in block copolymers. The theory applies to both the weak segregation (WS) and the strong segregation (SS) regimes. A new energy functional is proposed without appealing to the random phase approximation (RPA). We find new terms unaccounted for within RPA. We work out in detail transitions to the lamellar state and compare the method to other existing theories of ODT and numerical simulations. We find good agreements with recent experimental results and predict that the intermediate segregation regime may have more than one scaling behavior.Comment: 23 pages, 8 figure

Representing and extending ensembles of parsimonious evolutionary histories with a directed acyclic graph

In many situations, it would be useful to know not just the best phylogenetic tree for a given data set, but the collection of high-quality trees. This goal is typically addressed using Bayesian techniques, however, current Bayesian methods do not scale to large data sets. Furthermore, for large data sets with relatively low signal one cannot even store every good tree individually, especially when the trees are required to be bifurcating. In this paper, we develop a novel object called the "history subpartition directed acyclic graph" (or "history sDAG" for short) that compactly represents an ensemble of trees with labels (e.g. ancestral sequences) mapped onto the internal nodes. The history sDAG can be built efficiently and can also be efficiently trimmed to only represent maximally parsimonious trees. We show that the history sDAG allows us to find many additional equally parsimonious trees, extending combinatorially beyond the ensemble used to construct it. We argue that this object could be useful as the "skeleton" of a more complete uncertainty quantification.Comment: To appear in JM