Source author record

Monica Olvera de la Cruz

Monica Olvera de la Cruz 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

21works
12topics
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

21 published item(s)

preprint2023arXiv

Modulation of ionic conduction using polarizable surfaces

Hybrid ionic-electronic conductors have the potential to generate memory effects and neuronal behavior. The functionality of these mixed materials depends on ion motion through thin polarizable channels. Here, we explore different polarization models to show that the current and conductivity of electrolytes is higher when confined by conductors than by dielectrics. We find non-linear currents in both dielectrics and conductors, and we recover the known linear (Ohmic) result only in the two-dimensional limit between conductors. We show that the polarization charge location impacts electrolyte structure and transport properties. This work suggests a mechanism to induce memristor hysteresis loops using conductor-dielectric switchable materials.

preprint2022arXiv

Stokesian Dynamics with odd viscosity

Stokesian Dynamics is a well-established computational method for simulating dynamics of many particles suspended in a conventional passive fluid medium. Active fluids composed of self-propelling particles with broken time reversal symmetry permit the emergence of a so-called odd viscosity. In this work, we extended the conventional Stokesian Dynamics formalism to incorporate the additional hydrodynamic effects due to odd viscosity, which enables simulating collective behaviors of many particles suspended in an active fluid medium with both even viscosity and odd viscosity.

preprint2021arXiv

Polar state memory in active fluids

Spontaneous emergence of correlated states such as flocks and vortices are prime examples of remarkable collective dynamics and self-organization observed in active matter. The formation of globally correlated polar states in geometrically confined systems proceeds through the emergence of a macroscopic steadily rotating vortex that spontaneously selects a clockwise or counterclockwise global chiral state. Here, we reveal that a global vortex formed by colloidal rollers exhibits state memory. The information remains stored even when the energy injection is ceased and the activity is terminated. We show that a subsequent formation of the collective states upon re-energizing the system is not random. We combine experiments and simulations to elucidate how a combination of hydrodynamic and electrostatic interactions leads to hidden asymmetries in the local particle positional order encoding the chiral state of the system. The stored information can be accessed and exploited to systematically command subsequent polar states of active liquid through temporal control of the activity. With the chirality of the emergent collective states controlled on-demand, active liquids offer new possibilities for flow manipulation, transport, and mixing at the microscale.

preprint2020arXiv

Dynamics of a driven confined polyelectrolyte solution

The transport of polyelectrolytes confined by oppositely charged surfaces and driven by a constant electric field is of interest in studies of DNA separation according to size. Using molecular dynamics simulations that include surface polarization effect, we find that the mobilities of the polyelectrolytes and their counterions change non-monotonically with the confinement surface charge density. For an optimum value of the confinement charge density, efficient separation of polyelectrolytes can be achieved over a wide range of polyelectrolyte charge due to the differential friction imparted by the oppositely charged confinement on the polyelectrolyte chains. Furthermore, by altering the placement of the charged confinement counterions, enhanced polyelectrolyte separation can be achieved by utilizing surface polarization effect due to dielectric mismatch between the media inside and outside the confinement.

preprint2020arXiv

Multicanonical Monte Carlo Ensemble Growth Algorithm

We present a novel Ensemble Monte Carlo Growth method to sample the equilibrium thermodynamic properties of random chains. The method is based on the multicanonical technique of computing the density of states in the energy space. Such a quantity is temperature independent, and therefore microcanonical and canonical thermodynamic quantities, including the free energy, entropy, and thermal averages, can be obtained by re-weighting with a Boltzmann factor. The algorithm we present combines two approaches: the first is the Monte Carlo ensemble growth method, where a "population" of samples in the state space is considered, as opposed to traditional sampling by long random walks, or iterative single-chain growth. The second is the flat-histogram Monte Carlo, similar to the popular Wang-Landau sampling, or to multicanonical chain-growth sampling. We discuss the performance and relative simplicity of the proposed algorithm, and we apply it to known test cases.

preprint2020arXiv

Sublattice Melting in Binary Superionic Colloidal Crystals

In superionic compounds one component pre-melts providing high ionic conductivity to solid state electrolytes. Here, we find sublattice melting in colloidal crystals of oppositely charged particles that are highly asymmetric in size and charge in salt solutions. The small particles in ionic compounds melt when the temperature increases forming a superionic phase. These delocalized small particles in a crystal of large oppositely charged particles, in contrast to superionic phases in atomic systems, form crystals with non-electroneutral stoichiometric ratios. This generates structures with multiple domains of ionic crystals in percolated superionic phases with adjustable stoichiometries.

preprint2019arXiv

