In 1911, J\"uttner proposed the generalization, for a relativistic gas, of the Maxwell-Boltzmann distribution of velocities. Here we want to discuss, among others, J\"uttner probability density function (PDF). Both the velocity space and, consequently, the momentum space are not flat in special relativity. The velocity space corresponds to the Lobachevsky one, which has a negative curvature. This curvature induces a specific power for the Lorentz factor in the PDF, affecting J\"uttner normalization constant in one, two, and three dimensions. Furthermore, J\"uttner distribution, written in terms of a more convenient variable, the rapidity, presents a curvature change at the origin at sufficiently high energy, which does not agree with our computational dynamics simulations of a relativistic gas. However, in one dimension, the rapidity satisfies a simple additivity law. This allows us to obtain, through the Central Limit Theorem, a new, Lorentz-invariant, PDF whose curvature at the origin does not change for any energy value and which agrees with our computational dynamics simulations data. Also, we perform extensive first-principle simulations of a one-dimensional relativistic gas constituted by light and heavy particles.

Evaldo M. F. Curado Carlos E. Cedeño Ivano Damião Soares Constantino Tsallis

Modeling of the NICER X-ray waveform of the pulsar PSR J0030+0451, aimed to constrain the neutron star mass and radius, has inferred surface hot-spots (the magnetic polar caps) that imply significantly non-dipolar magnetic fields. To this end, we investigate magnetic field configurations that comprise offset dipole plus quadrupole components using static vacuum field and force-free global magnetosphere models. Taking into account the compactness and observer angle values provided by Miller et al. (2019) and Riley et al. (2019), we compute geodesics from the observer plane to the polar caps to compute the resulting X-ray light curve. We explore, through Markov chain Monte Carlo techniques, the detailed magnetic field configurations that can reproduce the observed X-ray light curve and have discovered degeneracies, i.e., diverse field configurations, which can provide sufficient descriptions to the NICER X-ray waveforms. Having obtained the force-free field structures, we then compute the corresponding synchronous gamma-ray light curves following Kalapotharakos et al. (2014) these we compare to those obtained by Fermi-LAT, to provide models consistent with both the X-ray and the gamma-ray data, thereby restricting further the multipole field parameters. An essential aspect of this approach is the proper computation of the relative phase between the synchronous X- and gamma-ray light curves. We conclude with a discussion of the broader implications of our study.

Constantinos Kalapotharakos Zorawar Wadiasingh Alice K. Harding Demosthenes Kazanas

Earlier work has shown that reusing experience from prior motion planning problems can improve the efficiency of similar, future motion planning queries. However, for robots with many degrees-of-freedom, these methods exhibit poor generalization across different environments and often require large datasets that are impractical to gather. We present SPARK and FLAME , two experience-based frameworks for sampling-based planning applicable to complex manipulators in 3 D environments. Both combine samplers associated with features from a workspace decomposition into a global biased sampling distribution. SPARK decomposes the environment based on exact geometry while FLAME is more general, and uses an octree-based decomposition obtained from sensor data. We demonstrate the effectiveness of SPARK and FLAME on a Fetch robot tasked with challenging pick-and-place manipulation problems. Our approaches can be trained incrementally and significantly improve performance with only a handful of examples, generalizing better over diverse tasks and environments as compared to prior approaches.

Constantinos Chamzas Zachary Kingston Carlos Quintero-Peña Anshumali Shrivastava Lydia E. Kavraki

Critiques presented by P. Grassberger in Phys. Rev. Lett. 95, 140601 (2005) are addressed.

Constantino Tsallis

The long-time dynamics of the critical contact process which is brought suddenly out of an uncorrelated initial state undergoes ageing in close analogy with quenched magnetic systems. In particular, we show through Monte Carlo simulations in one and two dimensions and through mean-field theory that time-translation invariance is broken and that dynamical scaling holds. We find that the autocorrelation and autoresponse exponents lambda_{Gamma} and lambda_R are equal but, in contrast to systems relaxing to equilibrium, the ageing exponents a and b are distinct. A recent proposal to define a non-equilibrium temperature through the short-time limit of the fluctuation-dissipation ratio is therefore not applicable.

Jose J. Ramasco Malte Henkel Maria Augusta Santos Constantino A. da Silva Santos

