X-ray and microwave cluster scaling relations are immensely valuable for cosmological analysis. However, their power is limited by astrophysical systematics that bias mass estimates and introduce additional scatter. Turbulence injected into the intracluster medium via mass assembly contributes substantially to cluster non-thermal pressure support, a significant source of such uncertainties. We use an analytical model to compute the assembly-driven non-thermal pressure profiles of haloes based on Monte Carlo-generated accretion histories. We introduce a fitting function for the average non-thermal pressure fraction profile, which exhibits minimal dependence on redshift at fixed peak height. Using the model, we predict deviations from self-similarity and the intrinsic scatter in the Sunyaev--Zel'dovich effect observable-mass scaling relation ($Y_\mathrm{SZ}-M$) due solely to inter-cluster variation in mass accretion histories. We study the dependence of $Y_\mathrm{SZ}-M$ on aperture radius, cosmology, redshift, and mass limit. The model predicts $5-9\%$ scatter in $Y_\mathrm{SZ}-M$ at $z=0$, increasing as the aperture used to compute $Y_\mathrm{SZ}$ increases from $R_\mathrm{500c}$ to $5R_\mathrm{500c}$. The predicted scatter lies slightly below that of studies based on non-radiative hydro-simulations, illustrating that assembly history variance is likely responsible for a substantial fraction of scatter in $Y_\mathrm{SZ}-M$. This should be regarded as a lower bound, which will likely increase with the use of an updated gas density model that incorporates a more realistic response to halo assembly. As redshift increases, $Y_\mathrm{SZ}-M$ deviates more from self-similarity and scatter increases. We show that the $Y_\mathrm{SZ}-M$ residuals correlate strongly with the recent halo mass accretion rate, potentially providing an opportunity to infer the latter.

We present a semi-analytic model of satellite galaxies, SatGen, which can generate large samples of satellite populations for a host halo of desired mass, redshift, and assembly history. The model combines dark-matter halo merger trees, empirical relations for the galaxy-halo connection, and analytic prescriptions for tidal effects, dynamical friction, and ram pressure stripping. SatGen emulates cosmological zoom-in hydro-simulations in certain aspects. Satellites can reside in cored or cuspy DM subhaloes, depending on the halo response to baryonic physics that can be formulated from hydro-simulations and physical modeling. The subhalo profile and the stellar mass and size of a satellite evolves depending on its tidal mass loss and initial structure. The host galaxy can include a baryonic disc and a stellar bulge, each described by a density profile that allows analytic orbit integration. SatGen complements simulations by propagating the effect of halo response found in simulated field galaxies to satellites (not properly resolved in simulations) and outperforms simulations by sampling the halo-to-halo variance of satellite statistics and overcoming artificial disruption due to insufficient resolution. As a first application, we use the model to study satellites of Milky Way sized hosts, making it emulate simulations of bursty star formation and of smooth star formation, respectively, and to experiment with a disc potential in the host halo. Our model reproduces the observed satellite statistics reasonably well. Different physical recipes make a difference in satellite abundance and spatial distribution at the 25% level, not large enough to be distinguished by current observations given the halo-to-halo variance. The MW disc depletes satellites by 20% and has a subtle effect of diversifying the internal structure of satellites, important for alleviating certain small-scale problems.

We introduce the use of autoregressive normalizing flows for rapid likelihood-free inference of binary black hole system parameters from gravitational-wave data with deep neural networks. A normalizing flow is an invertible mapping on a sample space that can be used to induce a transformation from a simple probability distribution to a more complex one: if the simple distribution can be rapidly sampled and its density evaluated, then so can the complex distribution. Our first application to gravitational waves uses an autoregressive flow, conditioned on detector strain data, to map a multivariate standard normal distribution into the posterior distribution over system parameters. We train the model on artificial strain data consisting of IMRPhenomPv2 waveforms drawn from a five-parameter $(m_1, m_2, \phi_0, t_c, d_L)$ prior and stationary Gaussian noise realizations with a fixed power spectral density. This gives performance comparable to current best deep-learning approaches to gravitational-wave parameter estimation. We then build a more powerful latent variable model by incorporating autoregressive flows within the variational autoencoder framework. This model has performance comparable to Markov chain Monte Carlo and, in particular, successfully models the multimodal $\phi_0$ posterior. Finally, we train the autoregressive latent variable model on an expanded parameter space, including also aligned spins $(\chi_{1z}, \chi_{2z})$ and binary inclination $\theta_{JN}$, and show that all parameters and degeneracies are well-recovered. In all cases, sampling is extremely fast, requiring less than two seconds to draw $10^4$ posterior samples.

We investigate dynamical self-friction, the process by which material that is stripped from a subhalo torques its remaining bound remnant, which causes it to lose orbital angular momentum. By running idealized simulations of a subhalo orbiting within an analytical host halo potential, we isolate the effect of self-friction from traditional dynamical friction due to the host halo. While at some points in a subhalo's orbit the torque of the stripped material can boost the orbital angular momentum of the remnant, the net effect over the long term is orbital decay regardless of the initial orbital parameters or subhalo mass. In order to quantify the strength of self-friction, we run a suite of simulations spanning typical host-to-subhalo mass ratios and orbital parameters. We find that the time-scale for self-friction, defined as the exponential decay time of the subhalo's orbital angular momentum, scales with mass ratio and orbital circularity similar to standard dynamical friction. The decay time due to self-friction is roughly an order of magnitude longer, suggesting that self-friction only contributes at the 10 percent level. However, along more radial orbits, self-friction can occasionally dominate over dynamical friction close to pericentric passage, where mass stripping is intense. This is also the epoch at which the self-friction torque undergoes large and rapid changes in both magnitude and direction, indicating that self-friction is an important process to consider when modeling pericentric passages of subhaloes and their associated satellite galaxies.

Accurately predicting the abundance and structural evolution of dark matter subhaloes is crucial for understanding galaxy formation, modeling galaxy clustering, and constraining the nature of dark matter. Due to the nonlinear nature of subhalo evolution, cosmological $N$-body simulations remain its primary method of investigation. However, it has recently been demonstrated that such simulations are still heavily impacted by artificial disruption, diminishing the information content on small scales and reducing the reliability of all simulation-calibrated semi-analytical models. In this paper, we utilize the recently released DASH library of high-resolution, idealized simulations of the tidal evolution of subhaloes, which are unhindered by numerical overmerging due to discreteness noise or force softening, to calibrate an improved, more-accurate model of the evolution of the density profiles of subhaloes that undergo tidal heating and stripping within their host halo. By testing previous findings that the structural evolution of a tidally truncated subhalo depends solely on the fraction of mass stripped, independent of the details of the stripping, we identify an additional dependence on the initial subhalo concentration. We provide significantly improved fitting functions for the subhalo density profiles and structural parameters ($V_\mathrm{max}$ and $r_\mathrm{max}$) that are unimpeded by numerical systematics and applicable to a wide range of parameter space. This model will be an integral component of a future semi-analytical treatment of substructure evolution, which can be used to predict key quantities, such as the evolved subhalo mass function and annihilation boost factors, and validate such calculations performed with cosmological simulations.

