Using the novel Magnetospheric Multiscale (MMS) mission data accumulated during the 2019 MMS Solar Wind Turbulence Campaign, we calculate the Taylor microscale $(\lambda_{\mathrm{T}})$ of the turbulent magnetic field in the solar wind. The Taylor microscale represents the onset of dissipative processes in classical turbulence theory. An accurate estimation of Taylor scale from spacecraft data is, however, usually difficult due to low time cadence, the effect of time decorrelation, and other factors. Previous reports were based either entirely on the Taylor frozen-in approximation, which conflates time dependence, or that were obtained using multiple datasets, which introduces sample-to-sample variation of plasma parameters, or where inter-spacecraft distance were larger than the present study. The unique configuration of linear formation with logarithmic spacing of the 4 MMS spacecraft, during the campaign, enables a direct evaluation of the $\lambda_{\mathrm{T}}$ from a single dataset, independent of the Taylor frozen-in approximation. A value of $\lambda_{\mathrm{T}} \approx 7000 \, \mathrm{km}$ is obtained, which is about 3 times larger than the previous estimates.

Both kinetic instabilities and strong turbulence have potential to impact the behavior of space plasmas. To assess effects of these two processes we compare results from a 3 dimensional particle-in-cell (PIC) simulation of collisionless plasma turbulence against observations by the MMS spacecraft in the terrestrial magnetosheath and by the Wind spacecraft in the solar wind. The simulation develops coherent structures and anisotropic ion velocity distributions that can drive micro-instabilities. Temperature-anisotropy driven instability growth rates are compared with inverse nonlinear turbulence time scales. Large growth rates occur near coherent structures; nevertheless linear growth rates are, on average, substantially less than the corresponding nonlinear rates. This result casts some doubt on the usual basis for employing linear instability theory, and raises questions as to why the linear theory appears to work in limiting plasma excursions in anisotropy and plasma beta.

To investigate the use of saliency-map analysis to aid in searches for transient signals, such as fast radio bursts and individual pulses from radio pulsars. We aim to demonstrate that saliency maps provide the means to understand predictions from machine learning algorithms and can be implemented in piplines used to search for transient events. We have implemented a new deep learning methodology to predict whether or not any segment of the data contains a transient event. The algorithm has been trained using real and simulated data sets. We demonstrate that the algorithm is able to identify such events. The output results are visually analysed via the use of saliency maps. We find that saliency maps can produce an enhanced image of any transient feature without the need for de-dispersion or removal of radio frequency interference. Such maps can be used to understand which features in the image were used in making the machine learning decision and to visualise the transient event. Even though the algorithm reported here was developed to demonstrate saliency-map analysis, we have detected, in archival data, a single burst event with dispersion measure of $41$\,cm$^{-3}$pc that is not associated with any currently known pulsar.

A familiar problem in space and astrophysical plasmas is to understand how dissipation and heating occurs. These effects are often attributed to the cascade of broadband turbulence which transports energy from large scale reservoirs to small scale kinetic degrees of freedom. When collisions are infrequent, local thermodynamic equilibrium is not established. In this case the final stage of energy conversion becomes more complex than in the fluid case, and both pressure-dilatation and pressure strain interactions (Pi-D $\equiv -\Pi_{ij} D_{ij}$) become relevant and potentially important. Pi-D in plasma turbulence has been studied so far primarily using simulations. The present study provides a statistical analysis of Pi-D in the Earth's magnetosheath using the unique measurement capabilities of the Magnetospheric Multiscale (MMS) mission. We find that the statistics of Pi-D in this naturally occurring plasma environment exhibit strong resemblance to previously established fully kinetic simulations results. The conversion of energy is concentrated in space and occurs near intense current sheets, but not within them. This supports recent suggestions that the chain of energy transfer channels involves regional, rather than pointwise, correlations.

Thanks to the VISTA Variables in the Via Lactea (VVV) ESO Public Survey it is now possible to explore a large number of objects in those regions. This paper addresses the variability analysis of all VVV point sources having more than 10 observations in VVVDR4 using a novel approach. In total, the near-IR light curves of 288,378,769 sources were analysed using methods developed in the New Insight Into Time Series Analysis project. As a result, we present a complete sample having 44, 998, 752 variable star candidates (VVV-CVSC), which include accurate individual coordinates, near-IR magnitudes (ZYJHKs), extinctions A(Ks), variability indices, periods, amplitudes, among other parameters to assess the science. Unfortunately, a side effect of having a highly complete sample, is also having a high level of contamination by non-variable (contamination ratio of non-variables to variables is slightly over 10:1). To deal with this, we also provide some flags and parameters that can be used by the community to de-crease the number of variable candidates without heavily decreasing the completeness of the sample. In particular, we cross-identified 339,601 of our sources with Simbad and AAVSO databases, which provide us with information for these objects at other wavelegths. This sub-sample constitutes a unique resource to study the corresponding near-IR variability of known sources as well as to assess the IR variability related with X-ray and Gamma-Ray sources. On the other hand, the other 99.5% sources in our sample constitutes a number of potentially new objects with variability information for the heavily crowded and reddened regions of the Galactic Plane and Bulge. The present results also provide an important queryable resource to perform variability analysis and to characterize ongoing and future surveys like TESS and LSST.

A large number of observations from the Parkes 64\,m-diameter radio telescope, recorded with high time resolution, are publicly available. We have re-processed all of the observations obtained during the first four years (from 1997 to 2001) of the Parkes Multibeam receiver system in order to identify transient events and have built a database that records the 568,736,756 pulse candidates generated during this search. We have discovered a new fast radio burst (FRB), FRB~010305, with a dispersion measure of 350$\pm$5\,\,cm$^{-3}\,$pc and explored why so few FRBs have been discovered in data prior to 2001. After accounting for the dispersion smearing across the channel bandwidth and the sky regions surveyed, the number of FRBs is found to be consistent with model predictions. We also present five single pulse candidates from unknown sources, but with Galactic dispersion measures. We extract a diverse range of sources from the database, which can be used, for example, as a training set of data for new software being developed to search for FRBs in the presence of radio frequency interference.

Using in situ data, accumulated in the turbulent magnetosheath by the Magnetospheric Multiscale (MMS) Mission, we report a statistical study of magnetic field curvature and discuss its role in the turbulent space plasmas. Consistent with previous simulation results, the Probability Distribution Function (PDF) of the curvature is shown to have distinct power-law tails for both high and low value limits. We find that the magnetic-field-line curvature is intermittently distributed in space. High curvature values reside near weak magnetic-field regions, while low curvature values are correlated with small magnitude of the force acting normal to the field lines. A simple statistical treatment provides an explanation for the observed curvature distribution. This novel statistical characterization of magnetic curvature in space plasma provides a starting point for assessing, in a turbulence context, the applicability and impact of particle energization processes, such as curvature drift, that rely on this fundamental quantity.

We describe 14 years of public data from the Parkes Pulsar Timing Array (PPTA), an ongoing project that is producing precise measurements of pulse times of arrival from 26 millisecond pulsars using the 64-m Parkes radio telescope with a cadence of approximately three weeks in three observing bands. A comprehensive description of the pulsar observing systems employed at the telescope since 2004 is provided, including the calibration methodology and an analysis of the stability of system components. We attempt to provide full accounting of the reduction from the raw measured Stokes parameters to pulse times of arrival to aid third parties in reproducing our results. This conversion is encapsulated in a processing pipeline designed to track provenance. Our data products include pulse times of arrival for each of the pulsars along with an initial set of pulsar parameters and noise models. The calibrated pulse profiles and timing template profiles are also available. These data represent almost 21,000 hrs of recorded data spanning over 14 years. After accounting for processes that induce time-correlated noise, 22 of the pulsars have weighted root-mean-square timing residuals of < 1 ${\mu}$s in at least one radio band. The data should allow end users to quickly undertake their own gravitational-wave analyses (for example) without having to understand the intricacies of pulsar polarisation calibration or attain a mastery of radio-frequency interference mitigation as is required when analysing raw data files.

The first-order Fermi acceleration of electrons requires an injection of electrons into a mildly relativistic energy range. However, the mechanism of injection has remained a puzzle both in theory and observation. We present direct evidence for a novel stochastic shock drift acceleration theory for the injection obtained with Magnetospheric Multiscale (MMS) observations at Earth's bow shock. The theoretical model can explain electron acceleration to mildly relativistic energies at high-speed astrophysical shocks, which may provide a solution to the long-standing issue of electron injection.

The inner parsec of our Galaxy contains tens of Wolf-Rayet stars whose powerful outflows are constantly interacting while filling the region with hot, diffuse plasma. Theoretical models have shown that, in some cases, the collision of stellar winds can generate cold, dense material in the form of clumps. However, their formation process and properties are not well understood yet. In this work we present, for the first time, a statistical study of the clump formation process in unstable wind collisions. We study systems with dense outflows (${\sim}10^{-5}\rm\ M_{\odot}\ yr^{-1}$), wind speeds of $500$-$1500\rm\ km\ s^{-1}$, and stellar separations of ${\sim}20$-$200\rm\ au$. We develop 3D high resolution hydrodynamical simulations of stellar wind collisions with the adaptive-mesh refinement grid-based code Ramses. We aim to characterise the initial properties of clumps that form through hydrodynamic instabilities, mostly via the non-linear thin shell instability (NTSI). Our results confirm that more massive clumps are formed in systems whose winds are close to the transition between the radiative and adiabatic regimes. Increasing either the wind speed or the degree of asymmetry increases the dispersion of the clump mass and ejection speed distributions. Nevertheless, the most massive clumps are very light (${\sim}10^{-3}$-$10^{-2}\rm\ M_{\oplus}$), about three orders of magnitude less massive than theoretical upper limits. Applying these results to the Galactic Centre we find that clumps formed through the NTSI should not be heavy enough either to affect the thermodynamic state of the region or to survive for long enough to fall onto the central super-massive black hole.

