It has been shown that a nanoliter chamber separated by a wall of asymmetric obstacles can lead to an inhomogeneous distribution of self-propelled microorganisms. Although it is well established that this rectification effect arises from the interaction between the swimmers and the noncentrosymmetric pillars, here we demonstrate numerically that its efficiency is strongly dependent on the detailed dynamics of the individual microorganism. In particular, for the case of run-and-tumble dynamics, the distribution of run lengths, the rotational diffusion, and the partial preservation of run orientation memory through a tumble are important factors when computing the rectification efficiency. In addition, we optimize the geometrical dimensions of the asymmetric pillars in order to maximize the swimmer concentration and we illustrate how it can be used for sorting by swimming strategy in a long array of parallel obstacles.

Self-organization of hydrogen bubbles generated by laser-treated areas of an aluminum plate etched in a basic aqueous solution of ammonia is studied experimentally and theoretically. The dynamics of the establishment of a stationary pattern of gas bubbles is experimentally is shown. In the theoretical model, the velocity field of liquid flows around an ensemble of several bubbles is obtained. Modeling of the process of self-organization of gas bubbles is performed on the basis of a continuum model of a bubble jet. Under certain assumptions, the pressure of a diluted system of bubbles is described by an equation similar to that for nonideal gas, which follows the van der Waals equation of state. The model predicts an alignment of gas bubbles along bisectors of the laser-treated area limited by a square, which is in good agreement with experimental observations. Further development of the model leads to an equation with a negative diffusion coefficient that may be responsible for symmetry breakdown and pattern formation.

We report large-scale computer simulations of the hard-disk system at high densities in the region of the melting transition. Our simulations reproduce the equation of state, previously obtained using the event-chain Monte Carlo algorithm, with a massively parallel implementation of the local Monte Carlo method and with event-driven molecular dynamics. We analyze the relative performance of these simulation methods to sample configuration space and approach equilibrium. Our results confirm the first-order nature of the melting phase transition in hard disks. Phase coexistence is visualized for individual configurations via the orientational order parameter field. The analysis of positional order confirms the existence of the hexatic phase.

The quantum entropy is usually defined using von Neumann's formula, which measures lack of information and vanishes for pure states. In contrast, we obtain a formula for the entropy of a pure state as it is measured from thermodynamic experiments, solely from the self-entanglement of the wave function, and find strong numerical evidence that the two are in agreement for nonintegrable systems, both for energy eigenstates and for states that are obtained at long times under the evolution of more general initial conditions. This is an extension of Boltzmann's hypothesis for classical systems, relating microscopic motion to thermodynamics.

We consider the area coverage of radial Lévy flights in a finite square area with periodic boundary conditions. From simulations we show how the fractal path dimension $d_f$ and thus the degree of area coverage depends on the number of steps of the trajectory, the size of the area, and the resolution of the applied box counting algorithm. For sufficiently long trajectories and not too high resolution, the fractal dimension returned by the box counting method equals two, and in that sense the Lévy flight fully covers the area. Otherwise, the determined fractal dimension equals the stable index of the distribution of jump lengths of the Lévy flight. We provide mathematical expressions for the turnover between these two scaling regimes. As complementary methods to analyze confined Lévy flights we investigate fractional order moments of the position for which we also provide scaling arguments. Finally, we study the time evolution of the probability density function and the first passage time density of Lévy flights in a square area. Our findings are of interest for a general understanding of Lévy flights as well as for the analysis of recorded trajectories of animals searching for food or for human motion patterns.

By using an achiral monomer and a photoinitiator, we introduced polymer networks into the Sm-$C^{*}$ phase of an antiferroelectric liquid crystal that forms ferro-, ferri-, and antiferroelectric phases. We then investigated the temperature dependence of Bragg wavelengths selectively reflected from these samples and found that the reflection bands shift to shorter wavelengths with increasing polymer concentration. The intermediate and the Sm-$C^{*}$${}_A$ phases dominate over the Sm-$C^{*}$ phase for polymer concentration ⩾4$\%$ by weight and the Sm-$C^{*}$ phase disappears completely for a 6$\%$ polymer sample.

The nature of coil-globule transition and scaling behavior of a strongly charged polyelectrolyte chain in a solution system with explicit neutralizing counterions and solvent molecules are studied using replica-exchange Monte Carlo simulations, focusing on the effects of finite chain length. The results reveal that at the thermodynamic limit of infinitely long chain length, the coil-globule transition may remain first order. Phase transition temperatures at various ion concentrations are obtained by extrapolating the values obtained at finite chain lengths. Furthermore, it is found that the exponent $\nu$ of the radius of gyration, $\langle R_g$${}^2\rangle$ ∼ $N^{2\nu }$, can be slightly larger than 1 under some conditions.

