Gaming self-consistent field theory: Generative block polymer phase discovery.

Block polymers are an attractive platform for uncovering the factors that give rise to self-assembly in soft matter owing to their relatively simple thermodynamic description, as captured in self-consistent field theory (SCFT). SCFT historically has found great success explaining experimental data, allowing one to construct phase diagrams from a set of candidate phases, and there is now strong interest in deploying SCFT as a screening tool to guide experimental design. However, using SCFT for phase discovery leads to a conundrum: How does one discover a new morphology if the set of candidate phases needs to be specified in advance? This long-standing challenge was surmounted by training a deep convolutional generative adversarial network (GAN) with trajectories from converged SCFT solutions, and then deploying the GAN to generate input fields for subsequent SCFT calculations. The power of this approach is demonstrated for network phase formation in neat diblock copolymer melts via SCFT. A training set of only five networks produced 349 candidate phases spanning known and previously unexplored morphologies, including a chiral network. This computational pipeline, constructed here entirely from open-source codes, should find widespread application in block polymer phase discovery and other forms of soft matter.

DNA in nanochannels: theory and applications.

Nanofluidic structures have over the last two decades emerged as a powerful platform for detailed analysis of DNA on the kilobase pair length scale. When DNA is confined to a nanochannel, the combination of excluded volume and DNA stiffness leads to the DNA being stretched to near its full contour length. Importantly, this stretching takes place at equilibrium, without any chemical modifications to the DNA. As a result, any DNA can be analyzed, such as DNA extracted from cells or circular DNA, and it is straight-forward to study reactions on the ends of linear DNA. In this comprehensive review, we first give a thorough description of the current understanding of the polymer physics of DNA and how that leads to stretching in nanochannels. We then describe how the versatility of nanofabrication can be used to design devices specifically tailored for the problem at hand, either by controlling the degree of confinement or enabling facile exchange of reagents to measure DNA-protein reaction kinetics. The remainder of the review focuses on two important applications of confining DNA in nanochannels. The first is optical DNA mapping, which provides the genomic sequence of intact DNA molecules in excess of 100 kilobase pairs in size, with kilobase pair resolution, through labeling strategies that are suitable for fluorescence microscopy. In this section, we highlight solutions to the technical aspects of genomic mapping, including the use of enzyme-based labeling and affinity-based labeling to produce the genomic maps, rather than recent applications in human genetics. The second is DNA-protein interactions, and several recent examples of such studies on DNA compaction, filamentous protein complexes, and reactions with DNA ends are presented. Taken together, these two applications demonstrate the power of DNA confinement and nanofluidics in genomics, molecular biology, and biophysics.

Diffusion of knots in nanochannel-confined DNA molecules.

We used Langevin dynamics simulations without hydrodynamic interactions to probe knot diffusion mechanisms and the time scales governing the evolution and the spontaneous untying of trefoil knots in nanochannel-confined DNA molecules in the extended de Gennes regime. The knot untying follows an "opening up process," wherein the initially tight knot continues growing and fluctuating in size as it moves toward the end of the DNA molecule before its annihilation at the chain end. The mean knot size increases significantly and sub-linearly with increasing chain contour length. The knot diffusion in nanochannel-confined DNA molecules is subdiffusive, with the unknotting time scaling with chain contour length with an exponent of 2.64 ± 0.23 to within a 95% confidence interval. The scaling exponent for the mean unknotting time vs chain contour length, along with visual inspection of the knot conformations, suggests that the knot diffusion mechanism is a combination of self-reptation and knot region breathing for the simulated parameters.

Thermodynamics and morphology of linear multiblock copolymers at homopolymer interfaces.

Block copolymers at homopolymer interfaces are poised to play a critical role in the compatibilization of mixed plastic waste, an area of growing importance as the rate of plastic accumulation rapidly increases. Using molecular dynamics simulations of Kremer-Grest polymer chains, we have investigated how the number of blocks and block degree of polymerization in a linear multiblock copolymer impacts the interface thermodynamics of strongly segregated homopolymer blends, which is key to effective compatibilization. The second virial coefficient reveals that interface thermodynamics are more sensitive to block degree of polymerization than to the number of blocks. Moreover, we identify a strong correlation between surface pressure (reduction of interfacial tension) and the spatial uniformity of block junctions on the interface, yielding a morphological framework for interpreting the role of compatibilizer architecture (number of blocks) and block degree of polymerization. These results imply that, especially at high interfacial loading, the choice of architecture of a linear multiblock copolymer compatibilizing surfactant does not greatly affect the modification of interfacial tension.

