We consider the Peano curve separating a spanning tree from its dual spanning tree on an embedded planar graph, where the tree and dual tree are weighted by $y$ to the number of active edges, and “active” is in the sense of the Tutte polynomial. When the graph is a portion of the square grid approximating a simply connected domain, it is known ($y=1$ and $y=1+\sqrt{2}$) or believed ($1<y<3$) that the Peano curve converges to a space-filling ${SLE}_{\kappa }$ loop, where $y=1-2cos(4\pi /\kappa )$, corresponding to $4<\kappa \le 8$. We argue that the same should hold for $0\le y<1$, which corresponds to $8<\kappa \le 12$.

Relaxation dynamics of complex quantum systems with strong interactions towards the steady state is a fundamental problem in statistical mechanics. The steady state of subsystems weakly interacting with their environment is described by the canonical ensemble which assumes the probability distribution for energy to be of the Boltzmann form. The emergence of this probability distribution is ensured by the detailed balance of the transitions induced by the interaction with the environment. Here we consider relaxation of an open correlated quantum system brought into contact with a reservoir in the vacuum state. We refer to such a system as emissive since particles irreversibly evaporate into the vacuum. The steady state of the system is a statistical mixture of the stable eigenstates. We found that, despite the absence of the detailed balance, the stationary probability distribution over these eigenstates is of the Boltzmann form in each $N$-particle sector. A quantum statistical ensemble corresponding to the steady state is characterized by different temperatures in the different sectors, in contrast to the Gibbs ensemble. We investigate the transition rates between the eigenstates to understand the emergence of the Boltzmann distribution and find their exponential dependence on the transition energy. We argue that this property of transition rates is generic for a wide class of emissive quantum many-body systems.

We introduce a generalized $d$-dimensional Fermi-Pasta-Ulam model in the presence of long-range interactions, and perform a first-principle study of its chaos for $d=1,2,3$ through large-scale numerical simulations. The nonlinear interaction is assumed to decay algebraically as $d_{ij}^{-\alpha }$ $(\alpha \ge 0)$, $\left\{d_{ij}\right\}$ being the distances between $N$ oscillator sites. Starting from random initial conditions we compute the maximal Lyapunov exponent ${\lambda }_{max}$ as a function of $N$. Our $N\gg 1$ results strongly indicate that ${\lambda }_{max}$ remains constant and positive for $\alpha /d>1$ (implying strong chaos, mixing, and ergodicity), and that it vanishes like $N^{-\kappa }$ for $0\le \alpha /d<1$ (thus approaching weak chaos and opening the possibility of breakdown of ergodicity). The suitably rescaled exponent $\kappa$ exhibits universal scaling, namely that $(d+2)\kappa$ depends only on $\alpha /d$ and, when $\alpha /d$ increases from zero to unity, it monotonically decreases from unity to zero, remaining so for all $\alpha /d>1$. The value $\alpha /d=1$ can therefore be seen as a critical point separating the ergodic regime from the anomalous one, $\kappa$ playing a role analogous to that of an order parameter. This scaling law is consistent with Boltzmann-Gibbs statistics for $\alpha /d>1$, and possibly with $q$ statistics for $0\le \alpha /d<1$.

When a droplet with a higher density falls in a miscible solution, the droplet deforms and breaks up. The instability of a vortex ring, formed by droplet deformation during the falling process, causes the breakup. To determine the origin of the instability, the wavelengths and thicknesses of the vortex rings are investigated at the time when the instability occurs. The experimental results are almost in agreement with the calculated results for the Rayleigh-Taylor instability using the thickness of a higher-density solution. Furthermore, we performed simulations considering the torus shapes and circulations of the vortex ring. The simulations provided patterns similar to those observed experimentally for the breakup process, and showed that the circulations suppress the instability of the vortex ring. These results imply that the Rayleigh-Taylor instability plays a dominant role in the instability of vortex rings.

It is commonly assumed in percolation theories that at most one percolating cluster can exist in a network. We show that several coexisting percolating clusters (CPCs) can emerge in networks due to limited mixing, i.e., a finite and sufficiently small number of interlinks between network modules. We develop an approach called modular message passing (MMP) to describe and verify these observations. We demonstrate that the appearance of CPCs is an important source of inaccuracy in previously introduced percolation theories, such as the message passing (MP) approach, which is a state-of-the-art theory based on the belief propagation method. Moreover, we show that the MMP theory improves significantly over the predictions of MP for percolation on synthetic networks with limited mixing and also on several real-world networks. These findings have important implications for understanding the robustness of networks and in quantifying epidemic outbreaks in the susceptible-infected-recovered (SIR) model of disease spread.

