The interaction between thin elastic films and soft-adhesive foundations has recently gained interest due to technological applications that require control over such objects. Motivated by these applications we investigate the equilibrium configuration of an open cylindrical shell with natural curvature $\kappa$ and bending modulus $B$ that is adhered to soft and adhesive foundation with stiffness $K$. We derive an analytical model that predicts the delamination criterion, i.e., the critical natural curvature, ${\kappa }_{\text{cr}}$, at which delamination first occurs, and the ultimate shape of the shell. While in the case of a rigid foundation, $K\to \infty$, our model recovers the known two-states solution at which the shell either remains completely attached to the substrate or completely detaches from it, on a soft foundation our model predicts the emergence of a new branch of solutions. This branch corresponds to partially adhered shells, where the contact zone between the shell and the substrate is finite and scales as ${\ell }_w\sim {(B/K)}^{1/4}$. In addition, we find that the criterion for delamination depends on the total length of the shell along the curved direction, $L$. While relatively short shells, $L\sim {\ell }_w$, transform continuously between adhered and delaminated solutions, long shells, $L\gg {\ell }_w$, transform discontinuously. Notably, our work provides insights into the detachment phenomena of thin elastic sheets from soft and adhesive foundations.

Nondilute transport occurs routinely in porous medium systems. Experimental observations have revealed effects that seemingly depend upon density, viscosity, velocity, and chemical activity. Macroscale models based upon averaged behavior over many pores have been relied upon to describe such systems to date, which require parametrization of important physical phenomena in material coefficients. To advance fundamental understanding of these complex systems, we examine nondilute transport from a fundamental microscale, or pore-scale, continuum modeling perspective. We approximate the solution of a model based upon the variable-density Navier-Stokes equations and a nondilute species transport equation. Known dependencies of the densities, viscosities, chemical activity, and diffusion for a salt solution on chemical composition are included in the model. Microscale model solutions are averaged to the macroscale and compared with extant experimental observations. Investigation of the effects of various physical phenomena on the microscale velocity distribution and the observed macroscale dispersion are considered using dimensional analysis and constrained simulations. Simulation results are used to explain observed experimental results in light of underlying mechanisms. Conditions under which the various physicochemical effects investigated are important are revealed.

We investigate the forcing strength needed to sustain a flow using linear forcing. A critical Reynolds number $R_c$ is determined, based on the longest wavelength allowed by the system, the forcing strength and the viscosity. A simple model is proposed for the dissipation rate, leading to a closed expression for the kinetic energy of the flow as a function of the Reynolds number. The dissipation model and the prediction for the kinetic energy are assessed using direct numerical simulations and two-point closure integrations. An analysis of the dissipation-rate equation and the triadic structure of the nonlinear transfer allows to refine the model in order to reproduce the low-Reynolds-number asymptotic behavior, where the kinetic energy is proportional to $R-R_c$.

We study the finite-time erasure of a one-bit memory consisting of a one-dimensional double-well potential, with each well encoding a memory macrostate. We focus on setups that provide full control over the form of the potential-energy landscape and derive protocols that minimize the average work needed to erase the bit over a fixed amount of time. We allow for cases where only some of the information encoded in the bit is erased. For systems required to end up in a local-equilibrium state, we calculate the minimum amount of work needed to erase a bit explicitly, in terms of the equilibrium Boltzmann distribution corresponding to the system's initial potential. The minimum work is inversely proportional to the duration of the protocol. The erasure cost may be further reduced by relaxing the requirement for a local-equilibrium final state and allowing for any final distribution compatible with constraints on the probability to be in each memory macrostate. We also derive upper and lower bounds on the erasure cost.

We develop the framework of classical observational entropy, which is a mathematically rigorous and precise framework for nonequilibrium thermodynamics, explicitly defined in terms of a set of observables. Observational entropy can be seen as a generalization of Boltzmann entropy to systems with indeterminate initial conditions, and it describes the knowledge achievable about the system by a macroscopic observer with limited measurement capabilities; it becomes Gibbs entropy in the limit of perfectly fine-grained measurements. This quantity, while previously mentioned in the literature, has been investigated in detail only in the quantum case. We describe this framework reasonably pedagogically, then show that in this framework, certain choices of coarse-graining lead to an entropy that is well-defined out of equilibrium, additive on independent systems, and that grows toward thermodynamic entropy as the system reaches equilibrium, even for systems that are genuinely isolated. Choosing certain macroscopic regions, this dynamical thermodynamic entropy measures how close these regions are to thermal equilibrium. We also show that in the given formalism, the correspondence between classical entropy (defined on classical phase space) and quantum entropy (defined on Hilbert space) becomes surprisingly direct and transparent, while manifesting differences stemming from noncommutativity of coarse-grainings and from nonexistence of a direct classical analog of quantum energy eigenstates.

