BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Akira Onuki

Akira Onuki appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.softcond-mat.stat-mechcond-mat.mtrl-scicond-mat.dis-nncond-mat.quant-gascond-mat.str-elphysics.chem-ph

Catalog footprint

What is connected

42works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

42 published item(s)

preprint2020arXiv

Theory of electrolytes including steric, attractive, and hydration interactions

We present a continuum theory of electrolytes composed of a waterlike solvent and univalent ions. First, we start with a density functional $\cal F$ for the coarse-grained solvent, cation, and anion densities, including the Debye-Hückel free energy, the Coulombic interaction, and the direct interactions among these three components. These densities fluctuate obeying the distribution $\propto \exp(- {\cal F}/k_BT)$. Eliminating the solvent density deviation in $\cal F$, we obtain the effective non-Coulombic interactions among the ions, which consist of the direct ones and the solvent-mediated ones. %where the latter are written in terms of the ion volumes %and are inversely proportional to the solvent compressibility. We then derive general expressions for the ion correlation, the apparent partial volume, and the activity and osmotic coefficients up to linear order in the average salt density $n_{\rm s}$. Secondly, we perform numerical analysis using the Mansoori-Carnahan-Starling-Leland model $[$J. Chem. Phys. {\bf 54}, 1523 (1971)$]$ for three-component hardspheres. The effective interactions sensitively depend on the cation and anion sizes due to competition between the steric and hydration effects, which are repulsive between small-large ion pairs and attractive between symmetric pairs. These agree with previous experiments and Collins' rule $[$Biophys. J. {\bf 72}, 65 (1997)$]$. We also give simple approximate expressions for the ionic interaction coefficients valid for any ion sizes.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Ionization at a solid-water interface in an applied electric field: Charge regulation

We investigate ionization at a solid-water interface in applied electric field. We attach an electrode to a dielectric film bearing silanol or carboxyl groups with an areal density $Γ_0$, where the degree of dissociation $α$ is determined by the proton density in water close to the film. We show how $α$ depends on the density $n_0$ of NaOH in water and the surface charge density $σ_m$ on the electrode. For $σ_m>0$, the protons are expelled away from the film, leading to an increase in $α$. In particular, in the range $0<σ_m<eΓ_0$, self-regulation occurs to realize $α\cong σ_m/eΓ_0 $ for $n_0\ll n_c$, where $n_c$ is $0.01$ mol$/$L for silica surfaces and is $2\times 10^{-5}$ mol$/$L for carboxyl-bearing surfaces. We also examine the charge regulation with decreasing the cell thickness $H$ below the Debye length $κ^{-1}$, where a crossover occurs at the Gouy-Chapman length. In particular, when $σ_m \sim eΓ_0$ and $H\ll κ^{-1}$, the surface charges remain only partially screened by ions, leading to an electric field in the interior.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Structure Formation due to Antagonistic Salts

Antagonistic salts are composed of hydrophilic and hydrophobic ions. In a mixture solvent (water-oil) such ion pairs are preferentially attracted to water or oil, giving rise to a coupling between the charge density and the composition. First, they form a large electric double layer at a water-oil interface, reducing the surface tension and producing mesophases. Here, the cations and anions are loosely bound by the Coulomb attraction across the interface on the scale of the Debye screening length. Second, on solid surfaces, hydrophilic (hydrophobic) ions are trapped in a water-rich (oil-rich) adsorption layer, while those of the other species are expelled from the layer. This yields a solvation mechanism of local charge separation near a solid. In particular, near the solvent criticality, disturbances around solid surfaces can become oscillatory in space. In mesophases, we calculate periodic structures, which resemble those in experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Bubble formation in water with addition of a hydrophobic solute

We show that phase separation can occur in a one-component liquid outside its coexistence curve (CX) with addition of a small amount of a solute. The solute concentration at the transition decreases with increasing the difference of the solvation chemical potential between liquid and gas. As a typical bubble-forming solute, we consider ${\rm O}_2$ in ambient liquid water, which exhibits mild hydrophobicity and its critical temperature is lower than that of water. Such a solute can be expelled from the liquid to form gaseous domains while the surrounding liquid pressure is higher than the saturated vapor pressure $p_{cx}$. This solute-induced bubble formation is a first-order transition in bulk and on a partially dried wall, while a gas film grows continuously on a completely dried wall. We set up a bubble free energy $ΔG$ for bulk and surface bubbles with a small volume fraction $ϕ$. It becomes a function of the bubble radius $R$ under the Laplace pressure balance. Then, for sufficiently large solute densities above a threshold, $ΔG$ exhibits a local maximum at a critical radius and a minimum at an equilibrium radius. We also examine solute-induced nucleation taking place outside CX, where bubbles larger than the critical radius grow until attainment of equilibrium.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Density functional theory of gas-liquid phase separation in dilute binary mixtures