The development of new models of biological tissues that consider cells in a discrete manner is becoming increasingly popular as an alternative to continuum methods based on partial differential equations, although formal relationships between the discrete and continuum frameworks remain to be established. For crystal mechanics, the discrete-to-continuum bridge is often made by assuming that local atom displacements can be mapped homogeneously from the mesoscale deformation gradient, an assumption known as the Cauchy-Born rule (CBR). Although the CBR does not hold exactly for noncrystalline materials, it may still be used as a first-order approximation for analytic calculations of effective stresses or strain energies. In this work, our goal is to investigate numerically the applicability of the CBR to two-dimensional cellular-scale models by assessing the mechanical behavior of model biological tissues, including crystalline (honeycomb) and noncrystalline reference states. The numerical procedure involves applying an affine deformation to the boundary cells and computing the quasistatic position of internal cells. The position of internal cells is then compared with the prediction of the CBR and an average deviation is calculated in the strain domain. For center-based cell models, we show that the CBR holds exactly when the deformation gradient is relatively small and the reference stress-free configuration is defined by a honeycomb lattice. We show further that the CBR may be used approximately when the reference state is perturbed from the honeycomb configuration. By contrast, for vertex-based cell models, a similar analysis reveals that the CBR does not provide a good representation of the tissue mechanics, even when the reference configuration is defined by a honeycomb lattice. The paper concludes with a discussion of the implications of these results for concurrent discrete and continuous modeling, adaptation of atom-to-continuum techniques to biological tissues, and model classification.

This paper deals with the theoretical modeling of a rotating solitary surface wave that was observed during water drainage from a cylindrical reservoir, when shallow water conditions were reached. It represents an improvement of our previous study, where the radial flow perturbation was neglected. This assumption led to the classical planar Korteweg–de Vries equation for the wall wave profile, which did not account for the rotational character of the base flow. The present formulation is based on a less restricting condition and consequently corrects the last shortcoming. Now the influence of the background flow appears in the wave characteristics. The theory provides a better physical depiction of the unique experiment by predicting fairly well the wave profile at least in the first half of its lifetime and estimating the speed of the observed wave with good accuracy.

We present two-dimensional numerical investigations of the temperature and velocity evolution of a pure near-critical fluid confined in microchannels. The fluid is subjected to two sides heating after it reached isothermal steady state. We focus on the abnormal behaviors of the near-critical fluid in response to the sudden imposed heat flux. New thermal-mechanical effects dominated by fluid instability originating from the boundary and local equilibrium process are reported. Near the microchannel boundaries, the instability grows very quickly and an unexpected vortex formation mode is identified when near-critical thermal-mechanical effect is interacting with the microchannel shear flow. The mechanism of the new kind of Kelvin-Helmholtz instability induced by boundary expansion and density stratification processes is also discussed in detail. This mechanism may bring about innovations in the field of microengineering.

Interactions between nearly bicritical modes in Taylor-Couette flow, which have been concerned with the framework of weakly nonlinear theory, are extended to fully nonlinear Navier-Stokes computation. For this purpose, a standard Newton solver for axially periodic flows is generalized to compute any mixed solutions having up to two phases, which typically arise from interactions of two spiral or Taylor vortex modes. Also, a simple theory is developed in order to classify the mixed solutions. With these methods, we elucidate pattern formation phenomena, which have been observed in a Taylor-Couette flow experiment. Focusing on the counter-rotating parameter range, all possible classes of interaction of various solutions with different azimuthal and axial wave numbers are considered within our computational restriction, and we observe numerous connection branches, e.g., footbridge solutions. Some of the mixed solutions result in a three-dimensional wavy spiral solution with axial relative periodicity or an axially doubly periodic toroidally closed vortex solution. The possible connection of the former solution family to spiral turbulence, which has been observed in highly counter-rotating Taylor-Couette flow, is discussed.

The theoretical model by Sorasio et al. [Phys. Rev. Lett. 96, 045005 (2006)] for the steady state Mach number of electrostatic shocks formed in the interaction of two plasma slabs of arbitrary density and temperature is generalized for relativistic electron and nonrelativistic ion temperatures. We find that the relativistic correction leads to lower Mach numbers and as a consequence ions are reflected with lower energies. The steady state bulk velocity of the downstream population is introduced as an additional parameter to describe the transition between the minimum and maximum Mach numbers that is dependent on the initial density and temperature ratios. In order to transform the solitonlike solution in the upstream region into a shock, a population of reflected ions is considered and differences from a zero-ion temperature model are discussed.

We combine the finite-size scaling method with the mesh-free spectral method to calculate quantum critical parameters for a given Hamiltonian. The basic idea is to expand the exact wave function in a finite exponential basis set and extrapolate the information about system criticality from a finite basis to the infinite basis set limit. The used exponential basis set, though chosen intuitively, allows handling a very wide range of exponential decay rates and calculating multiple eigenvalues simultaneously. As a benchmark system to illustrate the combined approach, we choose the Hulthen potential. The results show that the method is very accurate and converges faster when compared with other basis functions. The approach is general and can be extended to examine near-threshold phenomena for atomic and molecular systems based on even-tempered exponential and Gaussian basis functions.

