# Effective temperature concept evaluated in an active colloid mixture

^{a}Applied Physics Graduate Program, Northwestern University, Evanston, IL 60208;^{b}Department of Materials Science and Engineering, University of Illinois, Urbana, IL 61801;^{c}Institute for Basic Science Center for Soft and Living Matter, Ulsan National Institute of Science and Technology, Ulsan 689-798, South Korea;^{d}Department of Materials Science and Engineering, Northwestern University, Evanston, IL 60208;^{e}Department of Engineering Sciences and Applied Mathematics, Northwestern University, Evanston, IL 60208;^{f}Department of Physics and Astronomy, Northwestern University, Evanston, IL 60208

See allHide authors and affiliations

Contributed by Steve Granick, June 11, 2017 (sent for review April 24, 2017; reviewed by Denis Bartolo and Angelo Cacciuto)

## Significance

In the long-standing debate on whether the insights of statistical thermodynamics could apply to nonequilibrium systems, a key unsolved problem is whether it is useful to postulate an emergent effective temperature, different from the thermal temperature. Here, we introduce a two-component system of driven Janus colloids such that collisions induced by external energy sources play the role of temperature. We find that this nonequilibrium system quantitatively behaves as if at equilibrium, with collisions caused by differential rhythmic motion between the two components acting as a strict analog to thermal motion. We demonstrate that the effective temperature can serve as a common control parameter for both kinetics and phase behavior.

## Abstract

Thermal energy agitates all matter, and its competition with ordering tendencies is a fundamental organizing principle in the physical world; this observation suggests that an effective temperature might emerge when external energy input enhances agitation. However, despite the repeated proposal of this concept based on kinetics for various nonequilibrium systems, the value of an effective temperature as a thermodynamic control parameter has been unclear. Here, we introduce a two-component system of driven Janus colloids, such that collisions induced by external energy sources agitate the system, and we demonstrate quantitative agreement with hallmarks of statistical thermodynamics for binary phase behavior: the archetypal phase diagram with equilibrium critical exponents, Gaussian displacement distributions, and even capillarity. The significance is to demonstrate a class of dynamical conditions under which thermodynamic analysis extends quantitatively to systems that are decidedly nonequilibrium except that the effective temperature differs from the physical temperature.

The tug of war between randomization due to temperature and ordering from enthalpy governs whether components mix or demix at equilibrium. The even larger realm of nonequilibrium systems similarly shows empirically that particles with different dynamic modes may segregate spontaneously. Examples include bacteria (1), granular particles under shaking (2), and active colloids (3). This concept has inspired various efforts (4, 5) to extend the profound insights of statistical thermodynamics to nonequilibrium systems. In particular, the concept of defining an effective temperature based upon the kinetics of a system has repeatedly been proposed for synthetic active materials (6, 7), biological motors (8⇓–10), granular materials (11⇓–13), and turbulent fluids (14). However, the predictive value of this effective temperature has been explored only to a limited extent (15, 16). Here, without wishing to minimize important differences between specific nonequilibrium phenomena, we test the concept that some of them may be indeed governed by an emergent effective temperature.

## Results

Experimentally, we apply external electric and magnetic fields (Fig. 1*A*) to Janus silica spheres with a diameter of 3 μm, one hemisphere coated with a thin Ni/Ti cap and an SiO_{2} protective layer to render their outermost chemical makeup homogeneous. The particles settle to the bottom of the planar sample cell to form a quasi–two-dimensional (quasi-2D) system, but, as they remain dispersed in water, we avoid complications of solid–solid friction (17) often encountered in granular systems. The ac electric field applied perpendicular to the monolayer causes the particles to orient their equators vertically and self-propel (“swim”) with their nonmetallic hemisphere facing forward (18), while a magnetic field, rotating at low frequency in the plane, forces swimming particles to trace circular orbits with a tunable radius *R* (Movie S1). Spontaneous symmetry breaking results. The Janus colloids dynamically differentiate into two types that orbit 180 degrees out of phase, with their nonmetallic sides facing in opposite directions. Particles that are in phase never collide, whereas those that orbit out of phase have opportunities to collide with each other when close by (Fig. 1*B*). Brownian dynamics computer simulations, matching the experimental conditions in their essential respects, extend the experiments to otherwise inaccessible regimes. In these simulations, particles carrying an electric dipole moment are confined within a 2D plane, and experience a swimming force directed from the metal toward the silica hemisphere. The effect of the rotating magnetic field is modeled as an external torque driving the circular motion of each particle. Without this rotation, oppositely swimming particles would separate into lanes to reduce collisions, a known phenomenon (19). However, at finite *R*, lanes grow only to a finite length (Fig. S1) as they continually break and reform to follow the orbits (Movie S2). Viscous drag dissipates the injected energy, and the collisions disturbing the intrinsic circular orbits cause particle trajectories to fluctuate.

