The interplay between demixing and freezing in two-dimensional symmetrical binary mixtures of Lennard-Jones particles is studied using Monte Carlo simulation. It is demonstrated that different scenarios are possible. For example, the line of continuous liquid demixing transition can start at the liquid side of the vapor-liquid coexistence at the lower critical end point and then it can terminate at the liquid side of the liquid–demixed solid coexistence at the upper critical end point. Other situations are also possible. We distinguish four different scenarios depending on the interactions between unlike particles.

Recently, different approaches have been proposed for studying basic properties of time series from a complex network perspective. In this work, the corresponding potentials and limitations of networks based on recurrences in phase space are investigated in some detail. We discuss the main requirements that permit a feasible system-theoretic interpretation of network topology in terms of dynamically invariant phase-space properties. Possible artifacts induced by disregarding these requirements are pointed out and systematically studied. Finally, a rigorous interpretation of the clustering coefficient and the betweenness centrality in terms of invariant objects is proposed.

It is shown that the lattice Boltzmann equation (LBE) corresponds to an explicit Verlet time-marching scheme for a continuum generalized Boltzmann equation with a memory delay equal to a half time step. This proves second-order accuracy of LBE with respect to this generalized equation, with no need of resorting to any implicit time-marching procedure (Crank-Nicholson) and associated nonlinear variable transformations. It is also shown, and numerically demonstrated, that this equivalence is not only formal, but it also translates into a complete equivalence of the corresponding computational schemes with respect to the hydrodynamic equations. Second-order accuracy with respect to the continuum kinetic equation is also numerically demonstrated for the case of the Taylor-Green vortex. It is pointed out that the equivalence is however broken for the case in which mass and/or momentum are not conserved, such as for chemically reactive flows and mixtures. For such flows, the time-centered implicit formulation may indeed offer a better numerical accuracy.

It is well known that the constant injection rate flow in radial Hele-Shaw cells leads to the formation of highly branched patterns, where finger tip-splitting events are plentiful. Different kinds of patterns arise in the lifting Hele-Shaw flow problem, where the cell’s gap width grows linearly with time. In this case, the morphology of the emerging structures is characterized by the strong competition among inward moving fingers. By employing a mode-coupling theory we find that both finger tip-splitting and finger competition can be restrained by properly adjusting the injection rate and the time-dependent gap width, respectively. Our theoretical model approaches the problem analytically and is capable of capturing these important controlling mechanisms already at weakly nonlinear stages of the dynamics.

We investigate the particle and kinetic-energy densities for a system of $N$ fermions bound in a local (mean-field) potential $V\left(\mathbf{r}\right)$. We generalize a recently developed semiclassical theory [J. Roccia and M. Brack, Phys. Rev. Lett. 100, 200408 (2008)] in which the densities are calculated in terms of the closed orbits of the corresponding classical system to $D>1$ dimensions. We regularize the semiclassical results (i) for the $U\left(1\right)$ symmetry breaking occurring for spherical systems at $r=0$ and (ii) near the classical turning points where the Friedel oscillations are predominant and well reproduced by the shortest orbit going from $r$ to the closest turning point and back. For systems with spherical symmetry, we show that there exist two types of oscillations which can be attributed to radial and nonradial orbits, respectively. The semiclassical theory is tested against exact quantum-mechanical calculations for a variety of model potentials. We find a very good overall numerical agreement between semiclassical and exact numerical densities even for moderate particle numbers $N$. Using a “local virial theorem,” shown to be valid (except for a small region around the classical turning points) for arbitrary local potentials, we can prove that the Thomas-Fermi functional ${\tau }_{\text{TF}}\left[\rho \right]$ reproduces the oscillations in the quantum-mechanical densities to first order in the oscillating parts.

The generalized Langevin theory for a solution has been derived as the infinite dilution limit of the theory for a two component mixture. Following a similar formalism, the mode coupling approximations of the memory kernel have been also obtained. We have applied this method for one component bulk liquid of Lennard-Jones spheres and proved this approximation theoretically. The analysis of the space and time pair correlation proposed by Van Hove has been carried out as a function of solute particle sizes. It is found that the size of the solute particle is deeply related to the relaxation process of the solvation structure characterized around a solute particle at equilibrium. We have also investigated the relation between the different thermodynamic environments and relaxation process. From these studies, we have obtained the useful information about the rapidity of the relaxation of the solvation structure around a solute at equilibrium.

