Physical Review. E, Statistical, Nonlinear, and Soft Matter Physics

This corrects the article DOI: 10.1103/PhysRevE.91.040901.

This corrects the article DOI: 10.1103/PhysRevE.92.033003.

This corrects the article DOI: 10.1103/PhysRevE.90.022142.

This corrects the article DOI: 10.1103/PhysRevE.92.052715.

This corrects the article DOI: 10.1103/PhysRevE.92.043308.

In this Reply we respond to the criticism of Frewer et al. presented in the Comment on "Statistical symmetries of the Lundgren-Monin-Novikov hierarchy" by Wacławczyk et al. [Phys. Rev. E 90, 013022 (2014)]. We discuss physical interpretation of the statistical symmetries, and respond to criticism on the violation of the causality principle. We derive the Lundgren-Monin-Novikov equations for a flow with boundaries. Last, we stress that our work addressed the phenomenon of "external intermittency" (separation between laminar and turbulent flow), and not "internal intermittency" (strong fluctuations at small scales)...

The article by M. Wacławczyk et al. [Phys. Rev. E 90, 013022 (2014)] proposes two new statistical symmetries in the classical theory for turbulent hydrodynamic flows. In this Comment, however, we show that both symmetries are unphysical due to violating the principle of causality. In addition, they must get broken in order to be consistent with all physical constraints naturally arising in the statistical Lundgren-Monin-Novikov (LMN) description of turbulence. As a result, we state that besides the well-known classical symmetries of the LMN equations no new statistical symmetries exist...

In a recent paper [M. Cristelli, A. Zaccaria, and L. Pietronero, Phys. Rev. E 85, 066108 (2012)], the authors analyzed the relation between skewness and kurtosis for complex dynamical systems, and they identified two power-law regimes of non-Gaussianity, one of which scales with an exponent of 2 and the other with 4/3. They concluded that the observed relation is a universal fact in complex dynamical systems. In this Comment, we test the proposed universal relation between skewness and kurtosis with a large number of synthetic data, and we show that in fact it is not a universal relation and originates only due to the small number of data points in the datasets considered...

The equations of motion of a single particle subject to an arbitrary electric and a static magnetic field form a Poisson system. We present a second-order time integration method which preserves well the Poisson structure and compare it to commonly used algorithms, such as the Boris scheme. All the methods are represented in a general framework of splitting methods. We use the so-called ϕ functions, which give efficient ways for both analyzing and implementing the algorithms. Numerical experiments show an excellent long term stability for the method considered...

In this work, a numerical scheme based on multipoles and Stokeslet is proposed for calculating the radiation force acting on a single rigid sphere in a viscous fluid. First-order velocity and pressure are obtained from the multipole series solution, and the volumetric force in the acoustic streaming is subsequently calculated from the first-order velocity and pressure. The acoustic streaming equations are solved using the Stokeslet method within a finite domain descretized by tetrahedral elements. The boundary conditions for streaming are imposed using the weighted residue method to obtain the unknown coefficients in the multipole series expansion for the second-order velocity potentials...

Many physical systems can be modeled as large sets of domains "glued" together along boundaries-biological cells meet along cell membranes, soap bubbles meet along thin films, countries meet along geopolitical boundaries, and metallic crystals meet along grain interfaces. Each class of microstructures results from a complex interplay of initial conditions and particular evolutionary dynamics. The statistical steady-state microstructure resulting from isotropic grain growth of a polycrystalline material is canonical in that it is the simplest example of a cellular microstructure resulting from a gradient flow of an energy that is directly proportional to the total length or area of all cell boundaries...

Population annealing is an efficient sequential Monte Carlo algorithm for simulating equilibrium states of systems with rough free-energy landscapes. The theory of population annealing is presented, and systematic and statistical errors are discussed. The behavior of the algorithm is studied in the context of large-scale simulations of the three-dimensional Ising spin glass and the performance of the algorithm is compared to parallel tempering. It is found that the two algorithms are similar in efficiency though with different strengths and weaknesses...

