Source author record

Paul van der Schoot

Paul van der Schoot appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

24works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

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

24 published item(s)

preprint2026arXiv

Chirality Amplification and Chiral Segregation in Liquid Crystals

Liquid crystal mesophases of achiral molecules are normally achiral, yet in a few materials they spontaneously segregate and form right- and left-handed chiral domains. One mechanism that drives chiral segregation is molecular shape fluctuations between axial chiral conformations, where molecular interactions favor matching chirality and promote helical twist. Cooperative chiral ordering may also play a role in chirality amplification, as when a tiny fraction of chiral dopant drives a nematic phase to become cholesteric. We present a model of cooperative chiral ordering in liquid crystals using Maier-Saupe theory, and predict a phase diagram with a segregated cholesteric phase with alternating domains of left- and right-handed chiral twist, with opposite enantiomeric excess, in addition to racemic nematic and isotropic phases. Our model also demonstrates how chiral molecular fluctuations influence the helical twisting power of dopants in the nematic phase, which may be observed even in materials where the segregated cholesteric phase is preempted by a transition to another phase. We compare these results with Monte Carlo simulation studies of the switchable chiral Lebwohl-Lasher model, where each spin switches between right- and left-handed chiral states. Simulation results validate the predicted phase diagram, demonstrate chiral amplification in the racemic nematic phase, and reveal complex coarsening dynamics in the segregated cholesteric phase. These results suggest that molecular fluctuations between degenerate chiral configurations may be a common mechanism to produce cooperative chiral order in achiral liquid crystals.

preprint2022arXiv

Geometric percolation of colloids in shear flow

We combine a heuristic theory of geometric percolation and the Smoluchowski theory of colloid dynamics to predict the impact of shear flow on the percolation threshold of hard spherical colloidal particles, and verify our findings by means of molecular dynamics simulations. It appears that the impact of shear flow is subtle and highly non-trivial, even in the absence of hydrodynamic interactions between the particles. The presence of shear flow can both increase and decrease the percolation threshold, depending on the criterion used for determining whether or not two particles are connected and on the Péclet number. Our approach opens up a route to quantitatively predict the percolation threshold in nanocomposite materials that, as a rule, are produced under non-equilibrium conditions, making comparison with equilibrium percolation theory tenuous. Our theory can be adapted straightforwardly for application in other types of flow field, and particles of different shape or interacting via other than hard-core potentials.

preprint2022arXiv

The dynamics of viruslike capsid assembly and disassembly

Cowpea chlorotic mottle virus (CCMV) is a widely used model for virus replication studies. A major challenge lies in distinguishing between the roles of the interaction between coat proteins and that between the coat proteins and the viral RNA in assembly and disassembly processes. Here, we report on the spontaneous and reversible size conversion of the empty capsids of a CCMV capsid protein functionalized with a hydrophobic elastin-like polypeptide which occurs following a pH jump. We monitor the concentration of T = 3 and T = 1 capsids as a function of time and show that the time evolution of the conversion from one T number to another is not symmetric: the conversion from T = 1 to T = 3 is a factor of 10 slower than that of T = 3 to T = 1. We explain our experimental findings using a simple model based on classical nucleation theory applied to virus capsids, in which we account for the change in the free protein concentration, as the different types of shells assemble and disassemble by shedding or absorbing single protein subunits. As far as we are aware, this is the first study confirming that both the assembly and disassembly of viruslike shells can be explained through classical nucleation theory, reproducing quantitatively results from time-resolved experiments.

preprint2022arXiv

Transient response and domain formation in electrically deforming liquid crystal networks

