Starting from Landauer's slogan “information is physical,” we revise and modify Landauer's principle stating that the erasure of information has a minimal price in the form of a certain quantity of free energy. We establish a direct link between the erasure cost and the work value of a piece of information and show that the former is essentially the length of the string's best compression by a reversible computation. We generalize the principle by deriving bounds on the free energy to be invested for—or gained from, for that matter—a general computation. We then revisit the second law of thermodynamics and compactly rephrase it (assuming the Church-Turing-Deutsch hypothesis that physical reality can be simulated by a universal Turing machine): Time evolutions are logically reversible—“the future fully remembers the past (but not necessarily vice versa).” We link this view to previous formulations of the second law, and we argue that it has a particular feature that suggests its “logico-informational” nature, namely, simulation resilience: If a computation faithfully simulates a physical process violating the law, then that very computation procedure violates it as well.

We examine two-level systems which undergo an instantaneous change in the energy level differences due to the application of an external protocol. In particular the behavior of the work distribution as well as the probability of second-law violations as the thermodynamic limit, and separately the quasistatic limit, is approached are investigated.

This article discusses the numerical result predicted by the quantum Langevin equation of the generalized diffusion function of a Brownian particle immersed in an Ohmic quantum bath of harmonic oscillators. The time dependence of the standard deviation of the reduced Wigner function of the system, obtained by integrating the whole function in the momentum space, is determined as well. The complexity of the equations leads to resort to a much simpler calculation based in the position correlation function. They are done for the three possible regimes of the system, namely, periodic, aperiodic, and overdamped. It is found in the periodic case that the generalized diffusion is a discontinuous function exhibiting negative values during short time periods of time. This counterintuitive result, found theoretically in other systems and waiting for its experimental confirmation, can be perfectly explained in the framework of the quantum Langevin equation. Its oscillatory behavior is primordially due to the response to the external field while its quantum origin contributes only in its magnitude. The results are compared to those of the continuum limit which exhibits similar behavior.

In this paper, we explore the self-assembly behavior of disk-coil block copolymers (BCPs) confined within a cylinder using molecular dynamics simulations. As functions of the diameter of the confining cylinder and the number of coil beads, concentric lamellar structures are obtained with a different number of alternating disk-rich and coil-rich bilayers. Our paper focuses on the curvature-induced structural behavior in the disk-rich domain of a self-assembled structure, which is investigated by calculating the local density distribution $P\left(r\right)$ and the orientational distribution $G(r,\theta )$. In the inner layers of cylinder-confined disk-coil BCPs, both $P\left(r\right)$ and $G(r,\theta )$ show characteristic asymmetry within a bilayer which is directly contrasted with the bulk and slab-confined disk-coil BCPs. We successfully explain the structural frustration of disks arising from the curved structure due to packing frustration of disks and asymmetric stretching of coils to the regions with different curvatures in a bilayer. Our results are important to understand the self-assembly behavior of BCPs containing a rigid motif in a confined structure, such as a self-assembled structure of bacteriochlorophyll molecules confined by a lipid layer to form a chlorosome, the photosynthetic antennae complex found in nature.

It is well known that glass-forming liquids exhibit a number of anomalous dynamical phenomena, most notably a two-step relaxation in the self-intermediate scattering function and the breakdown of the Stokes-Einstein (SE) relation, as they are cooled toward the glass transition temperature. While these phenomena are generally ascribed to dynamic heterogeneity, specifically to the presence of slow- and fast-moving particles, a quantitative elucidation of the two-step relaxation and the violation of the SE relation in terms of these concepts has not been successful. In this work, we propose a classification of particles according to the rank order of their displacements (from an arbitrarily defined origin of time), and we divide the particles into long-distance (LD), medium-distance, and short-distance (SD) traveling particle groups. Using molecular-dynamics simulation data of the Kob-Andersen model, we show quantitatively that the LD group is responsible for the fast relaxation in the two-step relaxation process in the intermediate scattering function, while the SD group gives rise to the slow ($\alpha$) relaxation. Furthermore, our analysis reveals that ${\tau }_{\alpha }$ is controlled by the SD group, while the ensemble-averaged diffusion coefficient $D$ is controlled by both the LD and SD groups. The combination of these two features provides a natural explanation for the breakdown in the SE relation at low temperature. In addition, we find that the $\alpha$-relaxation time, ${\tau }_{\alpha }$, of the overall system is related to the relaxation time of the LD particles, ${\tau }_{\text{LD}}$, as ${\tau }_{\alpha }={\tau }_{0}exp(\mathrm{\Omega }{\tau }_{\text{LD}}/k_{\text{B}}T)$.

