Measurement of Indeterminacy in Packings of Perfectly Rigid Disks

Static packings of perfectly rigid particles are investigated theoretically and numerically. The problem of finding the contact forces in such packings is formulated mathematically. Letting the values of the contact forces define a vector in a high-dimensional space enable us to show that the set of all possible contact forces is convex, facilitating its numerical exploration. It is also found that the boundary of the set is connected with the presence of sliding contacts, suggesting that a stable packing should not have more than 2M-3N sliding contacts in two dimensions, where M is the number of contacts and N is the number of particles. These results were used to analyze packings generated in different ways by either molecular dynamics or contact dynamics simulations. The dimension of the set of possible forces and the number of sliding contacts agrees with the theoretical expectations. The indeterminacy of each component of the contact forces are found, as well as the an estimate for the diameter of the set of possible contact forces. We also show that contacts with high indeterminacy are located on force chains. The question of whether the simulation methods can represent a packing's memory of its formation is addressed.Comment: 12 pages, 13 figures, submitted to Phys Rev

Homogeneous cooling of rough, dissipative particles: Theory and simulations

We investigate freely cooling systems of rough spheres in two and three dimensions. Simulations using an event driven algorithm are compared with results of an approximate kinetic theory, based on the assumption of a generalized homogeneous cooling state. For short times $t$, translational and rotational energy are found to change linearly with $t$. For large times both energies decay like $t^{-2}$ with a ratio independent of time, but not corresponding to equipartition. Good agreement is found between theory and simulations, as long as no clustering instability is observed. System parameters, i.e. density, particle size, and particle mass can be absorbed in a rescaled time, so that the decay of translational and rotational energy is solely determined by normal restitution and surface roughness.Comment: 10 pages, 10 eps-figure

The powerful outburst in Hercules A

The radio source Hercules A resides at the center of a cooling flow cluster of galaxies at redshift z = 0.154. A Chandra X-ray image reveals a shock front in the intracluster medium (ICM) surrounding the radio source, about 160 kpc from the active galactic nucleus (AGN) that hosts it. The shock has a Mach number of 1.65, making it the strongest of the cluster-scale shocks driven by an AGN outburst found so far. The age of the outburst ~5.9e7 y, its energy about 3e61 erg and its mean power ~1.6e46 erg/s. As for the other large AGN outbursts in cooling flow clusters, this outburst overwhelms radiative losses from the ICM of the Hercules A cluster by a factor of ~100. It adds to the case that AGN outbursts are a significant source of preheating for the ICM. Unless the mechanical efficiency of the AGN in Hercules A exceeds 10%, the central black hole must have grown by more than 1.7e8 Msun to power this one outburst.Comment: 4 pages, 5 figures, accepted by ApJ

Using geophysical surveys to test tracer-based storage estimates in headwater catchments

Acknowledgements The authors are grateful to Stian Bradford, Chris Gabrielli, and Julie Timms for practical and logistical assistance. The provision of transport by Iain Malcolm and Ross Glover of Marine Scotland Science was greatly appreciated. We also thank the European Research Council ERC (project GA 335910 VEWA) for funding through the VeWa project and the Leverhulme Trust for funding through PLATO (RPG-2014-016).Peer reviewedPostprin

An Energetic AGN Outburst Powered by a Rapidly Spinning Supermassive Black Hole or an Accreting Ultramassive Black Hole

Powering the 10^62 erg nuclear outburst in the MS0735.6+7421 cluster central galaxy by accretion implies that its supermassive black hole (SMBH) grew by ~6x10^8 solar masses over the past 100 Myr. We place upper limits on the amount of cold gas and star formation near the nucleus of <10^9 solar masses and <2 solar masses per year, respectively. These limits imply that an implausibly large fraction of the preexisting cold gas in the bulge must have been consumed by its SMBH at the rate of ~3-5 solar masses per year while leaving no trace of star formation. Such a high accretion rate would be difficult to maintain by stellar accretion or the Bondi mechanism, unless the black hole mass approaches 10^11 solar masses. Its feeble nuclear luminosities in the UV, I, and X-ray bands compared to its enormous mechanical power are inconsistent with rapid accretion onto a ~5x10^9 solar mass black hole. We suggest instead that the AGN outburst is powered by a rapidly-spinning black hole. A maximally-spinning, 10^9 solar mass black hole contains enough rotational energy, ~10^62 erg, to quench a cooling flow over its lifetime and to contribute significantly to the excess entropy found in the hot atmospheres of groups and clusters. Two modes of AGN feedback may be quenching star formation in elliptical galaxies centered in cooling halos at late times. An accretion mode that operates in gas-rich systems, and a spin mode operating at modest accretion rates. The spin conjecture may be avoided in MS0735 by appealing to Bondi accretion onto a central black hole whose mass greatly exceeds 10^10 solar mass. The host galaxy's unusually large, 3.8 kpc stellar core radius (light deficit) may witness the presence of an ultramassive black hole.Comment: Accepted for publication in ApJ. Modifications: adopted slightly higher black hole mass using Lauer's M_SMBH vs L_bulge relation and adjusted related quantities; considered more seriously the consequences of a ultramassive black hole, motivated by new Kormendy & Bender paper published after our submission; other modifications per referee comments by Ruszkowsk

Jet Interactions with the Hot Halos of Clusters and Galaxies