We examine statics and dynamics of phase-separated states of dilute binary mixtures using density functional theory. In our systems, the difference in the salvation chemical potential $Δμ_s$ between liquid and gas is considerably larger than the thermal energy $k_BT$ for each solute particle and the attractive interaction among the solute particles is weaker than that among the solvent particles. In these conditions, the saturated vapor pressure increases by an amount equal to the solute density in liquid multiplied by the large factor $k_BT \exp(Δμ_s/k_BT)$. As a result, phase separation is induced at low solute densities in liquid and the new phase remains in gaseous states, while the liquid pressure is outside the coexistence curve of the solvent. This explains the widely observed formation of stable nanobubbles in ambient water with a dissolved gas. We calculate the density and stress profiles across planar and spherical interfaces, where the surface tension decreases with increasing the interfacial solute adsorption. We realize stable solute-rich bubbles with radius about 30 nm, which minimize the free energy functional. We then study dynamics around such a bubble after a decompression of the surrounding liquid, where the bubble undergoes a damped oscillation. In addition, we present some exact and approximate expressions for the surface tension and the interfacial stress tensor.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Fluctuations of local electric field and dipole moments in water between metal walls

We examine the thermal fluctuations of the local electric field $E_k^{\rm loc}$ and the dipole moment $μ_k$ in liquid water at $T=298$ K between metal walls in electric field applied in the perpendicular direction. We use analytic theory and molecular dynamics simulation. In this situation, there is a global electrostatic coupling between the surface charges on the walls and the polarization in the bulk. Then, the correlation function of the polarization density $p_z(r)$ along the applied field contains a homogeneous part inversely proportional to the cell volume $V$. Accounting for the long-range dipolar interaction, we derive the Kirkwood-Fr$\ddot{\rm{o}}$hlich formula for the polarization fluctuations when the specimen volume $v$ is much smaller than $V$. However, for not small $v/V$, the homogeneous part comes into play in dielectric relations. We also calculate the distribution of $E_k^{\rm loc}$ in applied field. As a unique feature of water, its magnitude $|E_k^{\rm loc}|$ obeys a Gaussian distribution with a large mean value $E_0 \cong 17~$V$/$nm, which arises mainly from the surrounding hydrogen-bonded molecules. Since $|μ_k|E_0\sim 30 k_{\rm B}T$, $μ_k$ becomes mostly parallel to $E_k^{\rm loc}$. As a result, the orientation distributions of these two vectors nearly coincide, assuming the classical exponential form. In dynamics, the component of $μ_k(t)$ parallel to $E_k^{\rm loc}(t)$ changes on the timescale of the hydrogen bonds $\sim 5$ ps, while its smaller perpendicular component undergoes librational motions on timescales of 0.01 ps.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Hydrodynamics in bridging and aggregation of two colloidal particles in a near-critical binary mixture

We investigate bridging and aggregation of two colloidal particles in a near-critical binary mixture when the fluid far from the particles is outside the coexistence (CX) curve and is rich in the component disfavored by the colloid surfaces. In such situations, the adsorption-induced interaction is enhanced, leading to bridging and aggregation of the particles. We realize bridging firstly by changing the temperature with a fixed interparticle separation and secondly by letting the two particles aggregate. The interparticle attractive force dramatically increases upon bridging. The dynamics is governed by hydrodynamic flow around the colloid surfaces. In aggregation, the adsorption layers move with the particles and squeezing occurs at narrow separation. We take into account the renormalization effect due to the critical fluctuations using the recent local functional theory [J. Chem. Phys. 136, 114704 (2012)].

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Molecular Dynamics Simulation of Water between Metal Walls under Electric Field: Dielectric Response and Dynamics after Field Reversal

We study water between parallel metal walls under applied electric field accounting for the image effect at $T=298$ K. The electric field due to the surface charges serves to attract and orient nearby water molecules, while it tends to a constant determined by the mean surface charge density away from the walls. We find Stern boundary layers with thickness about $5$ $\rm Å$ and a homogeneously polarized bulk region. The molecules in the layers more sensitively respond to the applied field than in the bulk. As a result, the potential drop in the layers is larger than that in the bulk unless the cell length exceeds 10 nm. We also examine the hydrogen bonds, which tend to make small angles with respect to the walls in the layers even without applied field. The average local field considerably deviates from the classical Lorentz field and the local field fluctuations are very large in the bulk. If we suppose a nanometer-size sphere around each molecule, the local field contribution from its exterior is nearly equal to that from the continuum electrostatics and that from its interior yields the deviation from the classical Lorentz field. As a nonequilibrium problem, we investigate the dynamics after a reversal of applied field, where the relaxation is mostly caused by large-angle rotational jumps after 1 ps due to the presence of the hydrogen bond network. The molecules undergoing these jumps themselves form hydrogen-bonded clusters heterogeneously distributed in space.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Applying electric field to charged and polar particles between metallic plates: Extension of the Ewald method