Recently, Van der Kooij and co-workers recognised three distinct, transient regimes in the dynamics of electrically-deforming liquid crystal networks [Van der Kooij et al., Nat. Commun. 10, 1 (2019)]. Based on a Landau-theoretical framework, which encompasses spatially resolved information, we interpret these regimes: initially, the response is dominated by thermal noise, then the top of the film expands, followed by a permeation of this response into the bulk. An important signature of this interpretation is a significant dependence of the regime time scales on film thickness, where we observe a clear thin-film-to-bulk transition. The point of transition coincides with the emergence of spatial inhomogeneities in the bulk, i.e., domain formation, and should be avoided due to the less predictable steady-state expansion it gives rise to. Finally, we show that this domain formation can be suppressed by decreasing the initial thickness of the film, and increasing the linear dimensions of the mesogens, or their orientational order when crosslinked into the network, though this comes at the cost of the deformation magnitude. Our results contribute to achieving finer control over how smart liquid crystal network coatings are activated.

preprint2021arXiv

Enhanced ordering in length-polydisperse carbon nanotube solutions at high concentrations as revealed by small angle X-ray scattering

Carbon nanotubes (CNTs) are stiff, all-carbon macromolecules with diameters as small as one nanometer and few microns long. Solutions of CNTs in chlorosulfonic acid (CSA) follow the phase behavior of rigid rod polymers interacting via a repulsive potential and display a liquid crystalline phase at sufficiently high concentration. Here, we show that small-angle X-ray scattering and polarized light microscopy data can be combined to characterize quantitatively the morphology of liquid crystalline phases formed in CNT solutions at concentrations from 3 to 6.5 % by volume. We find that upon increasing their concentration, CNTs self-assemble into a liquid crystalline phase with a pleated texture and with a large inter-particle spacing that could be indicative of a transition to higher-order liquid crystalline phases. We explain how thermal undulations of CNTs can enhance their electrostatic repulsion and increase their effective diameter by an order of magnitude. By calculating the critical concentration, where the mean amplitude of undulation of an unconstrained rod becomes comparable to the rod spacing, we find that thermal undulations start to affect steric forces at concentrations as low as the isotropic cloud point in CNT solutions.

preprint2019arXiv

Geometric percolation of hard nanorods: the interplay of spontaneous and externally induced uniaxial particle alignment

We present a numerical study on geometric percolation in liquid dispersions of hard slender colloidal particles subjected to an external orienting field. In the formulation and liquid-state processing of nanocomposite materials, the alignment of particles by external fields such as electric, magnetic or flow fields is practically inevitable, and often works against the emergence of large nanoparticle networks. Using continuum percolation theory in conjunction with Onsager theory, we investigate how the interplay between externally induced alignment and the spontaneous symmetry breaking of the uniaxial nematic phase affects cluster formation within nanoparticle dispersions. It is known that the enhancement of particle alignment by means of a density increase or an external field may result in the breakdown of an already percolating network. As a result, percolation can be limited to a small region of the phase diagram only. Here, we demonstrate that the existence and shape of such a "percolation island" in the phase diagram crucially depends on the connectivity length -- a critical distance defining direct connections between neighbouring particles. Deformations of this percolation island can lead to peculiar re-entrance effects, in which a system-spanning network forms and breaks down multiple times with increasing particle density.

preprint2019arXiv

Unusual geometric percolation of hard nanorods in the uniaxial nematic liquid crystalline phase

We investigate by means of continuum percolation theory and Monte Carlo simulations how spontaneous uniaxial symmetry breaking affects geometric percolation in dispersions of hard rod-like particles. If the particle aspect ratio exceeds about twenty, percolation in the nematic phase can be lost upon adding particles to the dispersion. This contrasts with percolation in the isotropic phase, where a minimum particle loading is always required to obtain system-spanning clusters. For sufficiently short rods, percolation in the uniaxial nematic mimics that of the isotropic phase, where the addition of particles always aids percolation. For aspect ratios between twenty and infinity, but not including infinity, we find re-entrance behavior: percolation in the low-density nematic may be lost upon increasing the amount of nanofillers but can be re-gained by the addition of even more particles to the suspension. Our simulation results for aspect ratios of 5, 10, 20, 50 and 100 strongly support our theoretical predictions, with almost quantitative agreement. We show that a new closure of the connectedness Ornstein-Zernike equation, inspired by Scaled Particle Theory, is more accurate than the Lee-Parsons closure that effectively describes the impact of many-body direct contacts.

