We develop a new numerical algorithm to model collisional cascades in debris disks. Because of the large dynamical range in particle masses, we solve the integro-differential equations describing erosive and catastrophic collisions in a particle-in-a-box approach, while treating the orbital dynamics of the particles in an approximate fashion. We employ a new scheme for describing erosive (cratering) collisions that yields a continuous set of outcomes as a function of colliding masses. We demonstrate the stability and convergence characteristics of our algorithm and compare it with other treatments. We show that incorporating the effects of erosive collisions results in a decay of the particle distribution that is significantly faster than with purely catastrophic collisions.

The observed flat rotation curves of galaxies require either the presence of dark matter in Newtonian gravitational potentials or a significant modification to the theory of gravity at galactic scales. Detecting relativistic Doppler shifts and gravitational effects in the rotation curves offers a tool for distinguishing between predictions of gravity theories that modify the inertia of particles and those that modify the field equations. These higher-order effects also allow us in principle, to test whether dark matter particles obey the equivalence principle. We calculate here the magnitudes of the relativistic Doppler and gravitational shifts expected in realistic models of galaxies in a general metric theory of gravity. We identify a number of observable quantities that measure independently the special- and general-relativistic effects in each galaxy and suggest that both effects might be detected in a statistical sense by combining appropriately the rotation curves of a large number of galaxies.

We study the structure of neutron stars in f(R) gravity theories with perturbative constraints. We derive the modified Tolman-Oppenheimer-Volkov equations and solve them for a polytropic equation of state. We investigate the resulting modifications to the masses and radii of neutron stars and show that observations of surface phenomena alone cannot break the degeneracy between altering the theory of gravity versus choosing a different equation of state of neutron-star matter. On the other hand, observations of neutron-star cooling, which depends on the density of matter at the stellar interior, can place significant constraints on the parameters of the theory.

Major observational efforts in the coming decade are designed to probe the equation of state of dark energy. Measuring a deviation of the equation-of-state parameter w from -1 would indicate a dark energy that cannot be represented solely by a cosmological constant. While it is commonly assumed that any implied modification to the LambdaCDM model amounts to the addition of new dynamical fields, we propose here a framework for investigating whether or not such new fields are required when cosmological observations are combined with a set of minimal assumptions about the nature of gravitational physics. In our approach, we treat the additional degrees of freedom as perturbatively constrained and calculate a number of observable quantities, such as the Hubble expansion rate and the cosmic acceleration, for a homogeneous Universe. We show that current observations place our Universe within the perturbative validity of our framework and allow for the presence of non-dynamical gravitational degrees of freedom at cosmological scales.