We present a new method, called x-cut cosmic shear, which optimally removes sensitivity to poorly modeled scales from the two-point cosmic shear signal. We show that the x-cut cosmic shear covariance matrix can be computed from the correlation function covariance matrix in a few minutes, enabling a likelihood analysis at virtually no additional computational cost. Further we show how to generalize x-cut cosmic shear to galaxy-galaxy lensing. Performing an x-cut cosmic shear analysis of the Dark Energy Survey Year 1 (DESY1) shear data, we reduce the error on S8 by 20% relative to a correlation function analysis, while showing our constraints are robust to different baryonic feedback models. Largely driven by information at small angular scales, our result, S8 = 0.734 plus/minus 0.038, yields a 2.6 sigma tension with the Planck Legacy analysis of the cosmic microwave background. As well as alleviating baryonic modelling uncertainties, our method can be used to optimally constrain a large number of theories of modified gravity where computational limitations make it infeasible to model the power spectrum down to extremely small scales. The key parts of our code are made publicly available.

We predict the mass fraction of oxygen lost from galaxies in a cosmological simulation as a function of stellar mass and environment at the present day. The distribution with stellar mass is bimodal, separating star-forming and quenched galaxies. The metallicity of gas and stars is self-consistently calculated using a chemical evolution model that includes supernovae type II and Ia, hypernovae, and asymptotic giant branch stars. The mass of oxygen lost from each galaxy is calculated by comparing the existing oxygen in gas and stars in the galaxy to the oxygen that should have been produced by the present-day population of stars. More massive galaxies are able to retain a greater fraction of their metals ($\sim 100$ per cent) than low-mass galaxies ($\sim 40 - 70$ per cent). As in the star formation main sequence, star-forming galaxies follow a tight relationship also in terms of oxygen mass lost -- a metal flow main sequence, ZFMS -- whereas massive quenched galaxies tend to have lost a greater fraction of oxygen (up to 20 per cent), due to AGN-driven winds. The amount of oxygen lost by satellite galaxies depends on the details of their interaction history, and those in richer groups tend to have lost a greater fraction of their oxygen. Observational estimates of metal retention in galaxies will provide a strong constraint on models of galaxy evolution.

Stage IV weak lensing experiments will offer more than an order of magnitude leap in precision. We must therefore ensure that our analyses remain accurate in this new era. Accordingly, previously ignored systematic effects must be addressed. In this work, we evaluate the impact of the reduced shear approximation and magnification bias, on the information obtained from the angular power spectrum. To first-order, the statistics of reduced shear, a combination of shear and convergence, are taken to be equal to those of shear. However, this approximation can induce a bias in the cosmological parameters that can no longer be neglected. A separate bias arises from the statistics of shear being altered by the preferential selection of galaxies and the dilution of their surface densities, in high-magnification regions. The corrections for these systematic effects take similar forms, allowing them to be treated together. We calculated the impact of neglecting these effects on the cosmological parameters that would be determined from Euclid, using cosmic shear tomography. To do so, we employed the Fisher matrix formalism, and included the impact of the super-sample covariance. We also demonstrate how the reduced shear correction can be calculated using a lognormal field forward modelling approach. These effects cause significant biases in Omega_m, sigma_8, n_s, Omega_DE, w_0, and w_a of -0.53 sigma, 0.43 sigma, -0.34 sigma, 1.36 sigma, -0.68 sigma, and 1.21 sigma, respectively. We then show that these lensing biases interact with another systematic: the intrinsic alignment of galaxies. Accordingly, we develop the formalism for an intrinsic alignment-enhanced lensing bias correction. Applying this to Euclid, we find that the additional terms introduced by this correction are sub-dominant.

During the first three observing runs of the Advanced gravitational-wave detector network, the LIGO/Virgo collaboration detected several black hole binary (BHBH) mergers. As the population of detected BHBH mergers grows, it will become possible to constrain different channels for their formation. Here we consider the chemically homogeneous evolution (CHE) channel in close binaries, by performing population synthesis simulations that combine realistic binary models with detailed cosmological calculations of the chemical and star-formation history of the Universe. This allows us to constrain population properties, as well as cosmological and aLIGO detection rates of BHBH mergers formed through this pathway. We predict a BHBH merger rate at redshift zero of $5.8 \hspace{1mm} \textrm{Gpc}^{-3} \textrm{yr}^{-1}$ through the CHE channel, to be compared with aLIGO's measured rate of ${53.2}_{-28.2}^{+55.8} \hspace{1mm} \text{Gpc}^{-3}\text{yr}^{-1}$, and find that eventual merger systems have BH masses in the range $17 - 43 \hspace{1mm} \textrm{M}_{\odot}$ below the pair-instability supernova (PISN) gap, and $>124 \hspace{1mm} \textrm{M}_{\odot}$ above the PISN gap. We further investigate the effects of momentum kicks during black hole formation, calculate cosmological and magnitude limited PISN rates and investigate the effects of high-redshift deviations in the star formation rate. We find that momentum kicks tend to increase delay times of BHBH systems, and our magnitude limited PISN rate estimates indicate that current deep surveys should be able to detect such events. Lastly, we find that our cosmological merger rate estimates change by at most $\sim 8\%$ for mild deviations of the star formation rate in the early Universe, and by up to $\sim 40\%$ for extreme deviations.

The Yarkovsky effect is a thermal process acting upon the orbits of small celestial bodies, which can cause these orbits to slowly expand or contract with time. The effect is subtle (da/dt ~ 10^-4 au/My for a 1 km diameter object) and is thus generally difficult to measure. We analyzed both optical and radar astrometry for 600 near-Earth asteroids (NEAs) for the purpose of detecting and quantifying the Yarkovsky effect. We present 247 NEAs with measured drift rates, which is the largest published set of Yarkovsky detections. This large sample size provides an opportunity to examine the Yarkovsky effect in a statistical manner. In particular, we describe two independent population-based tests that verify the measurement of Yarkovsky orbital drift. First, we provide observational confirmation for the Yarkovsky effect's theoretical size dependence of 1/D, where D is diameter. Second, we find that the observed ratio of negative to positive drift rates in our sample is 2.34, which, accounting for bias and sampling uncertainty, implies an actual ratio of $2.7^{+0.3}_{-0.7}$. This ratio has a vanishingly small probability of occurring due to chance or statistical noise. The observed ratio of retrograde to prograde rotators is two times lower than the ratio expected from numerical predictions from NEA population studies and traditional assumptions about the sense of rotation of NEAs originating from various main belt escape routes. We also examine the efficiency with which solar energy is converted into orbital energy and find a median efficiency in our sample of 12%. We interpret this efficiency in terms of NEA spin and thermal properties.