preprint2016arXiv

A Landau-de Gennes theory for hard colloidal rods: defects and tactoids

We construct a phenomenological Landau-de Gennes theory for hard colloidal rods by performing an order parameter expansion of the chemical-potential dependent grand potential. By fitting the coefficients to known results of Onsager theory, we are not only able to describe the isotropic-nematic phase transition as function of density, including the well-known density jump, but also the isotropic-nematic planar interface. The resulting theory is applied in calculations of the isotropic core size in a radial hedgehog defect, the density dependence of linear defects of hard rods in square confinement, and the formation of a nematic droplet in an isotropic background.

preprint2016arXiv

Dynamic Landau Theory for Supramolecular Self-Assembly

Although pathway-specific kinetic theories are fundamentally important to describe and under- stand reversible polymerisation kinetics, they come in principle at a cost of having a large number of system-specific parameters. Here, we construct a dynamical Landau theory to describe the kinetics of activated linear supramolecular self-assembly, which drastically reduces the number of parameters and still describes most of the interesting and generic behavior of the system in hand. This phenomenological approach hinges on the fact that if nucleated, the polymerisation transition resembles a phase transition. We are able to describe hysteresis, overshooting, undershooting and the existence of a lag time before polymerisation takes o?, and pinpoint the conditions required for observing these types of phenomenon in the assembly and disassembly kinetics. We argue that the phenomenological kinetic parameter in our theory is a pathway controller, i.e., it controls the relative weights of the molecular pathways through which self-assembly takes place.

preprint2016arXiv

Effects of RNA branching on the electrostatic stabilization of viruses

Many single-stranded (ss) RNA viruses self assemble from capsid protein subunits and the nucleic acid to form an infectious virion. It is believed that the electrostatic interactions between the negatively charged RNA and the positively charged viral capsid proteins drive the encapsidation, although there is growing evidence that the sequence of the viral RNA also plays a role in packaging. In particular the sequence will determine the possible secondary structures that the ssRNA will take in solution. In this work, we use a mean field theory to investigate how the secondary structure of the RNA combined with electrostatic interactions affects the efficiency of assembly and stability of the assembled virions. We show that the secondary structure of RNA may result in negative osmotic pressures while a linear polymer causes positive osmotic pressures for the same conditions. This may suggest that the branched structure makes the RNA more effectively packaged and the virion more stable.

preprint2016arXiv

RNA topology remolds electrostatic stabilization of viruses

Simple RNA viruses efficiently encapsulate their genome into a nano-sized protein shell: the capsid. Spontaneous co-assembly of the genome and the capsid proteins is driven predominantly by electrostatic interactions between the negatively charged RNA and the positively charged inner capsid wall. Using field theoretic formulation we show that the inherently branched RNA secondary structure allows viruses to maximize the amount of encapsulated genome and make assembly more efficient, allowing viral RNAs to out-compete cellular RNAs during replication in infected host cells.

preprint2016arXiv

Role of Genome in the Formation of Conical Retroviral Shells

Human immunodeficiency virus (HIV) capsid proteins spontaneously assemble around the genome into a protective protein shell called the capsid, which can take on a variety of shapes broadly classified as conical, cylindrical and irregular. The majority of capsids seen in in vivo studies are conical in shape, while in vitro experiments have shown a preference for cylindrical capsids. The factors involved in the selection of the unique shape of HIV capsids are not well understood, and in particular the impact of RNA on the formation of the capsid is not known. In this work, we study the role of the genome and its interaction with the capsid protein by modeling the genomic RNA through a mean-field theory. Our results show that the confinement free energy for a homopolymeric model genome confined in a conical capsid is lower than that in a cylindrical capsid, at least when the genome does not interact with the capsid, which seems to be the case in in vivo experiments. Conversely, the confinement free energy for the cylinder is lower than for a conical capsid if the genome is attracted to the capsid proteins as the in vitro experiments. Understanding the factors that contribute to the formation of conical capsids may shed light on the infectivity of HIV particles.