Hallmarks of classical Brownian motion are observed when we study the long-time behavior of particle dynamics by inspecting their positions once per rotation cycle (“stroboscopic” view), at symmetric (50:50) composition. First, the mean squared displacements are proportional to the number of cycles (Fig. 1*C*). Second, displacements trace a Gaussian distribution over three orders of magnitude of probability density that were inspected (Fig. 1*D*). The Gaussian distribution was further confirmed by inspection of the kurtosis (Fig. S2). We emphasize that our experiments are consistent with the simulations wherever comparison is possible, whereas the simulations offer better and more exhaustive statistics. Increasing diffusion with *R* suggests that *R* plays the role of effective temperature, but it remains to explain why.

Strong claims demand strong evidence. According to the textbook Langevin approach (20) that describes thermal fluctuation dynamics in equilibrium, random forces must satisfy the following requirements: (*i*) their (time or ensemble) average must vanish; (*ii*) orthogonal components of the random forces on a particle must be independent; and (*iii*) each component of the random force must decorrelate instantaneously and have a variance proportional to temperature. This very approach, mathematically equivalent to the Fokker–Planck equation (20), is used widely in molecular dynamics simulations to control temperature. We tested these points in simulations. Condition *i* is satisfied by symmetry of the circular motion. In this overdamped system, random forces originate from interparticle collisions, but are immediately balanced by the viscous forces and manifest themselves as perturbations around the circular orbits. Satisfying conditions *ii* and *iii*, orthogonal components of the perturbations tangential and normal to the circular orbit are independent, and both decorrelate at short times, within one tenth of a cycle, as shown in Fig. 2.

The resulting Smoluchowski equation (21) explains the emergence of Brownian motion. Consistently, the diffusivity *D* measured on long time scales and its value computed from the autocorrelation function of the perturbations agree quantitatively (Fig. 2*C*). Furthermore, the Gaussian displacement distribution that we observe is also explained, as the Smoluchowski equation yields such a distribution with variance determined by the variance of the random forces, regardless of whether these forces are Gaussian random variables (20) (Fig. 2*A*). The rapid decorrelation of the random forces (condition *iii* above; Fig. 2*B*) guarantees Markovian dynamics in the stroboscopic view. In combination with the steady-state nature of the system, this Markovian dynamics implies that the fluctuation–dissipation theorem applies even out of equilibrium (9). According to the Einstein relation, *D* is proportional to absolute temperature. Because we find that, in our system, *D* is linear with the radius *R* (Fig. 2*C*), it is physically meaningful to represent and control the effective temperature by the radius. Indeed, as the radius increases, more interparticle collisions occur within each cycle, and thus larger randomness is introduced into the system, analogous to increasing the temperature.

Another hallmark of thermodynamics is phase separation, as illustrated in Fig. 3. Investigating this, we note that phase separation between the two dynamical species arises because collisions between these particles of different types lead to effective attraction between particles of the same type. Despite the absence of enthalpy of interaction, there is a tendency to aggregate, akin to the concept of osmotic pressure. To quantify this, we extract the potential of mean force (22, 23) *V*_{PMF} from the radial distribution functions *g*(*r*) of the two species (Fig. 3*B*). This procedure yields the pair potential in the intrinsic “thermal” energy unit, even though, here, that is not *k*_{B}*T*. The fact that attraction strengthens with decreasing *R* suggests that phase segregation might ensue. Indeed, we observe this (Fig. 3*A* and Movie S3) experimentally.