Symmetry breaking in particle-forming diblock polymer/homopolymer blends.

Compositionally asymmetric diblock copolymers provide an attractive platform for understanding the emergence of tetragonally close-packed, Frank-Kasper phases in soft matter. Block-polymer phase behavior is governed by a straightforward competition between chain stretching and interfacial tension under the constraint of filling space at uniform density. Experiments have revealed that diblock copolymers with insufficient conformational asymmetry to form Frank-Kasper phases in the neat-melt state undergo an interconversion from body-centered cubic (bcc) close-packed micelles to a succession of Frank-Kasper phases (σ to C14 to C15) upon the addition of minority-block homopolymer in the dry-brush regime, accompanied by the expected transition from bcc to hexagonally packed cylinders in the wet-brush regime. Self-consistent field theory data presented here qualitatively reproduce the salient features of the experimental phase behavior. A particle-by-particle analysis of homopolymer partitioning furnishes a basis for understanding the symmetry breaking from the high-symmetry bcc phase to the lower-symmetry Frank-Kasper phases, wherein the reconfiguration of the system into polyhedra of increasing volume asymmetry delays the onset of macroscopic phase separation.

Tuning conformational asymmetry in particle-forming diblock copolymer alloys.

Self-consistent field theory is employed to compute the phase behavior of binary blends of conformationally asymmetric, micelle-forming diblock copolymers with miscible corona blocks and immiscible core blocks (a diblock copolymer "alloy"). The calculations focus on establishing conditions that promote the formation of Laves phases by tuning the relative softness of the cores of the two different Laves phase particles via independent control of their conformational asymmetries. Increasing the conformational asymmetry of the more spherical particles of the Laves structure has a stabilizing effect, consistent with the expectations of increased imprinting of the Wigner-Seitz cells on the core/corona interface as conformational asymmetry increases. The resulting phase diagram in the temperature-blend composition space features a more stable Laves phase field than that predicted for conformationally symmetric systems. The phase field closes at low temperatures in favor of macrophase separation between a hexagonally-packed cylinder (hex) phase and a body-centered cubic phase. Companion calculations, using an alloy whose components do not produce a hex phase in the neat melt state, suggest that the Laves phase field in such a blend will persist at strong segregation.

Hydrophobic catalysis and a potential biological role of DNA unstacking induced by environment effects.

Hydrophobic base stacking is a major contributor to DNA double-helix stability. We report the discovery of specific unstacking effects in certain semihydrophobic environments. Water-miscible ethylene glycol ethers are found to modify structure, dynamics, and reactivity of DNA by mechanisms possibly related to a biologically relevant hydrophobic catalysis. Spectroscopic data and optical tweezers experiments show that base-stacking energies are reduced while base-pair hydrogen bonds are strengthened. We propose that a modulated chemical potential of water can promote "longitudinal breathing" and the formation of unstacked holes while base unpairing is suppressed. Flow linear dichroism in 20% diglyme indicates a 20 to 30% decrease in persistence length of DNA, supported by an increased flexibility in single-molecule nanochannel experiments in poly(ethylene glycol). A limited (3 to 6%) hyperchromicity but unaffected circular dichroism is consistent with transient unstacking events while maintaining an overall average B-DNA conformation. Further information about unstacking dynamics is obtained from the binding kinetics of large thread-intercalating ruthenium complexes, indicating that the hydrophobic effect provides a 10 to 100 times increased DNA unstacking frequency and an "open hole" population on the order of 10-2 compared to 10-4 in normal aqueous solution. Spontaneous DNA strand exchange catalyzed by poly(ethylene glycol) makes us propose that hydrophobic residues in the L2 loop of recombination enzymes RecA and Rad51 may assist gene recombination via modulation of water activity near the DNA helix by hydrophobic interactions, in the manner described here. We speculate that such hydrophobic interactions may have catalytic roles also in other biological contexts, such as in polymerases.

The C36 Laves phase in diblock polymer melts.