IRS~13E is an enigmatic compact group of massive stars located in projection only 3.6 arcseconds away from Sgr A*. This group has been suggested to be bounded by an intermediate-mass black hole (IMBH). We present a multi-wavelength study of the group and its interplay with the environment. Based on Chandra observations, we find the X-ray spectrum of IRS~13E can be well characterized by an optically thin thermal plasma. The emission peaks between two strongly mass-losing Wolf-Rayet stars of the group. These properties can be reasonably well reproduced by simulated colliding winds of these two stars. However, this scenario under-predicts the X-ray intensity in outer regions. The residual emission likely results from the ram-pressure confinement of the IRS~13E group wind by the ambient medium and is apparently associated with a shell-like warm gas structure seen in Pa-alpha and in ALMA observations. These latter observations also show strongly peaked thermal emission with unusually large velocity spread between the two stars. These results indicate that the group is colliding with the bar of the dense cool gas mini-spiral around Sgr A*. The extended X-ray morphology of IRS~13E and its association with the bar further suggest that the group is physically much farther away than the projected distance from Sgr A*. The presence of an IMBH, while favorable to keep the stars bound together, is not necessary to explain the observed stellar and gas properties of IRS~13E.

It is widely accepted that dark matter contributes about a quarter of the critical mass-energy density in our Universe. The nature of dark matter is currently unknown, with the mass of possible constituents spanning nearly one hundred orders of magnitude. The ultralight scalar field dark matter, consisting of extremely light bosons with $m \sim 10^{-22}$ eV and often called "fuzzy" dark matter, provides intriguing solutions to some challenges at sub-Galactic scales for the standard cold dark matter model. As shown by Khmelnitsky and Rubakov, such a scalar field in the Galaxy would produce an oscillating gravitational potential with nanohertz frequencies, resulting in periodic variations in the times of arrival of radio pulses from pulsars. The Parkes Pulsar Timing Array (PPTA) has been monitoring 20 millisecond pulsars at two to three weeks intervals for more than a decade. In addition to the detection of nanohertz gravitational waves, PPTA offers the opportunity for direct searches for fuzzy dark matter in an astrophysically feasible range of masses. We analyze the latest PPTA data set which includes timing observations for 26 pulsars made between 2004 and 2016. We perform a search in this data set for evidence of ultralight dark matter in the Galaxy using Bayesian and Frequentist methods. No statistically significant detection has been made. We therefore place upper limits on the local dark matter density. Our limits, improving on previous searches by a factor of two to five, constrain the dark matter density of ultralight bosons with $m \leq 10^{-23}$ eV to be below $6\,\text{GeV}\,\text{cm}^{-3}$ with 95\% confidence in the Earth neighborhood. Finally, we discuss the prospect of probing the astrophysically favored mass range $m \gtrsim 10^{-22}$ eV with next-generation pulsar timing facilities.

The central super-massive black hole of the Milky Way, Sgr A*, accretes at a very low rate making it a very underluminous galactic nucleus. Despite the tens of Wolf-Rayet stars present within the inner parsec supplying ${\sim}10^{-3}\rm\ M_{\odot}\ yr^{-1}$ in stellar winds, only a negligible fraction of this material ($<10^{-4}$) ends up being accreted onto Sgr A*. The recent discovery of cold gas (${\sim}10^4\rm\ K$) in its vicinity raised questions about how such material could settle in the hostile (${\sim}10^7\rm\ K$) environment near Sgr A*. In this work we show that the system of mass-losing stars blowing winds can naturally account for both the hot, inefficient accretion flow, as well as the formation of a cold disk-like structure. We run hydrodynamical simulations using the grid-based code Ramses starting as early in the past as possible to observe the state of the system at the present time. Our results show that the system reaches a quasi-steady state in about ${\sim}500\rm\ yr$ with material being captured at a rate of ${\sim}10^{-6}\rm\ M_{\odot}\ yr^{-1}$ at scales of ${\sim}10^{-4}\rm\ pc$, consistent with the observations and previous models. However, on longer timescales ($\gtrsim3000\rm\ yr$) the material accumulates close to the black hole in the form of a disk. Considering the duration of the Wolf-Rayet phase (${\sim}10^5\rm\ yr$), we conclude that this scenario likely has already happened, and could be responsible for the more active past of Sgr A*, and/or its current outflow. We argue that the hypothesis of the mass-losing stars being the main regulator of the activity of the black hole deserves further consideration.

We present analysis of more than one hundred large-amplitude bipolar electrostatic structures in a quasi-perpendicular supercritical Earth's bow shock crossing, measured by the Magnetospheric Multiscale spacecraft. The occurrence of the bipolar structures is shown to be tightly correlated with magnetic field gradients in the shock transition region. The bipolar structures have negative electrostatic potentials and spatial scales of a few Debye lengths. The bipolar structures propagate highly oblique to the shock normal with velocities (in the plasma rest frame) of the order of the ion-acoustic velocity. We argue that the bipolar structures are ion phase space holes produced by the two-stream instability between incoming and reflected ions. This is the first identification of the ion two-stream instability in collisionless shocks. The implications for electron acceleration are discussed.

We have constructed a new timescale, TT(IPTA16), based on observations of radio pulsars presented in the first data release from the International Pulsar Timing Array (IPTA). We used two analysis techniques with independent estimates of the noise models for the pulsar observations and different algorithms for obtaining the pulsar timescale. The two analyses agree within the estimated uncertainties and both agree with TT(BIPM17), a post-corrected timescale produced by the Bureau International des Poids et Mesures (BIPM). We show that both methods could detect significant errors in TT(BIPM17) if they were present. We estimate the stability of the atomic clocks from which TT(BIPM17) is derived using observations of four rubidium fountain clocks at the US Naval Observatory. Comparing the power spectrum of TT(IPTA16) with that of these fountain clocks suggests that pulsar-based timescales are unlikely to contribute to the stability of the best timescales over the next decade, but they will remain a valuable independent check on atomic timescales. We also find that the stability of the pulsar-based timescale is likely to be limited by our knowledge of solar-system dynamics, and that errors in TT(BIPM17) will not be a limiting factor for the primary goal of the IPTA, which is to search for the signatures of nano-Hertz gravitational waves.

In this paper, we describe the International Pulsar Timing Array second data release, which includes recent pulsar timing data obtained by three regional consortia: the European Pulsar Timing Array, the North American Nanohertz Observatory for Gravitational Waves, and the Parkes Pulsar Timing Array. We analyse and where possible combine high-precision timing data for 65 millisecond pulsars which are regularly observed by these groups. A basic noise analysis, including the processes which are both correlated and uncorrelated in time, provides noise models and timing ephemerides for the pulsars. We find that the timing precisions of pulsars are generally improved compared to the previous data release, mainly due to the addition of new data in the combination. The main purpose of this work is to create the most up-to-date IPTA data release. These data are publicly available for searches for low-frequency gravitational waves and other pulsar science.

The Parkes Pulsar Timing Array (PPTA) project monitors two dozen millisecond pulsars (MSPs) in order to undertake a variety of fundamental physics experiments using the Parkes 64m radio telescope. Since June 2017 we have been undertaking commensal searches for fast radio bursts (FRBs) during the MSP observations. Here, we report the discovery of four FRBs (171209, 180309, 180311 and 180714). The detected events include an FRB with the highest signal-to-noise ratio ever detected at the Parkes observatory, which exhibits unusual spectral properties. All four FRBs are highly polarized. We discuss the future of commensal searches for FRBs at Parkes.

The degree of deviation from the thermodynamic equilibrium in the ion velocity distribution functions (VDFs), measured by the Magnetospheric Multiscale (MMS) mission in the Earth's turbulent magnetosheath, is quantitatively investigated. Taking advantage of MMS ion data, having a resolution never reached before in space missions, and of the comparison with Vlasov-Maxwell simulations, this analysis aims at relating any deviation from Maxwellian equilibrium to typical plasma parameters. Correlations of the non-Maxwellian features with plasma quantities such as electric fields, ion temperature, current density and ion vorticity are very similar in both magnetosheath data and numerical experiments, and suggest that distortions in the ion VDFs occur close to (but not exactly at) peaks in current density and ion temperature. Similar results have also been found during a magnetopause crossing by MMS. This work could help clarifying the origin of distortion of the ion VDFs in space plasmas.

Cosmic-ray acceleration has been a long-standing mystery and despite more than a century of study, we still do not have a complete census of acceleration mechanisms. The collision of strong stellar winds in massive binary systems creates powerful shocks, which have been expected to produce high-energy cosmic-rays through Fermi acceleration at the shock interface. The accelerated particles should collide with stellar photons or ambient material, producing non-thermal emission observable in X-rays and gamma-rays. The supermassive binary star eta Carinae drives the strongest colliding wind shock in the solar neighborhood. Observations with non-focusing high-energy observatories indicate a high energy source near eta Carinae, but have been unable to conclusively identify eta Carinae as the source because of their relatively poor angular resolution. Here we present the first direct focussing observations of the non-thermal source in the extremely hard X-ray band, which is found to be spatially coincident with the star within several arc-seconds. These observations show that the source of non-thermal X-rays varies with the orbital phase of the binary, and that the photon index of the emission is similar to that derived through analysis of the gamma-ray spectrum. This is conclusive evidence that the high-energy emission indeed originates from non-thermal particles accelerated at colliding wind shocks.

We report the discovery of a new fast radio burst (FRB), FRB~010312, in archival data from a 1.4\,GHz survey of the Magellanic Clouds using the multibeam receiver on the Parkes 64\,m-diameter radio telescope. These data sets include the Lorimer burst (FRB~010724), which it pre-dates and which we also re-detect. The new burst has a much higher dispersion measure of 1187\,cm$^{-3}$pc. The burst is one of the broadest found to date, the second earliest FRB known, and the ninth FRB discovered with a dispersion measure larger than 1000\,cm$^{-3}\,$pc. Our discovery indicates that there are likely to be more burst events still to be found in the existing Parkes data archive.