We apply the event-chain Monte Carlo algorithm to the three-dimensional ferromagnetic Heisenberg model. The algorithm is rejection-free and also realizes an irreversible Markov chain that satisfies global balance. The autocorrelation functions of the magnetic susceptibility and the energy indicate a dynamical critical exponent z≈1 at the critical temperature, while that of the magnetization does not measure the performance of the algorithm. We show that the event-chain Monte Carlo algorithm substantially reduces the dynamical critical exponent from the conventional value of z≃2...

In this paper a full Eulerian lattice Boltzmann model is proposed for conjugate heat transfer. A unified governing equation with a source term for the temperature field is derived. By introducing the source term, we prove that the continuity of temperature and its normal flux at the interface is satisfied automatically. The curved interface is assumed to be zigzag lines. All physical quantities are recorded and updated on a Cartesian grid. As a result, any complicated treatment near the interface is avoided, which makes the proposed model suitable to simulate the conjugate heat transfer with complex interfaces efficiently...

We validate the application of our recent orbital-free density functional theory (DFT) approach [Phys. Rev. Lett. 113, 155006 (2014);] for the calculation of ionic and electronic transport properties of dense plasmas. To this end, we calculate the self-diffusion coefficient, the viscosity coefficient, the electrical and thermal conductivities, and the reflectivity coefficient of hydrogen and aluminum plasmas. Very good agreement is found with orbital-based Kohn-Sham DFT calculations at lower temperatures. Because the computational costs of the method do not increase with temperature, we can produce results at much higher temperatures than is accessible by the Kohn-Sham method...

A reconstruction methodology based on different-phase-neighbor (DPN) pixel swapping and multigrid hierarchical annealing is presented. The method performs reconstructions by starting at a coarse image and successively refining it. The DPN information is used at each refinement stage to freeze interior pixels of preformed structures. This preserves the large-scale structures in refined images and also reduces the number of pixels to be swapped, thereby resulting in a decrease in the necessary computational time to reach a solution...

We review the momentum exchange method to compute the flow force and torque on a submerged body in lattice-Boltzmann methods by presenting an alternative derivation. Our derivation does not depend on a particular implementation of the boundary conditions at the body surface, and it relies on general principles. After the introduction of the momentum exchange method in lattice-Boltzmann methods, some formulations were introduced to compute the fluid force on static and moving bodies. These formulations were introduced in a rather intuitive, ad hoc way...

Density matrix perturbation theory [Niklasson and Challacombe, Phys. Rev. Lett. 92, 193001 (2004)] is generalized to canonical (NVT) free-energy ensembles in tight-binding, Hartree-Fock, or Kohn-Sham density-functional theory. The canonical density matrix perturbation theory can be used to calculate temperature-dependent response properties from the coupled perturbed self-consistent field equations as in density-functional perturbation theory. The method is well suited to take advantage of sparse matrix algebra to achieve linear scaling complexity in the computational cost as a function of system size for sufficiently large nonmetallic materials and metals at high temperatures...

We consider a harmonically driven acoustic medium in the form of a (finite length) highly nonlinear granular crystal with an amplitude- and frequency-dependent boundary drive. Despite the absence of a linear spectrum in the system, we identify resonant periodic propagation whereby the crystal responds at integer multiples of the drive period and observe that this can lead to local maxima of transmitted force at its fixed boundary. In addition, we identify and discuss minima of the transmitted force ("antiresonances") between these resonances...

By invoking Bogoliubov's spectrum, we show that for the nonlinear Schrödinger equation, the modulation instability (MI) of its n=1 Fourier mode on a finite background automatically triggers a further cascading instability, forcing all the higher modes to grow exponentially in locked step with the n=1 mode. This fundamental insight, the enslavement of all higher modes to the n=1 mode, explains the formation of a triangular-shaped spectrum that generates the Akhmediev breather, predicts its formation time analytically from the initial modulation amplitude, and shows that the Fermi-Pasta-Ulam (FPU) recurrence is just a matter of energy conservation with a period twice the breather's formation time...