We present Cassini VIMS observations of sun glitter -- wave-induced reflections from a liquid surface offset from a specular point -- on Kraken Mare. Sun glitter reveals rough sea surfaces around Kraken Mare, namely the coasts and narrow straits. The sun glitter observations indicate wave activity driven by the winds and tidal currents in Kraken Mare during northern summer. T104 Cassini VIMS observations show three sun glitter features in Bayta Fretum indicative of variegated wave fields. We cannot uniquely determine one source for the coastal Bayta waves, but we lean toward the interpretation of surface winds, because tidal currents should be too weak to generate capillary-gravity waves in Bayta Fretum. T105 and T110 observations reveal wave fields in the straits of Seldon Fretum, Lulworth Sinus, and Tunu Sinus that likely originate from the constriction of tidal currents. Coastlines of Bermoothes and Hufaidh Insulae adjoin rough sea surfaces, suggesting a complex interplay of wind-roughened seas and localized tidal currents. Bermoothes and Hufaidh Insulae may share characteristics of either the Torres Strait off Australia or the Aland region of Finland, summarized as an island-dense strait with shallow bathymetry that hosts complex surface circulation patterns. Hufaidh Insulae could host seafloor bedforms formed by tidal currents with an abundant sediment supply, similar to the Torres Strait. The coastlines of Hufaidh and Bermoothes Insulae likely host ria or flooded coastal inlets, suggesting the Insulae may be local peaks of primordial crust isolated by an episode of sea-level rise or tectonic uplift.

The nearest exoplanets to the Sun are our best possibilities for detailed characterization. We report the discovery of a compact multi-planet system of super-Earths orbiting the nearby red dwarf GJ 887, using radial velocity measurements. The planets have orbital periods of 9.3 and 21.8~days. Assuming an Earth-like albedo, the equilibrium temperature of the 21.8 day planet is approx 350 K; which is interior, but close to the inner edge, of the liquid-water habitable zone. We also detect a further unconfirmed signal with a period of 50 days which could correspond to a third super-Earth in a more temperate orbit. GJ 887 is an unusually magnetically quiet red dwarf with a photometric variability below 500 parts-per-million, making its planets amenable to phase-resolved photometric characterization.

Short gamma-ray bursts are thought to result from the mergers of two neutron stars or a neutron star and stellar mass black hole. The final stages of the merger are generally accompanied by the production of one or more tidal "tails" of ejecta, which fall back onto the remnant-disc system at late times. Using the results of a linear stability analysis, we show that if the material comprising these tails is modeled as adiabatic and the effective adiabatic index satisfies $\gamma \ge 5/3$, then the tails are gravitationally unstable and collapse to form small-scale knots. We analytically estimate the properties of these knots, including their spacing along the tidal tail and the total number produced, and their effect on the mass return rate to the merger remnant. We perform hydrodynamical simulations of the disruption of a polytropic (with the polytropic and adiabatic indices $\gamma$ equal), $\gamma =2$ neutron star by a black hole, and find agreement between the predictions of the linear stability analysis and the distribution of knots that collapse out of the instability. The return of these knots to the black hole induces variability in the fallback rate, which can manifest as variability in the lightcurve of the GRB and -- depending on how rapidly the instability operates -- the prompt emission. The late-time variability induced by the return of these knots is also consistent with the extended emission observed in some GRBs.

We analyse TESS photometry of DMPP-1 (HD 38677; TIC 66560666), a nearby F8V star hosting hot super-Earth planets and a warm Neptune. Using the Transit Least Squares algorithm and other methods we find a transit signal at $P=3.2854^{+0.0032}_{-0.0025}$ d with depth 87$^{+25}_{-30}$ ppm and false alarm probability 1.6%. This is shallower than hitherto published TESS discoveries. The 3.285 d signal is recovered for several, but not all, methods for detrending stellar astrophysical variability. Further observations are needed to improve the significance of the detection. If this transit were due to an Earth-like rocky planet it would have been detected in the RV data, but it is not. The TESS data cover seven individual transits, one of which is consistent with zero depth. The insolation of the putative planet is $990 {\rm S_{\oplus}}$, typical of fluxes experienced by the three known catastrophically disintegrating exoplanets (CDEs). The transits can be self-consistently attributed to a CDE with a mass below the RV detection threshold. We searched for transits of the known RV planets, finding null results and detection thresholds of $< 100$ ppm, which we quantify for each. The DMPP-1 planetary system was discovered as a consequence of circumstellar gas attributed to ablation of hot planets. The RV planets may have been ablated to near-pure iron cores. We place limits on the orbital inclinations of the RV planets where the expected transit depth exceeds the detection threshold. If the 3.2854 day transit detection is confirmed, e.g. with CHEOPS photometry, DMPP-1 would be a first-rate target for JWST spectroscopy.

We present new Doppler images of both components of the double-lined binary $\sigma^{2}$ CrB, based on the high-resolution spectroscopic data collected during 11 nights in 2015 March--April. The observed spectra form two independent data sets with sufficient phase coverage. We apply the least-squares deconvolution to all observed spectra to obtain high signal-to-noise mean profiles, from which we derive the Doppler images of both components of $\sigma^{2}$ CrB simultaneously. The surfaces of both F9 and G0 components are dominated by pronounced polar spots. The F9 component exhibits a weak spot at latitude 30$^{\circ}$ and its mid-to-low latitudes are relatively featureless. The G0 star shows an extended spot structure at latitude 30$^{\circ}$, and its surface spot coverage is larger than that of the F9 star, which suggests a higher level of magnetic activity. With the cross-correlation method, we derive a solar-like surface differential rotation on the G0 star of $\sigma^{2}$ CrB for the first time, and the surface shear rate is $\Delta \Omega = 0.180 \pm 0.004$ rad d$^{-1}$ and $\alpha = \Delta \Omega / \Omega_{eq} = 0.032 \pm 0.001$. We do not obtain a clear surface shear law for the F9 star due to the lack of mid-to-low latitude features, but detect a systematic longitude shift of high-latitude spots, which indicates a slower rotation with respect to the co-rotating frame.

KELT-9 b is an ultra hot Jupiter transiting a rapidly rotating, oblate early-A-type star in a polar orbit. We model the effect of rapid stellar rotation on KELT-9 b's transit light curve using photometry from the Transiting Exoplanet Survey Satellite (\tess) to constrain the planet's true spin-orbit angle and to explore how KELT-9 b may be influenced by stellar gravity darkening. We constrain the host star's equatorial radius to be $1.089\pm0.017$ times as large as its polar radius and its local surface brightness to vary by $\sim38$\% between its hot poles and cooler equator. We model the stellar oblateness and surface brightness gradient and find that it causes the transit light curve to lack the usual symmetry around the time of minimum light. We take advantage of the light curve asymmetry to constrain KELT-9 b's true spin orbit angle (${87^\circ}^{+10^\circ}_{-11^\circ}$), agreeing with \citet{gaudi2017giant} that KELT-9 b is in a nearly polar orbit. We also apply a gravity darkening correction to the spectral energy distribution model from \citet{gaudi2017giant} and find that accounting for rapid rotation gives a better fit to available spectroscopy and yields a more reliable estimate for the star's polar effective temperature.

Previous studies have shown that planets that rotate retrograde (backwards with respect to their orbital motion) generally experience less severe obliquity variations than those that rotate prograde (the same direction as their orbital motion). Here we examine retrograde-rotating planets on eccentric orbits and find a previously unknown secular spin-orbit resonance that can drive significant obliquity variations. This resonance occurs when the frequency of the planet's rotation axis precession becomes commensurate with an orbital eigenfrequency of the planetary system. The planet's eccentricity enables a participating orbital frequency through an interaction in which the apsidal precession of the planet's orbit causes a cyclic nutation of the planet's orbital angular momentum vector. The resulting orbital frequency follows the relationship $f = 2 \dot{\varpi} - \dot{\Omega}$, where $\dot{\varpi}$ and $\dot{\Omega}$ are the rates of the planet's changing longitude of periapsis and ascending node, respectively. We test this mechanism by simulating cases of a simple Earth-Jupiter system, and confirm the predicted resonance. Over the course of 100 Myr, the test Earths with rotation axis precession rates near the predicted resonant frequency experienced pronounced obliquity variations of order $10^\circ$-$30^\circ$. These variations can be significant, and suggest that while retrograde rotation is a stabilizing influence most of the time, retrograde rotators can experience large obliquity variations if they are on eccentric orbits and enter this spin-orbit resonance.