We report a comprehensive analysis of the global spectrophotometric properties of Ceres using Dawn Framing Camera images collected from April to June 2015 during the RC3 and Survey mission phases. The single-scattering albedo of Ceres at 555 nm is 0.14$\pm$0.04, the geometric albedo is 0.096$\pm$0.006, and the Bond albedo is 0.037$\pm$0.002. The asymmetry factors calculated from the best-fit two-term Henyey-Greenstein (HG) single-particle phase function (SPPF) show a wavelength dependence, suggesting that the phase reddening of Ceres is dominated by single-particle scattering rather than multiple scattering or small-scale surface roughness. The Hapke roughness parameter of Ceres is derived to be 20$^\circ\pm$6$^\circ$ with no wavelength dependence. The phase function of Ceres shows appreciably strong scattering around 90$^\circ$ phase angle that cannot be fitted with a single-term HG SPPF, suggesting possible stronger forward scattering than other asteroids previously analyzed with spacecraft data. We speculate that such a scattering characteristic of Ceres might be related to its unique surface composition. We grouped the reflectance data into a 1$^\circ$ latitude-longitude grid and fitted each grid independently to study the spatial variations of photometric properties. The albedo and color maps are consistent with previous studies. The SPPF over the surface of Ceres shows stronger backscattering associated with lower albedo and vice versa, consistent with the general trend among asteroids. The Hapke roughness parameter does not vary much across the surface of Ceres, except for the ancient Vendimia Planitia region that has a slightly higher roughness. Based on the wavelength dependence of the SPPF of Ceres, we hypothesize that its regolith grains either contain a considerable fraction of $\lessapprox\mu$m-sized particles, or are strongly affected by internal scatterers of this size.

The surface reflectance of planetary regoliths may increase dramatically towards zero phase angle, a phenomenon known as the opposition effect (OE). Two physical processes that are thought to be the dominant contributors to the brightness surge are shadow hiding (SH) and coherent backscatter (CB). The occurrence of shadow hiding in planetary regoliths is self-evident, but it has proved difficult to unambiguously demonstrate CB from remote sensing observations. One prediction of CB theory is the wavelength dependence of the OE angular width. The Dawn spacecraft observed the OE on the surface of dwarf planet Ceres. We characterize the OE over the resolved surface, including the bright Cerealia Facula, and to find evidence for SH and/or CB. We analyze images of the Dawn framing camera by means of photometric modeling of the phase curve. We find that the OE of most of the investigated surface has very similar characteristics, with an enhancement factor of 1.4 and a FWHM of 3{\deg} (broad OE). A notable exception are the fresh ejecta of the Azacca crater, which display a very narrow brightness enhancement that is restricted to phase angles $< 0.5${\deg} (narrow OE); suggestively, this is in the range in which CB is thought to dominate. We do not find a wavelength dependence for the width of the broad OE, and lack the data to investigate the dependence for the narrow OE. The prediction of a wavelength-dependent CB width is rather ambiguous. The zero-phase observations allow us to determine Ceres' visible geometric albedo as $p_V = 0.094 \pm 0.005$. A comparison with other asteroids suggests that Ceres' broad OE is typical for an asteroid of its spectral type, with characteristics that are primarily linked to surface albedo. Our analysis suggests that CB may occur on the dark surface of Ceres in a highly localized fashion.

We study spectral features of ion velocity and magnetic field correlations in the solar wind and in the magnetosheath using data from the Magnetospheric Multi-Scale (MMS) spacecraft. High resolution MMS observations enable the study of transition of these correlations between their magnetofluid character at larger scales into the sub-proton kinetic range, previously unstudied in spacecraft data. Cross-helicity, angular alignment and energy partitioning is examined over a suit- able range of scales, employing measurements based on the Taylor frozen-in approximation as well as direct two-spacecraft correlation measurements. The results demonstrate signatures of alignment at large scales. As kinetic scales are approached, the alignment between v and b is destroyed by demagnetization of protons.

Protons (ionized hydrogen) in the solar wind frequently exhibit distinct temperatures ($T_{\perp p}$ and $T_{\parallel p}$) perpendicular and parallel to the plasma's background magnetic-field. Numerous prior studies of the interplanetary solar-wind have shown that, as plasma beta ($\beta_{\parallel p}$) increases, a narrower range of temperature-anisotropy ($R_p\equiv T_{\perp p}\,/\,T_{\parallel p}$) values is observed. Conventionally, this effect has been ascribed to the actions of kinetic microinstabilities. This study is the first to use data from the Magnetospheric Multiscale Mission (MMS) to explore such $\beta_{\parallel p}$-dependent limits on $R_p$ in Earth's magnetosheath. The distribution of these data across the $(\beta_{\parallel p},R_p)$-plane reveals limits on both $R_p>1$ and $R_p<1$. Linear Vlasov theory is used to compute contours of constant growth-rate for the ion-cyclotron, mirror, parallel-firehose, and oblique-firehose instabilities. These instability thresholds closely align with the contours of the data distribution, which suggests a strong association of instabilities with extremes of ion temperature anisotropy in the magnetosheath. The potential for instabilities to regulate temperature anisotropy is discussed.

We have observed the remnant of supernova SN~1987A (SNR~1987A), located in the Large Magellanic Cloud (LMC), to search for periodic and/or transient radio emission with the Parkes 64\,m-diameter radio telescope. We found no evidence of a radio pulsar in our periodicity search and derived 8$\sigma$ upper bounds on the flux density of any such source of $31\,\mu$Jy at 1.4~GHz and $21\,\mu$Jy at 3~GHz. Four candidate transient events were detected with greater than $7\sigma$ significance, with dispersion measures (DMs) in the range 150 to 840\,cm$^{-3}\,$pc. For two of them, we found a second pulse at slightly lower significance. However, we cannot at present conclude that any of these are associated with a pulsar in SNR~1987A. As a check on the system, we also observed PSR~B0540$-$69, a young pulsar which also lies in the LMC. We found eight giant pulses at the DM of this pulsar. We discuss the implications of these results for models of the supernova remnant, neutron star formation and pulsar evolution.

Future remote sensing of exoplanets will be enhanced by a thorough investigation of our solar system Ice Giants (Neptune-size planets). What can the configuration of the magnetic field tell us (remotely) about the interior, and what implications does that field have for the structure of the magnetosphere; energy input into the atmosphere, and surface geophysics (for example surface weathering of satellites that might harbour sub-surface oceans). How can monitoring of auroral emission help inform future remote observations of emission from exoplanets? Our Solar System provides the only laboratory in which we can perform in-situ experiments to understand exoplanet formation, dynamos, systems and magnetospheres.

Data from the NASA Magnetospheric Multiscale (MMS) mission are used to investigate asymmetric magnetic reconnection at the dayside boundary between the Earth's magnetosphere and the solar wind (the magnetopause). High-resolution measurements of plasmas, electric and magnetic fields, and waves are used to identify highly localized (~15 electron Debye lengths) standing wave structures with large electric-field amplitudes (up to 100 mV/m). These wave structures are associated with spatially oscillatory dissipation, which appears as alternatingly positive and negative values of J dot E (dissipation). For small guide magnetic fields the wave structures occur in the electron stagnation region at the magnetosphere edge of the EDR. For larger guide fields the structures also occur near the reconnection x-line. This difference is explained in terms of channels for the out-of-plane current (agyrotropic electrons at the stagnation point and guide-field-aligned electrons at the x-line).

The Galactic Centre is a hotbed of astrophysical activity, with the injection of wind material from $\sim$30 massive Wolf-Rayet (WR) stars orbiting within 12 arcsec of the supermassive black hole (SMBH) playing an important role. Hydrodynamic simulations of such colliding and accreting winds produce a complex density and temperature structure of cold wind material shocking with the ambient medium, creating a large reservoir of hot, X-ray-emitting gas. This work aims to confront the 3 Ms of Chandra X-ray Visionary Program observations of this diffuse emission by computing the X-ray emission from these hydrodynamic simulations of the colliding WR winds, amid exploring a variety of SMBH feedback mechanisms. The major success of the model is that it reproduces the spectral shape from the 2-5 arcsec ring around the SMBH, where most of the stellar wind material that is ultimately captured by Sgr A* is shock-heated and thermalised. This naturally explains that the hot gas comes from colliding WR winds, and that the wind speeds of these stars are, in general, well constrained. The flux level of these spectra, as well as 12 $\times$ 12-arcsec$^2$ images of 4-9 keV, shows that the X-ray flux is tied to the SMBH feedback strength; stronger feedback clears out more hot gas, thereby decreasing the thermal X-ray emission. The model in which Sgr A* produced an intermediate-strength outflow during the last few centuries best matches the observations to within about 10 per cent, showing that SMBH feedback is required to interpret the X-ray emission in this region.

Plasma turbulence is investigated using high-resolution ion velocity distributions measured by the Magnetospheric Multiscale Mission (MMS) in the Earth's magnetosheath. The particle distribution is highly structured, suggesting a cascade-like process in velocity space. This complex velocity space structure is investigated using a three-dimensional Hermite transform that reveals a power law distribution of moments. In analogy to hydrodynamics, a Kolmogorov approach leads directly to a range of predictions for this phase-space cascade. The scaling theory is in agreement with observations, suggesting a new path for the study of plasma turbulence in weakly collisional space and astrophysical plasmas.

We compute the thermal X-ray emission from hydrodynamic simulations of the 30 Wolf-Rayet (WR) stars orbiting within a parsec of Sgr A$^*$, with the aim of interpreting the Chandra X-ray observations of this region. The model well reproduces the spectral shape of the observations, indicating that the shocked WR winds are the dominant source of this thermal emission. The model X-ray flux is tied to the strength of the Sgr A$^*$ outflow, which clears out hot gas from the vicinity of Sgr A$^*$. A moderate outflow best fits the present-day observations, even though this supermassive black hole (SMBH) outflow ended $\sim$100 yr ago.

360-degree videos are a new type of movie that renders over all 4$\pi$ steradian. Video sharing sites such as YouTube now allow this unique content to be shared via virtual reality (VR) goggles, hand-held smartphones/tablets, and computers. Creating 360$^\circ$ videos from astrophysical simulations is not only a new way to view these simulations as you are immersed in them, but is also a way to create engaging content for outreach to the public. We present what we believe is the first 360$^\circ$ video of an astrophysical simulation: a hydrodynamics calculation of the central parsec of the Galactic centre. We also describe how to create such movies, and briefly comment on what new science can be extracted from astrophysical simulations using 360$^\circ$ videos.

We report on the first multi-color precision light curve of the bright Wolf-Rayet binary $\gamma^2$ Velorum, obtained over six months with the nanosatellites in the BRITE- Constellation fleet. In parallel, we obtained 488 high-resolution optical spectra of the system. In this first report on the datasets, we revise the spectroscopic orbit and report on the bulk properties of the colliding winds. We find a dependence of both the light curve and excess emission properties that scales with the inverse of the binary separation. When analyzing the spectroscopic properties in combination with the photometry, we find that the phase dependence is caused only by excess emission in the lines, and not from a changing continuum. We also detect a narrow, high-velocity absorption component from the He I $\lambda$5876 transition, which appears twice in the orbit. We calculate smoothed-particle hydrodynamical simulations of the colliding winds and can accurately associate the absorption from He I to the leading and trailing arms of the wind shock cone passing tangentially through our line of sight. The simulations also explain the general strength and kinematics of the emission excess observed in wind lines such as C III $\lambda$5696 of the system. These results represent the first in a series of investigations into the winds and properties of $\gamma^2$ Velorum through multi-technique and multi-wavelength observational campaigns.