We develop an efficient Ewald method of molecular dynamics simulation for calculating the electrostatic interactions among charged and polar particles between parallel metallic plates, where we may apply an electric field with an arbitrary size. We use the fact that the potential from the surface charges is equivalent to the sum of those from image charges and dipoles located outside the cell. We present simulation results on boundary effects of charged and polar fluids, formation of ionic crystals, and formation of dipole chains, where the applied field and the image interaction are crucial. For polar fluids, we find a large deviation of the classical Lorentz-field relation between the local field and the applied field due to pair correlations along the applied field. As general aspects, we clarify the difference between the potential-fixed and the charge-fixed boundary conditions and examine the relationship between the discrete particle description and the continuum electrostatics.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Attractive interaction and bridging transition between neutral colloidal particles due to preferential adsorption in a near-critical binary mixture

We examine the solvent-mediated interaction between two neutral colloidal particles due to preferential adsorption in a near-critical binary mixture. We take into account the renormalization effect due to the critical fluctuations using the recent local functional theory $[$J. Chem. Phys. {\bf 136}, 114704 (2012)$]$. We calculate the free energy and the force between two colloidal particles as functions of the temperature $T$, the composition far from the colloidal particles $c_\infty$, and the colloid separation $\ell$. The interaction is much enhanced when the component favored by the colloid surfaces is poor in the reservoir. For such off-critical compositions, we find a surface of a first-order bridging transition $\ell= \ell_{\rm cx}(T,c_\infty)$ in the $T$-$c_\infty$-$\ell$ space in a universal, scaled form, across which a discontinuous change occurs between separated and bridged states. This surface starts from the bulk coexistence surface (CX) and ends at a bridging critical line $\ell= \ell_{c}(T)$. On approaching the critical line, the discontinuity vanishes and the derivatives of the force with respect to $T$ and $\ell$ both diverge. Furthermore, bridged states continuously change into separated states if $c_\infty$ (or $T$) is varied from a value on CX to value far from CX with $\ell$ kept smaller than $\ell_c(T)$.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Formation of double glass in binary mixtures of anisotropic particles

We study glass transitions in mixtures of elliptic and circular particles in two dimensions using an orientation-dependent Lennard-Jones potential. Changing anisotropic parameters of the potential, the size ratio, and the concentration, we realize double glass, where both the particle positions and orientations are disordered but still hold mesoscopic order. The ellipses are anchored around the circular impurities in the homeotropic or planar directions. We examine slowing-down of rotational and translational time-correlation functions. Turnover motions of the ellipses are activated more frequently than the configuration changes, where the latter cause the structural relaxation.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Formation of double glass in binary mixtures of anisotropic particles: Dynamic heterogeneities in rotations and displacements

We study glass behavior in a mixture of elliptic and circular particles in two dimensions at low temperatures using an orientation-dependent Lennard-Jones potential. The ellipses have a mild aspect ratio ($\sim 1.2$) and tend to align at low temperatures, while the circular particles play the role of impurities disturbing the ellipse orientations at a concentration of 20%. These impurities have a size smaller than that of the ellipses and attract them in the homeotropic alignment. As a result, the coordination number around each impurity is mostly five or four in glassy states. We realize double glass, where both the orientations and the positions are disordered but still hold mesoscopic order. We find a strong heterogeneity in the flip motions of the ellipses, which sensitively depends on the impurity clustering. The average rotational diffusion constant is determined by rapidly flipping ellipses, which still remain even at low temperatures in our model. In contrast, the non-flip rotations (with angle changes not close to $\pm π$) are mainly caused by the cooperative configuration changes involving many particles. Then, there arises a long-time heterogeneity in the non-flip rotations closely correlated with the dynamic heterogeneity in displacements.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Orientational glass in mixtures of elliptic and circular particles: Structural heterogeneities, rotational dynamics, and rheology

Using molecular dynamics simulation with an angle-dependent Lennard-Jones potential, we study orientational glass with quadrupolar symmetry in mixtures of elliptic particles and circular impurities in two dimensions. With a mild aspect ratio ($=1.2$) and a mild size ratio ($=1.2$), we realize a plastic crystal at relatively high temperature $T$. With further lowering $T$, we find a structural phase transition for very small impurity concentration $c$ and pinned disordered orientations for not small $c$. The ellipses are anchored by the impurities in the planar alignment. With increasing $c$, the orientation domains composed of isosceles triangles gradually become smaller, resulting in orientational glass with crystal order. In our simulation, the impurity distribution becomes heterogeneous during quenching from liquid, which then produces rotational dynamic heterogeneities. We also examine rheology in orientational glass to predict a shape memory effect and a superelasticity effect, where a large fraction of the strain is due to collective orientation changes.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Phase separation in a binary mixture confined between symmetric parallel plates: Capillary condensation transition near the bulk critical point