We use a quantitative phase-field approach to study directional solidification in various three-dimensional geometries for realistic parameters of a transparent binary alloy. The geometries are designed to study the steady-state growth of spatially extended hexagonal arrays, linear arrays in thin samples, and axisymmetric shapes constrained in a tube. As a basis to address issues of dynamical pattern selection, the phase-field simulations are specifically geared to identify ranges of primary spacings for the formation of the classically observed “fingers” (deep cells) with blunt tips and “needles” with parabolic tips. Three distinct growth regimes are identified that include a low-velocity regime with only fingers forming, a second intermediate-velocity regime characterized by coexistence of fingers and needles that exist on separate branches of steady-state growth solutions for small and large spacings, respectively, and a third high-velocity regime where those two branches merge into a single one. Along the latter, the growth shape changes continuously from fingerlike to needlelike with increasing spacing. These regimes are strongly influenced by crystalline anisotropy with the third regime extending to lower velocity for larger anisotropy. Remarkably, however, steady-state shapes and tip undercoolings are only weakly dependent on the growth geometry. Those results are used to test existing theories of directional finger growth as well as to interpret the hysteretic nature of the cell-to-dendrite transition.

Compositional lipid domains (lipid rafts) in plasma membranes are believed to be important components of many cellular processes. The mechanisms by which cells regulate the sizes, lifetimes, and spatial localization of these domains are rather poorly understood at the moment. We propose a robust mechanism for the formation of finite-sized lipid raft domains in plasma membranes, the competition between phase separation in an immiscible lipid system and active cellular lipid transport processes naturally leads to the formation of such domains. Simulations of a continuum model reveal that the raft size distribution is broad and the average raft size is strongly dependent on the rates of cellular and interlayer lipid transport processes. We demonstrate that spatiotemporal variations in the recycling may enable the cell to localize larger raft aggregates at specific parts along the membrane. Moreover, we show that membrane compartmentalization may further facilitate spatial localization of the raft domains. Finally, we demonstrate that local interactions with immobile membrane proteins can spatially localize the rafts and lead to further clustering.

Plasmids are extra chromosomal DNA that can confer to their hosts’ supplementary characteristics such as antibiotic resistance. Plasmids code for their copy number through their own replication frequency. Even though the biochemical networks underlying the plasmid copy number (PCN) regulation processes have been studied and modeled, no measurement of the heterogeneity in PCN within a whole population has been done. We have developed a fluorescent-based measurement system, which enables determination of the mean and noise in PCN within a monoclonal population of bacteria. Two different fluorescent protein reporters were inserted: one on the chromosome and the other on the plasmid. The fluorescence of these bacteria was measured with a microfluidic flow cytometry device. We show that our measurements are consistent with known plasmid characteristics. We find that the partitioning system lowers the PCN mean and standard deviation. Finally, bacterial populations were allowed to grow without selective pressure. In this case, we were able to determine the plasmid loss rate and growth inhibition effect.

We study regimes of strong and weak generalized synchronization in chaotically forced nonlinear flows. The transition between these dynamical states can occur via a number of different routes, and here we examine the onset of weak generalized synchrony through intermittency and blowout bifurcations. The quantitative characterization of this dynamical transition is facilitated by measures that have been developed for the study of strange nonchaotic motion. Weak and strong generalized synchronous motion show contrasting sensitivity to parametric variation and have distinct distributions of finite-time Lyapunov exponents.

The generation of arbitrarily shaped nonspherical laser-induced cavitation bubbles is demonstrated with a optical technique. The nonspherical bubbles are formed using laser intensity patterns shaped by a spatial light modulator using linear absorption inside a liquid gap with a thickness of $40\mu \text{m}$. In particular we demonstrate the dynamics of elliptic, toroidal, square, and V-shaped bubbles. The bubble dynamics is recorded with a high-speed camera at framing rates of up to $300000$ frames per second. The observed bubble evolution is compared to predictions from an axisymmetric boundary element simulation which provides good qualitative agreement. Interesting dynamic features that are observed in both the experiment and simulation include the inversion of the major and minor axis for elliptical bubbles, the rotation of the shape for square bubbles, and the formation of a unidirectional jet for V-shaped bubbles. Further we demonstrate that specific bubble shapes can either be formed directly through the intensity distribution of a single laser focus, or indirectly using secondary bubbles that either confine the central bubble or coalesce with the main bubble. The former approach provides the ability to generate in principle any complex bubble geometry.

We present a few-mode Galerkin model for convection in binary fluid layers subject to an approximation to realistic horizontal boundary conditions at positive separation ratios. The model exhibits convection patterns in form of rolls and squares. The stable squares at onset develop into stable rolls at higher thermal driving. In between, a regime of a so-called cross roll structure is found. The results of our few-mode model are in good agreement with both experiments and numerical multimode simulations.