We present analytical results for the distribution of shortest path lengths between random pairs of nodes in configuration model networks. The results, which are based on recursion equations, are shown to be in good agreement with numerical simulations for networks with degenerate, binomial, and power-law degree distributions. The mean, mode, and variance of the distribution of shortest path lengths are also evaluated. These results provide expressions for central measures and dispersion measures of the distribution of shortest path lengths in terms of moments of the degree distribution, illuminating the connection between the two distributions.

We analyze the connectivity of an $M$-layer network over a common set of nodes that are active only in a fraction of the layers. Each layer is assumed to be a subgraph (of an underlying connectivity graph $G$) induced by each node being active in any given layer with probability $q$. The $M$-layer network is formed by aggregating the edges over all $M$ layers. We show that when $q$ exceeds a threshold $q_c\left(M\right)$, a giant connected component appears in the $M$-layer network—thereby enabling far-away users to connect using “bridge” nodes that are active in multiple network layers—even though the individual layers may only have small disconnected islands of connectivity. We show that $q_c\left(M\right)\lesssim \sqrt{-ln(1-p_c)}/\sqrt{M}$, where $p_c$ is the bond percolation threshold of $G$, and $q_c\left(1\right)\equiv q_c$ is its site-percolation threshold. We find $q_c\left(M\right)$ exactly for when $G$ is a large random network with an arbitrary node-degree distribution. We find $q_c\left(M\right)$ numerically for various regular lattices and find an exact lower bound for the kagome lattice. Finally, we find an intriguingly close connection between this multilayer percolation model and the well-studied problem of site-bond percolation in the sense that both models provide a smooth transition between the traditional site- and bond-percolation models. Using this connection, we translate known analytical approximations of the site-bond critical region, which are functions only of $p_c$ and $q_c$ of the respective lattice, to excellent general approximations of the multilayer connectivity threshold $q_c\left(M\right)$.

We experimentally investigated the splashing of dense suspension droplets impacting a solid surface, extending prior work to the regime where the viscosity of the suspending liquid becomes a significant parameter. The overall behavior can be described by a combination of two trends. The first one is that the splashing becomes favored when the kinetic energy of individual particles at the surface of a droplet overcomes the confinement produced by surface tension. This is expressed by a particle-based Weber number ${\mathrm{We}}_p$. The second is that splashing is suppressed by increasing the viscosity of the solvent. This is expressed by the Stokes number St, which influences the effective coefficient of restitution of colliding particles. We developed a phase diagram where the splashing onset is delineated as a function of both ${\mathrm{We}}_p$ and St. A surprising result occurs at very small Stokes number, where not only splashing is suppressed but also plastic deformation of the droplet. This leads to a situation where droplets can bounce back after impact, an observation we are able to reproduce using discrete particle numerical simulations that take into account viscous interaction between particles and elastic energy.

The force network of a granular assembly, defined by the contact network and the corresponding contact forces, carries valuable information about the state of the packing. Simple analysis of these networks based on the distribution of force strengths is rather insensitive to the changes in preparation protocols or to the types of particles. In this and the companion paper [Kondic et al., Phys. Rev. E 93, 062903 (2016)], we consider two-dimensional simulations of tapped systems built from frictional disks and pentagons, and study the structure of the force networks of granular packings by considering network's topology as force thresholds are varied. We show that the number of clusters and loops observed in the force networks as a function of the force threshold are markedly different for disks and pentagons if the tangential contact forces are considered, whereas they are surprisingly similar for the network defined by the normal forces. In particular, the results indicate that, overall, the force network is more heterogeneous for disks than for pentagons. Such differences in network properties are expected to lead to different macroscale response of the considered systems, despite the fact that averaged measures (such as force probability density function) do not show any obvious differences. Additionally, we show that the states obtained by tapping with different intensities that display similar packing fraction are difficult to distinguish based on simple topological invariants.

