Custom flow in overdamped Brownian Dynamics

Daniel de las Heras,Johannes Renner,Matthias Schmidt
DOI: https://doi.org/10.1103/PhysRevE.99.023306
2018-12-06
Abstract:When an external field drives a colloidal system out of equilibrium, the ensuing colloidal response can be very complex and obtaining a detailed physical understanding often requires case-by-case considerations. In order to facilitate systematic analysis, here we present a general iterative scheme for the determination of the unique external force field that yields a prescribed inhomogeneous stationary or time-dependent flow in an overdamped Brownian many-body system. The computer simulation method is based on the exact one-body force balance equation and allows to specifically tailor both gradient and rotational velocity contributions, as well as to freely control the one-body density distribution. Hence compressibility of the flow field can be fully adjusted. The practical convergence to a unique external force field demonstrates the existence of a functional map from both velocity and density to external force field, as predicted by the power functional variational framework. In equilibrium, the method allows to find the conservative force field that generates a prescribed target density profile, and hence implements the Mermin-Evans classical density functional map from density distribution to external potential. The conceptual tools developed here enable one to gain detailed physical insight into complex flow behaviour, as we demonstrate in prototypical situations.
Soft Condensed Matter,Statistical Mechanics,Fluid Dynamics
What problem does this paper attempt to address?
The problem that this paper attempts to solve is how to determine the unique external force field that can generate a preset non - uniform steady - state or time - dependent flow in an overdamped Brownian dynamics system. Specifically, the researchers developed an iterative method to determine the specific external force field that can keep a given density distribution and velocity field stable in a non - equilibrium state. This method is based on the exact one - body force - balance equation and allows researchers to freely control the density distribution as well as the gradient and rotational components in the flow field, thereby achieving a full adjustment of the compressibility of the flow field. ### Background and Problem of the Paper When an external field drives a colloidal system away from equilibrium, the response of the colloids can be very complex, and obtaining a detailed physical understanding usually requires case - by - case analysis. To facilitate systematic analysis, the author proposes a general iterative scheme for determining the unique external force field that can generate a preset non - uniform steady - state or time - dependent flow. This method is based on computer simulations, utilizes the exact one - body force - balance equation, and can specifically customize the gradient and rotational components in the flow field and freely control the one - body density distribution. ### Main Contributions 1. **Iterative Method**: The author proposes an iterative method that can determine the unique external force field that generates a preset flow. This method is applicable not only to the steady - state case but also to the time - dependent non - equilibrium state. 2. **Functional Mapping**: Through this method, the author verifies the existence and uniqueness of the functional mapping from velocity and density to the external force field, which is consistent with the prediction of the power - functional variational framework. 3. **Application Examples**: The author shows how to customize the flow in a many - body system in several model cases, for example, generating a specific density distribution and velocity field in a two - dimensional system. ### Mathematical Description - **Density Distribution and Current Distribution**: \[ \rho(\mathbf{r}, t)=\left\langle\sum_{i} \delta(\mathbf{r}-\mathbf{r}_{i})\right\rangle \] \[ \mathbf{J}(\mathbf{r}, t)=\left\langle\sum_{i} \delta(\mathbf{r}-\mathbf{r}_{i}) \mathbf{v}_{i}\right\rangle \] - **Continuity Equation**: \[ \frac{\partial \rho(\mathbf{r}, t)}{\partial t}=-\nabla \cdot \mathbf{J}(\mathbf{r}, t) \] - **Force - Balance Equation**: \[ \gamma \mathbf{J}=-k_{B} T \nabla \rho+\mathbf{F}_{\text {int}}+\rho \mathbf{f}_{\text {ext}} \] where $\mathbf{F}_{\text {int}}$ is the internal force density field and $\mathbf{f}_{\text {ext}}$ is the external force field. - **Iterative Steps**: \[ \mathbf{f}_{\text {ext}}^{(k)}=k_{B} T \nabla \ln \rho-\mathbf{f}_{\text {int}}^{(k - 1)}+\gamma \mathbf{v} \] Initial Conditions: \[ \mathbf{f}_{\text {ext}}^{(0)}=k_{B} T \nabla \ln \rho+\gamma \mathbf{v} \] ### Results and Discussion The author verifies the effectiveness of this method through specific numerical simulations. For example, in a one - dimensional system, they show how to generate multiple steady - states with the same density distribution but different constant currents. In addition, they also generate a specific density distribution and velocity field in a two - dimensional system, demonstrating the wide applicability of this method. ### Conclusion This research provides a powerful tool for systematically studying non - equilibrium.