We consider the statics and dynamics of a stable, mobile three-dimensional (3D) spatiotemporal light bullet in a cubic-quintic nonlinear medium with a focusing cubic nonlinearity above a critical value and any defocusing quintic nonlinearity. The 3D light bullet can propagate with a constant velocity in any direction. Stability of the light bullet under a small perturbation is established numerically. We consider frontal collision between two light bullets with different relative velocities. At large velocities the collision is elastic with the bullets emerge after collision with practically no distortion. At small velocities two bullets coalesce to form a bullet molecule. At a small range of intermediate velocities the localized bullets could form a single entity which expands indefinitely, leading to a destruction of the bullets after collision. The present study is based on an analytic Lagrange variational approximation and a full numerical solution of the 3D nonlinear Schrödinger equation.

The paradigm of layered networks is used to describe many real-world systems, from biological networks to social organizations and transportation systems. While recently there has been much progress in understanding the general properties of multilayer networks, our understanding of how to control such systems remains limited. One fundamental aspect that makes this endeavor challenging is that each layer can operate at a different time scale; thus, we cannot directly apply standard ideas from structural control theory of individual networks. Here we address the problem of controlling multilayer and multi-time-scale networks focusing on two-layer multiplex networks with one-to-one interlayer coupling. We investigate the practically relevant case when the control signal is applied to the nodes of one layer. We develop a theory based on disjoint path covers to determine the minimum number of inputs $\left(N_{\mathrm{i}}\right)$ necessary for full control. We show that if both layers operate on the same time scale, then the network structure of both layers equally affect controllability. In the presence of time-scale separation, controllability is enhanced if the controller interacts with the faster layer: $N_{\mathrm{i}}$ decreases as the time-scale difference increases up to a critical time-scale difference, above which $N_{\mathrm{i}}$ remains constant and is completely determined by the faster layer. We show that the critical time-scale difference is large if layer I is easy and layer II is hard to control in isolation. In contrast, control becomes increasingly difficult if the controller interacts with the layer operating on the slower time scale and increasing time-scale separation leads to increased $N_{\mathrm{i}}$, again up to a critical value, above which $N_{\mathrm{i}}$ still depends on the structure of both layers. This critical value is largely determined by the longest path in the faster layer that does not involve cycles. By identifying the underlying mechanisms that connect time-scale difference and controllability for a simplified model, we provide crucial insight into disentangling how our ability to control real interacting complex systems is affected by a variety of sources of complexity.

When an entangled interpenetrating collection of long flexible polymer chains dispersed in a suitable solvent is cooled to low enough temperatures, thin lamellar crystals form. Remarkably, these lamellae are sectored, with several growth sectors that have differing melting temperatures and growth kinetics, eluding so far an understanding of their origins. We present a theoretical model to explain this six-decade-old challenge by addressing the elasticity of fold surfaces of finite-sized lamella in the presence of disclination-type topological defects arising from anisotropic line tension. Entrapment of a disclination defect in a lamella results in sectors separated by walls, which are soliton solutions of a two-dimensional elliptic sine-Gordon equation. For flat square morphologies, exact results show that sectored squares are more stable than plain squares if the dimensionless anisotropic line tension parameter $\alpha ={\gamma }_{an}/\sqrt{h_4{K}_{\varphi }}$ (${\gamma }_{an}$ = anisotropic line tension, $h_4$ = fold energy parameter, $K_{\varphi }$ = elastic constant for two-dimensional orientational deformation) is above a critical value, which depends on the size of the square.

Tubular crystals, two-dimensional lattices wrapped into cylindrical topologies, arise in many contexts, including botany and biofilaments, and in physical systems such as carbon nanotubes. The geometrical principles of botanical phyllotaxis, describing the spiral packings on cylinders commonly found in nature, have found application in all these systems. Several recent studies have examined defects in tubular crystals associated with crystalline packings that must accommodate a fixed tube radius. Here we study the mechanics of tubular crystals with variable tube radius, with dislocations interposed between regions of different phyllotactic packings. Unbinding and separation of dislocation pairs with equal and opposite Burgers vectors allow the growth of one phyllotactic domain at the expense of another. In particular, glide separation of dislocations offers a low-energy mode for plastic deformations of solid tubes in response to external stresses, reconfiguring the lattice step by step. Through theory and simulation, we examine how the tube's radius and helicity affects, and is in turn altered by, the mechanics of dislocation glide. We also discuss how a sufficiently strong bending rigidity can alter or arrest the deformations of tubes with small radii.