The development of equilibrium statistical thermodynamics was revolutionized when it succeeded to explain the critical point with the concomitant universal singular behavior (24), so now we likewise apply this stringent test to the nonequilibrium dynamical system introduced here. Symmetry dictates that this point must occur in a system with equal numbers of A and B particles. When *R* is decreased to a critical value *R*_{c} = 1.84*σ* at area fraction 0.37 in the simulations (between *σ* and 2*σ* in experiment, depending on particle density), the spatial distribution of the local order parameter *s* = *ϕ*_{A} − *ϕ*_{B} (*ϕ* is the composition fraction) shows self-similar patterns indicating a diverging correlation length, the hallmark of criticality (Fig. 4*A* and Movie S4). Consistently, we find that the susceptibility *χ*, the variance of the local order parameter *s*, also diverges around *R*_{c} in the limit of large system size (Fig. 4*B*). Applying finite-size scaling (25), we find that the measured susceptibility maximum *χ*_{max} diverges as a function of linear system size *L* (Fig. 2*C*, *Inset*) with exponent 1.71 ± 0.05, agreeing with the 2D Ising critical exponent ratio *γ*/*ν* = 7/4 for an equilibrium system with scalar order parameter and short-range interactions. Moreover, as shown in Fig. 4*C**,* the global order parameter exhibits singular behavior *M* ≈ |*R* − *R*_{c}|^{β} with *β* = 0.127 ± 0.003, again consistent with the pertinent equilibrium prediction 1/8. This agreement with expectations for equilibrium critical phenomena extends beyond static properties, encompassing even dynamics: The relaxation time of the system diverges as (*R* − *R*_{c})^{−zν} upon approaching the critical point from the supercritical regime (Fig. 4*D*), with *zν* = 3.7 ± 0.1, in agreement with *zν* = 4 − *η* = 3.75 for a conserved order parameter (26). The discovery of these highly nontrivial (24) power laws shows that the dynamical system quantitatively reflects all universal aspects of a continuous phase transition, including a consistent and complete set of critical exponents belonging to this known equilibrium universality class. We conjecture that a similar framework may apply to other steady-state nonequilibrium systems also, provided that the reasoning leading to designing this experimental system is satisfied, as, in an equilibrium system, all of these findings can be derived from a free-energy functional (25, 26). Although no quantitative derivation can be presented at this time, to discover this hint of such a functional in a nonequilibrium system is stimulating.

Furthermore, mapping a complete “phase diagram” parametrized by *R* and *ϕ*_{A}, Fig. 3*C* visually and quantitatively resembles the established phase diagram of binary fluids, sectioned by a binodal and a spinodal curve. Validating the simulations, experiments within the phase boundary display the expected interconnected patterns known as spinodal decomposition (Fig. 3*A* and Movie S5) at the critical composition, and nucleation-and-growth behavior (Movie S6) at off-critical compositions close to the phase boundary. At the initial stage when (in the case of thermal phase separation) most differences between these behaviors are expected, the differences here too are clear (Fig. 5*A*): Whereas the former displays patterns with a specific wavelength, the latter relies on random excitations to form a nucleus and hence lacks a characteristic length scale. Nascent nuclei fluctuate wildly, then merge and grow into spherical domains with time (Movies S6 and S7). Late-stage coarsening behavior of the spinodal pattern is also consistent with the equilibrium concept of surface tension driving domain growth with a power law exponent of 0.5 (Fig. 5 *B*–*D* and Fig. S3). The existence of an effective surface tension can be clearly seen from the spreading of a spherical droplet onto a larger domain (Fig. 5*B* and Movie S8). However, the capillary velocity (Fig. 5*D*, *Inset*) is seven orders of magnitude less than that of water in which these colloidal suspensions are immersed. Spinodal patterns at the critical composition, and nucleation patterns off the critical composition, present another pleasing analogy to equilibrium mixtures.

## Discussion

In conclusion, by designing this two-component system of driven Janus colloids, we have demonstrated the possibility of a nearly complete mapping between nonequilibrium phenomena and thermodynamics. Here, constant collisions powered by external fields play a role similar to that of thermal agitations, controlling both kinetics and steady states of the system. Akin to random forces in the Langevin approach, they cause the particles to diffuse at large time scales and display the fluctuation–dissipation relation. The diffusivity *D* measured from this long-term kinetics precisely matches the prediction based upon the short-term fluctuations (Fig. 2*C*). The linear dependence of *D* on the radius *R* of particle rotation further suggests that this radius plays the role of an effective temperature *T*, controlling the randomness of the system. More importantly, the effective temperature consistently governs the steady states: As *R* (or *T*) is decreased, the system undergoes a transition from a mixture to segregated domains similar to a thermal binary fluid, even with an archetypal phase diagram and 2D Ising critical exponents. This temperature dependence is different from that of the motility-induced phase behavior found in systems of linearly self-propelled particles (16, 27). There, with the ideal single-particle diffusivity *D*_{0} defined as effective temperature *T*, phase separation occurs as *D*_{0} (or *T*) is increased, the opposite of conventional thermal systems; this is because the diffusion of linearly self-propelled particles is primarily a combination of rotational fluctuations and particle motility. The motility neither introduces nor distributes randomness among the particles, making the diffusivity *D*_{0} an improper measure of the effective temperature. Recent work suggests that, instead, the motility gives rise to an effective pressure that causes particle aggregation (7). In contrast, in our system, randomness is induced solely by interparticle collisions and, in the one-phase regime, can be shared among particles, making it possible to define a universal effective temperature that governs both kinetics and phase behavior.

A devil’s advocate might argue that there is no need to go beyond the diffusivity *D* to define a temperature, as the Fokker–Planck or Smoluchowski equation that describes the microstructure in the suspension is completely specified by diffusivity and interparticle forces. However, we consider that defining an emergent temperature by the imposed radius of rotation *R* offers deeper understanding by connecting to the mature framework of statistical thermodynamics and enabling the thermodynamics-based predictions that our simulations and experiments validate. Furthermore, because no noise is introduced in our simulations, the effective temperature is not an arbitrary simulation input but instead emerges naturally from particle collisions.

We note limitations of the effective temperature concept based upon interparticle collisions. First, the concept depends on particle dynamics, so that particles moving at different radii do not share the same temperature. Second, the temperature also varies with local environment once phase separation sets in. Then, the effective temperature locally drops to zero and these systems crystallize (Movies S6–S8), unlike equilibrium systems where temperature is inherently uniform. This nonuniform temperature distribution has been also reported in other nonequilibrium systems (23, 27). Third, this effective temperature is only valid on long time scales, and is not necessarily reflected in instantaneous particle arrangements and kinetics. For instance, the mixed state at high effective temperature consists of lanes (Movie S2), unlike an equilibrium system. Moreover, during coarsening, collisions at the domain boundary drive particles to move in opposite directions parallel to the boundary (Fig. 6). Particle currents along domain boundaries were also reported in a recent work on the phase separation of oppositely spinning particles (23), but there the role of effective temperature was not discussed. The findings in this paper should generalize naturally to other dynamic components that perform three-dimensional periodic or chaotic motion that traces back to the origin with statistically significant numbers of collisions. However, these limitations underscore that care is needed to apply the effective temperature concept despite the nontrivial predictive capability demonstrated here.

## Materials and Methods

### Particle Synthesis.

Onto a planar submonolayer of monodisperse 3-µm silica particles (Tokuyama), obtained by drying 20 µL of 2 wt% particle suspension in deionized (DI) water onto a half glass slide (1.5 inch × 1 inch), a thin film of 4-nm Ni/10-nm Ti/5-nm SiO_{2} is sequentially deposited vertically using electron beam deposition in vacuum. The monolayer is washed thoroughly with DI water and isopropyl alcohol, and then sonicated in 20 mL of DI water to collect the particles in a 50-mL centrifuge tube.

### Setup.

A spatially homogeneous rotating magnetic field is generated by two orthogonal pairs of solenoids as described previously (28). Briefly, the two pairs of solenoids receive two sinusoidal voltages from a function generator, amplified by power amplifiers, with π/2 phase difference between each signal to produce an in-plane rotating field *B*.

The electric field is applied by sandwiching the particles between two coverslips coated with indium tin oxide (ITO) from SPI Supplies (29). The ITO-coated coverslips are further coated with 25 nm of SiO_{2} using electron beam deposition to electrostatically prevent the particles from surface adsorption. A thin strip on one side of the ITO-coated coverslip is left uncoated during the SiO_{2} deposition, and is used later to connect to a function generator (Agilent 33522A) via copper tapes. A square wave of 5 kHz and 7 V is used for all experiments. The two ITO-coated coverslips are separated by a spacer (SecureSeal; Grace Biolab) about 120 μm thick, with a 9-mm-diameter hole.

Movies are taken in a customized microscope using an LED light source (MCWHL2; Thorlabs) and a CMOS camera (5012M GigE; Edmund Optics). For a high-magnification view in which the Janus features can be resolved, a 50× long-working-distance objective (numerical aperture 0.55; Mitutoyo) is used. For large-scale views, a 10× long-working-distance objective (numerical aperture 0.28; Mitutoyo) is used. In this case, it is not possible to differentiate the silica/metal sides of each individual sphere. Instead, we use “tilted illumination” and shine light on the particle from an oblique angle, so that a particle appears dark or bright depending on whether the metal side is facing toward or away from the light source. We then use stroboscopic sampling, only analyzing one frame per cycle, such that one species is consistently brighter than the other. When we zoom in to take a magnified view of the system, the optics is well aligned, allowing us to clearly differentiate between the metal-coated and the silica hemisphere.

### Experimental Procedure.

Ten microliters of concentrated particle suspension is dropped onto an ITO-coated coverslip with spacer, and then sandwiched with another ITO coverslip. This procedure produces a monolayer of colloids with area fraction between 0.35 and 0.45. Repeated experiments can be done on one sample for 2 to 3 h, after which significant particle clustering or substrate adsorption occurs. Before each run, a constant DC magnetic field is applied in the image plane to cause particles to assemble into zigzag chains, which ensures a completely mixed initial state without memory of earlier states. A rotating magnetic field is then applied with a strength of 5 mT and a frequency ranging from 16 Hz to 0.05 Hz, followed immediately by an ac electric field (5 kHz, 58 V/mm); this we define as *t* = 0 s. Movies are taken up to 20 min for each condition, at a rate of 16 frames per second.

To change the ratio between the two species, a pulse of additional *z*-axis magnetic field is applied. Although the precise mechanism of the conversion is unclear, we are able to roughly control the composition by varying the strength and duration of this *z* field.

### Image Processing for Large View.

We first subtract the background light gradient due to the intentional optical misalignment, obtained by averaging over a whole movie. The image is then smoothed and binarized, with each pixel now having a local order parameter *s*(*r*) (defined in *Analysis of Phase Transition*) of value +1 or −1. Spatial correlation *C*(*r*), defined as 〈*s*(0)·*s*(*r*)〉, is calculated from the binarized image. Correlation length *ξ* is defined as the *r* value where *C*(*r*) first passes zero. Linear interpolation between the last positive value and first negative value is used to improve accuracy. Structure factor is obtained by applying fast Fourier transform to the image and taking a radial average.

### Simulation Procedure.

We model this overdamped and fully driven quasi-2D system via athermal Brownian dynamics simulations. Confined to the *x*–*y* plane, Janus particles of diameter *σ* = 3 μm with area fraction 0.37 (to match the experiment) move in a square, periodic domain of side length 200*σ*. Their positions ** r**(

*t*) and orientations Ω

**(**

*t*) follow the master equations

where ** F** and

**are, respectively, the force and the torque exerted on the particle, and**

*T**ζ*

_{t}and

*ζ*

_{r}refer to the translational and rotational drag coefficients. Here reduced units are used, with the particle mass

*m*= 3.74 × 10

^{−14}kg as mass unit, the particle diameter

*σ*as length unit, and

*τ*= 10

^{−3}s as time unit. Given the dynamic viscosity

*μ*of the water, those spherical particles have

*ζ*

_{t}= 3π

*μσ*= 672.8 and

*ζ*

_{r}= π

*μσ*

^{3}= 224.3. According to the reciprocal theory (30), the action of the electric and magnetic fields in experiment can be modeled as a driven circular motion caused by an external force

*F*

_{ex}=

*ζ*

_{t}

*v*= 4.7 (with

*ν*≈ 21 μm/s as the measured terminal swimming velocity) and an external torque

*T*

_{ex}=

*ζ*

_{r}

*ω*(with

*ω*as the tunable angular frequency). Two types of particles are considered, with oppositely directed driving forces at each instant.

In addition to the external driving, interactions between the particles also affect their motions. Excluded-volume effects are implemented as shifted−truncated Lennard-Jones interactions with *σ*_{LJ} = 1, *ε*_{LJ} = 1, and cutoff at *r*_{c} = 2^{1/6}. In experiment, due to the *z*-directional electric field, the particles experience an induced dipolar force, which is isotropic in the *x*–*y* plane and decays as *A*⋅*r*^{−4}. Here we set the strength *A* = 7.5 so that the dipolar interaction balances the driving force at *r*_{c}. Note that this is a simplified version of our previous model (29), where both sides of a Janus sphere carry an induced electric dipole shifted away from the particle center, and the resultant off-centered dipolar interaction can reorient the particle. In the current system, the magnetic field applied to control particle orientation overwhelms this interaction-induced reorientation mechanism. Therefore, we only consider the overall induced electric dipole moment on each particle and the resultant soft repulsion between particles.

All simulations start with a random configuration and proceed for at least 10^{9} steps with time step d*t* = *τ*. This duration corresponds to around 10^{5} rotation cycles (equivalent to 50 h in experiment), which we found to be sufficient for the system to reach a steady state.

### Analysis of Phase Transition.

To quantify the phase transition, we define the order parameter as the composition difference *s* = *ϕ*_{A} − *ϕ*_{B}, which is conserved globally but fluctuates locally in the system. Here we divide the entire system into relatively small subcells of side length *L* = 12.5*σ* and record their composition fluctuations, which mimic those present in the grand canonical ensemble. The histogram yields the order-parameter distribution, which changes from a Gaussian distribution centered at 0 to a symmetric bimodal distribution as *R* decreases at symmetric composition (Fig. S4). The peak position denotes the global order parameter *M* of the system, specifically,

It displays singular behavior at the transition point *R*_{c} (Fig. 4*C*). The susceptibility is defined as

To study the divergence of *χ*, we perform a finite-size scaling analysis by systematically varying the subcell size *L* from 5.6*σ* to 20*σ* and exploring the rise of the maximum susceptibility *χ*_{max} with *L* (Fig. 2*C*).

## Acknowledgments

At the Institute for Basic Science Center for Soft and Living Matter, S.G. acknowledges support by the Institute for Basic Science, project code IBS-R020-D1. This work was supported by the U.S. Department of Energy, Division of Basic Energy Sciences, under Award DE-FG02-07ER46471, through the Frederick Seitz Materials Research Laboratory at the University of Illinois at Urbana–Champaign (to J.Y. and S.G.) and by the National Science Foundation (NSF) under Award Nos. DMR-1121262 and DMR-1610796 (to M.H. and E.L.). We acknowledge support from NSF CBET-0853737 for equipment and from the Quest high-performance computing facility at Northwestern University (to M.H. and E.L.). J.Y. holds a Career Award at the Scientific Interface from the Burroughs Wellcome Fund.

## Footnotes

↵

^{1}M.H. and J.Y. contributed equally to this work.- ↵
^{2}To whom correspondence may be addressed. Email: sgranick{at}ibs.re.kr or luijten{at}northwestern.edu.

Author contributions: J.Y., S.G., and E.L. designed research; M.H. and J.Y. performed research; M.H., J.Y., S.G., and E.L. analyzed data; and M.H., J.Y., S.G., and E.L. wrote the paper.

Reviewers: D.B., Ecole Normale Supérieure de Lyon; and A.C., Columbia University.

The authors declare no conflict of interest.

This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1706702114/-/DCSupplemental.

Freely available online through the PNAS open access option.

## References

- ↵
- ↵
- ↵.
- Stenhammar J,
- Wittkowski R,
- Marenduzzo D,
- Cates ME

- ↵
- ↵.
- Dieterich E,
- Camunas-Soler J,
- Ribezzi-Crivellari M,
- Seifert U,
- Ritort F

- ↵
- ↵.
- Takatori SC,
- Brady JF

- ↵
- ↵
- ↵.
- Wang S,
- Wolynes PG

- ↵.
- Edwards SF

- ↵
- ↵.
- Song C,
- Wang P,
- Makse HA

- ↵
- ↵.
- Bechinger C, et al.

- ↵
- ↵
- ↵
- ↵.
- Vissers T, et al.

- ↵.
- Van Kampen NG

- ↵.
- von Smoluchowski M

- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵.
- Yan J, et al.

- ↵.
- Happel J,
- Brenner H

## Citation Manager Formats

## Article Classifications

- Physical Sciences
- Applied Physical Sciences