The Hubble Space Telescope (HST) has obtained multi-epoch observations providing the opportunity for a comprehensive variability search aiming to uncover new variables. We have therefore undertaken the task of creating a catalog of variable sources based on version 3 of the Hubble Source Catalog (HSC), which relies on publicly available images obtained with the WFPC2, ACS, and WFC3 instruments onboard the HST. We adopted magnitude-dependent thresholding in median absolute deviation (a robust measure of light curve scatter) combined with sophisticated preprocessing techniques and visual quality control to identify and validate variable sources observed by Hubble with the same instrument and filter combination five or more times. The Hubble Catalog of Variables (HCV) includes 84,428 candidate variable sources (out of 3.7 million HSC sources that were searched for variability) with $V \leq 27$ mag; for 11,115 of them the variability is detected in more than one filter. The data points in the light curves of the variables in the HCV catalog range from five to 120 points (typically having less than ten points); the time baseline ranges from under a day to over 15 years; while $\sim$8% of all variables have amplitudes in excess of 1 mag. Visual inspection performed on a subset of the candidate variables suggests that at least 80% of the candidate variables that passed our automated quality control are true variable sources rather than spurious detections resulting from blending, residual cosmic rays, and calibration errors. The HCV is the first, homogeneous catalog of variable sources created from archival HST data and currently is the deepest catalog of variables available. The catalog includes variable stars in our Galaxy and nearby galaxies, as well as transients and variable active galactic nuclei. (abbreviated)

Density-estimation likelihood-free inference (DELFI) has recently been proposed as an efficient method for simulation-based cosmological parameter inference. Compared to the standard likelihood-based Markov Chain Monte Carlo (MCMC) approach, DELFI has several advantages: it is highly parallelizable, there is no need to assume a possibly incorrect functional form for the likelihood and complicated effects (e.g the mask and detector systematics) are easier to handle with forward models. In light of this, we present two DELFI pipelines to perform weak lensing parameter inference with lognormal realizations of the tomographic shear field -- using the C_l summary statistic. The first pipeline accounts for the non-Gaussianities of the shear field, intrinsic alignments and photometric-redshift error. We validate that it is accurate enough for Stage III experiments and estimate that O(1000) simulations are needed to perform inference on Stage IV data. By comparing the second DELFI pipeline, which makes no assumption about the functional form of the likelihood, with the standard MCMC approach, which assumes a Gaussian likelihood, we test the impact of the Gaussian likelihood approximation in the MCMC analysis. We find it has a negligible impact on Stage IV parameter constraints. Our pipeline is a step towards seamlessly propagating all data-processing, instrumental, theoretical and astrophysical systematics through to the final parameter constraints.

The white paper discusses Arecibo Observatory's plan for facility improvements and activities over the next decade. The facility improvements include: (a) improving the telescope surface, pointing and focusing to achieve superb performance up to ~12.5 GHz; (b) equip the telescope with ultrawide-band feeds; (c) upgrade the instrumentation with a 4 GHz bandwidth high dynamic range digital link and a universal backend and (d) augment the VLBI facility by integrating the 12m telescope for phase referencing. These upgrades to the Arecibo telescope are critical to keep the national facility in the forefront of research in radio astronomy while maintaining its dominance in radar studies of near-Earth asteroids, planets and satellites. In the next decade, the Arecibo telescope will play a synergistic role with the upcoming facilities such as ngVLA, SKA and the now commissioned FAST telescope. Further, the observatory will be actively engaged in mentoring and training programs for students from a diverse background.

We present a method to measure the growth of structure and the background geometry of the Universe -- with no a priori assumption about the underlying cosmological model. Using Canada-France-Hawaii Lensing Survey (CFHTLenS) shear data we simultaneously reconstruct the lensing amplitude, the linear intrinsic alignment amplitude, the redshift evolving matter power spectrum, P(k,z), and the co-moving distance, r(z). We find that lensing predominately constrains a single global power spectrum amplitude and several co-moving distance bins. Our approach can localise precise scales and redshifts where Lambda-Cold Dark Matter (LCDM) fails -- if any. We find that below z = 0.4, the measured co-moving distance r (z) is higher than that expected from the Planck LCDM cosmology by ~1.5 sigma, while at higher redshifts, our reconstruction is fully consistent. To validate our reconstruction, we compare LCDM parameter constraints from the standard cosmic shear likelihood analysis to those found by fitting to the non-parametric information and we find good agreement.

Explaining the internal distribution and motions of stars and gas in galaxies is a key aspect in understanding their evolution. In previous work we identified five well resolved galaxies with atypical kinematics from a cosmological simulation; two had kinematically distinct cores (KDCs), and three had counter-rotating gas and stars (CRGD). In this paper, we show that i) the KDC galaxies have flattening of stellar [O/Fe] at large galacto-centric radii due to the minor mergers that gave rise to the KDCs, and ii) the CRGD galaxies have an abrupt transition in the gas metallicity maps, from high metallicity in the centre to very low metallicity further out. These galaxies are embedded in dark matter filaments where there is a ready supply of near-pristine gas to cause this effect. The non-linear increase in gas metallicity is also seen in the radial profiles, but when the metallicity gradients are measured, the difference is buried in the scatter of the relation. We also find that all five galaxies are fairly compact, with small effective radii given their stellar masses. This is because they have not experienced major mergers that kinematically heat the stars, and would have destroyed their unusual kinematics. In order to detect these signatures of minor mergers or accretion, the galaxy scaling relations or radial metallicity profiles are not enough, and it is necessary to obtain the 2D maps with integral field spectroscopy observations.

We report Arecibo S-band (2380 MHz; 12.6 cm) radar observations of near-Earth asteroid (3200) Phaethon during the December 2017 apparition when Phaethon passed within 0.07 au of Earth. Radar images with a resolution of 75 m per pixel reveal a roughly spheroidal shape more than 6 km in diameter at the equator with several discernible surface features hundreds of meters in extent. These include a possible crater more than 1 km across located below 30$^{\circ}$ latitude and a roughly 600-m radar-dark region near one of the poles. Overall, the radar images of Phaethon are reminiscent of those of (101955) Bennu, target of the OSIRIS-REx mission. As such, the shape of Phaethon is suspected to have an equatorial ridge similar to the top-shaped models of several other radar-observed near-Earth asteroids as well as the optical images of (162173) Ryugu returned by the Hayabusa2 spacecraft. Preliminary analysis of the radar data finds no satellites and gives no indication of a dusty coma at the time of these observations.

Pan-STARRS1 has carried out a set of distinct synoptic imaging sky surveys including the $3\pi$ Steradian Survey and the Medium Deep Survey in 5 bands ($grizy_{P1}$). The mean 5$\sigma$ point source limiting sensitivities in the stacked 3$\pi$ Steradian Survey in $grizy_{P1}$ are (23.3, 23.2, 23.1, 22.3, 21.4) respectively. The upper bound on the systematic uncertainty in the photometric calibration across the sky is 7-12 millimag depending on the bandpass. The systematic uncertainty of the astrometric calibration using the Gaia frame comes from a comparison of the results with Gaia: the standard deviation of the mean and median residuals ($ \Delta ra, \Delta dec $) are (2.3, 1.7) milliarcsec, and (3.1, 4.8) milliarcsec respectively. The Pan-STARRS system and the design of the PS1 surveys are described and an overview of the resulting image and catalog data products and their basic characteristics are described together with a summary of important results. The images, reduced data products, and derived data products from the Pan-STARRS1 surveys are available to the community from the Mikulski Archive for Space Telescopes (MAST) at STScI.

