Search

Machine-learned force fields have transformed the atomistic modelling of materials by enabling simulations of ab initio quality on unprecedented time and length scales. However, they are currently limited by: (i) the significant computational and human effort that must go into development and validation of potentials for each particular system of interest; and (ii) a general lack of transferability from one chemical system to the next. Here, using the state-of-the-art MACE architecture we introduce a single general-purpose ML model, trained on a public database of 150k inorganic crystals, that is capable of running stable molecular dynamics on molecules and materials. We demonstrate the power of the MACE-MP-0 model - and its qualitative and at times quantitative accuracy - on a diverse set problems in the physical sciences, including the properties of solids, liquids, gases, chemical reactions, interfaces and even the dynamics of a small protein. The model can be applied out of the box and as a starting or "foundation model" for any atomistic system of interest and is thus a step towards democratising the revolution of ML force fields by lowering the barriers to entry., Comment: 119 pages, 63 figures, 37MB PDF

The DArk Matter Particle Explorer (DAMPE) is well suitable for searching for monochromatic and sharp $\gamma$-ray structures in the GeV$-$TeV range thanks to its unprecedented high energy resolution. In this work, we search for $\gamma$-ray line structures using five years of DAMPE data. To improve the sensitivity, we develop two types of dedicated data sets (including the BgoOnly data which is the first time to be used in the data analysis for the calorimeter-based gamma-ray observatories) and adopt the signal-to-noise ratio optimized regions of interest (ROIs) for different DM density profiles. No line signals or candidates are found between 10 and 300 GeV in the Galaxy. The constraints on the velocity-averaged cross section for $\chi\chi \to \gamma\gamma$ and the decay lifetime for $\chi \to \gamma\nu$, both at 95% confidence level, have been calculated and the systematic uncertainties have been taken into account. Comparing to the previous Fermi-LAT results, though DAMPE has an acceptance smaller by a factor of $\sim 10$, similar constraints on the DM parameters are achieved and below 100 GeV the lower limits on the decay lifetime are even stronger by a factor of a few. Our results demonstrate the potential of high-energy-resolution observations on dark matter detection., Comment: 14 pages, 8 figures. Update the content to keep up with the published version

The measurement of Born cross section of $e^+e^-\to ZH$ process is one of the major goals of the future Circular Electron Positron Collider, which may reach a precision of 0.5%. Such unprecedented precision must be guaranteed by both theoretical and experimental sides, such as the calculations of high order corrections, the knowledge of the $\sigma^{\mathrm{B}}(e^+e^-\to ZH)$ line shape. The uncertainty of radiative correction factor due to lack of knowledge of the $\sigma^{\mathrm{B}}(e^+e^-\to ZH)$ line shape is evaluated in this work. Therefore, the dedicated data-taking schemes are proposed in order to precisely calculate the ISR correction factor., Comment: 9 pages, 5 figures

The discovery of the Higgs boson with its mass around 125 GeV by the ATLAS and CMS Collaborations marked the beginning of a new era in high energy physics. The Higgs boson will be the subject of extensive studies of the ongoing LHC program. At the same time, lepton collider based Higgs factories have been proposed as a possible next step beyond the LHC, with its main goal to precisely measure the properties of the Higgs boson and probe potential new physics associated with the Higgs boson. The Circular Electron Positron Collider~(CEPC) is one of such proposed Higgs factories. The CEPC is an $e^+e^-$ circular collider proposed by and to be hosted in China. Located in a tunnel of approximately 100~km in circumference, it will operate at a center-of-mass energy of 240~GeV as the Higgs factory. In this paper, we present the first estimates on the precision of the Higgs boson property measurements achievable at the CEPC and discuss implications of these measurements., Comment: 46 pages, 37 figures

The framework of statistical inference has been successfully used to detect the meso-scale structures in complex networks, such as community structure, core-periphery (CP) structure. The main principle is that the stochastic block model (SBM) is used to fit the observed network and the learnt parameters indicate the group assignment, in which the parameters of model are often calculated via an expectation-maximization (EM) algorithm and a belief propagation (BP) algorithm is implemented to calculate the decomposition itself. In the derivation process of the BP algorithm, some approximations were made by omitting the effects of node's neighbors, the approximations do not hold if networks are dense or some nodes holding large degrees. As a result, for example, the BP algorithm cannot well detect CP structure in networks and even yields wrong detection because the nodal degrees in core group are very large. In doing so, we propose an improved BP algorithm to solve the problem in the original BP algorithm without increasing any computational complexity. By comparing the improved BP algorithm with the original BP algorithm on community detection and CP detection, we find that the two algorithms yield the same performance on the community detection when the network is sparse, for the community structure in dense networks or CP structure in networks, our improved BP algorithm is much better and more stable. The improved BP algorithm may help us correctly partition different types of meso-scale structures in networks., Comment: 7 pages, 5 figures

It has been recognized that many complex dynamical systems in the real world require a description in terms of multiplex networks, where a set of common, mutually connected nodes belong to distinct network layers and play a different role in each layer. In spite of recent progress towards data based inference of single-layer networks, to reconstruct complex systems with a multiplex structure remains largely open. We articulate a mean-field based maximum likelihood estimation framework to solve this outstanding and challenging problem. We demonstrate the power of the reconstruction framework and characterize its performance using binary time series from a class of prototypical duplex network systems that host two distinct types of spreading dynamics. In addition to validating the framework using synthetic and real-world multiplex networks, we carry out a detailed analysis to elucidate the impacts of structural and dynamical parameters as well as noise on the reconstruction accuracy and robustness., Comment: One main context and a supplementary information