We investigate phase separation of near-critical binary mixtures between parallel symmetric walls in the strong adsorption regime. We take into account the renormalization effect due to the critical fluctuations using the recent local functional theory [J. Chem. Phys. 136, 114704 (2012)]. In statics, a van der Waals loop is obtained in the relation between the average order parameter $<psi>$ in the film and the chemical potential when the temperature $T$ is lower than the film critical temperature $T_c^{ca}$ (in the case of an upper critical solution temperature). In dynamics, we lower $T$ below the capillary condensation line from above $T_c^{ca}$. We calculate the subsequent time-development assuming no mass exchange between the film and the reservoir. In the early stage, the order parameter $psi$ changes only in the direction perpendicular to the walls. For sufficiently deep quenching, such one-dimensional profiles become unstable with respect to the fluctuations varying in the lateral directions. The late-stage coarsening is then accelerated by the hydrodynamic interaction. A pancake domain of the phase disfavored by the walls finally appears in the middle of the film.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Casimir amplitudes and capillary condensation of near-critical fluids between parallel plates: Renormalized local functional theory

We investigate the critical behavior of a near-critical fluid confined between two parallel plates in contact with a reservoir by calculating the order parameter profile and the Casimir amplitudes (for the force density and for the grand potential). Our results are applicable to one-component fluids and binary mixtures. We assume that the walls absorb one of the fluid components selectively for binary mixtures. We propose a renormalized local functional theory accounting for the fluctuation effects. Analysis is performed in the plane of the temperature T and the order parameter in the reservoir ψ_{\infty} . Our theory is universal if the physical quantities are scaled appropriately. If the component favored by the walls is slightly poor in the reservoir, there appears a line of first-order phase transition of capillary condensation outside the bulk coexistence curve. The excess adsorp- tion changes discontinuously between condensed and noncondensed states at the transition. With increasing T, the transition line ends at a capillary critical point T = T_c^{ca} slightly lower than the bulk critical temperature T_c. The Casimir amplitudes are larger than their critical-point values by 10-100 times between the transition line and the bulk coexistence curve and slightly above the capillary critical point.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Molecular dynamics simulation of orientational glass formation in anisotropic particle systems in three dimensions

We propose a simple microscopic model of molecular dynamics simulation to study orientational glass in three dimensions. We present simulation results for mixtures of mildly anisotropic particles and spherical impurities. We realize fcc solids without orientational order in a rotator phase. As the temperature $T$ is lowered, the disordered matrix is gradually replaced by four kinds of orientationally ordered, rhombohedral domains. Two-phase coexistence is realized in a temperature window. The impurities serve to anchor the orientations of the surrounding anisotropic particles, resulting in finely divided domains or medium long-range orientational order. We examine the rotational dynamics of the molecular orientations which is slowed down at low $T$. We predict the shape memory effect under a stretching cycle due to inter-variant transformation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Relationship between bond-breakage correlations and four-point correlations in heterogeneous glassy dynamics: Configuration changes and vibration modes

We investigate the dynamic heterogeneities of glassy particle systems in the theoretical schemes of bond breakage and four-point correlation functions. In the bond-breakage scheme, we introduce the structure factor S_b(q,t) and the susceptibility chi_b(t) to detect the spatial correlations of configuration changes. Here chi_b(t) attains a maximum at t =t_b^max as a function of time t, where the fraction of the particles with broken bonds phi_b(t)$ is about 1/2. In the four-point scheme, treating the structure factor S_4(q,t) and the susceptibility chi_4(t), we detect superpositions of the heterogeneity of bond breakage and that of thermal low-frequency vibration modes. While the former grows slowly, the latter emerges quickly to exhibit complex space-time behavior. In two dimensions, the vibration modes extending over the system yield significant contributions to the four-point correlations, which depend on the system size logarithmically. A maximum of chi_4(t) is attained at t= t_4^max, where these two contributions become of the same order. As a result, t_4^max is considerably shorter than t_b^max.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Selective solvation in aqueous mixtures: Interface deformations and instability

We briefly review the effects of selective solvation of ions in aqueous mixtures, where the ion densities and the composition fluctuations are strongly coupled. We then examine the surface tension γof a liquid-liquid interface in the presence of ions. We show that γcan be decreased drastically due to the electrostatic and solvation interactions near the interface. We calculate how the free energy is changed due to small surface undulations in the presence of an electric double layer. A surface instability occurs for negative γ, which can easily be realized for antagonistic ion pairs near the solvent criticality. Three-dimensional simulation shows how the surface instability is induced.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Slow relaxations and stringlike jump motions in fragile glass-forming liquids: Breakdown of the Stokes-Einstein relation