In this paper we address the problem of finding optimal cosmic shear tomographic bins. We generalise the definition of a cosmic shear tomographic bin to be a set of commonly labelled voxels in photometric colour space; rather than bins defined directly in redshift. We explore this approach by using a self-organising map to define the multi-dimensional colour space, and a we define a 'label space' of connected regions on the self-organising map using overlapping elliptical disks. This allows us to then find optimal labelling schemes by searching the label space. We use a metric that is the signal-to-noise ratio of a dark energy equation of state measurement, and in this case we find that for up to five tomographic bins the optimal colour-space labelling is an approximation of an equally-spaced binning in redshift; that is in all cases the best configuration. We also show that such a redefinition is more robust to photometric redshift outliers than a standard tomographic bin selection.

Most cosmological simulations of galaxy evolution include active galactic nucleus (AGN) feedback, typically seeding black holes with masses of $\geq 10^5\, h^{-1}\, \rm{M}_{\odot}$ when the dark matter halo exceeds a given threshold mass. Taylor & Kobayashi (2014) introduced a new model, which seeds black holes at $10^3\, h^{-1}\, \rm{M}_{\odot}$ based on gas properties alone, and motivated by the channel of black hole formation due to the collapse of the most massive first stars in the Universe. We compare the black hole mass when the dark matter halo mass is $10^{10}\, h^{-1}\, \rm{M}_{\odot}$ between the different seeding methods. We find that seeding based upon gas properties gives a distribution of black hole masses with $\langle \log M_{\rm{BH}} {/ \rm{M}_{\odot}} \rangle = (5.18 \pm 0.54)$ when dark matter halo mass is $10^{10}\, h^{-1}\, \rm{M}_{\odot}$, consistent with the {seeding criteria} used in other simulations. However, the evolution of individual galaxies can be strongly affected by the different seeding mechanisms. We also find that the mean value of the distribution of black hole masses at a given halo mass evolves over time, with higher masses at higher redshifts, indicative of downsizing. Our results can inform more physically motivated black hole and AGN feedback models in cosmological simulations and semi-analytic models.

Cosmic shear - the weak gravitational lensing effect generated by fluctuations of the gravitational tidal fields of the large-scale structure - is one of the most promising tools for current and future cosmological analyses. The spherical-Bessel decomposition of the cosmic shear field ("3D cosmic shear") is one way to maximise the amount of redshift information in a lensing analysis and therefore provides a powerful tool to investigate in particular the growth of cosmic structure that is crucial for dark energy studies. However, the computation of simulated 3D cosmic shear covariance matrices presents numerical difficulties, due to the required integrations over highly oscillatory functions. We present and compare two numerical methods and relative implementations to perform these integrations. We then show how to generate 3D Gaussian random fields on the sky in spherical coordinates, starting from the 3D cosmic shear covariances. To validate our field-generation procedure, we calculate the Minkowski functionals associated with our random fields, compare them with the known expectation values for the Gaussian case and demonstrate parameter inference from Minkowski functionals from a cosmic shear survey. This is a first step towards producing fully 3D Minkowski functionals for a lognormal field in 3D to extract Gaussian and non-Gaussian information from the cosmic shear field, as well as towards the use of Minkowski functionals as a probe of cosmology beyond the commonly used two-point statistics.

If left unchecked modeling uncertainties at small scales, due to poorly understood baryonic physics and non-linear structure formation, will significantly bias Stage IV cosmic shear two-point statistic parameter constraints. While it is perhaps possible to run N-body or hydrodynamical simulations to determine the impact of these effects this approach is computationally expensive; especially to test a large number of theories of gravity. Instead we propose directly removing sensitivity to small-scale structure from the lensing spectrum, creating a statistic that is robust to these uncertainties. We do this by taking a redshift-dependent l-cut after applying the Bernardeau-Nishimichi-Taruya (BNT) nulling scheme. This reorganizes the information in the lensing spectrum to make the relationship between the angular scale, l, and the structure scale, k, much clearer compared to standard cosmic shear power spectra -- for which no direct relationship exists. We quantify the effectiveness of this method at removing sensitivity to small scales and compute the predicted Fisher error on the dark energy equation of state, w0, for different k-cuts in the matter power spectrum.

Upcoming photometric lensing surveys will considerably tighten constraints on the neutrino mass and the dark energy equation of state. Nevertheless it remains an open question of how to optimally extract the information and how well the matter power spectrum must be known to obtain unbiased cosmological parameter estimates. By performing a Principal Component Analysis (PCA), we quantify the sensitivity of 3D cosmic shear and tomography with different binning strategies to different regions of the lensing kernel and matter power spectrum, and hence the background geometry and growth of structure in the Universe. We find that a large number of equally spaced tomographic bins in redshift can extract nearly all the cosmological information without the additional computational expense of 3D cosmic shear. Meanwhile a large fraction of the information comes from small poorly understood scales in the matter power spectrum, that can lead to biases on measurements of cosmological parameters. In light of this, we define and compute a cosmology-independent measure of the bias due to imperfect knowledge of the power spectrum. For a Euclid-like survey, we find that the power spectrum must be known to an accuracy of less than 1% on scales with k = 7 h /Mpc This requirement is not absolute since the bias depends on the magnitude of modelling errors, where they occur in k-z space, and the correlation between them, all of which are specific to any particular model. We therefore compute the bias in several of the most likely modelling scenarios and introduce a general formalism and public code, RequiSim, to compute the expected bias from any non-linear model.

We introduce the Generalised Lensing and Shear Spectra GLaSS code which is available for download from https://github.com/astro-informatics/GLaSS It is a fast and flexible public code, written in Python, that computes generalized spherical cosmic shear spectra. The commonly used tomographic and spherical Bessel lensing spectra come as built-in run-mode options. GLaSS is integrated into the Cosmosis modular cosmological pipeline package. We outline several computational choices that accelerate the computation of cosmic shear power spectra. Using GLaSS, we test whether the assumption that using the lensing and projection kernels for a spatially-flat universe -- in a universe with a small amount of spatial curvature -- negligibly impacts the lensing spectrum. We refer to this assumption as The Spatially-Flat Universe Approximation, that has been implicitly assumed in all cosmic shear studies to date. We confirm that The Spatially-Flat Universe Approximation has a negligible impact on Stage IV cosmic shear experiments.