Carefully injected noise can speed the average convergence of Markov chain Monte Carlo (MCMC) estimates and simulated annealing optimization. This includes quantum annealing and the MCMC special case of the Metropolis-Hastings algorithm. MCMC seeks the solution to a computational problem as the equilibrium probability density of a reversible Markov chain. The algorithm must cycle through a long burn-in phase until it reaches equilibrium because the Markov samples are statistically correlated. The special injected noise reduces this burn-in period in MCMC. A related theorem shows that it reduces the cooling time in simulated annealing. Simulations showed that optimal noise gave a 76% speed-up in finding the global minimum in the Schwefel optimization benchmark. The noise-boosted simulations found the global minimum in 99.8% of trials compared with only 95.4% of trials in noiseless simulated annealing. Simulations also showed that the noise boost is robust to accelerated cooling schedules and that noise decreased convergence times by more than 32% under aggressive geometric cooling. Molecular dynamics simulations showed that optimal noise gave a 42% speed-up in finding the minimum potential energy configuration of an eight-argon-atom gas system with a Lennard-Jones 12-6 potential. The annealing speed-up also extends to quantum Monte Carlo implementations of quantum annealing. Noise improved ground-state energy estimates in a 1024-spin simulated quantum annealing simulation by 25.6%. The quantum noise flips spins along a Trotter ring. The noisy MCMC algorithm brings each Markov step closer on average to equilibrium if an inequality holds between two expectations. Gaussian or Cauchy jump probabilities reduce the noise-benefit inequality to a simple quadratic inequality. Simulations show that noise-boosted simulated annealing is more likely than noiseless annealing to sample high probability regions of the search space and to accept solutions that increase the search breadth.

We propose a generalization of stochastic thermodynamics to systems of active particles, which move under the combined influence of stochastic internal self-propulsions (activity) and a heat bath. The main idea is to consider joint trajectories of particles' positions and self-propulsions. It is then possible to exploit formal similarity of an active system and a system consisting of two subsystems interacting with different heat reservoirs and coupled by a nonsymmetric interaction. The resulting thermodynamic description closely follows the standard stochastic thermodynamics. In particular, total entropy production, $\mathrm{\Delta }{s}_{\text{tot}}$, can be decomposed into housekeeping, $\mathrm{\Delta }{s}_{\text{hk}}$, and excess, $\mathrm{\Delta }{s}_{\text{ex}}$, parts. Both $\mathrm{\Delta }{s}_{\text{tot}}$ and $\mathrm{\Delta }{s}_{\text{hk}}$ satisfy fluctuation theorems. The average rate of the steady-state housekeeping entropy production can be related to the violation of the fluctuation-dissipation theorem via a Harada-Sasa relation. The excess entropy production enters into a Hatano-Sasa-like relation, which leads to a generalized Clausius inequality involving the change of the system's entropy and the excess entropy production. Interestingly, although the evolution of particles' self-propulsions is free and uncoupled from that of their positions, nontrivial steady-state correlations between these variables lead to the nonzero excess dissipation in the reservoir coupled to the self-propulsions.

An Ising-type Vicsek model is proposed for collective motion and sudden direction change in a population of self-propelled particles. Particles move on a linear lattice with velocity $+1$ or $-1$ in the one-dimensional model. The probability of the velocity of a particle at the next step is determined by the number difference of the right- and left-moving particles at the present lattice site and its nearest-neighboring sites. A solitary wave appears also in our model similarly to previous models. In some parameter range, the moving direction of the solitary wave sometimes changes rather suddenly, which is like the sudden change of direction of a flock of birds. We study the average reversal time of traveling direction numerically and compare the results with a mean-field theory. The one-dimensional model is generalized to a two-dimensional model. Flip motion of a bandlike soliton is observed in the two-dimensional model.

