A self-generation mechanism of magnetic dipoles and the anomalous energy dissipation of fast electrons in a magnetized beam-plasma system are presented. Based on two-dimensional particle-in-cell simulations, it is found that the magnetic dipoles are self-organized and play important roles in the beam electron energy dissipation. These dipoles drift slowly in the direction of the return flow with a quasisteady velocity, which depends upon the magnetic amplitude of the dipole and the imposed external magnetic field. This dipole formation provides a mechanism for the anomalous energy dissipation of a relativistic electron beam, which would play an important role in collisionless shock and ion shock acceleration.

We investigate coupled nonlinear Schrödinger equations (NLSEs) with variable coefficients and gain. The coupled NLSE is a model equation for optical soliton propagation and their interaction in a multimode fiber medium or in a fiber array. By using Hirota's bilinear method, we obtain the bright-bright, dark-bright combinations of a one-soliton solution (1SS) and two-soliton solutions (2SS) for an $n$-coupled NLSE with variable coefficients and gain. Crucial properties of two-soliton (dark-bright pair) interactions, such as elastic and inelastic interactions and the dynamics of soliton bound states, are studied using asymptotic analysis and graphical analysis. We show that a bright $2\text{-soliton}$, in addition to elastic interactions, also exhibits multiple inelastic interactions. A dark $2\text{-soliton}$, on the other hand, exhibits only elastic interactions. We also observe a breatherlike structure of a bright 2-soliton, a feature that become prominent with gain and disappears as the amplitude acquires a minimum value, and after that the solitons remain parallel. The dark $2\text{-soliton}$, however, remains parallel irrespective of the gain. The results found by us might be useful for applications in soliton control, a fiber amplifier, all optical switching, and optical computing.

Current implementations of fluctuating lattice Boltzmann equations (FLBEs) describe single component fluids. In this paper, a model based on the continuum kinetic Boltzmann equation for describing multicomponent fluids is extended to incorporate the effects of thermal fluctuations. The thus obtained fluctuating Boltzmann equation is first linearized to apply the theory of linear fluctuations, and expressions for the noise covariances are determined by invoking the fluctuation-dissipation theorem directly at the kinetic level. Crucial for our analysis is the projection of the Boltzmann equation onto the orthonormal Hermite basis. By integrating in space and time the fluctuating Boltzmann equation with a discrete number of velocities, the FLBE is obtained for both ideal and nonideal multicomponent fluids. Numerical simulations are specialized to the case where mean-field interactions are introduced on the lattice, indicating a proper thermalization of the system.

We demonstrate reservoir computing with a physical system using a single autonomous Boolean logic element with time-delay feedback. The system generates a chaotic transient with a window of consistency lasting between 30 and 300 ns, which we show is sufficient for reservoir computing. We then characterize the dependence of computational performance on system parameters to find the best operating point of the reservoir. When the best parameters are chosen, the reservoir is able to classify short input patterns with performance that decreases over time. In particular, we show that four distinct input patterns can be classified for 70 ns, even though the inputs are only provided to the reservoir for 7.5 ns.

The computation of the critical exponent $\eta$ characterizing the universal elastic behavior of crystalline membranes in the flat phase continues to represent challenges to theorists as well as computer simulators that manifest themselves in a considerable spread of numerical results for $\eta$ published in the literature. We present additional insight into this problem that results from combining Wilson's momentum shell renormalization-group method with the power of modern computer simulations based on the Fourier Monte Carlo algorithm. After discussing the ideas and difficulties underlying this combined scheme, we present a calculation of the renormalization-group flow of the effective two-dimensional Young modulus for momentum shells of different thickness. Extrapolation to infinite shell thickness allows us to produce results in reasonable agreement with those obtained by functional renormalization group or by Fourier Monte Carlo simulations in combination with finite-size scaling. Moreover, our method allows us to obtain a decent estimate for the value of the Wegner exponent $\omega$ that determines the leading correction to scaling, which in turn allows us to refine our numerical estimate for $\eta$ previously obtained from precise finite-size scaling data.