X-ray observations of cavities and shock fronts produced by jets streaming through hot halos have significantly advanced our understanding of the energetics and dynamics of extragalactic radio sources. Radio sources at the centers of clusters have dynamical ages between ten and several hundred million years. They liberate between 1E58-1E62 erg per outburst, which is enough energy to regulate cooling of hot halos from galaxies to the richest clusters. Jet power scales approximately with the radio synchrotron luminosity to the one half power. However, the synchrotron efficiency varies widely from nearly unity to one part in 10,000, such that relatively feeble radio source can have quasar-like mechanical power. The synchrotron ages of cluster radio sources are decoupled from their dynamical ages, which tend to be factors of several to orders of magnitude older. Magnetic fields and particles in the lobes tend to be out of equipartition. The lobes may be maintained by heavy particles (e.g., protons), low energy electrons, a hot, diffuse thermal gas, or possibly magnetic (Poynting) stresses. Sensitive X-ray images of shock fronts and cavities can be used to study the dynamics of extragalactic radio sources.Comment: 10 pages, 3 figures, invited review, "Extragalactic Jets: Theory and Observation from Radio to Gamma Ray, held in Girdwood, Alaska, U.S.A. 21-24 May, 2007, minor text changes; one added referenc

Energy flows in vibrated granular media

We study vibrated granular media, investigating each of the three components of the energy flow: particle-particle dissipation, energy input at the vibrating wall, and particle-wall dissipation. Energy dissipated by interparticle collisions is well estimated by existing theories when the granular material is dilute, and these theories are extended to include rotational kinetic energy. When the granular material is dense, the observed particle-particle dissipation rate decreases to as little as 2/5 of the theoretical prediction. We observe that the rate of energy input is the weight of the granular material times an average vibration velocity times a function of the ratio of particle to vibration velocity. `Particle-wall' dissipation has been neglected in all theories up to now, but can play an important role when the granular material is dilute. The ratio between gravitational potential energy and kinetic energy can vary by as much as a factor of 3. Previous simulations and experiments have shown that E ~ V^delta, with delta=2 for dilute granular material, and delta ~ 1.5 for dense granular material. We relate this change in exponent to the departure of particle-particle dissipation from its theoretical value.Comment: 19 pages revtex, 10 embedded eps figures, accepted by PR

Controls on the temporal and spatial variability of soil moisture in a mountainous landscape: the signature of snow and complex terrain

The controls on the spatial distribution of soil moisture include static and dynamic variables. The superposition of static and dynamic controls can lead to different soil moisture patterns for a given catchment during wetting, draining, and drying periods. These relationships can be further complicated in snow-dominated mountain regions where soil water input by precipitation is largely dictated by the spatial variability of snow accumulation and melt. In this study, we assess controls on spatial and temporal soil moisture variability in a small (0.02 km&lt;sup&gt;2&lt;/sup&gt;), snow-dominated, semi-arid catchment by evaluating spatial correlations between soil moisture and site characteristics through different hydrologic seasons. We assess the relative importance of snow with respect to other catchment properties on the spatial variability of soil moisture and track the temporal persistence of those controls. Spatial distribution of snow, distance from divide, soil texture, and soil depth exerted significant control on the spatial variability of moisture content throughout most of the hydrologic year. These relationships were strongest during the wettest period and degraded during the dry period. As the catchment cycled through wet and dry periods, the relative spatial variability of soil moisture tended to remain unchanged. We suggest that the static properties in complex terrain (slope, aspect, soils) impose first order controls on the spatial variability of snow and resulting soil moisture patterns, and that the interaction of dynamic (timing of water input) and static influences propagate that relative constant spatial variability through most of the hydrologic year. The results demonstrate that snow exerts significant influence on how water is retained within mid-elevation semi-arid catchments and suggest that reductions in annual snowpacks associated with changing climate regimes may strongly influence spatial and temporal soil moisture patterns and catchment physical and biological processes

Interference of stochastic resonances: Splitting of Kramers' rate

We consider the escape of particles located in the middle well of a symmetric triple well potential driven sinusoidally by two forces such that the potential wells roll as in stochastic resonance and the height of the potential barrier oscillates symmetrically about a mean as in resonant activation. It has been shown that depending on their phase difference the application of these two synchronized signals may lead to a splitting of time averaged Kramers' escape rate and a preferential product distribution in a parallel chemical reaction in the steady state

The Open Cluster NGC 7789: I. Radial Velocities for Giant Stars

A total of 597 radial-velocity observations for 112 stars in the ~1.6 Gyr old open cluster NGC 7789 have been obtained since 1979 with the radial velocity spectrometer at the Dominion Astrophysical Observatory. The mean cluster radial velocity is -54.9 +/- 0.12 km/s and the dispersion is 0.86 km/s, from 50 constant-velocity stars selected as members from this radial-velocity study and the proper motion study of McNamara and Solomon (1981). Twenty-five stars (32%) among 78 members are possible radial-velocity variable stars, but no orbits are determined because of the sparse sampling. Seventeen stars are radial-velocity non-members, while membership estimates of six stars are uncertain. There is a hint that the observed velocity dispersion falls off at large radius. This may due to the inclusion of long-period binaries preferentially in the central area of the cluster. The known radial-velocity variables also seem to be more concentrated toward the center than members with constant velocity. Although this is significant at only the 85% level, when combined with similar result of Raboud and Mermilliod (1994) for three other clusters, the data strongly support the conclusion that mass segregation is being detected.Comment: 16 pages (including 3 figures) and 3 table