We study how the predicted [$\alpha$/Fe]-[Fe/H] relations in early-type galaxies vary as functions of their stellar masses, ages and stellar velocity dispersions, by making use of cosmological chemodynamical simulations with feedback from active galactic nuclei. Our model includes a detailed treatment for the chemical enrichment from dying stars, core-collapse supernovae (both Type II and hypernovae) and Type Ia supernovae. At redshift $z=0$, we create a catalogue of $526$ galaxies, among which we determine $80$ early-type galaxies. From the analysis of our simulations, we find [$\alpha$/Fe]-[Fe/H] relations similar to the Galactic bulge. We also find that, in the oldest galaxies, Type Ia supernovae start to contribute at higher [Fe/H] than in the youngest ones. On the average, early-type galaxies with larger stellar masses (and, equivalently, higher stellar velocity dispersions) have higher [$\alpha$/Fe] ratios, at fixed [Fe/H]. This is qualitatively consistent with the recent observations of Sybilska et al., but quantitatively there are mismatches, which might require stronger feedback, sub-classes of Type Ia Supernovae, or a variable initial mass function to address.

Integral field spectroscopy surveys provide spatially resolved gas and stellar kinematics of galaxies. They have unveiled a range of atypical kinematic phenomena, which require detailed modelling to understand. We present results from a cosmological simulation that includes stellar and AGN feedback. We find that the distribution of angles between the gas and stellar angular momenta of galaxies is not affected by projection effects. We examine five galaxies ($\approx 6$ per cent of well resolved galaxies) that display atypical kinematics; two of the galaxies have kinematically distinct cores (KDC), while the other three have counter-rotating gas and stars. All five form the majority of their stars in the field, subsequently falling into cosmological filaments where the relative orientation of the stellar angular momentum and the bulk gas flow leads to the formation of a counter-rotating gas disc. The accreted gas exchanges angular momentum with pre-existing co-rotating gas causing it to fall to the centre of the galaxy. This triggers low-level AGN feedback, which reduces star formation. Later, two of the galaxies experience a minor merger (stellar mass ratio $\sim1/10$) with a galaxy on a retrograde orbit compared to the spin of the stellar component of the primary. This produces the KDCs, and is a different mechanism than suggested by other works. The role of minor mergers in the kinematic evolution of galaxies may have been under-appreciated in the past, and large, high-resolution cosmological simulations will be necessary to gain a better understanding in this area.

We present integral field spectroscopy of galaxy cluster Abell 3827, using ALMA and VLT/MUSE. It reveals an unusual configuration of strong gravitational lensing in the cluster core, with at least seven lensed images of a single background spiral galaxy. Lens modelling based on HST imaging had suggested that the dark matter associated with one of the cluster's central galaxies may be offset. The new spectroscopic data enable better subtraction of foreground light, and better identification of multiple background images. The inferred distribution of dark matter is consistent with being centered on the galaxies, as expected by LCDM. Each galaxy's dark matter also appears to be symmetric. Whilst we do not find an offset between mass and light (suggestive of self-interacting dark matter) as previously reported, the numerical simulations that have been performed to calibrate Abell 3827 indicate that offsets and asymmetry are still worth looking for in collisions with particular geometries. Meanwhile, ALMA proves exceptionally useful for strong lens image identifications.

The internal distribution of heavy elements, in particular the radial metallicity gradient, offers insight into the merging history of galaxies. Using our cosmological, chemodynamical simulations that include both detailed chemical enrichment and feedback from active galactic nuclei (AGN), we find that stellar metallicity gradients in the most massive galaxies ($\sim3\times10^{10}$M$_\odot$) are made flatter by mergers and are unable to regenerate due to the quenching of star formation by AGN feedback. The fitting range is chosen on a galaxy-by-galaxy basis in order to mask satellite galaxies. The evolutionary paths of the gradients can be summarised as follows; i) creation of initial steep gradients by gas-rich assembly, ii) passive evolution by star formation and/or stellar accretion at outskirts, iii) sudden flattening by mergers. There is a significant scatter in gradients at a given mass, which originates from the last path, and therefore from galaxy type. Some variation remains at given galaxy mass and type because of the complexity of merging events, and hence we find only a weak environmental dependence. Our early-type galaxies (ETGs), defined from the star formation main sequence rather than their morphology, are in excellent agreement with the observed stellar metallicity gradients of ETGs in the SAURON and ATLAS3D surveys. We find small positive [O/Fe] gradients of stars in our simulated galaxies, although they are smaller with AGN feedback. Gas-phase metallicity and [O/Fe] gradients also show variation, the origin of which is not as clear as for stellar populations.

Star formation in galaxies relies on the availability of cold, dense gas, which, in turn, relies on factors internal and external to the galaxies. In order to provide a simple model for how star formation is regulated by various physical processes in galaxies, we analyse data at redshift $z=0$ from a hydrodynamical cosmological simulation that includes prescriptions for star formation and stellar evolution, active galactic nuclei (AGN), and their associated feedback processes. This model can determine the star formation rate (SFR) as a function of galaxy stellar mass, gas mass, black hole mass, and environment. We find that gas mass is the most important quantity controlling star formation in low-mass galaxies, and star-forming galaxies in dense environments have higher SFR than their counterparts in the field. In high-mass galaxies, we find that black holes more massive than $\sim10^{7.5}$ M$_\odot$ can be triggered to quench star formation in their host; this mass scale is emergent in our simulations. Furthermore, this black hole mass corresponds to a galaxy bulge mass $\sim2\times10^{10}$ M$_\odot$, consistent with the mass at which galaxies start to become dominated by early types ($\sim3\times10^{10}$ M$_\odot$, as previously shown in observations by Kauffmann et al.). Finally, we demonstrate that our model can reproduce well the SFR measured from observations of galaxies in the GAMA and ALFALFA surveys.

Simulations of self-interacting dark matter (SIDM) predict that dark matter should lag behind galaxies during a collision. If the interaction is mediated by a high-mass force carrier, the distribution of dark matter can also develop asymmetric dark matter tails. To search for this asymmetry, we compute the gravitational lensing properties of a mass distribution with a free {\em skewness} parameter. We apply this to the dark matter around the four central galaxies in cluster Abell~3827. In the galaxy whose dark matter peak has previously been found to be offset, we tentatively measure a skewness $s=0.23^{+0.05}_{-0.22}$ in the same direction as the peak offset. Our method may be useful in future gravitational lensing analyses of colliding galaxy clusters and merging galaxies.

We present thermal properties and an improved shape model for potentially hazardous asteroid (162421) 2000 ET70. In addition to the radar data from 2000 ET70's apparition in 2012, our model incorporates optical lightcurves and infrared spectra that were not included in the analysis of Naidu et al. (2013, Icarus 226, 323-335). We confirm the general "clenched fist" appearance of the Naidu et al. model, but compared to their model, our best-fit model is about 10% longer along its long principal axis, nearly identical along the intermediate axis, and about 25% shorter along the short axis. We find the asteroid's dimensions to be 2.9 km $\times$ 2.2 km $\times$ 1.5 km (with relative uncertainties of about 10%, 15%, and 25%, respectively). With the available data, 2000 ET70's period and pole position are degenerate with each other. The radar and lightcurve data together constrain the pole direction to fall along an arc that is about twenty-three degrees long and eight degrees wide. Infrared spectra from the NASA InfraRed Telescope Facility (IRTF) provide an additional constraint on the pole. Thermophysical modeling, using our SHERMAN software, shows that only a subset of the pole directions, about twelve degrees of that arc, are compatible with the infrared data. Using all of the available data, we find that 2000 ET70 has a sidereal rotation period of 8.944 hours ($\pm$ 0.009 h) and a north pole direction of ecliptic coordinates $(52^{\circ}, -60^{\circ}) \pm 6^{\circ}$. The infrared data, acquired over several dates, require that the thermal properties (albedo, thermal inertia, surface roughness) must change across the asteroid's surface. By incorporating the detailed shape model and spin state into our thermal modeling, the multiple ground-based observations at different viewing geometries have allowed us to constrain the levels of the variations in the surface properties.

