We study the thermalization dynamics of one-dimensional diatomic lattices (which represents the simplest system possessing multi-branch phonons), exemplified by the famous Fermi-Pasta-Ulam-Tsingou (FPUT)-$\beta$ and the Toda models. Here we focus on how the system relaxes to the equilibrium state when part of highest-frequency optical modes are initially excited, which is called the anti-FPUT problem comparing with the original FPUT problem (low frequency excitations of the monatomic lattice). It is shown numerically that the final thermalization time $T_{\rm eq}$ of the diatomic FPUT-$\beta$ chain depends on whether its acoustic modes are thermalized, whereas the $T_{\rm eq}$ of the diatomic Toda chain depends on the optical ones; in addition, the metastable state of both models have different energy distributions and lifetimes. Despite these differences, in the near-integrable region, the $T_{\rm eq}$ of both models still follows the same scaling law, i.e., $T_{\rm eq}$ is inversely proportional to the square of the perturbation strength. Finally, comparisons of the thermalization behavior between different models under various initial conditions are briefly summarized.

This study analyzes the Collatz map through nonlinear dynamics. By embedding integers in Sharkovsky's ordering, we show that odd initial values suffice for full dynamical characterization. We introduce ``direction phases'' to partition iterations into upward and downward phases, and derive a recursive function family parameterized by upward phase counts. Consequently, a logarithmic scaling law between iteration steps and initial values is revealed, demonstrating finite-time convergence to the period-three orbit. Moreover, we establish the equivalence of the Collatz map to a binary shift map, whose ergodicity guarantees universal convergence to attractors, providing additional support for convergence. Furthermore, we identify that basins of attraction follow power-law distributions and find that odd numbers classified by upward phases follow Gamma statistics. These results offer valuable insights into the dynamics of discrete systems and their connections to number theory.

Nonlinear normal modes are periodic orbits that survive in nonlinear chains, whose instability plays a crucial role in the dynamics of many-body Hamiltonian systems toward thermalization. Here we focus on how the stability of nonlinear modes depends on the perturbation strength and the system size to observe whether they have the same behavior in different models. To this end, as illustrating examples, the instability dynamics of the ${N}/{2}$ mode in both the Fermi-Pasta-Ulam-Tsingou (FPUT) -$\alpha$ and -$\beta$ chains under fixed boundary conditions are studied systematically. Applying the Floquet theory, we show that for both models the stability time $T$ as a function of the perturbation strength $\lambda$ follows the same behavior; i.e., $T\propto(\lambda-\lambda_c)^{-\frac{1}{2}}$, where $\lambda_c$ is the instability threshold. The dependence of $\lambda_c$ on $N$ is also obtained. The results of $T$ and $\lambda_c$ agree well with those obtained by the direct molecular dynamics simulations. Finally, the effect of instability dynamics on the thermalization properties of a system is briefly discussed.

Exploring meaningful structural regularities embedded in networks is a key to understanding and analyzing the structure and function of a network. The node-attribute information can help improve such understanding and analysis. However, most of the existing methods focus on detecting traditional communities, i.e., groupings of nodes with dense internal connections and sparse external ones. In this paper, based on the connectivity behavior of nodes and homogeneity of attributes, we propose a principle model (named GNAN), which can generate both topology information and attribute information. The new model can detect not only community structure, but also a range of other types of structure in networks, such as bipartite structure, core-periphery structure, and their mixture structure, which are collectively referred to as generalized structure. The proposed model that combines topological information and node-attribute information can detect communities more accurately than the model that only uses topology information. The dependency between attributes and communities can be automatically learned by our model and thus we can ignore the attributes that do not contain useful information. The model parameters are inferred by using the expectation-maximization algorithm. And a case study is provided to show the ability of our model in the semantic interpretability of communities. Experiments on both synthetic and real-world networks show that the new model is competitive with other state-of-the-art models.

The symmetry of interparticle interaction potential (IIP) has a crucial influence on the thermodynamic and transport properties of solids. Here we focus on the effect of the asymmetry of IIP on thermalization properties. In general, asymmetry and nonlinearity interweave with each other. To explore the effects of asymmetry and nonlinearity on thermalization separately, we here introduce an asymmetric harmonic (AH) model, whose IIP has the only asymmetry but no nonlinearity in the sense that the frequency of a system is independent of input energy. Through extensive numerical simulations, a power-law relationship between the thermalization time $T_{\rm eq}$ and the perturbation strength (here the degree of asymmetry) is still confirmed, yet a larger exponent is observed instead of the previously found inverse square law in the thermodynamic limit. Then the quartic (symmetric) nonlinearity is added into the AH model, and the thermalization behavior under the coexistence of asymmetry and nonlinearity is systematically studied. It is found that Matthiessen's rule gives a better estimation of the $T_{\rm eq}$ of the system. This work further suggests that the asymmetry plays a distinctive role in the properties of relaxation of a system.