In order to explain the absence of hysteresis in ferromagnetic p-type (Cd,Mn)Te quantum wells (QWs), spin dynamics was previously investigated by Monte Carlo simulations combining the Metropolis algorithm with the determination of hole eigenfunctions at each Monte Carlo sweep. Short-range antiferromagnetic superexchange interactions between Mn spins - which compete with the hole-mediated long-range ferromagnetic coupling - were found to accelerate magnetization dynamics if the the layer containing Mn spins is wider than the vertical range of the hole wave function. Employing this approach it is shown here that appreciate magnitudes of remanence and coercivity can be obtained if Mn ions are introduced to the quantum well in a delta-like fashion.

C. Simserides A. Lipinska K. N. Trohidou T. Dietl

Markov Chain Monte Carlo (MCMC) techniques are now widely used for cosmological parameter estimation. Chains are generated to sample the posterior probability distribution obtained following the Bayesian approach. An important issue is how to optimize the efficiency of such sampling and how to diagnose whether a finite-length chain has adequately sampled the underlying posterior probability distribution. We show how the power spectrum of a single such finite chain may be used as a convergence diagnostic by means of a fitting function, and discuss strategies for optimizing the distribution for the proposed steps. The methods developed are applied to current CMB and LSS data interpreted using both a pure adiabatic cosmological model and a mixed adiabatic/isocurvature cosmological model including possible correlations between modes. For the latter application, because of the increased dimensionality and the presence of degeneracies, the need for tuning MCMC methods for maximum efficiency becomes particularly acute.

Joanna Dunkley Martin Bucher Pedro G. Ferreira Kavilan Moodley Constantinos Skordis

Recent measurements of the temperature fluctuations of the cosmic microwave background (CMB) radiation from the WMAP satellite provide indication of a non-Gaussian behavior. Although the observed feature is small, it is detectable and analyzable. Indeed, the temperature distribution P^{CMB}(Delta T) of these data can be quite well fitted by the anomalous probability distribution emerging within nonextensive statistical mechanics, based on the entropy S_q = k (1 - \int dx [P(x)]^q)/(q - 1) (where in the limit case q -> 1 we obtain the Boltzmann-Gibbs entropy S_1 = - k \int dx P(x) ln[P(x)]). For the CMB frequencies analysed, \nu= 40.7, 60.8, and 93.5 GHz, P^{CMB}(Delta T) is well described by P_q(Delta T) \propto 1/[1 + (q-1) B(\nu) (Delta T)^2]^{1/(q-1)}, with q = 1.04 \pm 0.01, the strongest non-Gaussian contribution coming from the South-East sector of the celestial sphere. Moreover, Monte Carlo simulations exclude, at the 99% confidence level, P_1(Delta T) \propto e^{- B(\nu) (Delta T)^2} to fit the three-year WMAP data.

Armando Bernui Constantino Tsallis Thyrso Villela

We address an iterative procedure that can be used to detect coarse-grained hyperbolic unstable equilibria (saddle points) of microscopic simulators when no equations at the macroscopic level are available. The scheme is based on the concept of coarse timestepping [Kevrekidis et al., 2003] incorporating an adaptive mechanism based on the chord method allowing the location of coarse-grained saddle points directly. Ultimately, it can be used in a consecutive manner to trace the coarse-grained open-loop saddle-node bifurcation diagrams. We illustrate the procedure through two indicatively examples including (i) a kinetic Monte Carlo simulation (kMC) of simple surface catalytic reactions and (ii) a simple agent-based model, a financial caricature which is used to simulate the dynamics of buying and selling of a large population of interacting individuals in the presence of mimesis. Both models exhibit coarse-grained regular turning points which give rise to branches of saddle points.

A. C. Tsoumanis C. I. Siettos

This article reports a measurement of the production cross section of prompt isolated photon pairs in proton-antiproton collisions at \sqrt{s} = 1.96 TeV using the CDF II detector at the Fermilab Tevatron collider. The data correspond to an integrated luminosity of 5.36/fb. The cross section is presented as a function of kinematic variables sensitive to the reaction mechanisms. The results are compared with three perturbative QCD calculations: (1) a leading order parton shower Monte Carlo, (2) a fixed next-to-leading order calculation and (3) a next-to-leading order/next-to-next-to-leading-log resummed calculation. The comparisons show that, within their known limitations, all calculations predict the main features of the data, but no calculation adequately describes all aspects of the data.

CDF Collaboration T. Aaltonen