# Monte Carlo computer simulation of sedimentation of charged hard spherocylinders

## Abstract

In this article we present a NVT Monte Carlo computer simulation study of sedimentation of an electroneutral mixture of oppositely charged hard spherocylinders (CHSC) with aspect ratio L/σ = 5, where L and σ are the length and diameter of the cylinder and hemispherical caps, respectively, for each particle. This system is an extension of the restricted primitive model for spherical particles, where L/σ = 0, and it is assumed that the ions are immersed in an structureless solvent, i.e., a continuum with dielectric constant D. The system consisted of N = 2000 particles and the Wolf method was implemented to handle the coulombic interactions of the inhomogeneous system. Results are presented for different values of the strength ratio between the gravitational and electrostatic interactions, Γ = (mgσ)/(e{sup 2}/Dσ), where m is the mass per particle, e is the electron's charge and g is the gravitational acceleration value. A semi-infinite simulation cell was used with dimensions L{sub x} ≈ L{sub y} and L{sub z} = 5L{sub x}, where L{sub x}, L{sub y}, and L{sub z} are the box dimensions in Cartesian coordinates, and the gravitational force acts along the z-direction. Sedimentation effects were studied by looking at every layer formedmore »

- Authors:

- Unidad Académica de Física, Universidad Autónoma de Zacatecas, Calzada Solidaridad esq. Paseo, La Bufa s/n, 98060 Zacatecas, Zacatecas, México (Mexico)
- Departamento de Ingeniería Física, División de Ciencias e Ingenierías, Campus León, Universidad de Guanajuato, Loma del Bosque 103, Lomas del Campestre, 37150 León, Guanajuato, México (Mexico)

- Publication Date:

- OSTI Identifier:
- 22419966

- Resource Type:
- Journal Article

- Resource Relation:
- Journal Name: Journal of Chemical Physics; Journal Volume: 141; Journal Issue: 4; Other Information: (c) 2014 AIP Publishing LLC; Country of input: International Atomic Energy Agency (IAEA)

- Country of Publication:
- United States

- Language:
- English

- Subject:
- 37 INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL CHEMISTRY; COMPUTERIZED SIMULATION; ELECTRONS; GRAVITATIONAL FIELDS; INTERACTIONS; MONTE CARLO METHOD; PARTICLES; PERMITTIVITY; SEDIMENTATION; SOLVENTS

### Citation Formats

```
Viveros-Méndez, P. X., E-mail: xviveros@fisica.uaz.edu.mx, Aranda-Espinoza, S., and Gil-Villegas, Alejandro.
```*Monte Carlo computer simulation of sedimentation of charged hard spherocylinders*. United States: N. p., 2014.
Web. doi:10.1063/1.4890819.

```
Viveros-Méndez, P. X., E-mail: xviveros@fisica.uaz.edu.mx, Aranda-Espinoza, S., & Gil-Villegas, Alejandro.
```*Monte Carlo computer simulation of sedimentation of charged hard spherocylinders*. United States. doi:10.1063/1.4890819.

```
Viveros-Méndez, P. X., E-mail: xviveros@fisica.uaz.edu.mx, Aranda-Espinoza, S., and Gil-Villegas, Alejandro. Mon .
"Monte Carlo computer simulation of sedimentation of charged hard spherocylinders". United States.
doi:10.1063/1.4890819.
```

```
@article{osti_22419966,
```

title = {Monte Carlo computer simulation of sedimentation of charged hard spherocylinders},

author = {Viveros-Méndez, P. X., E-mail: xviveros@fisica.uaz.edu.mx and Aranda-Espinoza, S. and Gil-Villegas, Alejandro},

abstractNote = {In this article we present a NVT Monte Carlo computer simulation study of sedimentation of an electroneutral mixture of oppositely charged hard spherocylinders (CHSC) with aspect ratio L/σ = 5, where L and σ are the length and diameter of the cylinder and hemispherical caps, respectively, for each particle. This system is an extension of the restricted primitive model for spherical particles, where L/σ = 0, and it is assumed that the ions are immersed in an structureless solvent, i.e., a continuum with dielectric constant D. The system consisted of N = 2000 particles and the Wolf method was implemented to handle the coulombic interactions of the inhomogeneous system. Results are presented for different values of the strength ratio between the gravitational and electrostatic interactions, Γ = (mgσ)/(e{sup 2}/Dσ), where m is the mass per particle, e is the electron's charge and g is the gravitational acceleration value. A semi-infinite simulation cell was used with dimensions L{sub x} ≈ L{sub y} and L{sub z} = 5L{sub x}, where L{sub x}, L{sub y}, and L{sub z} are the box dimensions in Cartesian coordinates, and the gravitational force acts along the z-direction. Sedimentation effects were studied by looking at every layer formed by the CHSC along the gravitational field. By increasing Γ, particles tend to get more packed at each layer and to arrange in local domains with an orientational ordering along two perpendicular axis, a feature not observed in the uncharged system with the same hard-body geometry. This type of arrangement, known as tetratic phase, has been observed in two-dimensional systems of hard-rectangles and rounded hard-squares. In this way, the coupling of gravitational and electric interactions in the CHSC system induces the arrangement of particles in layers, with the formation of quasi-two dimensional tetratic phases near the surface.},

doi = {10.1063/1.4890819},

journal = {Journal of Chemical Physics},

number = 4,

volume = 141,

place = {United States},

year = {Mon Jul 28 00:00:00 EDT 2014},

month = {Mon Jul 28 00:00:00 EDT 2014}

}