We perform molecular dynamics simulation on a glass-forming liquid binary mixture with the soft-core potential in three dimensions. We investigate crossover of the configuration changes caused by stringlike jump motions. With lowering the temperature $T$, the motions of the particles composing strings become larger in sizes and displacements, while those of the particles surrounding strings become smaller. Then, the contribution of the latter to time-correlation functions tends to be long-lived as $T$ is lowered. As a result, the relaxation time $τ_α$ and the viscosity $η$ grow more steeply than the inverse diffusion constant $D^{-1}$ at low $T$, leading to breakdown of the Stokes-Einstein relation. At low $T$, the diffusion occurs as activation processes and may well be described by short-time analysis of rare jump motions with broken bonds and large displacements. Some characteristic features of the van Hove self-correlation function arise from escape jumps over high potential barriers.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Structural phase transition and orientation-strain glass formation in anisotropic particle systems with impurities in two dimensions

Using a modified Lennard-Jones model for elliptic particles and spherical impurities, we present results of molecular dynamics simulation in two dimensions. In one-component systems of elliptic particles, we find an orientation phase transition on a hexagonal lattice as the temperature $T$ is lowered. It is also a structural one because of spontaneous strain. At low $T$, there arise three martensitic variants due to the underlying lattice, leading to a shape memory effect without dislocation formation. Thermal hysteresis, a minimum of the shear modulus, and a maximum of the specific heat are also found with varying $T$. With increasing the composition $c$ of impurities, the three kinds of orientation domains are finely divided, yielding orientation-strain glass with mesoscopically ordered regions still surviving. If the impurities are large and repulsive, planar anchoring of the elliptic particles occurs around the impurity surfaces. If they are small and attractive, homeotropic anchoring occurs. Clustering of impurities is conspicuous. With increasing the anchoring power and/or the composition of the impurities, positional disorder can also be enhanced. We also investigate the rotational dynamics of the molecular orientations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Charged colloids in an aqueous mixture with a salt

We calculate the ion and composition distributions around colloid particles in an aqueous mixture, accounting for preferential adsorption, electrostatic interaction, selective solvation among ions and polar molecules, and composition-dependent ionization. On the colloid surface, we predict a precipitation transition induced by strong preference of hydrophilic ions to water and a prewetting transition between weak and strong adsorption and ionization. These transition lines extend far from the solvent coexistence curve in the plane of the interaction parameter χ(or the temperature) and the average solvent composition. The colloid interaction is drastically altered by these phase transitions on the surface. In particular, the interaction is much amplified upon bridging of wetting layers formed above the precipitation line. Such wetting layers can either completely or partially cover the colloid surface depending on the average solvent composition.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Phase Transitions in Soft Matter Induced by Selective Solvation

We review our recent studies on selective solvation effects in phase separation in polar binary mixtures with a small amount of solutes. Such hydrophilic or hydrophobic particles are preferentially attracted to one of the solvent components. We discuss the role of antagonistic salt composed of hydrophilic and hydrophobic ions, which undergo microphase separation at water-oil interfaces leading to mesophases. We then discuss phase separation induced by a strong selective solvent above a critical solute density np, which occurs far from the solvent coexistence curve. We also give theories of ionic surfactant systems and weakly ionized polyelectrolytes including solvation among charged particles and polar molecules. We point out that the Gibbs formula for the surface tension needs to include an electrostatic contribution in the presence of an electric double layer.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Pre-dewetting transition on a hydrophobic wall: Statics and dynamics

For one-component fluids, we predict a pre-dewetting phase transition between a thin and thick low-density layer in liquid on a wall repelling the fluid. This is the case of a hydrophobic wall for water. A pre-dewetting line starts from the coexistence curve and ends at a surface critical point in the phase diagram. We calculate this line numerically using the van der Waals model and analytically using the free energy expansion up to the quartic order. We also examine the pre-dewetting dynamics of a layer created on a hydrophobic spot on a heterogeneous wall. It is from a thin to thick layer during decompression and from a thick to thin layer during compression. Upon the transition, a liquid region above the film is cooled for decompression and heated for compression due to latent heat convection and a small pressure pulse is emitted from the film into the liquid.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Structural phase transition of anisotropic particles and formation of orientation-strain glass with addition of impurities

Using a modified Lennard-Jones model for anisotropic particles, we present results of molecular dynamics simulation in two dimensions. In one-component systems, we find crystallization, a Berezinskii-Kosterlitz-Thouless phase, and a structural phase transition, as the temperatures is lowered. In the lowest temperature range, the crystal is composed of three martensitic variants on a hexagonal lattice, exhibiting the shape memory effect. With addition of larger spherical particles (impurities), these domains are finely divided, yielding glass with slow time evolution. With increasing the impurity size, the structural or translational disorder is also proliferated.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Visualization of stacking faults in fcc crystals in plastic deformations