The C14 and C15 Laves phases form as micelle packing structures in many types of soft matter, but the related C36 phase, which consists of alternating C14-type and C15-type layers, has not been observed in any such system. To understand this absence in the context of diblock polymers, we used self-consistent field theory to relate the morphology and energetics of C36 to other known mesophases. Two case studies were conducted: blends of AB diblock polymers with A homopolymers (where A forms the micelle core), in which C14 and C15 have stability windows, and neat AB diblock melts, in which Laves phases are metastable. Laves phases exhibit nearly identical micelle morphologies and nearly degenerate free energies, with the free energy of C36 being a near-perfect bisector of the C14 and C15 free energies in all cases, revealing an intrinsic symmetry in free energy that is attributed solely to the structural relationship between the phases in which the packing of C36 is intermediate between C14 and C15. Based on this connection between structure and free energy, C36 is thus not expected to form in flexible diblock polymers, since C14 and C15 can always form instead via facile mass transfer.

Open-source platform for block polymer formulation design using particle swarm optimization.

Facile exploration of large design spaces is critical to the development of new functional soft materials, including self-assembling block polymers, and computational inverse design methodologies are a promising route to initialize this task. We present here an open-source software package coupling particle swarm optimization (PSO) with an existing open-source self-consistent field theory (SCFT) software for the inverse design of self-assembling block polymers to target bulk morphologies. To lower the barrier to use of the software and facilitate exploration of novel design spaces, the underlying SCFT calculations are seeded with algorithmically generated initial fields for four typical morphologies: lamellae, network phases, cylindrical phases, and spherical phases. In addition to its utility within PSO, the initial guess tool also finds generic applicability for stand-alone SCFT calculations. The robustness of the software is demonstrated with two searches for classical phases in the conformationally symmetric diblock system, as well as one search for the Frank-Kasper [Formula: see text] phase in conformationally asymmetric diblocks. The source code for both the initial guess generation and the PSO wrapper is publicly available.

Interactions between two knots in nanochannel-confined DNA molecules.

Experimental data on the interaction between two knots in deoxyribonucleic acid (DNA) confined in nanochannels produced two particular behaviors of knot pairs along the DNA molecules: (i) widely separated knots experience an attractive interaction but only remain in close proximity for several seconds and (ii) knots tend to remain separated until one of the knots unravels at the chain end. The associated free energy profile of the knot-knot separation distance for an ensemble of DNA knots exhibits a global minimum when knots are separated, indicating that the separated knot state is more stable than the intertwined knot state, with dynamics in the separated knot state that are consistent with independent diffusion. The experimental observations of knot-knot interactions under nanochannel confinement are inconsistent with previous simulation-based and experimental results for stretched polymers under tension wherein the knots attract and then stay close to each other. This inconsistency is postulated to result from a weaker fluctuation-induced attractive force between knots under confinement when compared to the knots under tension, the latter of which experience larger fluctuations in transverse directions.

Stable Frank-Kasper phases of self-assembled, soft matter spheres.

Single molecular species can self-assemble into Frank-Kasper (FK) phases, finite approximants of dodecagonal quasicrystals, defying intuitive notions that thermodynamic ground states are maximally symmetric. FK phases are speculated to emerge as the minimal-distortional packings of space-filling spherical domains, but a precise measure of this distortion and how it affects assembly thermodynamics remains ambiguous. We use two complementary approaches to demonstrate that the principles driving FK lattice formation in diblock copolymers emerge directly from the strong-stretching theory of spherical domains, in which a minimal interblock area competes with a minimal stretching of space-filling chains. The relative stability of FK lattices is studied first using a diblock foam model with unconstrained particle volumes and shapes, which correctly predicts not only the equilibrium σ lattice but also the unequal volumes of the equilibrium domains. We then provide a molecular interpretation for these results via self-consistent field theory, illuminating how molecular stiffness increases the sensitivity of the intradomain chain configurations and the asymmetry of local domain packing. These findings shed light on the role of volume exchange on the formation of distinct FK phases in copolymers and suggest a paradigm for formation of FK phases in soft matter systems in which unequal domain volumes are selected by the thermodynamic competition between distinct measures of shape asymmetry.

Origins of low-symmetry phases in asymmetric diblock copolymer melts.