In this work, existing methods for plasma crystal analysis are investigated using artificial and simulated calibration data, which reproduce a multiphase system consisting of fcc, hcp, and bcc crystal data. Disturbances of the artificial data including Gaussian noise, stretching, and randomly missing particles are used to investigate the methods thoroughly. A popular method, called bond order parameter, has been repeatedly criticized as a structure analysis tool and will be improved with the help of a recent development. The method is called the bcc-sensitive Minkowski structure metric. It enhances robustness and consistency, while remaining compatible with previous bond-order-based results. Also, a promising method rooted in the molecular dynamics community is tested, yielding detailed insight of bond-order-specific drawbacks. With this investigation, the state of three-dimensional plasma crystal analysis will be significantly improved.

We report on experiments irradiating isolated plastic spheres with a peak laser intensity of $2\text{–}3\times {10}^{20}{\mathrm{W}\mathrm{cm}}^{-2}$. With a laser focal spot size of $10 \mu \mathrm{m}$ full width half maximum (FWHM) the sphere diameter was varied between 520 nm and $19.3 \mu \mathrm{m}$. Maximum proton energies of $\sim 25$ MeV are achieved for targets matching the focal spot size of $10 \mu \mathrm{m}$ in diameter or being slightly smaller. For smaller spheres the kinetic energy distributions of protons become nonmonotonic, indicating a change in the accelerating mechanism from ambipolar expansion towards a regime dominated by effects caused by Coulomb repulsion of ions. The energy conversion efficiency from laser energy to proton kinetic energy is optimized when the target diameter matches the laser focal spot size with efficiencies reaching the percent level. The change of proton acceleration efficiency with target size can be attributed to the reduced cross-sectional overlap of subfocus targets with the laser. Reported experimental observations are in line with 3D3V particle in cell simulations. They make use of well-defined targets and point out pathways for future applications and experiments.

We consider the diffusion-advection problem in two simple cellular flow models (often invoked as examples of subdiffusive tracer motion) and concentrate on the intermediate time range, in which the tracer motion indeed may show subdiffusion. We perform extensive numerical simulations of the systems under different initial conditions and show that the pure intermediate-time subdiffusion regime is only evident when the particles start at the border between different cells, i.e., at the separatrix, and is less pronounced or absent for other initial conditions. The motion moreover shows quite peculiar aging properties, which are also mirrored in the behavior of the time-averaged mean squared displacement for single trajectories. This kind of behavior is due to the complex motion of tracers trapped inside the cell and is absent in classical models based on continuous-time random walks with no dynamics in the trapped state.

We consider the transfer of lattice wave packets through a tilted discrete breather (TDB) in opposite directions in the discrete nonlinear Schrödinger model with asymmetric defects, which may be realized as a Bose-Einstein condensate trapped in a deep optical lattice, or as optical beams in a waveguide array. A unidirectional transport mode is found, in which the incident wave packets, whose energy belongs to a certain interval between full reflection and full passage regions, pass the TDB only in one direction, while in the absence of the TDB, the same lattice admits bidirectional propagation. The operation of this mode is accurately explained by an analytical consideration of the respective energy barriers. The results suggest that the TDB may emulate the unidirectional propagation of atomic and optical beams in various settings. In the case of the passage of the incident wave packet, the scattering TDB typically shifts by one lattice unit in the direction from which the wave packet arrives, which is an example of the tractor-beam effect, provided by the same system, in addition to the rectification of incident waves.

We consider a homogeneous DNA molecule and investigate the effect of random force applied on the unzipping profile of the molecule. How the critical force varies as a function of the chain length or number of base pairs is the objective of this study. In general, the ratio of the critical forces that is applied on the middle of the chain to that which is applied on one of the ends is two. Our study shows that this ratio depends on the length of the chain. This means that the force which is applied to a point can be experienced by a section of the chain. Beyond a length, the base pairs have no information about the applied force. In the case when the chain length is shorter than this length, this ratio may vary. Only in the case when the chain length exceeds a critical length, this ratio is found to be two. Based on the de Gennes formulation, we developed a method to calculate these forces at zero temperature. The exact results at zero temperature match numerical calculations.

