The fractional operator is introduced to an activator-inhibitor system to describe species anomalous superdiffusion. The effects of the superdiffusive exponent on pattern formation and pattern selection are studied. Our linear stability analysis shows that the wave number of the Turing pattern increases with the superdiffusive exponent. A weakly nonlinear analysis yields a system of amplitude equations and the analysis of these amplitude equations predicts parameter regimes where hexagons, stripes, and their coexistence are expected. Numerical simulations of the activator-inhibitor model near the stability boundaries confirm our analytical results. Since diffusion rate manifests in both diffusion constant and diffusion exponent, we numerically explore their interactions on the emergence of Turing patterns. When the activator and inhibitor have different superdiffusive exponents, we find that the critical ratio of the diffusion rate of the inhibitor to the activator, required for the formation of the Turing pattern, increases monotonically with the superdiffusive exponent. We conclude that small ratio (than unity) of anomalous diffusion exponent between the inhibitor and activator is more likely to promote the emergence of the Turing pattern, relative to the normal diffusion.

The topological complexity inherent to all porous media imparts persistent chaotic advection under steady flow conditions, which, in concert with the no-slip boundary condition, generates anomalous transport. We explore the impact of this mechanism upon longitudinal dispersion via a model random porous network and develop a continuous-time random walk that predicts both preasymptotic and asymptotic transport. In the absence of diffusion, the ergodicity of chaotic fluid orbits acts to suppress longitudinal dispersion from ballistic to superdiffusive transport, with asymptotic variance scaling as ${\sigma }_L^2\left(t\right)\sim t^2/{(lnt)}^3$. These results demonstrate that anomalous transport is inherent to homogeneous porous media and has significant implications for macrodispersion.

We perform a Leslie-type perturbative treatment on stably stratified turbulence with Boussinesq approximation, where the buoyancy terms in the corresponding dynamical equations are treated as perturbations against the isotropic background fields. Thus we calculate the anisotropic corrections to various correlation functions, namely, velocity-velocity, temperature-temperature, and velocity-temperature correlations, up to second order in this scheme. We find that the prefactors associated with the anisotropic corrections depend on the energy flux, scalar flux, Kolmogorov constant, Batchelor constant, and the eddy-damping amplitudes. The correlation functions further yield the anisotropic parts of the energy and mean-square temperature spectra as $k^{-3}$ and the anisotropic buoyancy spectrum as $k^{-7/3}$. The resulting angle-dependent energy density is found to be concentrated predominantly around the vertical wave vector signifying layered structures in the physical space.

We discuss implementation of the lattice-switching Monte Carlo method (LSMC) as a binary sampling between two synchronized Markov chains exploring separated minima in the potential energy landscape. When expressed in this fashion, the generalization to more complex crystals is straightforward. We extend the LSMC method to a flexible model of linear alkanes, incorporating bond length and angle constraints. Within this model, we accurately locate a transition between two polymorphs of $n$-butane with increasing density, and suggest this as a benchmark problem for other free-energy methods.

Most of the signals recorded in experiments are inevitably contaminated by measurement noise. Hence, it is important to understand the effect of such noise on estimating causal relations between such signals. A primary tool for estimating causality is Granger causality. Granger causality can be computed by modeling the signal using a bivariate autoregressive (AR) process. In this paper, we greatly extend the previous analysis of the effect of noise by considering a bivariate AR process of general order $p$. From this analysis, we analytically obtain the dependence of Granger causality on various noise-dependent system parameters. In particular, we show that measurement noise can lead to spurious Granger causality and can suppress true Granger causality. These results are verified numerically. Finally, we show how true causality can be recovered numerically using the Kalman expectation maximization algorithm.

The shear viscosity in the dilute regime of a model for confined granular matter is studied by simulations and kinetic theory. The model consists on projecting into two dimensions the motion of vibrofluidized granular matter in shallow boxes by modifying the collision rule: besides the restitution coefficient that accounts for the energy dissipation, there is a separation velocity that is added in each collision in the normal direction. The two mechanisms balance on average, producing stationary homogeneous states. Molecular dynamics simulations show that in the steady state the distribution function departs from a Maxwellian, with cumulants that remain small in the whole range of inelasticities. The shear viscosity normalized with stationary temperature presents a clear dependence with the inelasticity, taking smaller values compared to the elastic case. A Boltzmann-like equation is built and analyzed using linear response theory. It is found that the predictions show an excellent agreement with the simulations when the correct stationary distribution is used but a Maxwellian approximation fails in predicting the inelasticity dependence of the viscosity. These results confirm that transport coefficients depend strongly on the mechanisms that drive them to stationary states.

