In the present Reply we show that a Comment casting doubts on the results of our recent paper [Fronczak, Fronczak, and Siudem, Phys. Rev. E 101, 022111 (2020)] results from a misunderstanding of the assumptions of our model and from overinterpretation based on this misunderstanding.

Imitation and aspiration update rules are frequently observed in human and animal populations. While the imitation process entails payoff comparisons with surroundings, the aspiration process refers to self-evaluation. This work explores the evolution of cooperation in dyadic games under the coexistence of these two dynamics in an infinitely large well-mixed population. Two situations have been explored: (i) individuals adopt either an imitation or aspiration update rule with a certain probability, and (ii) the entire population is divided into two groups where one group only uses imitative rules and the other obeys aspiration updating alone. Both premises have been modeled by taking an infinite approximation of the finite population. In particular, the second mixing principle follows an additive property: the outcome of the whole population is the weighted average of outcomes from imitators and aspiration-driven individuals. Our work progressively investigates several variants of aspiration dynamics under strong selection, encompassing symmetric, asymmetric, and adaptive aspirations, which then coalesce with imitative dynamics. We also demonstrate which of the update rules performs better, under different social dilemmas, by allowing the evolution of the preference of update rules besides strategies. Aspiration dynamics always outperform imitation dynamics in the prisoner's dilemma, however, in the chicken and stag-hunt games the predominance of either update rule depends on the level of aspirations as well as on the extent of greed and fear present in the system. Finally, we examine the coevolution of strategies, aspirations, and update rules which leads to a binary state of obeying either imitation or aspiration dynamics. In such a circumstance, when aspiration dynamics prevail over imitation dynamics, cooperators and defectors coexist to an equal extent.

We focus here on the thermodynamic properties of adsorbates formed by two-species $A+B\to \oslash$ reactions on a one-dimensional infinite lattice with heterogeneous “catalytic” properties. In our model hard-core $A$ and $B$ particles undergo continuous exchanges with their reservoirs and react when dissimilar species appear at neighboring lattice sites in presence of a “catalyst.” The latter is modeled by supposing either that randomly chosen bonds in the lattice promote reactions (Model I) or that reactions are activated by randomly chosen lattice sites (Model II). In the case of annealed disorder in spatial distribution of a catalyst we calculate the pressure of the adsorbate by solving three-site (Model I) or four-site (Model II) recursions obeyed by the corresponding averaged grand-canonical partition functions. In the case of quenched disorder, we use two complementary approaches to find exact expressions for the pressure. The first approach is based on direct combinatorial arguments. In the second approach, we frame the model in terms of random matrices; the pressure is then represented as an averaged logarithm of the trace of a product of random $3\times 3$ matrices—either uncorrelated (Model I) or sequentially correlated (Model II).

For a classical system with long-range interactions, a soft mode exists whenever a stationary state spontaneously breaks a continuous symmetry of the Hamiltonian. Besides that, if the corresponding coordinate associated to the symmetry breaking is periodic, then the same energy of the different stationary states and finite $N$ thermal fluctuations result in a superdiffusive motion of the center of mass for total zero momentum, that tends to a normal diffusion for very long times. As examples of this, we provide a two-dimensional self-gravitating system, a free electron laser, and the Hamiltonian mean-field (HMF) model. For the latter, a detailed theory for the motion of the center of mass is given. We also discuss how the coupling of the soft mode to the mean-field motion of individual particles may lead to strong chaotic behavior for a finite particle number, as illustrated by the HMF model.