A novel type of waves is examined in the context of non-Hermitian photonics. We can identify a class of complex guided structures that support localized paraxial solutions whose intensity distribution is exactly the same as the intensity of a corresponding solution in homogeneous media (free or bulk space). In other words, intensity-wise the two solutions are identical and their phase is different by a factor $exp\left[i\theta \right(x,y\left)\right]$. The non-Hermitian potential is determined by the phase $\theta$, as well as the amplitude and phase of the bulk space solution that contributes to the imaginary and real part of the potential, respectively. That way we can connect the plane waves and Gaussian beams of free space to constant-intensity waves and what we call the equal-intensity waves (EI waves) in non-Hermitian media. Such a relation allows us to study three different physical problems: Propagating EI waves inside random media, interface lattice solitons, and moving solitons in photonic waveguide structures with free-space characteristics. The relation of EI waves to unidirectional invisibility and Bohmian photonics is also examined.

We investigate the heterogeneity of outcomes of repeated instances of percolation experiments in complex networks using a message-passing approach to evaluate heterogeneous, node-dependent probabilities of belonging to the giant or percolating cluster, i.e., the set of mutually connected nodes whose size scales linearly with the size of the system. We evaluate these both for large finite single instances and for synthetic networks in the configuration model class in the thermodynamic limit. For the latter, we consider both Erdős-Rényi and scale-free networks as examples of networks with narrow and broad degree distributions, respectively. For real-world networks we use an undirected version of a Gnutella peer-to-peer file-sharing network with $N=62568$ nodes as an example. We derive the theory for multiple instances of both uncorrelated and correlated percolation processes. For the uncorrelated case, we also obtain a closed-form approximation for the large mean degree limit of Erdős-Rényi networks.

We study both numerically and experimentally the breakup of a viscoelastic liquid bridge formed between two parallel electrodes. The polymer solutions and applied voltages are those commonly used in electrospinning and near-field electrospinning. We solve the leaky-dielectric finitely extensible nonlinear elastic-Peterlin (FENE-P) model to describe the dynamical response of the liquid bridge under isothermal conditions. The results show that the surface charge screens the inner electric field perpendicular to the free surface over the entire dynamical process. The liquid bridge deformation produces a normal electric field on the outer side of the free surface that is commensurate with the axial one. The surface conduction does not significantly affect the current intensity in the time interval analyzed in the experiments. The force due to the shear electric stress becomes comparable to both the viscoelastic and surface tension forces in the last stage of the filament. However, it does not alter the elastocapillary balance in the filament. As a consequence, the extensional relaxation times measured from the filament exponential thinning approximately coincides with the stress relaxation time prescribed in the FENE-P model. The above results allow us to interpret correctly the experiments. In the experiments, we measure the filament electrical conductivity and extensional relaxation time for polyethylene oxide (PEO) dissolved in deionized water and in a mixture of water and glycerine. We compare the filament electrical conductivity with the value measured in hydrostatic conditions for the same estimated temperature. Good agreement was found for PEO dissolved in water $+$ glycerine, which indicates that the change in the filament microscopic structure due to the presence of stretched polymeric chains does not significantly alter the ion mobility in the stretching direction. Significant deviations are found for PEO dissolved in deionized water. These deviations may be attributed to the heat transferred to the ambient, which is neglected in the calculation of the filament temperature. We measure the extensional relaxation time from the images acquired during the filament thinning. The relaxation times obtained in the first stage of the exponential thinning hardly depend on the applied voltage. Little but measurable influence of the applied voltage is found in the last phase of the filament thinning.

Magnetic reconnection in a relativistic electron magnetization regime was observed in a laboratory plasma produced by a high-intensity, large energy, picoseconds laser pulse. Magnetic reconnection conditions realized with a laser-driven several kilotesla magnetic field is comparable to that in the accretion disk corona of black hole systems, i.e., Cygnus X-1. We observed particle energy distributions of reconnection outflow jets, which possess a power-law component in a high-energy range. The hardness of the observed spectra could explain the hard-state x-ray emission from accreting black hole systems.

In this paper we study thermodynamic properties of uniform electron gas (UEG) over wide density and temperature range, using the improved fermionic-path-integral Monte Carlo (FPIMC) method. This method demonstrates a significant reduction of the “fermionic sign problem,” which takes place in standard path-integral Monte Carlo simulations of degenerate fermionic systems. We introduce three basic improvements. The first one is the improved treatment of exchange interaction, achieved by the proper change of variables in the path-integral measure. The second improvement is the inclusion of long-range Coulomb effects into an angle-averaged effective potential, as proposed by Yakub and Ronchi [J. Chem. Phys. 119, 11556 (2003)]. The third improvement is the angle-averaging of an exchange determinant, describing the fermionic exchange interaction not only between particles in the main Monte Carlo cell, but also with electrons in the nearest periodic images. The FPIMC shows very good agreement with analytical data for ideal Fermi gas. For strongly coupled UEG under warm dense matter conditions we compare our total and exchange-correlation energy results with other Monte Carlo approaches.