Pulsar positions can be measured with high precision using both pulsar timing methods and very-long-baseline interferometry (VLBI). Pulsar timing positions are referenced to a solar-system ephemeris, whereas VLBI positions are referenced to distant quasars. Here we compare pulsar positions from published VLBI measurements with those obtained from pulsar timing data from the Nanshan and Parkes radio telescopes in order to relate the two reference frames. We find that the timing positions differ significantly from the VLBI positions (and also differ between different ephemerides). A statistically significant change in the obliquity of the ecliptic of $2.16\pm0.33$\,mas is found for the JPL ephemeris DE405, but no significant rotation is found in subsequent JPL ephemerides. The accuracy with which we can relate the two frames is limited by the current uncertainties in the VLBI reference source positions and in matching the pulsars to their reference source. Not only do the timing positions depend on the ephemeris used in computing them, but also different segments of the timing data lead to varying position estimates. These variations are mostly common to all ephemerides, but slight changes are seen at the 10$\mu$as level between ephemerides.

We study the spectrophotometric properties of dwarf planet Ceres in the VIS-IR spectral range by means of hyper-spectral images acquired by the VIR imaging spectrometer on board the NASA Dawn mission. Disk-resolved observations with a phase angle within the $7^{\circ}<\alpha<132^{\circ}$ interval were used to characterize Ceres' phase curve in the 0.465-4.05 $\mu$m spectral range. Hapke's model was applied to perform the photometric correction of the dataset, allowing us to produce albedo and color maps of the surface. The $V$-band magnitude phase function of Ceres was fitted with both the classical linear model and H-G formalism. The single-scattering albedo and the asymmetry parameter at 0.55$\mu$m are $w=0.14\pm0.02$ and $\xi=-0.11\pm0.08$, respectively (two-lobe Henyey-Greenstein phase function); the modeled geometric albedo is $0.094\pm0.007$; the roughness parameter is $\bar{\theta}=29^{\circ}\pm6^{\circ}$. Albedo maps indicate small variability on a global scale with an average reflectance of $0.034 \pm 0.003$. Isolated areas such as the Occator bright spots, Haulani, and Oxo show an albedo much higher than average. We measure a significant spectral phase reddening, and the average spectral slope of Ceres' surface after photometric correction is $1.1\%k\AA^{-1}$ and $0.85\%k\AA^{-1}$ at VIS and IR wavelengths, respectively. Broadband color indices are $V-R=0.38\pm0.01$ and $R-I=0.33\pm0.02$. H-G modeling of the $V$-band magnitude phase curve for $\alpha<30^{\circ}$ gives $H=3.14\pm0.04$ and $G=0.10\pm0.04$, while the classical linear model provides $V(1,1,0^{\circ})=3.48\pm0.03$ and $\beta=0.036\pm0.002$. The comparison with spectrophotometric properties of other minor bodies indicates that Ceres has a less back-scattering phase function and a slightly higher albedo than comets and C-type objects. However, the latter represents the closest match in the usual asteroid taxonomy.

MMS observations recently confirmed that crescent-shaped electron velocity distributions in the plane perpendicular to the magnetic field occur in the electron diffusion region near reconnection sites at Earth's magnetopause. In this paper, we re-examine the origin of the crescent-shaped distributions in the light of our new finding that ions and electrons are drifting in opposite directions when displayed in magnetopause boundary-normal coordinates. Therefore, ExB drifts cannot cause the crescent shapes. We performed a high-resolution multi-scale simulation capturing sub-electron skin depth scales. The results suggest that the crescent-shaped distributions are caused by meandering orbits without necessarily requiring any additional processes found at the magnetopause such as the highly asymmetric magnetopause ambipolar electric field. We use an adiabatic Hamiltonian model of particle motion to confirm that conservation of canonical momentum in the presence of magnetic field gradients causes the formation of crescent shapes without invoking asymmetries or the presence of an ExB drift. An important consequence of this finding is that we expect crescent-shaped distributions also to be observed in the magnetotail, a prediction that MMS will soon be able to test.

We present a global spectrophotometric characterization of the Ceres surface using Dawn Framing Camera (FC) images. We identify the photometric model that yields the best results for photometrically correcting images. Corrected FC images acquired on approach to Ceres were assembled into global maps of albedo and color. Generally, albedo and color variations on Ceres are muted. The albedo map is dominated by a large, circular feature in Vendimia Planitia, known from HST images (Li et al., 2006), and dotted by smaller bright features mostly associated with fresh-looking craters. The dominant color variation over the surface is represented by the presence of "blue" material in and around such craters, which has a negative spectral slope over the visible wavelength range when compared to average terrain. We also mapped variations of the phase curve by employing an exponential photometric model, a technique previously applied to asteroid Vesta (Schr\"oder et al., 2013b). The surface of Ceres scatters light differently from Vesta in the sense that the ejecta of several fresh-looking craters may be physically smooth rather than rough. High albedo, blue color, and physical smoothness all appear to be indicators of youth. The blue color may result from the desiccation of ejected material that is similar to the phyllosilicates/water ice mixtures in the experiments of Poch et al. (2016). The physical smoothness of some blue terrains would be consistent with an initially liquid condition, perhaps as a consequence of impact melting of subsurface water ice. We find red terrain (positive spectral slope) near Ernutet crater, where De Sanctis et al. (2017) detected organic material. The spectrophotometric properties of the large Vendimia Planitia feature suggest it is a palimpsest, consistent with the Marchi et al. (2016) impact basin hypothesis. The central bright area in Occator crater, Cerealia...

During its year-long orbital mission, the Dawn spacecraft has mapped the surface of main-belt asteroid Vesta multiple times at different spatial resolutions and illumination and viewing angles. The onboard Framing Camera has acquired thousands of clear filter and narrow band images, which, with the availability of high-resolution global shape models, allows for a photometric characterization of the surface in unprecedented detail. We analyze clear filter images to retrieve the photometric properties of the regolith. In the first part of the paper we evaluate different photometric models for the global average. In the second part we use these results to study variations in albedo and steepness of the phase curve over the surface. Maps of these two photometric parameters show large scale albedo variations, which appear to be associated with compositional differences. They also reveal the location of photometrically extreme terrains, where the phase curve is unusually shallow or steep. We find that shallow phase curves are associated with steep slopes on crater walls and faults, as calculated from a shape model. On the other hand, the phase curve of ejecta associated with young impact craters is steep. We interpret these variations in phase curve slope in terms of physical roughness of the regolith. The lack of rough ejecta around older craters suggests that initially rough ejecta associated with impact craters on Vesta are smoothed over a relatively short time of several tens of Myr. We propose that this process is the result of impact gardening, and as such represents a previously unrecognized aspect of Vesta space weathering (Pieters et al., 2012). If this type of space weathering is common, we may expect to encounter this photometric phenomenon on other main belt asteroids.

We report the observations of an electron vortex magnetic hole corresponding to a new type of coherent structures in the magnetosheath turbulent plasma using the Magnetospheric Multiscale (MMS) mission data. The magnetic hole is characterized by a magnetic depression, a density peak, a total electron temperature increase (with a parallel temperature decrease but a perpendicular temperature increase), and strong currents carried by the electrons. The current has a dip in the center of the magnetic hole and a peak in the outer region of the magnetic hole. The estimated size of the magnetic hole is about 0.23 \r{ho}i (~ 30 \r{ho}e) in the circular cross-section perpendicular to its axis, where \r{ho}i and \r{ho}e are respectively the proton and electron gyroradius. There are no clear enhancement seen in high energy electron fluxes, but an enhancement in the perpendicular electron fluxes at ~ 90{\deg} pitch angles inside the magnetic hole is seen, implying that the electron are trapped within it. The variations of the electron velocity components Vem and Ven suggest that an electron vortex is formed by trapping electrons inside the magnetic hole in the circular cross-section (in the M-N plane). These observations demonstrate the existence of a new type of coherent structures behaving as an electron vortex magnetic hole in turbulent space plasmas as predicted by recent kinetic simulations.

The Galactic centre is a hotbed of astrophysical activity, with the injection of wind material from $\sim$30 massive Wolf-Rayet (WR) stars orbiting within 12" of the super-massive black hole (SMBH) playing an important role. Hydrodynamic simulations of such colliding and accreting winds produce a complex density and temperature structure of cold wind material shocking with the ambient medium, creating a large reservoir of hot, X-ray-emitting gas. This work aims to confront the 3Ms of Chandra X-ray Visionary Program (XVP) observations of this diffuse emission by computing the X-ray emission from these hydrodynamic simulations of the colliding WR winds, amid exploring a variety of SMBH feedback mechanisms. The major success of the model is that it reproduces the spectral shape from the 2"--5" ring around the SMBH, where most of the stellar wind material that is ultimately captured by Sgr A* is shock-heated and thermalised. This naturally explains that the hot gas comes from colliding WR winds, and that the wind speeds of these stars are in general well constrained. The flux level of these spectra, as well as 12"$\times$12" images of 4--9 keV, show the X-ray flux is tied to the SMBH feedback strength; stronger feedback clears out more hot gas, thereby decreasing the thermal X-ray emission. The model in which Sgr A* produced an intermediate-strength outflow during the last few centuries best matches the observations to within about 10%, showing SMBH feedback is required to interpret the X-ray emission in this region.

Gamma Cassiopeiae is an enigmatic Be star with unusually strong hard X-ray emission. The Suzaku observatory detected six rapid X-ray spectral hardening events called "softness dips" in a ~100 ksec duration observation in 2011. All the softness dip events show symmetric softness ratio variations, and some of them have flat bottoms apparently due to saturation. The softness dip spectra are best described by either ~40% or ~70% partial covering absorption to kT ~12 keV plasma emission by matter with a neutral hydrogen column density of ~2-8e21 cm-2, while the spectrum outside of these dips is almost free of absorption. This result suggests the presence of two distinct X-ray emitting spots in the gamma Cas system, perhaps on a white dwarf companion with dipole mass accretion. The partial covering absorbers may be blobs in the Be stellar wind, the Be disk, or rotating around the white dwarf companion. Weak correlations of the softness ratios to the hard X-ray flux suggest the presence of stable plasmas at kT ~0.9 and 5 keV, which may originate from the Be or white dwarf winds. The formation of a Be star and white dwarf binary system requires mass transfer between two stars; gamma Cas may have experienced such activity in the past.

We have monitored the massive binary star $\eta$ Carinae with the CTIO/SMARTS 1.5~m telescope and CHIRON spectrograph from the previous apastron passage of the system through the recent 2014.6 periastron passage. Our monitoring has resulted in a large, homogeneous data set with an unprecedented time-sampling, spectral resolving power, and signal-to-noise. This allowed us to investigate temporal variability previously unexplored in the system and discover a kinematic structure in the P Cygni absorption troughs of neutral helium wind lines. The features observed occurred prior to the periastron passage and are seen as we look through the trailing arm of the wind-wind collision shock cone. We show that the bulk of the variability is repeatable across the last five periastron passages, and that the absorption occurs in the inner 230 AU of the system. In addition, we found an additional, high-velocity absorption component super-imposed on the P Cygni absorption troughs that has been previously un-observed in these lines, but which bears resemblance to the observations of the He~I $\lambda$10830 \AA\ feature across previous cycles. Through a comparison of the current smoothed particle hydrodynamical simulations, we show that the observed variations are likely caused by instabilities in the wind-wind collision region in our line-of-sight, coupled with stochastic variability related to clumping in the winds.

The magnetospheric cusps are important sites of the coupling of a magnetosphere with the solar wind. The combination of both ground- and space-based observations at Earth have enabled considerable progress to be made in understanding the terrestrial cusp and its role in the coupling of the magnetosphere to the solar wind via the polar magnetosphere. Voyager 2 fully explored Neptune's cusp in 1989 but highly inclined orbits of the Cassini spacecraft at Saturn present the most recent opportunity to repeatedly studying the polar magnetosphere of a rapidly rotating planet. In this paper we discuss observations made by Cassini during two passes through Saturn's southern polar magnetosphere. Our main findings are that i) Cassini directly encounters the southern polar cusp with evidence for the entry of magnetosheath plasma into the cusp via magnetopause reconnection, ii) magnetopause reconnection and entry of plasma into the cusp can occur over a range of solar wind conditions, and iii) double cusp morphologies are consistent with the position of the cusp oscillating in phase with Saturn's global magnetospheric periodicities.