Currently, effective phonons (renormalized or interacting phonons) rather than solitary waves (for short, solitons) are regarded as the energy carriers in nonlinear lattices. In this work, by using the approximate soliton solutions of the corresponding equations of motion and adopting the Boltzmann distribution for these solitons, the average velocities of solitons are obtained and are compared with the sound velocities of energy transfer. Excellent agreements with the numerical results and the predictions of other existing theories are shown in both the symmetric Fermi-Pasta-Ulam-$\beta$ lattices and the asymmetric Fermi-Pasta-Ulam-$\alpha \beta$ lattices. These clearly indicate that solitons are suitable candidates for energy carriers in Fermi-Pasta-Ulam lattices. In addition, the root-mean-square velocity of solitons can be obtained from the effective phonons theory., Comment: 18 pages, 4 figures

Complex networks hosting binary-state dynamics arise in a variety of contexts. In spite of previous works, to fully reconstruct the network structure from observed binary data remains to be challenging. We articulate a statistical inference based approach to this problem. In particular, exploiting the expectation-maximization (EM) algorithm, we develop a method to ascertain the neighbors of any node in the network based solely on binary data, thereby recovering the full topology of the network. A key ingredient of our method is the maximum likelihood estimation of the probabilities associated with actual or non-existent links, and we show that the EM algorithm can distinguish the two kinds of probability values without any ambiguity, insofar as the length of the available binary time series is reasonably long. Our method does not require any {\em a priori} knowledge of the detailed dynamical processes, is parameter free, and is capable of accurate reconstruction even in the presence of noise. We demonstrate the method using combinations of distinct types of binary dynamical processes and network topologies, and provide a physical understanding of the underlying reconstruction mechanism. Our statistical inference based reconstruction method contributes an additional piece to the rapidly expanding "toolbox" of data based reverse engineering of complex networked systems., Comment: 5 figures,accepted by Physical Review E

Recently, the core-periphery (CP) structure of networks as one type of meso-scale structure has received attention. The CP structure is composed of a dense core and a sparse connected periphery. In this paper, we propose an algorithm to detect the CP structure based on the 3-tuple motif, which is inspired by the idea of motif. In this algorithm, we first define a 3-tuple motif by considering the property of nodes, and then a motif adjacency matrix is formed based on the defined motif, finally, the detection of the CP structure is converted to find a cluster that minimizes the smallest motif conductance. Our algorithm can detect different CP structures: including single or multiple CP structure; and local or global CP structures. Results in the synthetic and the empirical networks indicate that the method is efficient and can apply to large-scale networks. More importantly, our algorithm is parameter free, where the core and periphery are detected without the need for any predefined parameters., Comment: 10 pages, 9 figures, submitted

Core-periphery structure and community structure are two typical meso-scale structures in complex networks. Though the community detection has been extensively investigated from different perspectives, the definition and the detection of core-periphery structure have not received much attention. Furthermore, the detection problems of the core-periphery and community structure were separately investigated. In this paper, we develop a unified framework to simultaneously detect core-periphery structure and community structure in complex networks. Moreover, there are several extra advantages of our algorithm: our method can detect not only single but also multiple pairs of core-periphery structures; the overlapping nodes belonging to different communities can be identified; different scales of core-periphery structures can be detected by adjusting the size of core. The good performance of the method has been validated on synthetic and real complex networks. So we provide a basic framework to detect the two typical meso-scale structures: core-periphery structure and community structure., Comment: 10 figures. Accepted by Chaos

Studies on how to model the interplay between diseases and behavioral responses (so-called coupled disease-behavior interaction) have attracted increasing attention. Owing to the lack of obvious clinical evidence of diseases, or the incomplete information related to the disease, the risks of infection cannot be perceived and may lead to inappropriate behavioral responses. Therefore, how to quantitatively analyze the impacts of asymptomatic infection on the interplay between diseases and behavioral responses is of particular importance. In this Letter, under the complex network framework, we study the coupled disease-behavior interaction model by dividing infectious individuals into two states: U-state (without evident clinical symptoms, labelled as U) and I-state (with evident clinical symptoms, labelled as I). A susceptible individual can be infected by U- or I-nodes, however, since the U-nodes cannot be easily observed, susceptible individuals take behavioral responses \emph{only} when they contact I-nodes. The mechanism is considered in the improved Susceptible-Infected-Susceptible (SIS) model and the improved Susceptible-Infected-Recovered (SIR) model, respectively. Then, one of the most concerned problems in spreading dynamics: the epidemic thresholds for the two models are given by two methods. The analytic results \emph{quantitatively} describe the influence of different factors, such as asymptomatic infection, the awareness rate, the network structure, and so forth, on the epidemic thresholds. Moreover, because of the irreversible process of the SIR model, the suppression effect of the improved SIR model is weaker than the improved SIS model., Comment: 6 pages, 3 figures

The spreading of an infectious disease can trigger human behavior responses to the disease, which in turn plays a crucial role on the spreading of epidemic. In this study, to illustrate the impacts of the human behavioral responses, a new class of individuals, $S^F$, is introduced to the classical susceptible-infected-recovered ($SIR$) model. In the model, $S^F$ state represents that susceptible individuals who take self-initiate protective measures to lower the probability of being infected, and a susceptible individual may go to $S^F$ state with a response rate when contacting an infectious neighbor. Via the percolation method, the theoretical formulas for the epidemic threshold as well as the prevalence of epidemic are derived. Our finding indicates that, with the increasing of the response rate, the epidemic threshold is enhanced and the prevalence of epidemic is reduced. The analytical results are also verified by the numerical simulations. In addition, we demonstrate that, because the mean field method neglects the dynamic correlations, a wrong result based on the mean field method is obtained---the epidemic threshold is not related to the response rate, i.e., the additional $S^F$ state has no impact on the epidemic threshold., Comment: 6 figures. accepted by chaos