Motivated by the recently observed phenomenon of topology trivialization of potential energy landscapes (PELs) for several statistical mechanics models, we perform a numerical study of the finite-size 2-spin spherical model using both numerical polynomial homotopy continuation and a reformulation via non-Hermitian matrices. The continuation approach computes all of the complex stationary points of this model while the matrix approach computes the real stationary points. Using these methods, we compute the average number of stationary points while changing the topology of the PEL as well as the variance. Histograms of these stationary points are presented along with an analysis regarding the complex stationary points. This work connects topology trivialization to two different branches of mathematics: algebraic geometry and catastrophe theory, which is fertile ground for further interdisciplinary research.

The effect of a uniform (nonspatially varying) external field on the liquid-vapor interface of the Stockmayer fluid (Lennard-Jones particles embedded with a point dipole) has been investigated by molecular-dynamics simulations. The long-ranged parts of both the dipole and Lennard-Jones interactions are treated using an Ewald summation, which removes the effects of the cutoff. The direction of the field shifts the critical point and interfacial properties in different directions. For an external field parallel to the interface, the critical temperature increases, while for a field applied perpendicular to the interface, it decreases. The effects of the field on surface tension and interfacial width are also investigated. For zero field, dipoles near the liquid-vapor interface show a weak orientation parallel to the interface. For fields parallel to the interface, ordering in the liquid phase is greater than the vapor, while for fields perpendicular to the interface, the opposite is true.

In directional solidification of binary eutectics, it is often observed that two-phase lamellar growth patterns grow tilted with respect to the direction z of the imposed temperature gradient. This crystallographic effect depends on the orientation of the two crystal phases $\alpha$ and $\beta$ with respect to z. Recently, an approximate theory was formulated that predicts the lamellar tilt angle as a function of the anisotropy of the free energy of the solid($\alpha$)-solid($\beta$) interphase boundary. We use two different numerical methods—phase field (PF) and dynamic boundary integral (BI)—to simulate the growth of steady periodic patterns in two dimensions as a function of the angle ${\theta }_R$ between z and a reference crystallographic axis for a fixed relative orientation of $\alpha$ and $\beta$ crystals, that is, for a given anisotropy function (Wulff plot) of the interphase boundary. For Wulff plots without unstable interphase-boundary orientations, the two simulation methods are in excellent agreement with each other and confirm the general validity of the previously proposed theory. In addition, a crystallographic “locking” of the lamellae onto a facet plane is well reproduced in the simulations. When unstable orientations are present in the Wulff plot, it is expected that two distinct values of the tilt angle can appear for the same crystal orientation over a finite ${\theta }_R$ range. This bistable behavior, which has been observed experimentally, is well reproduced by BI simulations but not by the PF model. Possible reasons for this discrepancy are discussed.