The X-ray emission of $\eta$ Carinae shows multiple features at various spatial and temporal scales. The central constant emission (CCE) component is centred on the binary and arises from spatial scales much smaller than the bipolar Homunculus nebula, but likely larger than the central wind--wind collision region between the stars as it does not vary over the $\sim$2-3 month X-ray minimum when it can be observed. Using large-scale 3D smoothed particle hydrodynamics (SPH) simulations, we model both the colliding-wind region between the stars, and the region where the secondary wind collides with primary wind ejected from the previous periastron passage. The simulations extend out to one hundred semimajor axes and make two limiting assumptions (strong coupling and no coupling) about the influence of the primary radiation field on the secondary wind. We perform 3D radiative transfer calculations on the SPH output to synthesize the X-ray emission, with the aim of reproducing the CCE spectrum. For the preferred primary mass-loss rate $\dot{M}_A\approx8.5\times10^{-4}$ M$_\odot$ yr$^{-1}$, the model spectra well reproduce the observation as the strong- and no-coupling spectra bound the CCE observation for longitude of periastron $\omega\approx252^\circ$, and bound/converge on the observation for $\omega\approx90^\circ$. This suggests that $\eta$ Carinae has moderate coupling between the primary radiation and secondary wind, that both the region between the stars and the comoving collision on the backside of the secondary generate the CCE, and that the CCE cannot place constraints on the binary's line of sight. We also discuss comparisons with common X-ray fitting parameters.

The evolved, massive highly eccentric binary system, eta Carinae, underwent a periastron passage in the summer of 2014. We obtained two coordinated X-ray observations with XMM-Newton and NuSTAR during the elevated X-ray flux state and just before the X-ray minimum flux state around this passage. These NuSTAR observations clearly detected X-ray emission associated with eta Car extending up to ~50 keV for the first time. The NuSTAR spectrum above 10 keV can be fit with the bremsstrahlung tail from a kT ~6 keV plasma. This temperature is Delta kT ~2 keV higher than those measured from the iron K emission line complex, if the shocked gas is in collisional ionization equilibrium. This result may suggest that the companion star's pre-shock wind velocity is underestimated. The NuSTAR observation near the X-ray minimum state showed a gradual decline in the X-ray emission by 40% at energies above 5 keV in a day, the largest rate of change of the X-ray flux yet observed in individual eta Car observations. The column density to the hardest emission component, NH ~1e24 cm-2, marked one of the highest values ever observed for eta Car, strongly suggesting increased obscuration of the wind-wind colliding X-ray emission by the thick primary stellar wind prior to superior conjunction. Neither observation detected the power-law component in the extremely hard band that INTEGRAL and Suzaku observed prior to 2011. If the non-detection by NuSTAR is caused by absorption, the power-law source must be small and located very near the WWC apex. Alternatively, it may be that the power-law source is not related to either eta Car or the GeV gamma-ray source.

Previous observations suggested that Ceres has active but possibly sporadic water outgassing, and possibly varying spectral characteristics in a time scale of months. We used all available data of Ceres collected in the past three decades from the ground and the Hubble Space Telescope, and the newly acquired images by Dawn Framing Camera to search for spectral and albedo variability on Ceres, in both a global scale and local regions, particularly the bright spots inside Occator crater, over time scales of a few months to decades. Our analysis has placed an upper limit on the possible temporal albedo variation on Ceres. Sporadic water vapor venting, or any possibly ongoing activity on Ceres, is not significant enough to change the albedo or the area of the bright features in Occator crater by >15%, or the global albedo by >3% over various time scales that we searched. Recently reported spectral slope variations can be explained by changing Sun-Ceres-Earth geometry. The active area on Ceres is less than 1 km$^2$, too small to cause global albedo and spectral variations detectable in our data. Impact ejecta due to impacting projectiles of tens of meters in size like those known to cause observable changes to the surface albedo on Asteroid Scheila cannot cause detectable albedo change on Ceres due to its relatively large size and strong gravity. The water vapor activity on Ceres is independent of Ceres' heliocentric distance, rulling out the possibility of comet-like sublimation process as a possible mechanism driving the activity.

{\eta} Carinae is an extremely massive binary system in which rapid spectrum variations occur near periastron. Most notably, near periastron the He II $\lambda 4686$ line increases rapidly in strength, drops to a minimum value, then increases briefly before fading away. To understand this behavior, we conducted an intense spectroscopic monitoring of the He II $\lambda 4686$ emission line across the 2014.6 periastron passage using ground- and space-based telescopes. Comparison with previous data confirmed the overall repeatability of EW(He II $\lambda 4686$), the line radial velocities, and the timing of the minimum, though the strongest peak was systematically larger in 2014 than in 2009 by 26%. The EW(He II $\lambda 4686$) variations, combined with other measurements, yield an orbital period $2022.7\pm0.3$ d. The observed variability of the EW(He II $\lambda 4686$) was reproduced by a model in which the line flux primarily arises at the apex of the wind-wind collision and scales inversely with the square of the stellar separation, if we account for the excess emission as the companion star plunges into the hot inner layers of the primary's atmosphere, and including absorption from the disturbed primary wind between the source and the observer. This model constrains the orbital inclination to $135^\circ$-$153^\circ$, and the longitude of periastron to $234^\circ$-$252^\circ$. It also suggests that periastron passage occurred on $T_0 = 2456874.4\pm1.3$ d. Our model also reproduced EW(He II $\lambda 4686$) variations from a polar view of the primary star as determined from the observed He II $\lambda 4686$ emission scattered off the Homunculus nebula.

Colliding Wolf-Rayet (WR) winds produce thermal X-ray emission widely observed by X-ray telescopes. In wide WR+O binaries, such as WR 140, the X-ray flux is tied to the orbital phase, and is a direct probe of the winds' properties. In the Galactic center, $\sim$30 WRs orbit the super massive black hole (SMBH) within $\sim$10", leading to a smorgasbord of wind-wind collisions. To model the X-ray emission of WR 140 and the Galactic center, we perform 3D hydrodynamic simulations to trace the complex gaseous flows, and then carry out 3D radiative transfer calculations to compute the variable X-ray spectra. The model WR 140 RXTE light curve matches the data well for all phases except the X-ray minimum associated with periastron, while the model spectra agree with the RXTE hardness ratio and the shape of the Suzaku observations throughout the orbit. The Galactic center model of the Chandra flux and spectral shape match well in the region r$<$3", but the model flux falls off too rapidly beyond this radius.

Eta Carinae is the nearest example of a supermassive, superluminous, unstable star. Mass loss from the system is critical in shaping its circumstellar medium and in determining its ultimate fate. Eta Car currently loses mass via a dense, slow stellar wind and possesses one of the largest mass loss rates known. It is prone to episodes of extreme mass ejection via eruptions from some as-yet unspecified cause; the best examples of this are the large-scale eruptions which occurred in 19th century. Eta Car is a colliding wind binary in which strong variations in X-ray emission and in other wavebands are driven by the violent collision of the wind of eta Car-A and the fast, less dense wind of an otherwise hidden companion star. X-ray variations are the simplest diagnostic we have to study the wind-wind collision and allow us to measure the state of the stellar mass loss from both stars. We present the X-ray lightcurve over the last 20 years from ROSAT observations and monitoring with the Rossi X-ray Timing Explorer and the X-ray Telescope on the Swift satellite. We compare and contrast the behavior of the X-ray emission from the system over that timespan, including surprising variations during the 2014 X-ray minimum.

