Source author record

Margarida M. Telo da Gama

Margarida M. Telo da Gama 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

18works
5topics
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

18 published item(s)

preprint2022arXiv

The role of disorder in the motion of chiral swimmers in the presence of obstacles

The presence of obstacles is intuitively expected to hinder the diffusive transport of micro-swimmers. However, for chiral micro-swimmers, a low density of obstacles near a surface can enhance their diffusive behavior, due to the rectification of the chiral motion by the obstacles. Here, we study numerically the role that disorder plays in determining the transport dynamics of chiral micro-swimmers on surfaces with obstacles. We consider different densities of regularly spaced obstacles and distinct types of disorder: noise in the dynamics of the micro-swimmer, quenched noise in the positions of the obstacles as well as obstacle size polydispersity. We show that, depending on the type and strength of the disorder, the presence of obstacles can either enhance or hinder transport, and discuss implications for the control of active transport in disordered media.

preprint2021arXiv

Dynamics of flowing 2D skyrmions

We investigate, numerically, the effects of externally imposed material flows on the structure and temporal evolution of liquid crystal skyrmions. The dynamics of a 2D system of skyrmions is modeled using the Ericksen-Leslie theory, which is based on two coupled equations, one for material flow and the other for the director field. As the time scales of the velocity and director fields differ by several orders of magnitude for realistic values of the system parameters, we have simplified the calculations by assuming that the velocity relaxes instantaneously when compared to the relaxation of the director field. Thus, we have used a finite-differences method known as artificial compressibility with adaptive time step to solve the velocity field and a fourth-order Runge-Kutta method for the director field. We characterized the skyrmion shape or configuration as a function of the time and the average velocity of the flow field. We found that for velocities above a certain threshold, the skyrmions stretch in the direction perpendicular to the flow, by contrast to the regime of weak flows where the skyrmions stretch along the streamlines of the flow field. These two regimes are separated by an abrupt (first-order) dynamical transition, which is robust with respect to e.g., the liquid crystal elastic anisotropy. Additionally, we have found how the presence of a second skyrmion affects the evolution of the shape of the skyrmions, by comparing the evolution of pairs of skyrmions to the evolution of a single-skyrmion.

preprint2021arXiv

Effect of particle deformability on shear thinning in a 3D channel

Soft particles suspended in fluids flowing through microchannels are often encountered in biological flows such as cells in vessels. They can deform under flow or when subject to mechanical stresses as they interact with themselves. Deformability adds distinct characteristics to the nature of the flow of these particles. We simulate deformable particles suspended in a fluid at a high volume fraction flowing through a 3D wetting channel. We find a discontinuous shear thinning transition, which depends on the particle deformability. The capillary number is the main parameter that controls the transition. We also notice significant differences in the velocity profiles when compared to their 2D counterpart. To obtain these results, we improved and extended to 3D a multicomponent Lattice Boltzmann method which prevents the coalescence between the droplets.

preprint2019arXiv

Active nematic-isotropic interfaces in channels

We use numerical simulations to investigate the hydrodynamic behavior of the interface between nematic (N) and isotropic (I) phases of a confined active liquid crystal. At low activities, a stable interface with constant shape and velocity is observed separating the two phases. For nematics in homeotropic channels, the velocity of the interface at the NI transition increases from zero (i) linearly with the activity for contractile systems and (ii) quadratically for extensile ones. Interestingly, the nematic phase expands for contractile systems while it contracts for extensile ones, as a result of the active forces at the interface. Since both activity and temperature affect the stability of the nematic, for active nematics in the stable regime the temperature can be tuned to observe static interfaces, providing an operational definition for the coexistence of active nematic and isotropic phases. At higher activities, beyond the stable regime, an interfacial instability is observed for extensile nematics. In this regime defects are nucleated at the interface and move away from it. The dynamics of these defects is regular and persists asymptotically for a finite range of activities. We used an improved hybrid model of finite differences and lattice Boltzmann method with multi-relaxation-time collision operator, the accuracy of which allowed us to characterize the dynamics of the distinct interfacial regimes.

preprint2019arXiv

Lattice-Boltzmann simulation of free nematic-isotropic interfaces