Cooling disordered compositionally asymmetric diblock copolymers leads to the formation of nearly spherical particles, each containing hundreds of molecules, which crystallize upon cooling below the order-disorder transition temperature (TODT). Self-consistent field theory (SCFT) reveals that dispersity in the block degrees of polymerization stabilizes various Frank-Kasper phases, including the C14 and C15 Laves phases, which have been accessed experimentally in low-molar-mass poly(isoprene)-b-poly(lactide) (PI-PLA) diblock copolymers using thermal processing strategies. Heating and cooling a specimen containing 15% PLA above and below the TODT from the body-centered cubic (BCC) or C14 states regenerates the same crystalline order established at lower temperatures. This memory effect is also demonstrated with a specimen containing 20% PLA, which recrystallizes to either C15 or hexagonally ordered cylinders (HEXC) upon heating and cooling. The process-path-dependent formation of crystalline order shapes the number of particles per unit volume, n/V, which is retained in the highly structured disordered liquid as revealed by small-angle X-ray scattering (SAXS) experiments. We hypothesize that symmetry breaking during crystallization is governed by the particle number density imprinted in the liquid during ordering at lower temperature, and this metastable liquid is kinetically constrained from equilibrating due to prohibitively large free energy barriers for micelle fusion and fission. Ordering at fixed n/V is enabled by facile chain exchange, which redistributes mass as required to meet the multiple particle sizes and packing associated with specific low-symmetry Frank-Kasper phases. This discovery exposes universal concepts related to order and disorder in self-assembled soft materials.

Open-source code for self-consistent field theory calculations of block polymer phase behavior on graphics processing units.

Self-consistent field theory (SCFT) is a powerful approach for computing the phase behavior of block polymers. We describe a fast version of the open-source Polymer Self-Consistent Field (PSCF) code that takes advantage of the massive parallelization provided by a graphical processing unit (GPU). Benchmarking double-precision calculations indicate up to 30× reduction in time to converge SCFT calculations of various diblock copolymer phases when compared to the Fortran CPU version of PSCF using the same algorithms, with the speed-up increasing with increasing unit cell size for the diblock polymer problems examined here. Where double-precision accuracy is not needed, single-precision calculations can provide speed-up of up to 60× in convergence time. These improvements in speed within an open-source format open up new vistas for SCFT-driven block polymer materials discovery by the community at large.

Influence of charge sequence on the adsorption of polyelectrolytes to oppositely-charged polyelectrolyte brushes.

When a solution of polyanionic chains is placed in contact with a polycationic brush, the polyanions adsorb into the brush. We investigate the influence of the charge sequences of the free and bound species on the thermodynamics of polyelectrolyte adsorption. As model systems, we consider free and brush polyelectrolytes with either block or alternating charge sequences, and study the adsorption process using coarse-grained Langevin dynamics with implicit solvent, explicit counterions, and excess salt. Free energy, internal energy, and entropy of adsorption are computed using umbrella sampling methods. When the number of polyanions exceed the number of polycations, the brush becomes overcharged. Free chains adsorb most strongly when both free and tethered chains have a block charge sequence, and most weakly when both species have an alternating sequence. Adsorption is stronger when the free polyanion has a block sequence and the tethered polycation is alternating than in the reverse case of an alternating free polymer and a tethered block copolymer. Sequence-dependent effects are shown to be largely energetic, rather than entropic, in origin.

Extension distribution for DNA confined in a nanochannel near the Odijk regime.

DNA confinement in a nanochannel typically is understood via mapping to the confinement of an equivalent neutral polymer by hard walls. This model has proven to be effective for confinement in relatively large channels where hairpin formation is frequent. An analysis of existing experimental data for Escherichia coli DNA extension in channels smaller than the persistence length, combined with an additional dataset for λ-DNA confined in a 34 nm wide channel, reveals a breakdown in this approach as the channel size approaches the Odijk regime of strong confinement. In particular, the predicted extension distribution obtained from the asymptotic solution to the weakly correlated telegraph model for a confined wormlike chain deviates significantly from the experimental distribution obtained for DNA confinement in the 34 nm channel, and the discrepancy cannot be resolved by treating the alignment fluctuations or the effective channel size as fitting parameters. We posit that the DNA-wall electrostatic interactions, which are sensible throughout a significant fraction of the channel cross section in the Odijk regime, are the source of the disagreement between theory and experiment. Dimensional analysis of the wormlike chain propagator in channel confinement reveals the importance of a dimensionless parameter, reflecting the magnitude of the DNA-wall electrostatic interactions relative to thermal energy, which has not been considered explicitly in the prevailing theories for DNA confinement in a nanochannel.