We present a machine learning approach for estimating galaxy cluster masses, trained using both Chandra and eROSITA mock X-ray observations of 2,041 clusters from the Magneticum simulations. We train a random forest regressor, an ensemble learning method based on decision tree regression, to predict cluster masses using an input feature set. The feature set uses core-excised X-ray luminosity and a variety of morphological parameters, including surface brightness concentration, smoothness, asymmetry, power ratios, and ellipticity. The regressor is cross-validated and calibrated on a training sample of 1,615 clusters (80% of sample), and then results are reported as applied to a test sample of 426 clusters (20% of sample). This procedure is performed for two different mock observation series in an effort to bracket the potential enhancement in mass predictions that can be made possible by including dynamical state information. The first series is computed from idealized Chandra-like mock cluster observations, with high spatial resolution, long exposure time (1 Ms), and the absence of background. The second series is computed from realistic-condition eROSITA mocks with lower spatial resolution, short exposures (2 ks), instrument effects, and background photons modeled. We report a 20% reduction in the mass estimation scatter when either series is used in our random forest model compared to a standard regression model that only employs core-excised luminosity. The morphological parameters that hold the highest feature importance are smoothness, asymmetry, and surface brightness concentration. Hence, these parameters, which encode the dynamical state of the cluster, can be used to make more accurate predictions of cluster masses in upcoming surveys, offering a crucial step forward for cosmological analyses.

We present an online application, the SuperNova Interactive Fitter (SNIF), which can be used to interactively fit one-zone, analytical models to optical transient light curves powered by (1) the radioactive decay of 56Ni, (2) the spin-down of newly born magnetars, or (3) the interaction between supernova ejecta and pre-existing circumstellar material. This application utilizes the Open Astronomy Catalogs API to connect the user to thousands of publicly available transient light curves from the literature.

In keeping with the Luxembourg government's initiative to support the future use of space resources, ASIME 2018 was held in Belval, Luxembourg on April 16-17, 2018. The goal of ASIME 2018: Asteroid Intersections with Mine Engineering, was to focus on asteroid composition for advancing the asteroid in-space resource utilisation domain. What do we know about asteroid composition from remote-sensing observations? What are the potential caveats in the interpretation of Earth-based spectral observations? What are the next steps to improve our knowledge on asteroid composition by means of ground-based and space-based observations and asteroid rendez-vous and sample return missions? How can asteroid mining companies use this knowledge? ASIME 2018 was a two-day workshop of almost 70 scientists and engineers in the context of the engineering needs of space missions with in-space asteroid utilisation. The 21 Questions from the asteroid mining companies were sorted into the four asteroid science themes: 1) Potential Targets, 2) Asteroid-Meteorite Links, 3) In-Situ Measurements and 4) Laboratory Measurements. The Answers to those Questions were provided by the scientists with their conference presentations and collected by A. Graps or edited directly into an open-access collaborative Google document or inserted by A. Graps using additional reference materials. During the ASIME 2018, first day and second day Wrap-Ups, the answers to the questions were discussed further. New readers to the asteroid mining topic may find the Conversation boxes and the Mission Design discussions especially interesting.

(abridged) We present the Significant Cosmic Holes in Universe (SCHU) method for identifying cosmic voids and loops of filaments in cosmological datasets and assigning their statistical significance using techniques from topological data analysis. Persistent homology is used to find different dimensional holes. For dark matter halo catalogs and galaxy surveys, the 0-, 1-, and 2-dimensional holes can be identified with clusters, loops of filaments, and voids. The procedure overlays halos/galaxies on a 3D grid, and a distance-to-measure (DTM) function is calculated at each point of the grid. A filtration is generated over the lower-level sets of the DTM across increasing threshold values. The filtered simplicial complex can be used to summarize the birth/death times of the different dimension homology group generators (i.e., the holes). Persistence diagrams are produced from the dimension and birth/death times of each homology group generator. Using the persistence diagrams and bootstrap sampling, we explain how $p$-values can be assigned to each homology group generator. The homology group generators on a persistence diagram are not, in general, uniquely located back in the original dataset volume so we propose a method for finding a representation of the homology group generators. This method provides a novel, statistically rigorous approach for locating informative generators in cosmological datasets, which may be useful for providing complementary cosmological constraints on the effects of, for example, the sum of the neutrino masses. The method is tested on a Voronoi foam simulation, and then applied to a subset of the SDSS galaxy survey and a cosmological simulation. Lastly, we calculate Betti functions for two of the MassiveNuS simulations and discuss implications for using the persistent homology of the density field to help break degeneracy in the cosmological parameters.

The Bubble Nebula (or NGC 7635) is a parsec-scale seemingly spherical wind-blown bubble around the relatively unevolved O star BD+60$^\circ$2522. The small dynamical age of the nebula and significant space velocity of the star suggest that the Bubble Nebula might be a bow shock. We have run 2D hydrodynamic simulations to model the interaction of the central star's wind with the interstellar medium (ISM). The models cover a range of possible ISM number densities of $n=50-200 {\rm cm}^{-3}$ and stellar velocities of $v_{\star}=20-40$ km s$^{-1}$. Synthetic H$\alpha$ and 24 $\mu$m emission maps predict the same apparent spherical bubble shape with quantitative properties similar to observations. The synthetic maps also predict a maximum brightness similar to that from the observations and agree that the maximum brightness is at the apex of the bow shock. The best-matching simulation had $v_{\star}\approx20$ km s$^{-1}$ into an ISM with $n\sim100 {\rm cm}^{-3}$, at an angle of 60$^\circ$ with respect to the line of sight. Synthetic maps of soft ($0.3-2$ keV) and hard ($2-10$ keV) X-ray emission show that the brightest region is in the wake behind the star and not at the bow shock itself. The unabsorbed soft X-rays have luminosity $\sim10^{32}-10^{33}$ erg s$^{-1}$. The hard X-rays are fainter, luminosity $\sim 10^{30} - 10^{31}$ erg s$^{-1}$, and may be too faint for current X-ray instruments to successfully observe. Our results imply that the O star creates a bow shock as it moves through the ISM and in turn creates an asymmetric bubble visible at optical and infrared wavelengths, and predicted to be visible in X-rays. The Bubble Nebula does not appear to be unique, it could be just a favourably oriented very dense bow shock. The dense ISM surrounding BD+60$^\circ$2522 and its strong wind suggest that it could be a good candidate for detecting non-thermal emission.