preprint2016arXiv

Stochastic Lag Time in Nucleated Linear Self-Assembly

Protein aggregation is of great importance in biology, e.g., in amyloid fibrillation. The aggregation processes that occur at the cellular scale must be highly stochastic in nature because of the statistical number fluctuations that arise on account of the small system size at the cellular scale. We study the nucleated reversible self-assembly of monomeric building blocks into polymer-like aggregates using the method of kinetic Monte Carlo. Kinetic Monte Carlo, being inherently stochastic, allows us to study the impact of fluctuations on the polymerisation reactions. One of the most important characteristic features in this kind of problem is the existence of a lag phase before self-assembly takes off, which is what we focus attention on. We study the associated lag time as a function of the system size and kinetic pathway. We find that the leading order stochastic contribution to the lag time before polymerisation commences is inversely proportional to the system volume for large-enough system size for all nine reaction pathways tested. Finite-size corrections to this do depend on the kinetic pathway.

preprint2015arXiv

Percolation in suspensions of hard nanoparticles: From spheres to needles

We investigate geometric percolation and scaling relations in suspensions of nanorods, covering the entire range of aspect ratios from spheres to extremely slender needles. A new version of connectedness percolation theory is introduced and tested against specialized Monte Carlo simulations. The theory accurately predicts percolation thresholds for aspect ratios as low as 10. The percolation threshold for rod-like particles of aspect ratios below 1000 deviates significantly from the inverse aspect ratio scaling prediction, thought to be valid in the limit of infinitely slender rods and often used as a rule of thumb for nano-fibers in composite materials. Hence, most fibers that are currently used as fillers in composite materials cannot be regarded as practically infinitely slender for the purposes of percolation theory. Comparing percolation thresholds of hard rods and new benchmark results for ideal rods, we find that (i) for large aspect ratios, they differ by a factor that is inversely proportional to the connectivity distance between the hard cores, and (ii) they approach the slender rod limit differently.

preprint2015arXiv

Percolation in suspensions of polydisperse hard rods : quasi-universality and finite-size effects

We present a study of connectivity percolation in suspensions of hard spherocylinders by means of Monte Carlo simulation and connectedness percolation theory. We focus attention on polydispersity in the length, the diameter and the connectedness criterion, and invoke bimodal, Gaussian and Weibull distributions for these. The main finding from our simulations is that the percolation threshold shows quasi universal behaviour, i.e., to a good approximation it depends only on certain cumulants of the full size and connectivity distribution. Our connectedness percolation theory hinges on a Lee-Parsons type of closure recently put forward that improves upon the often-used second virial approximation [ArXiv e-prints, May 2015, 1505.07660]. The theory predicts exact universality. Theory and simulation agree quantitatively for aspect ratios in excess of 20, if we include the connectivity range in our definition of the aspect ratio of the particles. We further discuss the mechanism of cluster growth that, remarkably, differs between systems that are polydisperse in length and in width, and exhibits non-universal aspects.

preprint2013arXiv

Fractional hopping-type motion in columnar mesophase of semi-flexible rod-like particles

We report on single-particle dynamics of strongly interacting filamentous fd virus particles in the liquid-crystalline columnar state in aqueous solution. From fluorescence microscopy we find that rare, discrete events take place, in which individual particles engage in sudden, jump-like motion along the main rod axis. The jump length distribution is bimodal and centered at half and full particle lengths. Our Brownian dynamics simulations of hard semi-flexible particles mimic our experiments, and indicate that full-length jumps must be due to collective dynamics in which particles move in string-like fashion in and between neighboring columns, while half jumps arise as a result of particles moving into defects. We find that the finite domain structure of the columnar phase strongly influences the observed dynamics.

