Efficient Monte Carlo algorithm and high-precision results for percolation

We present a new Monte Carlo algorithm for studying site or bond percolation on any lattice. The algorithm allows us to calculate quantities such as the cluster size distribution or spanning probability over the entire range of site or bond occupation probabilities from zero to one in a single run which takes an amount of time scaling linearly with the number of sites on the lattice. We use our algorithm to determine that the percolation transition occurs at occupation probability 0.59274621(13) for site percolation on the square lattice and to provide clear numerical confirmation of the conjectured 4/3-power stretched-exponential tails in the spanning probability functions.Comment: 8 pages, including 3 postscript figures, minor corrections in this version, plus updated figures for the position of the percolation transitio

Transistors

Contains reports on eight research projects.Lincoln Laboratory under Contract AF19(122)-45

Universality of finite-size corrections to the number of critical percolation clusters

Monte-Carlo simulations on a variety of 2d percolating systems at criticality suggest that the excess number of clusters in finite systems over the bulk value of nc is a universal quantity, dependent upon the system shape but independent of the lattice and percolation type. Values of nc are found to high accuracy, and for bond percolation confirm the theoretical predictions of Temperley and Lieb, and Baxter, Temperley, and Ashley, which we have evaluated explicitly in terms of simple algebraic numbers. Predictions for the fluctuations are also verified for the first time.Comment: 13 pages, 2 figs., Latex, submitted to Phys. Rev. Let

A fast Monte Carlo algorithm for site or bond percolation

We describe in detail a new and highly efficient algorithm for studying site or bond percolation on any lattice. The algorithm can measure an observable quantity in a percolation system for all values of the site or bond occupation probability from zero to one in an amount of time which scales linearly with the size of the system. We demonstrate our algorithm by using it to investigate a number of issues in percolation theory, including the position of the percolation transition for site percolation on the square lattice, the stretched exponential behavior of spanning probabilities away from the critical point, and the size of the giant component for site percolation on random graphs.Comment: 17 pages, 13 figures. Corrections and some additional material in this version. Accompanying material can be found on the web at http://www.santafe.edu/~mark/percolation

Transistors

Contains reports on six research projects

Exact results at the 2-D percolation point

We derive exact expressions for the excess number of clusters b and the excess cumulants b_n of a related quantity at the 2-D percolation point. High-accuracy computer simulations are in accord with our predictions. b is a finite-size correction to the Temperley-Lieb or Baxter-Temperley-Ashley formula for the number of clusters per site n_c in the infinite system limit; the bn correct bulk cumulants. b and b_n are universal, and thus depend only on the system's shape. Higher-order corrections show no apparent dependence on fractional powers of the system size.Comment: 12 pages, 2 figures, LaTeX, submitted to Physical Review Letter

Comparing apples and oranges: assessment of the relative video quality in the presence of different types of distortions

<p>Abstract</p> <p>Video quality assessment is essential for the performance analysis of visual communication applications. Objective metrics can be used for estimating the relative quality differences, but they typically give reliable results only if the compared videos contain similar types of quality distortion. However, video compression typically produces different kinds of visual artifacts than transmission errors. In this article, we focus on a novel subjective quality assessment method that is suitable for comparing different types of quality distortions. The proposed method has been used to evaluate how well different objective quality metrics estimate the relative subjective quality levels for content with different types of quality distortions. Our conclusion is that none of the studied objective metrics works reliably for assessing the co-impact of compression artifacts and transmission errors on the subjective quality. Nevertheless, we have observed that the objective metrics' tendency to either over- or underestimate the perceived impact of transmission errors has a high correlation with the spatial and temporal activity levels of the content. Therefore, our results can be useful for improving the performance of objective metrics in the presence of both source and channel distortions.</p

Multi-step self-guided pathways for shape-changing metamaterials