In the companion paper [Pugnaloni et al., Phys. Rev. E 93, 062902 (2016)], we use classical measures based on force probability density functions (PDFs), as well as Betti numbers (quantifying the number of components, related to force chains, and loops), to describe the force networks in tapped systems of disks and pentagons. In the present work, we focus on the use of persistence analysis, which allows us to describe these networks in much more detail. This approach allows us not only to describe but also to quantify the differences between the force networks in different realizations of a system, in different parts of the considered domain, or in different systems. We show that persistence analysis clearly distinguishes the systems that are very difficult or impossible to differentiate using other means. One important finding is that the differences in force networks between disks and pentagons are most apparent when loops are considered: the quantities describing properties of the loops may differ significantly even if other measures (properties of components, Betti numbers, force PDFs, or the stress tensor) do not distinguish clearly or at all the investigated systems.

We use an explicitly isochoric two-temperature theory to analyze recent x-ray laser scattering data for aluminum in the ultra-fast-matter (UFM) regime up to 6 eV. The observed surprisingly low conductivities are explained by including strong electron-ion scattering effects using the phase shifts calculated via the neutral-pseudo-atom model. The difference between the static conductivity for UFM-Al and equilibrium aluminum in the warm-dense matter state is clearly brought out by comparisons with available density-fucntional+molecular-dynamics simulations. Thus the applicability of the Mermin model to UFM is questioned. The static and dynamic conductivity, collision frequency, and the plasmon line shape, evaluated within the simplest Born approximation for UFM aluminum, are in good agreement with experiment.

We present an approximation for calculating the equation of state (EOS) of warm and hot dense matter that is built on the previously published pseudoatom molecular dynamics (PAMD) model of dense plasmas [Starrett et al., Phys. Rev. E 91, 013104 (2015)]. While the EOS calculation with PAMD was previously limited to orbital-free density functional theory (DFT), the new approximation presented here allows a Kohn-Sham DFT treatment of the electrons. The resulting EOS thus includes a quantum mechanical treatment of the electrons with a self-consistent model of the ionic structure, while remaining tractable at high temperatures. The method is validated by comparisons with pressures from ab initio simulations of Be, Al, Si, and Fe. The EOS in the Thomas-Fermi approximation shows remarkable thermodynamic consistency over a wide range of temperatures for aluminum. We calculate the principal Hugoniots of aluminum and silicon up to 500 eV. We find that the ionic structure of the plasma has a modest effect that peaks at temperatures of a few eV and that the features arising from the electronic structure agree well with ab initio simulations.

Nanomechanical electron shuttles can work as ratchets for radio-frequency rectification. We develop a full stochastic model of coupled shuttles, where the mechanical motion of nanopillars and the incoherent electronic tunneling are modeled by a Markov chain. In particular, the interaction of their randomness is taken into account, so that a linear master equation is constructed. Numerical solutions from our fast approximate method and analytical derivation reveal the symmetry breaking, which results in the direct current observed in earlier measurements [Phys. Rev. Lett. 105, 067204 (2010)]. Additionally, the method can facilitate device simulation of more complex designs such as shuttle arrays.

Mean-field theory (MFT) is one of the main available tools for analytical calculations entailed in investigations regarding many-body systems. Recently, there has been a surge of interest in ameliorating this kind of method, mainly with the aim of incorporating geometric and correlation properties of these systems. The correlated cluster MFT (CCMFT) is an improvement that succeeded quite well in doing that for classical spin systems. Nevertheless, even the CCMFT presents some deficiencies when applied to quantum systems. In this article, we address this issue by proposing the quantum CCMFT (QCCMFT), which, in contrast to its former approach, uses general quantum states in its self-consistent mean-field equations. We apply the introduced QCCMFT to the transverse Ising model in honeycomb, square, and simple cubic lattices and obtain fairly good results both for the Curie temperature of thermal phase transition and for the critical field of quantum phase transition. Actually, our results match those obtained via exact solutions, series expansions or Monte Carlo simulations.

The fermionic and Majorana edge mode dynamics of various topological systems are compared, after a sudden global quench of the Hamiltonian parameters takes place. Attention is focused on the regimes where the survival probability of an edge state has oscillations either due to critical or off-critical quenches. The nature of the wave functions and the overlaps between the eigenstates of different points in parameter space determine the various types of behaviors, and the distinction due to the Majorana nature of the excitations plays a lesser role. Performing a sequence of quenches, it is shown that the edge states, including Majorana modes, may be switched off and on. Also, the generation of Majoranas due to quenching from a trivial phase is discussed.