preprint2013arXiv

Nematic droplets in aqueous dispersions of carbon nanotubes

Aqueous dispersions of exfoliated, bile-salt stabilized single-wall carbon nanotubes exhibit a first order transition to a nematic liquid-crystalline phase. The nematic phase presents itself in the form of micron-sized nematic droplets also known as tactoids, freely floating in the isotropic host dispersion. The nematic droplets are spindle shaped and have an aspect ratio of about four, irrespective of their size. We attribute this to a director field that is uniform rather than bipolar, which is confirmed by polarization microscopy. It follows that the ratio of the anchoring strength and the surface tension must be about four, which is quite larger than predicted theoretically but in line with earlier observations of bipolar tactoids. From the scatter in the data we deduce that the surface tension of the coexisting isotropic and nematic phases must be extremely low, that is, of the order of nN/m.

preprint2013arXiv

Orientational order of carbon nanotube guests in a nematic host suspension of colloidal viral rods

In order to investigate the coupling between the degrees of alignment of elongated particles in binary nematic dispersions, surfactant stabilized single-wall carbon nanotubes (CNTs) have been added to nematic suspensions of colloidal rodlike viruses in aqueous solution.We have independently measured the orientational order parameter of both components of the guest-host system by means of polarized Raman spectroscopy and by optical birefringence, respectively. Our system allows us therefore to probe the regime where the guest particles (CNTs) are shorter and thinner than the fd virus host particles. We show that the degree of order of the CNTs is systematically smaller than that of the fd virus particles for the whole nematic range. These measurements are in good agreement with predictions of an Onsager-type second-viral theory, which explicitly includes the flexibility of the virus particles, and the polydispersity of the CNTs.

preprint2013arXiv

Quasiuniversal connectedness percolation of polydisperse rod systems

The connectedness percolation threshold (eta_c) and critical coordination number (Z_c) of systems of penetrable spherocylinders characterized by a length polydispersity are studied by way of Monte Carlo simulations for several aspect ratio distributions. We find that (i) η_c is a nearly universal function of the weight-averaged aspect ratio, with an approximate inverse dependence that extends to aspect ratios that are well below the slender rod limit and (ii) that percolation of impenetrable spherocylinders displays a similar quasiuniversal behavior. For systems with a sufficiently high degree of polydispersity, we find that Z_c can become smaller than unity, in analogy with observations reported for generalized and complex networks.

preprint2012arXiv

Impact of Charge Variation on the Encapsulation of Nanoparticles by Virus Coat Proteins

Electrostatic interaction is the driving force for the encapsulation by virus coat proteins of nanoparticles such as quantum dots, gold particles and magnetic beads for, e.g., imaging and therapeutic purposes. In recent experimental work, Daniel et al. [ACS Nano 4 (2010), 3853-3860] found the encapsulation efficiency to sensitively depend on the interplay between the surface charge density of negatively charged gold nanoparticles and the number of positive charges on the RNA binding domains of the proteins. Surprisingly, these experiments reveal that despite the highly cooperative nature of the co-assembly at low pH, the efficiency of encapsulation is a gradual function of their surface charge density. We present a simple all-or-nothing mass action law combined with an electrostatic interaction model to explain the experiments. We find quantitative agreement with experimental observations, supporting the existence of a natural statistical charge distribution between nanoparticles.

preprint2012arXiv

Random, blocky and alternating ordering in supramolecular polymers of chemically bidisperse monomers