Many of the most interesting complex media are non-Newtonian and exhibit time-dependent behavior of thixotropy and rheopecty. They may also have temporal responses described by power laws. The material behavior is represented by the relaxation modulus and the creep compliance. On the one hand, it is shown that in the special case of a Maxwell model characterized by a linearly time-varying viscosity, the medium's relaxation modulus is a power law which is similar to that of a fractional derivative element often called a springpot. On the other hand, the creep compliance of the time-varying Maxwell model is identified as Lomnitz's logarithmic creep law, making this possibly its first direct derivation. In this way both fractional derivatives and Lomnitz's creep law are linked to time-varying viscosity. A mechanism which yields fractional viscoelasticity and logarithmic creep behavior has therefore been found. Further, as a result of this linking, the curve-fitting parameters involved in the fractional viscoelastic modeling, and the Lomnitz law gain physical interpretation.

As a granular material is compressed, the particles and forces within the system arrange to form complex and heterogeneous collective structures. Force chains are a prime example of such structures, and are thought to constrain bulk properties such as mechanical stability and acoustic transmission. However, capturing and characterizing the evolving nature of the intrinsic inhomogeneity and mesoscale architecture of granular systems can be challenging. A growing body of work has shown that graph theoretic approaches may provide a useful foundation for tackling these problems. Here, we extend the current approaches by utilizing multilayer networks as a framework for directly quantifying the progression of mesoscale architecture in a compressed granular system. We examine a quasi-two-dimensional aggregate of photoelastic disks, subject to biaxial compressions through a series of small, quasistatic steps. Treating particles as network nodes and interparticle forces as network edges, we construct a multilayer network for the system by linking together the series of static force networks that exist at each strain step. We then extract the inherent mesoscale structure from the system by using a generalization of community detection methods to multilayer networks, and we define quantitative measures to characterize the changes in this structure throughout the compression process. We separately consider the network of normal and tangential forces, and find that they display a different progression throughout compression. To test the sensitivity of the network model to particle properties, we examine whether the method can distinguish a subsystem of low-friction particles within a bath of higher-friction particles. We find that this can be achieved by considering the network of tangential forces, and that the community structure is better able to separate the subsystem than a purely local measure of interparticle forces alone. The results discussed throughout this study suggest that these network science techniques may provide a direct way to compare and classify data from systems under different external conditions or with different physical makeup.

In this paper, a comparative study of the lattice Boltzmann (LB) models for the Allen-Cahn (A-C) and Cahn-Hilliard (C-H) equations is conducted. To this end, a new LB model for the A-C equation is first proposed, where the equilibrium distribution function and the source term distribution function are delicately designed to recover the A-C equation correctly. The gradient term in this model can be computed by the nonequilibrium part of the distribution function such that the collision process can be implemented locally. Then a detailed numerical study on several classical problems is performed to give a comparison between the present model for the A-C equation and the recently developed LB model [H. Liang et al., Phys. Rev. E 89, 053320 (2014)] for the C-H equation in terms of tracking the interface of two-phase flow. The results show that the present LB model for the A-C equation is more accurate and more stable, and also has a second-order convergence rate in space, while the convergence rate of the previous LB model for the C-H equation is only about 1.5.

We present an approach for modeling nanoscale wetting and dewetting of textured solid surfaces that exploits recently developed, sophisticated techniques for computing exact long-range dispersive van der Waals (vdW) or (more generally) Casimir forces in arbitrary geometries. We apply these techniques to solve the variational formulation of the Young-Laplace equation and predict the equilibrium shapes of liquid-vacuum interfaces near solid gratings. We show that commonly employed methods of computing vdW interactions based on additive Hamaker or Derjaguin approximations, which neglect important electromagnetic boundary effects, can result in large discrepancies in the shapes and behaviors of liquid surfaces compared to exact methods.

In this paper, we generalize the theory of Brownian motion and the Onsager–Machlup theory of fluctuations for spatially symmetric systems to equilibrium and nonequilibrium steady-state systems with a preferred spatial direction, due to an external force. To do this, we extend the Langevin equation to include a bias, which is introduced by an external force and alters the Gaussian structure of the system's fluctuations. In addition, by solving this extended equation, we provide a physical interpretation for the statistical properties of the fluctuations in these systems. Connections of the extended Langevin equation with the theory of active Brownian motion are discussed as well.