The abundance and demographics of dark matter substructure is important for many areas in astrophysics and cosmological $N$-body simulations have been the primary tool used to investigate them. However, it has recently become clear that the simulations are subject to numerical artefacts, which hampers a proper treatment of the tidal evolution of subhaloes. Unfortunately, no analytical models that accurately describe subhalo evolution exist either. We therefore present a library of idealized, high resolution $N$-body simulations of the tidal evolution of individual subhaloes that can be used to calibrate semi-analytical models and to complement cosmological simulations. The simulations focus on minor mergers, i.e., the mass of the subhalo is much smaller than that of the host halo, such that the impact of dynamical friction is negligible. This setup allows the adoption of a fixed analytical potential for modelling the host halo. The dynamical evolution of subhaloes is followed with $N$-body computations. In the library, four parameters, two of which characterize the subhalo orbit with respect to the host halo, and the two concentrations of the host- and subhalo, are varied over the ranges encountered in cosmological simulations. We show several representative examples from the library that illustrate the evolution of the subhalo mass and velocity dispersion profiles. Additionally, we make publicly available a pre-trained non-parametric model of the subhalo mass evolution based on random forest regression. This model is able to interpolate the simulation data at the 0.1\,dex level and provides efficient access to the data for further use in modelling.

Rotational spin-up due to outgassing of comet nuclei has been identified as a possible mechanism for considerable mass-loss and splitting. We report a search for spin changes for three large Jupiter-family comets (JFCs): 14P/Wolf, 143P/Kowal-Mrkos, and 162P/Siding Spring. None of the three comets has detectable period changes, and we set conservative upper limits of 4.2 (14P), 6.6 (143P) and 25 (162P) minutes per orbit. Comparing these results with all eight other JFCs with measured rotational changes, we deduce that none of the observed large JFCs experiences significant spin changes. This suggests that large comet nuclei are less likely to undergo rotationally-driven splitting, and therefore more likely to survive more perihelion passages than smaller nuclei. We find supporting evidence for this hypothesis in the cumulative size distributions of JFCs and dormant comets, as well as in recent numerical studies of cometary orbital dynamics. We added 143P to the sample of 13 other JFCs with known albedos and phase-function slopes. This sample shows a possible correlation of increasing phase-function slopes for larger geometric albedos. Partly based on findings from recent space missions to JFCs, we hypothesise that this correlation corresponds to an evolutionary trend for JFCs. We propose that newly activated JFCs have larger albedos and steeper phase functions, which gradually decrease due to sublimation-driven erosion. If confirmed, this could be used to analyse surface erosion from ground and to distinguish between dormant comets and asteroids.

Observations of gravitational waves from inspiralling neutron star binaries---such as GW170817---can be used to constrain the nuclear equation of state by placing bounds on stellar tidal deformability. For slowly rotating neutron stars, the response to a weak quadrupolar tidal field is characterized by four internal-structure-dependent constants called "Love numbers." The tidal Love numbers $k_2^\text{el}$ and $k_2^\text{mag}$ measure the tides raised by the gravitoelectric and gravitomagnetic components of the applied field, and the rotational-tidal Love numbers $\mathfrak{f}^\text{o}$ and $\mathfrak{k}^\text{o}$ measure those raised by couplings between the applied field and the neutron star spin. In this work we compute these four Love numbers for perfect fluid neutron stars with realistic equations of state. We discover (nearly) equation-of-state independent relations between the rotational-tidal Love numbers and the moment of inertia, thereby extending the scope of I-Love-Q universality. We find that similar relations hold among the tidal and rotational-tidal Love numbers. These relations extend the applications of I-Love universality in gravitational-wave astronomy. As our findings differ from those reported in the literature, we derive general formulas for the rotational-tidal Love numbers in post-Newtonian theory and confirm numerically that they agree with our general-relativistic computations in the weak-field limit.

The interplanetary space probe Cassini/Huygens reached Saturn in July 2004 after seven years of cruise phase. The Cosmic Dust Analyzer (CDA) measures the interplanetary, interstellar and planetary dust in our solar system since 1999 and provided unique discoveries. In 1999, CDA detected interstellar dust in the inner solar system followed by the detection of electrical charges of interplanetary dust grains during the cruise phase between Earth and Jupiter. The instrument determined the composition of interplanetary dust and the nanometre sized dust streams originating from Jupiter's moon Io. During the approach to Saturn in 2004, similar streams of submicron grains with speeds in the order of 100 km/s were detected from Saturn's inner and outer ring system and are released to the interplanetary magnetic field. Since 2004 CDA measured more than one million dust impacts characterizing the dust environment of Saturn. The instrument is one of three experiments which discovered the active ice geysers located at the south pole of Saturn's moon Enceladus in 2005. Later, a detailed compositional analysis of the water ice grains in Saturn's E ring system lead to the discovery of large reservoirs of liquid water (oceans) below the icy crust of Enceladus. Finally, the determination of the dust- magnetosphere interaction and the discovery of the extended E ring (at least twice as large as predicted) allowed the definition of a dynamical dust model of Saturn's E ring describing the observed properties. This paper summarizes the discoveries of a ten year story of success based on reliable measurements with the most advanced dust detector flown in space until today. This paper focuses on cruise results and findings achieved at Saturn with a focus on flux and density measurements.

CASTAway is a mission concept to explore our Solar System's main asteroid belt. Asteroids and comets provide a window into the formation and evolution of our Solar System and the composition of these objects can be inferred from space-based remote sensing using spectroscopic techniques. Variations in composition across the asteroid populations provide a tracer for the dynamical evolution of the Solar System. The mission combines a long-range (point source) telescopic survey of over 10,000 objects, targeted close encounters with 10 to 20 asteroids and serendipitous searches to constrain the distribution of smaller (e.g. 10 m) size objects into a single concept. With a carefully targeted trajectory that loops through the asteroid belt, CASTAway would provide a comprehensive survey of the main belt at multiple scales. The scientific payload comprises a 50 cm diameter telescope that includes an integrated low-resolution (R = 30 to 100) spectrometer and visible context imager, a thermal (e.g. 6 to 16 microns) imager for use during the flybys, and modified star tracker cameras to detect small (approx. 10 m) asteroids. The CASTAway spacecraft and payload have high levels of technology readiness and are designed to fit within the programmatic and cost caps for a European Space Agency medium class mission, whilst delivering a significant increase in knowledge of our Solar System.