Multi-step pathways, constituted of a sequence of reconfigurations, are central to a wide variety of natural and man-made systems. Such pathways autonomously execute in self-guided processes such as protein folding and self-assembly, but require external control in macroscopic mechanical systems, provided by, e.g., actuators in robotics or manual folding in origami. Here we introduce shape-changing mechanical metamaterials, that exhibit self-guided multi-step pathways in response to global uniform compression. Their design combines strongly nonlinear mechanical elements with a multimodal architecture that allows for a sequence of topological reconfigurations, i.e., modifications of the topology caused by the formation of internal self-contacts. We realized such metamaterials by digital manufacturing, and show that the pathway and final configuration can be controlled by rational design of the nonlinear mechanical elements. We furthermore demonstrate that self-contacts suppress pathway errors. Finally, we demonstrate how hierarchical architectures allow to extend the number of distinct reconfiguration steps. Our work establishes general principles for designing mechanical pathways, opening new avenues for self-folding media, pluripotent materials, and pliable devices in, e.g., stretchable electronics and soft robotics.Comment: 16 pages, 3 main figures, 10 extended data figures. See https://youtu.be/8m1QfkMFL0I for an explanatory vide

ϕ\phi meson production in dd++Au collisions at sNN=200\sqrt{s_{_{NN}}}=200 GeV

The PHENIX experiment has measured $\phi$ meson production in $d$$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV using the dimuon and dielectron decay channels. The $\phi$ meson is measured in the forward (backward) $d$-going (Au-going) direction, $1.2<y<2.2$ ($-2.2<y<-1.2$) in the transverse-momentum ($p_T$) range from 1--7 GeV/$c$, and at midrapidity $|y|<0.35$ in the $p_T$ range below 7 GeV/$c$. The $\phi$ meson invariant yields and nuclear-modification factors as a function of $p_T$, rapidity, and centrality are reported. An enhancement of $\phi$ meson production is observed in the Au-going direction, while suppression is seen in the $d$-going direction, and no modification is observed at midrapidity relative to the yield in $p$$+$$p$ collisions scaled by the number of binary collisions. Similar behavior was previously observed for inclusive charged hadrons and open heavy flavor indicating similar cold-nuclear-matter effects.Comment: 484 authors, 16 pages, 12 figures, 6 tables. v1 is the version accepted for publication in Phys. Rev. C. Data tables for the points plotted in the figures are given in the paper itsel

Measurements of double-helicity asymmetries in inclusive J/ψJ/\psi production in longitudinally polarized p+pp+p collisions at s=510\sqrt{s}=510 GeV

We report the double helicity asymmetry, $A_{LL}^{J/\psi}$, in inclusive $J/\psi$ production at forward rapidity as a function of transverse momentum $p_T$ and rapidity $|y|$. The data analyzed were taken during $\sqrt{s}=510$ GeV longitudinally polarized $p$$+$$p$ collisions at the Relativistic Heavy Ion Collider (RHIC) in the 2013 run using the PHENIX detector. At this collision energy, $J/\psi$ particles are predominantly produced through gluon-gluon scatterings, thus $A_{LL}^{J/\psi}$ is sensitive to the gluon polarization inside the proton. We measured $A_{LL}^{J/\psi}$ by detecting the decay daughter muon pairs $\mu^+ \mu^-$ within the PHENIX muon spectrometers in the rapidity range $1.2<|y|<2.2$. In this kinematic range, we measured the $A_{LL}^{J/\psi}$ to be $0.012 \pm 0.010$~(stat)~$\pm$~$0.003$(syst). The $A_{LL}^{J/\psi}$ can be expressed to be proportional to the product of the gluon polarization distributions at two distinct ranges of Bjorken $x$: one at moderate range $x \approx 0.05$ where recent RHIC data of jet and $\pi^0$ double helicity spin asymmetries have shown evidence for significant gluon polarization, and the other one covering the poorly known small-$x$ region $x \approx 2\times 10^{-3}$. Thus our new results could be used to further constrain the gluon polarization for $x< 0.05$.Comment: 335 authors, 10 pages, 4 figures, 3 tables, 2013 data. Version accepted for publication by Phys. Rev. D. Plain text data tables for the points plotted in figures for this and previous PHENIX publications are (or will be) publicly available at http://www.phenix.bnl.gov/papers.htm