Irreversible adsorption of objects of different shapes and sizes on Euclidean, fractal, and random lattices is studied. The adsorption process is modeled by using random sequential adsorption algorithm. Objects are adsorbed on one-, two-, and three-dimensional Euclidean lattices, on Sierpinski carpets having dimension $d$ between 1 and 2, and on Erdős-Rényi random graphs. The number of sites is $M=L^d$ for Euclidean and fractal lattices, where $L$ is a characteristic length of the system. In the case of random graphs, such a characteristic length does not exist, and the substrate can be characterized by a fixed set of $M$ vertices (sites) and an average connectivity (or degree) $g$. This paper concentrates on measuring (i) the probability $W_{L\left(M\right)}\left(\theta \right)$ that a lattice composed of $L^d\left(M\right)$ elements reaches a coverage $\theta$ and (ii) the exponent ${\nu }_j$ characterizing the so-called jamming transition. The results obtained for Euclidean, fractal, and random lattices indicate that the quantities derived from the jamming probability $W_{L\left(M\right)}\left(\theta \right)$, such as ${(d{W}_L/d\theta )}_{\text{max}}$ and the inverse of the standard deviation ${\mathrm{\Delta }}_L$, behave asymptotically as $M^{1/2}$. In the case of Euclidean and fractal lattices, where $L$ and $d$ can be defined, the asymptotic behavior can be written as $M^{1/2}=L^{d/2}=L^{1/{\nu }_j}$, with ${\nu }_j=2/d$.

The effect of the Kolmogorov-Arnold-Moser (KAM) theorem in quantum systems is manifested in dividing eigenstates into regular and irregular states. We propose an effective method based on the Wannier basis in phase space to illustrate this division of eigenstates. The quantum kicked-rotor model is used to illustrate this method, which allows us to define the area and effective dimension of each eigenstate to distinguish quantitatively regular and irregular eigenstates. This Wannier basis method also allows us to define the length of a Planck cell in the spectrum that measures how many Planck cells the system will traverse if it starts at the given Planck cell. Moreover, with this Wannier approach, we are able to clarify the distinction between the KAM effect and Anderson localization.

We explore the consequences of introducing higher-order interactions in a geometric complex network of Morris-Lecar neurons. We focus on the regime where traveling synchronization waves are observed from a first-neighbors-based coupling to evaluate the changes induced when higher-order dynamical interactions are included. We observe that the traveling-wave phenomenon gets enhanced by these interactions, allowing the activity to travel further in the system without generating pathological full synchronization states. This scheme could be a step toward a simple phenomenological modelization of neuroglial networks.

Receptor-mediated endocytosis requires that the energy of adhesion overcomes the deformation energy of the plasma membrane. The resulting driving force is balanced by dissipative forces, leading to deterministic dynamical equations. While the shape of the free membrane does not play an important role for tensed and loose membranes, in the intermediate regime it leads to an important energy barrier. Here we show that this barrier is similar to but different from an effective line tension and suggest a simple analytical approximation for it. We then explore the rich dynamics of uptake for particles of different shapes and present the corresponding dynamical state diagrams. We also extend our model to include stochastic fluctuations, which facilitate uptake and lead to corresponding changes in the phase diagrams.

Motivated by the development of cellulose-based functional materials, we investigate the microscopic dynamics of suspensions of cellulose nanocrystals (CNCs) at different ionic strengths, both in the absence and in the presence of AC electric fields and for various temperatures. A concentration of 5 wt % of the CNCs is chosen for which the dispersions are in the full chiral-nematic state at low ionic strengths. Dynamic light scattering is used to characterize the wave vector-dependent decay rates of number-density fluctuations. Contrary to an isotropic suspension, the dispersion relations (the wave vector dependence of the correlation-function decay rates) as obtained by means of depolarized light scattering are found to exhibit anomalous behavior. The dispersion relations, both without and with an external field, exhibit minima at small wave vectors, which is attributed to coupling of translational motion to the orientation of the CNCs, shown in the chiral-nematic state. The location of the minima is found to weakly depend on ionic strength and shifts significantly towards larger wave vectors upon applying an external electric field for sufficiently high ionic strengths. Finally, preliminary results are presented for smaller length-scale density fluctuations (at larger wave vectors) as a function of temperature, revealing the anisotropic mobilities in the chiral-nematic state of CNCs.

Magnetic nanorods rotating in a viscous liquid are very sensitive to any ambient magnetic field. We theoretically predicted and experimentally validated the conditions for two-dimensional synchronous and asynchronous rotation as well as three-dimensional precession and tumbling of nanorods in an ambient field superimposed on a planar rotating magnetic field. We discovered that any ambient field stabilizes the synchronous precession of the nanorod so that the nanorod precession can be completely controlled. This effect opens up different applications of magnetic nanorods as sensors of weak magnetic fields, for microrheology, and generally for magnetic levitation.