The processes that led to the formation of the planetary bodies in the Solar System are still not fully understood. Using the results obtained with the comprehensive suite of instruments on-board ESA's Rosetta mission, we present evidence that comet 67P/Churyumov-Gerasimenko likely formed through the gentle gravitational collapse of a bound clump of mm-sized dust aggregates ("pebbles"), intermixed with microscopic ice particles. This formation scenario leads to a cometary make-up that is simultaneously compatible with the global porosity, homogeneity, tensile strength, thermal inertia, vertical temperature profiles, sizes and porosities of emitted dust, and the steep increase in water-vapour production rate with decreasing heliocentric distance, measured by the instruments on-board the Rosetta spacecraft and the Philae lander. Our findings suggest that the pebbles observed to be abundant in protoplanetary discs around young stars provide the building material for comets and other minor bodies.

We report new lightcurves and phase functions for nine Jupiter-family comets (JFCs). They were observed in the period 2004-2015 with various ground telescopes as part of the Survey of Ensemble Physical Properties of Cometary Nuclei (SEPPCoN) as well as during devoted observing campaigns. We add to this a review of the properties of 35 JFCs with previously published rotation properties. The photometric time-series were obtained in Bessel R, Harris R and SDSS r' filters and were absolutely calibrated using stars from the Pan-STARRS survey. This specially-developed method allowed us to combine data sets taken at different epochs and instruments with absolute-calibration uncertainty down to 0.02 mag. We used the resulting time series to improve the rotation periods for comets 14P/Wolf, 47P/Ashbrook-Jackson, 94P/Russell, and 110P/Hartley 3 and to determine the rotation rates of comets 93P/Lovas and 162P/Siding-Spring for the first time. In addition to this, we determined the phase functions for seven of the examined comets and derived geometric albedos for eight of them. We confirm the known cut-off in bulk densities at $\sim$0.6 g $\mathrm{cm^{-3}}$ if JFCs are strengthless. Using the model of Davidsson (2001) for prolate ellipsoids with typical density and elongations, we conclude that none of the known JFCs require tensile strength larger than 10-25 Pa to remain stable against rotational instabilities. We find evidence for an increasing linear phase function coefficient with increasing geometric albedo. The median linear phase function coefficient for JFCs is 0.046 mag/deg and the median geometric albedo is 4.2 per cent.