We establish a stochastic thermodynamics for a Fermionic level driven by a time-dependent force and interacting with initially thermalized levels playing the role of a reservoir. The driving induces consecutive avoided crossings between system and reservoir levels described within Landau-Zener theory. We derive the resulting system dynamics and thermodynamics and identify energy, work, heat, entropy, and dissipation. Our theory perfectly reproduces the numerically exact quantum work statistics obtained using a two point measurements approach of the total energy and provides an explicit expression for the dissipation in terms of diabatic transitions.

We take a central spin model (CSM), consisting of a one-dimensional environmental Ising spin chain and a single qubit connected globally to all the spins of the environment, to study the excess energy (EE) of the environment and the logarithm of decoherence factor namely, generalized fidelity susceptibility per site (GFSS), associated with the qubit under a periodic driving of the transverse field term of environment across its critical point using the Floquet theory. The coupling to the qubit, prepared in a pure state, with the transverse field of the spin chain yields two sets of EE corresponding to the two species of Floquet operators. In the limit of weak coupling, we derive an approximated expression of GFSS after an infinite number of driving period which can successfully estimate the low- and intermediate-frequency behavior of GFSS obtained numerically with a large number of time periods. Our main focus is to analytically investigate the effect of system-environment coupling strength on the EEs and GFSS and relate the behavior of GFSS to EEs as a function of frequency by plausible analytical arguments. We explicitly show that the low-frequency beatinglike pattern of GFSS is an outcome of two frequencies, causing the oscillations in the two branches of EEs, that are dependent on the coupling strength. In the intermediate frequency regime, dip structure observed in GFSS can be justified by the resonance peaks of EEs at those coupling parameter-dependent frequencies; high-frequency saturation behavior of EEs and GFSS are controlled by the same static Hamiltonian and the associated saturation values are related to the coupling strength.

For a long time, extreme events happening in complex systems, such as financial markets, earthquakes, and neurological networks, were thought to follow power-law size distributions. More recently, evidence suggests that in many systems the largest and rarest events differ from the other ones. They are dragon kings, outliers that make the distribution deviate from a power law in the tail. Understanding the processes of formation of extreme events and what circumstances lead to dragon kings or to a power-law distribution is an open question and it is a very important one to assess whether extreme events will occur too often in a specific system. In the particular system studied in this paper, we show that the rate of occurrence of dragon kings is controlled by the value of a parameter. The system under study here is composed of two nearly identical chaotic oscillators which fail to remain in a permanently synchronized state when coupled. We analyze the statistics of the desynchronization events in this specific example of two coupled chaotic electronic circuits and find that modifying a parameter associated to the local instability responsible for the loss of synchronization reduces the occurrence of dragon kings, while preserving the power-law distribution of small- to intermediate-size events with the same scaling exponent. Our results support the hypothesis that the dragon kings are caused by local instabilities in the phase space.

The binding energies of complexes of DNA nucleobase pairs are evaluated using quantum mechanical calculations at the level of dispersion corrected density functional theory. We begin with Watson-Crick base pairs of singlets, duplets, and triplets and calculate their binding energies. At a second step, mismatches are incorporated into the Watson-Crick complexes in order to evaluate the variation in the binding energy with respect to the canonical Watson-Crick pairs. A linear variation of this binding energy with the degree of mismatching is observed. The binding energies for the duplets and triplets containing mismatches are further compared to the energies of the respective singlets in order to assess the degree of collectivity in these complexes. This study also suggests that mismatches do not considerably affect the energetics of canonical base pairs. Our work is highly relevant to the recognition process in DNA promoted through the RecA protein and suggests a clear distinction between recognition in singlets, and recognition in duplets or triplets. Our work assesses the importance of collectivity in the homology recognition of DNA.

We performed parallel-tempering Monte Carlo simulations to investigate the formation and stability of helical tertiary structures for flexible and semiflexible polymers, employing a generic coarse-grained model. Structural conformations exhibit helical order with tertiary ordering into single helices, multiple helical segments organized into bundles, and disorganized helical arrangements. For both bending-restrained semiflexible and bending-unrestrained flexible helical polymers, the stability of the structural phases is discussed systematically by means of hyperphase diagrams parametrized by suitable order parameters, temperature, and torsion strength. This exploration lends insight into the restricted flexibility of biological polymers such as double-stranded DNA and proteins.