Role of Chain Length in the Formation of Frank-Kasper Phases in Diblock Copolymers.

The phase behavior of poly(styrene)-b-poly(1,4-butadiene) diblock copolymers with a polymer block invariant degree of polymerization N[over ¯]_{b}≈800 shows no evidence of Frank-Kasper phases, in contrast to low molar mass diblock copolymers (N[over ¯]_{b}<100) with the same conformational asymmetry. A universal self-concentration crossover parameter N[over ¯]_{x}≈400 is identified, directly related to the crossover to entanglement dynamics in polymer melts. Mean-field behavior is recovered when N[over ¯]_{b}>N[over ¯]_{x}, while complex low symmetry phase formation is attributed to fluctuations and space-filling constraints, which dominate when N[over ¯]_{b}

Alpha-Synuclein Modulates the Physical Properties of DNA.

Fundamental research on Parkinson's disease (PD) most often focuses on the ability of α-synuclein (aS) to form oligomers and amyloids, and how such species promote brain cell death. However, there are indications that aS also plays a gene-regulatory role in the cell nucleus. Here, the interaction between monomeric aS and DNA in vitro has been investigated with single-molecule techniques. Using a nanofluidic channel system, it was discovered that aS binds to DNA and by studying the DNA-protein complexes at different confinements we determined that aS binding increases the persistence length of DNA from 70 to 90 nm at high coverage. By atomic force microscopy it was revealed that at low protein-to-DNA ratio, the aS binding occurs as small protein clusters scattered along the DNA; at high protein-to-DNA ratio, the DNA is fully covered by protein. As DNA-aS interactions may play roles in PD, it is of importance to characterize biophysical properties of such complexes in detail.

Measuring the wall depletion length of nanoconfined DNA.

Efforts to study the polymer physics of DNA confined in nanochannels have been stymied by a lack of consensus regarding its wall depletion length. We have measured this quantity in 38 nm wide, square silicon dioxide nanochannels for five different ionic strengths between 15 mM and 75 mM. Experiments used the Bionano Genomics Irys platform for massively parallel data acquisition, attenuating the effect of the sequence-dependent persistence length and finite-length effects by using nick-labeled E. coli genomic DNA with contour length separations of at least 30 µm (88 325 base pairs) between nick pairs. Over 5 × 106 measurements of the fractional extension were obtained from 39 291 labeled DNA molecules. Analyzing the stretching via Odijk's theory for a strongly confined wormlike chain yielded a linear relationship between the depletion length and the Debye length. This simple linear fit to the experimental data exhibits the same qualitative trend as previously defined analytical models for the depletion length but now quantitatively captures the experimental data.

The Statistical Segment Length of DNA: Opportunities for Biomechanical Modeling in Polymer Physics and Next-Generation Genomics.

The development of bright bisintercalating dyes for deoxyribonucleic acid (DNA) in the 1990s, most notably YOYO-1, revolutionized the field of polymer physics in the ensuing years. These dyes, in conjunction with modern molecular biology techniques, permit the facile observation of polymer dynamics via fluorescence microscopy and thus direct tests of different theories of polymer dynamics. At the same time, they have played a key role in advancing an emerging next-generation method known as genome mapping in nanochannels. The effect of intercalation on the bending energy of DNA as embodied by a change in its statistical segment length (or, alternatively, its persistence length) has been the subject of significant controversy. The precise value of the statistical segment length is critical for the proper interpretation of polymer physics experiments and controls the phenomena underlying the aforementioned genomics technology. In this perspective, we briefly review the model of DNA as a wormlike chain and a trio of methods (light scattering, optical or magnetic tweezers, and atomic force microscopy (AFM)) that have been used to determine the statistical segment length of DNA. We then outline the disagreement in the literature over the role of bisintercalation on the bending energy of DNA, and how a multiscale biomechanical approach could provide an important model for this scientifically and technologically relevant problem.