Using molecular dynamics simulation, we investigate the dynamics of stacking faults in fcc crystals in uniaxial stretching in a Lennard-Jones binary mixture composed of 4096 particles in three dimensions. We visualize stacking faults using a disorder variable $D_j(t)$ for each particle $j$ constructed from local bond order parameters based on spherical harmonics (Steinhardt order parameters). Also introducing a method of bond breakage, we examine how stacking faults are formed and removed by collective particle motions. These processes are relevant in plasticity of fcc crystals.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Jammed Particle Configurations and Dynamics in High-Density Lennard-Jones Binary Mixtures in Two Dimensions

We examine the changeover in the particle configurations and the dynamics in dense Lennard-Jones binary mixtures composed of small and large particles. By varying the composition at a low temperature, we realize crystal with defects, polycrystal with small grains, and glass with various degrees of disorder. In particular, we show configurations where small crystalline regions composed of the majority species are enclosed by percolated amorphous layers composed of the two species. We visualize the dynamics of configuration changes using the method of bond breakage and following the particle displacements. In quiescent jammed states, the dynamics is severely slowed down and is highly heterogeneous at any compositions. In shear, plastic deformations multiply occur in relatively fragile regions, growing into large-scale shear bands where the strain is highly localized. Such bands appear on short time scales and change on long time scales with finite life times.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Phase field model of solid-liquid and liquid-liquid phase transitions in flow and elastic fields in one-component systems

We construct a phase field model including hydrodynamics and elasticity in one-component systems. It can be used to investigate solid-liquid and liquid-liquid phase transitions. Upon first-order phase transition, a velocity field is induced around interfaces in the presence of a density difference between the two phases even without applied shear flow. As applications, we present simulation results on two cases of melting, where a solid domain is placed on a heated wall in one case and is suspended in a warmer liquid under shear flow in the other. We find that the solid domain moves or rotates as a whole due to elasticity, releasing latent heat. We also examine the liquid-liquid phase transition of a highly viscous domain into a less viscous liquid on a heated wall, where an inhomogeneous velocity field is induced within a projected part of the domain. In these phase transitions, the interface temperature is nearly equal to the coexisting temperature $T_{\rm cx}(p)$ away from the heated wall in the presence of heat flow in the surrounding liquid.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Plastic deformations in crystal, polycrystal, and glass in binary mixtures under shear: Collective yielding

Using molecular dynamics simulation, we examine the dynamics of crystal, polycrystal, and glass in a Lennard-Jones binary mixture composed of small and large particles in two dimensions. The crossovers occur among these states as the composition c is varied at fixed size ratio. Shear is applied to a system of 9000 particles in contact with moving boundary layers composed of 1800 particles. The particle configurations are visualized with a sixfold orientation angle alpha_j(t) and a disorder variable D_j(t) defined for particle j, where the latter represents the deviation from hexagonal order. Fundamental plastic elements are classified into dislocation gliding and grain boundary sliding. At any c, large-scale yielding events occur on the acoustic time scale. Moreover, they multiply occur in narrow fragile areas, forming shear bands. The dynamics of plastic flow is highly hierarchical with a wide range of time scales for slow shearing. We also clarify the relationship between the shear stress averaged in the bulk region and the wall stress applied at the boundaries.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Precipitation in aqueous mixtures with addition of strongly hydrophilic or hydrophobic solute

We examine phase separation in aqueous mixtures due to preferential solvation with a low-density solute (hydrophilic ions or hydrophobic particles). For hydrophilic ions, preferential solvation can stabilize water domains enriched with ions. This precipitation occurs in wide ranges of the temperature and the average composition above a critical solute density $n_p$, where the mixture solvent would be in a one-phase state without solute. The volume fraction of precipitated domains tends to zero as the average solute density $\bar n$ is decreased to $n_p$ or as the interaction parameter $χ$ is decreased to a critical value $χ_p$. If we start with one-phase states with ${\bar n}>n_p$ or $χ>χ_p$, precipitation proceeds via homogeneous nucleation or via heterogeneous nucleation, for example, around suspended colloids. In the latter case, colliod particles are wrapped by thick wetting layers. We also predict a first-order prewetting transition for $\bar n$ or $χ$ slightly below $n_p$ or $χ_p$.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Precipitation of water from aqueous mixtures with addition of hydrophilic ions

We examine phase separation in aqueous mixtures at fixed amounts of hydrophilic monovalent ions. When water is the minority component, preferential solvation can stabilize water domains enriched with ions. This ion-induced precipitation occurs in wide ranges of the temperature and the average composition where the solvent would be in one-phase states without ions. The volume fraction of such water domains is decreased to zero as the interaction parameter $χ$ (dependent on the temperature) is decreased toward a critical value for each average composition.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Self-organization in He4 near the superfluid transition in heat flow and gravity