The Dispersed Matter Planet Project targets stars with anomalously low Ca II H&K chromospheric emission. High precision, high cadence radial velocity measurements of the F8V star HD 38677 / DMPP-1 reveal four short period planets. DMPP-1 has log(R'HK) = -5.16 which probably indicates the presence of circumstellar absorbing gas arising from an ablating hot planet. The planets have P$_{\rm orb}$ ~ 2.9 - 19 d, i.e., a compact planetary system with super-Earth (~3 M$_\oplus$) to Neptune-mass (~24 M$_\oplus$) planets. These irradiated planets may be chthonian: remnant cores of giant planets after mass-loss while crossing the Neptune desert. Modelling the possible long-term activity indicators while searching for Keplerian signals modifies the recovered planetary signals. A priori inferences about the presence of short period planets allowed the efficient discovery of the DMPP-1 planets. They have great potential for novel and informative follow-up characterisation studies.

Synoptic searches for the optical counterpart to a binary neutron star (BNS) or neutron star-black hole (NSBH) merger can pose significant challenges towards the discovery of kilonovae and performing multi-messenger science. In this work, we describe the advantage of a global multi-telescope network towards this end, with a particular focus on the key and complementary role the Dark Energy Camera (DECam) plays in multi-facility follow-up. We describe the Global Relay of Observatories Watching Transients Happen (GROWTH) Target-of-Opportunity (ToO) Marshal, a common web application we built to ingest events, plan observations, search for transient candidates, and retrieve performance summary statistics for all of the telescopes in our network. Our infrastructure enabled us to conduct observations of two events during O3a, S190426c and S190510g. Furthermore, our analysis of deep DECam observations of S190814bv conducted by the DESGW team, and access to a variety of global follow-up facilities allowed us to place meaningful constraints on the parameters of the kilonova and the merging binary. We emphasize the importance of a global telescope network in conjunction with a power telescope like DECam in performing searches for the counterparts to gravitational-wave sources.

Conditions on Saturn's moon Titan suggest dust devils, which are convective, dust-laden plumes, may be active. Although the exact nature of dust on Titan is unclear, previous observations confirm an active aeolian cycle, and dust devils may play an important role in Titan's aeolian cycle, possibly contributing to regional transport of dust and even production of sand grains. The Dragonfly mission to Titan will document dust devil and convective vortex activity and thereby provide a new window into these features, and our analysis shows that associated winds are likely to be modest and pose no hazard to the mission.

The nearby ultra-compact multiplanetary system YZ Ceti consists of at least three planets. The orbital period of each planet is the subject of discussion in the literature due to strong aliasing in the radial velocity data. The stellar activity of this M dwarf also hampers significantly the derivation of the planetary parameters. With an additional 229 radial velocity measurements obtained since the discovery publication, we reanalyze the YZ Ceti system and resolve the alias issues. We use model comparison in the framework of Bayesian statistics and periodogram simulations based on a method by Dawson and Fabrycky to resolve the aliases. We discuss additional signals in the RV data, and derive the planetary parameters by simultaneously modeling the stellar activity with a Gaussian process regression model. To constrain the planetary parameters further we apply a stability analysis on our ensemble of Keplerian fits. We resolve the aliases: the three planets orbit the star with periods of $2.02$ d, $3.06$ d, and $4.66$ d. We also investigate an effect of the stellar rotational signal on the derivation of the planetary parameters, in particular the eccentricity of the innermost planet. Using photometry we determine the stellar rotational period to be close to $68$ d. From the absence of a transit event with TESS, we derive an upper limit of the inclination of $i_\mathrm{max} = 87.43$ deg. YZ Ceti is a prime example of a system where strong aliasing hindered the determination of the orbital periods of exoplanets. Additionally, stellar activity influences the derivation of planetary parameters and modeling them correctly is important for the reliable estimation of the orbital parameters in this specific compact system. Stability considerations then allow additional constraints to be placed on the planetary parameters.

We report on a multi-band variability and correlation study of the TeV blazar Mrk 421 during an exceptional flaring activity observed from 2013 April 11 to 2013 April 19. The study uses, among others, data from GASP-WEBT, Swift, NuSTAR, Fermi-LAT, VERITAS, and MAGIC. The large blazar activity, and the 43 hours of simultaneous NuSTAR and MAGIC/VERITAS observations, permitted variability studies on 15 minute time bins, and over three X-ray bands (3-7 keV, 7-30 keV and 30-80 keV) and three very-high-energy (>0.1 TeV, hereafter VHE) gamma-ray bands (0.2-0.4 TeV, 0.4-0.8 TeV and >0.8 TeV). We detected substantial flux variations on multi-hour and sub-hour timescales in all the X-ray and VHE gamma-ray bands. The characteristics of the sub-hour flux variations are essentially energy-independent, while the multi-hour flux variations can have a strong dependence on the energy of the X-ray and the VHE gamma rays. The three VHE bands and the three X-ray bands are positively correlated with no time-lag, but the strength and the characteristics of the correlation changes substantially over time and across energy bands. Our findings favour multi-zone scenarios for explaining the achromatic/chromatic variability of the fast/slow components of the light curves, as well as the changes in the flux-flux correlation on day-long timescales. We interpret these results within a magnetic reconnection scenario, where the multi-hour flux variations are dominated by the combined emission from various plasmoids of different sizes and velocities, while the sub-hour flux variations are dominated by the emission from a single small plasmoid moving across the magnetic reconnection layer.

Some highly irradiated close-in exoplanets orbit stars showing anomalously low stellar chromospheric emission. We attribute this to absorption by circumstellar gas replenished by mass loss from ablating planets. Here we report statistics validating this hypothesis. Among ~3000 nearby, bright, main sequence stars ~40 show depressed chromospheric emission indicative of undiscovered mass-losing planets. The Dispersed Matter Planet Project uses high precision, high cadence radial velocity measurements to detect these planets. We summarise results for two planetary systems (DMPP-1 and DMPP-3) and fully present observations revealing a Mp sin i = 0.469 M$_{\rm J}$ planet in a 5.207 d orbit around the $\gamma$-Doradus pulsator HD 11231 (DMPP-2). We have detected short period planets wherever we have made more than 60 RV measurements, demonstrating that we have originated a very efficient method for detecting nearby compact planetary systems. Our shrouded, ablating planetary systems may be a short-lived phase related to the Neptunian desert: i.e. the dearth of intermediate-mass planets at short orbital periods. The circumstellar gas facilitates compositional analysis; allowing empirical exogeology in the cases of sublimating rocky planets. Dispersed Matter Planet Project discoveries will be important for establishing the empirical mass-radius-composition relationship(s) for low mass planets.

We present the first Doppler images of the prototypical active binary star RS CVn, derived from high-resolution spectra observed in 2004, 2016 and 2017, using three different telescopes and observing sites. We apply the least-squares deconvolution technique to all observed spectra to obtain high signal-to-noise line profiles, which are used to derive the surface images of the active K-type component. Our images show a complex spot pattern on the K star, distributed widely in longitude. All starspots revealed by our Doppler images are located below a latitude of about 70$^{\circ}$. In accordance with previous light-curve modeling studies, we find no indication of a polar spot on the K star. Using Doppler images derived from two consecutive rotational cycles, we estimate a surface differential rotation rate of $\Delta\Omega = -0.039 \pm 0.003 ~rad~d^{-1}$ and $\alpha = \Delta\Omega/\Omega_{eq} = -0.030 \pm 0.002$ for the K star. Given the limited phase coverage during those two rotations, the uncertainty of our differential rotation estimate is presumably higher.

Earth mass exoplanets are difficult to detect. The Dispersed Matter Planet Project (DMPP) identifies stars which are likely to host the most detectable low mass exoplanets. The star DMPP-3 (HD 42936) shows signs of circumstellar absorption, indicative of mass loss from ablating planets. Here we report the radial velocity (RV) discovery of a highly eccentric 507 d binary companion and a hot super-Earth planet in a 6.67 d orbit around the primary star. DMPP-3A is a solar type star while DMPP-3B is just massive enough to fuse hydrogen. The binary, with semi-major axis 1.22 $\pm$ 0.02 AU, is considerably tighter than others known to host planets orbiting only one of the component stars. The configuration of the DMPP-3 planetary system is rare and indicates dynamical interactions, though the evolutionary history is not entirely clear. DMPP-3Ab is possibly the residual core of a giant planet precursor, consistent with the inferred circumstellar gas shroud.

We present a sample of 446 galaxy pairs constructed using the cosmological simulation IllustrisTNG-100 at z = 0, with M$_{FoF, dm}$ = 10$^{11}$-10$^{13.5}$ M$_{\odot}$. We produce ideal mock SDSS g-band images of all pairs to test the reliability of visual classification schema employed to produce samples of interacting galaxies. We visually classify each image as interacting or not based on the presence of a close neighbour, the presence of stellar debris fields, disturbed discs, and/or tidal features. By inspecting the trajectories of the pairs, we determine that these indicators correctly identify interacting galaxies $\sim45\%$ of the time. We subsequently split the sample into the visually identified interacting pairs (VIP; 38 pairs) and those which are interacting but are not visually identified (nonVIP; 47 pairs). We find that VIP have undergone a close passage nearly twice as recently as the nonVIP, and typically have higher stellar masses. Further, the VIP sit in dark matter haloes that are approximately 2.5 times as massive, in environments nearly 2 times as dense, and are almost a factor of 10 more affected by the tidal forces of their surroundings than the nonVIP. These factors conspire to increase the observability of tidal features and disturbed morphologies, making the VIP more likely to be identified. Thus, merger rate calculations which rely on stellar morphologies are likely to be significantly biased toward massive galaxy pairs which have recently undergone a close passage.

Small low-mass stars are favourable targets for the detection of rocky habitable planets. In particular, planetary systems in the solar neighbourhood are interesting and suitable for precise characterisation. The Red Dots campaigns seek to discover rocky planets orbiting nearby low-mass stars. The 2018 campaign targeted GJ 1061, which is the 20$^{\rm th}$ nearest star to the Sun. For three consecutive months we obtained nightly, high-precision radial velocity measurements with the HARPS spectrograph. We analysed these data together with archival HARPS data. We report the detection of three planet candidates with periods of $3.204\pm 0.001$, $6.689\pm 0.005$ and $13.03\pm 0.03$ days, which is close to 1:2:4 period commensurability. After several considerations related to the properties of the noise and sampling, we conclude that a 4$^{\rm th}$ signal is most likely explained by stellar rotation, although it may be due to a planet. The proposed three-planet system (and the potential four-planet solution) is long-term dynamically stable. Planet-planet gravitational interactions are below our current detection threshold. The minimum masses of the three planets range from $1.4\pm 0.2$ to $1.8\pm 0.3$ M$_\oplus$. Planet d, with $m \sin i = 1.68\pm0.25$ M$_\oplus$, receives a similar amount of energy as Earth receives from the Sun. Consequently it lies within the liquid-water habitable zone of the star and has a similar equilibrium temperature to Earth. GJ 1061 has very similar properties to Proxima Centauri but activity indices point to lower levels of stellar activity.

We measure the bulk system parameters of the seismically active, rapidly-rotating $\delta$-Scuti KOI-976 and constrain the orbit geometry of its transiting binary companion using a combined approach of asteroseismology and gravity-darkening light curve analysis. KOI-976 is a $1.62\pm0.2~\mathrm{M_\odot}$ star with a measured $v\sin(i)$ of $120\pm2$ km/s and seismically-induced variable signal that varies by $\sim$ 0.6\% of the star's total photometric brightness. We take advantage of the star's oblate shape and seismic activity to perform three measurements of its obliquity angle relative to the plane of the sky. We first apply rotational splitting theory to the star's variable signal observed in short-cadence \emph{Kepler} photometry to constrain KOI-976's obliquity angle, and then subtract off variability from that dataset using the linear algorithm for significance reduction software {\tt LASR}. We perform gravity-darkened fits to \emph{Kepler} variability-subtracted short-cadence photometry and to \emph{Kepler's} phase-folded long-cadence photometry to obtain two more measurements of the star's obliquity. We find that the binary system transits in a grazing configuration with measured obliquity values of $36^\circ\pm17^\circ$, $46^\circ\pm16^\circ$, and $43^\circ\pm20^\circ$ respectively for the three measurements. We perform these analyses as a way to demonstrate overcoming the challenges high-mass stars can present to transit light curve fitting and to prepare for the large number of exoplanets \emph{TESS} will discover orbiting A/F stars.

The most abundant stars in the Galaxy, M dwarfs, are very commonly hosts to diverse systems of low-mass planets. Their abundancy implies that the general occurrence rate of planets is dominated by their occurrence rate around such M dwarfs. In this article, we combine the M dwarf surveys conducted with the HIRES/Keck, PFS/Magellan, HARPS/ESO, and UVES/VLT instruments supported with data from several other instruments. We analyse the radial velocities of an approximately volume- and brightness-limited sample of 426 nearby M dwarfs in order to search for Doppler signals of cadidate planets. In addition, we analyse spectroscopic activity indicators and ASAS photometry to rule out radial velocity signals corresponding to stellar activity as Doppler signals of planets. We calculate estimates for the occurrence rate of planets around the sample stars and study the properties of this occurrence rate as a function of stellar properties. Our analyses reveal a total of 118 candidate planets orbiting nearby M dwarfs. Based on our results accounting for selection effects and sample detection threshold, we estimate that M dwarfs have on average at least 2.39$^{+4.58}_{-1.36}$ planets per star orbiting them. Accounting for the different sensitivities of radial velocity surveys and Kepler transit photometry implies that there are at least 3.0 planets per star orbiting M dwarfs. We also present evidence for a population of cool mini-Neptunes and Neptunes with indications that they are found an order of magnitude more frequently orbiting the least massive M dwarfs in our sample.

The American Astronomical Society (AAS) Journals are a vital asset of our professional society. With the push towards open access, page charges are a viable and sustainable option for continuing to effectively fund and publish the AAS Journals. However, the existing page charge model, which requires individual authors to pay page charges out of their grants or even out of pocket, is already challenging to some researchers and could be exacerbated in the Open Access (OA) era if charges increase. A discussion of alternative models for funding page charges and publishing costs should be part of the Astro2020 decadal survey if we wish to continue supporting the sustainable and accessible publication of US research in AAS journals in the rapidly-shifting publication landscape. The AAS Publications Committee recommends that the National Academy of Sciences form a task force to develop solutions and recommendations with respect to the urgent concerns and considerations highlighted in this White Paper.

Variations in the axial tilt, or obliquity, of terrestrial planets can affect their climates and therefore their habitability. Kepler-62f is a 1.4 R$_\oplus$ planet orbiting within the habitable zone of its K2 dwarf host star (Borucki et al. 2013). We perform N-body simulations that monitor the evolution of obliquity of Kepler-62f for 10 million year timescales to explore the effects on model assumptions, such as the masses of the Kepler-62 planets and the possibility of outer bodies. Significant obliquity variation occurs when the rotational precession frequency overlaps with one or more of the secular orbital frequencies, but most variations are limited to $\lesssim$10$^\circ$. Moderate variations ($\sim$10$^\circ - 20^\circ$) can occur over a broader range of initial obliquities when the relative nodal longitude ($\Delta\Omega$) overlaps with the frequency and phase of a given secular mode. However, we find that adding outer gas giants on long period orbits ($\gtrsim$ 1000 days) can produce large ($\sim$60$^\circ$) variations in obliquity if Kepler-62f has a very rapid (4 hr) rotation period. The possibility of giant planets on long period orbits impacts the climate and habitability of Kepler-62f through variations in the latitudinal surface flux, where the timescale for large variation can occur on million year timescales.

The first two months of the third Advanced LIGO and Virgo observing run (2019 April-May) showed that distant gravitational wave (GW) events can now be readily detected. Three candidate mergers containing neutron stars (NS) were reported in a span of 15 days, all likely located more than 100 Mpc away. However, distant events such as the three new NS mergers are likely to be coarsely localized, which highlights the importance of facilities and scheduling systems that enable deep observations over hundreds to thousands of square degrees to detect the electromagnetic counterparts. On 2019-05-10 02:59:39.292 UT the GW candidate S190510g was discovered and initially classified as a BNS merger with 98% probability. The GW event was localized within an area of 3462 deg2, later refined to 1166 deg2 (90%) at a distance of 227 +- 92 Mpc. We triggered Target of Opportunity observations with the Dark Energy Camera (DECam), a wide-field optical imager mounted at the prime focus of the 4m Blanco Telescope at CTIO in Chile. This Letter describes our DECam observations and our real-time analysis results, focusing in particular on the design and implementation of the observing strategy. Within 24 hours of the merger time, we observed 65% of the total enclosed probability of the final skymap with an observing efficiency of 94%. We identified and publicly announced 13 candidate counterparts. S190510g was re-classified 1.7 days after the merger, after our observations were completed, with a "binary neutron star merger" probability reduced from 98% to 42% in favor of a "terrestrial" classification.

The discovery of a binary neutron star merger (NSM) through both its gravitational wave and electromagnetic emission has revealed these events to be key sites of r-process nucleosynthesis. Here, we evaluate the prospects of finding the remnants of Galactic NSMs by detecting the gamma-ray decay lines from their radioactive r-process ejecta. We find that $^{126}$Sn, which has several lines in the energy range 415-695 keV and resides close to the second r-process peak, is the most promising isotope, because of its half-life $t_{1/2}=2.30(14)\times 10^{5}$ yr being comparable to the ages of recent NSMs. Using a Monte Carlo procedure, we predict that multiple remnants are detectable as individual sources by next-generation gamma-ray telescopes which achieve sub-MeV line sensitivities of $\sim 10^{-8}$-$10^{-6}$ $\gamma$ cm$^{-2}$ s$^{-1}$. However, given the unknown locations of the remnants, the most promising search strategy is a systematic survey of the Galactic plane and bulge extending to high Galactic latitudes. Individual known supernova remnants which may be mis-classified NSM remnants could also be targeted, especially those located outside the Galactic plane. Detection of a moderate sample of Galactic NSM remnants would provide important clues to unresolved issues such as the production of actinides in NSMs, properties of merging NS binaries, and even help distinguish them from rare supernovae as current Galactic r-process sources. We also investigate the diffuse flux from longer-lived nuclei (e.g. $^{182}$Hf) that could in principle trace the Galactic spatial distribution of NSMs over longer timescales, but find that the detection of the diffuse flux appears challenging even with next-generation telescopes.

As of today, we have directly detected exactly one source in both gravitational waves (GWs) and electromagnetic (EM) radiation, the binary neutron star merger GW170817, its associated gamma-ray burst GRB170817A, and the subsequent kilonova SSS17a/AT 2017gfo. Within ten years, we will detect hundreds of events, including new classes of events such as neutron-star-black-hole mergers, core-collapse supernovae, and almost certainly something completely unexpected. As we build this sample, we will explore exotic astrophysical topics ranging from nucleosynthesis, stellar evolution, general relativity, high-energy astrophysics, nuclear matter, to cosmology. The discovery potential is extraordinary, and investments in this area will yield major scientific breakthroughs. Here we outline some of the most exciting scientific questions that can be answered by combining GW and EM observations.

The field of time-domain astrophysics has entered the era of Multi-messenger Astronomy (MMA). One key science goal for the next decade (and beyond) will be to characterize gravitational wave (GW) and neutrino sources using the next generation of Extremely Large Telescopes (ELTs). These studies will have a broad impact across astrophysics, informing our knowledge of the production and enrichment history of the heaviest chemical elements, constrain the dense matter equation of state, provide independent constraints on cosmology, increase our understanding of particle acceleration in shocks and jets, and study the lives of black holes in the universe. Future GW detectors will greatly improve their sensitivity during the coming decade, as will near-infrared telescopes capable of independently finding kilonovae from neutron star mergers. However, the electromagnetic counterparts to high-frequency (LIGO/Virgo band) GW sources will be distant and faint and thus demand ELT capabilities for characterization. ELTs will be important and necessary contributors to an advanced and complete multi-messenger network.

The kilonova emission observed following the binary neutron star merger event GW170817 provided the first direct evidence for the synthesis of heavy nuclei through the r-process. The late-time transition in the spectral energy distribution to near-infrared wavelengths was interpreted as indicating the production of lanthanide nuclei, with atomic mass number A>140. However, compelling evidence for the presence of heavier third-peak r-process elements (e.g., gold, platinum) or translead nuclei remains elusive. At early times (~days) most of the r-process heating arises from a large statistical ensemble of beta-decays, which thermalize efficiently while the ejecta is still dense, generating a heating rate that is reasonably approximated by a single power-law. However, at later times (weeks to months), the decay energy input can possibly be dominated by a discrete number of alpha-decays, 223Ra (half-life t_{1/2}=11.43d), 225Ac (t_{1/2}=10.0d, following the beta-decay of 225Ra with t_{1/2}=14.9d), and the fissioning isotope 254Cf (t_{1/2}=60.5d), which liberate more energy per decay and thermalize with greater efficiency than beta-decay products. Late-time nebular observations of kilonovae which constrain the radioactive power provide the potential to identify signatures of these individual isotopes, thus confirming the production of heavy nuclei. In order to constrain the bolometric light to the required accuracy, multi-epoch & wide-band observations are required with sensitive instruments like the James Webb Space Telescope. In addition, by comparing the nuclear heating rate obtained with an abundance distribution following the Solar r abundance pattern, to the bolometric lightcurve of AT2017gfo, we find that the yet-uncertain r abundance of 72Ge plays a decisive role in powering the lightcurve, if one assumes that GW170817 has produced a full range of the Solar r abundances down to A~70.

At a distance of 1.8 parsecs, Barnard's star (Gl 699) is a red dwarf with the largest apparent motion of any known stellar object. It is the closest single star to the Sun, second only to the alpha Centauri triple stellar system. Barnard's star is also among the least magnetically active red dwarfs known and has an estimated age older than our Solar System. Its properties have made it a prime target for planet searches employing techniques such as radial velocity, astrometry, and direct imaging, all with different sensitivity limits but ultimately leading to disproved or null results. Here we report that the combination of numerous measurements from high-precision radial velocity instruments reveals the presence of a low-amplitude but significant periodic signal at 233 days. Independent photometric and spectroscopic monitoring, as well as the analysis of instrumental systematic effects, show that this signal is best explained as arising from a planetary companion. The candidate planet around Barnard's star is a cold super-Earth with a minimum mass of 3.2 Earth masses orbiting near its snow-line. The combination of all radial velocity datasets spanning 20 years additionally reveals a long-term modulation that could arise from a magnetic activity cycle or from a more distant planetary object. Because of its proximity to the Sun, the proposed planet has a maximum angular separation of 220 milli-arcseconds from Barnard's star, making it an excellent target for complementary direct imaging and astrometric observations.

We present an analysis of gas densities in the central R=300 parsecs of the Milky Way, focusing on three clouds: GCM-0.02-0.07 (the 50 km/s cloud), GCM-0.13-0.08 (the 20 km/s cloud), and GCM0.25+0.01 (the "Brick"). Densities are determined using observations of the J=(3-2), (4-3), (5-4), (10-9), (18-17), (19-18), (21-20), and (24-23) transitions of the molecule HC3N. We find evidence of at least two excitation regimes for HC3N and constrain the low-excitation component to have a density less than 10^4 cm^-3 and the high-excitation component to have a density between 10^5 and 10^6 cm^-3. This is much less than densities of 10^7 cm^-3 that are found in Sgr B2, the most actively star-forming cloud in the Galactic center. This is consistent with the requirement of a higher density threshold for star formation in the Galactic center than is typical in the Galactic disk. We are also able to constrain the column density of each component in order to determine the mass fraction of 'dense' (n>10^5 cm^-3) gas for these clouds. We find that this is ~15% for all three clouds. Applying the results of our models to ratios of the (10-9) and (3-2) line across the entire central R=300 pc, we find that the fraction of gas with n>10^4 cm^-3 increases inward of a radius of ~140 pc, consistent with the predictions of recent models for the gas dynamics in this region. Our observations show that HC3N is an excellent molecule for probing the density structure of clouds in the Galactic center.

The binary neutron star merger GW170817 detected by Advanced LIGO/Virgo was followed by a short gamma-ray burst (GRB) and thermal radiation ('kilonova') powered by the radioactive decay of heavy nuclei synthesized in the merger ejecta by the rapid neutron capture process ('r-process'). The large inferred quantity of ejecta is best understood as originating in an outflow from the accretion disk surrounding the newly-formed black hole, the same engine that was likely responsible for the GRB jet. Similar accretion flows accompany the collapse of rotating massive stars ('collapsars'), powering the class of long GRBs and their associated supernovae. Here we show that collapsar accretion disks also produce neutron-rich outflows that synthesize heavy r-process nuclei, despite the comparatively proton-rich composition of the infalling star. Though occurring less frequently than mergers, the much greater accreted mass in collapsars---and their correspondingly larger disk wind ejecta---implicate them as dominant contributors to the Galactic r-process. Collapsars provide the rare r-process source needed to promptly enrich a small fraction of ultra-faint dwarf galaxies early in the Universe which is also compatible with the Galactic chemical enrichment of europium relative to iron over longer timescales. We predict an excess in the late-time near-infrared emission from GRB supernovae, testable by future observations.

We have analyzed the complete Visual and Infrared Mapping Spectrometer (VIMS) data archive of Titan. Our objective is to build global surface cartographic products, by combining all the data gathered during the 127 targeted flybys of Titan into synthetic global maps interpolated on a grid at 32 pixels per degree (~1.4 km/pixel at the equator), in seven infrared spectral atmospheric windows. Multispectral summary images have been computed for each single VIMS cube in order to rapidly identify their scientific content and assess their quality. These summary images are made available to the community on a public website (vims.univ-nantes.fr). The global mapping work faced several challenges due to the strong absorbing and scattering effects of the atmosphere coupled to the changing observing conditions linked to the orbital tour of the Cassini mission. We determined a surface photometric function which accounts for variations in incidence, emergence and phase angles, and which is able to mitigate brightness variations linked to the viewing geometry of the flybys. The atmospheric contribution has been reduced using the subtraction of the methane absorption band wings, considered as proxies for atmospheric haze scattering. We present a new global three color composite map of band ratios (red: 1.59/1.27 {\mu}m; green: 2.03/1.27 {\mu}m; blue: 1.27/1.08 {\mu}m), which has also been empirically corrected from an airmass (the solar photon path length through the atmosphere) dependence. This map provides a detailed global color view of Titan's surface partially corrected from the atmosphere and gives a global insight of the spectral variability, with the equatorial dunes fields appearing in brownish tones, and several occurrences of bluish tones localized in areas such as Sinlap, Menvra and Selk craters. (...)

In July 2018 an FRIB Theory Alliance program was held on the implications of GW170817 and its associated kilonova for r-process nucleosynthesis. Topics of discussion included the astrophysical and nuclear physics uncertainties in the interpretation of the GW170817 kilonova, what we can learn about the astrophysical site or sites of the r process from this event, and the advances in nuclear experiment and theory most crucial to pursue in light of the new data. Here we compile a selection of scientific contributions to the workshop, broadly representative of progress in r-process studies since the GW170817 event.

Compact object mergers can produce a thermal electromagnetic counterpart (a "kilonova") powered by the decay of freshly synthesized radioactive isotopes. The luminosity of kilonova light curves depends on the efficiency with which beta-decay electrons are thermalized in the ejecta. Here we derive a simple analytic solution for thermalization by calculating how electrons accumulate in the ejecta and lose energy adiabatically and via plasma losses. We find that the time-dependent thermalization efficiency is well described by $f(t) \approx (1 + t/t_e)^{-n}$ where $n \approx 1$ and the timescale $t_e$ is a function of the ejecta mass and velocity. For a statistical distribution of r-process isotopes with radioactive power $\dot{Q} \propto t^{-4/3}$, the late time kilonova luminosity asymptotes to $L \propto t^{-7/3}$ and depends super-linearly on the ejecta mass, $L \propto M^{5/3}$. If a kilonova is instead powered by a single dominate isotope, we show that the late time luminosity can deviate substantially from the underlying exponential decay and eventually become brighter than the instantaneous radioactivity due to the accumulation of trapped electrons. Applied to the kilonova associated with the gravitational wave source GW170817, these results imply that a possible steepening of the observed light curve at $\gtrsim 7$ days is unrelated to thermalization effects and instead could mark the onset of translucency in a high opacity component of ejecta. The analytic results should be convenient for estimating the properties of observed kilonovae and assessing the potential late time detectability of future events.

Neutron star mergers offer unique conditions for the creation of the heavy elements and additionally provide a testbed for our understanding of this synthesis known as the $r$-process. We have performed dynamical nucleosynthesis calculations and identified a single isotope, $^{254}$Cf, which has a particularly high impact on the brightness of electromagnetic transients associated with mergers on the order of 15 to 250 days. This is due to the anomalously long half-life of this isotope and the efficiency of fission thermalization compared to other nuclear channels. We estimate the fission fragment yield of this nucleus and outline the astrophysical conditions under which $^{254}$Cf has the greatest impact to the light curve. Future observations in the middle-IR which are bright during this regime could indicate the production of actinide nucleosynthesis.

Dynamical interactions between colliding spiral galaxies strongly affect the state and distribution of their interstellar gas. Observations indicate that interactions funnel gas toward the nuclei, fueling bursts of star formation and nuclear activity. To date, most numerical simulations of galaxy mergers have assumed that the gaseous and stellar disks initially have the same distribution and size. However, observations of isolated disk galaxies show that this is seldom the case; in fact, most spirals have as much or more gas beyond their optical radii as they do within. Can gas in such extended disks be efficiently transported to the nuclei during interactions? To address this question, we examine the effect of various parameters on the transport of gas to the nuclei of interacting galaxies. In addition to the relative radii of the gaseous and stellar disks, these parameters include the pericentric separation, disk orientation, fractional gas mass, presence of a bulge, treatment of gas thermodynamics, and the spatial resolution of the numerical simulation. We found that gas accumulates in most of our simulated nuclei, but the efficiency of inflow is largely dependent upon the encounter geometry. Dissipation alone is not enough to produce inflows; an efficient mechanism for extracting angular momentum from the gas is necessary. Several different mechanisms are seen in these experiments. Aside from mode-driven inflows (such as, but not limited to, bars) and ram-pressure sweeping, both of which have been previously described and well studied, we supply the first quantitative study of an often-seen process: the formation of massive gas clumps in Jeans-unstable tidal shocks, and their subsequent delivery to the nuclei via dynamical friction.

One of the aims of the BCool programme is to search for cycles in other stars and to understand how similar they are to the Sun. In this paper we aim to monitor the evolution of $\tau$ Boo's large-scale magnetic field using high-cadence observations covering its chromospheric activity maximum. For the first time, we detect a polarity switch that is in phase with $\tau$ Boo's 120 day chromospheric activity maximum and its inferred X-ray activity cycle maximum. This means that $\tau$ Boo has a very fast magnetic cycle of only 240 days. At activity maximum $\tau$ Boo's large-scale field geometry is very similar to the Sun at activity maximum: it is complex and there is a weak dipolar component. In contrast, we also see the emergence of a strong toroidal component which has not been observed on the Sun, and a potentially overlapping butterfly pattern where the next cycle begins before the previous one has finished.

We have analyzed the complete archive of the Visual and Infrared Mapping Spectrometer (VIMS) data in order to monitor and analyze the evolution of the clouds and haze coverage at both poles of Titan during the entire Cassini mission. Our objective is to give a cartographic synopsis from a VIMS perspective, to provide a global view of the seasonal evolution of Titan's atmosphere over the poles. We leave the detailed comparison with the Imaging Science Subsystem (ISS) and the Composite Infrared Spectrometer (CIRS) data sets to further studies. We have computed global hyperspectral mosaics for each of the 127 targeted flybys of Titan to produce synthetic color maps emphasizing the main atmospheric features. The north pole appears fully covered by a huge cloud as soon as the first observations in 2004 and up to the equinox in 2009 (Le Mou\'elic et al. 2012). The northern skies then became progressively clearer, after the circulation turnover in 2009, revealing the underlying lakes and seas to the optical instruments up to 2017. The reverse situation is observed over the south pole, which was mostly clear of such a high obscuring cloud during the first years of the mission, but started to develop a polar cloud in 2012. This feature grew up month after month until the end of the mission in 2017, with a poleward latitudinal extent of 75$^\circ$S in 2013 up to 58$^\circ$S in April 2017. Thanks to the spectral capabilities of VIMS, we have detected HCN spectral signatures over the north pole in almost all flybys between 2004 and 2008. These HCN signatures started then to show up over the south pole in almost all flybys between 2012 and 2017, so perfectly matching the timing and spatial extent of the northern and southern polar atmospheric features.

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.

We develop a technique for removing stellar variability in the light curves of $\delta$-Scuti and similar stars. Our technique, which we name the Linear Algorithm for Significance Reduction (LASR), subtracts oscillations from a time series by minimizing their statistical significance in frequency space. We demonstrate that LASR can subtract variable signals of near-arbitrary complexity and can robustly handle close frequency pairs and overtone frequencies. We demonstrate that our algorithm performs an equivalent fit as prewhitening to the straightforward variable signal of KIC 9700322. We also show that LASR provides a better fit to seismic activity than prewhitening in the case of the complex $\delta$-Scuti KOI-976.

We explore the effect of using different kinematic tracers (HI and H$\alpha$) on reconstructing the encounter parameters of the Mice major galaxy merger (NGC 4676A/B). We observed the Mice using the SparsePak Integral Field Unit (IFU) on the WIYN telescope, and compared the H$\alpha$ velocity map with VLA HI observations. The relatively high spectral resolution of our data (R $\approx$ 5000) allows us to resolve more than one kinematic component in the emission lines of some fibers. We separate the H$\alpha$-[N II] emission of the star-forming regions from shocks using their [N II]/H$\alpha$ line ratio and velocity dispersion. We show that the velocity of star-forming regions agree with that of the cold gas (HI), particularly, in the tidal tails of the system. We reconstruct the morphology and kinematics of these tidal tails utilizing an automated modeling method based on the Identikit software package. We quantify the goodness of fit and the uncertainties of the derived encounter parameters. Most of the initial conditions reconstructed using H$\alpha$ and HI are consistent with each other, and qualitatively agree with the results of previous works. For example, we find 210$\pm^{50}_{40}$ Myrs, and 180$\pm^{50}_{40}$ Myrs for the time since pericenter, when modeling H$\alpha$ and HI kinematics, respectively. This confirms that in some cases, H$\alpha$ kinematics can be used instead of HI kinematics for reconstructing the initial conditions of galaxy mergers, and our automated modeling method is applicable to some merging systems.

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

The young and magnetically active K dwarf Epsilon Eridani exhibits a chromospheric activity cycle of about 3 years. Previous reconstructions of its large-scale magnetic field show strong variations at yearly epochs. To understand how Epsilon Eridani's large-scale magnetic field geometry evolves over its activity cycle we focus on high cadence observations spanning 5 months at its activity minimum. Over this timespan we reconstruct 3 maps of Epsilon Eridani's large-scale magnetic field using the tomographic technique of Zeeman Doppler Imaging. The results show that at the minimum of its cycle, Epsilon Eridani's large-scale field is more complex than the simple dipolar structure of the Sun and 61 Cyg A at minimum. Additionally we observe a surprisingly rapid regeneration of a strong axisymmetric toroidal field as Epsilon Eridani emerges from its S-index activity minimum. Our results show that all stars do not exhibit the same field geometry as the Sun and this will be an important constraint for the dynamo models of active solar-type stars.

We present an implementation of the Gehrels et al. (2016) galaxy-targeted strategy for gravitational-wave (GW) follow-up using the Las Cumbres Observatory global network of telescopes. We use the Galaxy List for the Advanced Detector Era (GLADE) galaxy catalog, which we show is complete (with respect to a Schechter function) out to ~300 Mpc for galaxies brighter than the median Schechter function galaxy luminosity. We use a prioritization algorithm to select the galaxies with the highest chance of containing the counterpart given their luminosity, their position, and their distance relative to a GW localization, and in which we are most likely to detect a counterpart given its expected brightness compared to the limiting magnitude of our telescopes. This algorithm can be easily adapted to any expected transient parameters and telescopes. We implemented this strategy during the second Advanced Detector Observing Run (O2) and followed the black hole merger GW170814 and the neutron star merger GW170817. For the latter, we identified an optical kilonova/macronova counterpart thanks to our algorithm selecting the correct host galaxy fifth in its ranked list among 182 galaxies we identified in the Laser Interferometer Gravitational-wave Observatory LIGO-Virgo localization. This also allowed us to obtain some of the earliest observations of the first optical transient ever triggered by a GW detection (as presented in a companion paper).

We present the spectroscopic evolution of AT 2017gfo, the optical counterpart of the first binary neutron star (BNS) merger detected by LIGO and Virgo, GW170817. While models have long predicted that a BNS merger could produce a kilonova (KN), we have not been able to definitively test these models until now. From one day to four days after the merger, we took five spectra of AT 2017gfo before it faded away, which was possible because it was at a distance of only 39.5 Mpc in the galaxy NGC 4993. The spectra evolve from blue ($\sim6400$K) to red ($\sim3500$K) over the three days we observed. The spectra are relatively featureless --- some weak features exist in our latest spectrum, but they are likely due to the host galaxy. However, a simple blackbody is not sufficient to explain our data: another source of luminosity or opacity is necessary. Predictions from simulations of KNe qualitatively match the observed spectroscopic evolution after two days past the merger, but underpredict the blue flux in our earliest spectrum. From our best-fit models, we infer that AT 2017gfo had an ejecta mass of $0.03M_\odot$, high ejecta velocities of $0.3c$, and a low mass fraction $\sim10^{-4}$ of high-opacity lanthanides and actinides. One possible explanation for the early excess of blue flux is that the outer ejecta is lanthanide-poor, while the inner ejecta has a higher abundance of high-opacity material. With the discovery and follow-up of this unique transient, combining gravitational-wave and electromagnetic astronomy, we have arrived in the multi-messenger era.

The merger of two neutron stars has been predicted to produce an optical-infrared transient (lasting a few days) known as a 'kilonova', powered by the radioactive decay of neutron-rich species synthesized in the merger. Evidence that short gamma-ray bursts also arise from neutron-star mergers has been accumulating. In models of such mergers a small amount of mass ($10^{-4}$-$10^{-2}$ solar masses) with a low electron fraction is ejected at high velocities (0.1-0.3 times light speed) and/or carried out by winds from an accretion disk formed around the newly merged object. This mass is expected to undergo rapid neutron capture (r-process) nucleosynthesis, leading to the formation of radioactive elements that release energy as they decay, powering an electromagnetic transient. A large uncertainty in the composition of the newly synthesized material leads to various expected colours, durations and luminosities for such transients. Observational evidence for kilonovae has so far been inconclusive as it was based on cases of moderate excess emission detected in the afterglows of gamma-ray bursts. Here we report optical to near-infrared observations of a transient coincident with the detection of the gravitational-wave signature of a binary neutron-star merger and of a low-luminosity short-duration gamma-ray burst. Our observations, taken roughly every eight hours over a few days following the gravitational-wave trigger, reveal an initial blue excess, with fast optical fading and reddening. Using numerical models, we conclude that our data are broadly consistent with a light curve powered by a few hundredths of a solar mass of low-opacity material corresponding to lanthanide-poor (a fraction of $10^{-4.5}$ by mass) ejecta.

The cosmic origin of the elements heavier than iron has long been uncertain. Theoretical modelling shows that the matter that is expelled in the violent merger of two neutron stars can assemble into heavy elements such as gold and platinum in a process known as rapid neutron capture (r-process) nucleosynthesis. The radioactive decay of isotopes of the heavy elements is predicted to power a distinctive thermal glow (a 'kilonova'). The discovery of an electromagnetic counterpart to the gravitational-wave source GW170817 represents the first opportunity to detect and scrutinize a sample of freshly synthesized r-process elements. Here we report models that predict the detailed electromagnetic emission of kilonovae and enable the mass, velocity and composition of ejecta to be derived from the observations. We compare the models to the optical and infrared radiation associated with GW170817 event to argue that the observed source is a kilonova. We infer the presence of two distinct components of ejecta, one composed primarily of light (atomic mass number less than 140) and one of heavy (atomic mass number greater than 140) r-process elements. Inferring the ejected mass and a merger rate from GW170817 implies that such mergers are a dominant mode of r-process production in the Universe.

Merging neutron stars offer an exquisite laboratory for simultaneously studying strong-field gravity and matter in extreme environments. We establish the physical association of an electromagnetic counterpart EM170817 to gravitational waves (GW170817) detected from merging neutron stars. By synthesizing a panchromatic dataset, we demonstrate that merging neutron stars are a long-sought production site forging heavy elements by r-process nucleosynthesis. The weak gamma-rays seen in EM170817 are dissimilar to classical short gamma-ray bursts with ultra-relativistic jets. Instead, we suggest that breakout of a wide-angle, mildly-relativistic cocoon engulfing the jet elegantly explains the low-luminosity gamma-rays, the high-luminosity ultraviolet-optical-infrared and the delayed radio/X-ray emission. We posit that all merging neutron stars may lead to a wide-angle cocoon breakout; sometimes accompanied by a successful jet and sometimes a choked jet.

Unusually high velocities (< ~0.1c) and correspondingly high kinetic energies have been observed in a subset of Type Ic supernovae (so-called "broad-lined Ic" supernovae; SNe Ic-BL), prompting a search for a central engine model capable of generating such energetic explosions. A clue to the explosion mechanism may lie in the fact that all supernovae that accompany long-duration gamma-ray bursts belong to the SN Ic-BL class. Using a combination of two-dimensional relativistic hydrodynamics and radiation transport calculations, we demonstrate that the central engine responsible for long gamma-ray bursts can also trigger a SN Ic-BL. We find that a reasonable gamma-ray burst engine injected into a stripped Wolf-Rayet progenitor produces a relativistic jet with energy ~10^51 ergs, as well as a SN whose synthetic light curves and spectra are fully consistent with observed SNe Ic-BL during the photospheric phase. As a result of the jet's asymmetric energy injection, the SN spectra and light curves depend on viewing angle. The impact of viewing angle on the spectrum is particularly pronounced at early times, while the viewing angle dependence for the light curves (~10% variation in bolometric luminosity) persists throughout the photospheric phase.

The removal of noise typically correlated in time and wavelength is one of the main challenges for using the radial velocity method to detect Earth analogues. We analyze radial velocity data of tau Ceti and find robust evidence for wavelength dependent noise. We find this noise can be modeled by a combination of moving average models and "differential radial velocities". We apply this noise model to various radial velocity data sets for tau Ceti, and find four periodic signals at 20.0, 49.3, 160 and 642 d which we interpret as planets. We identify two new signals with orbital periods of 20.0 and 49.3 d while the other two previously suspected signals around 160 and 600 d are quantified to a higher precision. The 20.0 d candidate is independently detected in KECK data. All planets detected in this work have minimum masses less than 4$M_\oplus$ with the two long period ones located around the inner and outer edges of the habitable zone, respectively. We find that the instrumental noise gives rise to a precision limit of the HARPS around 0.2 m/s. We also find correlation between the HARPS data and the central moments of the spectral line profile at around 0.5 m/s level, although these central moments may contain both noise and signals. The signals detected in this work have semi-amplitudes as low as 0.3 m/s, demonstrating the ability of the radial velocity technique to detect relatively weak signals.

We aim to understand how stellar parameters such as mass and rotation impact the distribution of starspots on the stellar surface. To this purpose, we have used Doppler imaging to reconstruct the surface brightness distributions of three fully convective M dwarfs with similar rotation rates. We secured high cadence spectral time series observations of the 5.5 AU separation binary GJ 65, comprising GJ 65A (M5.5V, Prot = 0.24 d) and GJ 65B (M6V, Prot = 0.23 d). We also present new observations of GJ 791.2A (M4.5V, Prot = 0.31 d). Observations of each star were made on two nights with UVES, covering a wavelength range from 0.64 - 1.03 microns. The time series spectra reveal multiple line distortions, which we interpret as cool starspots and which are persistent on both nights suggesting stability on the timescale of 3 days. Spots are recovered with resolutions down to 8.3 degs at the equator. The global spot distributions for GJ 791.2A are similar to observations made a year earlier. Similar high latitude and circumpolar spot structure is seen on GJ791.2A and GJ 65A. However, they are surprisingly absent on GJ 65B, which instead reveals more extensive, larger, spots concentrated at intermediate latitudes. All three stars show small amplitude latitude-dependent rotation that is consistent with solid body rotation. We compare our measurements of differential rotation with previous Doppler imaging studies and discuss the results in the wider context of other observational estimates and recent theoretical predictions.

Stellar oscillations appear all across the Hertzsprung-Russell diagram. Recent theoretical studies support their existence also in the atmospheres of M dwarfs. These studies predict for them short periodicities ranging from 20~min to 3~h. Our Cool Tiny Beats (CTB) programme aims at finding these oscillations for the very first time. With this goal, CTB explores the short time domain of M dwarfs using radial velocity data from the HARPS-ESO and HARPS-N high-precision spectrographs. Here we present the results for the two most long-term stable targets observed to date with CTB, GJ~588 and GJ~699 (i.e. Barnard's star). In the first part of this work we detail the correction of several instrumental effects. These corrections are specially relevant when searching for sub-night signals. Results show no significant signals in the range where M dwarfs pulsations were predicted. However, we estimate that stellar pulsations with amplitudes larger than $\sim0.5\,\mathrm{m\,s}^{-1}$ can be detected with a 90% completeness with our observations. This result, along with the excess of power regions detected in the periodograms, open the possibility of non-resolved very low amplitude pulsation signals. Next generation more precise instrumentation would be required to detect such oscillations. However, the possibility of detecting pulsating M-dwarf stars with larger amplitudes is feasible due to the short size of the analysed sample. This motivates the need for completeness of the CTB survey.

We present 90 mas (37 pc) resolution ALMA imaging of Arp 220 in the CO (1-0) line and continuum at $\lambda = 2.6$ mm. The internal gas distribution and kinematics of both galactic nuclei are well-resolved for the first time. In the West nucleus, the major gas and dust emission extends out to 0.2\arcsec radius (74 pc); the central resolution element shows a strong peak in the dust emission but a factor 3 dip in the CO line emission. In this nucleus, the dust is apparently optically thick ($\tau_{\rm 2.6mm} \sim1$) at $\lambda = 2.6$ mm with a dust brightness temperature $\sim147$ K. The column of ISM at this nucleus is $\rm N_{H2} \geq 2\times10^{26}$ cm$^{-2}$, corresponding to $\sim$900 gr cm$^{-2}$. The East nucleus is more elongated with radial extent 0.3\arcsec or $\sim111$ pc. The derived kinematics of the nuclear disks provide a good fit to the line profiles, yielding the emissivity distributions, the rotation curves and velocity dispersions. In the West nucleus, there is evidence of a central Keplerian component requiring a central mass of $8\times10^8$ \msun. The intrinsic widths of the emission lines are $\Delta \rm v (FWHM)$ = 250 (West) and 120 (East) \kms. Given the very short dissipation timescales for turbulence ($\lesssim10^5$ yrs), we suggest that the line widths may be due to semi-coherent motions within the nuclear disks. The symmetry of the nuclear disk structures is impressive -- implying the merger timescale is significantly longer than the rotation period of the disks.

We measured the chromospheric activity of the four hot Jupiter hosts WASP-43, WASP-51/HAT-P-30, WASP-72 & WASP-103 to search for anomalous values caused by the close-in companions. The Mount Wilson Ca II H&K S-index was calculated for each star using observations taken with the Robert Stobie Spectrograph at the Southern African Large Telescope. The activity level of WASP-43 is anomalously high relative to its age and falls among the highest values of all known main sequence stars. We found marginal evidence that the activity of WASP-103 is also higher than expected from the system age. We suggest that for WASP-43 and WASP-103 star-planet interactions (SPI) may enhance the Ca II H&K core emission. The activity levels of WASP-51/HAT-P-30 and WASP-72 are anomalously low, with the latter falling below the basal envelope for both main sequence and evolved stars. This can be attributed to circumstellar absorption due to planetary mass loss, though absorption in the ISM may contribute. A quarter of known short period planet hosts exhibit anomalously low activity levels, including systems with hot Jupiters and low mass companions. Since SPI can elevate and absorption can suppress the observed chromospheric activity of stars with close-in planets, their Ca II H&K activity levels are an unreliable age indicator. Systems where the activity is depressed by absorption from planetary mass loss are key targets for examining planet compositions through transmission spectroscopy.

Accounting for stellar activity is a crucial component of the search for ever-smaller planets orbiting stars of all spectral types. We use Doppler imaging methods to demonstrate that starspot induced radial velocity variability can be effectively reduced for moderately rotating, fully convective stars. Using starspot distributions extrapolated from sunspot observations, we adopt typical M dwarf starspot distributions with low contrast spots to synthesise line profile distortions. The distortions are recovered using maximum entropy regularised fitting and the corresponding stellar radial velocities are measured. The procedure is demonstrated for a late-M star harbouring an orbiting planet in the habitable zone. The technique is effective for stars with vsini = 1-10 km/s, reducing the stellar noise contribution by factors of nearly an order of magnitude. With a carefully chosen observing strategy, the technique can be used to determine the stellar rotation period and is robust to uncertainties such as unknown stellar inclination. While demonstrated for late-type M stars, the procedure is applicable to all spectral types.

We report the discovery of eight new giant planets, and updated orbits for four known planets, orbiting dwarf and subgiant stars using the CORALIE, HARPS, and MIKE instruments as part of the Calan-Hertfordshire Extrasolar Planet Search. The planets have masses in the range 1.1-5.4MJs, orbital periods from 40-2900 days, and eccentricities from 0.0-0.6. They include a double-planet system orbiting the most massive star in our sample (HD147873), two eccentric giant planets (HD128356b and HD154672b), and a rare 14 Herculis analogue (HD224538b). We highlight some population correlations from the sample of radial velocity detected planets orbiting nearby stars, including the mass function exponential distribution, confirmation of the growing body of evidence that low-mass planets tend to be found orbiting more metal-poor stars than giant planets, and a possible period-metallicity correlation for planets with masses >0.1MJ, based on a metallicity difference of 0.16 dex between the population of planets with orbital periods less than 100 days and those with orbital periods greater than 100 days.

At a distance of 1.295 parsecs, the red-dwarf Proxima Centauri ($\alpha$ Centauri C, GL 551, HIP 70890, or simply Proxima) is the Sun's closest stellar neighbor and one of the best studied low-mass stars. It has an effective temperature of only $\sim$ 3050 K, a luminosity of $\sim$0.1 per cent solar, a measured radius of 0.14 R$_\odot$ and a mass of about 12 per cent the mass of the Sun. Although Proxima is considered a moderately active star, its rotation period is $\sim$ 83 days, and its quiescent activity levels and X-ray luminosity are comparable to the Sun's. New observations reveal the presence of a small planet orbiting Proxima with a minimum mass of 1.3~Earth masses and an orbital period of $\sim$11.2 days. Its orbital semi-major axis is $\sim0.05$ AU, with an equilibrium temperature in the range where water could be liquid on its surface.

We investigate retrieval of the stellar rotation signal for Proxima Centauri. We make use of high-resolution spectra taken with uves and harps of Proxima Centauri over a 13-year period as well as photometric observations of Proxima Centauri from asas and hst. We measure the H{\alpha} equivalent width and H{\alpha} index, skewness and kurtosis and introduce a method that investigates the symmetry of the line, the Peak Ratio, which appears to return better results than the other measurements. Our investigations return a most significant period of 82.6 $\pm$ 0.1 days, confirming earlier photometric results and ruling out a more recent result of 116.6 days which we conclude to be an alias induced by the specific harps observation times. We conclude that whilst spectroscopic H{\alpha} measurements can be used for period recovery, in the case of Proxima Centauri the available photometric measurements are more reliable. We make 2D models of Proxima Centauri to generate simulated H{\alpha}, finding that reasonable distributions of plage and chromospheric features are able to reproduce the equivalent width variations in observed data and recover the rotation period, including after the addition of simulated noise and flares. However the 2D models used fail to generate the observed variety of line shapes measured by the peak ratio. We conclude that only 3D models which incorporate vertical motions in the chromosphere can achieve this.

We present a geomorphologic map of Titan's polar terrains. The map was generated from a combination of Cassini Synthetic Aperture Radar (SAR) and Imaging Science Subsystem imaging products, as well as altimetry, SARTopo and radargrammetry topographic datasets. In combining imagery with topographic data, our geomorphologic map reveals a stratigraphic sequence from which we infer process interactions between units. In mapping both polar regions with the same geomorphologic units, we conclude that processes that formed the terrains of the north polar region also acted to form the landscape we observe at the south. Uniform, SAR-dark plains are interpreted as sedimentary deposits, and are bounded by moderately dissected uplands. These plains contain the highest density of filled and empty lake depressions, and canyons. These units unconformably overlay a basement rock that outcrops as mountains and SAR-bright dissected terrains at various elevations across both poles. All these units are then superposed by surficial units that slope towards the seas, suggestive of subsequent overland transport of sediment. From estimates of the depths of the embedded empty depressions and canyons that drain into the seas, the SAR-dark plains must be >600 m thick in places, though the thickness may vary across the poles. At the lowest elevations of each polar region, there are large seas, which are currently liquid methane/ethane filled at the north and empty at the south. The large plains deposits and the surrounding hillslopes may represent remnant landforms that are a result of previously vast polar oceans, where larger liquid bodies may have allowed for a sustained accumulation of soluble and insoluble sediments, potentially forming layered sedimentary deposits. Coupled with vertical crustal movements, the resulting layers would be of varying solubilities and erosional resistances.

We have detected excess absorption in the emission cores of Ca II H & K during transits of HD 189733b for the first time. Using observations of three transits we investigate the origin of the absorption, which is also seen in H{\alpha} and the Na I D lines. Applying differential spectrophotometry methods to the Ca II H and Ca II K lines combined, using respective passband widths of $\Delta\lambda$ = 0.4 & 0.6 \AA\ yields excess absorption of t$_d$ = 0.0074 $\pm$ 0.0044 (1.7$\sigma$; Transit 1) and 0.0214 +/- 0.0022 (9.8$\sigma$; Transit 2). Similarly, we detect excess H{\alpha} absorption in a passband of width $\Delta\lambda$ = 0.7 \AA, with t$_d$ = 0.0084 $\pm$ 0.0016 (5.2$\sigma$) and 0.0121 $\pm$ 0.0012 (9.9$\sigma$). For both lines, Transit 2 is thus significantly deeper. Combining all three transits for the Na I D lines yields excess absorption of t$_d$ = 0.0041 $\pm$ 0.0006 (6.5$\sigma$). By considering the time series observations of each line, we find that the excess apparent absorption is best recovered in the stellar reference frame. These findings lead us to postulate that the main contribution to the excess transit absorption in the differential light curves arises because the normalising continuum bands form in the photosphere, whereas the line cores contain a chromospheric component. We can not rule out that part of the excess absorption signature arises from the planetary atmosphere, but we present evidence which casts doubt on recent claims to have detected wind motions in the planet's atmosphere in these data.

One of the most promising electromagnetic signatures of compact object mergers are kilonovae: approximately isotropic radioactively-powered transients that peak days to weeks post-merger. Key uncertainties in modeling kilonovae include the emission profiles of the radioactive decay products---non-thermal beta- and alpha-particles, fission fragments, and gamma-rays---and the efficiency with which they deposit their energy in the ejecta. The total radioactive energy and the efficiency of its thermalization sets the luminosity budget and is therefore necessary for predicting kilonova light curves. We outline the uncertainties in r-process decay, describe the physical processes by which the energy of the decay products is absorbed in the ejecta, and present time-dependent thermalization efficiencies for each particle type. We determine the net heating efficiency and explore its dependence on r-process yields---in particular, the production of translead nuclei that undergo alpha-decay---and on the ejecta's mass, velocity, composition, and magnetic field configuration. We incorporate our results into new time-dependent, multi-wavelength radiation transport simulations, and calculate updated predictions of kilonova light curves. Thermalization has a substantial effect on kilonova photometry, reducing the luminosity by a factor of roughly 2 at peak, and by an order of magnitude or more at later times (15 days or more after explosion). We present simple analytic fits to time-dependent net thermalization efficiencies, which can easily be used to improve light curve models. We briefly revisit the putative kilonova that accompanied gamma ray burst 130603B, and offer new estimates of the mass ejected in that event. We find that later-time kilonova light curves can be significantly impacted by alpha-decay from translead isotopes; data at these times may therefore be diagnostic of ejecta abundances.

We present a new, detailed analysis of the morphologies and molecular gas fractions for a complete sample of 65 local luminous infrared galaxies (LIRGs) from the Great Observatories All-Sky LIRG Survey (GOALS) using high resolution $I$-band images from The Hubble Space Telescope, the University of Hawaii 2.2m Telescope and the Pan-STARRS1 Survey. Our classification scheme includes single undisturbed galaxies, minor mergers, and major mergers, with the latter divided into five distinct stages from pre-first pericenter passage to final nuclear coalescence. We find that major mergers of molecular gas-rich spirals clearly play a major role for all sources with $L_{\rm IR} > 10^{11.5} L_\odot $; however, below this luminosity threshold, minor mergers and secular processes dominate. Additionally, galaxies do not reach $L_{\rm IR} > 10^{12.0} L_\odot $ until late in the merger process when both disks are near final coalescence. The mean molecular gas fraction (MGF $= M_{\rm H_2} / (M_* + M_{\rm H_2})$) for non-interacting and early-stage major merger LIRGs is 18$\pm 2$%, which increases to 33$\pm 3$%, for intermediate stage major merger LIRGs, consistent with the hypothesis that, during the early-mid stages of major mergers, most of the initial large reservoir of atomic gas (HI) at large galactocentric radii is swept inward where it is converted into molecular gas (H$_2$).

Understanding the sources of instrumental systematic noise is a must to improve the design of future spectrographs. In this study, we alternated observations of the well-suited pair of M-stars GJ 725A+B to delve into the sub-night HARPS-N response. Besides the possible presence of a low-mass planet orbiting GJ 725B, our observations reveal changes in the spectral energy distribution (SED) correlated with measurements of the width of the instrumental line profile and, to a lower degree, with the Doppler measurements. To study the origin of these effects, we searched for correlations among several quantities defined and measured on the spectra and on the acquisition images. We find that the changes in apparent SED are very likely related to flux losses at the fibre input. Further tests indicate that such flux losses do not seriously affect the shape of the instrumental point spread function of HARPS-N, but identify an inefficient fitting of the continuum as the most likely source of the systematic variability observed in the FWHM. This index, accounting for the HARPS-N cross-correlation profiles width, is often used to decorrelate Doppler time-series. We show that the Doppler measurement obtained by a parametric least-squares fitting of the spectrum accounting for continuum variability is insensitive to changes in the slope of the SED, suggesting that forward modeling techniques to measure moments of the line profile are the optimal way to achieve higher accuracy. Remaining residual variability at ~1 m/s suggests that for M-stars Doppler surveys the current noise floor still has an instrumental origin.

Venus currently rotates slowly, with its spin controlled by solid-body and atmospheric thermal tides. However, conditions may have been far different 4 billion years ago, when the Sun was fainter and most of the carbon within Venus could have been in solid form, implying a low-mass atmosphere. We investigate how the obliquity would have varied for a hypothetical rapidly rotating Early Venus. The obliquity variation structure of an ensemble of hypothetical Early Venuses is simpler than that Earth would have if it lacked its large Moon (Lissauer et al., 2012), having just one primary chaotic regime at high prograde obliquities. We note an unexpected long-term variability of up to $\pm7^\circ$ for retrograde Venuses. Low-obliquity Venuses show very low total obliquity variability over billion-year timescales -- comparable to that of the real Moon-influenced Earth.

We document the similarities in composition between the equatorial basins Tui Regio, Hotei Regio, and other 5-$\mu$m-bright materials, notably the north polar evaporites, by investigating the presence and extent of an absorption feature at 4.92 $\mu$m. In most observations, Woytchugga Lacuna, Ontario Lacus, MacKay Lacus, deposits near Fensal, some of the lakes and dry lake beds south of Ligeia, and the southern shores of Kraken Mare share the absorption feature at 4.92 \um observed in the spectra of Tui and Hotei. Besides Woytchugga and at Fensal, these 5-$\mu$m-bright deposits are geomorphologically-substantiated evaporites. Thus, the similarity in composition strengthens the hypothesis that Tui and Hotei once contained liquid. Other evaporite deposits, however, do not show the 4.92 \um absorption, notably Muggel Lacus and the shores of Ligeia Mare at the north pole. This difference in composition suggests that there are more than one kind of soluble material in Titan's lakes that can create evaporite and/or that the surface properties at the VIMS wavelength scale are not uniform between the different deposits (crystal size, abundance, etc). Our results indicate that the surface structure, composition, and formation history of Titan's evaporites may be at least as dynamic and complex as their Earth counterparts.

Numerous geological features that could be evaporitic in origin have been identified on the surface of Titan. Although they seem to be water-ice poor, their main properties -chemical composition, thickness, stratification- are essentially unknown. In this paper, which follows on a previous one focusing on the surface composition (Cordier et al., 2013), we provide some answers to these questions derived from a new model. This model, based on the up-to-date thermodynamic theory known as "PC-SAFT", has been validated with available laboratory measurements and specifically developed for our purpose. 1-D models confirm the possibility of an acetylene and/or butane enriched central layer of evaporitic deposit. The estimated thickness of this acetylene-butane layer could explain the strong RADAR brightness of the evaporites. The 2-D computations indicate an accumulation of poorly soluble species at the deposit's margin. Among these species, HCN or aerosols similar to tholins could play a dominant role. Our model predicts the existence of chemically trimodal "bathtub rings" which is consistent with what it is observed at the south polar lake Ontario Lacus. This work also provides plausible explanations to the lack of evaporites in the south polar region and to the high radar reflectivity of dry lakebeds.

We use rotational gravity darkening in the disk of \emph{Kepler} star KOI-2138 to show that the orbit of $2.1-R_\oplus$ transiting planet candidate KOI-2138.01 has a low projected spin-orbit alignment of $\lambda=1^\circ\pm13$. KOI-2138.01 is just the second super-Earth with a measured spin-orbit alignment after 55 Cancri e, and the first to be aligned. With a 23.55-day orbital period, KOI-2138.01 may represent the tip of a future iceberg of solar-system-like terrestrial planets having intermediate periods and low-inclination circular orbits.

We present coordinated multiwavelength observations of the bright, nearby BL Lac object Mrk 421 taken in 2013 January-March, involving GASP-WEBT, Swift, NuSTAR, Fermi-LAT, MAGIC, VERITAS, and other collaborations and instruments, providing data from radio to very-high-energy (VHE) gamma-ray bands. NuSTAR yielded previously unattainable sensitivity in the 3-79 keV range, revealing that the spectrum softens when the source is dimmer until the X-ray spectral shape saturates into a steep power law with a photon index of approximately 3, with no evidence for an exponential cutoff or additional hard components up to about 80 keV. For the first time, we observed both the synchrotron and the inverse-Compton peaks of the spectral energy distribution (SED) simultaneously shifted to frequencies below the typical quiescent state by an order of magnitude. The fractional variability as a function of photon energy shows a double-bump structure which relates to the two bumps of the broadband SED. In each bump, the variability increases with energy which, in the framework of the synchrotron self-Compton model, implies that the electrons with higher energies are more variable. The measured multi-band variability, the significant X-ray-to-VHE correlation down to some of the lowest fluxes ever observed in both bands, the lack of correlation between optical/UV and X-ray flux, the low degree of polarization and its significant (random) variations, the short estimated electron cooling time, and the significantly longer variability timescale observed in the NuSTAR light curves point toward in-situ electron acceleration, and suggest that there are multiple compact regions contributing to the broadband emission of Mrk 421 during low-activity states.

We present the first Doppler images of the active eclipsing binary system SZ Psc, based on the high-resolution spectral data sets obtained in 2004 November and 2006 September--December. The least-squares deconvolution technique was applied to derive high signal-to-noise profiles from the observed spectra of SZ Psc. Absorption features contributed by a third component of the system were detected in the LSD profiles at all observed phases. We estimated the mass and period of the third component to be about $0.9 M_{\odot}$ and $1283 \pm 10$ d, respectively. After removing the contribution of the third body from the LSD profiles, we derived the surface maps of SZ Psc. The resulting Doppler images indicate significant starspot activities on the surface of the K subgiant component. The distributions of starspots are more complex than that revealed by previous photometric studies. The cooler K component exhibited pronounced high-latitude spots as well as numerous low- and intermediate-latitude spot groups during the entire observing seasons, but did not show any large, stable polar cap, different from many other active RS CVn-type binaries.

We constrain the true spin-orbit alignment of the KOI-89 system by numerically fitting the two \emph{Kepler} photometric lightcurves produced by transiting planets KOI-89.01 and KOI-89.02. The two planets have periods of 84.69 days and 207.58 days, respectively. We find that the two bodies are low-density giant planets with radii $0.45 \pm 0.03~\mathrm{R_{jup}}$ and $0.43 \pm 0.05~\mathrm{R_{jup}}$ and spin-orbit misalignments $72^{\circ} \pm 3^\circ$ and $73^{\circ+11}_{-5}$, respectively. Via dynamic stability tests we demonstrate the general trend of higher system stability with the two planets close to mutual alignment and estimate their coalignment angle to $20^\circ \pm 20^\circ$ -- i.e. the planets are misaligned with the star but may be aligned with each other. From these results, we limit KOI-89's misalignment mechanisms to star-disk-binary interactions, disk warping via planet-disk interactions, planet-planet scattering, Kozai resonance, or internal gravity waves.

Tidal interactions between disc galaxies depend on galaxy structure, but the details of this relationship are incompletely understood. I have constructed a three-parameter grid of bulge/disc/halo models broadly consistent with $\Lambda$CDM, and simulated an extensive series of encounters using these models. Halo mass and extent strongly influence the dynamics of orbit evolution. In close encounters, the transfer of angular momentum mediated by the dynamical response of massive, extended haloes can reverse the direction of orbital motion of the central galaxies after their first passage. Tidal response is strongly correlated with the ratio $v_\mathrm{e} / v_\mathrm{c}$ of escape to circular velocity within the participating discs. Moreover, the same ratio also correlates with the rate at which tidal tails are reaccreted by their galaxies of origin; consequently, merger remnants with `twin tails', such as NGC 7252, may prove hard to reproduce unless $(v_\mathrm{e} / v_\mathrm{c})^2 \lesssim 5.5$. The tidal morphology of an interacting system can provide useful constraints on progenitor structure. In particular, encounters in which halo dynamics reverses orbital motion exhibit a distinctive morphology which may be recognized observationally. Detailed models attempting to reproduce observations of interacting galaxies should explore the likely range of progenitor structures along with other encounter parameters.

We develop and test an automated technique to model the dynamics of interacting galaxy pairs. We use Identikit (Barnes & Hibbard 2009, Barnes 2011) as a tool for modeling and matching the morphology and kinematics of the interacting pairs of equal-mass galaxies. In order to reduce the effect of subjective human judgement, we automate the selection of phase-space regions used to match simulations to data, and we explore how selection of these regions affects the random uncertainties of parameters in the best-fit model. In this work, we use an independent set of GADGET SPH simulations as input data to determine the systematic bias in the measured encounter parameters based on the known initial conditions of these simulations. We test both cold gas and young stellar components in the GADGET simulations to explore the effect of choosing HI vs. H$\alpha$ as the line of sight velocity tracer. We find that we can group the results into tests with good, fair, and poor convergence based on the distribution of parameters of models close to the best-fit model. For tests with good and fair convergence, we rule out large fractions of parameter space and recover merger stage, eccentricity, pericentric distance, viewing angle, and initial disc orientations within 3$\sigma$ of the correct value. All of tests on prograde-prograde systems have either good or fair convergence. The results of tests on edge-on discs are less biased than face-on tests. Retrograde and polar systems do not converge and may require constraints from regions other than the tidal tails and bridges.

(shortened for arXiv) We aim to progress towards more efficient exoplanet detection around active stars by optimizing the use of Doppler Imaging in radial velocity measurements. We propose a simple method to simultaneously extract a brightness map and a set of orbital parameters through a tomographic inversion technique derived from classical Doppler mapping. Based on the maximum entropy principle, the underlying idea is to determine the set of orbital parameters that minimizes the information content of the resulting Doppler map. We carry out a set of numerical simulations to perform a preliminary assessment of the robustness of our method, using an actual Doppler map of the very active star HR 1099 to produce a realistic synthetic data set for various sets of orbital parameters of a single planet in a circular orbit. Using a simulated time-series of 50 line profiles affected by a peak-to-peak activity jitter of 2.5 km/s, we are able in most cases to recover the radial velocity amplitude, orbital phase and orbital period of an artificial planet down to a radial velocity semi-amplitude of the order of the radial velocity scatter due to the photon noise alone (about 50 m/s in our case). One noticeable exception occurs when the planetary orbit is close to co-rotation, in which case significant biases are observed in the reconstructed radial velocity amplitude, while the orbital period and phase remain robustly recovered. The present method constitutes a very simple way to extract orbital parameters from heavily distorted line profiles of active stars, when more classical radial velocity detection methods generally fail. It is easily adaptable to most existing Doppler Imaging codes, paving the way towards a systematic search for close-in planets orbiting young, rapidly-rotating stars.

The spatial distribution of the tropospheric methane on Titan was measured using near-infrared spectroscopy. Ground-based observations at 1.5$\mu{\rm m}$ (H-band) were performed during the same night using instruments with adaptive optics at both the W. M. Keck Observatory and at the Paranal Observatory on 17 July 2014 UT. The integral field observations with SINFONI on the VLT covered the entire H-band at moderate resolving power, $R=\lambda/\Delta\lambda\approx1,500$, while the Keck observations were performed with NIRSPAO near 1.55254$\mu{\rm m}$ at higher resolution, $R\approx25,000$. The moderate resolution observations are used for flux calibration and for the determination of model parameters that can be degenerate in the interpretation of high resolution spectra. Line-by-line calculations of CH$_4$ and CH$_3$D correlated $k$ distributions from the HITRAN 2012 database were used, which incorporate revised line assignments near 1.5$\mu{\rm m}$. We fit the surface albedo and aerosol distributions in the VLT SINFONI observations that cover the entire H-band window and used these quantities to constrain the models of the high-resolution Keck NIRSPAO spectra when retrieving the methane abundances. Cassini VIMS images of the polar regions, acquired on 20 July 2014 UT, are used to validate the assumption that the opacity of tropospheric aerosol is relatively uniform below 10 km. We retrieved methane abundances at latitudes between 42$^{\circ}$ S and 80$^{\circ}$ N. The tropospheric methane in the Southern mid-latitudes was enhanced by a factor of $\sim$10--40% over the nominal profile that was measured using the GCMS on Huygens. The Northern hemisphere had $\sim$90% of the nominal methane abundance up to polar latitudes (80$^{\circ}$N). These measurements suggest that a source of saturated polar air is equilibrating with dryer conditions at lower latitudes.

Since M4.5 - M9 dwarfs exhibit equatorial rotation velocities of order 10 km/s on average, radial velocity surveys targeting this stellar population will likely need to find methods to effectively remove starspot jitter. We present the first high resolution Doppler images of the M4.5 dwarf, GJ 791.2A, and the M9 dwarf, LP 944-20. The time series spectra of both objects reveal numerous line profile distortions over the rotation period of each star which we interpret as starspots. The transient distortions are modelled with spot/photosphere contrast ratios that correspond to model atmosphere temperature differences of Tphot-Tspot = 200 and 300 K. GJ 791.2A is a fully convective star with vsini = 35.1 km/s. Although we find more starspot structure at high latitudes, we reconstruct spots at a range of phases and latitudes with a mean spot filling of ~3%. LP 944-20 is one of the brightest known late-M dwarfs, with spectral type M9V and vsini = 30.8 km/s. Its spectral time series exhibits two dominant transient line distortions that are reconstructed as high latitude spots, while a mean spot filling factor of only 1.5% is found. The occurrence of low-contrast spots at predominantly high latitudes, which we see in both targets here, is in general likely to be responsible for the low amplitude photometric variability seen in late-M dwarfs. For GJ 791.2A, the radial velocities induced by the starspot features yield an rms velocity variability of 138 m/s, which can be reduced by a factor of 1.9 using our reconstructed surface brightness distributions.

We present Spitzer 4.5\micron\ light curve observations, Keck NIRSPEC radial velocity observations, and LCOGT optical light curve observations of PTFO~8-8695, which may host a Jupiter-sized planet in a very short orbital period (0.45 days). Previous work by \citet{vaneyken12} and \citet{barnes13} predicts that the stellar rotation axis and the planetary orbital plane should precess with a period of $300 - 600$ days. As a consequence, the observed transits should change shape and depth, disappear, and reappear with the precession. Our observations indicate the long-term presence of the transit events ($>3$ years), and that the transits indeed do change depth, disappear and reappear. The Spitzer observations and the NIRSPEC radial velocity observations (with contemporaneous LCOGT optical light curve data) are consistent with the predicted transit times and depths for the $M_\star = 0.34\ M_\odot$ precession model and demonstrate the disappearance of the transits. An LCOGT optical light curve shows that the transits do reappear approximately 1 year later. The observed transits occur at the times predicted by a straight-forward propagation of the transit ephemeris. The precession model correctly predicts the depth and time of the Spitzer transit and the lack of a transit at the time of the NIRSPEC radial velocity observations. However, the precession model predicts the return of the transits approximately 1 month later than observed by LCOGT. Overall, the data are suggestive that the planetary interpretation of the observed transit events may indeed be correct, but the precession model and data are currently insufficient to confirm firmly the planetary status of PTFO~8-8695b.

Exoplanet science is now in its full expansion, particularly after the CoRoT and Kepler space missions that led us to the discovery of thousands of extra-solar planets. The last decade has taught us that UV observations play a major role in advancing our understanding of planets and of their host stars, but the necessary UV observations can be carried out only by HST, and this is going to be the case for many years to come. It is therefore crucial to build a treasury data archive of UV exoplanet observations formed by a dozen "golden systems" for which observations will be available from the UV to the infrared. Only in this way we will be able to fully exploit JWST observations for exoplanet science, one of the key JWST science case.

This paper presents the simultaneous Bayesian analysis of the radial velocity and activity indices for the K5Ve active star BD+20 1790, that was proposed to host a close-in massive planet by Hern\'an-Obispo et al. (2010). The Bayesian analysis supports the hypothesis of a planetary companion. We present a new orbital solution after removing the two main contributions to the stellar jitter, one that varies with the photometric period and another that varies with the synodic period of the star+planet system. We use a new method to determine these jitter components, considering them as a second and third signals in the system. A discussion on possible magnetic star-planet-interaction is included, based on the Bayesian analysis of the activity indices. We propose two possible sources for flare events of this star: one related to the geometry of the system and the relative movement of the star+planet, and a second one purely stochastic source that is related to the evolution of stellar active regions. Also, we observe for the first time the magnetic field of the star, from spectropolarimetric data.

We report ALMA Band 7 (350 GHz) imaging at 0.4 - 0.6arcsec resolution and Band 9 (696 GHz) at ~0.25arcsec resolution of the luminous IR galaxies Arp 220 and NGC 6240. The long wavelength dust continuum is used to estimate ISM masses for Arp 220 East, West and NGC 6240 of 1.9, 4.2 and 1.6x10^9 msun within radii of 69, 65 and 190 pc. The HCN emission was modeled to derive the emissivity distribution as a function of radius and the kinematics of each nuclear disk, yielding dynamical masses consistent with the masses and sizes derived from the dust emission. In Arp 220, the major dust and gas concentrations are at radii less than 50 pc in both counter-rotating nuclear disks. The thickness of the disks in Arp 220estimated from the velocity dispersion and rotation velocities are 10-20 pc and the mean gas densities are n_H2 ~10^5 cm^-3 at R < 50 pc. We develop an analytic treatment for the molecular excitation (including photon trapping), yielding volume densities for both the HCN and CS emission with n_H2 ~2x10^5 cm^-3. The agreement of the mean density from the total mass and size with that required for excitation suggests that the volume is essentially filled with dense gas, i.e. it is not cloudy or like swiss cheese.

Using three dimensional Monte Carlo radiation transfer models of photoionisation and dust scattering, we explore different components of the widespread diffuse H$\alpha$ emission observed in the interstellar medium of the Milky Way and other galaxies. We investigate the relative contributions of H$\alpha$ from recombination emission in ionised gas and H$\alpha$ that originates in H{\sc ii} regions near the Galactic midplane and scatters off high altitude dust in the diffuse interstellar medium. For the radiation transfer simulations we consider two geometries for the interstellar medium: a three dimensional fractal geometry that reproduces the average density structure inferred for hydrogen in the Milky Way, and a density structure from a magneto hydrodynamic simulation of a supernovae driven turbulent interstellar medium. Although some sight lines that are close to H{\sc ii} regions can be dominated by scattered light, overall we find that less than $\sim 20\%$ of the total H$\alpha$ intensity in our simulations can be attributed to dust scattering. Our findings on the relative contribution of scattered H$\alpha$ are consistent with previous observational and theoretical analyses. We also investigate the relative contributions of dust scattering and in situ ionisation of high density dust clouds in the diffuse gas. Dust scattering in these partially ionised clouds contribute $\sim 40\%$ to the total intensity of H$\alpha$.

We present Doppler images of both components of the eclipsing binary system ER Vul, based on the spectra obtained in 2004 November, 2006 September and 2008 November. The least-squares deconvolution technique is used for enhancing the signal-to-noise ratios of the observed profiles. The new surface images reveal that both stars of ER Vul show strong starspot activities and the starspots appear at various latitudes. The surface maps of 2006 and 2008 both show the presence of large high-latitude starspots on each component of ER Vul. We find no obvious phase shift of the active regions during our observations. The longitude distributions of starspots are non-uniform on both stars. At low-to-mid latitudes, the active regions are almost exclusively found in the hemisphere facing the other star. However, we find no pronounced concentration of spots at the sub-stellar points.

We report the discovery of WTS-1b, the first extrasolar planet found by the WFCAM Transit Survey, which began observations at the 3.8-m United Kingdom Infrared Telescope (UKIRT) in August 2007. Light curves comprising almost 1200 epochs with a photometric precision of better than 1 per cent to J ~ 16 were constructed for ~60000 stars and searched for periodic transit signals. For one of the most promising transiting candidates, high-resolution spectra taken at the Hobby-Eberly Telescope (HET) allowed us to estimate the spectroscopic parameters of the host star, a late-F main sequence dwarf (V=16.13) with possibly slightly subsolar metallicity, and to measure its radial velocity variations. The combined analysis of the light curves and spectroscopic data resulted in an orbital period of the substellar companion of 3.35 days, a planetary mass of 4.01 +- 0.35 Mj and a planetary radius of 1.49+0.16-0.18 Rj. WTS-1b has one of the largest radius anomalies among the known hot Jupiters in the mass range 3-5 Mj. The high irradiation from the host star ranks the planet in the pM class.

Water-ice-poor, 5-$\mu$m-bright material on Saturn's moon Titan has previously been geomorphologically identified as evaporitic. Here we present a global distribution of the occurrences of the 5-$\mu$m-bright spectral unit, identified with Cassini's Visual Infrared Mapping Spectrometer (VIMS) and examined with RADAR when possible. We explore the possibility that each of these occurrences are evaporite deposits. The 5-$\mu$m-bright material covers 1\% of Titan's surface and is not limited to the poles (the only regions with extensive, long-lived surface liquid). We find the greatest areal concentration to be in the equatorial basins Tui Regio and Hotei Regio. Our interpretations, based on the correlation between 5-$\mu$m-bright material and lakebeds, imply that there was enough liquid present at some time to create the observed 5-$\mu$m-bright material. We address the climate implications surrounding a lack of evaporitic material at the south polar basins: if the south pole basins were filled at some point in the past, then where is the evaporite?

Ultraluminous and luminous infrared galaxies (ULIRGs and LIRGs) are the most extreme star-forming galaxies in the universe, and dominate the total star formation rate density at z>1. In the local universe (z<0.3), the majority of ULIRGs and a significant portion of LIRGs are triggered by interactions between gas-rich spiral galaxies, yet it is unclear if this is still the case at high-z. To investigate the relative importance of galaxy interactions in infrared luminous galaxies, we carry out a comparison of optical morphological properties between local (U)LIRGs and (U)LIRGs at z=0.5-1.5 based on the same sample selection, morphology classification scheme, and optical morphology at similar rest-frame wavelengths. In addition, we quantify the systematics in comparing local and high-z datasets by constructing a redshifted dataset from local (U)LIRGs, in which its data quality mimics the high-z dataset. Based on the Gini-M20 classification scheme, we find that the fraction of interacting systems decreases by ~8% from local to z<~1, and it is consistent with the reduction between local and redshifted datasets (6(+14-6)%). Based on visual classifications, the merger fraction of local ULIRGs is found to be ~20% lower compared to published results, and the reduction due to redshifiting is 15(+10-8)%. Consequently, the differences of merger fractions between local and z<~1 (U)LIRGs is only ~17%. These results demonstrate that there is no strong evolution in the fraction of (U)LIRGs classified as mergers at least out to z~1. At z>1, the morphology types of ~30% of (U)LIRGs can not be determined due to their faintness in the F814W-band, and thus the merger fraction measured at z>1 suffers from large uncertainties.

We fit the \emph{Kepler} photometric light curve of the KOI-368 system using an oblate, gravity-darkened stellar model in order to constrain its spin-orbit alignment. We find that the system is relatively well-aligned with a sky-projected spin-orbit alignment of $\lambda=10^{\circ}\pm2^{\circ}$, a stellar obliquity of $\psi=3^{\circ}\pm7^{\circ}$, and a true spin-orbit alignment of $\varphi=11^{\circ}\pm3^{\circ}$. Although our measurement differs significantly from zero, the low value for $\varphi$ is consistent with spin-orbit alignment. We also measure various transit parameters of the KOI-368 system: $R_{\star}=2.28\pm0.02 R_{\odot}$, $R_{p}=1.83\pm0.02R_{jup}$, and $i=89.221^{\circ}\pm0.013^{\circ}$. This work shows that our gravity-darkened model can constrain long-period, well-aligned planets and M-class stars orbiting fast-rotators, allowing for measurement of a new subcategory of transiting bodies.

The Diffuse Ionised Gas (DIG) in galaxies traces photoionisation feedback from massive stars. Through three dimensional photoionisation simulations, we study the propagation of ionising photons, photoionisation heating and the resulting distribution of ionised and neutral gas within snapshots of magnetohydrodynamic simulations of a supernova driven turbulent interstellar medium. We also investigate the impact of non-photoionisation heating on observed optical emission line ratios. Inclusion of a heating term which scales less steeply with electron density than photoionisation is required to produce diagnostic emission line ratios similar to those observed with the Wisconsin H{\alpha} Mapper. Once such heating terms have been included, we are also able to produce temperatures similar to those inferred from observations of the DIG, with temperatures increasing to above 15000 K at heights |z| > 1 kpc. We find that ionising photons travel through low density regions close to the midplane of the simulations, while travelling through diffuse low density regions at large heights. The majority of photons travel small distances (< 100pc); however some travel kiloparsecs and ionise the DIG.

Due to their higher planet-star mass-ratios, M dwarfs are the easiest targets for detection of low-mass planets orbiting nearby stars using Doppler spectroscopy. Furthermore, because of their low masses and luminosities, Doppler measurements enable the detection of low-mass planets in their habitable zones that correspond to closer orbits than for Solar-type stars. We re-analyse literature UVES radial velocities of 41 nearby M dwarfs in a combination with new velocities obtained from publicly available spectra from the HARPS-ESO spectrograph of these stars in an attempt to constrain any low-amplitude Keplerian signals. We apply Bayesian signal detection criteria, together with posterior sampling techniques, in combination with noise models that take into account correlations in the data and obtain estimates for the number of planet candidates in the sample. More generally, we use the estimated detection probability function to calculate the occurrence rate of low-mass planets around nearby M dwarfs. We report eight new planet candidates in the sample (orbiting GJ 27.1, GJ 160.2, GJ 180, GJ 229, GJ 422, and GJ 682), including two new multiplanet systems, and confirm two previously known candidates in the GJ 433 system based on detections of Keplerian signals in the combined UVES and HARPS radial velocity data that cannot be explained by periodic and/or quasiperiodic phenomena related to stellar activities. Finally, we use the estimated detection probability function to calculate the occurrence rate of low-mass planets around nearby M dwarfs. According to our results, M dwarfs are hosts to an abundance of low-mass planets and the occurrence rate of planets less massive than 10 M$_{\oplus}$ is of the order of one planet per star, possibly even greater. ...

We report the discovery of WTS-2 b, an unusually close-in 1.02-day hot Jupiter (Mp=1.12MJ, Rp=1.363RJ) orbiting a K2V star, which has a possible gravitationally-bound M-dwarf companion at 0.6 arcsec separation contributing ~20 percent of the total flux in the observed J-band light curve. The planet is only 1.5 times the separation from its host star at which it would be destroyed by Roche lobe overflow, and has a predicted remaining lifetime of just ~40 Myr, assuming a tidal dissipation quality factor of Q'*=10^6. Q'* is a key factor in determining how frictional processes within a host star affect the orbital evolution of its companion giant planets, but it is currently poorly constrained by observations. We calculate that the orbital decay of WTS-2 b would correspond to a shift in its transit arrival time of T_shift~17 seconds after 15 years assuming Q'*=10^6. A shift less than this would place a direct observational constraint on the lower limit of Q'* in this system. We also report a correction to the previously published expected T_shift for WASP-18 b, finding that T_shift=356 seconds after 10 years for Q'*=10^6, which is much larger than the estimated 28 seconds quoted in WASP-18 b discovery paper. We attempted to constrain Q'* via a study of the entire population of known transiting hot Jupiters, but our results were inconclusive, requiring a more detailed treatment of transit survey sensitivities at long periods. We conclude that the most informative and straight-forward constraints on Q'* will be obtained by direct observational measurements of the shift in transit arrival times in individual hot Jupiter systems. We show that this is achievable across the mass spectrum of exoplanet host stars within a decade, and will directly probe the effects of stellar interior structure on tidal dissipation.

We present a study of age-related spectral signatures observed in 25 young low-mass objects that we have previously determined as possible kinematic members of five young moving groups: the Local Association (Pleiades moving group, age=20 - 150 Myr), the Ursa Major group (Sirius supercluster, age=300 Myr), the Hyades supercluster (age=600 Myr), IC 2391 supercluster (age=35--55 Myr) and the Castor moving group (age=200 Myr). In this paper we characterize the spectral properties of observed high or low resolution spectra of our kinematic members by fitting theoretical spectral distributions. We study signatures of youth, such as lithium {\sc i} 6708 \AA, H$\alpha$ emission and other age-sensitive spectroscopic signatures in order to confirm the kinematic memberships through age constraints. We find that 21 ($84\%$) targets show spectroscopic signatures of youth in agreement with the age ranges of the moving group to which membership is implied. For two further objects, age-related constraints remain difficult to determine from our analysis. In addition, we confirm two moving group kinematic candidates as brown dwarfs.

We present radial velocity measurements of a sample of M5V-M9V stars from our Red-Optical Planet Survey, ROPS, operating at 0.65-1.025 micron. Radial velocities for 15 stars, with r.m.s. precision down to 2.5 m/s over a week long time scale are achieved using Thorium-Argon reference spectra. We are sensitive to planets with m_psin(i) >= 1.5 MEarth (3 MEarth at 2-sigma) in the classical habitable zone and our observations currently rule out planets with m_psin(i) > 0.5 MJup at 0.03 AU for all our targets. A total of 9 of the 15 targets exhibit r.m.s. < 16 m/s, which enables us to rule out the presence of planets with m_psin(i) > 10 MEarth in 0.03 AU orbits. Since the mean rotation velocity is of order 8 km/s for an M6V star and 15 km/s by M9V, we avoid observing only slow rotators that would introduce a bias towards low axial inclination i << 90 deg systems, which are unfavourable for planet detection. Our targets with the highest vsini values exhibit radial velocities significantly above the photon-noise limited precision, even after accounting for vsini. We therefore monitored stellar activity via chromospheric emission from the Halpha and Ca II infrared triplet lines. A clear trend of log_10(Halpha/Lbol) with radial velocity r.m.s. is seen, implying that significant starspot activity is responsible for the observed radial velocity precision floor. The implication that most late M dwarfs are significantly spotted, and hence exhibit time varying line distortions, indicates that observations to detect orbiting planets need strategies to reliably mitigate against the effects of activity induced radial velocity variations.

Theoretical predictions suggest that the distribution of planets in very young stars could be very different to that typically observed in Gyr old systems that are the current focus of radial velocity surveys. However, the detection of planets around young stars is hampered by the increased stellar activity associated with young stars, the signatures of which can bias the detection of planets. In this paper we place realistic limitations on the possibilities for detecting planets around young active G and K dwarfs. The models of stellar activity based on tomographic imaging of the G dwarf HD 141943 and the K1 dwarf AB Dor and also include contributions from plage and many small random starspots. Our results show that the increased stellar activity levels present on young Solar-type stars strongly impacts the detection of Earth-mass and Jupiter mass planets and that the degree of activity jitter is directly correlated with stellar \vsinis. We also show that for G and K dwarfs, the distribution of activity in individual stars is more important than the differences in induced radial velocities as a function of spectral type. We conclude that Jupiter mass planets can be detected close-in around fast-rotating young active stars, Neptune-mass planets around moderate rotators and that Super-Earths are only detectable around very slowly rotating stars. The effects of an increase in stellar activity jitter by observing younger stars can be compensated for by extending the observational base-line to at least 100 epochs.

Galaxy interactions/mergers have been shown to dominate the population of IR luminous galaxies (log(LIR)>11.6Lsun) in the local Universe (z<0.25). Recent studies based on the relation between galaxies' star formation rates and stellar mass (the SFR-M relation or the galaxy main sequence (MS)) have suggested that galaxy interaction/mergers may only become significant when galaxies fall well above the galaxy MS. Since the typical SFR at given M increases with redshift, the existence of galaxy MS implies that massive, IR-luminous galaxies at high-z may not necessarily be driven by galaxy interactions. We examine the role of galaxy interactions in the SFR-M relation by carrying out a morphological analysis of 2084 Herschel-selected galaxies at 0.2 < z < 1.5 in the COSMOS field. Herschel-PACS and -SPIRE observations covering the full 2-deg^2 COSMOS field provide one of the largest far-IR selected samples of high-redshift galaxies with well-determined redshifts to date, with sufficient sensitivity at z ~ 1, to sample objects lying on and above the galaxy MS. Using a detailed visual classification scheme, we show that the fraction of "disk galaxies" decreases and the fraction of "irregular" galaxies increases systematically with increasing LIR out to z ~ 1.5 and z ~ 1.0, respectively. At log(LIR) > 11.5 Lsun, >50% of the objects show evident features of strongly interacting/merger systems, where this percentage is similar to the studies of local IR-luminous galaxies. The fraction of interacting/merger systems also systematically increases with the deviation from the SFR-M relation, supporting the view that galaxies fall above the MS are more dominated by mergers than the MS galaxies. Meanwhile, we find that ~18% of massive IR-luminous MS galaxies are classified as interacting systems, where this population may not evolve through the evolutionary track predicted by a simple gas exhaustion model.

Titan, the main satellite of Saturn, has an active cycle of methane in its troposphere. Among other evidence for a mechanism of evaporation at work on the ground, dry lakebeds have been discovered. Recent Cassini infrared observations of these empty lakes have revealed a surface composition poor in water ice compared to that of the surrounding terrains --- suggesting the existence of organic evaporites deposits. The chemical composition of these possible evaporites is unknown. In this paper, we study evaporite composition using a model that treats both organic solids dissolution and solvent evaporation. Our results suggest the possibility of large abundances of butane and acetylene in the lake evaporites. However, due to uncertainties of the employed theory, these determinations have to be confirmed by laboratory experiments.

PTFO 8-8695b represents the first transiting exoplanet candidate orbiting a pre-main-sequence star. We find that the unusual lightcurve shapes of PTFO 8-8695 can be explained by transits of a planet across an oblate, gravity-darkened stellar disk. We simultaneously and self-consistently fit two separate lightcurves observed in 2009 December and 2010 December. Our two self-consistent fits yield M_p = 3.0 M_Jup and M_p = 3.6 M_Jup for assumed stellar masses of M_* = 0.34 M_Sun and M_* = 0.44 M_Sun respectively. The two fits have precession periods of 293 days and 581 days. These mass determinations (consistent with previous upper limits) along with the strength of the gravity-darkened precessing model together validate PTFO 8-8695b as just the second Hot Jupiter known to orbit an M-dwarf. Our fits show a high degree of spin-orbit misalignment in the PTFO 8-8695 system: 69 +/- 2 or 73.1 +/- 0.5 degrees, in the two cases. The large misalignment is consistent with the hypothesis that planets become Hot Jupiters with random orbital plane alignments early in a system's lifetime. We predict that as a result of the highly misaligned, precessing system, the transits should disappear for months at a time over the course of the system's precession period. The precessing, gravity-darkened model also predicts other observable effects: changing orbit inclination that could be detected by radial velocity observations, changing stellar inclination that would manifest as varying v sin i, changing projected spin-orbit alignment that could be seen by the Rossiter-McLaughlin effect, changing transit shapes over the course of the precession, and differing lightcurves as a function of wavelength. Our measured planet radii of 1.64 R_Jup and 1.68 R_Jup in each case are consistent with a young, hydrogen-dominated planet that results from a hot-start formation mechanism.

Using the Position and Proper Motion Extended-L (PPMXL) catalogue, we have used optical and near-infrared colour cuts together with a reduced proper motion cut to find bright M dwarfs for future exoplanet transit studies. PPMXL's low proper motion uncertainties allow us to probe down to smaller proper motions than previous similar studies. We have combined unique objects found with this method to that of previous work to produce 8479 K<9 M dwarfs. Low resolution spectroscopy was obtained of a sample of the objects found using this selection method to gain statistics on their spectral type and physical properties. Results show a spectral type range of K7-M4V. This catalogue is the most complete collection of K<9 M dwarfs currently available and is made available here.

We present dynamical models of four interacting systems: NGC 5257/8, The Mice, the Antennae, and NGC 2623. The parameter space of the encounters are constrained using the Identikit model-matching and visualization tool. Identikit utilizes hybrid N-body and test particle simulations to enable rapid exploration of the parameter space of galaxy mergers. The Identikit-derived matches of these systems are reproduced with self-consistent collisionless simulations which show very similar results. The models generally reproduce the observed morphology and \HI\ kinematics of the tidal tails in these systems with reasonable properties inferred for the progenitor galaxies. The models presented here are the first to appear in the literature for NGC 5257/8 and NGC 2623, and The Mice and the Antennae are compared with previously published models. Based on the assumed mass model and our derived initial conditions, the models indicate the four systems are currently being viewed 175-260 Myr after first passage and cover a wide range of merger stages. In some instances there are mismatches between the models and the data (e.g., in the length of a tail); these are likely due to our adoption of a single mass model for all galaxies. Despite the use of a single mass model, these results demonstrate the utility of Identikit in constraining the parameter space for galaxy mergers when applied to real data.

We compare emission line intensities from photoionization models of smooth and fractal shell geometries for low density H II regions, with particular focus on the low-ionization diagnostic diagram [N II]/H-alpha vs H-alpha. Building on previously published models and observations of Barnard's Loop, we show that the observed range of intensities and variations in the line intensity ratios may be reproduced with a three dimensional shell geometry. Our models adopt solar abundances throughout the model nebula, in contrast with previous one dimensional modeling which suggested the variations in line intensity ratios could only be reproduced if the heavy element abundances were increased by a factor of 1.4. For spatially resolved H II regions, the multiple sightlines that pierce and sample different ionization and temperature conditions within smooth and fractal shells produce a range of line intensities that are easily overlooked if only the total integrated intensities from the entire nebula model are computed. Our conclusion is that inference of H II region properties, such as elemental abundances, via photoionization models of one dimensional geometries must be treated with caution and further tested through three dimensional modeling.

The coalescence of compact objects are a promising astrophysical sources of gravitational wave (GW) signals. The ejection of r-process material from such mergers may lead to a radioactively-powered electromagnetic counterpart which, if discovered, would enhance the science return of a GW detection. As very little is known about the optical properties of heavy r-process elements, previous light curve models have adopted opacities similar to those of iron group elements. Here we report that the presence of heavier elements, particularly the lanthanides, increase the ejecta opacity by several orders of magnitude. We include these higher opacities in time dependent, multi-wavelength radiative transport calculations to predict the broadband light curves of one-dimensional models over a range of parameters (ejecta masses from 0.001 to 0.1 solar masses and velocities from 0.1 to 0.3c). We find that the higher opacities lead to much longer duration light curves which can last a week or more. The emission is shifted toward the infrared bands due to strong optical line blanketing, and the colors at later times are representative of a blackbody near the recombination temperature of the lanthanides (T ~ 2500K). We further consider the case in which a second mass outflow, composed of 56Ni, is ejected from a disk wind, and show that the net result is a distinctive two component spectral energy distribution, with a bright optical peak due to 56Ni and an infrared peak due to r-process ejecta. We briefly consider the prospects for detection and identification of these transients.

Material ejected during (or immediately following) the merger of two neutron stars may assemble into heavy elements by the r-process. The subsequent radioactive decay of the nuclei can power electromagnetic emission similar to, but significantly dimmer than, an ordinary supernova. Identifying such events is an important goal of future transient surveys, offering new perspectives on the origin of r-process nuclei and the astrophysical sources of gravitational waves. Predictions of the transient light curves and spectra, however, have suffered from the uncertain optical properties of heavy ions. Here we consider the opacity of expanding r-process material and argue that it is dominated by line transitions from those ions with the most complex valence electron structure, namely the lanthanides. For a few representative ions, we run atomic structure models to calculate radiative data for tens of millions of lines. We find that the resulting r-process opacities are orders of magnitude larger than that of ordinary (e.g., iron-rich) supernova ejecta. Radiative transport calculations using these new opacities indicate that the transient emission should be dimmer and redder than previously thought. The spectra appear pseudo-blackbody, with broad absorption features, and peak in the infrared (~1 micron). We discuss uncertainties in the opacities and attempt to quantify their impact on the spectral predictions. The results have important implications for observational strategies to find and study the radioactively powered electromagnetic counterparts to compact object mergers.

With time-series spectroscopic observations taken with the Near Infrared Spectrometer (NIRSPEC) at Keck II, we investigated the atmosphere of the close orbiting transiting extrasolar giant planet, HD 189733b. In particular, we intended to measure the dense absorption line forest around 2.3 micron, which is produced by carbon monoxide (CO). CO is expected to be present in the planetary atmosphere, although no detection of this molecule has been claimed yet. To identify the best suited data analysis method, we created artificial spectra of planetary atmospheres and analyzed them by three approaches found in the literature, the deconvolution method, data modeling via chi^2-minimization, and cross-correlation. As a result, we found that cross-correlation and chi^2-data modeling show systematically a higher sensitivity than the deconvolution method. We analyzed the NIRSPEC data with cross-correlation and detect CO absorption in the day-side spectrum of HD 189733b at the known planetary radial velocity semi-amplitude with 3.4 sigma confidence.

In these proceedings we give a status update of the Calan-Hertfordshire Extrasolar Planet Search, an international collaboration led from Chile that aims to discover more planets around super metal-rich and Sun-like stars, and then follow these up with precision photometry to hunt for new bright transit planets. We highlight some results from this program, including exoplanet and brown dwarf discoveries, and a possible correlation between metallicity and planetary minimum mass at the lowest planetary masses detectable. Finally we discuss the short-term and long-term future pathways this program can take.

The abilities of radial velocity exoplanet surveys to detect the lowest-mass extra-solar planets are currently limited by a combination of instrument precision, lack of data, and "jitter". Jitter is a general term for any unknown features in the noise, and reflects a lack of detailed knowledge of stellar physics (asteroseismology, starspots, magnetic cycles, granulation, and other stellar surface phenomena), as well as the possible underestimation of instrument noise. We study an extensive set of radial velocities for the star HD 10700 ($\tau$ Ceti) to determine the properties of the jitter arising from stellar surface inhomogeneities, activity, and telescope-instrument systems, and perform a comprehensive search for planetary signals in the radial velocities. We perform Bayesian comparisons of statistical models describing the radial velocity data to quantify the number of significant signals and the magnitude and properties of the excess noise in the data. We reach our goal by adding artificial signals to the "flat" radial velocity data of HD 10700 and by seeing which one of our statistical noise models receives the greatest posterior probabilities while still being able to extract the artificial signals correctly from the data. We utilise various noise components to assess properties of the noise in the data and analyse the HARPS, AAPS, and HIRES data for HD 10700 to quantify these properties and search for previously unknown low-amplitude Keplerian signals. ...

We announce the discovery of a low-mass planet orbiting the super metal-rich K0V star HD77338 as part of our on-going Calan-Hertfordshire Extrasolar Planet Search. The best fit planet solution has an orbital period of 5.7361\pm0.0015 days and with a radial velocity semi-amplitude of only 5.96\pm1.74 m/s, we find a minimum mass of 15.9+4.7-5.3 Me. The best fit eccentricity from this solution is 0.09+0.25-0.09, and we find agreement for this data set using a Bayesian analysis and a periodogram analysis. We measure a metallicity for the star of +0.35\pm0.06 dex, whereas another recent work (Trevisan et al. 2011) finds +0.47\pm0.05 dex. Thus HD77338b is one of the most metal-rich planet host stars known and the most metal-rich star hosting a sub-Neptune mass planet. We searched for a transit signature of HD77338b but none was detected. We also highlight an emerging trend where metallicity and mass seem to correlate at very low masses, a discovery that would be in agreement with the core accretion model of planet formation. The trend appears to show that for Neptune-mass planets and below, higher masses are preferred when the host star is more metal-rich. Also a lower boundary is apparent in the super metal-rich regime where there are no very low-mass planets yet discovered in comparison to the sub-solar metallicity regime. A Monte Carlo analysis shows that this, low-mass planet desert, is statistically significant with the current sample of 36 planets at around the 4.5\sigma\ level. In addition, results from Kepler strengthen the claim for this paucity of the lowest-mass planets in super metal-rich systems. Finally, this discovery adds to the growing population of low-mass planets around low-mass and metal-rich stars and shows that very low-mass planets can now be discovered with a relatively small number of data points using stable instrumentation.

We present an analysis of the photometric variability of M dwarfs in the WFCAM Transit Survey. Although periodic lightcurve variability in low mass stars is generally dominated by photospheric star spot activity, M dwarf variability in the J band has not been as thoroughly investigated as at visible wavelengths. Spectral type estimates for a sample of over 200,000 objects are made using spectral type-colour relations, and over 9600 dwarfs (J<17) with spectral types later than K7 were found. The light curves of the late-type sample are searched for periodicity using a Lomb-Scargle periodogram analysis. A total of 68 periodic variable M dwarfs are found in the sample with periods ranging from 0.16 days to 90.33 days, with amplitudes in the range of ~0.009 to ~0.115 in the J band. We simulate active M dwarfs with a range of latitude-independent spot coverages and estimate a periodically variable fractions of 1-3 per cent for stars where spots cover more than 10 per cent of the star's surface. Our simulated spot distributions indicate that operating in the J band, where spot contrast ratios are minimised, enables variability in only the most active of stars to be detected. These findings affirm the benefits of using the $J$ band for planetary transit searches compared to visible bands. We also serendipitously find a \Delta J>0.2 mag flaring event from an M4V star in our sample.

We present observations of the full January 2010 transit of HD80606b from the Canadian microsatellite, Microvariability and Oscillations of Stars (MOST). By employing a space-based telescope, we monitor the entire transit thus limiting systematic errors that result from ground observations. We determine measurements for the planetary radius (R_{p}=0.987\pm0.061R_{Jup}) and inclination (i=89.283^{o}\pm0.024) by constraining our fits with the observed parameters of different groups. Our measured mid-transit time of 2455210.6449\pm0.0034 HJD is consistant with the 2010 Spitzer results and is 20 minutes earlier than predicted by groups who observed the June 2009 transit.

Observatories constructed solely for photometric monitoring make it possible to understand the temporal nature of objects over time scales that historically have been difficult to achieve. We report on one such observatory, the Remote Observatory for Variable Object Research (ROVOR) which was constructed to enable both long-term and rapid cadence observations of brighter objects. ROVOR is an 0.4m optical telescope located in central Utah and commissioned for scientific observations in 2008. Principle research has been monitoring blazars, x-ray binaries, AGN, and an occasional gamma-ray burst afterglow. We describe the observatory, the control system, and its unique roof.

We report the discovery of WTS-1b, the first extrasolar planet found by the WFCAM Transit Survey, which began observations at the 3.8-m United Kingdom Infrared Telescope. Light curves comprising almost 1200 epochs with a photometric precision of better than 1 per cent to J=16 were constructed for 60000 stars and searched for periodic transit signals. For one of the most promising transiting candidates, high-resolution spectra taken at the Hobby-Eberly Telescope allowed us to estimate the spectroscopic parameters of the host star, a late-F main sequence dwarf (V=16.13) with possibly slightly subsolar metallicity, and to measure its radial velocity variations. The combined analysis of the light curves and spectroscopic data resulted in an orbital period of the substellar companion of 3.35 days, a planetary mass of 4.01+-0.35 Mj and a planetary radius of 1.49+-0.17 Rj. WTS-1b has one of the largest radius anomalies among the known hot Jupiters in the mass range 3-5 Mj.

Luminous and ultraluminous infrared galaxies ((U)LIRGs) are the most extreme star forming galaxies in the universe. The local (U)LIRGs provide a unique opportunity to study their multi-wavelength properties in detail for comparison to their more numerous counterparts at high redshifts. We present common large aperture photometry at radio through X-ray wavelengths, and spectral energy distributions (SEDs) for a sample of 53 nearby LIRGs and 11 ULIRGs spanning log (LIR/Lsun) = 11.14-12.57 from the flux-limited Great Observatories All-sky LIRG Survey (GOALS). The SEDs for all objects are similar in that they show a broad, thermal stellar peak and a dominant FIR thermal dust peak, where nuLnu(60um) / nuLnu(V) increases from ~2-30 with increasing LIR. When normalized at IRAS-60um, the largest range in the luminosity ratio, R(lambda)=log[nuLnu(lambda)/nuLnu(60um)] observed over the full sample is seen in the Hard X-rays (HX=2-10 keV). A small range is found in the Radio (1.4GHz), where the mean ratio is largest. Total infrared luminosities, LIR(8-1000um), dust temperatures, and dust masses were computed from fitting thermal dust emission modified blackbodies to the mid-infrared (MIR) through submillimeter SEDs. The new results reflect an overall ~0.02 dex lower luminosity than the original IRAS values. Total stellar masses were computed by fitting stellar population synthesis models to the observed near-infrared (NIR) through ultraviolet (UV) SEDs. Mean stellar masses are found to be log(M/Msun) = 10.79+/-0.40. Star formation rates have been determined from the infrared (SFR_IR~45Msun/yr) and from the monochromatic UV luminosities (SFR_UV~1.3Msun/yr), respectively. Multiwavelength AGN indicators have be used to select putative AGN: about 60% of the ULIRGs would have been classified as an AGN by at least one of the selection criteria.

In self-consistent N-body simulations of collisionless systems, gravitational interactions are modified on small scales to remove singularities and simplify the task of numerically integrating the equations of motion. This `gravitational softening' is sometimes presented as an ad-hoc departure from Newtonian gravity. However, softening can also be described as a smoothing operation applied to the mass distribution; the gravitational potential and the smoothed density obey Poisson's equation precisely. While `softening' and `smoothing' are mathematically equivalent descriptions, the latter has some advantages. For example, the smoothing description suggests a way to set up N-body initial conditions in almost perfect dynamical equilibrium.

We formulated tidal decay lifetimes for hypothetical moons orbiting extrasolar planets with both lunar and stellar tides. Previous work neglected the effect of lunar tides on planet rotation, and are therefore applicable only to systems in which the moon's mass is much less than that of the planet. This work, in contrast, can be applied to the relatively large moons that might be detected around newly-discovered Neptune-mass and super-Earth planets. We conclude that moons are more stable when the planet/moon systems are further from the parent star, the planets are heavier, or the parent stars are lighter. Inclusion of lunar tides allows for significantly longer lifetimes for a massive moon relative to prior formulations. We expect that the semi-major axis of the planet hosting the first detected exomoon around a G-type star is 0.4-0.6 AU and is 0.2-0.4 AU for an M-type star.

We present the first results from our Red Optical Planet Survey (ROPS) to search for low mass planets orbiting late type dwarfs (M5.5V - M9V) in their habitable zones (HZ). Our observations, with the red arm of the MIKE spectrograph (0.5 - 0.9 microns) at the 6.5 m Magellan Clay telescope at Las Campanas Observatory indicate that >= 92 per cent of the flux lies beyond 0.7 microns. We use a novel approach that is essentially a hybrid of the simultaneous iodine and ThAr methods for determining precision radial velocities. We apply least squares deconvolution to obtain a single high S/N ratio stellar line for each spectrum and cross correlate against the simultaneously observed telluric line profile, which we derive in the same way. Utilising the 0.62 - 0.90 micron region, we have achieved an r.m.s. precision of 10 m/s for an M5.5V spectral type star with spectral S/N ~160 on 5 minute timescales. By M8V spectral type, a precision of ~30 m/s at S/N = 25 is suggested, although more observations are needed. An assessment of our errors and scatter in the radial velocity points hints at the presence of stellar radial velocity variations. Of our sample of 7 stars, 2 show radial velocity signals at 6-sigma and 10-sigma of the cross correlation uncertainties. If the signals are planetary in origin, our findings are consistent with estimates of Neptune mass planets that predict a frequency of 13 - 27 per cent for early M dwarfs.Our current analysis indicates the we can achieve a sensitivity that is equivalent to the amplitude induced by a 6 M_Earth planet orbiting in the habitable zone. Based on simulations, we estimate that <10 M_Earth habitable zone planets will be detected in a new stellar mass regime, with <=20 epochs of observations.

We present a new model for Ellipsoidal Variations Induced by a Low-Mass Companion, the EVIL-MC model. We employ several approximations appropriate for planetary systems to substantially increase the computational efficiency of our model relative to more general ellipsoidal variation models and improve upon the accuracy of simpler models. This new approach gives us a unique ability to rapidly and accurately determine planetary system parameters. We use the EVIL-MC model to analyze Kepler Quarter 0-2 (Q0-2) observations of the HAT-P-7 system, an F-type star orbited by a nearly Jupiter-mass companion. Our analysis corroborates previous estimates of the planet-star mass ratio q = (1.10 +/- 0.06) x 10^(-3), and we have revised the planet's dayside brightness temperature to 2680 +10/-20 K. We also find a large difference between the day- and nightside planetary flux, with little nightside emission. Preliminary dynamical+radiative modeling of the atmosphere indicates this result is qualitatively consistent with high altitude absorption of stellar heating. Similar analyses of Kepler and CoRoT photometry of other planets using EVIL-MC will play a key role in providing constraints on the properties of many extrasolar systems, especially given the limited resources for follow-up and characterization of these systems. However, as we highlight, there are important degeneracies between the contributions from ellipsoidal variations and planetary emission and reflection. Consequently, for many of the hottest and brightest Kepler and CoRoT planets, accurate estimates of the planetary emission and reflection, diagnostic of atmospheric heat budgets, will require accurate modeling of the photometric contribution from the stellar ellipsoidal variation.

Planetary satellites are an integral part of the heirarchy of planetary systems. Here we make two predictions concerning their formation. First, primordial satellites, which have an array of distinguishing characteristics, form only around giant planets. If true, the size and duration of a planetary system's protostellar nebula, as well as the location of its snow line, can be constrained by knowing which of its planets possess primordial satellites and which do not. Second, all satellites around terrestrial planets form by impacts. If true, this greatly enhances the constraints that can be placed on the history of terrestrial planets by their satellites' compositions, sizes, and dynamics.

IDENTIKIT was originally developed as a fast approximate scheme for modeling the tidal morphology and kinematics of disk galaxy encounters and mergers. In this form, it was first used to implement an interactive modeling tool for galaxy collisions; tests with artificial data showed that the morphology and kinematics of merging galaxies strongly constrain their initial conditions. This tool is now being applied to real galaxies. More recently, IDENTIKIT has been used to develop a mapping from the present state of a tidal encounter back to the initial conditions; this offers a way to partly automate the search for dynamical models of galaxy encounters. Finally, IDENTIKIT's theoretical applications include a comprehensive way to evaluate the mass and extent of tidal features as functions of halo structure.

Here we present the study of FR Cnc, a young, active and spotted star. We performed analysis of ASAS-3 (The All Sky Automated Survey) data for the years 2002-2008 and amended the value of the rotational period to be 0.826518 d. The amplitude of photometric variations decreased abruptly in the year 2005, while the mean brightness remained the same, which was interpreted as a quick redistribution of spots. BVRc and Ic broad band photometric calibration was performed for 166 stars in FR Cnc vicinity. The photometry at Terskol Observatory shows two brightening episodes, one of which occurred at the same phase as the flare of 2006 November 23. Polarimetric BVR observations indicate the probable presence of a supplementary source of polarization. We monitored FR Cnc spectroscopically during the years 2004-2008. We concluded that the RV changes cannot be explained by the binary nature of FR Cnc. We determined the spectral type of FR Cnc as K7V. Calculated galactic space-velocity components (U, V, W) indicate that FR Cnc belongs to the young disc population and might also belong to the IC 2391 moving group. Based on LiI 6707.8 measurement, we estimated the age of FR Cnc to be between 10-120 Myr. Doppler Tomography was applied to create a starspot image of FR Cnc. We optimized the goodness of fit to the deconvolved profiles for axial inclination, equivalent width and v sin i, finding v sin i=46.2 km s^-1 and i=55 degrees. The starspot distribution of FR Cnc is also of interest since it is one of the latest spectral types to have been imaged. No polar spot was detected on FR Cnc.

New optical HST, Spitzer, GALEX, and Chandra observations of the single-nucleus, luminous infrared galaxy (LIRG) merger IC 883 are presented. The galaxy is a member of the Great Observatories All-sky LIRG Survey (GOALS), and is of particular interest for a detailed examination of a luminous late-stage merger due to the richness of the optically-visible star clusters and the extended nature of the nuclear X-ray, mid-IR, CO and radio emission. In the HST ACS images, the galaxy is shown to contain 156 optically visible star clusters distributed throughout the nuclear regions and tidal tails of the merger, with a majority of visible clusters residing in an arc ~ 3-7 kpc from the position of the mid-infrared core of the galaxy. The luminosity functions of the clusters have an alpha_F435W ~ -2.17+/-0.22 and alpha_F814W ~ -2.01+/-0.21. Further, the colors and absolute magnitudes of the majority of the clusters are consistent with instantaneous burst population synthesis model ages in the range of a few x10^7 - 10^8 yrs (for 10^5 M_sun clusters), but may be as low as few x10^6 yrs with extinction factored in. The X-ray and mid-IR spectroscopy are indicative of predominantly starburst-produced nuclear emission, and the star formation rate is ~ 80 M_sun / yr. The kinematics of the CO emission and the morphology of both the CO and radio emission are consistent with the nuclear starburst being situated in a highly inclined disk 2 kpc in diameter with an infrared surface brightness mu_IR ~ 2x10^11 L_sun kpc^-2, a factor of 10 less than that of the Orion star-forming region. Finally, the detection of the [Ne V] 14.32 um emission line is evidence that an AGN is present. The faintness of the line (i.e., [Ne V] / [Ne II] um ~ 0.01) and the small equivalent width of the 6.2 um PAH feature ($= 0.39\mu$m) are both indicative of a relatively weak AGN. (abridged)

We model the asymmetry of the KOI-13.01 transit lightcurve assuming a gravity-darkened rapidly-rotating host star in order to constrain the system's spin-orbit alignment and transit parameters. We find that our model can reproduce the Kepler lightcurve for KOI-13.01 with a sky-projected alignment of {\lambda} = $23 \pm 4$ degrees and with the star's north pole tilted away from the observer by $48 \pm 4$ degrees (assuming $M_* = 2.05 M_\odot$). With both these determinations, we calculate that the net misalignment between this planet's orbit normal and its star's rotational pole is $56 \pm 4$ degrees. Degeneracies in our geometric interpretation also allow a retrograde spin-orbit angle of $124 \pm 4$ degrees. This is the first spin-orbit measurement to come from gravity darkening, and is one of only a few measurements of the full (not just the sky-projected) spin-orbit misalignment of an extrasolar planet. We also measure accurate transit parameters incorporating stellar oblateness and gravity darkening: $R_* = 1.756 \pm 0.014 R_\odot$, $R_p = 1.445 \pm 0.016 R_{Jup}$, and $i = 85.9 \pm 0.4$ degrees. The new lower planetary radius falls within the planetary mass regime for plausible interior models for the transiting body. A simple initial calculation shows that KOI-13.01's circular orbit is apparently inconsistent with the Kozai mechanism having driven its spin-orbit misalignment; planet-planet scattering and stellar spin migration remain viable mechanisms. Future Kepler data will improve the precision of the KOI-13.01 transit lightcurve, allowing more precise determination of transit parameters and the opportunity to use the Photometric Rossiter-McLaughlin effect to resolve the prograde/retrograde orbit determination degeneracy.

Using a combination of self-consistent and test-particle techniques, Identikit 1 provided a way to vary the initial geometry of a galactic collision and instantly visualize the outcome. Identikit 2 uses the same techniques to define a mapping from the current morphology and kinematics of a tidal encounter back to the initial conditions. By requiring that various regions along a tidal feature all originate from a single disc with a unique orientation, this mapping can be used to derive the initial collision geometry. In addition, Identikit 2 offers a robust way to measure how well a particular model reproduces the morphology and kinematics of a pair of interacting galaxies. A set of eight self-consistent simulations is used to demonstrate the algorithm's ability to search a ten-dimensional parameter space and find near-optimal matches; all eight systems are successfully reconstructed.

In light of the growing interest in searching for low mass, rocky planets, we investigate the impact of starspots on radial velocity searches for earth-mass planets in orbit about M dwarf stars. Since new surveys targeting M dwarfs will likely be carried out at infrared wavelengths, a comparison between V and Y band starspot induced jitter is made, indicating a reduction of up to an order of magnitude when observing in the Y band. The exact reduction in jitter is dependent on the photosphere to spot contrast ratio, with greater improvements at smaller contrasts. We extrapolate a model used to describe solar spot distributions to simulate the spot patterns that we expect to find on M dwarfs. Under the assumption that M dwarfs are near or fully convective, we randomly place starspots on the stellar surface, simulating different levels of spot coverage. Line profiles, distorted by spots are derived and are used to investigate the starspot induced jitter. By making assumptions about the degree of spot activity, detection limits for earth-mass planets in habitable zones are simulated for between 10 and 500 observation epochs. We find that <= 50 epochs are required to detect 1 - 2 MEarth planets (with < 1 per cent false alarm probability) orbiting slowly rotating 0.1 and 0.2 MSun stars. This sensitivity decreases when typical rotation velocities and activity levels for each stellar mass/spectral type are considered. No detections of below 20 MEarth planets are expected for <= 500 observations for the most active stars with vsini >= 20 km/s and dark spots.

We present VLA observations of the neutral hydrogen and radio continuum of NGC 34 (= NGC 17 = Mrk 938). This object is an ideal candidate to study the fate of gas in mergers, since, as shown by an optical study done by Schweizer & Seitzer (2007), it is a gas-rich ("wet") merger remnant of two disk galaxies of unequal mass hosting a strong central starburst and a weak AGN. We detect HI emission from both tidal tails and from nearby galaxies, suggesting that NGC 34 is actually part of a gas-rich group and might have recently interacted with one of its companions. The kinematics of the gas suggests this remnant is forming an outer disk of neutral hydrogen from the gas of the northern tail. We also detect broad HI absorption (514 +/- 21 km/s wide) at both negative and positive velocities with respect to the systemic velocity. This absorption could be explained by the motions of the tidal tails or by the presence of a circumnuclear disk. In addition, we present radio-continuum images that show both nuclear (62.4 +/- 0.3 mJy) and extra-nuclear emission (26.5 +/- 3.0 mJy). The extra-nuclear component is very diffuse and in the shape of two radio lobes, spanning 390 kpc overall. This emission could be a signature of an AGN that has turned off, or it could originate from a starburst-driven superwind. We discuss the possible scenarios that explain our observations, and what they tell us about the location of the gas and the future evolution of NGC 34.

We study a target sample of 68 low-mass objects (with spectral types in the range M4.5-L1) previously selected via photometric and astrometric criteria, as possible members of five young moving groups: the Local Association (Pleiades moving group, age=20 - 150 Myr), the Ursa Mayor group (Sirius supercluster, age=300 Myr), the Hyades supercluster (age=600 Myr), IC 2391 supercluster (age=35 - 55 Myr) and the Castor moving group (age=200 Myr). In this paper we assess their membership by using different kinematic and spectroscopic criteria. We use high resolution echelle spectroscopic observations of the sample to measure accurate radial velocities (RVs). Distances are calculated and compared to those of the moving group from the literature, we also calculate the kinematic Galactic components (U,V,W) of the candidate members and apply kinematic criterion of membership to each group. In addition we measure rotational velocities (v sin i) to place further constraints on membership of kinematic members. We find that 49 targets have young disk kinematics and that 36 of them possibly belong to one of our five moving groups. From the young disk target ob jects, 31 have rotational velocities in agreement with them belonging to the young disk population. We also find that one of our moving group candidates, 2MASS0123- 3610, is a low-mass double lined spectroscopic binary, with probable spectral types around M7.

We present new dynamical models of the merger remnant NGC 7252 which include star formation simulated according to various phenomenological rules. By using interactive software to match our model with the observed morphology and gas velocity field, we obtain a consistent dynamical model for NGC 7252. In our models, this proto-elliptical galaxy formed by the merger of two similar gas-rich disk galaxies which fell together with an initial pericentric separation of ~2 disk scale lengths approximately 620 Myr ago. Results from two different star formation rules--- density-dependent and shock-induced--- show significant differences in star formation during and after the first passage. Shock-induced star formation yields a prompt and wide-spread starburst at the time of first passage, while density-dependent star formation predicts a more slowly rising and centrally concentrated starburst. A comparison of the distributions and ages of observed clusters with results of our simulations favors shock-induced mechanism of star formation in NGC 7252. We also present simulated color images of our model of NGC 7252, constructed by incorporating population synthesis with radiative transfer and dust attenuation. Overall the predicted magnitudes and colors of the models are consistent with observations, although the simulated tails are fainter and redder than observed. We suggest that a lack of star formation in the tails, reflected by the redder colors, is due to an incomplete description of star formation in our models rather than insufficient gas in the tails.

An analysis of data from the Spitzer Space Telescope, Hubble Space Telescope, Chandra X-ray Observatory, and AKARI Infrared Astronomy Satellite is presented for the z=0.036 merging galaxy system II Zw 096 (CGCG 448-020). Because II Zw 096 has an infrared luminosity of log(L_IR/L_sun) = 11.94, it is classified as a Luminous Infrared Galaxy (LIRG), and was observed as part of the Great Observatories All-sky LIRG Survey (GOALS). The Spitzer data suggest that 80% of the total infrared luminosity comes from an extremely compact, red source not associated with the nuclei of the merging galaxies. The Spitzer mid-infrared spectra indicate no high-ionization lines from a buried active galactic nucleus in this source. The strong detection of the 3.3 micron and 6.2 micron PAH emission features in the AKARI and Spitzer spectra also implies that the energy source of II Zw 096 is a starburst. Based on Spitzer infrared imaging and AKARI near-infrared spectroscopy, the star formation rate is estimated to be 120 M_sun/yr and > 45 M_sun/yr, respectively. Finally, the high-resolution B, I, and H-band images show many star clusters in the interacting system. The colors of these clusters suggest at least two populations - one with an age of 1-5 Myr and one with an age of 20-500 Myr, reddened by 0-2 magnitudes of visual extinction. The masses of these clusters span a range between 10^6-10^8 M_sun. This starburst source is reminiscent of the extra-nuclear starburst seen in NGC 4038/9 (the Antennae Galaxies) and Arp 299 but approximately an order of magnitude more luminous than the Antennae. The source is remarkable in that the off-nuclear infrared luminosity dominates the enitre system.

We report the observation of a cloud system on Titan that remained localized near 40S latitude and 60W longitude for at least 34 hours. Ground-based observations obtained with the SINFONI imaging spectrograph at the Very Large Telescope over 4 consecutive nights recorded the lifetime and altitude of the unresolved cloud system. Concomitant measurements made by Cassini/VIMS over 3 hours resolved changes in the altitude and opacity of individual regions within the system during this time. Clouds are measured from 13 to 37 km altitude with optical depths per pixel ranging from $\tau$=0.13 to 7. Short timescale rise times are consistent with previous measurements of the evolution of mid-latitude clouds; however the long timescale localization of the cloud structure is unexplained. We speculate about the role of meso-scale circulation in relation to cloud formation.

Recent results comparing interacting galaxies to the mass-metallicity relation show that their nuclear oxygen abundances are unexpectedly low. We present analysis of N-body/SPH numerical simulations of equal-mass mergers that confirm the hypothesis that these underabundances are accounted for by radial inflow of low-metallicity gas from the outskirts of the two merging galaxies. The underabundances arise between first and second pericenter, and the simulated abundance dilution is in good agreement with observations. The simulations further predict that the radial metallicity gradients of the disk galaxies flatten shortly after first passage, due to radial mixing of gas. These predictions will be tested by future observations of the radial metallicity distributions in interacting galaxies.

The detailed study of the exoplanetary systems HD189733 and HD209458 has given rise to a wealth of exciting information on the physics of exoplanetary atmospheres. To further our understanding of the make-up and processes within these atmospheres we require a larger sample of bright transiting planets. We have began a project to detect more bright transiting planets in the southern hemisphere by utilising precision radial-velocity measurements. We have observed a constrained sample of bright, inactive and metal-rich stars using the HARPS instrument and here we present the current status of this project, along with our first discoveries which include a brown dwarf/extreme-Jovian exoplanet found in the brown dwarf desert region around the star HD191760 and improved orbits for three other exoplanetary systems HD48265, HD143361 and HD154672. Finally, we briefly discuss the future of this project and the current prospects we have for discovering more bright transiting planets.

BD+20 1790 is a young active, metal-rich, late-type K5Ve star. We have undertaken a study of stellar activity and kinematics for this star over the past few years. Previous results show a high level of stellar activity, with the presence of prominence-like structures, spots on the surface and strong flare events, despite the moderate rotational velocity of the star. In addition, radial velocity variations with a semi-amplitude of up to 1 km/s were detected. We investigated the nature of these radial velocity variations, in order to determine whether they are due to stellar activity or the reflex motion of the star induced by a companion. We have analysed high-resolution echelle spectra and also two-band photometry was obtained to produce the light curve and determine the photometric period. Based upon the analysis of the bisector velocity span, as well as spectroscopic indices of chromospheric indicators and taking into account the photometric analysis, we report that the best explanation for the RV variation is the presence of a sub-stellar companion. The Keplerian fit of the RV data yields a solution for a close-in massive planet with an orbital period of 7.78 days. The presence of the close-in massive planet could also be an interpretation for the high level of stellar activity detected. Since the RV data are not part of a planet search program, we can consider our results as a serendipitous evidence of a planetary companion. To date, this is the youngest main sequence star for which a planetary candidate has been reported.

We associate 132 low-mass ultracool dwarfs in the southern hemisphere as candidate members of five moving groups using photometric and astrometric selection techniques. Of these objects, we present high resolution spectroscopy for seven candidates and combine these with previous measurements from the literature to determine spectral types and radial velocities. We thus constrain distance and space motion spectroscopically, allowing the kinematic membership of the moving groups to be assessed. Possible membership of moving groups has allowed ages and metallicities to be constrained for these objects and evolutionary models have been used to estimate their mass. We estimate that up to ~75 of our candidate moving group members should be genuine, and discuss future work that will confirm and exploit this major new sample.

Several processes can cause the shape of an extrasolar giant planet's shadow, as viewed in transit, to depart from circular. In addition to rotational effects, cloud formation, non-homogenous haze production and movement, and dynamical effects (winds) could also be important. When such a planet transits its host star as seen from Earth, the asphericity will introduce a deviation in the transit lightcurve relative to the transit of a perfectly spherical (or perfectly oblate) planet. We develop a theoretical framework to interpret planet shapes. We then generate predictions for transiting planet shapes based on a published theoretical dynamical model of HD189733b. Using these shape models we show that planet shapes are unlikely to introduce detectable lightcurve deviations (those >~1e-5 of the host star), but that the shapes may lead to astrophysical sources of systematic error when measuring planetary oblateness, transit time, and impact parameter.

We use signal enhancement techniques and a matched filter analysis to search for the K band spectroscopic absorption signature of the close orbiting extrasolar giant planet, HD 189733b. With timeseries observations taken with NIRSPEC at the Keck II telescope, we investigate the relative abundances of H2O and carbon bearing molecules, which have now been identified in the dayside spectrum of HD 189733b. We detect a candidate planet signature with a low level of significance, close to the ~153 km/s velocity amplitude of HD 189733b. However, some systematic variations, mainly due to imperfect telluric line removal, remain in the residual spectral timeseries in which we search for the planetary signal. The robustness of our candidate signature is assessed, enabling us to conclude that it is not possible to confirm the presence of any planetary signal which appears at Fp/F* contrasts deeper than the 95.4 per cent confidence level. Our search does not enable us to detect the planet at a contrast ratio of Fp/F* = 1/1920 with 99.9 per cent confidence. We also investigate the effect of model uncertainties on our ability to reliably recover a planetary signal. The use of incorrect temperature, model opacity wavelengths and model temperature-pressure profiles have important consequences for the least squares deconvolution procedure that we use to boost the S/N ratio in our spectral timeseries observations. We find that mismatches between the empirical and model planetary spectrum may weaken the significance of a detection by ~30-60 per cent, thereby potentially impairing our ability to recover a planetary signal with high confidence.

Main-sequence stars earlier than spectral type ~F6 or so are expected to rotate rapidly due to their radiative exteriors. This rapid rotation leads to an oblate stellar figure. It also induces the photosphere to be hotter (by up to several thousand Kelvin) at the pole than at the equator as a result of a process called gravity darkening that was first predicted by von Zeipel (1924). Transits of extrasolar planets across such a non-uniform, oblate disk yield unusual and distinctive lightcurves that can be used to determine the relative alignment of the stellar rotation pole and the planet orbit normal. This spin-orbit alignment can be used to constrain models of planet formation and evolution. Orderly planet formation and migration within a disk that is coplanar with the stellar equator will result in spin-orbit alignment. More violent planet-planet scattering events should yield spin-orbit misaligned planets. Rossiter-McLaughlin measurements of transits of lower-mass stars show that some planets are spin-orbit aligned, and some are not. Since Rossiter-McLaughlin measurements are difficult around rapid rotators, lightcurve photometry may be the best way to determine the spin-orbit alignment of planets around massive stars. The Kepler mission will monitor ~10^4 of these stars within its sample. The lightcurves of any detected planets will allow us to probe the planet formation process around high-mass stars for the first time.

We present spectroscopic rotation velocities (v sin i) for 56 M dwarf stars using high resolution HET HRS red spectroscopy. In addition we have also determined photometric effective temperatures, masses and metallicities ([Fe/H]) for some stars observed here and in the literature where we could acquire accurate parallax measurements and relevant photometry. We have increased the number of known v sin is for mid M stars by around 80% and can confirm a weakly increasing rotation velocity with decreasing effective temperature. Our sample of v sin is peak at low velocities (~3 km/s). We find a change in the rotational velocity distribution between early M and late M stars, which is likely due to the changing field topology between partially and fully convective stars. There is also a possible further change in the rotational distribution towards the late M dwarfs where dust begins to play a role in the stellar atmospheres. We also link v sin i to age and show how it can be used to provide mid-M star age limits. When all literature velocities for M dwarfs are added to our sample there are 198 with v sin i less than or equal to 10 km/s and 124 in the mid-to-late M star regime (M3.0-M9.5) where measuring precision optical radial-velocities is difficult. In addition we also search the spectra for any significant Halpha emission or absorption. 43% were found to exhibit such emission and could represent young, active objects with high levels of radial-velocity noise. The tables presented here will aid any future M star planet search target selection to extract stars that will exhibit low radial-velocity jitter.

We present Doppler imaging and a Balmer line analysis of the weak-line T Tauri star TWA 17. Spectra were taken in 2006 with the UCL Echelle Spectrograph on the Anglo-Australian Telescope. Using least-squares deconvolution to improve the effective signal-to-noise ratio we produced a Doppler map of the surface spot distribution. This shows similar features to maps of other rapidly rotating T Tauri stars, i.e. a polar spot with more spots extending out of it down to the equator. In addition to the photospheric variability, the chromospheric variability was studied using the Balmer emission. The mean H-alpha profile has a narrow component consistent with rotational broadening and a broad component extending out to +/-220 km/s. The variability in H-alpha suggests that the chromosphere has at least one slingshot prominence 3 stellar radii above the surface.

We present evidence for condensed phase methane precipitation near Xanadu using nine nights of observations from the SINFONI integral-field spectrograph at the Very Large Telescope and imaging analysis with empirical surface subtraction. Radiative transfer models are used to support the imaging technique by simulating the spectrometer datacubes and testing for variations in both the surface reflectivity spectrum and atmospheric opacity. We use the models and observations together to argue against artifacts that may arise in the image analysis. High phase angle observations from Cassini/VIMS are used to test against surface scattering artifacts that may be confused with sources of atmospheric opacity. Although changes in the surface reflectivity spectrum can reproduce observations from a particular viewing geometry on a given night, multiple observations are best modeled by condensed-phase methane opacity near the surface. These observations and modeling indicate that the condensed-phase methane opacity observed with this technique occurs predominantly near Xanadu and is most likely due to precipitation.

Clouds on Titan result from the condensation of methane and ethane and, as on other planets, are primarily structured by circulation of the atmosphere. At present, cloud activity mainly occurs in the southern (summer) hemisphere, arising near the pole and at mid-latitudes from cumulus updrafts triggered by surface heating and/or local methane sources, and at the north (winter) pole, resulting from the subsidence and condensation of ethane-rich air into the colder troposphere. General circulation models predict that this distribution should change with the seasons on a 15-year timescale, and that clouds should develop under certain circumstances at temperate latitudes (~40\degree) in the winter hemisphere. The models, however, have hitherto been poorly constrained and their long-term predictions have not yet been observationally verified. Here we report that the global spatial cloud coverage on Titan is in general agreement with the models, confirming that cloud activity is mainly controlled by the global circulation. The non-detection of clouds at latitude ~40\degree N and the persistence of the southern clouds while the southern summer is ending are, however, both contrary to predictions. This suggests that Titan's equator-to-pole thermal contrast is overestimated in the models and that its atmosphere responds to the seasonal forcing with a greater inertia than expected.

The Great Observatories All-sky LIRG Survey (GOALS) combines data from NASA's Spitzer, Chandra, Hubble and GALEX observatories, together with ground-based data into a comprehensive imaging and spectroscopic survey of over 200 low redshift Luminous Infrared Galaxies (LIRGs). The LIRGs are a complete subset of the IRAS Revised Bright Galaxy Sample (RBGS). The LIRGs targeted in GOALS span the full range of nuclear spectral types defined via traditional optical line-ratio diagrams as well as interaction stages. They provide an unbiased picture of the processes responsible for enhanced infrared emission in galaxies in the local Universe. As an example of the analytic power of the multi-wavelength GOALS dataset, we present data for the interacting system VV 340 (IRAS F14547+2449). Between 80-95% of the total far-infrared emission (or about 5E11 solar luminosities) originates in VV 340 North. While the IRAC colors of VV 340 North and South are consistent with star-forming galaxies, both the Spitzer IRS and Chandra ACIS data indicate the presence of a buried AGN in VV 340 North. The GALEX far and near-UV fluxes imply a extremely large infrared "excess" (IRX) for the system (IR/FUV = 81) which is well above the correlation seen in starburst galaxies. Most of this excess is driven by VV 340 N, which alone has an IR excess of nearly 400. The VV 340 system seems to be comprised of two very different galaxies - an infrared luminous edge-on galaxy (VV 340 North) that dominates the long-wavelength emission from the system and which hosts a buried AGN, and a face-on starburst (VV 340 South) that dominates the short-wavelength emission.

We report the discovery of a brown dwarf on an eccentric orbit and with a semimajor axis that places it in the brown dwarf desert region around the star HD191760. The star has a spectral type of G3IV/V and a metallicity ([Fe/H]) of 0.29 dex. HD191760 adds to the small number of metal-rich stars with brown dwarf companions. The brown dwarf (HD191760b) is found to have an orbital period of 505.57+/-0.40 days and semimajor axis of 1.35+/-0.01 AU, placing it firmly in the brown dwarf desert. The eccentricity of HD191760b is found to be 0.63+/-0.01, meaning it reaches as close as 0.5 AU from the host star. Dynamical simulations indicate that no inner planets could reside at separations beyond ~0.17 AU due to the disastrous gravity imposed by HD191760b. In addition to these first results we also refine the orbits found for the exoplanets around the stars HD48265, HD143361 and HD154672. All 1-planet solutions are in agreement with those previously published by the Magellan Planet Search.

We report the results from a comprehensive study of 74 ultraluminous infrared galaxies (ULIRGs) and 34 Palomar-Green (PG) quasars within z ~ 0.3$ observed with the Spitzer Infrared Spectrograph (IRS). The contribution of nuclear activity to the bolometric luminosity in these systems is quantified using six independent methods that span a range in wavelength and give consistent results within ~ +/-10-15% on average. The average derived AGN contribution in ULIRGs is ~35-40%, ranging from ~15-35% among "cool" (f_25/f_60 =< 0.2) optically classified HII-like and LINER ULIRGs to ~50 and ~75% among warm Seyfert 2 and Seyfert 1 ULIRGs, respectively. This number exceeds ~80% in PG QSOs. ULIRGs fall in one of three distinct AGN classes: (1) objects with small extinctions and large PAH equivalent widths are highly starburst-dominated; (2) systems with large extinctions and modest PAH equivalent widths have larger AGN contributions, but still tend to be starburst-dominated; and (3) ULIRGs with both small extinctions and small PAH equivalent widths host AGN that are at least as powerful as the starbursts. The AGN contributions in class 2 ULIRGs are more uncertain than in the other objects, and we cannot formally rule out the possibility that these objects represent a physically distinct type of ULIRGs. A morphological trend is seen along the sequence (1)-(2)-(3), in general agreement with the standard ULIRG - QSO evolution scenario and suggestive of a broad peak in extinction during the intermediate stages of merger evolution. However, the scatter in this sequence, implies that black hole accretion, in addition to depending on the merger phase, also has a strong chaotic/random component, as in local AGN. (abridged)

The paper was withdrawn due to another possible solution to the dataset that is significantly different in nature. This issue will be addressed shortly and clarified with an additional data point.

By combining test-particle and self-consistent techniques, we have developed a method to rapidly explore the parameter space of galactic encounters. Our method, implemented in an interactive graphics program, can be used to find the parameters required to reproduce the observed morphology and kinematics of interacting disk galaxies. We test this system on an artificial data-set of 36 equal-mass merging encounters, and show that it is usually possible to reproduce the morphology and kinematics of these encounters and that a good match strongly constrains the encounter parameters.

We have begun a metal-rich planet search project using the HARPS instrument in La Silla, Chile to target planets with a high potential to transit their host star and add to the number of bright benchmark transiting planets. The sample currently consists of 100, bright (7.5 </= V </= 9.5) solar-type stars (0.5 </= B-V </= 0.9) in the southern hemisphere which are both inactive (logR'HK </= -4.5) and metal-rich ([Fe/H] >/= 0.1 dex). We determined the chromospheric activity and metallicity status of our sample using high resolution FEROS spectra. We also introduce the first result from our HARPS planet search and show that the radial-velocity amplitude of this star is consistent with an orbiting planetary-mass companion (i.e. Msini < 0.5MJ) with a period of ~5 days. We are currently engaged in follow-up to confirm this signal as a bonafide orbiting planet.

We have designed and tested an in-plane echelle spectrograph configured to investigate precision radial velocities from ground-based near-infrared observations. The spectrograph operates across the spectral range of 0.9-1.7 mm at a spectral resolution of R = 50,000, and uses a liquid nitrogen-cooled HAWAII 1K detector. Repeated measurements of the Earth's rotation via integrated Sunlight with two different instrument arrangements in the near infrared Y band have produced radial velocities with ~10 m/s RMS over a period of several hours. The most recent instrument configuration has achieved an unbinned RMS of 7 m/s and suggests that infrared radial velocity precisions may be able to approach those achieved at optical wavelengths.

We have carried out a search for the 2.14 micron spectroscopic signature of the close orbiting extrasolar giant planet, HD 179949b. High cadence time series spectra were obtained with the CRIRES spectrograph at VLT1 on two closely separated nights. Deconvolution yielded spectroscopic profiles with mean S/N ratios of several thousand, enabling the near infrared contrast ratios predicted for the HD 179949 system to be achieved. Recent models have predicted that the hottest planets may exhibit spectral signatures in emission due to the presence of TiO and VO which may be responsible for a temperature inversion high in the atmosphere. We have used our phase dependent orbital model and tomographic techniques to search for the planetary signature under the assumption of an absorption line dominated atmospheric spectrum, where T and V are depleted from the atmospheric model, and an emission line dominated spectrum, where TiO and VO are present. We do not detect a planet in either case, but the 2.120 - 2.174 micron wavelength region covered by our observations enables the deepest near infrared limits yet to be placed on the planet/star contrast ratio of any close orbiting extrasolar giant planet system. We are able to rule out the presence of an atmosphere dominated by absorption opacities in the case of HD 179949b at a contrast ratio of F_p/F_* ~ 1/3350, with 99 per cent confidence.

We present the first measurements of surface differential rotation on a pre-main sequence binary system. Using intensity (Stokes I) and circularly polarised (Stokes V) timeseries spectra, taken over eleven nights at the Anglo-Australian Telescope (AAT), we incorporate a solar-like differential rotation law into the surface imaging process. We find that both components of the young, 18 Myr, HD 155555 (V824 Ara, G5IV + K0IV) binary system show significant differential rotation. The equator-pole laptimes as determined from the intensity spectra are 80 days for the primary star and 163 days for the secondary. Similarly for the magnetic spectra we obtain equator-pole laptimes of 44 and 71 days respectively, showing that the shearing timescale of magnetic regions is approximately half that found for stellar spots. Both components are therefore found to have rates of differential rotation similar to those of the same spectral type main sequence single stars. The results for HD 155555 are therefore in contrast to those found in other, more evolved, binary systems where negligible or weak differential rotation has been discovered. We discuss two possible explanations for this; firstly that at the age of HD 155555 binary tidal forces have not yet had time to suppress differential rotation, secondly that the weak differential rotation previously observed on evolved binaries is a consequence of their large convection zone depths. We suggest that the latter is the more likely solution and show that both temperature and convection zone depth (from evolutionary models) are good predictors of differential rotation strength. Finally, we also examine the possible consequences of the measured differential rotation on the interaction of binary star coronae.

We present the results from high-resolution spectroscopic measurements to determine metallicities and activities of bright stars in the southern hemisphere. We have measured the iron abundances ([Fe/H]) and chromospheric emission indices (logR'HK) of 353 solar-type stars with V=7.5-9.5. [Fe/H] abundances are determined using a custom chi^2 fitting procedure within a large grid of Kurucz model atmospheres. The chromospheric activities were determined by measuring the amount of emission in the cores of the strong CaII HK lines. The sample of metallicities has been compared to other [Fe/H] determinations and was found to agree with these at the +/-0.05 dex level for spectroscopic values and at the +/-0.1 dex level for photometric values. The distribution of chromospheric activities is found to be described by a bimodal distribution, agreeing well with the conclusions from other works. Also an analysis of Maunder Minimum status was attempted and it was found that 6+/-4 stars in the sample could be in a Maunder Minimum phase of their evolution and hence the Sun should only spend a few per cent of its main sequence lifetime in Maunder Minimum.

We present Doppler imaging and Balmer line analysis of the weak-line T Tauri star TWA 6. Using this data we have made one of the first attempts to measure differential rotation in a T Tauri star, and the first detection of a slingshot prominence in such a star. We also show the most direct evidence to date of the existence of solar-type plages in a star other than the Sun. Observations were made over six nights: 11-13th February 2006 and 18-20th February 2006, when spectra were taken with the UCL Echelle Spectrograph on the 3.9-m Anglo-Australian Telescope. Using least-squares deconvolution to improve the effective signal--to--noise ratio we produced two Doppler maps. These show similar features to maps of other rapidly rotating T Tauri stars, i.e. a polar spot with more spots extending out of it down to equator. Comparison of the two maps was carried out to measure the differential rotation. Cross-correlation and parameter fitting indicates that TWA 6 does not have detectable differential rotation. The Balmer emission of the star was studied. The mean H-alpha profile has a narrow component consistent with rotational broadening and a broad component extending out to 250km/s. The variability in H-alpha suggests that the chromosphere has active regions that are cospatial with the spots in the photosphere, similar to the 'plages' observed on the Sun. In addition the star has at least one slingshot prominence 3 stellar radii above the surface - the first such detection in a T Tauri star.

We obtained 238 spectra of the close orbiting extrasolar giant planet HD 189733b with resolution R ~ 15,000 during one night of observations with the near infrared spectrograph, NIRSPEC, at the Keck II Telescope. We have searched for planetary absorption signatures in the 2.0 - 2.4 micron region where H_2O and CO are expected to be the dominant atmospheric opacities. We employ a phase dependent orbital model and tomographic techniques to search for the planetary absorption signatures in the combined stellar and planetary spectra. Because potential absorption signatures are hidden in the noise of each single exposure, we use a model list of lines to apply a spectral deconvolution. The resulting mean profile possesses a S/N ratio that is 20 times greater than that found in individual lines. Our spectral timeseries thus yields spectral signatures with a mean S/N = 2720. We are unable to detect a planetary signature at a contrast ratio of log_10(F_p/F_*) = -3.40, with 63.8 per cent confidence. Our findings are not consistent with model predictions which nevertheless give a good fit to mid-infrared observations of HD 189733. The 1-sigma result is a factor of 1.7 times less than the predicted 2.185 micron planet/star flux ratio of log_10(F_p/F_*) ~ -3.16.

It is shown herein that planets with eccentric orbits are more likely to transit than circularly orbiting planets with the same semimajor axis by a factor of (1-e^2)^{-1}. If the orbital parameters of discovered transiting planets are known, as from follow-up radial velocity observations, then the transit-detected planet population is easily debiased of this effect. The duration of a planet's transit depends upon of its eccentricity and longitude of periastron; transits near periastron are shorter, and those near apoastron last longer, for a given impact parameter. If fitting for the stellar radius with the other transit parameters, this effect causes a systematic error in the resulting measurements. If the stellar radius is instead held fixed at a value measured independently, then it is possible to place a lower limit on the planet's eccentricity using photometry alone. Orbital accelerations cause a difference in the planet's ingress and egress durations that lead to an asymmetry in the transit lightcurve that could be used along with the transit velocity measurement to uniquely measure the planet's eccentricity and longitude of periapsis. However, the effect is too small to be measured with current technology. The habitability of transiting terrestrial planets found by Kepler depends on those planets' orbital eccentricities. While Kepler will be able to place lower limits on those planets' orbital eccentricity, the actual value for any given planet will likely remain unknown.

To aid in the physical interpretation of planetary radii constrained through observations of transiting planets, or eventually direct detections, we compute model radii of pure hydrogen-helium, water, rock, and iron planets, along with various mixtures. Masses ranging from 0.01 Earth masses to 10 Jupiter masses at orbital distances of 0.02 to 10 AU are considered. For hydrogen-helium rich planets, our models are the first to couple planetary evolution to stellar irradiation over a wide range of orbital separations (0.02 to 10 AU) through a non-gray radiative-convective equilibrium atmosphere model. Stellar irradiation retards the contraction of giant planets, but its effect is not a simple function of the irradiation level: a planet at 1 AU contracts as slowly as a planet at 0.1 AU. For hydrogen-helium planets, we consider cores up to 90% of the total planet mass, comparable to those of Uranus and Neptune. If "hot Neptunes" have maintained their original masses and are not remnants of more massive planets, radii of 0.30-0.45 times Jupiter's radius are expected. Water planets are ~40-50% larger than rocky planets, independent of mass. Finally, we provide tables of planetary radii at various ages and compositions, and for ice-rock-iron planets we fit our results to analytic functions, which will allow for quick composition estimates, given masses and radii, or mass estimates, given only planetary radii. These results will assist in the interpretation of observations for both the current transiting planet surveys as well as upcoming space missions, including CoRoT and Kepler.

The SuperWASP-I instrument observed 6.7 million stars between 8-15 mag from La Palma during the 2004 May-September season. Our transit-hunting algorithm selected 11,626 objects from the 184,442 stars within the range RA 18hr-21hr. We describe our thorough selection procedure whereby catalogue information is exploited along with careful study of the SuperWASP data to filter out, as far as possible, transit mimics. We have identified 35 candidates which we recommend for follow-up observations.

We have performed photometric observations of nearly 7 million stars with 8 < V < 15 with the SuperWASP-North instrument from La Palma between 2004 May-September. Fields in the RA range 17-18hr, yielding over 185,000 stars with sufficient quality data, have been searched for transits using a modified box least-squares (BLS) algorithm. We find a total of 58 initial transiting candidates which have high S/N in the BLS, show multiple transit-like dips and have passed visual inspection. Analysis of the blending and inferred planetary radii for these candidates leaves a total of 7 transiting planet candidates which pass all the tests plus 4 which pass the majority. We discuss the derived parameters for these candidates and their properties and comment on the implications for future transit searches.

We present a search for the near infrared spectroscopic signature of the close orbiting extrasolar giant planet HD 75289b. We obtained ~230 spectra in the wavelength range 2.18 - 2.19 microns using the Phoenix spectrograph at Gemini South. By considering the direct spectrum, derived from irradiated model atmospheres, we search for the absorption profile signature present in the combined star and planet light. Since the planetary spectrum is separated from the stellar spectrum at most phases, we apply a phase dependent orbital model and tomographic techniques to search for absorption signatures. Because the absorption signature lies buried in the noise of a single exposure we apply a multiline deconvolution to the spectral lines available in order to boost the effective S/N ratio of the data. The wavelength coverage of 80 angstroms is expected to contain ~100 planetary lines, enabling a mean line with S/N ratio of ~800 to be achieved after deconvolution. We are nevertheless unable to detect the presence of the planet in the data and carry out further simulations to show that broader wavelength coverage should enable a planet like HD 75289b to be detected with 99.9 per cent (4 sigma) confidence. We investigate the sensitivity of our method and estimate detection tolerances for mismatches between observed and model planetary atmospheres.

Galaxy interactions are known to trigger starbursts. The young star clusters formed in mergers may be young globular clusters. The ages of these young star clusters yield the timing of interaction-triggered star formation and provide an important way to reconstruct the history of merging galaxies. Here we present the first results from our investigation into age and metallicity of twelve young clusters in the merging galaxy pair NGC 4676, using spectra from the multi-object Low Resolution Imaging Spectrometer (LRIS) on Keck. For ten clusters, comparison of the Balmer emission lines with model equivalent widths (EWs) yields ages less than 10 Myr. Two spectra display Balmer absorption lines typical of star clusters dominated by A-type stars, with estimated ages of about 170 Myr. These ages are comparable to the dynamical age of the tidal tails and are consistent with star formation triggered during the first passage of the pair. The locations of these two clusters in the tidal tails are generally consistent with predictions of shock-induced star formation models. One of these older objects appears unresolved on the image and is luminous enough to qualify as a young globular cluster. Using EWs of the diagnostic lines [OII] and [OIII], we obtain oxygen abundances in the range 7.3 < 12+log(O/H) < 9.0. These values show a nearly flat distribution along the northern tail, suggesting efficient gas mixing in the tail.

The Spitzer-COSMOS survey (S-COSMOS) is a Legacy program (Cycles 2+3) designed to carry out a uniform deep survey of the full 2 sq deg COSMOS field in all seven Spitzer bands (3.6, 4.5, 5.6, 8.0, 24.0, 70.0, 160.0 u). This paper describes the survey parameters, mapping strategy, data reduction procedures, achieved sensitivities to date, and the complete data set for future reference. We show that the observed infrared backgrounds in the S-COSMOS field are within 10% of the predicted background levels. The fluctuations in the background at 24u have been measured and do not show any significant contribution from cirrus, as expected. In addition, we report on the number of asteroid detections in the low galactic latitude COSMOS field. We use the Cycle 2 S-COSMOS data to determine preliminary number counts, and compare our results with those from previous Spitzer Legacy surveys (e.g. SWIRE, GOODS). The results from this "first analysis" confirm that the S-COSMOS survey will have sufficient sensitivity with IRAC to detect ~ L* disks and spheroids out to z ~ 3, and with MIPS to detect ultraluminous starbursts and AGN out to z ~3 at 24u and out to z ~1.5-2 at 70u and 160u.

Observations of stellar prominences on young rapidly rotating stars provide unique probes of their magnetic fields out to many stellar radii. We compare two independently obtained datasets of the K3 dwarf Speedy Mic using the Anglo-Australian Telescope (AAT) and the ESO Very Large Telescope (VLT). Taken more than a fortnight apart they provide the first insight into the evolution of the prominence system on such a young rapidly rotating star. The largest prominences observed transiting the stellar disc are found at very similar rotational phases between the epochs. This suggests that the magnetic structures supporting the prominences retain their identity on a two to three week timescale. We identify prominences as transient absorption features in all lines of the Hydrogen Balmer series down to H_10. We use the ratios of the prominence EWs in these lines to determine their column densities in the first excited state of hydrogen. We determine the optical depths, finding prominences to be rather optically thick (\tau \approx 20) in the H_alpha line. The total hydrogen column density and thus the prominence masses are determined via observations of the CaII H&K lines. We find typical masses for four of the largest prominences to be in the range 0.5 - 2.3 x10^{14} kg, slightly larger than giant solar prominence masses. Rotationally modulated emission is seen outside of the H_alpha line. These loops of emission are shown to be caused by prominences seen off the stellar disc. We find that all of the large emission loops can be associated with prominences we see transiting the stellar disc. This combined with the fact that many prominences appear to eclipse the off-disc emission of others, strongly suggests that the prominence system is highly flattened and likely confined to low stellar latitudes.

We report the discovery of a medium-strength (~0.5kG) magnetic field on the young, massive star tauSco (B0.2V), which becomes the third-hottest magnetic star known. Circularly polarized Zeeman signatures are clearly detected in observations collected mostly with the ESPaDOnS spectropolarimeter, recently installed on the 3.6-m Canada-France-Hawaii Telescope; temporal variability is also clearly established in the polarimetry, and can be unambiguously attributed to rotational modulation with a period close to 41d. Archival UV spectra confirm that this modulation repeats over timescales of decades. By reconstructing the large-scale structure of its magnetic topology, we find that the magnetic structure is unusually complex for a hot star. The surface topology is dominated by a potential field, although a moderate toroidal component is probably present. We fail to detect intrinsic temporal variability of the magnetic structure over the 1.5-yr period of our spectropolarimetric observations (in agreement with the stable temporal variations of the UV spectra), and infer that any differential surface rotation must be very small. The topology of the extended magnetic field that we derive from the photospheric magnetic maps is also more complex than a global dipole, and features in particular a significantly warped torus of closed magnetic loops encircling the star (tilted at about 90deg to the rotation axis), with additional, smaller, networks of closed field lines. This topology appears to be consistent with the exceptional Xray properties of tauSco and also provides a natural explanation of the variability observed in wind-formed UV lines. We conclude that its magnetic field is most probably a fossil remnant from the star-formation stage.

Understanding how cool stars produce magnetic fields within their interiors is crucial for predicting the impact of such fields, such as the activity cycle of the Sun. In this respect, studying fully convective stars enables us to investigate the role of convective zones in magnetic field generation. We produced a magnetic map of a rapidly rotating, very-low-mass, fully convective dwarf through tomographic imaging from time series of spectropolarimetric data. Our results, which demonstrate that fully convective stars are able to trigger axisymmetric large-scale poloidal fields without differential rotation, challenge existing theoretical models of field generation in cool stars.

We have used maximum entropy eclipse mapping to recover images of the visual surface brightness distribution of the primary component of the RS CVn eclipsing binary SV Cam, using high-precision photometry data obtained during three primary eclipses with STIS aboard the Hubble Space Telescope. These were augmented by contemporaneous ground-based photometry secured around the rest of the orbit. The goal of these observations was to determine the filling factor and size distribution of starspots too small to be resolved by Doppler imaging. The information content of the final image and the fit to the data were optimised with respect to various system parameters using the chi^2 landscape method, using an eclipse mapping code that solves for large-scale spot coverage. It is only with the unprecedented photometric precision of the HST data (0.00015 mag) that it is possible to see strong discontinuities at the four contact points in the residuals of the fit to the lightcurve. These features can only be removed from the residual lightcurve by the reduction of the photospheric temperature, to synthesise high unresolvable spot coverage, and the inclusion of a polar spot. We show that this spottedness of the stellar surface can have a significant impact on the determination of the stellar binary parameters and the fit to the lightcurve by reducing the secondary radius from 0.794 +/- 0.009 Rsun to 0.727 +/- 0.009 Rsun. This new technique can also be applied to other binary systems with high precision spectrophotometric observations.

We present new observations of the prominence system on the K3 dwarf Speedy Mic (BO Mic, HD 197890). Using an improved technique to track the absorption features in H_alpha we find a very active prominence system with approximately ten prominences on the observable hemisphere per rotation. From a total of 25 prominences we find an average axial distance of 2.85+/-0.54 R_* which is twice the co-rotation radius above the stellar surface. We discuss the consequences of these observations on the nature of the supporting magnetic structures. Two consecutive nights, with complete phase coverage, combined with a further night after a three night gap allow us to study the evolution of the prominence system on two different timescales. Several of the prominences have counterparts at similar phases on consecutive nights. During this interval many prominences show evidence for evolution in their heights and phases of observation. Five nights (13 rotation cycles) later we recover many prominences at approximately the same phases. Whilst individual prominences change axial distances or appear/re-appear from night-to-night the underlying prominence supporting structures appear to be stable over as many as 13 stellar rotations.

We present high resolution images of the young rapidly rotating K3 dwarf Speedy Mic (BO Mic, HD 197890). The photospheric spot maps reveal a heavily and uniformly spotted surface from equatorial to high latitude regions. Contrary to many images of similar objects, Speedy Mic does not possess a uniform filling at high latitudes, but exhibits structure in the polar regions showing greatest concentration in a particular longitude range. The asymmetric rotation profile of Speedy Mic indicates the presence of a companion or nearby star which shows radial velocity shifts over a timescale of several years. Using a simple dynamical argument, we show that Speedy Mic is unlikely to be a binary system, and conclude that the feature must be the result of a chance alignment with a background binary. Complete phase coverage on two consecutive nights in addition to 60 per cent phase coverage after a three night gap has enabled us to track the evolution of spots with time. By incorporating a solar-like differential rotation model into the image reconstruction process, we find that the equator laps the polar regions once every $191 \pm 17$ d. This finding is in close agreement with measurements for other late-type rapid rotators.

We have used spectrophotometric data from the Hubble Space Telescope to eclipse-map the primary component of the RS CVn binary SV Cam over 9 HST orbits. We find from these observations and the HIPPARCOS parallax that the surface flux in the eclipsed low-latitude region of the primary is about 30% lower than computed from a PHOENIX model atmosphere at the effective temperature that best fits the spectral energy distribution of the eclipsed flux. This can only be accounted for if about a third of the primary's surface is covered with unresolved dark star-spots. Extending this to the full surface of the primary, we find that even taking into account this spot filling factor there is an additional flux deficit on the primary star. This can only be explained if there is a large polar spot on the primary star extending from the pole to latitude 48 degrees.

We use Doppler imaging techniques to determine the dependence of starspot rotation rates on latitude in an homogeneous sample of young, rapidly-rotating solar analogues. A solar-like differential rotation law is used, where the rotation depends on sin$^2$($\theta$), where $\theta$ is the stellar latitude. By including this term in the image reconstruction process, using starspots as tracers, we are able to determine the magnitude of the shear over more than one rotation cycle. We also consider results from matched filter starspot tracking techniques, where individual starspot rotation rates are determined. In addition we have re-analysed published results and present a new measurement for the K3 dwarf, Speedy Mic. A total of 10 stars of spectral type G2 - M2 are considered. We find a trend towards decreasing surface differential rotation with decreasing effective temperature. The implied approach to solid body rotation with increasing relative convection zone depth implies that the dynamo mechanism operating in low-mass stars may be substantially different from that in the Sun.

We present Doppler images of the young K5V - K7V rapid rotator LO Peg from seven nights of continuous spectroscopy obtained in 1998 from July 04 to July 10. The images reveal the presence of a strong polar cap with appendages extending to mid-latitudes, but no starspots are seen below 15\degr. We briefly discuss the distribution of spots in light of recent flux transport simulations which are able to reproduce the observed latitude dependence. With the full timeseries of spectra, of which 314 are useful, many phases are observed three times over the seven nights of observations. Using starspots as tracers of a solar-like latitudinal differential rotation in our image reconstructions, we find that the equatorial regions complete one more rotation than the polar regions every {181 $\pm$ 35} d. LO Peg is the second coolest star for which such a measurement has been made using indirect imaging methods. The degree of latitudinal shear is less than that seen in G and early K dwarfs, suggesting a trend in which differential rotation decreases with stellar mass in (pre-)main sequence objects.

We investigate whether rings around extrasolar planets could be detected from those planets' transit lightcurves. To this end we develop a basic theoretical framework for calculating and interpreting the lightcurves of ringed planet transits based on the existing framework used for stellar occultations, a technique which has been effective for discovering and probing ring systems in the solar system. We find that the detectability of large, Saturn-like ring systems is largest during ingress and egress, and that reasonable photometric precisions of $\sim 1-3 \times 10^{-4}$ with 15-minute time resolution should be sufficient to discover such ring systems. For some ring particle sizes, diffraction around individual particles leads to a detectable level of forward-scattering that can be used to measure modal ring particle diameters. An initial census of large ring systems can be carried out using high-precision follow-up observations of detected transits and by the upcoming NASA \emph{Kepler} mission. The distribution of ring systems as a function of stellar age and as a function of planetary semimajor axis will provide empirical evidence to help constrain how rings form and how long rings last.

Star formation plays an important role in the fate of interacting galaxies. To date, most galactic simulations including star formation have used a density-dependent star formation rule designed to approximate a Schmidt law. Here, I present a new star formation rule which is governed by the local rate of energy dissipation in shocks. The new and old rules are compared using self-consistent simulations of NGC 4676; shock-induced star formation provides a better match to the observations of this system.

We present results from an HST study of the morphology and kinematics of NGC 6240. This merging galaxy with a double nucleus is one of the nearest and best-studied ultraluminous infrared galaxies. HST resolves both nuclei into seperate components. The distance between the northern and southern optical/near-infrared components is greater than that observed in radio and X-ray studies, arguing that even in K-band we may not be seeing all the way through the dust to the true nuclei. The ionized gas does not display rotation around either of the nuclei, or equilibrium motion in general. There is a strong velocity gradient between the nuclei, similar to what is seen in CO data. There is no such gradient in our stellar kinematics. The velocity dispersion of the gas is larger than expected for a cold disk. We also map and model the emission-line velocity field at an off-nuclear position where a steep velocity gradient was previously detected in ground-based data. Overall, the data indicate that line-of-sight projection effects, dust absorption, non-equilibrium merger dynamics, and the possible influence of a wind may be playing an important role in the observed kinematics. Chandra observations of hard X-rays have shown that both of the nuclei contain an Active Galactic Nucleus (AGN). The HST data show no clear sign of the two AGNs: neither continuum nor narrow-band imaging shows evidence for unresolved components in the nuclei, and there are no increased emission line widths or rapid rotation near the nuclei. This underscores the importance of X-ray data for identifying AGNs in highly dust-enshrouded environments.

We present Doppler Images of the RS CVn binary system HR 1099 (V711 Tau) from spectra taken in two different sites, KPNO and MSO, during the MUSICO S 1998 campaign. Contemporaneous APT photometry is used to constrain the Doppler Images. The resulting maximum entropy reconstructions based on the least-square s deconvolved profiles, derived from $\sim$2000 photospheric absorption lines, r eveal the presence of starspots at medium-high latitudes. We have obtained maps for both components of the binary system for the first time. The predominant str ucture in the primary component is an off-centered polar spot, confirming previo us works on the same target by using independent codes. The result is verified b y using both data sets independently. The lower spectral resolution data set giv es a less detailed map for the MSO data set. The images obtained for the seconda ry component show a low latitude spot around orbital phase 0.7. This spot seems to mirror the structure seen on the primary. It might suggest that tidal forces may influence the spot distribution on this binary system.

We investigate the prospects for characterizing extrasolar giant planets by measuring planetary oblateness from transit photometry and inferring planetary rotational periods. The rotation rates of planets in the solar system vary widely, reflecting the planets' diverse formational and evolutionary histories. A measured oblateness, assumed composition, and equation of state yields a rotation rate from the Darwin-Radau relation. The lightcurve of a transiting oblate planet should differ significantly from that of a spherical one with the same cross-sectional area under identical stellar and orbital conditions. However, if the stellar and orbital parameters are not known a priori, fitting for them allows changes in the stellar radius, planetary radius, impact parameter, and stellar limb darkening parameters to mimic the transit signature of an oblate planet, diminishing the oblateness signature. Thus even if HD209458b had an oblateness of 0.1 instead of our predicted 0.003, it would introduce a detectable departure from a model spherical lightcurve at the level of only one part in 10^5. Planets with nonzero obliquity break this degeneracy because their ingress lightcurve is asymmetric relative to that from egress, and their best-case detectability is of order 10^{-4}. However, the measured rotation rate for these objects is non-unique due to degeneracy between obliquity and oblateness and the unknown component of obliquity along the line of sight. Detectability of oblateness is maximized for planets transiting near an impact parameter of 0.7 regardless of obliquity. Future measurements of oblateness will be challenging because the signal is near the photometric limits of current hardware and inherent stellar noise levels.

We have obtained echelle spectroscopy of 14 Population II objects selected from those previously observed by Bonifacio & Molaro (1997). For one object, HD 140283, we obtained exquisite data with the High Dispersion Spectrograph on the Subaru Telescope, with S/N exceeding 1000 per 0.018A pixel. Li abundances have been determined by spectral synthesis from both the 6708A resonance line and also from 6104A subordinate feature. Firm detections of the weak line have been made in seven objects, and upper limits are reported for the remainder. Our 6708A abundances agree with those reported by Bonifacio & Molaro at the 99%-confidence level. Abundances from the 6104A line hint at a higher Li abundance than that determined from the resonance feature, but this evidence is mixed; the weakness of the 6104A line and the large number of upper limits make it difficult to draw firm conclusions. NLTE-corrections increase (rather than eliminate) the size of the (potential) discrepancy, and binarity appears unlikely to affect any abundance difference. The effect of multi-dimensional atmospheres on the line abundances was also considered, although it appears that use of 3-D (LTE) models could again act to increase the discrepancy, if one is indeed present.

We investigate the long-term dynamical stability of hypothetical moons orbiting extrasolar giant planets. Stellar tides brake a planet's rotation and, together with tidal migration, act to remove satellites; this process limits the lifetimes of larger moons in extrasolar planetary systems. Because more massive satellites are removed more quickly than less massive ones, we are able to derive an upper mass limit for those satellites that might have survived to the present day. For example, we estimate that no primordial satellites with masses greater than 7 x 10^{-7} M_Earth$ (~ \70km radius for rho=3 g cm^{-3}) could have survived around the transiting planet HD209458b for the age of the system. No meaningful mass limits can be placed on moons orbiting Jovian planets more than ~0.6 AU from their parent stars. Earth-like moons of Jovian planets could exist for 5 Gyr in systems where the stellar mass is greater than 0.15 M_Sun. Transits show the most promise for the discovery of extrasolar moons -- we discuss prospects for satellite detection via transits using space-based photometric surveys and the limits on the planetary tidal dissipation factor Q that a discovery would imply.

LPL dominance over Steward is demonstrated through several lines of evidence including observations and modelling of previous April Fools day behavior. The decrease in Steward coolness is attributed to the departure of interesting graduate students from Steward that enrolled at LPL.

Observations indicate that much of the interstellar gas in merging galaxies may settle into extended gaseous disks. Here, I present simulations of disk formation in mergers of gas-rich galaxies. Up to half of the total gas settles into embedded disks; the most massive instances result from encounters in which both galaxies are inclined to the orbital plane. These disks are often warped, many have rather complex kinematics, and roughly a quarter have counter-rotating or otherwise decoupled central components. Disks typically grow from the inside out; infall from tidal tails may continue disk formation over long periods of time.

We present new VLA C+D-array HI observations and optical and NIR imaging of the well known interacting system NGC 4038/9, ``The Antennae''. The radio data reveal a wealth of gaseous sub-structure both within the main bodies of the galaxies and along the tidal tails. In agreement with previous HI studies, we find that the northern tail has HI along its outer length, but none along its base. We suggest that the HI at the base of this tail has been ionized by massive stars in the disk of NGC 4038. The gas in the southern tail has a bifurcated structure, with one filament lying along the optical tail and another running parallel to it but with no optical counterpart. The two filaments join just before the location of several star forming regions near the end of the tail. The HI velocity field at the end of the tail is dominated by strong velocity gradients which suggest that at this location the tail is bending away from us. We delineate and examine two regions within the tail previously identified as possible sites of a so-called ``tidal dwarf galaxy'' condensing out of the expanding tidal material. The tail velocity gradients mask any clear kinematic signature of a self-gravitating condensation in this region. A dynamical analysis suggest that there is not enough mass in gas alone for either of these regions to be self-gravitating. Conversely, if they are bound they require a significant contribution to their dynamical mass from evolved stars or dark matter. (Abridged)

WFPC2 images and STIS spectroscopic observations are presented of the double nucleus in the merger system NGC 6240. We find that: (a) the kinematics of the ionized gas is similar to that of the molecular gas, despite a different morphology; (b) the gaseous and stellar kinematics are quite different, suggesting an early merger stage; (c) neither the gaseous nor the stellar kinematics show an obvious sign of the supermassive black hole believed to be responsible for the X-ray emission of NGC 6240; and (d) the steep off-nuclear velocity gradient is not due to a 10E11 solar mass black hole, in contrast to earlier suggestions.

We present ground-based near-infrared (H-band) imaging of the circumstellar disk around the nearby classical T Tauri star TW Hydrae. The scattered light image shows a face-on disk with radius 4 arcseconds (corresponding to 225 AU) and a morphology which agrees with recent images from the Hubble Space Telescope and the Very Large Array. The best fit power-law for the disk's radial surface brightness profile obeys the law r^{-3.3\pm0.3}. We use our image and published continuum flux densities to derive properties of the disk with a simple model of emission from an optically thick, flat disk. The best-fit values for disk mass and inner radius are 0.03 M_{\odot} and 0.3 AU; best-fit values for temperature, density, and grain opacity power law exponents (q, p, and \beta) are 0.7, 1.3, and 0.9, respectively. These properties are similar to those of disks around classical T Tauri stars located in more distant molecular clouds. Because of TW Hydrae's nearby location and pole-on orientation, it is a uniquely favorable object for future studies of radial disk structure at the classical T Tauri stage.

Lithium abundances, radial velocities and rotational velocities are reported for 20 candidate F, G and K stars of the sparse Coma Berenices open cluster. All the stars are proper-motion selected, and our radial velocities support the membership credentials of at least 12 of the candidates. Combining our data with that in the literature, we have a close-to-complete census of Li abundances for late-type stars with Mv < 5.8. These data show that the Li-depletion pattern in Coma Ber is similar, but not identical, to that in the Hyades cluster, which has a similar age but higher metallicity. Several Coma Ber F stars have suffered significantly more Li depletion than their counterparts in the Hyades. The G and early-K stars of Coma Ber have undergone less Li depletion than those of the Hyades, but much more than that predicted by standard evolutionary models featuring only convective mixing. This provides strong evidence for additional mixing and Li depletion operating in these stars during their first 400Myr on the main sequence, amounting to 0.3 dex at 6000K and rising to 0.8-1.2 dex at 5400K. We find that 4 of the radial-velocity non-members are among a small number of low-mass stars which were previously reported as part of an extra-tidal moving group associated with the Coma Ber cluster. As a result we now find that the luminosity function of this moving group is indistinguishable from that of the central cluster. It is uncertain whether there may be a significant number of stars with even lower masses among the moving group.

I compare gas-dynamical and stellar-dynamical models of collisions. These two models have distinctly different physics; for example, shocks introduce irreversibility in gas systems, while stellar systems evolve in a completely reversible fashion. Nonetheless, both models yield broadly similar results, suggesting that analogies between gas and stellar dynamics have some heuristic validity even applied to collisions.

In interacting and merging galaxies, gas is subject to direct hydrodynamic effects as well as tidal forces. One consequence of interactions is the rapid inflows of gas which may fuel starbursts and AGN. But gas dynamics is not limited to inflows; a small survey of equal-mass and unequal-mass encounters produces a wide variety of features, including plumes between galaxies, extended disks formed by infall of tidal debris, and counterrotating nuclear disks. An even richer spectrum of behavior awaits better thermodynamic models for gas in merging galaxies.

We describe a simple algorithm for classifying orbits into orbit families. This algorithm works by finding patterns in the sign changes of the principal coordinates. Orbits in the logarithmic potential are studied as an application; we classify orbits into boxlet families, and examine the influence of the core radius on the set of stable orbit families.

We use Gauss-Hermite functions to study the line of sight velocity distributions in simulated merger remnants. Our sample contains sixteen remnants; eight produced by mergers between disk galaxies of equal mass and eight produced by mergers between disk galaxies with mass ratios of 3:1. The equal-mass mergers display a wide range of kinematic features, including counterrotation at large radii, orthogonally rotating cores, and misaligned rotational axes. Most of the unequal-mass remnants exhibit fairly regular disk-like kinematics, although two have kinematics more typical of the equal-mass remnants. Our results may be compared to observations of early type objects, including ellipticals with misaligned kinematic axes, counterrotating systems, and S0 galaxies.

Theoretical considerations and observational data support the idea that mergers were more frequent in the past. At high redshifts, violent interactions and mergers may be implicated in the origin of Lyman-break galaxies, sub-mm starbursts, and active galactic nuclei. Most stars in cluster ellipticals probably formed at redshifts z > 2, as did most of the halo and globular clusters of the Milky Way; these events may all be connected with mergers. But what kind of galaxies merged at high redshifts, and how are these early events connected to present-epoch mergers? I will approach these questions by describing ideas for the formation of the Milky Way, elliptical galaxies, and populations of globular clusters.

Theoretical considerations and observational data support the idea that mergers were more frequent in the past. At redshifts z = 2 to 5, violent interactions and mergers may be implicated by observations of Lyman-break galaxies, sub-mm starbursts, and active galactic nuclei. Most stars in cluster ellipticals probably formed at such redshifts, as did most of the halo and globular clusters of the Milky Way; these events may all be connected with mergers. But what kind of galaxies merged at high redshifts, and are present-epoch mergers useful guides to these early collisions? I will approach these questions by describing ideas for the formation of the Milky Way, elliptical galaxies, and systems of globular clusters.

Interacting galaxies are a natural arena for studies of non-equilibrium stellar dynamics, gas dynamics, and thermodynamics. Only galaxy formation itself is as deeply concerned with as many different aspects of dynamics, and the connection between interactions and the formation of galaxies is probably no coincidence. This review discusses tidal interactions, halos and orbit decay, dissipative effects in galaxy interactions, properties of merger remnants, and origins of starbursts.

This review focus on some issues which seem relevant to recent discussions: (1) how halo structure influences tail length, (2) the fate of power-law density cusps, (3) the results of unequal-mass disk galaxy mergers, and (4) the behavior of hot and cold gas in merging disk galaxies.