We present a summary of the campaign of remote observations that supported the European Space Agency's Rosetta mission. Telescopes across the globe (and in space) followed comet 67P/Churyumov-Gerasimenko from before Rosetta's arrival until nearly the end of mission in September 2016. These provided essential data for mission planning, large-scale context information for the coma and tails beyond the spacecraft, and a way to directly compare 67P with other comets. The observations revealed 67P to be a relatively `well behaved' comet, typical of Jupiter family comets and with activity patterns that repeat from orbit-to-orbit. Comparison between this large collection of telescopic observations and the in situ results from Rosetta will allow us to better understand comet coma chemistry and structure. This work is just beginning as the mission ends -- in this paper we present a summary of the ground-based observations and early results, and point to many questions that will be addressed in future studies.

The aim of the Asteroid Science Intersections with In-Space Mine Engineering (ASIME) 2016 conference on September 21-22, 2016 in Luxembourg City was to provide an environment for the detailed discussion of the specific properties of asteroids, with the engineering needs of space missions that utilize asteroids. The ASIME 2016 Conference produced a layered record of discussions from the asteroid scientists and the asteroid miners to understand each other's key concerns and to address key scientific questions from the asteroid mining companies: Planetary Resources, Deep Space Industries and TransAstra. These Questions were the focus of the two day conference, were addressed by scientists inside and outside of the ASIME Conference and are the focus of this White Paper. The Questions from the asteroid mining companies have been sorted into the three asteroid science themes: 1) survey, 2) surface and 3) subsurface and 4) Other. The answers to those Questions have been provided by the scientists with their conference presentations or edited directly into an early open-access collaborative Google document (August 2016-October 2016), or inserted by A. Graps using additional reference materials. During the ASIME 2016 last two-hours, the scientists turned the Questions from the Asteroid Miners around by presenting their own key concerns: Questions from the Asteroid Scientists . These answers in this White Paper will point to the Science Knowledge Gaps (SKGs) for advancing the asteroid in-space resource utilisation domain.

We describe the full nonlinear development of the superradiant instability for a charged massless scalar field, coupled to general relativity and electromagnetism, in the vicinity of a Reissner-Nordstrom-AdS black hole. The presence of the negative cosmological constant provides a natural context for considering perfectly reflecting boundary conditions and studying the dynamics as the scalar field interacts repeatedly with the black hole. At early times, small superradiant perturbations grow as expected from linearized studies. Backreaction then causes the black hole to lose charge and mass until the perturbation becomes nonsuperradiant, with the final state described by a stable hairy black hole. For large gauge coupling, the instability extracts a large amount of charge per unit mass, resulting in greater entropy increase. We discuss the implications of the observed behavior for the general problem of superradiance in black hole spacetimes.

We provide a simple discussion of our results on the backreaction effects of density inhomogeneities in cosmology, without mentioning one-parameter families or weak limits. Emphasis is placed on the manner in which "averaging" is done and the fact that one is solving Einstein's equation. The key assumptions and results that we rigorously derived within our original mathematical framework are thereby explained in a heuristic way.

We respond to the criticisms of a recent paper of Buchert et al. [arXiv:1505.07800]

We present THERMAP, a mid-infrared (8-16 {\mu}m) spectro-imager for space missions to small bodies in the inner solar system, developed in the framework of the MarcoPolo-R asteroid sample return mission. THERMAP is very well suited to characterize the surface thermal environment of a NEO and to map its surface composition. The instrument has two channels, one for imaging and one for spectroscopy: it is both a thermal camera with full 2D imaging capabilities and a slit spectrometer. THERMAP takes advantage of the recent technological developments of uncooled microbolometers detectors, sensitive in the mid-infrared spectral range. THERMAP can acquire thermal images (8-18 {\mu}m) of the surface and perform absolute temperature measurements with a precision better than 3.5 K above 200 K. THERMAP can acquire mid-infrared spectra (8-16 {\mu}m) of the surface with a spectral resolution {\Delta}{\lambda} of 0.3 {\mu}m. For surface temperatures above 350 K, spectra have a signal-to-noise ratio >60 in the spectral range 9-13 {\mu}m where most emission features occur.

Extremely well! In the $\Lambda$CDM model, the spacetime metric, $g_{ab}$, of our universe is approximated by an FLRW metric, $g_{ab}^{(0)}$, to about 1 part in $10^4$ or better on both large and small scales, except in the immediate vicinity of very strong field objects, such as black holes. However, derivatives of $g_{ab}$ are not close to derivatives of $g_{ab}^{(0)}$, so there can be significant differences in the behavior of geodesics and huge differences in curvature. Consequently, observable quantities in the actual universe may differ significantly from the corresponding observables in the FLRW model. Nevertheless, as we shall review here, we have proven general results showing that---within the framework of our approach to treating backreaction---the large matter inhomogeneities that occur on small scales cannot produce significant effects on large scales, so $g_{ab}^{(0)}$ satisfies Einstein's equation with the averaged stress-energy tensor of matter as its source. We discuss the flaws in some other approaches that have suggested that large backreaction effects may occur. As we also will review here, with a suitable "dictionary," Newtonian cosmologies provide excellent approximations to cosmological solutions to Einstein's equation (with dust and a cosmological constant) on all scales. Our results thereby provide strong justification for the mathematical consistency and validity of the $\Lambda$CDM model within the context of general relativistic cosmology.

The Yarkovsky (orbital drift) and YORP (spin state change) effects play important roles in the dynamical and physical evolution of asteroids. Thermophysical modelling of these observed effects, and of thermal-infrared observations, allows a detailed physical characterisation of an individual asteroid to be performed. We perform a detailed physical characterisation of near-Earth asteroid (1620) Geographos, a potential meteor stream source and former spacecraft target, using the same techniques as previously used in Rozitis et al. (2013) for (1862) Apollo. We use the advanced thermophysical model (ATPM) on published light-curve, radar, and thermal-infrared observations to constrain the thermophysical properties of Geographos. The derived properties are used to make detailed predictions of the Yarkovsky orbital drift and YORP rotational acceleration, which are then compared against published measurements to determine Geographos's bulk density. We find that Geographos has a thermal inertia of 340 +140/-100 J m-2 K-1 s-1/2, a roughness fraction of >50%, and a bulk density of 2100 +550/-450 kg m-3 when using the light-curve-derived shape model with the radar-derived maximum equatorial diameter of 5.04 +/- 0.07 km. It is also found that the radar observations had overestimated the z-axis in Geographos's shape model because of their near-equatorial view. This results in a poor fit to the thermal-infrared observations if its effective diameter is kept fixed in the model fitting. The thermal inertia derived for Geographos is slightly higher than the typical values for a near-Earth asteroid of its size, and its derived bulk density suggests a rubble-pile interior structure. Large uncertainties in shape model z-axes are likely to explain why radar and thermal-infrared observations sometimes give inconsistent diameter determinations for other asteroids.

We consider perfect fluid bodies (stars) in general relativity, characterized by particle number density and entropy per particle. A star is said to be in dynamic equilibrium if it is a stationary, axisymmetric solution to the Einstein-fluid equations with circular flow. We prove that for a star in dynamic equilibrium, the necessary and sufficient condition for thermodynamic equilibrium (extremum of total entropy S) is constancy of angular velocity (\Omega), redshifted temperature, and redshifted chemical potential. Friedman previously identified positivity of canonical energy, E, as a criterion for dynamic stability and argued that all rotating stars are dynamically unstable to sufficiently nonaxisymmetric perturbations (the CFS instability), so our main focus is on axisymmetric stability. We show that for a star in dynamic equilibrium, mode stability holds with respect to all axisymmetric perturbations if E is positive on a certain subspace, V, of axisymmetric Lagrangian perturbations that, in particular, have no Lagrangian change in angular momentum density. Conversely, if E fails to be positive on V, then there exist perturbations that can't become asymptotically stationary at late times. We further show that for a star in thermodynamic equilibrium, the canonical energy in the rotating frame, E_r, is related to second order changes in ADM mass, M, and angular momentum, J, by E_r = \delta^2 M - \Omega \delta^2 J. Thus, positivity of E_r for perturbations with \delta J = 0 is a necessary condition for thermodynamic stability (local maximum of S). For axisymmetric perturbations we have E = E_r, so a necessary condition for thermodynamic stability with respect to axisymmetric perturbations is positivity of E on all perturbations with \delta J = 0, not merely on the perturbations in V. Many of our results are in close parallel to results of Hollands and Wald for the theory of black holes.

We study the dynamics of a 2+1 dimensional relativistic viscous conformal fluid in Minkowski spacetime. Such fluid solutions arise as duals, under the "gravity/fluid correspondence", to 3+1 dimensional asymptotically anti-de Sitter (AAdS) black brane solutions to the Einstein equation. We examine stability properties of shear flows, which correspond to hydrodynamic quasinormal modes of the black brane. We find that, for sufficiently high Reynolds number, the solution undergoes an inverse turbulent cascade to long wavelength modes. We then map this fluid solution, via the gravity/fluid duality, into a bulk metric. This suggests a new and interesting feature of the behavior of perturbed AAdS black holes and black branes, which is not readily captured by a standard quasinormal mode analysis. Namely, for sufficiently large perturbed black objects (with long-lived quasinormal modes), nonlinear effects transfer energy from short to long wavelength modes via a turbulent cascade within the metric perturbation. As long wavelength modes have slower decay, this lengthens the overall lifetime of the perturbation. We also discuss various implications of this behavior, including expectations for higher dimensions, and the possibility of predicting turbulence in more general gravitational scenarios.

In previous work, we introduced a new framework to treat large scale backreaction effects due to small scale inhomogeneities in general relativity. We considered one-parameter families of spacetimes for which such backreaction effects can occur, and we proved that, provided the weak energy condition on matter is satisfied, the leading effect of small scale inhomogeneities on large scale dynamics is to produce a traceless effective stress-energy tensor that itself satisfies the weak energy condition. In this work, we illustrate the nature of our framework by providing two explicit examples of one-parameter families with backreaction. The first, based on previous work of Berger, is a family of polarized vacuum Gowdy spacetimes on a torus, which satisfies all of the assumptions of our framework. As the parameter approaches its limiting value, the metric uniformly approaches a smooth background metric, but spacetime derivatives of the deviation of the metric from the background metric do not converge uniformly to zero. The limiting metric has nontrivial backreaction from the small scale inhomogeneities, with an effective stress-energy that is traceless and satisfies the weak energy condition, in accord with our theorems. Our second one-parameter family consists of metrics which have a uniform Friedmann-Lemaitre-Robertson-Walker limit. This family satisfies all of our assumptions with the exception of the weak energy condition for matter. In this case, the limiting metric has an effective stress-energy tensor which is not traceless. We emphasize the importance of imposing energy conditions on matter in studies of backreaction.

Surfaces of planets and small bodies of our Solar System are often covered by a layer of granular material that can range from a fine regolith to a gravel-like structure of varying depths. Therefore, the dynamics of granular materials are involved in many events occurring during planetary and small-body evolution thus contributing to their geological properties. We demonstrate that the new adaptation of the parallel N-body hard-sphere code pkdgrav has the capability to model accurately the key features of the collective motion of bidisperse granular materials in a dense regime as a result of shaking. As a stringent test of the numerical code we investigate the complex collective ordering and motion of granular material by direct comparison with laboratory experiments. We demonstrate that, as experimentally observed, the scale of the collective motion increases with increasing small-particle additive concentration. We then extend our investigations to assess how self-gravity and external gravity affect collective motion. In our reduced-gravity simulations both the gravitational conditions and the frequency of the vibrations roughly match the conditions on asteroids subjected to seismic shaking, though real regolith is likely to be much more heterogeneous and less ordered than in our idealised simulations. We also show that collective motion can occur in a granular material under a wide range of inter-particle gravity conditions and in the absence of an external gravitational field. These investigations demonstrate the great interest of being able to simulate conditions that are to relevant planetary science yet unreachable by Earth-based laboratory experiments.

Despite their very low surface gravities, the surfaces of asteroids and comets are covered by granular materials - regolith - that can range from a fine dust to a gravel-like structure of varying depths. Understanding the dynamics of granular materials is, therefore, vital for the interpretation of the surface geology of these small bodies and is also critical for the design and/or operations of any device planned to interact with their surfaces. We present the first measurements of transient weakening of granular material after shear reversal in microgravity as well as a summary of experimental results recently published in other journals, which may have important implications for small-body surfaces. Our results suggest that the force contact network within a granular material may be weaker in microgravity, although the influence of any change in the contact network is felt by the granular material over much larger distances. This could mean that small body surfaces are even more unstable than previously imagined. However, our results also indicate that the consequences of, e.g., a meteorite impact or a spacecraft landing, may be very different depending on the impact angle and location, and depending on the prior history of the small body surface.

We investigate the role of gravity on convection in a dense granular shear flow. Using a microgravity modified Taylor-Couette shear cell under the conditions of parabolic flight microgravity, we demonstrate experimentally that secondary, convective-like flows in a sheared granular material are close to zero in microgravity and enhanced under high-gravity conditions, though the primary flow fields are unaffected by gravity. We suggest that gravity tunes the frictional particle-particle and particle-wall interactions, which have been proposed to drive the secondary flow. In addition, the degree of plastic deformation increases with increasing gravitational forces, supporting the notion that friction is the ultimate cause.

Despite their very low surface gravities, asteroids exhibit a number of different geological processes involving granular matter. Understanding the response of this granular material subject to external forces in microgravity conditions is vital to the design of a successful asteroid sub-surface sampling mechanism, and in the interpretation of the fascinating geology on an asteroid. We have designed and flown a Taylor-Couette shear cell to investigate granular flow due to rotational shear forces under the conditions of parabolic flight microgravity. The experiments occur under weak compression. First, we present the technical details of the experimental design with particular emphasis on how the equipment has been specifically designed for the parabolic flight environment. Then, we investigate how a steady state granular flow induced by rotational shear forces differs in varying gravitational environments. We find that the effect of constant shearing on the granular material, in a direction perpendicular to the effective acceleration, does not seem to be strongly influenced by gravity. This means that shear bands can form in the presence of a weak gravitational field just as on Earth.

Near-Earth asteroid (1862) Apollo has strong detections of both orbital semimajor axis drift and rotational acceleration. We produce a unified model that can accurately match both observed effects using a single set of thermophysical properties derived from ground-based observations, and we determine Apollo's long term evolution. We use light-curve shape inversion techniques and the ATPM on published light-curve, thermal-infrared, and radar observations to constrain Apollo's thermophysical properties. The derived properties are used to make detailed predictions of Apollo's Yarkovsky and YORP effects, which are then compared with published measurements of orbital drift and rotational acceleration. The ATPM explicitly incorporates 1D heat conduction, shadowing, multiple scattering of sunlight, global self-heating, and rough surface thermal-infrared beaming in the model predictions. We find that ATPM can accurately reproduce the light-curve, thermal-infrared, and radar observations of Apollo, and simultaneously match the observed orbital drift and rotational acceleration using: a shape model with axis ratios of 1.94:1.65:1.00, an effective diameter of 1.55 +/- 0.07 km, a geometric albedo of 0.20 +/- 0.02, a thermal inertia of 140 +140/-100 J m-2 K-1 s-1/2, a highly rough surface, and a bulk density of 2850 +480/-680 kg m-3. Using these properties we predict that Apollo's obliquity is increasing towards the 180 degree YORP asymptotic state at a rate of 1.5 +0.3/-0.5 degrees per 10^5 yr. The derived thermal inertia suggests that Apollo has loose regolith material resting on its surface, which is consistent with Apollo undergoing a recent resurfacing event based on its observed Q-type spectrum. The inferred bulk density is consistent with those determined for other S-type asteroids, and suggests that Apollo has a fractured interior.

We present an adaptation of the ATPM to simultaneously predict the Yarkovsky and YORP effects in the presence of global self-heating that occurs within the large concavities of irregularly shaped asteroids, which has been neglected or dismissed in all previous models. It is also combined with rough surface thermal-infrared beaming effects, which have been previously shown to enhance the Yarkovsky-orbital-drift and dampen on average the YORP-rotational-acceleration by orders of several tens of per cent. Tests on all published concave shape models of near-Earth asteroids, and also on one hundred Gaussian-random-spheres, show that the Yarkovsky effect is sensitive to shadowing and global self-heating effects at the few per cent level or less. For simplicity, Yarkovsky models can neglect these effects if the level of accuracy desired is of this order. Unlike the Yarkovsky effect, the YORP effect can be very sensitive to shadowing and global self-heating effects. Its sensitivity increases with decreasing relative strength of the YORP-rotational-acceleration, and doesn't appear to depend greatly on the degree of asteroid concavity. Global self-heating tends to produce a vertical offset in an asteroid's YORP-rotational-acceleration versus obliquity curve which is in opposite direction to that produced by shadowing effects. It also ensures that at least one critical obliquity angle exists at which zero YORP-rotational-acceleration occurs. Global self-heating must be included for accurate predictions of the YORP effect if an asteroid exhibits a large shadowing effect. If global self-heating effects are not included then it is found in ~75 per cent of cases that better predictions are produced when shadowing is also not included. Furthermore, global self-heating has implications for reducing the sensitivity of the YORP effect predictions to detailed variations in an asteroid's shape model.

In addition to collisions and gravitational forces, it is now becoming widely acknowledged that photon recoil forces and torques from the asymmetric reflection and thermal re-radiation of sunlight are primary mechanisms that govern the rotational evolution of an asteroid. The Yarkovsky-O'Keefe-Radzievskii-Paddack (YORP) effect causes changes in the rotation rate and pole direction of an irregularly shaped asteroid. We present a simple Monte Carlo method to estimate the range of YORP-rotational-accelerations acting on a near-Earth asteroid (NEA) without knowledge of its detailed shape, and to estimate its detectability using light-curve observations. The method requires knowledge of an asteroid's orbital properties and size, and assumes that the future observational circumstances of an asteroid have already been thought through. It is verified by application to the observational circumstances of the seven YORP-investigated asteroids, and is then applied to 540 NEAs with NEOWISE and/or other diameter measurements, and to all NEAs using MPCORB absolute magnitudes. The YORP-detectability is found to be a strong function of the combined asteroid orbital and diameter properties, and is independent of the rotation period for NEAs that don't have very fast or slow rotation rates. The median and 1-sigma spread of YORP-rotational-acceleration expected to be acting on a particular NEA can be estimated from its semimajor axis, eccentricity, and diameter by using (eq. 1) and/or by using (eq. 2) if the diameter is instead estimated from the absolute magnitude by assuming a geometric albedo of 0.1. The length of a light-curve observational campaign required to achieve a 50 per cent probability of detecting the YORP effect in a particular NEA can be estimated by using (eq. 3) and/or by using (eq. 4) for an absolute-magnitude-estimated diameter.

We present a new rough-surface thermophysical model (Advanced Thermophysical Model or ATPM) that describes the observed directional thermal emission from any atmosphereless planetary surface. It explicitly incorporates partial shadowing, scattering of sunlight, selfheating and thermal-infrared beaming (re-radiation of absorbed sunlight back towards the Sun as a result of surface roughness). The model is verified by accurately reproducing ground-based directional thermal emission measurements of the lunar surface using surface properties that are consistent with the findings of the Apollo missions and roughness characterised by an RMS slope of ~32 degrees. By considering the wide range of potential asteroid surface properties, the model implies a beaming effect that cannot be described by a simple parameter or function. It is highly dependent on the illumination and viewing angles as well as surface thermal properties and is predominantly caused by macroscopic rather than microscopic roughness. Roughness alters the effective Bond albedo and thermal inertia of the surface as well as moving the mean emission away from the surface normal. For accurate determination of surface properties from thermal-infrared observations of unresolved bodies or resolved surface elements, roughness must be explicitly modelled, preferably aided with thermal measurements at different emission angles and wavelengths.

It is now becoming widely accepted that photon recoil forces from the asymmetric reflection and thermal re-radiation of absorbed sunlight are, together with collisions and gravitational forces, primary mechanisms governing the dynamical and physical evolution of asteroids. The Yarkovsky effect causes orbital semi-major axis drift and the YORP effect causes changes in the rotation rate and pole orientation. We present an adaptation of the Advanced Thermophysical Model (ATPM) to simultaneously predict the Yarkovsky and YORP effects in the presence of thermal-infrared beaming caused by surface roughness, which has been neglected or dismissed in all previous models. Tests on Gaussian random sphere shaped asteroids, and on the real shapes of asteroids (1620) Geographos and (6489) Golevka, show that rough surface thermal-infrared beaming enhances the Yarkovsky orbital drift by typically tens of percent but it can be as much as a factor of two. The YORP rotational acceleration is on average dampened by up to a third typically but can be as much as one half. We find that the Yarkovsky orbital drift is only sensitive to the average degree, and not to the spatial distribution, of roughness across an asteroid surface. However, the YORP rotational acceleration is sensitive to the surface roughness spatial distribution, and can add significant uncertainties to the predictions for asteroids with relatively weak YORP effects. To accurately predict either effect the degree and spatial distribution of roughness across an asteroid surface must be known.

String theory contains sources like orientifold planes that support higher derivative interactions. These interactions make possible static flux compactifications which are forbidden in supergravity. They can also lead to violations of the strong energy condition (SEC) which is needed for an accelerating universe. We examine how large a violation is possible in the context of the heterotic string compactified to four dimensions. We find that de Sitter solutions are still not possible but that classically forbidden anti-de Sitter solutions are possible.

Cosmological N-body simulations are now being performed using Newtonian gravity on scales larger than the Hubble radius. It is well known that a uniformly expanding, homogeneous ball of dust in Newtonian gravity satisfies the same equations as arise in relativistic FLRW cosmology, and it also is known that a correspondence between Newtonian and relativistic dust cosmologies continues to hold in linearized perturbation theory in the marginally bound/spatially flat case. Nevertheless, it is far from obvious that Newtonian gravity can provide a good global description of an inhomogeneous cosmology when there is significant nonlinear dynamical behavior at small scales. We investigate this issue in the light of a perturbative framework that we have recently developed, which allows for such nonlinearity at small scales. We propose a relatively straightforward "dictionary"---which is exact at the linearized level---that maps Newtonian dust cosmologies into general relativistic dust cosmologies, and we use our "ordering scheme" to determine the degree to which the resulting metric and matter distribution solve Einstein's equation. We find that Einstein's equation fails to hold at "order 1" at small scales and at "order $\epsilon$" at large scales. We then find the additional corrections to the metric and matter distribution needed to satisfy Einstein's equation to these orders. While these corrections are of some interest in their own right, our main purpose in calculating them is that their smallness should provide a criterion for the validity of the original dictionary (as well as simplified versions of this dictionary). We expect that, in realistic Newtonian cosmologies, these additional corrections will be very small; if so, this should provide strong justification for the use of Newtonian simulations to describe relativistic cosmologies, even on scales larger than the Hubble radius.

Asteroid (175706) 1996 FG3 is a binary asteroid and the baseline target for the proposed MarcoPolo-R sample-return mission. We present thermal IR photometry obtained with the ESO VLT+VISIR together with optical photometry obtained with the ESO NTT+EFOSC2 . An absolute visual magnitude H_V = 17.833 \pm 0.024 and phase parameter G = -0.041 \pm 0.005 is derived. The Near-Earth Asteroid Thermal Model (NEATM) has been fitted to the measured fluxes to derive a geometric visual albedo p_v = 0.046 \pm 0.014, effective diameter at the observed aspect D_eff = 1.68 \pm 0.25 km, and beaming parameter {\eta} = 1.15 for phase angle {\alpha} = 11.7{\deg}. The Advanced Thermophysical Model (ATPM) has been fitted to the measured fluxes to derive a more accurate effective diameter D_eff = 1.71 \pm 0.07 km and albedo p_v = 0.044 \pm 0.004. Based on the ATPM results, assuming the same albedo for primary and secondary, we derive a primary mean spherical diameter D_p = (1.69 +0.18/-0.12) km, secondary diameter D_s = 0.51 \pm 0.03 km, and a secondary orbital semi-major axis a = (2.8 +1.7/-0.7) km. A low surface thermal inertia {\Gamma} = 120 \pm 50 J m^-2 s^-1/2 K^-1 was also derived, suggesting a dusty surface and raising questions as to the binary formation mechanism of this asteroid. These physical properties are used to predict a Yarkovsky drift in semi-major axis of (-60 +31/-45) m yr-1.

We develop a new, mathematically precise framework for treating the effects of nonlinear phenomena occurring on small scales in general relativity. Our approach is an adaptation of Burnett's formulation of the "shortwave approximation", which we generalize to analyze the effects of matter inhomogeneities as well as gravitational radiation. Our framework requires the metric to be close to a "background metric", but allows arbitrarily large stress-energy fluctuations on small scales. We prove that, within our framework, if the matter stress-energy tensor satisfies the weak energy condition (i.e., positivity of energy density in all frames), then the only effect that small scale inhomogeneities can have on the dynamics of the background metric is to provide an "effective stress-energy tensor" that is traceless and has positive energy density---corresponding to the presence of gravitational radiation. In particular, nonlinear effects produced by small scale inhomogeneities cannot mimic the effects of dark energy. We also develop "perturbation theory" off of the background metric. We derive an equation for the "long-wavelength part" of the leading order deviation of the metric from the background metric, which contains the usual terms occurring in linearized perturbation theory plus additional contributions from the small-scale inhomogeneities. Under various assumptions concerning the absence of gravitational radiation and the non-relativistic behavior of the matter, we argue that the "short wavelength" deviations of the metric from the background metric near a point $x$ should be accurately described by Newtonian gravity, taking into account only the matter lying within a "homogeneity lengthscale" of $x$. Finally, we argue that our framework should provide an accurate description of the actual universe.

We conduct a survey of numerical simulations to probe the structure and appearance of non-radiative black hole accretion flows like the supermassive black hole at the Galactic centre. We find a generic set of solutions, and make specific predictions for currently feasible rotation measure (RM) observations, which are accessible to current instruments including the EVLA, GMRT and ALMA. The slow time variability of the RM is a key quantitative signature of this accretion flow. The time variability of RM can be used to quantitatively measure the nature of the accretion flow, and to differentiate models. Sensitive measurements of RM can be achieved using RM synthesis or using pulsars. Our energy conserving ideal magneto-hydrodynamical simulations, which achieve high dynamical range by means of a deformed-mesh algorithm, stretch from several Bondi radii to about one thousandth of that radius, and continue for tens of Bondi times. Magnetized flows which lack outward convection possess density slopes around -1, almost independent of physical parameters, and are more consistent with observational constraints than are strongly convective flows We observe no tendency for the flows to become rotationally supported in their centres, or to develop steady outflow. We support these conclusions with formulae which encapsulate our findings in terms of physical and numerical parameters. We discuss the relation of these solutions to other approaches. The main potential uncertainties are the validity of ideal MHD and the absence of a fully relativistic inner boundary condition. The RM variability predictions are testable with current and future telescopes.

The immense volume of data generated by the suite of instruments on SDO requires new tools for efficient identifying and accessing data that is most relevant to research investigations. We have developed the Heliophysics Events Knowledgebase (HEK) to fill this need. The HEK system combines automated data mining using feature-detection methods and high-performance visualization systems for data markup. In addition, web services and clients are provided for searching the resulting metadata, reviewing results, and efficiently accessing the data. We review these components and present examples of their use with SDO data.

The Near-Earth Asteroid Thermal Model (NEATM, Harris, 1998) has proven to be a reliable simple thermal model for radiometric diameter determination. However NEATM assumes zero thermal emission on the night side of an asteroid. We investigate how this assumption affects the best-fit beaming parameter, overestimates the effective diameter and underestimates the albedo at large phase angles, by testing NEATM on thermal IR fluxes generated from simulated asteroid surfaces with different thermal inertia. We compare NEATM to radar diameters and find that NEATM overestimates the diameter when the beaming parameter is fitted to multi-wavelength observations and underestimates the diameter when the default beaming parameter is used. The Night Emission Simulated Thermal Model (NESTM) is introduced. NESTM models the night side temperature as an iso-latitudinal fraction (f) of the maximum day side temperature (Maximum temperature calculated for NEATM with beaming parameter = 1). A range of f is found for different thermal parameters, which depend on the thermal inertia. NESTM diameters are compared with NEATM and radar diameters, and it is shown that NESTM may reduce the systematic bias in overestimating diameters. It is suggested that a version of the NESTM which assumes the thermal inertia = 200 S.I. units is adopted as a default model when the solar phase angle is greater than 45 degrees.

In September 1996, a dust/debris detector: GORID was launched into the geostationary (GEO) region as a piggyback instrument on the Russian Express-2 telecommunications spacecraft. The instrument began its normal operation in April 1997 and ended its mission in July 2002. The goal of this work was to use GORID's particle data to identify and separate the space debris to interplanetary dust particles (IDPs) in GEO, to more finely determine the instrument's measurement characteristics and to derive impact fluxes. While the physical characteristics of the GORID impacts alone are insufficient for a reliable distinction between debris and interplanetary dust, the temporal behavior of the impacts are strong enough indicators to separate the populations based on clustering. Non-cluster events are predominantly interplanetary, while cluster events are debris. The GORID mean flux distributions (at mass thresholds which are impact speed dependent) for IDPs, corrected for dead time, are 1.35x10^{-4} m^{-2} s^{-1} using a mean detection rate: 0.54 d^{-1}, and for space debris are 6.1x10^{-4} m^{-2} s^{-1} using a mean detection rate: 2.5 d^{-1}. Beta-meteoroids were not detected. Clusters could be a closely-packed debris cloud or a particle breaking up due to electrostatic fragmentation after high charging.

Gamma-ray bursts (GRBs) are thought to result from the interaction of an extremely relativistic outflow interacting with a small amount of material surrounding the site of the explosion. Multi-wavelength observations covering the gamma-ray to radio wavebands allow investigations of this `fireball' model. On 23 January 1999 optical emission was detected while the gamma-ray burst was still underway. Here we report the results of gamma-ray, optical/infra-red, sub-mm, mm and radio observations of this burst and its afterglow, which indicate that the prompt and afterglow emissions from GRB 990123 are associated with three distinct regions in the fireball. The afterglow one day after the burst has a much lower peak frequency than those of previous bursts; this explains the short-lived nature of the radio emission, which is not expected to reappear. We suggest that such differences reflect variations in the magnetic-field strengths in the afterglow emitting regions.