Numerical simulations and finite-size scaling analysis have been performed to study the jamming and percolation behavior of elongated objects deposited on two-dimensional honeycomb lattices. The depositing particle is modeled as a linear array of length $k$ (so-called $k$-mer), maximizing the distance between first and last monomers in the chain. The separation between $k$-mer units is equal to the lattice constant. Hence, $k$ sites are occupied by a $k$-mer when adsorbed onto the surface. The adsorption process starts with an initial configuration, where all lattice sites are empty. Then, the sites are occupied following a random sequential adsorption mechanism. The process finishes when the jamming state is reached and no more objects can be deposited due to the absence of empty site clusters of appropriate size and shape. Jamming coverage ${\theta }_{j,k}$ and percolation threshold ${\theta }_{c,k}$ were determined for a wide range of values of $k$ ($2\le k\le 128$). The obtained results shows that $\left(i\right) {\theta }_{j,k}$ is a decreasing function with increasing $k$, being ${\theta }_{j,k\to \infty }=0.6007\left(6\right)$ the limit value for infinitely long $k$-mers; and $\left(ii\right) {\theta }_{c,k}$ has a strong dependence on $k$. It decreases in the range $2\le k<48$, goes through a minimum around $k=48$, and increases smoothly from $k=48$ up to the largest studied value of $k=128$. Finally, the precise determination of the critical exponents $\nu , \beta ,$ and $\gamma$ indicates that the model belongs to the same universality class as 2D standard percolation regardless of the value of $k$ considered.

We carry out Monte Carlo simulations of a colloidal fluid membrane with a free edge and composed of chiral rodlike viruses. The membrane is modeled by a triangular mesh of beads connected by bonds in which the bonds and beads are free to move at each Monte Carlo step. Since the constituent viruses are experimentally observed to twist only near the membrane edge, we use an effective energy that favors a particular sign of the geodesic torsion of the edge. The effective energy also includes the membrane bending stiffness, edge bending stiffness, and edge tension. We find three classes of membrane shapes resulting from the competition of the various terms in the free energy: branched shapes, chiral disks, and vesicles. Increasing the edge bending stiffness smooths the membrane edge, leading to correlations among the membrane normals at different points along the edge. The normalized power spectrum for edge displacements shows a peak with increasing preferred geodesic torsion. We also consider membrane shapes under an external force by fixing the distance between two ends of the membrane and finding the shape for increasing values of the distance between the two ends. As the distance increases, the membrane twists into a ribbon, with the force eventually reaching a plateau.

Collapse modes in compressed simple cubic (SC) and body-centered cubic (BCC) periodic arrangements of elastic frictionless beads were studied numerically using the discrete element method. Under pure hydrostatic compression, the SC arrangement tends to transform into a defective hexagonal close-packed or amorphous structure. The BCC assembly exhibits several modes of collapse, one of which, identified as cI16 structure, is consistent with the behavior of BCC metals Li and Na under high pressure. The presence of a deviatoric stress leads to the transformation of the BCC structure into face-centered cubic (FCC) one via the Bain path. The observed effects expand the knowledge on possible packings of soft elastic spheres and transformations between them, while providing an unexpected link with the mechanical behavior of certain atomic systems.

Nowadays metamaterials are at the focus of an intense research as promising for thermal and acoustic engineering. However, the computational cost associated to the large system size required for correctly simulating them imposes the use of finite-elements simulations, developing continuum models, able to grasp the physics at play without entering in the atomistic details. Still, a correct description should be able to reproduce not only the extrinsic scattering sources on waves propagation, as introduced by the metamaterial microstructure, but also the intrinsic wave attenuation of the material itself. This becomes dramatically important when the metamaterial is made out of a glass, which is intrinsically highly dissipative and with a wave attenuation strongly dependent on frequency. Here we propose a continuum mechanical model for a viscoelastic medium, able to bridge atomic and macroscopic scale in amorphous materials and describe phonon attenuation due to atomistic mechanisms, characterized by a defined frequency dependence. This represents a first decisive step for investigating the effect of a complex nano- or microstructure on acoustic attenuation, while including the atomistic contribution as well.

We show that similarly to the logarithmic mean-velocity profile in wall-bounded turbulence, the landscape topography presents an intermediate region with a logarithmic mean-elevation profile. Such profiles are present in complex topographies with channel branching and fractal river networks resulting from model simulation, controlled laboratory experiments, and natural landscapes. Dimensional and self-similarity arguments are used to corroborate this finding. We also tested the presence of logarithmic profiles in discrete, minimalist models of networks obtained from optimality principles (optimal channel networks) and directed percolation. The emergence of self-similar scaling appears as a robust outcome in dynamically different, but spatially bounded, complex systems, as a dimensional consequence of length-scale independence.