Near-Earth asteroid (NEA) 1566 Icarus ($a=1.08$ au, $e=0.83$, $i=22.8^{\circ}$) made a close approach to Earth in June 2015 at 22 lunar distances (LD). Its detection during the 1968 approach (16 LD) was the first in the history of asteroid radar astronomy. A subsequent approach in 1996 (40 LD) did not yield radar images. We describe analyses of our 2015 radar observations of Icarus obtained at the Arecibo Observatory and the DSS-14 antenna at Goldstone. These data show that the asteroid is a moderately flattened spheroid with an equivalent diameter of 1.44 km with 18% uncertainties, resolving long-standing questions about the asteroid size. We also solve for Icarus' spin axis orientation ($\lambda=270^{\circ}\pm10^{\circ}, \beta=-81^{\circ}\pm10^{\circ}$), which is not consistent with the estimates based on the 1968 lightcurve observations. Icarus has a strongly specular scattering behavior, among the highest ever measured in asteroid radar observations, and a radar albedo of $\sim$2\%, among the lowest ever measured in asteroid radar observations. The low cross-section suggests a high-porosity surface, presumably related to Icarus' cratering, spin, and thermal histories. Finally, we present the first use of our orbit determination software for the generation of observational ephemerides, and we demonstrate its ability to determine subtle perturbations on NEA orbits by measuring Icarus' orbit-averaged drift in semi-major axis ($(-4.62\pm0.48) \times 10^{-4}$ au/My, or $\sim$60 m per revolution). Our Yarkovsky rate measurement resolves a discrepancy between two published rates that did not include the 2015 radar astrometry.

Small near-Earth asteroids (>20 meters) are interesting because they are progenitors for meteorites in our terrestrial collection. Crucial to our understanding of the effectiveness of our atmosphere in filtering low-strength impactors is the physical characteristics of these small near-Earth asteroids (NEAs). In the past, characterization of small NEAs has been a challenge because of the difficulty in detecting them prior to close Earth flyby. In this study we physically characterized the 2-meter diameter near-Earth asteroid 2015 TC25 using ground-based optical, near-infrared and radar assets during a close flyby of the Earth (distance 69,000 miles) in Oct. 2015. Our observations suggest that its surface composition is similar to aubrites, a rare class of high albedo differentiated meteorites. Aubrites make up only 0.14 % of all know meteorites in our terrestrial meteorite collection. 2015 TC25 is also a very fast rotator with a rotation period of 133 seconds. We compared spectral and dynamical properties of 2015 TC25 and found the best candidate source body in the inner main belt to be the 70-km diameter E-type asteroid (44) Nysa. We attribute difference in spectral slope between the two objects to the lack of regolith on the surface of 2015 TC25. Using the albedo of E-type asteroids (50-60%) we refine the diameter of 2015 TC25 to 2-meters making it one of the smallest NEA ever to be characterized.

We present the first observation of cluster-scale radial metallicity gradients from star-forming galaxies. We use the DEIMOS spectrograph on the Keck II telescope to observe two CLASH clusters at z~0.35: MACS1115+0129 and RXJ1532+3021. Based on our measured interstellar medium (ISM) properties of star-forming galaxies out to a radius of 2.5 Mpc from the cluster centre, we find that the galaxy metallicity decreases as a function of projected cluster-centric distance (-0.15+/-0.08 dex/Mpc) in MACS1115+01. On the mass-metallicity relation (MZR), star-forming galaxies in MACS1115+01 are offset to higher metallicity (~0.2 dex) than the local SDSS galaxies at a fixed mass range. In contrast, the MZR of RXJ1532+30 is consistent with the local comparison sample. RXJ1532+30 exhibits a bimodal radial metallicity distribution, with one branch showing a similar negative gradient as MACS1115+01 (-0.14+/-0.05 dex/Mpc) and the other branch showing a positive radial gradient. The positive gradient branch in RXJ1532+30 is likely caused by either interloper galaxies or an in-plane merger, indicating that cluster-scale abundance gradients probe cluster substructures and thus the dynamical state of a cluster. Most strikingly, we discover that neither the radial metallicity gradient nor the offset from the MZR is driven by the stellar mass. We compare our observations with Rhapsody-G cosmological hydrodynamical zoom-in simulations of relaxed galaxy clusters and find that the simulated galaxy cluster also exhibits a negative abundance gradient, albeit with a shallower slope (-0.04+/-0.03 dex/Mpc). Our observations suggest that the negative radial gradient originates from ram-pressure stripping and/or strangulation processes in the cluster environments.

We predict the evolution of galaxy scaling relationships from cosmological, hydrodynamical simulations, that reproduce the scaling relations of present-day galaxies. Although we do not assume co-evolution between galaxies and black holes a priori, we are able to reproduce the black hole mass--velocity dispersion relation. This relation does not evolve, and black holes actually grow along the relation from significantly less massive seeds than have previously been used. AGN feedback does not very much affect the chemical evolution of our galaxies. In our predictions, the stellar mass--metallicity relation does not change its shape, but the metallicity significantly increases from $z\sim2$ to $z\sim1$, while the gas-phase mass-metallicity relation does change shape, having a steeper slope at higher redshifts ($z\lesssim3$). Furthermore, AGN feedback is required to reproduce observations of the most massive galaxies at $z\lesssim1$, specifically their positions on the star formation main sequence and galaxy mass--size relation.

We evaluated the planetary radar capabilities at Arecibo, the Goldstone 70-m DSS-14 and 34-m DSS-13 antennas, the 70-m DSS-43 antenna at Canberra, the Green Bank Telescope, and the Parkes Radio Telescope in terms of their relative sensitivities and the number of known near-Earth asteroids (NEAs) detectable per year in monostatic and bistatic configurations. In the 2015 calendar year, monostatic observations with Arecibo and DSS-14 were capable of detecting 253 and 131 NEAs respectively, with signal-to-noise ratios (SNRs) greater than 30/track. Combined, the two observatories were capable of detecting 276 NEAs. Of these, Arecibo detected 77 and Goldstone detected 32, or 30% and 24% the numbers that were possible. The two observatories detected an additional 18 and 7 NEAs respectively, with SNRs of less than 30/track. This indicates that a substantial number of potential targets are not being observed. The bistatic configuration with DSS-14 transmitting and the Green Bank Telescope receiving was capable of detecting about 195 NEAs, or ~50% more than with monostatic observations at DSS-14. Most of the detectable asteroids were targets of opportunity that were discovered less than 15 days before the end of their observing windows. About 50% of the detectable asteroids have absolute magnitudes > 25, which corresponds diameters < ~30 m.