We use a hybrid method of lattice Boltzmann and finite differences to simulate flat and curved interfaces between the nematic and isotropic phases of a liquid crystal described by the Landau-de Gennes theory. For the flat interface, we measure the interfacial velocity at different temperatures around the coexistence. We show that the interface is completely static at the coexistence temperature and that the profile width is in line with the theoretical predictions. The interface is stable in a range of temperatures around coexistence and disappears when one of the two phases becomes mechanically unstable. We stabilize circular nematic domains by a shift in temperature, related to the Laplace pressure, and estimate the spurious velocities of these lattice Boltzmann simulations.

preprint2019arXiv

Propagation of active nematic-isotropic interfaces on substrates

Motivated by results for the propagation of active-passive interfaces of bacterial Serratia marcescens swarms [Nat. Comm., 9, 5373 (2018)] we use a hydrodynamic multiphase model to investigate the propagation of interfaces of active nematics on substrates. We characterize the active nematic phase of the model and discuss its description of the statistical dynamics of the swarms. We show that the velocicity correlation functions and the energy spectrum of the active turbulent phase of the model scale with the active length, for a range of activities. In addition, the energy spectrum exhibits two power-law regimes with exponents close to those reported for other models and bacterial swarms. Although the exponent of the rising branch of the spectrum (small wavevector) appears to be independent of the activity, the exponent of the decay changes with activity systematically, albeit slowly. We characterize also the propagation of circular and flat active-passive interfaces. We find that the closing time of the circular passive domain decays quadratically with the activity and that the structure factor of the flat interface is similar to that reported for swarms, with an activity dependent exponent. Finally, the effect of the substrate friction was investigated. We found an activity dependent threshold, above which the turbulent active nematic forms isolated islands that shrink until the system becomes isotropic and below which the active nematic expands, with a well defined propagating interface. The interface may be stopped by a fricion gradient.

preprint2016arXiv

Nematic liquid crystals on sinusoidal channels: the zigzag instability

Substrates which are chemically or topographically patterned induce a variety of liquid crystal textures. The response of the liquid crystal to competing surface orientations, typical of patterned substrates, is determined by the anisotropy of the elastic constants and the interplay of the relevant lengths scales, such as the correlation length and the surface geometrical parameters. Transitions between different textures, usually with different symmetries, may occur under a wide range of conditions. We use the Landau-de Gennes free energy to investigate the texture of nematics in sinusoidal channels with parallel anchoring bounded by nematic-air interfaces that favour perpendicular (hometropic) anchoring. In micron size channels 5CB was observed to exhibit a non-trivial texture characterized by a disclination line, within the channel, which is broken into a zigzag pattern. Our calculations reveal that when the elastic anisotropy of the nematic does not favour twist distortions the defect is a straight disclination line that runs along the channel, which breaks into a zigzag pattern with a characteristic period, when the twist elastic constant becomes sufficiently small when compared to the splay and bend constants. The transition occurs through a twist instability that drives the defect line to rotate from its original position. The interplay between the energetically favourable twist distortions that induce the defect rotation and the liquid crystal anchoring at the surfaces leads to the zigzag pattern. We investigate in detail the dependence of the periodicity of the zigzag pattern on the geometrical parameters of the sinusoidal channels, which in line with the experimental results is found to be non-linear.

preprint2016arXiv

Pattern-induced anchoring transitions in nematic liquid crystals

In this paper we revisit the problem of a nematic liquid crystal in contact with patterned substrates. The substrate is modelled as a periodic array of parallel infinite grooves of well-defined cross section sculpted on a chemically homogeneous substrate which favors local homeotropic anchoring of the nematic. We consider three cases: a sawtooth, a crenellated and a sinusoidal substrate. We analyse this problem within the modified Frank-Oseen formalism. We argue that, for substrate periodicities much larger than the extrapolation length, the existence of different nematic textures with distinct far-field orientations, as well as the anchoring transitions between them, are associated with the presence of topological defects either on or close to the substrate. For the sawtooth and sinusoidal case, we observe a homeotropic to planar anchoring transition as the substrate roughness is increased. On the other hand, a homeotropic to oblique anchoring transition is observed for crenellated substrates. In this case, the anchoring phase diagram shows a complex dependence on the substrate roughness and substrate anchoring strength.

preprint2016arXiv

Percolation in binary and ternary mixtures of patchy colloids