We propose an expanded spin-glass model, called the quantum Ghatak-Sherrington model, which considers spin-1 quantum spin operators in a crystal field and in a transverse field. The analytic solutions and phase diagrams of this model are obtained by using the one-step replica symmetry-breaking ansatz under the static approximation. Our results represent the splitting within one spin-glass (SG) phase depending on the values of crystal and transverse fields. The two separated SG phases, characterized by a density of filled states, show certain differences in their shapes and phase boundaries. Such SG splitting becomes more distinctive when the degeneracy of the empty states of spins is larger than one of their filled states.

We obtain exact analytic expressions for a class of functions expressed as integrals over the Haar measure of the unitary group in $d$ dimensions. Based on these general mathematical results, we investigate generic dynamical properties of complex open quantum systems, employing arguments from ensemble theory. We further generalize these results to arbitrary eigenvalue distributions, allowing a detailed comparison of typical regular and chaotic systems with the help of concepts from random matrix theory. To illustrate the physical relevance and the general applicability of our results we present a series of examples related to the fields of open quantum systems and nonequilibrium quantum thermodynamics. These include the effect of initial correlations, the average quantum dynamical maps, the generic dynamics of system-environment pure state entanglement and, finally, the equilibration of generic open and closed quantum systems.

Mechanically nonlinear energy harvesters driven by broadband vibrations modeled as white noise are investigated. We derive an upper bound on output power versus load resistance and show that, subject to mild restrictions that we make precise, the upper-bound performance can be obtained by a linear harvester with appropriate stiffness. Despite this, nonlinear harvesters can have implementation-related advantages. Based on the Kramers equation, we numerically obtain the output power at weak coupling for a selection of phenomenological elastic potentials and discuss their merits.

We study and solve some variations of the random $K$-satisfiability ($K$-SAT) problem—balanced $K$-SAT and biased random $K$-SAT—on a regular tree, using techniques we have developed earlier. In both these problems as well as variations of these that we have looked at, we find that the transition from the satisfiable to the unsatisfiable regime obtained on the Bethe lattice matches the exact threshold for the same model on a random graph for $K=2$ and is very close to the numerical value obtained for $K=3$. For higher $K$, it deviates from the numerical estimates of the solvability threshold on random graphs but is very close to the dynamical one-step-replica-symmetry-breaking threshold as obtained from the first nontrivial fixed point of the survey propagation algorithm.

An implementation of quantum absorption chillers with three qubits has been recently proposed that is ideally able to reach the Carnot performance regime. Here we study the working efficiency of such self-contained refrigerators, adopting a consistent treatment of dissipation effects. We demonstrate that the coefficient of performance at maximum cooling power is upper bounded by $3/4$ of the Carnot performance. The result is independent of the details of the system and the equilibrium temperatures of the external baths. We provide design prescriptions that saturate the bound in the limit of a large difference between the operating temperatures. Our study suggests that delocalized dissipation, which must be taken into account for a proper modeling of the machine-baths interaction, is a fundamental source of irreversibility which prevents the refrigerator from approaching the Carnot performance arbitrarily closely in practice. The potential role of quantum correlations in the operation of these machines is also investigated.

I extend a previous work to susceptible-infected-susceptible (SIS) models on weighted Barabási-Albert scale-free networks. Numerical evidence is provided that phases with slow, power-law dynamics emerge as the consequence of quenched disorder and tree topologies studied previously with the contact process. I compare simulation results with spectral analysis of the networks and show that the quenched mean-field (QMF) approximation provides a reliable, relatively fast method to explore activity clustering. This suggests that QMF can be used for describing rare-region effects due to network inhomogeneities. Finite-size study of the QMF shows the expected disappearance of the epidemic threshold ${\lambda }_c$ in the thermodynamic limit and an inverse participation ratio $\sim 0.25$, meaning localization in case of disassortative weight scheme. Contrarily, for the multiplicative weights and the unweighted trees, this value vanishes in the thermodynamic limit, suggesting only weak rare-region effects in agreement with the dynamical simulations. Strong corrections to the mean-field behavior in case of disassortative weights explains the concave shape of the order parameter $\rho \left(\lambda \right)$ at the transition point. Application of this method to other models may reveal interesting rare-region effects, Griffiths phases as the consequence of quenched topological heterogeneities.

Processes that are far both from equilibrium and from phase transition are studied. It is shown that a process with mean velocity that exhibits power-law growth in time can be analyzed using the Langevin equation with multiplicative noise. The solution to the corresponding Fokker-Planck equation is derived. Results of the numerical solution of the Langevin equation and simulation of the motion of particles in a billiard system with a time-dependent boundary are presented.