We investigate the nonlinear dynamics of He4 slightly below the superffluid transition by integrating model F equations in three dimensions. When heated from above under gravity, a vortex tangle and a sheetlike phase slip are generated near the bottom plate. Then a self-organized superfluid containing high-density vortices and phase slips grows upward into an ordinary superfluid. The thermal resistance due to these defects produces a constant temperature gradient equal to the gradient of the pressure-dependent transition temperature $T_λ(p)$. In this self-organized region, the temperature deviation $T-T_λ(p)$ consists of a negative constant independent of the height and time-dependent fluctuations. Its time-average is calculated in good agreement with the experimental value (W.A. Moeur {\it et al.}, Phys. Rev. Lett. 78, 2421 (1997)).

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Solvation Effects in Phase Transitions in Soft Matter

Phase transitions in polar binary mixtures can be drastically altered even by a small amount of salt. This is because the preferential solvation strongly depends on the ambient composition. Together with a summary of our research in this problem, we present some detailed results on the role of antagonistic salt composed of hydrophilic and hydrophobic ions. These ions tend to segregate at liquid-liquid interfaces and selectively couple to water-rich and oil-rich composition fluctuations, leading to mesophase formation. In our two-dimensional simulation, the corasening of the domain structures can be stopped or slowed down, depending on the interaction parameter (or the temperature) and the salt density. We realize stripe patterns at the critical composition and droplet patterns at off-critical compositions. In the latter case, charged droplets emerge with considerable size dispersity in a percolated region. We also give the structure factors among the ions, accounting for the Coulomb interaction and the solvation interaction mediated by the composition fluctuations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Spreading with evaporation and condensation in one-component fluids

We investigate the dynamics of spreading of a small liquid droplet in gas in a one-component simple fluid, where the temperature is inhomogeneous around 0.9Tc and latent heat is released or generated at the interface upon evaporation or condensation (with Tc being the critical temperature). In the scheme of the dynamic van der Waals theory, the hydrodynamic equations containing the gradient stress are solved in the axisymmetric geometry. We assume that the substrate has a finite thickness and its temperature obeys the thermal diffusion equation. A precursor film then spreads ahead of the bulk droplet itself in the complete wetting condition. Cooling the substrate enhances condensation of gas onto the advancing film, which mostly takes place near the film edge and can be the dominant mechanism of the film growth in a late stage. The generated latent heat produces a temperature peak or a hot spot in the gas region near the film edge. On the other hand, heating the substrate induces evaporation all over the interface. For weak heating, a steady-state circular thin film can be formed on the substrate. For stronger heating, evaporation dominates over condensation, leading to eventual disappearance of the liquid region.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Bubble and droplet motion in binary mixtures: Evaporation-condensation mechanism and Marangoni effect

Bubble and droplet motion in binary mixtures is studied in weak heat and diffusion fluxes and in gravity by solving the linearized hydrodynamic equations supplemented with appropriate surface boundary conditions. Without gravity, the velocity field is induced by evaporation and condensation at the interface and by the Marangoni effect due to a surface tension gradient. In pure fluids, the latter nearly vanishes since the interface temperature tends to the coexistence temperature $T_{\rm cx}(p)$ even in heat flow. In binary mixtures, the velocity field can be much enhanced by the Marangoni effect above a crossover concentration $c^*$ inversely proportional to the radius $R$ of the bubble or droplet. Here $c^*$ is usually very small for large $R$ for non-azeotropic mixtures. The temperature and concentration deviations are also calculated.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Defect structures in nematic liquid crystals around charged particles

We numerically study the orientation deformations in nematic liquid crystals around charged particles. We set up a Ginzburg-Landau theory with inhomogeneous electric field. If the dielectric anisotropy varepsilon_1 is positive, Saturn ring defects are formed around the particles. For varepsilon_1<0, novel "ansa" defects appear, which are disclination lines with their ends on the particle surface. We find unique defect structures around two charged particles. To lower the free energy, oppositely charged particle pairs tend to be aligned in the parallel direction for varepsilon_1>0 and in the perpendicular plane for varepsilon_1<0 with respect to the background director . For identically charged pairs the preferred directions for varepsilon_1>0 and varepsilon_1<0 are exchanged. We also examie competition between the charge-induced anchoring and the short-range anchoring. If the short-range anchoring is sufficiently strong, it can be effective in the vicinity of the surface, while the director orientation is governed by the long-range electrostatic interaction far from the surface.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Dynamics of Binary Mixtures with Ions: Dynamic Structure Factor and Mesophase Formation

Dynamic equations are presented for polar binary mixtures containing ions in the presence of the preferential solvation. In one-phase states, we calculate the dynamic structure factor of the composition accounting for the ion motions. Microphase separation can take place for sufficiently large solvation asymmetry of the cations and the anions. We show two-dimensional simulation results of the mesophase formation with an antagonistic salt, where the cations are hydrophilic and the anions are hydrophobic. The structure factor S(q) in the resultant mesophase has a sharp peak at an intermediate wave number on the order of the Debye-Huckel wave number. As the quench depth is increased, the surface tension nearly vanishes in mesophases due to an electric double layer.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Ion distribution around a charged rod in one and two component solvents: Preferential solvation and first order ionization phase transition

