Computational modeling of passive transport of functionalized nanoparticles.

Functionalized nanoparticles (NPs) are complex objects present in a variety of systems ranging from synthetic grafted nanoparticles to viruses. The morphology and number of the decorating groups can vary widely between systems. Thus, the modeling of functionalized NPs typically considers simplified spherical objects as a first-order approximation.

A simple catch: Fluctuations enable hydrodynamic trapping of microrollers by obstacles.

It is known that obstacles can hydrodynamically trap bacteria and synthetic microswimmers in orbits, where the trapping time heavily depends on the swimmer flow field and noise is needed to escape the trap. Here, we use experiments and simulations to investigate the trapping of microrollers by obstacles. Microrollers are rotating particles close to a bottom surface, which have a prescribed propulsion direction imposed by an external rotating magnetic field.

Hydrodynamics of spike proteins dictate a transport-affinity competition for SARS-CoV-2 and other enveloped viruses.

Many viruses, such as SARS-CoV-2 or Influenza, possess envelopes decorated with surface proteins (a.k.a.

Metallic microswimmers driven up the wall by gravity.

Experiments on autophoretic bimetallic nanorods propelling within a fuel of hydrogen peroxide show that tail-heavy swimmers preferentially orient upwards and ascend along inclined planes. We show that such gravitaxis is strongly facilitated by interactions with solid boundaries, allowing even ultraheavy microswimmers to climb nearly vertical surfaces. Theory and simulations show that the buoyancy or gravitational torque that tends to align the rods is reinforced by a fore-aft drag asymmetry induced by hydrodynamic interactions with the wall.

Relating Rheotaxis and Hydrodynamic Actuation using Asymmetric Gold-Platinum Phoretic Rods.

We explore the behavior of micron-scale autophoretic Janus (Au/Pt) rods, having various Au/Pt length ratios, swimming near a wall in an imposed background flow. We find that their ability to robustly orient and move upstream, i.e.

Large scale Brownian dynamics of confined suspensions of rigid particles.

We introduce methods for large-scale Brownian Dynamics (BD) simulation of many rigid particles of arbitrary shape suspended in a fluctuating fluid. Our method adds Brownian motion to the rigid multiblob method [F. Balboa Usabiaga et al.

Brownian dynamics of confined suspensions of active microrollers.

We develop efficient numerical methods for performing many-body Brownian dynamics simulations of a recently observed fingering instability in an active suspension of colloidal rollers sedimented above a wall [M. Driscoll, B. Delmotte, M.

Rapid sampling of stochastic displacements in Brownian dynamics simulations.

We present a new method for sampling stochastic displacements in Brownian Dynamics (BD) simulations of colloidal scale particles. The method relies on a new formulation for Ewald summation of the Rotne-Prager-Yamakawa (RPY) tensor, which guarantees that the real-space and wave-space contributions to the tensor are independently symmetric and positive-definite for all possible particle configurations. Brownian displacements are drawn from a superposition of two independent samples: a wave-space (far-field or long-ranged) contribution, computed using techniques from fluctuating hydrodynamics and non-uniform fast Fourier transforms; and a real-space (near-field or short-ranged) correction, computed using a Krylov subspace method.

Brownian dynamics of confined rigid bodies.

We introduce numerical methods for simulating the diffusive motion of rigid bodies of arbitrary shape immersed in a viscous fluid. We parameterize the orientation of the bodies using normalized quaternions, which are numerically robust, space efficient, and easy to accumulate. We construct a system of overdamped Langevin equations in the quaternion representation that accounts for hydrodynamic effects, preserves the unit-norm constraint on the quaternion, and is time reversible with respect to the Gibbs-Boltzmann distribution at equilibrium.

A multiblob approach to colloidal hydrodynamics with inherent lubrication.

This work presents an intermediate resolution model of the hydrodynamics of colloidal particles based on a mixed Eulerian-Lagrangian formulation. The particle is constructed with a small set of overlapping Peskin's Immersed Boundary kernels (blobs) which are held together by springs to build up a particle impenetrable core. Here, we used 12 blobs placed in the vertexes of an icosahedron with an extra one in its center.

Brownian dynamics without Green's functions.

We develop a Fluctuating Immersed Boundary (FIB) method for performing Brownian dynamics simulations of confined particle suspensions. Unlike traditional methods which employ analytical Green's functions for Stokes flow in the confined geometry, the FIB method uses a fluctuating finite-volume Stokes solver to generate the action of the response functions "on the fly." Importantly, we demonstrate that both the deterministic terms necessary to capture the hydrodynamic interactions among the suspended particles, as well as the stochastic terms necessary to generate the hydrodynamically correlated Brownian motion, can be generated by solving the steady Stokes equations numerically only once per time step.

The Stokes-Einstein relation at moderate Schmidt number.

The Stokes-Einstein relation for the self-diffusion coefficient of a spherical particle suspended in an incompressible fluid is an asymptotic result in the limit of large Schmidt number, that is, when momentum diffuses much faster than the particle. When the Schmidt number is moderate, which happens in most particle methods for hydrodynamics, deviations from the Stokes-Einstein prediction are expected. We study these corrections computationally using a recently developed minimally resolved method for coupling particles to an incompressible fluctuating fluid in both two and three dimensions.

Applications of computational geometry to the molecular simulation of interfaces.

The identification of the interfacial molecules in fluid-fluid equilibrium is a long-standing problem in the area of simulation. We here propose an alternative point of view, making use of concepts taken from the field of computational geometry, where the definition of the "shape" of a set of points is a well-known problem. In particular, we employ the alpha -shape construction which, applied to the positions of the molecules, selects a shape and identifies its boundary points, which we will take to define our interfacial molecules.