We study a lattice-gas model of penetrable particles on a square-lattice substrate with same-site and nearest-neighbor interactions. Penetrability implies that the number of particles occupying a single lattice site is unlimited and the model itself is intended as a simple representation of penetrable particles encountered in realistic soft-matter systems. Our specific focus is on a binary mixture, where particles of the same species repel and those of the opposite species attract each other. As a consequence of penetrability and the unlimited occupation of each site, the system exhibits thermodynamic collapse, which in simulations is manifested by an emergence of extremely dense clusters scattered throughout the system with energy of a cluster $E\propto -n^2$, where $n$ is the number of particles in a cluster. After transforming a particle system into a spin system, in the large density limit the Hamiltonian recovers a simple harmonic form, resulting in the discrete Gaussian model used in the past to model the roughening transition of interfaces. For finite densities, due to the presence of a nonharmonic term, the system is approximated using a variational Gaussian model.

Classical arcsine law states that the fraction of occupation time on the positive or the negative side in Brownian motion does not converge to a constant but converges in distribution to the arcsine distribution. Here we consider how a preparation of the system affects the arcsine law, i.e., aging of the arcsine law. We derive an aging distributional theorem for occupation time statistics in Brownian motion, where the ratio of time when measurements start to the measurement time plays an important role in determining the shape of the distribution. Furthermore, we show that this result can be generalized as an aging distributional limit theorem in renewal processes.

A small tagged particle immersed in a fluid exhibits Brownian motion and diffuses on a long timescale. Meanwhile, on a short timescale, the dynamics of the tagged particle cannot be simply described by the usual generalized Langevin equation with Gaussian noise, since the number of collisions between the tagged particle and fluid particles is rather small. On such a timescale, we should explicitly consider individual collision events between the tagged particle and the surrounding fluid particles. In this study we analyze the short-time dynamics of a tagged particle in an ideal gas, where we do not have static or hydrodynamic correlations between fluid particles. We perform event-driven hard-sphere simulations and show that the short-time dynamics of the tagged particle is correlated even under such an idealized situation. Namely, the velocity autocorrelation function becomes negative when the tagged particle is relatively light and the fluid density is relatively high. This result can be attributed to the dynamical correlation between collision events. To investigate the physical mechanism which causes the dynamical correlation, we analyze the correlation between successive collision events. We find that the tagged particle can collide with the same ideal-gas particle several times and such collisions cause a strong dynamical correlation for the velocity.

A paradigmatic framework to study the phenomenon of spontaneous collective synchronization is provided by the Kuramoto model comprising a large collection of phase oscillators of distributed frequencies that are globally coupled through the sine of their phase differences. We study here a variation of the model by including nearest-neighbor interactions on a one-dimensional lattice. While the mean-field interaction resulting from the global coupling favors global synchrony, the nearest-neighbor interaction may have cooperative or competitive effects depending on the sign and the magnitude of the nearest-neighbor coupling. For unimodal and symmetric frequency distributions, we demonstrate that as a result, the model in the stationary state exhibits in contrast to the usual Kuramoto model both continuous and first-order transitions between synchronized and incoherent phases, with the transition lines meeting at a tricritical point. Our results are based on numerical integration of the dynamics as well as an approximate theory involving appropriate averaging of fluctuations in the stationary state.

Three-dimensional (3D) flow structures around a wall-mounted short cylinder of height-to-diameter ratio 1 were instantaneously measured by a high-resolution tomographic particle image velocimetry (Tomo-PIV) at Reynolds number of 10 720 in a water tunnel. 3D velocity fields, 3D vorticity, the $Q$ criterion, the rear separation region, and the characteristic of arch type vortex and tip vortices were first discussed. We found a strong 3D W-type arch vortex behind the short cylinder, which was originated by the interaction between upwash and downwash flows. This W-type arch vortex was reshaped to the M-shaped arch vortex downstream. It indicated that the head shape of the arch vortex structure depended on the aspect ratio of the cylinder. The large-scale streamwise vortices were originated by the downwash and upwash flows near the center location of W-type arch vortex. Then the 3D orthogonal wavelet multiresolution technique was developed to analyze instantaneous 3D velocity fields of Tomo-PIV in order to clarify 3D multiscale wake flow structures. The W-type shape arch vortex was extracted in the time-averaged intermediate-scale structure, while an M-shaped arch vortex was identified in the time-averaged large-scale structure. The tip vortices distributed in the time-averaged large- and intermediate-scale structures. The instantaneous intermediate-scale upwash vortices played an essential role in producing W-type head of arch structure. It was also observed that strong small-scale vortices appeared in the shear layer or near the bottom plate and most of them were contained in the intermediate-scale structures.