Correlations between intrinsic dynamics and local topology have become a new trend in the study of synchronization in complex networks. In this paper, we investigate the influence of topology on the dynamics of networks made up of second-order Kuramoto oscillators. In particular, based on mean-field calculations, we provide a detailed investigation of cluster explosive synchronization (CES) [Phys. Rev. Lett. 110, 218701 (2013)] in scale-free networks as a function of several topological properties. Moreover, we investigate the robustness of discontinuous transitions by including an additional quenched disorder, and we show that the phase coherence decreases with increasing strength of the quenched disorder. These results complement the previous findings regarding CES and also fundamentally deepen the understanding of the interplay between topology and dynamics under the constraint of correlating natural frequencies and local structure.

We identify the mechanism behind the existence of intensity-induced chimera states in globally coupled oscillators. We find that the effect of intensity in the system is to cause multistability by increasing the number of fixed points. This in turn increases the number of multistable attractors, and we find that their stability is determined by the strength of coupling. This causes the coexistence of different collective states in the system depending upon the initial state. We demonstrate that intensity-induced chimera is generic to both periodic and chaotic systems. We discuss possible applications of our results to real-world systems like the brain and spin torque nano-oscillators.

In a recent publication [Schmidt et al., Phys. Plasmas 21, 080701 (2014)], we reported results on beam-target interaction experiments that have been carried out at the CERN HiRadMat (High Radiation to Materials) facility using extended solid copper cylindrical targets that were irradiated with a 440-GeV proton beam delivered by the Super Proton Synchrotron (SPS). On the one hand, these experiments confirmed the existence of hydrodynamic tunneling of the protons that leads to substantial increase in the range of the protons and the corresponding hadron shower in the target, a phenomenon predicted by our previous theoretical investigations [Tahir et al., Phys. Rev. ST Accel. Beams 25, 051003 (2012)]. On the other hand, these experiments demonstrated that the beam heated part of the target is severely damaged and is converted into different phases of high energy density (HED) matter, as suggested by our previous theoretical studies [Tahir et al., Phys. Rev. E 79, 046410 (2009)]. The latter confirms that the HiRadMat facility can be used to study HED physics. In the present paper, we give details of the numerical simulations carried out to understand the experimental measurements. These include the evolution of the physical parameters, for example, density, temperature, pressure, and the internal energy in the target, during and after the irradiation. This information is important in order to determine the region of the HED phase diagram that can be accessed in such experiments. These simulations have been done using the energy deposition code fluka and a two–dimensional hydrodynamic code, big2, iteratively.

Following on our previous work [S. Delong, B. E. Griffith, E. Vanden-Eijnden, and A. Donev, Phys. Rev. E 87, 033302 (2013)], we develop temporal integrators for solving Langevin stochastic differential equations that arise in fluctuating hydrodynamics. Our simple predictor-corrector schemes add fluctuations to standard second-order deterministic solvers in a way that maintains second-order weak accuracy for linearized fluctuating hydrodynamics. We construct a general class of schemes and recommend two specific schemes: an explicit midpoint method and an implicit trapezoidal method. We also construct predictor-corrector methods for integrating the overdamped limit of systems of equations with a fast and slow variable in the limit of infinite separation of the fast and slow time scales. We propose using random finite differences to approximate some of the stochastic drift terms that arise because of the kinetic multiplicative noise in the limiting dynamics. We illustrate our integrators on two applications involving the development of giant nonequilibrium concentration fluctuations in diffusively mixing fluids. We first study the development of giant fluctuations in recent experiments performed in microgravity using an overdamped integrator. We then include the effects of gravity and find that we also need to include the effects of fluid inertia, which affects the dynamics of the concentration fluctuations greatly at small wave numbers.

To be useful, steady-state thermodynamics (SST) must be self-consistent and have predictive value. Consistency of SST was recently verified for driven lattice gases under global weak exchange. Here I verify consistency of SST under local (pointwise) exchange, but only in the limit of a vanishing exchange rate; for a finite exchange rate the coexisting systems have different chemical potentials. I consider the lattice gas with nearest-neighbor exclusion on the square lattice, with nearest-neighbor hopping, and with hopping to both nearest and next-nearest neighbors. I show that SST does not predict the coexisting densities under a nonuniform drive or in the presence of a nonuniform density provoked by a hard wall or nonuniform transition rates. The steady-state chemical potential profile is, moreover, nonuniform at coexistence, contrary to the basic principles of thermodynamics. Finally, I discuss examples of a pair of systems possessing identical steady states but which do not coexist when placed in contact. The results of these studies confirm the validity of SST for coexistence between spatially uniform systems but cast serious doubt on its consistency and predictive value in systems with a finite rate of particle exchange between coexisting regions exhibiting a nonuniform particle density.