We investigate percolation in binary and ternary mixtures of patchy colloidal particles theoretically and using Monte Carlo simulations. Each particle has three identical patches, with distinct species having different types of patch. Theoretically we assume tree-like clusters and calculate the bonding probabilities using Wertheim's first-order perturbation theory for association. For ternary mixtures we find up to eight fundamentally different percolated states. The states differ in terms of the species and pairs of species that have percolated. The strongest gel is a trigel or tricontinuous gel, in which each of the three species has percolated. The weakest gel is a mixed gel in which all of the particles have percolated, but none of the species percolates by itself. The competition between entropy of mixing and internal energy of bonding determines the stability of each state. Theoretical and simulation results are in very good agreement. The only significant difference is the temperature at the percolation threshold, which is overestimated by the theory due to the absence of closed loops in the theoretical description.

preprint2016arXiv

Temperature (de)activated patchy colloidal particles

We present a new model of patchy particles in which the interaction sites can be activated or deactivated by varying the temperature of the system. We study the thermodynamics of the system by means of Wertheim's first order perturbation theory, and use Flory-Stockmayer theory of polymerization to analyse the percolation threshold. We find a very rich phase behaviour including lower critical points and reentrant percolation.

preprint2016arXiv

The effect of curvature on cholesteric doughnuts

The confinement of liquid crystals inside curved geometries leads to exotic structures, with applications ranging from bio-sensors to optical switches and privacy windows. Here we study how curvature affects the alignment of a cholesteric liquid crystal. We model the system on the mesoscale using the Landau-de Gennes model. Our study was performed in three stages, analysing different curved geometries from cylindrical walls and pores, to toroidal domains, in order to isolate the curvature effects. Our results show that the stresses introduced by the curvature influence the orientation of the liquid crystal molecules, and cause distortions in the natural periodicity of the cholesteric that depend on the radius of curvature and on the commensurability between the pitch and the dimensions of the system. In particular, the cholesteric layers of toroidal droplets exhibit a symmetry breaking not seen in cylindrical pores and that is driven by the additional curvature.

preprint2015arXiv

Wetting of cholesteric liquid crystals

We investigate theoretically the wetting properties of cholesteric liquid crystals at a planar substrate. If the properties of substrate and of the interface are such that the cholesteric layers are not distorted the wetting properties are similar to those of a nematic liquid crystal. If, on the other hand, the anchoring conditions force the distortion of the liquid crystal layers the wetting properties are altered, the free cholesteric-isotropic interface is non-planar and there is a layer of topological defects close to the substrate. These deformations can either promote or hinder the wetting of the substrate by a cholesteric, depending on the properties of the cholesteric liquid crystal.

preprint2013arXiv

Computing the phase diagram of binary mixtures: a patchy particle case study

We investigate the phase behaviour of 2D mixtures of bi-functional and three-functional patchy particles and 3D mixtures of bi-functional and tetra-functional patchy particles by means of Monte Carlo simulations and Wertheim theory. We start by computing the critical points of the pure systems and then we investigate how the critical parameters change upon lowering the temperature. We extend the Successive Umbrella Sampling method to mixtures to make it possible to extract information about the phase behaviour of the system at a fixed temperature for the whole range of densities and compositions of interest.

preprint2013arXiv

Interfacial motion in flexo- and order-electric switching between nematic filled states

We consider a nematic liquid crystal, in coexistence with its isotropic phase, in contact with a substrate patterned with rectangular grooves. In such a system, the nematic phase may fill the grooves without the occurrence of complete wetting. There may exist multiple (meta)stable filled states, each characterised by the type of distortion (bend or splay) in each corner of the groove and by the shape of the nematic-isotropic interface, and additionally the plateaux that separate the grooves may be either dry or wet with a thin layer of nematic. Using numerical simulations, we analyse the dynamical response of the system to an externally- applied electric field, with the aim of identifying switching transitions between these filled states. We find that order-electric coupling between the fluid and the field provides a means of switching between states where the plateaux between grooves are dry and states where they are wet by a nematic layer, without affecting the configuration of the nematic within the groove. We find that flexoelectric coupling may change the nematic texture in the groove, provided that the flexoelectric coupling differentiates between the types of distortion at the corners of the substrate. We identify intermediate stages of the transitions, and the role played by the motion of the nematic-isotropic interface. We determine quantitatively the field magnitudes and orientations required to effect each type of transition.