We analyze the isotropic component of turbulent flows spanning a broad range or Reynolds numbers. The aim is to identify scaling laws and their Reynolds number dependence in flows under different mechanical forcings. To this end, we applied an SO(3) decomposition to data stemming from direct numerical simulations with spatial resolutions ranging from ${64}^3$ to ${1024}^3$ grid points, and studied the scaling of high order moments of the velocity field. The study was carried out for two different flows obtained forcing the system with a Taylor-Green vortex or the Arn’old-Beltrami-Childress flow. Our results indicate that helicity has no significant impact on the scaling exponents as obtained from the generalized structure functions. Intermittency effects increase with the Reynolds number in the range of parameters studied, and in some cases are larger than what can be expected from several models of intermittency in the literature. The observed dependence of intermittency with the Reynolds number decreases if extended self-similarity is used to estimate the exponents.

A self-consistent kinetic model of a low-pressure dc glow discharge with dust particles based on Boltzmann equation for the electron energy distribution function is presented. The ions and electrons production in ionizing processes as well as their recombination on the dust-particle surface and on the discharge tube wall were taken into account. The influence of dust-particle concentration $N_d$ on gas discharge and dust particles parameters was investigated. It is shown that the increase of $N_d$ leads to the increase of an averaged electric field and ion density, and to the decrease of a dust-particle charge and electron density in the dusty cloud. The results were obtained in a wide region of different discharge and dusty plasma parameters: dust particles density ${10}^2–{10}^8{\text{cm}}^{-3}$, discharge current density ${10}^{-1}–{10}^1\text{mA}/{\text{cm}}^2$, and dust particles radius 1, 2, and $5\mu \text{m}$. The scaling laws for dust-particle surface potential and electric filed dependencies on dust-particle density, particle radius and discharge currents were revealed. It is shown that the absorption of electrons and ions on the dust particles surface does not lead to the electron energy distribution function depletion due to a self-consistent adjustment of dust particles and discharge parameters.

Through an appropriate change in variables, we find that the three-dimensional acoustic wave equation is subject to the transformation media interpretation. In particular, we determine that this interpretation can be extended beyond the pressure difference to also account for the momentum transported by the wave. The suitability of momentum transport is especially interesting as it is an example where the field of interest is not governed by a wave equation. We examine how both fields behave in the case of cloaking. Explicit consideration of the boundary conditions shows that perfect cloaking is preserved, even when the incoming momentum is nonzero at the surface of the cloak.

We present a computational scheme based on classical molecular dynamics to study chaotic billiards in static external magnetic fields. The method allows us to treat arbitrary geometries and several interacting particles. We test the scheme for rectangular single-particle billiards in magnetic fields and find a sequence of regularity islands at integer aspect ratios. In the case of two Coulomb-interacting particles the dynamics is dominated by chaotic behavior. However, signatures of quasiperiodicity can be identified at weak interactions, as well as regular trajectories at strong magnetic fields. Our scheme provides a promising tool to monitor the classical limit of many-electron semiconductor nanostructures and transport systems up to high magnetic fields.

We present two numerical schemes of high accuracy for one-dimensional dynamical systems. They are modifications of the discrete gradient method and keep its advantages, including stability and conservation of the energy integral. However, their accuracy is higher by several orders of magnitude.

An all-chemical analog of clock-wave-front model for somitogenesis is proposed. The spatial periodicity can be obtained by arresting the homogeneous oscillations in a typical two-component reaction-diffusion system in the Hopf region by interacting with a chemical wave front. The patterns can be controlled by tuning the wave speed of the front.

We study, using exact numerical simulations, the statistics of the longest excursion $l_{\text{max}}\left(t\right)$ up to time $t$ for the fractional Brownian motion with Hurst exponent $0<H<1$. We show that in the large $t$ limit, $⟨l_{\text{max}}\left(t\right)⟩\propto Q_{\infty }t$, where $Q_{\infty }\equiv Q_{\infty }\left(H\right)$ depends continuously on $H$. These results are compared with exact analytical results for a renewal process with an associated persistence exponent $\theta =1-H$. This comparison shows that $Q_{\infty }\left(H\right)$ carries the clear signature of non-Markovian effects for $H\ne 1/2$. The preasymptotic behavior of $⟨l_{\text{max}}\left(t\right)⟩$ is also discussed.

The non-Markovian diffusion of dispersed particles in a semi-infinite cell of an isotropic fluid limited by an adsorbing-desorbing surface is theoretically investigated. The density of dispersed particles in the bulk is a time dependent function and the time dependent density of surface particles is governed by a modified kinetic equation with a time dependent kernel. In this framework, the densities of bulk and surface particles are analytically determined, taking into account the conservation of the number of particles immersed in the sample. This system exhibits anomalous diffusion behavior as well as memory effects in the adsorption-desorption process. The results obtained here are expected to be useful to investigate the adsorption-desorption phenomena of neutral as well as charged particles in an isotropic fluid in contact with a solid substrate when the anomalous diffusion is present.