When excited at sufficiently high acoustic pressures, a wall-attached bubble may exhibit asymmetric nonspherical modes. These vibration modes can be decomposed over the set of spherical harmonics $Y_{nm}(\theta ,\varphi )$ for a degree $n$ and order $m$. We experimentally capture the time-resolved dynamics of asymmetric bubble oscillations in a top-view configuration. A spatiotemporal modal analysis is performed and allowed recovering the set of zonal ($m=0$), tesseral ($0<m<n$), and sectoral ($m=n$) spherical harmonics that develop at the bubble interface. The analysis of the surface instability thresholds reveals that the frequencies of asymmetric modes differ from the standard Lamb spectrum. In addition, the nondegeneracy of asymmetric modes for a given degree $n$ is evidenced by noncompletely overlapping resonance bands. Finally, the coexistence between zonal and sectoral modes is analyzed through their modal interaction, amplitude interplay and relation of phase, as well as their geometric compatibility.

Recent experimental results indicate that mixing is enhanced by a reciprocal flow induced inside a levitated droplet with an oscillatory deformation [T. Watanabe et al., Sci. Rep. 8, 10221 (2018)]. Generally, reciprocal flow cannot convect the solutes in time average, and agitation cannot take place. In the present paper, we focus on the diffusion process coupled with the reciprocal flow. We theoretically derive that the diffusion process can be enhanced by the reciprocal flow, and the results are confirmed via numerical calculation of the over-damped Langevin equation with a reciprocal flow.

The resources needed for particle-in-cell simulations of laser wakefield acceleration can be greatly reduced in many cases of interest using an envelope model. However, the inclusion of tunneling ionization in this time-averaged treatment of laser-plasma acceleration is not straightforward, since the statistical features of the electron beams obtained through ionization should ideally be reproduced without resolving the high-frequency laser oscillations. In this context, an extension of an already known envelope ionization procedure is proposed, valid also for laser pulses with higher intensities, which consists in adding the initial longitudinal drift to the newly created electrons within the laser pulse ionizing the medium. The accuracy of the proposed procedure is shown with both linear and circular polarization in a simple benchmark where a nitrogen slab is ionized by a laser pulse and in a more complex benchmark of laser plasma acceleration with ionization injection in the nonlinear regime. With this addition to the envelope ionization algorithm, the main phase space properties of the bunches injected in a plasma wakefield with ionization by a laser (charge, average energy, energy spread, rms sizes, and normalized emittance) can be estimated with accuracy comparable to a nonenvelope simulation with significantly reduced resources, even in cylindrical geometry. Through this extended algorithm, preliminary studies of ionization injection in laser wakefield acceleration can be easily carried out even on a laptop.

A kinetic inviscid flux (KIF) is proposed for simulating incompressible and compressible flows. It is constructed based on the direct modeling of multiscale flow behaviors, which is used in the gas-kinetic scheme (GKS), the unified gas-kinetic scheme (UGKS), the discrete unified gas-kinetic scheme (DUGKS), etc. In KIF, the discontinuities (such as the shock wave) that cannot be well resolved by mesh cells are mainly solved by the kinetic flux vector splitting (KFVS) method representing the free transport mechanism (or microscale mechanism), while in other flow regions that are smooth, the flow behavior is solved mainly by the central-scheme-like totally thermalized transport (TTT). The weights of KFVS and TTT in KIF is automatically determined by those in the theory of direct modeling. Two ways of choosing the weights in KIF are proposed, which are actually the weights adopted in the UGKS and the DUGKS, respectively. By using the test cases of the Sod shock tube, the rarefaction wave, the boundary layer of a flat plate, the cavity flow, hypersonic flow over a circular cylinder, the shock and turbulent boundary iteration, and transonic flow over a three-dimensional M6 wing, the validity and accuracy of the present method are examined. The KIF does not suffer from the carbuncle phenomenon, and does not introduce extra numerical viscosity in smooth regions. Especially in the case of hypersonic cylinder, it gives quite sharp and clear density and temperature contours. The KIF can be viewed as an inviscid-viscous splitting version of the GKS. By doing this splitting, it is easy to be used in traditional computational fluid dynamics frameworks. It can also be classified as a type in the numerical schemes based on the kinetic theory that are represented by the works of Sun et al. [Adv. Appl. Math. Mech. 8, 703 (2016)] and Ohwada et al. [J. Comput. Phys. 362, 131 (2018)], except the weights are determined by the weights of direct modeling.