We present an overview of four phase-constrained Chandra HETGS X-ray observations of Delta Ori A. Delta Ori A is actually a triple system which includes the nearest massive eclipsing spectroscopic binary, Delta Ori Aa, the only such object which can be observed with little phase-smearing with the Chandra gratings. Since the fainter star, Delta Ori Aa2, has a much lower X-ray luminosity than the brighter primary, Delta Ori A provides a unique system with which to test the spatial distribution of the X-ray emitting gas around Delta Ori Aa1 via occultation by the photosphere of and wind cavity around the X-ray dark secondary. Here we discuss the X-ray spectrum and X-ray line profiles for the combined observation, having an exposure time of nearly 500 ksec and covering nearly the entire binary orbit. Companion papers discuss the X-ray variability seen in the Chandra spectra, present new space-based photometry and ground-based radial velocities simultaneous with the X-ray data to better constrain the system parameters, and model the effects of X-rays on the optical and UV spectrum. We find that the X-ray emission is dominated by embedded wind shock emission from star Aa1, with little contribution from the tertiary star Ab or the shocked gas produced by the collision of the wind of Aa1 against the surface of Aa2. We find a similar temperature distribution to previous X-ray spectrum analyses. We also show that the line half-widths are about $0.3-0.5\times$ the terminal velocity of the wind of star Aa1. We find a strong anti-correlation between line widths and the line excitation energy, which suggests that longer-wavelength, lower-temperature lines form farther out in the wind. Our analysis also indicates that the ratio of the intensities of the strong and weak lines of \ion{Fe}{17} and \ion{Ne}{10} are inconsistent with model predictions, which may be an effect of resonance scattering

We present time-resolved and phase-resolved variability studies of an extensive X-ray high-resolution spectral dataset of the $\delta$ Orionis Aa binary system. The four observations, obtained with Chandra ACIS HETGS, have a total exposure time of ~479 ks and provide nearly complete binary phase coverage. Variability of the total X-ray flux in the range 5-25 $\AA$ is confirmed, with maximum amplitude of about +/-15% within a single ~125 ks observation. Periods of 4.76d and 2.04d are found in the total X-ray flux, as well as an apparent overall increase in flux level throughout the 9-day observational campaign. Using 40 ks contiguous spectra derived from the original observations, we investigate variability of emission line parameters and ratios. Several emission lines are shown to be variable, including S XV, Si XIII, and Ne IX. For the first time, variations of the X-ray emission line widths as a function of the binary phase are found in a binary system, with the smallest widths at phase=0.0 when the secondary $\delta$ Orionis Aa2 is at inferior conjunction. Using 3D hydrodynamic modeling of the interacting winds, we relate the emission line width variability to the presence of a wind cavity created by a wind-wind collision, which is effectively void of embedded wind shocks and is carved out of the X-ray-producing primary wind, thus producing phase-locked X-ray variability.

We report on both high-precision photometry from the MOST space telescope and ground-based spectroscopy of the triple system delta Ori A consisting of a binary O9.5II+early-B (Aa1 and Aa2) with P = 5.7d, and a more distant tertiary (O9 IV P > 400 yrs). This data was collected in concert with X-ray spectroscopy from the Chandra X-ray Observatory. Thanks to continuous coverage for 3 weeks, the MOST light curve reveals clear eclipses between Aa1 and Aa2 for the first time in non-phased data. From the spectroscopy we have a well constrained radial velocity curve of Aa1. While we are unable to recover radial velocity variations of the secondary star, we are able to constrain several fundamental parameters of this system and determine an approximate mass of the primary using apsidal motion. We also detected second order modulations at 12 separate frequencies with spacings indicative of tidally influenced oscillations. These spacings have never been seen in a massive binary, making this system one of only a handful of such binaries which show evidence for tidally induced pulsations.

The Suzaku X-ray observatory monitored the supermassive binary system Eta Carinae 10 times during the whole 5.5 year orbital cycle between 2005-2011. This series of observations presents the first long-term monitoring of this enigmatic system in the extremely hard X-ray band between 15-40 keV. During most of the orbit, the 15-25 keV emission varied similarly to the 2-10 keV emission, indicating an origin in the hard energy tail of the kT ~4 keV wind-wind collision (WWC) plasma. However, the 15-25 keV emission declined only by a factor of 3 around periastron when the 2-10 keV emission dropped by two orders of magnitude due probably to an eclipse of the WWC plasma. The observed minimum in the 15-25 keV emission occurred after the 2-10 keV flux had already recovered by a factor of ~3. This may mean that the WWC activity was strong, but hidden behind the thick primary stellar wind during the eclipse. The 25-40 keV flux was rather constant through the orbital cycle, at the level measured with INTEGRAL in 2004. This result may suggest a connection of this flux component to the gamma-ray source detected in this field. The Helium-like Fe Kalpha line complex at ~6.7 keV became strongly distorted toward periastron as seen in the previous cycle. The 5-9 keV spectra can be reproduced well with a two-component spectral model, which includes plasma in collision equilibrium (CE) and a plasma in non-equilibrium ionization (NEI) with tau ~1e11 cm-3 s-1. The NEI plasma increases in importance toward periastron.

We present results from a study of the eclipsing, colliding-wind binary V444 Cyg that uses a combination of X-ray and optical spectropolarimetric methods to describe the 3-D nature of the shock and wind structure within the system. We have created the most complete X-ray light curve of V444 Cyg to date using 40 ksec of new data from Swift, and 200 ksec of new and archived XMM-Newton observations. In addition, we have characterized the intrinsic, polarimetric phase-dependent behavior of the strongest optical emission lines using data obtained with the University of Wisconsin's Half-Wave Spectropolarimeter. We have detected evidence of the Coriolis distortion of the wind-wind collision in the X-ray regime, which manifests itself through asymmetric behavior around the eclipses in the system's X-ray light curves. The large opening angle of the X-ray emitting region, as well as its location (i.e. the WN wind does not collide with the O star, but rather its wind) are evidence of radiative braking/inhibition occurring within the system. Additionally, the polarimetric results show evidence of the cavity the wind-wind collision region carves out of the Wolf-Rayet star's wind.

We present an analysis of the olivine-rich exposures at Bellicia and Arruntia craters using Dawn Framing Camera (FC) color data. Our results confirm the existence of olivine-rich materials at these localities as described by Ammannito et al. (2013a) using Visual Infrared Spectrometer (VIR) data. Analyzing laboratory spectra of various Howardite-Eucrite-Diogenite meteorites, high-Ca pyroxenes, olivines and olivine-orthopyroxene mixtures, we derive three FC spectral band parameters that are indicators of olivine-rich materials. Combining the three band parameters allows us, for the first time, to reliably identify sites showing modal olivine contents >40%. The olivine-rich exposures at Bellicia and Arruntia are mapped using higher spatial resolution FC data. The exposures are located on the slopes of outer/inner crater walls, on the floor of Arruntia, in the ejecta, as well as in nearby fresh small impact craters. The spatial extent of the exposures ranges from a few hundred meters to few kilometers. The olivine-rich exposures are in accordance with both the magma ocean and the serial magmatism model (e.g., Righter and Drake 1997; Yamaguchi et al. 1997). However, it remains unsolved why the olivine-rich materials are mainly concentrated in the northern hemisphere (~36-42{\deg} N, 46-74{\deg} E) and are almost absent in the Rheasilvia basin.

Vesta has a complex cratering history, with ancient terrains as well as recent large impacts that have led to regional resurfacing. Crater counts can help constrain the relative ages of different units on Vesta's surface, but converting those crater counts to absolute ages requires a chronology function. We present a cratering chronology based on the best current models for the dynamical evolution of asteroid belt, and calibrate it to Vesta using the record of large craters on its surface. While uncertainties remain, our chronology function is broadly consistent with an ancient surface of Vesta as well as other constraints such as the bombardment history of the rest of the inner Solar System and the Ar-Ar age distribution of howardite, eucrite and diogenite (HED) meteorites from Vesta.

The Dawn spacecraft observed the presence of dark material, which in turn proved to be associated with OH and H-rich material, on the surface of Vesta. The source of this dark material has been identified with the low albedo asteroids, but it is still a matter of debate whether the delivery of the dark material is associated with a few large impact events, to micrometeorites or to the continuous, secular flux of impactors on Vesta. The continuous flux scenario predicts that a significant fraction of the exogenous material accreted by Vesta should be due to non-dark impactors likely analogous to ordinary chondrites, which instead represent only a minor contaminant in the HED meteorites. We explored the continuous flux scenario and its implications for the composition of the vestan regolith, taking advantage of the data from the Dawn mission and the HED meteorites. We used our model to show that the stochastic events scenario and the micrometeoritic flux scenario are natural consequences of the continuous flux scenario. We then used the model to estimate the amounts of dark and hydroxylate materials delivered on Vesta since the LHB and we showed how our results match well with the values estimated by the Dawn mission. We used our model to assess the amount of Fe and siderophile elements that the continuous flux of impactors would mix in the vestan regolith: concerning the siderophile elements, we focused our attention on the role of Ni. The results are in agreement with the data available on the Fe and Ni content of the HED meteorites and can be used as a reference frame in future studies of the data from the Dawn mission and of the HED meteorites. Our model cannot yet provide an answer to the fate of the missing non-carbonaceous contaminants, but we discuss possible reasons for this discrepancy.

Cyg X-3 is a high mass X-ray binary consisting of a Wolf-Rayet star and a compact object in a very short orbital period of 4.8h. The only confirmed microquasar with high energy gamma-ray emission, Cyg X-3 provides a unique opportunity to study the relationship between the accretion power and the power in high energy emission. Because of a compact orbit and a slow Wolf-Rayet wind, the flow structure around the compact object is thought to be strongly affected by the orbital motion, details of which can be obtained only by numerical simulations. In this paper, we report on the results from 3D hydrodynamic simulations of the wind accretion in Cyg X-3. For simplicity we adopt an anti-gravity-like force that emulates the radiative acceleration consistent with the beta-velocity wind. Due to the rapid orbital motion, the flow around the compact object has large density gradients. As a result, the accretion rate onto the compact object is significantly lower than that of the Bondi-Hoyle-Lyttleton rate. We also calculate the model X-ray light curve. Although it roughly agrees with the observed light curve, more detailed modeling is needed for detailed comparison.