In this letter, a multi-wave quasi-resonance framework is established to analyze energy diffusion in classical lattices, uncovering that it is fundamentally determined by the characteristics of eigenmodes. Namely, based on the presence and the absence of extended modes, lattices fall into two universality classes with qualitatively different thermalization behavior. In particular, we find that while the one with extended modes can be thermalized under arbitrarily weak perturbations in the thermodynamic limit, the other class can be thermalized only when perturbations exceed a certain threshold, revealing for the first time the possibility that a lattice cannot be thermalized, violating the hypothesis of statistical mechanics. Our study addresses conclusively the renowned Fermi-Pasta-Ulam-Tsingou problem for large systems under weak perturbations, underscoring the pivotal roles of both extended and localized modes in facilitating energy diffusion and thermalization processes.

Nonlinear normal modes are periodic orbits that survive in nonlinear many-body Hamiltonian systems, and their instability is crucial for relaxation dynamics. Here, we study the instability process of the $\pi/3$-mode in the Fermi-Pasta-Ulam-Tsingou-$\alpha$ chain with fixed boundary conditions. We find that three types of bifurcations -- period-doubling, tangent, and Hopf -- coexist in this system, each driving instability at specific reduced wave-number $\tilde{k}$. Our analysis reveals a universal scaling law for the instability time $\mathcal{T} \propto (\lambda - \lambda_{\rm c})^{-1/2}$, independent of bifurcation types and models, where the critical perturbation strength $\lambda_{\rm c}$ scales as $\lambda_{\rm c} \propto (\tilde{k} - \tilde{k}_{\rm c})$, with$\tilde{k}_{\rm c}$ varying across bifurcations. We also observe a double instability phenomenon for certain system sizes, meaning that larger perturbations do not always lead to faster thermalization. These results provide new insights into the relaxation and thermalization dynamics in many-body systems.

On-site potentials are ubiquitous in physical systems and strongly influence their heat transport and energy localization. These potentials will inevitably affect the dynamical properties of $q$-breathers (QBs), defined as periodic orbits exponentially localized in normal mode space. By integrating on-site terms into the Fermi-Pasta-Ulam-Tsingou-$\beta$ system, this work utilizes numerical simulations and Floquet analysis to systematically explore the influence of on-site potentials on QB stability. For most QBs, except those at the phonon band edges, the instability is primarily governed by parametric resonance, and effectively described by coupled Mathieu equations. This approach provides a theoretical expression for the instability thresholds, which aligns well with numerical results. We demonstrate that the instability thresholds can be controlled through the strength of on-site potentials, and for a strong enough quadratic on-site potential, the QBs are always stable. Furthermore, the instability threshold is highly sensitive to the seed mode, in stark contrast to systems without on-site potentials. In addition, the instability phase diagrams exhibit joint interplay between different terms in the Hamiltonian, such as the quadratic on-site and quartic inter-site interaction terms, in regulating the QB dynamics. These findings offer valuable insights into QB stability and the manipulation of localized excitations in diverse physical systems with on-site potentials.

The Land$\acute{e}$ $g$ factors of Ba$^+$ are very important in high-precision measurement physics. The wave functions, energy levels, and Land$\acute{e}$ $g$ factors for the $6s$ $^{2}S_{1/2}$ and $5d$ $^{2}D_{3/2,5/2}$ states of Ba$^{+}$ ions were calculated using the multi-configuration Dirac-Hartree-Fock (MCDHF) method and the Model-QED method. The contributions of the electron correlation effects and quantum electrodynamics (QED) effects were discussed in detail. The transition energies are in excellent agreement with the experimental results, with differences of approximately 5 cm$^{-1}$. The presently calculated $g$ factor of 2.0024905(16) for the $6S_{1/2}$ agrees very well with the available experimental and theoretical results, with a difference at a level of 10$^{-6}$. For the $5D_{3/2, 5/2}$ states, the present results of 0.7993961(126) and 1.2003942(190) agree with the experimental results of 0.7993278(3) [\textcolor{blue}{Phys. Rev. A 54, 1199(1996)}] and 1.20036739(14) [\textcolor{blue}{Phys. Rev. Lett. 124, 193001 (2020)}] very well, with differences at the level of 10$^{-5}$.