Typical ferroelectric $\text{Sm-}{C}^{*}$ liquid crystal (FLC) cells exhibit a voltage threshold for switching from one stable state to another despite the FLC's response being inherently continuous and thresholdless (free rotation of the director around the tilt cone). This switching threshold is due to FLC-surface interactions and to the chevron smectic structure commonly formed in cells. It is shown here that the FLC electrostatic energy contribution $\sim P_S{}^2$ responsible for thresholdless switching of high-$P_S$ FLCs also plays a key role in the bistable switching of lower-$P_S$ FLCs. Among the consequences are that it can be difficult to lower a cell's threshold below $V_{TB}\sim 3.4{(B/{\varepsilon }_F)}^{1/2}$ ($B$ and ${\varepsilon }_F$ are the FLC's elastic and dielectric constants), that a cell's threshold becomes independent of cell thickness once it substantially exceeds the characteristic length ${\xi }_P={\left({\varepsilon }_{F}B\right)}^{1/2}/P_S$, and that there are conditions under which alignment layer capacitance can decrease rather than increase the threshold (i.e., transition to thresholdless switching). A model that predicts and explains these behaviors is presented along with threshold measurements of representative FLC cells.

Potassium accumulation and diffusion during neuronal epileptiform activity have been observed experimentally, and potassium lateral diffusion has been suggested to play an important role in nonsynaptic neuron networks. We adopt a hippocampal CA1 pyramidal neuron network in a zero-calcium condition to better understand the influence of extracellular potassium dynamics on the stimulus-induced activity. The potassium concentration in the interstitial space for each neuron is regulated by potassium currents, ${\mathrm{Na}}^{+}\text{−}{\mathrm{K}}^{+}$ pumps, glial buffering, and ion diffusion. In addition to potassium diffusion, nearby neurons are also coupled through gap junctions. Our results reveal that the latency of the first spike responding to stimulus monotonically decreases with increasing gap-junction conductance but is insensitive to potassium diffusive coupling. The duration of network oscillations shows a bell-like shape with increasing potassium diffusive coupling at weak gap-junction coupling. For modest electrical coupling, there is an optimal ${\text{K}}^{+}$ diffusion strength, at which the flow of potassium ions among the network neurons appropriately modulates interstitial potassium concentrations in a degree that provides the most favorable environment for the generation and continuance of the action potential waves in the network.

The line tension or edge energy of bilayer membranes self-assembled from binary amphiphilic molecules is studied using self-consistent-field theory (SCFT). Specifically, solutions of the SCFT equations corresponding to an infinite membrane with a circular pore, or an open membrane, are obtained for a coarse-grained model in which the amphiphilic species and hydrophilic solvents are represented by $AB\mathrm{and}ED$ diblock copolymers and $C$ homopolymers, respectively. The edge energy of the membrane is extracted from the free energy of the open membranes. Results for membranes composed of mixtures of symmetric and cone- or inverse cone-shaped amphiphilic molecules with neutral and/or repulsive interactions are obtained and analyzed. It is observed that an increase in the concentration of the cone-shaped species leads to a decrease of the line tension. In contrast, adding inverse cone-shaped copolymers results in an increase of the line tension. Furthermore, the density profile of the copolymers reveals that the line tension is regulated by the distribution of the amphiphiles at the bilayer edge.

We analyze the dynamics of a distribution circuit loaded with many induction motors and subjected to sudden changes in voltage at the beginning of the circuit. As opposed to earlier work by Duclut et al. [Phys. Rev. E 87, 062802 (2013)], the motors are disordered, i.e., the mechanical torque applied to the motors varies in a random manner along the circuit. In spite of the disorder, many of the qualitative features of a homogeneous circuit persist, e.g., long-range motor-motor interactions mediated by circuit voltage and electrical power flows result in coexistence of the spatially extended and propagating normal and stalled phases. We also observed a new phenomenon absent in the case without inhomogeneity or disorder. Specifically, the transition front between the normal and stalled phases becomes somewhat random, even when the front is moving very slowly or is even stationary. Motors within the blurred domain appear in a normal or stalled state depending on the local configuration of the disorder. We quantify the effects of the disorder and discuss the statistics of distribution dynamics, e.g., the front position and width, total active and reactive consumption of the feeder, and maximum clearing time.

This paper studies the Sznajd model for opinion formation in a population connected through a general network. A master equation describing the time evolution of opinions is presented and solved in a mean-field approximation. Although quite simple, this approximation allows us to capture the most important features regarding the steady states of the model. When spontaneous opinion changes are included, a discontinuous transition from consensus to polarization can be found as the rate of spontaneous change is increased. In this case we show that a hybrid mean-field approach including interactions between second nearest neighbors is necessary to estimate correctly the critical point of the transition. The analytical prediction of the critical point is also compared with numerical simulations in a wide variety of networks, in particular Barabási-Albert networks, finding reasonable agreement despite the strong approximations involved. The same hybrid approach that made it possible to deal with second-order neighbors could just as well be adapted to treat other problems such as epidemic spreading or predator-prey systems.

Mutually connected components (MCCs) play an important role as a measure of resilience in the study of interdependent networks. Despite their importance, an efficient algorithm to obtain the statistics of all MCCs during the removal of links has thus far been absent. Here, using a well-known fully dynamic graph algorithm, we propose an efficient algorithm to accomplish this task. We show that the time complexity of this algorithm is approximately $O\left(N^{1.2}\right)$ for random graphs, which is more efficient than $O\left(N^2\right)$ of the brute-force algorithm. We confirm the correctness of our algorithm by comparing the behavior of the order parameter as links are removed with existing results for three types of double-layer multiplex networks. We anticipate that this algorithm will be used for simulations of large-size systems that have been previously inaccessible.

In this paper we propose a quantum algorithm to measure the similarity between a pair of unattributed graphs. We design an experiment where the two graphs are merged by establishing a complete set of connections between their nodes and the resulting structure is probed through the evolution of continuous-time quantum walks. In order to analyze the behavior of the walks without causing wave function collapse, we base our analysis on the recently introduced quantum Jensen-Shannon divergence. In particular, we show that the divergence between the evolution of two suitably initialized quantum walks over this structure is maximum when the original pair of graphs is isomorphic. We also prove that under special conditions the divergence is minimum when the sets of eigenvalues of the Hamiltonians associated with the two original graphs have an empty intersection.

Chimera states are complex spatio-temporal patterns that consist of coexisting domains of spatially coherent and incoherent dynamics. This counterintuitive phenomenon was first observed in systems of identical oscillators with symmetric coupling topology. Can one overcome these limitations? To address this question, we discuss the robustness of chimera states in networks of FitzHugh-Nagumo oscillators. Considering networks of inhomogeneous elements with regular coupling topology, and networks of identical elements with irregular coupling topologies, we demonstrate that chimera states are robust with respect to these perturbations and analyze their properties as the inhomogeneities increase. We find that modifications of coupling topologies cause qualitative changes of chimera states: additional random links induce a shift of the stability regions in the system parameter plane, gaps in the connectivity matrix result in a change of the multiplicity of incoherent regions of the chimera state, and hierarchical geometry in the connectivity matrix induces nested coherent and incoherent regions.

We introduce a concept of noise robustness in dynamical systems with noise and argue that this concept is essential to guarantee the reliability of physical random bit generators (RBGs). As an example of promising physical RBGs we consider a chaotic laser system and show that it has the property of noise robustness with respect to changes in the temporal correlation of the noise source. Moreover, employing a simple model of tangent space dynamics, we give a theoretical interpretation of the numerical results and in particular show that the Lyapunov exponent determines a theoretical boundary of a noise-robust region in parameter space. These theoretical results are expected to be significant not only for chaotic lasers, but also for a broad class of chaotic dynamical systems with correlated noise.

In this paper, we use molecular dynamics (MD) simulations to study the mean free path distribution of nonequilibrium gases in micronanochannel and to model the Knudsen (Kn) layer effect. It is found that the mean free path is significantly reduced near the wall and rather insensitive to flow types (Poiseuille or Couette). The Cercignani relation between the mean free path and the viscosity is adopted to capture the velocity behavior of the special zone in the framework of the extended Navier-Stokes (NS) equations. MD simulations of flows are carried out at different Kn numbers. Results are then compared with the theoretical model.

The mathematical model of a steadily propagating Saffman-Taylor finger in a Hele-Shaw channel has applications to two-dimensional interacting streamer discharges which are aligned in a periodic array. In the streamer context, the relevant regularization on the interface is not provided by surface tension but instead has been postulated to involve a mechanism equivalent to kinetic undercooling, which acts to penalize high velocities and prevent blow-up of the unregularized solution. Previous asymptotic results for the Hele-Shaw finger problem with kinetic undercooling suggest that for a given value of the kinetic undercooling parameter, there is a discrete set of possible finger shapes, each analytic at the nose and occupying a different fraction of the channel width. In the limit in which the kinetic undercooling parameter vanishes, the fraction for each family approaches $1/2$, suggesting that this “selection” of $1/2$ by kinetic undercooling is qualitatively similar to the well-known analog with surface tension. We treat the numerical problem of computing these Saffman-Taylor fingers with kinetic undercooling, which turns out to be more subtle than the analog with surface tension, since kinetic undercooling permits finger shapes which are corner-free but not analytic. We provide numerical evidence for the selection mechanism by setting up a problem with both kinetic undercooling and surface tension and numerically taking the limit that the surface tension vanishes.

Reflectivity of shocked compressed xenon plasma is calculated within the framework of the density functional theory approach. Dependencies on the frequency of incident radiation and on the plasma density are analyzed. The Fresnel formula for the reflectivity is used. The longitudinal expression in the long-wavelength limit is applied for the calculation of the imaginary part of the dielectric function. The real part of the dielectric function is calculated by means of the Kramers-Kronig transformation. The results are compared with experimental data. The approach for the calculation of plasma frequency is developed.