We present a joint optical/X-ray analysis of the massive galaxy cluster Abell 2744 (z=0.308). Our strong- and weak-lensing analysis within the central region of the cluster, i.e., at R<1Mpc from the brightest cluster galaxy, reveals eight substructures, including the main core. All of these dark-matter halos are detected with a significance of at least 5sigma and feature masses ranging from 0.5 to 1.4x10^{14}Msun within R<150kpc. Merten et al. (2011) and Medezinski et al. (2016) substructures are also detected by us. We measure a slightly higher mass for the main core component than reported previously and attribute the discrepancy to the inclusion of our tightly constrained strong-lensing mass model built on Hubble Frontier Fields data. X-ray data obtained by XMM-Newton reveal four remnant cores, one of them a new detection, and three shocks. Unlike Merten et al. (2011), we find all cores to have both dark and luminous counterparts. A comparison with clusters of similar mass in the MXXL simulations yields no objects with as many massive substructures as observed in Abell 2744, confirming that Abell 2744 is an extreme system. We stress that these properties still do not constitute a challenge to $\Lambda$CDM, as caveats apply to both the simulation and the observations: for instance, the projected mass measurements from gravitational lensing and the limited resolution of the sub-haloes finders. We discuss implications of Abell 2744 for the plausibility of different dark-matter candidates and, finally, measure a new upper limit on the self-interaction cross-section of dark matter of sigma_{DM}<1.28cm2/g(68\% CL), in good agreement with previous results from Harvey et al. (2015).

In the past decade, the number of known binary near-Earth asteroids has more than quadrupled and the number of known large main belt asteroids with satellites has doubled. Half a dozen triple asteroids have been discovered, and the previously unrecognized populations of asteroid pairs and small main belt binaries have been identified. The current observational evidence confirms that small (<20 km) binaries form by rotational fission and establishes that the YORP effect powers the spin-up process. A unifying paradigm based on rotational fission and post-fission dynamics can explain the formation of small binaries, triples, and pairs. Large (>20 km) binaries with small satellites are most likely created during large collisions.

We show the effects of AGN-driven outflows on the ejection of heavy elements using our cosmological simulations, where super-massive black holes originate from the first stars. In the most massive galaxy, we have identified two strong outflows unambiguously driven by AGN feedback. These outflows have a speed greater than $\sim 8000$ km\,s$^{-1}$ near the AGN, and travel out to a half Mpc with $\sim 3000$ km\,s$^{-1}$. These outflows remove the remaining gas ($\sim 3$ per cent of baryons) and significant amounts of metals ($\sim 2$ per cent of total produced metals) from the host galaxy, chemically enriching the circumgalactic medium (CGM) and the intergalactic medium (IGM). 17.6 per cent of metals from this galaxy, and 18.4 per cent of total produced metals in the simulation, end up in the CGM and IGM, respectively. The metallicities of the CGM and IGM are higher with AGN feedback, while the mass--metallicity relation of galaxies is not affected very much. We also find `selective' mass-loss where iron is more effectively ejected than oxygen because of the time-delay of Type Ia Supernovae. AGN-driven outflows play an essential role not only in quenching of star formation in massive galaxies to match with observed down-sizing phenomena, but also in a large-scale chemical enrichment in the Universe. Observational constraints of metallicities and elemental abundance ratios in outflows are important to test the modelling of AGN feedback in galaxy formation.

Potentially hazardous asteroid (185851) 2000 DP107 was the first binary near-Earth asteroid to be imaged. Radar observations in 2000 provided images at 75 m resolution that revealed the shape, orbit, and spin-up formation mechanism of the binary. The asteroid made a more favorable flyby of the Earth in 2008, yielding images at 30 m resolution. We used these data to obtain shape models for the two components and to improve the estimates of the mutual orbit, component masses, and spin periods. The primary has a sidereal spin period of 2.7745 +/- 0.0007 h and is roughly spheroidal with an equivalent diameter of 863 m +/- 5 %. It has a mass of 4.656 +/- 0.43 x 10^11 kg and a density of 1381 +/- 244 kg m^{-3}. It exhibits an equatorial ridge similar to the (66391) 1999 KW4 primary, however the equatorial ridge in this case is not as regular and has a ~300 m diameter concavity on one side. The secondary has a sidereal spin period of 1.77 +/- 0.02 days commensurate with the orbital period. The secondary is slightly elongated and has overall dimensions of 377 x 314 x 268 m (6 % uncertainties). Its mass is 0.178 +/- 0.021 x 10^{11} kg and its density is 1047 +/- 230 kg m^{-3}. The mutual orbit has a semi-major axis of 2.659 +/- 0.08 km, an eccentricity of 0.019 +/- 0.01, and a period of 1.7556 +/- 0.0015 days. The normalized total angular momentum of this system exceeds the amount required for the expected spin-up formation mechanism. An increase of angular momentum from non-gravitational forces after binary formation is a possible explanation. The two components have similar radar reflectivity, suggesting a similar composition consistent with formation by spin-up. The secondary appears to exhibit a larger circular polarization ratio than the primary, suggesting a rougher surface or subsurface at radar wavelength scales.

We investigate the environmental dependence of the mass-metallicity relation at z=2 with MOSFIRE/Keck as part of the ZFIRE survey. Here, we present the chemical abundance of a Virgo-like progenitor at z=2.095 that has an established red sequence. We identified 43 cluster ($<z>=2.095\pm0.004$) and 74 field galaxies ($<z>=2.195\pm0.083$) for which we can measure metallicities. For the first time, we show that there is no discernible difference between the mass-metallicity relation of field and cluster galaxies to within 0.02dex. Both our field and cluster galaxy mass-metallicity relations are consistent with recent field galaxy studies at z~2. We present hydrodynamical simulations for which we derive mass-metallicity relations for field and cluster galaxies. We find at most a 0.1dex offset towards more metal-rich simulated cluster galaxies. Our results from both simulations and observations are suggestive that environmental effects, if present, are small and are secondary to the ongoing inflow and outflow processes that are governed by galaxy halo mass.

We show that feedback from active galactic nuclei (AGN) plays an essential role in reproducing the down-sizing phenomena, namely: the colour-magnitude relation; specific star formation rates; and the $\alpha$ enhancement of early type galaxies. In our AGN model, black holes originate from Population III stars, in contrast to the merging scenario of previous works. In this paper, we show how the properties of present-day galaxies in cosmological chemo-hydrodynamical simulations change when we include our model for AGN feedback. Massive galaxies become redder, older, less massive, less compact, and show greater $\alpha$ enhancement than their counterparts without AGN. Since we reproduce the black hole mass and galaxy mass relation, smaller galaxies do not host a supermassive black hole and their star formation history is affected very little, but they can get external enrichment from nearby AGN depending on their environment. Nonetheless, the metallicity change is negligible, and the mass--metallicity relations, which are mainly generated by supernova feedback at the first star burst, are preserved.

We present a new model for the formation of black holes in cosmological simulations, motivated by the first star formation. Black holes form from high density peaks of primordial gas, and grow via both gas accretion and mergers. Massive black holes heat the surrounding material, suppressing star formation at the centres of galaxies, and driving galactic winds. We perform an investigation into the physical effects of the model parameters, and obtain a `best' set of these parameters by comparing the outcome of simulations to observations. With this best set, we successfully reproduce the cosmic star formation rate history, black hole mass -- velocity dispersion relation, and the size -- velocity dispersion relation of galaxies. The black hole seed mass is 10^3Msun, which is orders of magnitude smaller than has been used in previous cosmological simulations with active galactic nuclei, but suggests that the origin of the seed black holes is the death of Population III stars.