preprint2012arXiv

Properties of patchy colloidal particles close to a surface: a Monte Carlo and density functional study

We investigate the behavior of a patchy particle model close to a hard-wall via Monte Carlo simulation and density functional theory (DFT). Two DFT approaches, based on the homogeneous and inhomogeneous versions of Wertheim's first order perturbation theory for the association free energy are used. We evaluate, by simulation and theory, the equilibrium bulk phase diagram of the fluid and analyze the surface properties for two isochores, one of which is close to the liquid side of the gas-liquid coexistence curve. We find that the density profile near the wall crosses over from a typical high-temperature adsorption profile to a low-temperature desorption one, for the isochore close to coexistence. We relate this behavior to the properties of the bulk network liquid and find that the theoretical descriptions are reasonably accurate in this regime. At very low temperatures, however, an almost fully bonded network is formed, and the simulations reveal a second adsorption regime which is not captured by DFT. We trace this failure to the neglect of orientational correlations of the particles, which are found to exhibit surface induced orientational order in this regime.

preprint2011arXiv

Bicontinuous and mixed gels in binary mixtures of patchy colloidal particles

We investigate the thermodynamics and percolation regimes of model binary mixtures of patchy colloidal particles. The particles of each species have three sites of two types, one of which promotes bonding of particles of the same species while the other promotes bonding of different species. We find up to four percolated structures at low temperatures and densities: two gels where only one species percolates, a mixed gel where particles of both species percolate but neither species percolates separately, and a bicontinuous gel where particles of both species percolate separately forming two interconnected networks. The competition between the entropy and the energy of bonding drives the stability of the different percolating structures. Appropriate mixtures exhibit one or more connectivity transitions between the mixed and bicontinuous gels, as the temperature and/or the composition changes.

preprint2011arXiv

Phase diagrams of binary mixtures of patchy colloids with distinct numbers and types of patches: The empty fluid regime

We investigate the effect of distinct bonding energies on the onset of criticality of low functionality fluid mixtures. We focus on mixtures of particles with two and three patches as this includes the mixture where 'empty' fluids were originally reported. In addition to the number of patches, the species differ in the type of patches or bonding sites. For simplicity, we consider that the patches on each species are identical: one species has 3 patches of type $A$ and the other 2 patches of type $B$. We have found a rich phase behaviour with closed miscibility gaps, liquid-liquid demixing and negative azeotropes. Liquid-liquid demixing was found to preempt the 'empty' fluid regime, of these mixtures, when the $AB$ bonds are weaker than the $AA$ or $BB$ bonds. By contrast, mixtures in this class exhibit 'empty' fluid behaviour when the $AB$ bonds are stronger than at least one of the other two. Mixtures with bonding energies $ε_{BB}=ε_{AB}$ and $ε_{AA}<ε_{BB}$, were found to exhibit an unusual negative azeotrope.

preprint2011arXiv

Phase diagrams of binary mixtures of patchy colloids with distinct numbers of patches: The network fluid regime

We calculate the network fluid regime and phase diagrams of binary mixtures of patchy colloids, using Wertheim's first order perturbation theory and a generalization of Flory-Stockmayer's theory of polymerization. The colloids are modelled as hard spheres with the same diameter and surface patches of the same type, $A$. The only difference between species is the number of their patches -or functionality-, $f_A^{(1)}$ and $f_A^{(2)}$ (with $f_A^{(2)}>f_A^{(1)}$). We have found that the difference in functionality is the key factor controlling the behaviour of the mixture in the network (percolated) fluid regime. In particular, when $f_A^{(2)}\ge2f_A^{(1)}$ the entropy of bonding drives the phase separation of two network fluids which is absent in other mixtures. This changes drastically the critical properties of the system and drives a change in the topology of the phase diagram (from type I to type V) when $f_A^{(1)}>2$. The difference in functionality also determines the miscibility at high (osmotic) pressures. If $f_A^{(2)}-f_A^{(1)}=1$ the mixture is completely miscible at high pressures, while closed miscibility gaps at pressures above the highest critical pressure of the pure fluids are present if $f_A^{(2)}-f_A^{(1)}>1$. We argue that this phase behaviour is driven by a competition between the entropy of mixing and the entropy of bonding, as the latter dominates in the network fluid regime.