Whether and how a system approaches equilibrium is central in nonequilibrium statistical physics, crucial to understanding thermalization and transport. Bogoliubov's three-stage (initial, kinetic, and hydrodynamic) evolution hypothesis offers a qualitative framework, but quantitative progress has focused on near-integrable systems like dilute gases. In this work, we investigate the relaxation dynamics of a one-dimensional diatomic hard-point (DHP) gas, presenting a phase diagram that characterizes relaxation behavior across the full parameter space, from near-integrable to far-from-integrable regimes. We analyze thermalization (local energy relaxation in nonequilibrium states) and identify three universal dynamical regimes: (i) In the near-integrable regime, kinetic processes dominate, local energy relaxation decays exponentially, and the thermalization time $\tau$ scales as $\tau \propto \delta^{-2}$. (ii) In the far-from-integrable regime, hydrodynamic effects dominate, energy relaxation decays power-law, and thermalization time scales linearly with system size $N$. (iii) In the intermediate regime, the Bogoliubov phase emerges, characterized by the transition from kinetic to hydrodynamic relaxation. The phase diagram also shows that hydrodynamic behavior can emerge in small systems when sufficiently far from the integrable regime, challenging the view that such effects occur only in large systems. In the thermodynamic limit, the system's relaxation depends on the order in which the limits ($N \to \infty$ or $\delta \to 0$) are taken. We then analyze heat transport (decay of heat-current fluctuations in equilibrium), demonstrating its consistency with thermalization, leading to a unified theoretical description of thermalization and transport. Our approach provides a pathway for studying relaxation dynamics in many-body systems, including quantum systems.

In this paper, we bring forward a completely perturbed nonconvex Schatten $p$-minimization to address a model of completely perturbed low-rank matrix recovery. The paper that based on the restricted isometry property generalizes the investigation to a complete perturbation model thinking over not only noise but also perturbation, gives the restricted isometry property condition that guarantees the recovery of low-rank matrix and the corresponding reconstruction error bound. In particular, the analysis of the result reveals that in the case that $p$ decreases $0$ and a>1 for the complete perturbation and low-rank matrix, the condition is the optimal sufficient condition \delta_{2r}<1 \cite{Recht et al 2010}. The numerical experiments are conducted to show better performance, and provides outperformance of the nonconvex Schatten $p$-minimization method comparing with the convex nuclear norm minimization approach in the completely perturbed scenario.

Compared to the the classical first-order Grünwald-Letnikov formula at time $t_{k+1} (\textmd{or}\, t_{k})$, we firstly propose a second-order numerical approximate scheme for discretizing the Riemann-Liouvile derivative at time $t_{k+\frac{1}{2}}$, which is very suitable for constructing the Crank-Niclson technique applied to the time-fractional differential equations. The established formula has the following form $$\begin{array}{lll} \displaystyle \,_{\mathrm{RL}}{\mathrm{D}}_{0,t}^{\alpha}u\left(t\right)\left|\right._{t=t_{k+\frac{1}{2}}}= \tau^{-\alpha}\sum\limits_{\ell=0}^{k} \varpi_{\ell}^{(\alpha)}u\left(t_k-\ell\tau\right) +\mathcal{O}(\tau^2),\,\,k=0,1,\ldots, \alpha\in(0,1), \end{array}$$ where the coefficients $\varpi_{\ell}^{(\alpha)}$ $(\ell=0,1,\ldots,k)$ can be determined via the following generating function $$\begin{array}{lll} \displaystyle G(z)=\left(\frac{3\alpha+1}{2\alpha}-\frac{2\alpha+1}{\alpha}z+\frac{\alpha+1}{2\alpha}z^2\right)^{\alpha},\;|z|<1. \end{array}$$ Applying this formula to the time fractional Cable equations with Riemann-liouville derivative in one or two space dimensions. Then the high-order compact finite difference schemes are obtained. The solvability, stability and convergence with orders $\mathcal{O}(\tau^2+h^4)$ and $\mathcal{O}(\tau^2+h_x^4+h_y^4)$ are shown, where $\tau$ is the temporal stepsize and $h$, $h_x$, $h_y$ are the spatial stepsizes, respectively. Finally, numerical experiments are provided to support the theoretical analysis.

This paper concentrates on the recovery of block-sparse signals, which is not only sparse but also nonzero elements are arrayed into some blocks (clusters) rather than being arbitrary distributed all over the vector, from linear measurements. We establish high-order sufficient conditions based on block RIP to ensure the exact recovery of every block $s$-sparse signal in the noiseless case via mixed $l_2/l_p$ minimization method, and the stable and robust recovery in the case that signals are not accurately block-sparse in the presence of noise. Additionally, a lower bound on necessary number of random Gaussian measurements is gained for the condition to be true with overwhelming probability. Furthermore, the numerical experiments conducted demonstrate the performance of the proposed algorithm.

Pressure plays a vital role in changing the transport properties of matter. To understand this phenomenon at a microscopic level, we here focus on a more fundamental problem, i.e., how pressure affects the thermalization properties of solids. As illustrating examples, we study the thermalization behavior of the monatomic chain and the mass-disordered chain of Fermi-Pasta-Ulam-Tsingou-$\beta$ under different strains in the thermodynamic limit. It is found that the pressure-induced change in nonintegrability results in qualitatively different thermalization processes for the two kinds of chains. However, for both cases, the thermalization time follows the same law -- it is inversely proportional to the square of the nonintegrability strength. This result suggests that pressure can significantly change the integrability of a system, which provides a new perspective for understanding the pressure-dependent thermal transport behavior.