We investigate how confinement may drastically change both the probability density of the first-encounter time and the associated survival probability in the case of two diffusing particles. To obtain analytical insights into this problem, we focus on two one-dimensional settings: a half-line and an interval. We first consider the case with equal particle diffusivities, for which exact results can be obtained for the survival probability and the associated first-encounter time density valid over the full time domain. We also evaluate the moments of the first-encounter time when they exist. We then turn to the case with unequal diffusivities and focus on the long-time behavior of the survival probability. Our results highlight the great impact of boundary effects in diffusion-controlled kinetics even for simple one-dimensional settings, as well as the difficulty of obtaining analytic results as soon as the translational invariance of such systems is broken.

The Bonabeau model of self-organized hierarchy formation is studied by using a piecewise linear approximation to the sigmoid function. Simulations of the piecewise-linear agent model show that there exist two-level and three-level hierarchical solutions and that each agent exhibits a transition from nonergodic to ergodic behaviors. Furthermore, by using a mean-field approximation to the agent model, it is analytically shown that there are asymmetric two-level solutions, even though the model equation is symmetric (asymmetry is introduced only through the initial conditions) and that linearly stable and unstable three-level solutions coexist. It is also shown that some of these solutions emerge through supercritical-pitchfork-like bifurcations in invariant subspaces. Existence and stability of the linear hierarchy solution in the mean-field model are also elucidated.

We consider a red-versus-blue coupled synchronization and spatial swarming (i.e., swarmalator) model that incorporates attraction and repulsion terms and an adversarial game of phases. The model exhibits behaviors such as spontaneous emergence of tactical manoeuvres of envelopment (e.g., flanking, pincer, and envelopment) that are often proposed in military theory or observed in nature. We classify these states based on a large set of features such as spatial densities, synchronization between clusters, and measures of cluster distances. These features are used to study the influence of coupling parameters on the expected presence of these states and the—sometimes sharp—transitions between them.

We present a method for evolving the projected Gross-Pitaevskii equation in an infinite rotating Bose-Einstein condensate, the ground state of which is a vortex lattice. We use quasiperiodic boundary conditions to investigate the behavior of the bulk superfluid in this system, in the absence of boundaries and edge effects. We also give the Landau gauge expression for the phase of a BEC subjected to these boundary conditions. Our spectral representation uses the eigenfunctions of the one-body Hamiltonian as basis functions. Since there is no known exact quadrature rule for these basis functions we approximately implement the projection associated with the energy cutoff, but we show that by choosing a suitably fine spatial grid the resulting error can be made negligible. We show how the convergence of this model is affected by simulation parameters such as the size of the spatial grid and the number of Landau levels. Adding dissipation, we use our method to find the lattice ground state for $N$ vortices. We can then perturb the ground-state, to investigate the melting of the lattice.

The discrete circle map is the archetypical example of a driven periodic system, showing a complex resonance structure under a change of the forcing frequency known as the devil's staircase. Adler's equation can be seen as the direct continuous equivalent of the circle map, describing locking effects in periodic systems with continuous forcing. This type of locking produces a single fundamental resonance tongue without higher-order resonances, and a devil's staircase is not observed. We show that, with harmonically modulated forcing, nonlinear oscillations close to a Hopf bifurcation generically reproduce the devil's staircase even in the continuous case. Experimental results on a semiconductor laser driven by a modulated optical signal show excellent agreement with our theoretical predictions. The locking appears as a modulation of the oscillation amplitude as well as the angular oscillation frequency. Our results show that by proper implementation of an external drive, additional regions of stable frequency locking can be introduced in systems which originally show only a single Adler-type resonance tongue. The induced resonances can be precisely controlled via the modulation parameters.

We investigate a simple forced harmonic oscillator with a natural frequency varying with time. It is shown that the time evolution of such a system can be written in a simplified form with Fresnel integrals, as long as the variation of the natural frequency is sufficiently slow compared to the time period of oscillation. Thanks to such a simple formulation, we found that a forced harmonic oscillator with a slowly varying natural frequency is essentially equivalent to diffraction of light.