Epilepsy is characterized by substantial network rearrangements leading to spontaneous seizures and little is known on how an epileptogenic focus impacts on neural activity in the contralateral hemisphere. Here, we used a model of unilateral epilepsy induced by injection of the synaptic blocker tetanus neurotoxin (TeNT) in the mouse primary visual cortex (V1). Local field potential (LFP) signals were simultaneously recorded from both hemispheres of each mouse in acute phase (peak of toxin action) and chronic condition (completion of TeNT effects). To characterize the neural electrical activities the corresponding LFP signals were analyzed with several methods of time series analysis. For the epileptic mice, the spectral analysis showed that TeNT determines a power redistribution among the different neurophysiological bands in both acute and chronic phases. Using linear and nonlinear interdependence measures in both time and frequency domains, it was found in the acute phase that TeNT injection promotes a reduction of the interhemispheric coupling for high frequencies ($12\text{–}30$ Hz) and small time lag ($<20$ ms), whereas an increase of the coupling is present for low frequencies ($0.5\text{–}4$ Hz) and long time lag ($>40$ ms). On the other hand, the chronic period is characterized by a partial or complete recovery of the interhemispheric interdependence level. Granger causality test and symbolic transfer entropy indicate a greater driving influence of the TeNT-injected side on activity in the contralateral hemisphere in the chronic phase. Lastly, based on experimental observations, we built a computational model of LFPs to investigate the role of the ipsilateral inhibition and exicitatory interhemispheric connections in the dampening of the interhemispheric coupling. The time evolution of the interhemispheric coupling in such a relevant model of epilepsy has been addressed here.

The growing sluggishness of glass-forming liquids is thought to be accompanied by growing structural order. The nature of such order, however, remains hotly debated. A decade ago, point-to-set (PTS) correlation lengths were proposed as measures of amorphous order in glass formers, but recent results raise doubts as to their generality. Here, we extend the definition of PTS correlations to agnostically capture any type of growing order in liquids, be it local or amorphous. This advance enables the formulation of a clear distinction between slowing down due to conventional critical ordering and that due to glassiness, and provides a unified framework to assess the relative importance of specific local order and generic amorphous order in glass formation.

In this paper general dynamic equations describing the time evolution of isothermal quasi-incompressible multicomponent liquids are derived in the framework of the classical Ginzburg-Landau theory of first order phase transformations. Based on the fundamental equations of continuum mechanics, a general convection-diffusion dynamics is set up first for compressible liquids. The constitutive relations for the diffusion fluxes and the capillary stress are determined in the framework of gradient theories. Next the general definition of incompressibility is given, which is taken into account in the derivation by using the Lagrange multiplier method. To validate the theory, the dynamic equations are solved numerically for the quaternary quasi-incompressible Cahn-Hilliard system. It is demonstrated that variable density (i) has no effect on equilibrium (in case of a suitably constructed free energy functional) and (ii) can influence nonequilibrium pattern formation significantly.

We present a rigorous efficient event-chain Monte Carlo algorithm for long-range interacting particle systems. Using a cell-veto scheme within the factorized Metropolis algorithm, we compute each single-particle move with a fixed number of operations. For slowly decaying potentials such as Coulomb interactions, screening line charges allow us to take into account periodic boundary conditions. We discuss the performance of the cell-veto Monte Carlo algorithm for general inverse-power-law potentials, and illustrate how it provides a new outlook on one of the prominent bottlenecks in large-scale atomistic Monte Carlo simulations.

In this work we study a one-dimensional lattice of Lipkin-Meshkov-Glick models with alternating couplings between nearest-neighbors sites, which resembles the Su-Schrieffer-Heeger model. Typical properties of the underlying models are present in our semiclassical-topological hybrid system, allowing us to investigate an interplay between semiclassical bifurcations at mean-field level and topological phases. Our results show that bifurcations of the energy landscape lead to diverse ordered quantum phases. Furthermore, the study of the quantum fluctuations around the mean-field solution reveals the existence of nontrivial topological phases. These are characterized by the emergence of localized states at the edges of a chain with free open-boundary conditions.

We study spin glass behavior in a random Ising Coulomb antiferromagnet in two and three dimensions using Monte Carlo simulations. In two dimensions, we find a transition at zero temperature with critical exponents consistent with those of the Edwards-Anderson model, though with large uncertainties. In three dimensions, evidence for a finite-temperature transition, as occurs in the Edwards-Anderson model, is rather weak. This may indicate that the sizes are too small to probe the asymptotic critical behavior, or possibly that the universality class is different from that of the Edwards-Anderson model and has a lower critical dimension equal to three.