Colliding wind binaries (CWBs) are unique laboratories for X-ray astrophysics. The massive stars in these systems possess powerful stellar winds with speeds up to $\sim$3000 km s$^{-1}$, and their collision leads to hot plasma (up to $\sim10^8$K) that emit thermal X-rays (up to $\sim$10 keV). Many X-ray telescopes have observed CWBs, including Suzaku, and our work aims to model these X-ray observations. We use 3D smoothed particle hydrodynamics (SPH) to model the wind-wind interaction, and then perform 3D radiative transfer to compute the emergent X-ray flux, which is folded through X-ray telescopes' response functions to compare directly with observations. In these proceedings, we present our models of Suzaku observations of the multi-year-period, highly eccentric systems $\eta$ Carinae and WR 140. The models reproduce the observations well away from periastron passage, but only $\eta$ Carinae's X-ray spectrum is reproduced at periastron; the WR 140 model produces too much flux during this more complicated phase.

X-ray emission from the supermassive binary system Eta Carinae declines sharply around periastron. This X-ray minimum has two distinct phases - the lowest flux phase in the first ~3 weeks and a brighter phase thereafter. In 2009, the Chandra X-ray Observatory monitored the first phase five times and found the lowest observed flux at ~1.9e-12 ergs cm-2 s-1 (3-8 keV). The spectral shape changed such that the hard band above ~4 keV dropped quickly at the beginning and the soft band flux gradually decreased to its lowest observed value in ~2 weeks. The hard band spectrum had begun to recover by that time. This spectral variation suggests that the shocked gas producing the hottest X-ray gas near the apex of the wind-wind collision (WWC) is blocked behind the dense inner wind of the primary star, which later occults slightly cooler gas downstream. Shocked gas previously produced by the system at earlier orbital phases is suggested to produce the faint residual X-ray emission seen when the emission near the apex is completely blocked by the primary wind. The brighter phase is probably caused by the re-appearance of the WWC plasma, whose emissivity significantly declined during the occultation. We interpret this to mean that the X-ray minimum is produced by a hybrid mechanism of an occultation and a decline in emissivity of the WWC shock.

Recent work suggests that the mass-loss rate of the primary star (Eta A) in the massive colliding wind binary Eta Carinae dropped by a factor of 2-3 between 1999 and 2010. We present results from large- (r=1545au) and small- (r=155au) domain, 3D smoothed particle hydrodynamic (SPH) simulations of Eta Car's colliding winds for 3 Eta A mass-loss rates (2.4, 4.8, and 8.5 x 10^-4 M_sun/yr), investigating the effects on the dynamics of the binary wind-wind collision (WWC). These simulations include orbital motion, optically thin radiative cooling, and radiative forces. We find that Eta A's mass-loss rate greatly affects the time-dependent hydrodynamics at all spatial scales investigated. The simulations also show that the post-shock wind of the companion star (Eta B) switches from the adiabatic to the radiative-cooling regime during periastron passage. The SPH simulations together with 1D radiative transfer models of Eta A's spectra reveal that a factor of 2 or more drop in Eta A's mass-loss rate should lead to substantial changes in numerous multiwavelength observables. Recent observations are not fully consistent with the model predictions, indicating that any drop in Eta A's mass-loss rate was likely by a factor < 2 and occurred after 2004. We speculate that most of the recent observed changes in Eta Car are due to a small increase in the WWC opening angle that produces significant effects because our line-of-sight to the system lies close to the dense walls of the WWC zone. A modest decrease in Eta A's mass-loss rate may be responsible, but changes in the wind/stellar parameters of Eta B cannot yet be fully ruled out. We suggest observations during Eta Car's next periastron in 2014 to further test for decreases in Eta A's mass-loss rate. If Eta A's mass-loss rate is declining and continues to do so, the 2014 X-ray minimum should be even shorter than that of 2009.

Dawn spacecraft orbited Vesta for more than one year and collected a huge volume of multispectral, high-resolution data in the visible wavelengths with the Framing Camera. We present a detailed disk-integrated and disk-resolved photometric analysis using the Framing Camera images with the Minnaert model and the Hapke model, and report our results about the global photometric properties of Vesta. The photometric properties of Vesta show weak or no dependence on wavelengths, except for the albedo. At 554 nm, the global average geometric albedo of Vesta is 0.38+/-0.04, and the Bond albedo range is 0.20+/-0.02. The bolometric Bond albedo is 0.18+/-0.01. The phase function of Vesta is similar to those of S-type asteroids. Vesta's surface shows a single-peaked albedo distribution with a full-width-half-max ~17% relative to the global average. This width is much smaller than the full range of albedos (from ~0.55x to >2x global average) in localized bright and dark areas of a few tens of km in sizes, and is probably a consequence of significant regolith mixing on the global scale. Rheasilvia basin is about 10% brighter than the global average. The phase reddening of Vesta measured from Dawn Framing Camera images is comparable or slightly stronger than that of Eros as measured by the Near Earth Asteroid Rendezvous mission, but weaker than previous measurements based on ground-based observations of Vesta and laboratory measurements of HED meteorites. The photometric behaviors of Vesta are best described by the Hapke model and the Akimov disk- function, when compared with the Minnaert model, Lommel-Seeliger model, and Lommel- Seeliger-Lambertian model. The traditional approach for photometric correction is validated for Vesta for >99% of its surface where reflectance is within +/-30% of global average.

NASA's Dawn mission observed a great variety of colored terrains on asteroid (4) Vesta during its survey with the Framing Camera (FC). Here we present a detailed study of the orange material on Vesta, which was first observed in color ratio images obtained by the FC and presents a red spectral slope. The orange material deposits can be classified into three types, a) diffuse ejecta deposited by recent medium-size impact craters (such as Oppia), b) lobate patches with well-defined edges, and c) ejecta rays from fresh-looking impact craters. The location of the orange diffuse ejecta from Oppia corresponds to the olivine spot nicknamed "Leslie feature" first identified by Gaffey (1997) from ground-based spectral observations. The distribution of the orange material in the FC mosaic is concentrated on the equatorial region and almost exclusively outside the Rheasilvia basin. Our in-depth analysis of the composition of this material uses complementary observations from FC, the visible and infrared spectrometer (VIR), and the Gamma Ray and Neutron Detector (GRaND). Combining the interpretations from the topography, geomorphology, color and spectral parameters, and elemental abundances, the most probable analog for the orange material on Vesta is impact melt.

Observations of asteroid 4 Vesta by NASA's Dawn spacecraft are interesting because its surface has the largest range of albedo, color and composition of any other asteroid visited by spacecraft to date. These hemispherical and rotational variations in surface brightness and composition have been attributed to impact processes since Vesta's formation. Prior to Dawn's arrival at Vesta, its surface properties were the focus of intense telescopic investigations for nearly a hundred years. Ground-based photometric and spectroscopic observations first revealed these variations followed later by those using Hubble Space Telescope. Here we compare interpretations of Vesta's rotation period, pole, albedo, topographic, color, and compositional properties from ground-based telescopes and HST with those from Dawn. Rotational spectral variations observed from ground-based studies are also consistent with those observed by Dawn. While the interpretation of some of these features was tenuous from past data, the interpretations were reasonable given the limitations set by spatial resolution and our knowledge of Vesta and HED meteorites at that time. Our analysis shows that ground-based and HST observations are critical for our understanding of small bodies and provide valuable support for ongoing and future spacecraft missions.

The NASA Dawn mission has extensively examined the surface of asteroid Vesta, the second most massive body in the main belt. The high quality of the gathered data provides us with an unique opportunity to determine the surface and internal properties of one of the most important and intriguing main belt asteroids (MBAs). In this paper, we focus on the size frequency distributions (SFDs) of sub-kilometer impact craters observed at high spatial resolution on several selected young terrains on Vesta. These small crater populations offer an excellent opportunity to determine the nature of their asteroidal precursors (namely MBAs) at sizes that are not directly observable from ground-based telescopes (i.e., below ~100 m diameter). Moreover, unlike many other MBA surfaces observed by spacecraft thus far, the young terrains examined had crater spatial densities that were far from empirical saturation. Overall, we find that the cumulative power-law index (slope) of small crater SFDs on Vesta is fairly consistent with predictions derived from current collisional and dynamical models down to a projectile size of ~10 m diameter (Bottke et al., 2005a,b). The shape of the impactor SFD for small projectile sizes does not appear to have changed over the last several billions of years, and an argument can be made that the absolute number of small MBAs has remained roughly constant (within a factor of 2) over the same time period. The apparent steady state nature of the main belt population potentially provides us with a set of intriguing constraints that can be used to glean insights into the physical evolution of individual MBAs as well as the main belt as an ensemble.

The surface composition of Vesta, the most massive intact basaltic object in the asteroid belt, is interesting because it provides us with an insight into magmatic differentiation of planetesimals that eventually coalesced to form the terrestrial planets. The distribution of lithologic and compositional units on the surface of Vesta provides important constraints on its petrologic evolution, impact history and its relationship with Vestoids and howardite-eucrite-diogenite (HED) meteorites. Using color parameters (band tilt and band curvature) originally developed for analyzing lunar data, we have identified and mapped HED terrains on Vesta in Dawn Framing Camera (FC) color data. The average color spectrum of Vesta is identical to that of howardite regions, suggesting an extensive mixing of surface regolith due to impact gardening over the course of solar system history. Our results confirm the hemispherical dichotomy (east-west and north-south) in albedo/color/composition that has been observed by earlier studies. The presence of diogenite-rich material in the southern hemisphere suggests that it was excavated during the formation of the Rheasilvia and Veneneia basins. Our lithologic mapping of HED regions provides direct evidence for magmatic evolution of Vesta with diogenite units in Rheasilvia forming the lower crust of a differentiated object.

NASA's Dawn spacecraft observations of asteroid (4) Vesta reveal a surface with the highest albedo and color variation of any asteroid we have observed so far. Terrains rich in low albedo dark material (DM) have been identified using Dawn Framing Camera (FC) 0.75 {\mu}m filter images in several geologic settings: associated with impact craters (in the ejecta blanket material and/or on the crater walls and rims); as flow-like deposits or rays commonly associated with topographic highs; and as dark spots (likely secondary impacts) nearby impact craters. This DM could be a relic of ancient volcanic activity or exogenic in origin. We report that the majority of the spectra of DM are similar to carbonaceous chondrite meteorites mixed with materials indigenous to Vesta. Using high-resolution seven color images we compared DM color properties (albedo, band depth) with laboratory measurements of possible analog materials. Band depth and albedo of DM are identical to those of carbonaceous chondrite xenolith-rich howardite Mt. Pratt (PRA) 04401. Laboratory mixtures of Murchison CM2 carbonaceous chondrite and basaltic eucrite Millbillillie also show band depth and albedo affinity to DM. Modeling of carbonaceous chondrite abundance in DM (1-6 vol%) is consistent with howardite meteorites. We find no evidence for large-scale volcanism (exposed dikes/pyroclastic falls) as the source of DM. Our modeling efforts using impact crater scaling laws and numerical models of ejecta reaccretion suggest the delivery and emplacement of this DM on Vesta during the formation of the ~400 km Veneneia basin by a low-velocity (<2 km/sec) carbonaceous impactor. This discovery is important because it strengthens the long-held idea that primitive bodies are the source of carbon and probably volatiles in the early Solar System.