In a steady state, the linear scaling laws are confirmed between the intensity characteristics of electroconvective (EC) vortex (including the vortex height and electroosmotic slip velocity) and the applied voltage for the nonshear EC flow with finite vortex height near permselective membranes. This finding in the nonshear EC flow is different from the shear EC flow [Kwak et al., Phys. Rev. Lett. 110, 114501 (2013)] and indicates that the local concentration gradient has a significant improvement in the analysis of slip velocity. Further, our study reveals that the EC vortex is mainly driven by the second peak effect of the Coulomb thrust in the extended space-charge layer, and the linear scaling law exhibited by the Coulomb thrust is an essential reason for the linear scaling laws of vortex intensity. The scaling laws proposed in this paper are supported by our direct numerical simulation data and previous experimental observations [Rubinstein et al., Phys. Rev. Lett. 101, 236101 (2008)].

Large-amplitude fluctuations of the solar wind magnetic field can scatter energetic ions. One of the main contributions to these fluctuations is provided by solar wind discontinuities, i.e., rapid rotations of the magnetic field. This study shows that the internal configuration of such discontinuities plays a crucial role in energetic ion scattering in pitch angles. Kinetic-scale discontinuities accomplish very fast ion pitch-angle scattering. The main mechanism of such pitch-angle scattering is the adiabatic invariant destruction due to separatrix crossings in the phase space. We demonstrate that efficiency of this scattering does not depend on the magnetic field component across the discontinuity surface, i.e., both rotational and almost tangential discontinuities scatter energetic ions with the same efficiency. We also examine how the strong scattering effect depends on the deviations of the discontinuity magnetic field from the force-free one.

Unsupervised learning requiring only raw data is not only a fundamental function of the cerebral cortex, but also a foundation for a next generation of artificial neural networks. However, a unified theoretical framework to treat sensory inputs, synapses, and neural activity together is still lacking. The computational obstacle originates from the discrete nature of synapses, and complex interactions among these three essential elements of learning. Here, we propose a variational mean-field theory in which the distribution of synaptic weights is considered. The unsupervised learning can then be decomposed into two intertwined steps: A maximization step is carried out as a gradient ascent of the lower bound on the data log-likelihood, in which the synaptic weight distribution is determined by updating variational parameters, and an expectation step is carried out as a message passing procedure on an equivalent or dual neural network whose parameter is specified by the variational parameters of the weight distribution. Therefore, our framework provides insights on how data (or sensory inputs), synapses, and neural activities interact with each other to achieve the goal of extracting statistical regularities in sensory inputs. This variational framework is verified in restricted Boltzmann machines with planted synaptic weights and handwritten-digits learning.

In this paper we analyze the entropy and entropy production of a nonisolated quantum system described within the quantum Brownian motion framework. This is a very general and paradigmatic framework for describing nonisolated quantum systems and can be used in any kind of coupling regime. We start by considering the application of von Neumann entropy to an arbitrarily damped quantum system making use of its reduced density operator. We argue that this application is formally valid and develop a path-integral method to evaluate that quantity analytically. We apply this technique to a harmonic oscillator in contact with a heat bath and obtain an exact form for its entropy. Then we study the entropy production of this system and enlighten important characteristics of its thermodynamical behavior on the pure quantum realm and also address their transition to the classical limit.

The general set of nonlocal $M$-component nonlinear Schrödinger (nonlocal $M$-NLS) equations obeying the $\mathcal{PT}$-symmetry and featuring focusing, defocusing, and mixed (focusing-defocusing) nonlinearities that has applications in nonlinear optics settings, is considered. First, the multisoliton solutions of this set of nonlocal $M$-NLS equations in the presence and in the absence of a background, particularly a periodic line wave background, are constructed. Then, we study the intriguing soliton collision dynamics as well as the interesting positon solutions on zero background and on a periodic line wave background. In particular, we reveal the fascinating shape-changing collision behavior similar to that of in the Manakov system but with fewer soliton parameters in the present setting. The standard elastic soliton collision also occurs for particular parameter choices. More interestingly, we show the possibility of such elastic soliton collisions even for defocusing nonlinearities. Furthermore, for the nonlocal $M$-NLS equations, the dependence of the collision characteristics on the speed of the solitons is analyzed. In the presence of a periodic line wave background, we notice that the soliton amplitude can be enhanced significantly, even for infinitesimal amplitude of the periodic line waves. In addition to these solutions, by considering the long-wavelength limit of the obtained soliton solutions with proper parameter constraints, higher-order positon solutions of the nonlocal $M$-NLS equations are derived. The background of periodic line waves also influences the wave profiles and amplitudes of the positons. Specifically, the positon amplitude can not only be enhanced but also be suppressed on the periodic line wave background of infinitesimal amplitude.