The target asteroid of the OSIRIS-REx asteroid sample return mission, (101955) Bennu (formerly 1999 RQ$_{36}$), is a half-kilometer near-Earth asteroid with an extraordinarily well constrained orbit. An extensive data set of optical astrometry from 1999--2013 and high-quality radar delay measurements to Bennu in 1999, 2005, and 2011 reveal the action of the Yarkovsky effect, with a mean semimajor axis drift rate $da/dt = (-19.0 \pm 0.1)\times 10^{-4}$ au/Myr or $284\pm 1.5\;\rm{m/yr}$. The accuracy of this result depends critically on the fidelity of the observational and dynamical model. As an example, neglecting the relativistic perturbations of the Earth during close approaches affects the orbit with $3\sigma$ significance in $da/dt$. The orbital deviations from purely gravitational dynamics allow us to deduce the acceleration of the Yarkovsky effect, while the known physical characterization of Bennu allows us to independently model the force due to thermal emissions. The combination of these two analyses yields a bulk density of $\rho = 1260\pm70\,\rm{kg/m^3}$, which indicates a macroporosity in the range $40\pm10$% for the bulk densities of likely analog meteorites, suggesting a rubble-pile internal structure. The associated mass estimate is $(7.8\pm0.9)\times 10^{10}\, \rm{kg}$ and $GM = 5.2\pm0.6\,\rm{m^3/s^2}$. Bennu's Earth close approaches are deterministic over the interval 1654--2135, beyond which the predictions are statistical in nature. In particular, the 2135 close approach is likely within the lunar distance and leads to strong scattering and therefore numerous potential impacts in subsequent years, from 2175--2196. The highest individual impact probability is $9.5\times 10^{-5}$ in 2196, and the cumulative impact probability is $3.7\times 10^{-4}$, leading to a cumulative Palermo Scale of -1.70.

We derive the locations of the fully synchronous end states of tidal evolution for binary asteroid systems having one spherical component and one oblate- or prolate-spheroid component. Departures from a spherical shape, at levels observed among binary asteroids, can result in the lack of a stable tidal end state for particular combinations of the system mass fraction and angular momentum, in which case the binary must collapse to contact. We illustrate our analytical results with near-Earth asteroids (8567) 1996 HW1, (66391) 1999 KW4, and 69230 Hermes.

We observed near-Earth asteroid (162421) 2000 ET70 using the Arecibo and Goldstone radar systems over a period of 12 days during its close approach to the Earth in February 2012. We obtained continuous wave spectra and range-Doppler images with range resolutions as fine as 15 m. Inversion of the radar images yields a detailed shape model with an effective spatial resolution of 100 m. The asteroid has overall dimensions of 2.6 km X 2.2 km X 2.1 km (5% uncertainties) and a surface rich with kilometer-scale ridges and concavities. This size, combined with absolute magnitude measurements, implies an extremely low albedo (~2%). It is a principal axis rotator and spins in a retrograde manner with a sidereal spin period of 8.96 +/- 0.01 hours. In terms of gravitational slopes evaluated at scales of 100 m, the surface seems mostly relaxed with over 99% of the surface having slopes less than 30 degrees, but there are some outcrops at the north pole that may have steeper slopes. Our precise measurements of the range and velocity of the asteroid, combined with optical astrometry, enables reliable trajectory predictions for this potentially hazardous asteroid in the interval 460-2813.

We present a relativistic model for the stationary axisymmetric accretion flow of a rotating cloud of non-interacting particles falling onto a Kerr black hole. Based on a ballistic approximation, streamlines are described analytically in terms of timelike geodesics, while a simple numerical scheme is introduced for calculating the density field. A novel approach is presented for describing all of the possible types of orbit by means of a single analytic expression. This model is a useful tool for highlighting purely relativistic signatures in the accretion flow dynamics coming from a strong gravitational field with frame-dragging. In particular, we explore the coupling due to this between the spin of the black hole and the angular momentum of the infalling matter. Moreover, we demonstrate how this analytic solution may be used for benchmarking general relativistic numerical hydrodynamics codes by comparing it against results of smoothed particle hydrodynamics simulations for a collapsar-like setup. These simulations are performed first for a ballistic flow (with zero pressure) and then for a hydrodynamical one where we measure the effects of pressure gradients on the infall, thus exploring the extent of applicability of the ballistic approximation.

A commonly cited drawback of SPH is the introduction of spurious shear viscosity by the artificial viscosity term in situations involving rotation. Existing approaches for quantifying its effect include approximate analytic formulae and disc-averaged be- haviour in specific ring-spreading simulations, based on the kinematic effects produced by the artificial viscosity. These methods have disadvantages, in that they typically are applicable to a very small range of physical scenarios, have a large number of simplifying assumptions, and often are tied to specific SPH formulations which do not include corrective (e.g., Balsara) or time-dependent viscosity terms. In this study we have developed a simple, generally applicable and practical technique for evaluating the local effect of artificial viscosity directly from the creation of specific entropy for each SPH particle. This local approach is simple and quick to implement, and it al- lows a detailed characterization of viscous effects as a function of position. Several advantages of this method are discussed, including its ease in evaluation, its greater accuracy and its broad applicability. In order to compare this new method with ex- isting ones, simple disc flow examples are used. Even in these basic cases, the very roughly approximate nature of the previous methods is shown. Our local method pro- vides a detailed description of the effects of the artificial viscosity throughout the disc, even for extended examples which implement Balsara corrections. As a further use of this approach, explicit dependencies of the effective viscosity in terms of SPH and flow parameters are estimated from the example cases. In an appendix, a method for the initial placement of SPH particles is discussed which is very effective in reducing numerical fluctuations.