Video images are being used with increased frequency in science, supplanting current methods such as light scattering by statistical evaluation of the images. In this study we use light turbidity data due to density-induced refractive index fluctuations to obtain critical amplitudes from image analysis. In order to bring hydrogen (${\mathrm{H}}_2$) very close to its critical point, we place the sample cell under weightlessness using a magnetic levitation device. Images of an ${\mathrm{H}}_2$-filled cell are taken near its critical temperature of 33 K by illuminating the cell with three different filters. We fit the turbidity data to a theoretical expression that allows us to estimate the critical amplitudes of isothermal compressibility and fluctuation correlation length. The values of isothermal compressibility and correlation length obtained from turbidity fitting are compared against literature values. Our data analysis shows a large sensitivity of the fitting parameters to the refractive index value and to even minute density deviations from criticality.

Economic studies have shown that there are two types of regulation schemes which can be considered as a vital part of today's global economy: self-regulation enforced by self-regulation organizations to govern industry practices and government regulation which is considered as another scheme to sustain corporate adherence. An outstanding problem of particular interest is to understand quantitatively the role of these regulation schemes in evolutionary dynamics. Typically, punishment usually occurs for enforcement of regulations. Taking into account both types of punishments to influence the regulations, we develop a game model where six evolutionary situations with corresponding combinations of strategies are considered. Furthermore, a semianalytical method is developed to allow us to give accurate estimations of the boundaries between the phases of full defection and nondefection. We find that, associated with the evolutionary dynamics, for an infinite well-mixed population, the mix of both punishments performs better than one punishment alone in promoting public cooperation, but for a networked population the cooperator-driven punishment turns out to be a better choice. We also reveal the monotonic facilitating effects of the synergy effect, punishment fine, and feedback sensitivity on the public cooperation for an infinite well-mixed population. Conversely, for a networked population an optimal intermediate range of feedback sensitivity is needed to best promote punishers' populations. Overall, a networked structure is overall more favorable for punishers and further for public cooperation, because of both network reciprocity and mutualism between punishers and cooperators who do not punish defectors. We provide physical understandings of the observed phenomena, through a detailed statistical analysis of frequencies of different strategies and spatial pattern formations in different evolution situations. These results provide valuable insights into how to select and enforce suitable regulation measures to let public cooperation remain prevalent, which has potential implications not only for self-regulation, but also for other topics in economics and social science.

A numerical method is described for producing stationary solutions of the Vlasov-Poisson system describing a relativistic charged-particle beam in a constant focusing accelerator channel, confined transversely by a general (linear or nonlinear) focusing potential. The method utilizes a variant of the spectral Galerkin algorithm to solve a nonlinear partial differential equation (PDE) in two degrees of freedom for the beam space charge potential in equilibrium. Numerical convergence with an increasing number of computed spectral modes is investigated for several benchmark problems. Preservation of the stationary phase space density is verified using a strongly nonlinear focusing channel based on the Integrable Optics Test Accelerator at Fermi National Accelerator Laboratory.

Dynamical universality classes are distinguished by their dynamical exponent $z$ and unique scaling functions encoding space-time asymmetry for, e.g., slow-relaxation modes or the distribution of time-integrated currents. So far the universality class of the Nagel-Schreckenberg (NaSch) model, which is a paradigmatic model for traffic flow on highways, was not known. Only the special case $v_{\text{max}}=1$, where the model corresponds to the totally asymmetric simple exclusion process, is known to belong to the superdiffusive Kardar-Parisi-Zhang (KPZ) class with $z=3/2$. In this paper, we show that the NaSch model also belongs to the KPZ class for general maximum velocities $v_{\text{max}}>1$. Using nonlinear fluctuating hydrodynamics theory we calculate the nonuniversal coefficients, fixing the exact asymptotic solutions for the dynamical structure function and the distribution of time-integrated currents. The results of large-scale Monte Carlo simulations match the exact asymptotic KPZ solutions without any fitting parameter left. Additionally, we find that nonuniversal early-time effects or the choice of initial conditions might have a strong impact on the numerical determination of the dynamical exponent and therefore lead to inconclusive results. We also show that the universality class is not changed by extending the model to a two-lane NaSch model with lane-changing rules.