Control of Ionic Mobility via Charge Size Asymmetry in Random Ionomers

Solid polymer electrolytes are considered a promising alternative to traditional liquid electrolytes in energy storage applications because of their good mechanical properties, and excellent thermal and chemical stability. A gap, however, still exists in understanding ion transport mechanisms and improving ion transport in solid polymer electrolytes. Therefore, it is crucial to bridge composition--structure and structure--property relationships. Here we demonstrate that size asymmetry, $λ$, represented by the ratio of counterion to charged monomer size, plays a key role in both the nanostructure and in the ionic dynamics. More specifically, when the nanostructure is modified by the external electric field such that the mobility cannot be described by linear response theory, two situations arise. The ionic mobility increases as $λ$ decreases (small counterions) in the weak electrostatics (high dielectric constant) regime. Whereas in systems with strong electrostatic interactions, ionomers with higher size symmetry ($λ\approx 1$) display higher ionic mobility. Moreover, ion transport is found to be dominated by the hopping of the ions and not by moving ionic clusters (also known as "vehicular'' charge transport). These results serve as a guide for designing ion-containing polymers for ion transport related applications.

preprint2015arXiv

Coulomb energy of uniformly-charged spheroidal shell systems

We provide exact expressions for the electrostatic energy of uniformly-charged prolate and oblate spheroidal shells. We find that uniformly-charged prolate spheroids of eccentricity greater than 0.9 have lower Coulomb energy than a sphere of the same area. For the volume-constrained case, we find that a sphere has the highest Coulomb energy among all spheroidal shells. Further, we derive the change in the Coulomb energy of a uniformly-charged shell due to small, area-conserving perturbations on the spherical shape. Our perturbation calculations show that buckling-type deformations on a sphere can lower the Coulomb energy. Finally, we consider the possibility of counterion condensation on the spheroidal shell surface. We employ a Manning-Oosawa two-state model approximation to evaluate the renormalized charge and analyze the behavior of the equilibrium free energy as a function of the shell's aspect ratio for both area-constrained and volume-constrained cases. Counterion condensation is seen to favor the formation of spheroidal structures over a sphere of equal area for high values of shell volume fractions.

preprint2015arXiv

Effective charges and virial pressure of concentrated macroion solutions

The stability of colloidal suspensions is crucial in a wide variety of processes including the fabrication of photonic materials and scaffolds for biological assemblies. The ionic strength of the electrolyte that suspends charged colloids is widely used to control the physical properties of colloidal suspensions. The extensively used two-body Derjaguin-Landau-Verwey-Overbeek (DLVO) approach allows for a quantitative analysis of the effective electrostatic forces between colloidal particles. DLVO relates the ionic double-layers, which enclose the particles, to their effective electrostatic repulsion. Nevertheless, the double layer is distorted at high macroion volume fractions. Therefore, DLVO cannot describe the many-body effects that arise in concentrated suspensions. We show that this problem can be largely resolved by identifying effective point charges for the macroions using cell theory. This extrapolated point charge (EPC) method assigns effective point charges in a consistent way, taking into account the excluded volume of highly charged macroions at any concentration, and thereby naturally accounting for high volume fractions in both salt-free and added-salt conditions. We provide an analytical expression for the effective pair potential and validate the EPC method by comparing molecular dynamics simulations of macroions and monovalent microions that interact via Coulombic potentials to simulations of macroions interacting via the derived EPC effective potential. The simulations reproduce the macroion-macroion spatial correlation and the virial pressure obtained with the EPC model. Our findings provide a route to relate the physical properties such as pressure in systems of screened-Coulomb particles to experimental measurements.

preprint2014arXiv

Electrostatics-driven shape transitions in soft shells

Manipulating the shape of nanoscale objects in a controllable fashion is at the heart of designing materials that act as building blocks for self-assembly or serve as targeted drug delivery carriers. Inducing shape deformations by controlling external parameters is also an important way of designing biomimetic membranes. In this paper, we demonstrate that electrostatics can be used as a tool to manipulate the shape of soft, closed membranes by tuning environmental conditions such as the electrolyte concentration in the medium. Using a molecular dynamics-based simulated annealing procedure, we investigate charged elastic shells that do not exchange material with their environment, such as elastic membranes formed in emulsions or synthetic nanocontainers. We find that by decreasing the salt concentration or increasing the total charge on the shell's surface, the spherical symmetry is broken, leading to the formation of ellipsoids, discs, and bowls. Shape changes are accompanied by a significant lowering of the electrostatic energy and a rise in the surface area of the shell. To substantiate our simulation findings, we show analytically that a uniformly charged disc has a lower Coulomb energy than a sphere of the same volume. Further, we test the robustness of our results by including the effects of charge renormalization in the analysis of the shape transitions and find the latter to be feasible for a wide range of shell volume fractions.

preprint2014arXiv

Free-energy functional of the electronic potential for Schrödinger-Poisson theory

In the study of model electronic device systems where electrons are typically under confinement, a key obstacle is the need to iteratively solve the coupled Schrödinger-Poisson (SP) equation. It is possible to bypass this obstacle by adopting a variational approach and obtaining the solution of the SP equation by minimizing a functional. Further, using molecular dynamics methods that treat the electronic potential as a dynamical variable, the functional can be minimized on the fly in conjunction with the update of other dynamical degrees of freedom leading to considerable reduction in computational costs. But such approaches require access to a true free-energy functional, one that evaluates to the equilibrium free energy at its minimum. In this paper, we present a variational formulation of the Schrödinger-Poisson (SP) theory with the needed free-energy functional of the electronic potential. We apply our formulation to semiconducting nanostructures and provide the expression of the free-energy functional for narrow channel quantum wells where the local density approximation yields accurate physics and for the case of wider channels where Thomas-Fermi approximation is valid.

preprint2013arXiv

A variational formulation of electrostatics in a medium with spatially varying dielectric permittivity

In biological and synthetic materials, many important processes involve charges that are present in a medium with spatially varying dielectric permittivity. To accurately understand the role of electrostatic interactions in such systems, it is important to take into account the spatial dependence of the permittivity of the medium. However, due to the ensuing theoretical and computational challenges, this inhomogeneous dielectric response of the medium is often ignored or excessively simplified. We develop a variational formulation of electrostatics to accurately investigate systems that exhibit this inhomogeneous dielectric response. Our formulation is based on a true energy functional of the polarization charge density. The defining characteristic of a true energy functional is that at its minimum it evaluates to the actual value of the energy; this is a feature not found in many commonly used electrostatic functionals. We explore in detail the charged systems that exhibit sharp discontinuous change in dielectric permittivity, and we show that for this case our functional reduces to a functional of only the surface polarization charge density. We apply this reduced functional to study model problems for which analytical solutions are well known. We demonstrate, in addition, that the functional has many properties that make it ideal for use in molecular dynamics simulations.

preprint2013arXiv

Electrostatic repulsion-driven crystallization model arising from filament networks

The crystallization of bundles in filament networks interacting via long-range repulsions in confinement is described by a phenomenological model. The model demonstrates the formation of the hexagonal crystalline order via the interplay of the confinement potential and the filament-filament repulsion. Two distinct crystallization mechanisms in the short- and large- screening length regimes are discussed, and the phase diagram is obtained. Simulation of large bundles predicts the existence of topological defects within the bundled filaments. This electrostatic repulsion-driven crystallization model arising from studying filament networks can even find a more general context extending to charged colloidal systems.

preprint2013arXiv

Free-energy functionals of the electrostatic potential for Poisson-Boltzmann theory

In simulating charged systems, it is often useful to treat some ionic components of the system at the mean-field level and solve the Poisson-Boltzmann (PB) equation to get their respective density profiles. The numerically intensive task of solving the PB equation at each step of the simulation can be bypassed using variational methods that treat the electrostatic potential as a dynamic variable. But such approaches require the access to a true free-energy functional; a functional that not only provides the correct solution of the PB equation upon extremization, it also evaluates to the true free energy of the system at its minimum. Moreover, the numerical efficiency of such procedures is further enhanced if the free-energy functional is local and is expressed in terms of the electrostatic potential. Existing PB functionals of the electrostatic potential, while possessing the local structure, are not free-energy functionals. We present a variational formulation with a local free-energy functional of the potential. In addition, we also construct a nonlocal free-energy functional of the electrostatic potential. These functionals are suited for employment in simulation schemes based on the ideas of dynamical optimization.

preprint2013arXiv

Generating true minima in constrained variational formulations via modified Lagrange multipliers

Variational principles are important in the investigation of large classes of physical systems. They can be used both as analytical methods as well as starting points for the formulation of powerful computational techniques such as dynamical optimization methods. Systems with charged objects in dielectric media and systems with magnetically active particles are important examples. In these examples and other important cases, the variational principles describing the system are required to obey a number of constraints. These constraints are implemented within the variational formulation by means of Lagrange multipliers. Such constrained variational formulations are in general not unique. For the application of efficient simulation methods, one must find specific formulations that satisfy a number of important conditions. An often required condition is that the functional be positive-definite, in other words, its extrema be actual minima. In this article, we present a general approach to attack the problem of finding, among equivalent variational functionals, those that generate true minima. The method is based on the modification of the Lagrange multiplier which allows us to generate large families of effective variational formulations associated with a single original constrained variational principle. We demonstrate its application to different examples and, in particular, to the important cases of Poisson and Poisson-Boltzmann equations. We show how to obtain variational formulations for these systems with extrema that are always minima.

preprint2013arXiv

Shapes of pored membranes

We study the shapes of pored membranes within the framework of the Helfrich theory under the constraints of fixed area and pore size. We show that the mean curvature term leads to a budding- like structure, while the Gaussian curvature term tends to flatten the membrane near the pore; this is corroborated by simulation. We propose a scheme to deduce the ratio of the Gaussian rigidity to the bending rigidity simply by observing the shape of the pored membrane. This ratio is usually difficult to measure experimentally. In addition, we briefly discuss the stability of a pore by relaxing the constraint of a fixed pore size and adding the line tension. Finally, the flattening effect due to the Gaussian curvature as found in studying pored membranes is extended to two-component membranes. We find that sufficiently high contrast between the components' Gaussian rigidities leads to budding which is distinct from that due to the line tension.

preprint2013arXiv

Simulation of Charged Systems in Heterogeneous Dielectric Media via a True Energy Functional

For charged systems in heterogeneous dielectric media, a key obstacle for molecular dynamics (MD) simulations is the need to solve the Poisson equation in the media. This obstacle can be bypassed using MD methods that treat the local polarization charge density as a dynamic variable, but such approaches require access to a true free energy functional; one that evaluates to the equilibrium electrostatic energy at its minimum. In this letter, we derive the needed functional. As an application, we develop a Car-Parrinello MD method for the simulation of free charges present near a spherical emulsion droplet separating two immiscible liquids with different dielectric constants. Our results show the presence of non-monotonic ionic profiles in the dielectric with lower dielectric constant.

preprint2012arXiv

A curvature-driven effective attraction in multicomponent membranes

We study closed liquid membranes that segregate into three phases due to differences in the chemical and physical properties of its components. The shape and in-plane membrane arrangement of the phases are coupled through phase-specific bending energies and line tensions. We use simulated annealing Monte Carlo simulations to find low-energy structures, allowing both phase arrangement and membrane shape to relax. The three-phase system is the simplest one in which there are multiple interface pairs, allowing us to analyze interfacial preferences and pairwise distinct line tensions. We observe the system's preference for interface pairs that maximize differences in spontaneous curvature. From a pattern selection perspective, this acts as an effective attraction between phases of most disparate spontaneous curvature. We show that this effective attraction is robust enough to persist even when the interface between these phases is the most penalized by line tension. This effect is driven by geometry and not by any explicit component-component interaction.

preprint2012arXiv

Nonlinear elastic model for faceting of vesicles with soft grain boundaries

We use an elastic model to explore faceting of solid-wall vesicles with elastic heterogeneities. We show that faceting occurs in regions where the vesicle wall is softer, such as areas of reduced wall thicknesses or concentrated in crystalline defects. The elastic heterogeneities are modeled as a second component with reduced elastic parameters. Using simulated annealing Monte Carlo simulations we obtain the vesicle shape by optimizing the distributions of facets and boundaries. Our model allows us to reduce the effects of the residual stress generated by crystalline defects, and reveals a robust faceting mechanism into polyhedra other than the icosahedron.

preprint2005arXiv

Phase equilibria and charge fractionation in polydisperse polyelectrolyte solutions

A new type of phase separation in the polyelectrolyte solutions consisting of several types of charged macromolecules differing in their degree of ionization is predicted via a general thermodynamic consideration. We show that even a small difference in the degree of ionization of otherwise equivalent components results in their spatial separation occurring upon decreasing the temperature much earlier than precipitation of the components from the solution. Some implications of charge fractionation in biological processes are discussed.

preprint2002arXiv

Structure and thermodynamics of associating rods solutions

Thermoreversible sol-gel transitions in solutions of rod-like associating polymers are analyzed by computer simulations and by mean field models. The sol-gel transition is determined by the divergence of the cluster weight average. The analytically determined sol-gel transition is in good agreement with the simulation results. At low temperatures we observe a peak in the heat capacity, which maximum is associated with the precipitation transition. The gelation transition is sensitive to the number of associating groups per rod but nearly insensitive to the spatial distribution of associating groups around the rod. The precipitation is strongly dependent on both the number and distribution of associating groups per rod. We find negligible nematic orientational order at the gelation and precipitation transitions.