The locations of the fully despun, double synchronous end states of tidal evolution are derived for spherical components. With the exception of nearly equal-mass binaries, binary asteroid systems are in the midst of lengthy tidal evolutions, far from their fully synchronous tidal end states. Calculations of material strength indicate that binaries in the main belt with 100-km-scale primary components are consistent with being made of monolithic or fractured rock as expected for binaries likely formed from sub-catastrophic impacts in the early solar system. To tidally evolve in their dynamical lifetime, near-Earth binaries with km-scale primaries or smaller must be much weaker mechanically than their main-belt counterparts even if formed in the main belt prior to injection into the near-Earth region. Small main-belt binaries with primary components less than 10 km in diameter, depending on their ages, could either be as strong as large main-belt binaries or as weak as near-Earth binaries because the inherent uncertainty in the age of a binary system can affect the calculation of material strength by orders of magnitude. Several other issues are considered, though these typically affect the calculation of material strength by no more than a factor of two. We also find indirect evidence within all three groups of binary asteroids that the inter-component separation may evolve via another mechanism(s) with the binary YORP effect being a likely candidate.

Three-body model fits to Arecibo and Goldstone radar data reveal the nature of two near-Earth asteroid triples. Triple-asteroid system 2001 SN263 is characterized by a primary of ~10^13 kg, an inner satellite ~1% as massive orbiting at ~3 primary radii in ~0.7 days, and an outer satellite ~2.5% as massive orbiting at ~13 primary radii in ~6.2 days. 1994 CC is a smaller system with a primary of mass ~2.6 \times 10^11 kg and two satellites ~2% and ~1% as massive orbiting at distances of ~5.5 and ~19.5 primary radii. Their orbital periods are ~1.2 and ~8.4 days. Examination of resonant arguments shows that the satellites are not currently in a mean-motion resonance. Precession of the apses and nodes are detected in both systems (2001 SN263 inner body: d{\varpi}/dt ~1.1 deg/day, 1994 CC inner body: d{\varpi}/dt ~ -0.2 deg/day), which is in agreement with analytical predictions of the secular evolution due to mutually interacting orbits and primary oblateness. Nonzero mutual inclinations between the orbital planes of the satellites provide the best fits to the data in both systems (2001 SN263: ~14 degrees, 1994 CC: ~16 degrees). Our best-fit orbits are consistent with nearly circular motion, except for 1994 CC's outer satellite which has an eccentric orbit of e ~ 0.19. We examine several processes that can generate the observed eccentricity and inclinations, including the Kozai and evection resonances, past mean-motion resonance crossings, and close encounters with terrestrial planets. In particular, we find that close planetary encounters can easily excite the eccentricities and mutual inclinations of the satellites' orbits to the currently observed values.

We observed the near-Earth asteroid 2008 EV5 with the Arecibo and Goldstone planetary radars and the Very Long Baseline Array during December 2008. EV5 rotates retrograde and its overall shape is a 400 /pm 50 m oblate spheroid. The most prominent surface feature is a ridge parallel to the asteroid's equator that is broken by a concavity 150 m in diameter. Otherwise the asteroid's surface is notably smooth on decameter scales. EV5's radar and optical albedos are consistent with either rocky or stony-iron composition. The equatorial ridge is similar to structure seen on the rubble-pile near-Earth asteroid (66391) 1999 KW4 and is consistent with YORP spin-up reconfiguring the asteroid in the past. We interpret the concavity as an impact crater. Shaking during the impact and later regolith redistribution may have erased smaller features, explaining the general lack of decameter-scale surface structure.

We provide a generalized discussion of tidal evolution to arbitrary order in the expansion of the gravitational potential between two spherical bodies of any mass ratio. To accurately reproduce the tidal evolution of a system at separations less than five times the radius of the larger primary component, the tidal potential due to the presence of a smaller secondary component is expanded in terms of Legendre polynomials to arbitrary order rather than truncated at leading order as is typically done in studies of well-separated system like the Earth and Moon. The equations of tidal evolution including tidal torques, the changes in spin rates of the components, and the change in semimajor axis (orbital separation) are then derived for binary asteroid systems with circular and equatorial mutual orbits. Accounting for higher-order terms in the tidal potential serves to speed up the tidal evolution of the system leading to underestimates in the time rates of change of the spin rates, semimajor axis, and mean motion in the mutual orbit if such corrections are ignored. Special attention is given to the effect of close orbits on the calculation of material properties of the components, in terms of the rigidity and tidal dissipation function, based on the tidal evolution of the system. It is found that accurate determinations of the physical parameters of the system, e.g., densities, sizes, and current separation, are typically more important than accounting for higher-order terms in the potential when calculating material properties. In the scope of the long-term tidal evolution of the semimajor axis and the component spin rates, correcting for close orbits is a small effect, but for an instantaneous rate of change in spin rate, semimajor axis, or mean motion, the close-orbit correction can be on the order of tens of percent.

We present 3D numerical simulations of the early evolution of long-duration gamma-ray bursts in the collapsar scenario. Starting from the core-collapse of a realistic progenitor model, we follow the formation and evolution of a central black hole and centrifugally balanced disk. The dense, hot accretion disk produces freely-escaping neutrinos and is hydrodynamically unstable to clumping and to forming non-axisymmetric (m=1, 2) modes. We show that these spiral structures, which form on dynamical timescales, can efficiently transfer angular momentum outward and can drive the high required accretion rates (>=0.1-1 M_sun) for producing a jet. We utilise the smoothed particle hydrodynamics code, Gadget-2, modified to implement relevant microphysics, such as cooling by neutrinos, a plausible treatment approximating the central object and relativistic effects. Finally, we discuss implications of this scenario as a source of energy to produce relativistically beamed gamma-ray jets.

We present the first results from a Chandra survey of the central region of the Small Magellanic Cloud. We detect a total of 122 sources down to a limiting luminosity of ~4.3x10^{33} erg\s (corrected for Galactic N_{H}), which is ~10 times lower than in any previous survey of the SMC. We identify 20 candidate transient sources: eighteen previously known sources in this area which are not detected in our observations, and two new bright sources. The spectral parameters of the brightest sources indicate that they are X-ray binary pulsars. The high spatial resolution of Chandra allows us to initially identify optical counterparts for 35 sources, 13 of which are new identifications.

We numerically investigate the migration of dust particles with initial orbits close to those of the numbered asteroids, observed trans-Neptunian objects, and Comet Encke. The fraction of silicate asteroidal particles that collided with the Earth during their lifetime varied from 1.1% for 100 micron particles to 0.008% for 1 micron particles. Almost all asteroidal particles with diameter d>4 microns collided with the Sun. The peaks in the migrating asteroidal dust particles' semi-major axis distribution at the n:(n+1) resonances with Earth and Venus and the gaps associated with the 1:1 resonances with these planets are more pronounced for larger particles. The probability of collisions of cometary particles with the Earth is smaller than for asteroidal particles, and this difference is greater for larger particles.

We numerically investigate the migration of dust particles with initial orbits close to those of the numbered asteroids. The fraction of silicate particles that collided with the Earth during their lifetimes varied from 0.2% for 40 micron particles to 0.008% for 1 micron particles. Almost all particles with diameter d>4 microns collided with the Sun. The peaks in the migrating asteroidal dust particles' semi-major axis distribution at the n/(n+1) resonances with Earth and Venus and the gaps associated with the 1:1 resonances with these planets are more pronounced for larger particles.