The derivation of general performance benchmarks is important in the design of highly optimized heat engines and refrigerators. To obtain them, one may model phenomenologically the leading sources of irreversibility ending up with results that are model independent, but limited in scope. Alternatively, one can take a simple physical system realizing a thermodynamic cycle and assess its optimal operation from a complete microscopic description. We follow this approach in order to derive the coefficient of performance at maximum cooling rate for any endoreversible quantum refrigerator. At striking variance with the universality of the optimal efficiency of heat engines, we find that the cooling performance at maximum power is crucially determined by the details of the specific system-bath interaction mechanism. A closed analytical benchmark is found for endoreversible refrigerators weakly coupled to unstructured bosonic heat baths: an ubiquitous case study in quantum thermodynamics.

The canonical statistics describes the statistical properties of an open system by assuming its coupling with the heat bath is infinitesimal in comparison with the total energy in thermodynamic limit. In this paper, we generally derive a noncanonical density matrix for the open system with a finite coupling to the heat bath, which deforms the energy shell to effectively modify the conventional canonical way. The obtained noncanonical distribution reflects the back action of system on the bath and thus depicts the statistical correlations between two subsystems by the mutual information as a result of energy conservation.

We demonstrate how the matrix-valued mode-coupling theory of the glass transition and glassy dynamics in planar confinement converges to the corresponding theory for two-dimensional (2D) planar and the three-dimensional bulk liquid, provided the wall potential satisfies certain conditions. Since the mode-coupling theory relies on the static properties as input, the emergence of a homogeneous limit for the matrix-valued intermediate scattering functions is directly connected to the convergence of the corresponding static quantities to their conventional counterparts. We show that the 2D limit is more subtle than the bulk limit, in particular, the in-planar dynamics decouples from the motion perpendicular to the walls. We investigate the frozen-in parts of the intermediate scattering function in the glass state and find that the limits time $t\to \infty$ and effective wall separation $L\to 0$ do not commute due to the mutual coupling of the residual transversal and lateral force kernels.

Low-frequency nonlinear wave dynamics is investigated in a two-dimensional inhomogeneous electron magnetohydrodynamic (EMHD) plasma in the presence of electron viscosity. In the long-wavelength limit, the dynamics of the wave is found to be governed by a novel nonlinear equation. The result of the moving-frame nonlinear analysis is noteworthy, which shows that this nonlinear equation does have a breather solution and electron viscosity is responsible for the breather. A breather is a nonlinear wave in which energy accumulates in a localized and oscillatory manner. Analytical solution and time-dependent numerical simulation of this novel equation reveal the collapse of a soliton (localized pulse) into a weak noise shelf and formation of shocklike structures.

We develop an approach to solving numerically the time-dependent Schrödinger equation when it includes source terms and time-dependent potentials. The approach is based on the generalized Crank-Nicolson method supplemented with an Euler-MacLaurin expansion for the time-integrated nonhomogeneous term. By comparing the numerical results with exact solutions of analytically solvable models, we find that the method leads to precision comparable to that of the generalized Crank-Nicolson method applied to homogeneous equations. Furthermore, the systematic increase in precision generally permits making estimates of the error.

We provide an algorithm for the construction of orthonormal multivariate polynomials that are symmetric with respect to the interchange of any two coordinates on the unit hypercube and are constrained to the hyperplane where the sum of the coordinates is one. These polynomials form a basis for the expansion of bosonic light-front momentum-space wave functions, as functions of longitudinal momentum, where momentum conservation guarantees that the fractions are on the interval $[0,1]$ and sum to one. This generalizes earlier work on three-boson wave functions to wave functions for arbitrarily many identical bosons. A simple application in two-dimensional ${\varphi }^4$ theory illustrates the use of these polynomials.

In this paper, a phase-field-based lattice Boltzmann (LB) model is proposed for axisymmetric multiphase flows. Modified equilibrium distribution functions and some source terms are properly added into the evolution equations such that multiphase flows in the axisymmetric coordinate system can be described. Different from previous axisymmetric LB multiphase models, the added source terms that arise from the axisymmetric effect contain no additional gradients, and therefore the present model is much simpler. Furthermore, through the Chapmann-Enskog analysis, the axisymmetric Chan-Hilliard equation and Navier-Stokes equations can be exactly derived from the present model. The model is also able to deal with flows with density contrast. A variety of numerical experiments, including planar and curve interfaces, an elongation field, a static droplet, a droplet oscillation, breakup of a liquid thread, and dripping of a liquid droplet under gravity, have been conducted to test the proposed model. It is found that the present model can capture accurate interface and the numerical results of multiphase flows also agree well with the analytical solutions and/or available experimental data.