As a first step to understanding the role of molecular or chemical polydispersity in self-assembly, we put forward a coarse-grained model that describes the spontaneous formation of quasi-linear polymers in solutions containing two self-assembling species. Our theoretical framework is based on a two-component self-assembled Ising model in which the bidispersity is parameterized in terms of the strengths of the binding free energies that depend on the monomer species involved in the pairing interaction. Depending upon the relative values of the binding free energies involved, different morphologies of assemblies that include both components are formed, exhibiting paramagnetic-, ferromagnetic- or anti ferromagnetic-like order, i.e., random, blocky or alternating ordering of the two components in the assemblies. Analyzing the model for the case of ferromagnetic ordering, which is of most practical interest, we find that the transition from conditions of minimal assembly to those characterized by strong polymerization can be described by a critical concentration that depends on the concentration ratio of the two species. Interestingly, the distribution of monomers in the assemblies is different from that in the original distribution, i.e., the ratio of the concentrations of the two components put into the system. The monomers with a smaller binding free energy are more abundant in short assemblies and monomers with a larger binding affinity are more abundant in longer assemblies. Under certain conditions the two components congregate into separate supramolecular polymeric species and in that sense phase separate. We find strong deviations from the expected growth law for supramolecular polymers even for modest amounts of a second component, provided it is chemically sufficiently distinct from the main one.

preprint2012arXiv

Random, blocky and alternating ordering in supramolecular polymers of chemically bidisperse monomers

As a first step to understanding the role of molecular or chemical polydispersity in self-assembly, we put forward a coarse-grained model that describes the spontaneous formation of quasi-linear polymers in solutions containing two self-assembling species. Our theoretical framework is based on a two-component self-assembled Ising model in which the bidispersity is parameterized in terms of the strengths of the binding free energies that depend on the monomer species involved in the pairing interaction. Depending upon the relative values of the binding free energies involved, different morphologies of assemblies that include both components are formed, exhibiting paramagnetic-, ferromagnetic- or anti ferromagnetic-like order,i.e., random, blocky or alternating ordering of the two components in the assemblies. Analyzing the model for the case of ferromagnetic ordering, which is of most practical interest, we find that the transition from conditions of minimal assembly to those characterized by strong polymerization can be described by a critical concentration that depends on the concentration ratio of the two species. Interestingly, the distribution of monomers in the assemblies is different from that in the original distribution, i.e., the ratio of the concentrations of the two components put into the system. The monomers with a smaller binding free energy are more abundant in short assemblies and monomers with a larger binding affinity are more abundant in longer assemblies. Under certain conditions the two components congregate into separate supramolecular polymeric species and in that sense phase separate. We find strong deviations from the expected growth law for supramolecular polymers even for modest amounts of a second component, provided it is chemically sufficiently distinct from the main one.

preprint2008arXiv

Self diffusion of particles in complex fluids: temporary cages and permanent barriers

We study the self diffusion of individual particles in dense (non-)uniform complex fluids within dynamic density functional theory and explicitly account for their coupling to the temporally fluctuating background particles. Applying the formalism to rod-like particles in uniaxial nematic and smectic liquid crystals, we find correlated diffusion in different directions: The temporary cage formed by the neighboring particles competes with permanent barriers in periodic inhomogeneous systems such as the lamellar smectic state and delays self diffusion of particles even in uniform systems. We compare our theory with recent experimental data on the self diffusion of fluorescently labelled filamentous virus particles in aqueous dispersions in the smectic phase and find qualitative agreement. This demonstrates the importance of explicitly dealing with the time-dependent self-consistent molecular field that every particle experiences.

preprint2004arXiv

Comment on "Large Difference in the Elastic Properties of fcc and hcp Hard-Sphere Crystals"

As is well known, hard-sphere crystals of the fcc and hcp type differ very little in their thermodynamic properties. Nonetheless, recent computer simulations by Pronk and Frenkel indicate that the elastic response to mechanical deformation of the two types of crystal should be quite different. By invoking a geometrical argument put forward by R. Martin some time ago, we suggest that this is largely due to the different symmetries of the fcc and hcp crystal structures. Indeed, we find that elastic constants obtained by means of computer simulations for the fcc hard-sphere crystal can be mapped onto the equivalent ones of the hcp crystal to very high accuracy. The same procedure applied to density functional theoretical predictions for the elastic properties of the fcc hard-sphere crystal also produces remarkably accurate predictions for those of the hcp hard-sphere crystal.