In one and two component solvents, we calculate the counterion distribution around a charged rod treating the degree of ionization $α$ as an annealed variable dependent on its local environment. In the one component case, $α$ is determined under various conditions without and with salt. In the two component case, we take into account the preferential solvation of the counterions and the ionized monomers and the short-range interaction between the rod and the solvent without salt. It then follows a composition-dependent mass action law. Mesoscopic variations of the composition and the counterions are produced around a chraged rod, which sensitively depend on various parameters of the molecular interactions. Furthermore, we predict a first order phase transition of weak-to-strong dissociation for strong preferential solvation. It can occur in expanded states of a polymer chain. This transition line starts from a point on the solvent coexistence curve and ends at a critical point in the plane of the temperature and the solvent composition. The composition change around a charged rod is long-ranged near the solvent critical point.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Onion structures induced by hydrophilic and hydrophobic ions in a binary mixture

Phase transition is observed between one-phase disordered phase and an ordered phase with multi-lamellar (onion) structures in an off-critical mixture of D$_2$O and 3-methylpyridine (3MP) containing a salt at 85mM. The salt consists of hydrophilic cations and hydrophobic anions, which interact asymmetrically with the solvent composition fluctuations inducing mesophases. The structure factor of the composition distribution obtained from small-angle neutron scattering has a peak at an intermediate wave number in the disordered phase and multiple peaks in the ordered phase. Lamellar layers forming onions are composed of solvation-induced charged membranes swollen by D$_2$O. The onion phase is realized only for small volume fractions of 3MP (in D$_2$O-rich solvent).

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Structural and dynamical heterogeneities in two-dimensional melting

Using molecular dynamics simulation, we study structural and dynamical heterogeneities at melting in two-dimensional one-component systems with 36000 particles. Between crystal and liquid we find intermediate hexatic states, where the density fluctuations are enhanced at small wave number k as well as those of the six-fold orientational order parameter. Their structure factors both grow up to the smallest wave number equal to the inverse system length. The intermediate scattering function of the density S(k,t) is found to relax exponentially with decay rate Gamma_k ~ k^z with z~2.6 at small k in the hexatic phase.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Dislocation Formation and Work-Hardening in Two-Phase Alloys

A phase field model is presented to investigate dislocation formation (coherency loss) and workhardening in two-phase binary alloys. In our model the elastic energy density is a periodic function of the shear and tetragonal strains, which allows multiple formation of slips (dislocation dipoles). By numerically integrating the dynamic equations in two dimensions, we find that dislocations appear in pairs in the interface region and grow into slips. One end of each slip glides preferentially into the softer region, while the other end remains trapped at the interface. Under uniaxial stretching at deep quenching, slips appear in the softer regions and do not penetrate into the harder domains, giving rise to an increase of the stress with increasing applied strain in plastic flow.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Plastic Flow in Two-Dimensional Solids

A time-dependent Ginzburg-Landau model of plastic deformation in two-dimensional solids is presented. The fundamental dynamic variables are the displacement field $\bi u$ and the lattice velocity ${\bi v}=\p {\bi u}/\p t$. Damping is assumed to arise from the shear viscosity in the momentum equation. The elastic energy density is a periodic function of the shear and tetragonal strains, which enables formation of slips at large strains. In this work we neglect defects such as vacancies, interstitials, or grain boundaries. The simplest slip consists of two edge dislocations with opposite Burgers vectors. The formation energy of a slip is minimized if its orientation is parallel or perpendicular to the flow in simple shear deformation and if it makes angles of $\pm π/4$ with respect to the stretched direction in uniaxial stretching. High-density dislocations produced in plastic flow do not disappear even if the flow is stopped. Thus large applied strains give rise to metastable, structurally disordered states. We divide the elastic energy into an elastic part due to affine deformation and a defect part. The latter represents degree of disorder and is nearly constant in plastic flow under cyclic straining.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

Kinetic Heterogeneities in a Highly Supercooled Liquid

We study a highly supercooled two-dimensional fluid mixture via molecular dynamics simulation. We follow bond breakage events among particle pairs, which occur on the scale of the $α$ relaxation time $τ_α$. Large scale heterogeneities analogous to the critical fluctuations in Ising systems are found in the spatial distribution of bonds which are broken in a time interval with a width of order $0.05τ_α$. The structure factor of the broken bond density is well approximated by the Ornstein-Zernike form. The correlation length is of order $100 σ_1$ at the lowest temperature studied, $σ_1$ being the particle size. The weakly bonded regions thus identified evolve in time with strong spatial correlations.

0 citations0 reviewsNo code linkedOpen access