Eta Carinae is considered to be a massive colliding wind binary system with a highly eccentric (e \sim 0.9), 5.54-yr orbit. However, the companion star continues to evade direct detection as the primary dwarfs its emission at most wavelengths. Using three-dimensional (3-D) SPH simulations of Eta Car's colliding winds and radiative transfer codes, we are able to compute synthetic observables across multiple wavebands for comparison to the observations. The models show that the presence of a companion star has a profound influence on the observed HST/STIS UV spectrum and H-alpha line profiles, as well as the ground-based photometric monitoring. Here, we focus on the Bore Hole effect, wherein the fast wind from the hot secondary star carves a cavity in the dense primary wind, allowing increased escape of radiation from the hotter/deeper layers of the primary's extended wind photosphere. The results have important implications for interpretations of Eta Car's observables at multiple wavelengths.

We present a three-dimensional (3-D) dynamical model for the broad [Fe III] emission observed in Eta Carinae using the Hubble Space Telescope/Space Telescope Imaging Spectrograph (HST/STIS). This model is based on full 3-D Smoothed Particle Hydrodynamics (SPH) simulations of Eta Car's binary colliding winds. Radiative transfer codes are used to generate synthetic spectro-images of [Fe III] emission line structures at various observed orbital phases and STIS slit position angles (PAs). Through a parameter study that varies the orbital inclination i, the PA {\theta} that the orbital plane projection of the line-of-sight makes with the apastron side of the semi-major axis, and the PA on the sky of the orbital axis, we are able, for the first time, to tightly constrain the absolute 3-D orientation of the binary orbit. To simultaneously reproduce the blue-shifted emission arcs observed at orbital phase 0.976, STIS slit PA = +38 degrees, and the temporal variations in emission seen at negative slit PAs, the binary needs to have an i \approx 130 to 145 degrees, {\theta} \approx -15 to +30 degrees, and an orbital axis projected on the sky at a PA \approx 302 to 327 degrees east of north. This represents a system with an orbital axis that is closely aligned with the inferred polar axis of the Homunculus nebula, in 3-D. The companion star, Eta B, thus orbits clockwise on the sky and is on the observer's side of the system at apastron. This orientation has important implications for theories for the formation of the Homunculus and helps lay the groundwork for orbital modeling to determine the stellar masses.

The colliding wind binary (CWB) systems \eta\ Carinae and WR140 provide unique laboratories for X-ray astrophysics. Their wind-wind collisions produce hard X-rays that have been monitored extensively by several X-ray telescopes, including RXTE. To interpret these X-ray light curves and spectra, we apply 3D hydrodynamic simulations of the wind-wind collision using smoothed particle hydrodynamics (SPH), with the recent improvements of radiative cooling and the acceleration of the stellar winds according to a \beta\ law. For both systems, the 2-10 keV RXTE light curves are well-reproduced in absolute units for most phases, but the light curve dips associated with the periastron passages are not well matched. In WR140, the dip is too weak, and in \eta\ Carinae, the large difference in wind speeds of the two stars leads to a hot, post-periastron bubble that produces excess emission toward the end of the X-ray minimum.

We report a comprehensive review of the UV-visible spectrum and rotational lightcurve of Vesta combining new observations by Hubble Space Telescope and Swift Gamma-ray Burst Observatory with archival International Ultraviolet Explorer observations. The geometric albedos of Vesta from 220 nm to 953 nm are derived by carefully comparing these observations from various instruments at different times and observing geometries. Vesta has a rotationally averaged geometric albedo of 0.09 at 250 nm, 0.14 at 300 nm, 0.26 at 373 nm, 0.38 at 673 nm, and 0.30 at 950 nm. The linear spectral slope as measured between 240 and 320 nm in the ultraviolet displays a sharp minimum near a sub-Earth longitude of 20^{\circ}, and maximum in the eastern hemisphere. This is consistent with the longitudinal distribution of the spectral slope in the visible wavelength. The photometric uncertainty in the ultraviolet is ~20%, and in the visible wavelengths it is better than 10%. The amplitude of Vesta's rotational lightcurves is ~10% throughout the range of wavelengths we observed, but is smaller at 950 nm (~6%) near the 1-\mum band center. Contrary to earlier reports, we found no evidence for any difference between the phasing of the ultraviolet and visible/near-infrared lightcurves with respect to sub-Earth longitude. Vesta's average spectrum between 220 and 950 nm can well be described by measured reflectance spectra of fine particle howardite-like materials of basaltic achondrite meteorites. Combining this with the in-phase behavior of the ultraviolet, visible, and near-infrared lightcurves, and the spectral slopes with respect to the rotational phase, we conclude that there is no global ultraviolet/visible reversal on Vesta. Consequently, this implies a lack of global space weathering on Vesta, as previously inferred from visible-near-infrared data.

The colliding wind binary (CWB) systems \eta\ Carinae and WR140 provide unique laboratories for X-ray astrophysics. Their wind-wind collisions produce hard X-rays that have been monitored extensively by several X-ray telescopes, including RXTE. To interpret these RXTE X-ray light curves, we model the wind-wind collision using 3D smoothed particle hydrodynamics (SPH) simulations. Adiabatic simulations that account for the absorption of X-rays from an assumed point source at the apex of the wind-collision shock cone by the distorted winds can closely match the observed 2-10keV RXTE light curves of both \eta\ Car and WR140. This point-source model can also explain the early recovery of \eta\ Car's X-ray light curve from the 2009.0 minimum by a factor of 2-4 reduction in the mass loss rate of \eta\ Car. Our more recent models relax the point-source approximation and account for the spatially extended emission along the wind-wind interaction shock front. For WR140, the computed X-ray light curve again matches the RXTE observations quite well. But for \eta\ Car, a hot, post-periastron bubble leads to an emission level that does not match the extended X-ray minimum observed by RXTE. Initial results from incorporating radiative cooling and radiatively-driven wind acceleration via a new anti-gravity approach into the SPH code are also discussed.

We conducted a satellite search around the dwarf planet 1 Ceres using Hubble Space Telescope and ground-based Palomar data. No candidate objects were found orbiting Ceres in its entire stability region down to ~500km from the surface of Ceres. Assuming a satellite would have the same albedo as Ceres, which has a visual geometric albedo of 0.07-0.10, our detection limit is sensitive to satellites larger than 1-2 km in diameter.

WR 140 is a canonical massive "colliding wind" binary system in which periodically-varying X-ray emission is produced by the collision between the wind of the WC7 and O4-5 star components in the space between the two stars. We have obtained X-ray observations using the RXTE satellite observatory through almost one complete orbital cycle including two consecutive periastron passages. We discuss the results of this observing campaign, and the implications of the X-ray data for our understanding of the orbital dynamics and the stellar mass loss.

This study examines the tail disconnection event on April 20, 2007 on comet 2P/Encke, caused by a coronal mass ejection (CME) at a heliocentric distance of 0.34 AU. During their interaction, both the CME and the comet are visible with high temporal and spatial resolution by the STEREO-A spacecraft. Previously, only current sheets or shocks have been accepted as possible reasons for comet tail disconnections, so it is puzzling that the CME caused this event. The MHD simulation presented in this work reproduces the interaction process and demonstrates how the CME triggered a tail disconnection in the April 20 event. It is found that the CME disturbs the comet with a combination of a $180^\circ$ sudden rotation of the interplanetary magnetic field (IMF), followed by a $90^\circ$ gradual rotation. Such an interpretation applies our understanding of solar wind-comet interactions to determine the \textit{in situ} IMF orientation of the CME encountering Encke.

The very massive star system $\eta$ Carinae exhibits regular 5.54-year (2024-day) period disruptive events in wavebands ranging from the radio to X-ray. There is a growing consensus that these events likely stem from periastron passage of an (as yet) unseen companion in a highly eccentric ($\epsilon \sim 0.9$) orbit. This paper presents three-dimensional (3-D) Smoothed Particle Hydrodynamics (SPH) simulations of the orbital variation of the binary wind-wind collision, and applies these to modeling the X-ray light curve observed by the Rossi X-ray Timing Explorer (RXTE). By providing a global 3-D model of the phase variation of the density of the interacting winds, the simulations allow computation of the associated variation in X-ray absorption, presumed here to originate from near the apex of the wind-wind interaction cone. We find that the observed RXTE light curve can be readily fit if the observer's line of sight is within this cone along the general direction of apastron. Specifically, the data are well fit by an assumed inclination $i = 45^{\circ}$ for the orbit's polar axis, which is thus consistent with orbital angular momentum being along the inferred polar axis of the Homunculus nebula. The fits also constrain the position angle $\phi$ that an orbital-plane projection makes with the apastron side of the semi-major axis, strongly excluding positions $\phi < 9^{\circ}$ along or to the retrograde side of the axis, with the best fit position given by $\phi = 27^{\circ}$. Overall the results demonstrate the utility of a fully 3-D dynamical model for constraining the geometric and physical properties of this complex colliding-wind binary system.

We study colliding winds in the superluminous binary eta Carinae by performing three-dimensional, Smoothed Particle Hydrodynamics (SPH) simulations. For simplicity, we assume both winds to be isothermal. We also assume that wind particles coast without any net external forces. We find that the lower density, faster wind from the secondary carves out a spiral cavity in the higher density, slower wind from the primary. Because of the phase-dependent orbital motion, the cavity is very thin on the periastron side, whereas it occupies a large volume on the apastron side. The model X-ray light curve using the simulated density structure fits very well with the observed light curve for a viewing angle of i=54 degrees and phi=36 degrees, where i is the inclination angle and phi is the azimuth from apastron.