We present C and X-band radio observations of the famous utraluminous X-ray source (ULX) Holmberg IX X-1, previously discovered to be associated with an optical emission line nebula several hundred pc in extent. Our recent infrared study of the ULX suggested that a jet could be responsible for the infrared excess detected at the ULX position. The new radio observations, performed using the Karl G. Jansky Very Large Array (VLA) in B-configuration, reveal the presence of a radio counterpart to the nebula with a spectral slope of -0.56 similar to other ULXs. Importantly, we find no evidence for an unresolved radio source associated with the ULX itself, and we set an upper limit on any 5 GHz radio core emission of 6.6 $\mu$Jy ($4.1\times10^{32}$ erg s$^{-1}$). This is 20 times fainter than what we expect if the bubble is energized by a jet. If a jet exists its core component is unlikely to be responsible for the infrared excess unless it is variable. Strong winds which are expected in super-Eddington sources could also play an important role in inflating the radio bubble. We discuss possible interpretations of the radio/optical bubble and we prefer the jet+winds-blown bubble scenario similar to the microquasar SS 433.

We present results of a systematic search for gravitationally lensed quasars in Pan-STARRS1. Our final sample of candidates comprises of 91 systems, not including 25 rediscovered lensed quasars and quasar pairs. In the absence of spectroscopy to verify the lensing nature of the candidates, the main sources of contaminants are likely to be quasar pairs, which we consider to be a byproduct of our work, and a smaller number of quasar$+$star associations. Amongst the independently discovered quads is 2M1134$-$2103, for which we obtained spectroscopy for the first time, finding a redshift of 2.77 for the quasar. There is evidence for microlensing in at least one image. We perform detailed mass modeling of this system using archival imaging data, and find that the unusually large shear responsible for the diamond-like configuration can be attributed mainly to a faint companion $\sim 4''$ away, and to a galaxy group/cluster $\sim 30''$ away. We also set limits of $z\sim0.5-1.5$ on the redshift of the lensing galaxy, based on its brightness, the image separation of the lensed images, and an analysis of the observed photometric flux ratios.

We investigate a sample of 3412 {\it International Celestial Reference Frame} (ICRF3) extragalactic radio-loud sources with accurate positions determined by VLBI in the S/X band, mostly active galactic nuclei (AGN) and quasars, which are cross-matched with optical sources in the second Gaia data release (Gaia DR2). The main goal of this study is to determine a core sample of astrometric objects that define the mutual orientation of the two fundamental reference frames, the Gaia (optical) and the ICRF3 (radio) frames. The distribution of normalized offsets between the VLBI sources and their optical counterparts is non-Rayleigh, with a deficit around the modal value and a tail extending beyond the 3-sigma confidence level. A few filters are applied to the sample in order to discard double cross-matches, confusion sources, and Gaia astrometric solutions of doubtful quality. {\it Panoramic Survey Telescope and Rapid Response System} (Pan-STARRS) and {\it Dark Energy Survey} (DES) stacked multi-color images are used to further deselect objects that are less suitable for precision astrometry, such as extended galaxies, double and multiple sources, and obvious misidentifications. After this cleaning, 2642 quasars remain, of which 20\% still have normalized offset magnitudes exceeding 3, or 99\% confidence level. We publish a list of 2118 radio-loud quasars of prime astrometric quality. The observed dependence of binned median offset on redshift shows the expected decline at small redshifts, but also an unexpected rise at $z\sim 1.6$, which may be attributed to the emergence of the C IV emission line in the Gaia's $G$ band. The Gaia DR2 parallax zero-point is found to be color-dependent, suggesting an uncorrected instrumental calibration effect.

We perform numerical simulations of the TRAPPIST-1 system of seven exoplanets orbiting a nearby M dwarf, starting with a previously suggested stable configuration. The long-term stability of this configuration is confirmed, but the motion of planets is found to be chaotic. The eccentricity values are found to vary within finite ranges. The rates of tidal dissipation and tidal evolution of orbits are estimated, assuming an Earth-like rheology for the planets. We find that under this assumption the planets b, d, e were captured in the 3:2 or higher spin-orbit resonances during the initial spin-down but slipped further down into the 1:1 resonance. Dependent on its rheology, the innermost planet b may be captured in a stable pseudosynchronous rotation. Non-synchronous rotation ensures higher levels of tidal dissipation and internal heating. The positive feedback between the viscosity and the dissipation rate -- and the ensuing runaway heating -- are terminated by a few self-regulation processes. When the temperature is high and the viscosity is low enough, the planet spontaneously leaves the 3:2 resonance. Further heating is stopped either by passing the peak dissipation or by the emergence of partial melt in the mantle. In the post-solidus state, the tidal dissipation is limited to the levels supported by the heat transfer efficiency. The tides on the host star are unlikely to have had a significant dynamical impact. The tides on the synchronized inner planets tend to reduce these planets' orbital eccentricity, possibly contributing thereby to the system's stability.

Highly volcanic exoplanets, which can be variously characterized as 'lava worlds', 'magma ocean worlds', or 'super-Ios' are high priority targets for investigation. The term 'lava world' may refer to any planet with extensive surface lava lakes, while the term 'magma ocean world' refers to planets with global or hemispherical magma oceans at their surface. 'Highly volcanic planets', including super-Ios, may simply have large, or large numbers of, active explosive or extrusive volcanoes of any form. They are plausibly highly diverse, with magmatic processes across a wide range of compositions, temperatures, activity rates, volcanic eruption styles, and background gravitational force magnitudes. Worlds in all these classes are likely to be the most characterizable rocky exoplanets in the near future due to observational advantages that stem from their preferential occurrence in short orbital periods and their bright day-side flux in the infrared. Transit techniques should enable a level of characterization of these worlds analogous to hot Jupiters. Understanding processes on highly volcanic worlds is critical to interpret imminent observations. The physical states of these worlds are likely to inform not just geodynamic processes, but also planet formation, and phenomena crucial to habitability. Volcanic and magmatic activity uniquely allows chemical investigation of otherwise spectroscopically inaccessible interior compositions. These worlds will be vital to assess the degree to which planetary interior element abundances compare to their stellar hosts, and may also offer pathways to study both the very young Earth, and the very early form of many silicate planets where magma oceans and surface lava lakes are expected to be more prevalent. We suggest that highly volcanic worlds may become second only to habitable worlds in terms of both scientific and public long-term interest.

The geophysics of extrasolar planets is a scientific topic often regarded as standing largely beyond the reach of near-term observations. This reality in no way diminishes the central role of geophysical phenomena in shaping planetary outcomes, from formation, to thermal and chemical evolution, to numerous issues of surface and near-surface habitability. We emphasize that for a balanced understanding of extrasolar planets, it is important to look beyond the natural biases of current observing tools, and actively seek unique pathways to understand exoplanet interiors as best as possible during the long interim prior to a time when internal components are more directly accessible. Such pathways include but are not limited to: (a) enhanced theoretical and numerical modeling, (b) laboratory research on critical material properties, (c) measurement of geophysical properties by indirect inference from imprints left on atmospheric and orbital properties, and (d) the purpose-driven use of Solar System object exploration expressly for its value in comparative planetology toward exoplanet-analogs. Breaking down barriers that envision local Solar System exploration, including the study of Earth's own deep interior, as separate from and in financial competition with extrasolar planet research, may greatly improve the rate of needed scientific progress for exoplanet geophysics. As the number of known rocky and icy exoplanets grows in the years ahead, we expect demand for expertise in 'exogeoscience' will expand at a commensurately intense pace. We highlight key topics, including: how water oceans below ice shells may dominate the total habitability of our galaxy by volume, how free-floating nomad planets may often attain habitable subsurface oceans supported by radionuclide decay, and how deep interiors may critically interact with atmospheric mass loss via dynamo-driven magnetic fields.

We report the serendipitous discovery of the first gravitationally lensed quasar candidate from Pan-STARRS. The grizy images reveal four point-like images with magnitudes between 14.9 mag and 18.1 mag. The colors of the point sources are similar, and they are more consistent with quasars than with stars or galaxies. The lensing galaxy is detected in the izy bands, with an inferred photometric redshift of ~0.6, lower than that of the point sources. We successfully model the system with a singular isothermal ellipsoid with shear, using the relative positions of the five objects as constraints. While the brightness ranking of the point sources is consistent with that of the model, we find discrepancies between the model-predicted and observed fluxes, likely due to microlensing by stars and millilensing due to the dark matter substructure. In order to fully confirm the gravitational lens nature of this system and add it to the small but growing number of the powerful probes of cosmology and astrophysics represented by quadruply lensed quasars, we require further spectroscopy and high-resolution imaging.

We investigate a sample of 2293 ICRF2 extragalactic radio-loud sources with accurate positions determined by VLBI, mostly active galactic nuclei (AGN) and quasars, which are cross-matched with optical sources in the first Gaia release (Gaia DR1). The distribution of offsets between the VLBI sources and their optical counterparts is strongly non-Gaussian, with powerful wings extending beyond 1 arcsecond. Limiting our analysis to only high-confidence difference detections, we find (and publish) a list of 188 objects with normalized variances above 12 and offsets below 1 arcsecond. Pan-STARRS stacked and monochromatic images resolve some of these sources indicating the presence of double sources, confusion sources, or pronounced extended structures. Some 89 high-quality objects, however, do not show any perturbations and appear to be star-like single sources, yet displaced by multiples of the expected error from the radio-loud AGN. We conclude that a fraction of luminous AGN (more than 4%) can be physically dislodged from the optical centers of their parent galaxies.

We present Spitzer Infrared Array Camera (IRAC) photometric observations of the Ultra-luminous X-ray Source (ULX, X-1) in Holmberg IX. We construct a spectral energy distribution (SED) for Holmberg IX X-1 based on published optical, UV and X-ray data combined with the IR data from this analysis. We modeled the X-ray and optical data with disk and stellar models, however we find a clear IR excess in the ULX SED that cannot be explained by fits or extrapolations of any of these models. Instead, further analysis suggests that the IR excess results either from dust emission, possibly from a circumbinary disk or from a variable jet.

We describe development and application of a Global Astrometric Solution (GAS) to the problem of Pan-STARRS1 (PS1) astrometry. Current PS1 astrometry is based on differential astrometric measurements using 2MASS reference stars, thus PS1 astrometry inherits the errors of the 2MASS catalog. The GAS, based on a single, least squares adjustment to approximately 750k grid stars using over 3000 extragalactic objects as reference objects, avoids this catalog-to-catalog propagation of errors to a great extent. The GAS uses a relatively small number of Quasi-Stellar Objects (QSOs, or distant AGN) with very accurate (<1 mas) radio positions, referenced to the ICRF2. These QSOs provide a hard constraint in the global least squares adjustment. Solving such a system provides absolute astrometry for all the stars simultaneously. The concept is much cleaner than conventional astrometry but is not easy to perform for large catalogs. In this paper we describe our method and its application to Pan-STARRS1 data. We show that large-scale systematic errors are easily corrected but our solution residuals for position (~60 mas) are still larger than expected based on simulations (~10 mas). We provide a likely explanation for the reason the small-scale residual errors are not corrected in our solution as would be expected.

We have used XMM-Newton's Optical Monitor (OM) images to study the local environment of a sample of 27 Ultraluminous X-ray Sources (ULXs) in nearby galaxies. UVW1 fluxes were extracted from 100 pc regions centered on the ULX positions. We find that at least 4 ULXs (out of 10 published) have spectral types that are consistent with previous literature values. In addition the colors are similar to those of young stars. For the highest-luminosity ULXs, the UVW1 fluxes may have an important contribution from the accretion disk. We find that the majority of ULXs are associated with recent star-formation. Many of the ULXs in our sample are located inside young OB associations or star-forming regions (SFRs). Based on their colors, we estimated ages and masses for star-forming regions located within 1 kpc from the ULXs in our sample. The resolution of the OM was insufficient to detect young dense super-clusters, but some of these star-forming regions are massive enough to contain such clusters. Only three ULXs have no associated SFRs younger than 50 Myr. The age and mass estimates for clusters were used to test runaway scenarios. The data are in general compatible with stellar-mass binaries accreting at super-Eddington rates and ejected by natal kicks. We also tested the hypothesis that ULXs are sub-Eddington accreting IMBHs ejected by three-body interactions, however this is not supported well by the data.

We investigate the dynamical evolution of the potentially habitable super-earth GJ 667Cc in the multiple system of at least two exoplanets orbiting a nearby M dwarf, paying special attention to its spin-orbital state. The published radial velocities for this star are re-analyzed and evidence is found for additional periodic signals, which could be taken for two additional planets on eccentric orbits. Such systems are not dynamically viable and break up quickly in numerical integrations. Limiting the scope to the two originally detected planets, we assess the dynamical stability of the system and find no evidence for bounded chaos in the orbital motion. The orbital eccentricity of the planets b and c is found to change cyclicly in the range 0.06 - 0.28 and 0.05 - 0.25, respectively, with a period of approximately 0.46 yr. Taking the eccentricity variation into account, numerical integrations are performed of the differential equations modeling the spin-orbit interaction of the planet GJ 667Cc with its host star assuming a terrestrial composition of its mantle. Depending on the interior temperature of the planet, it is likely to be entrapped in the 3:2 (probability 0.51) or even higher spin-orbit resonance. It is less likely to reach the 1:1 resonance (probability 0.24). The estimated characteristic spin-down times are quite short for the two planets, i.e., within 1 Myr for planet c and even shorter for planet b. Both planet arrived at their current and, most likely, ultimate spin-orbit states a long time ago. The rate of tidal dissipation of energy is formidably high in the planets of GJ 667, estimated at $10^{23.7}$ and $10^{26.7}$ J yr$^{-1}$ for c and b, respectively. This raises a question of how such relatively massive, close super-Earths could survive overheating and destruction.

GJ 581d is a potentially habitable super-Earth in the multiple system of exoplanets orbiting a nearby M dwarf. We investigate this planet's long-term dynamics, with an emphasis on its probable final rotation states acquired via tidal interaction with the host. The published radial velocities for the star are re-analysed with a benchmark planet detection algorithm, to confirm that there is no evidence for the recently proposed two additional planets (f and g). Limiting the scope to the four originally detected planets, we assess the dynamical stability of the system and find bounded chaos in the orbital motion. For the planet d, the characteristic Lyapunov time is 38 yr. Long-term numerical integration reveals that the system of four planets is stable, with the eccentricity of the planet d changing quasi-periodically in a tight range around 0.27, and with its semimajor axis varying only a little. The spin-orbit interaction of GJ 581d with its host star is dominated by the tides exerted by the star on the planet. We model this interaction, assuming a terrestrial composition of the mantle. Besides the customarily included secular parts of the triaxiality-caused and tidal torques, we also include these torques' oscillating components. It turns out that, dependent on the mantle temperature, the planet gets trapped into the 2:1 or an even higher spin-orbit resonance. It is very improbable that the planet could have reached the 1:1 resonance. This enhances the possibility of the planet being suitable for sustained life.

We propose a computer-efficient and accurate method of estimation of spatially correlated errors in astrometric positions, parallaxes and proper motions obtained by space and ground-based astrometry missions. In our method, the simulated observational equations are set up and solved for the coefficients of scalar and vector spherical harmonics representing the output errors, rather than for individual objects in the output catalog. Both accidental and systematic correlated errors of astrometric parameters can be accurately estimated. The method is demonstrated on the example of the JMAPS mission, but can be used for other projects of space astrometry, such as SIM or JASMINE.

We present Spitzer Infrared Spectrograph (IRS) observations of the ultra-luminous X-ray source (ULX) NGC 6946 X-1 and its associated nebula MF 16. This ULX has very similar properties to the famous Holmberg II ULX, the first ULX to show a prominent infrared [O IV] emission line comparable to those found in AGN. This paper attempts to constrain the ULX Spectral Energy Distribution (SED) given the optical/UV photometric fluxes and high-resolution X-ray observations. Specifically, Chandra X-ray data and published Hubble optical/UV data are extrapolated to produce a model for the full optical to X-ray SED. The photoionization modeling of the IR lines and ratios is then used to test different accretion spectral models. While either an irradiated disk model or an O-supergiant plus accretion disk model fit the data very well, we prefer the latter because it fits the nebular parameters slightly better. In this second case the accretion disk alone dominates the extreme-UV and X-ray emission, while an O-supergiant is responsible for most of the far-UV emission.

Quasometry is precision measurement of celestial positions and apparent motion of very distant extragalactic objects, such as quasars, galactic nuclei, and QSOs. We use this term to identify a specific area of research, the methodology of which differs from that of general astrometry. The main purpose of quasometry is to link the sub-milliarcsecond radio frame (ICRF) with the existing and emerging optical reference frames of similar accuracy, constructed by astrometric satellites. Some of the main difficulties in achieving this goal are discussed, e.g., the extended structures of quasar hosts, apparent motion on the sky, optical variability, galactic companions, faintness. Besides the strategic purpose, quasometry is undoubtedly useful for global astrometric surveys, as it helps to verify or even correct the resulting reference frames. There are two options of using measurements of distant quasars in a global astrometric solution: 1) hard constraints embedded in the fabric of observational equations; 2) {\it a posteriori} fitting of zonal errors. The relative benefits and shortcoming of the two options are reviewed. A relatively small set of about 200 carefully selected reference quasars can go a long way in improving the astrometric value of a space mission, if they are sufficiently bright, stable, fairly uniformly distributed on the sky, and are defining sources in the ICRF. We present an ongoing program at the USNO to construct a quality set of optical quasars with the required properties and to enhance the ICRF with new sources in the areas where known, well-observed quasars are scarce.

We compare mid-infrared emission-line properties, from high-resolution Spitzer spectra of a hard X-ray (14 -- 195 keV) selected sample of nearby (z < 0.05) AGN detected by the Burst Alert Telescope (BAT) aboard Swift. The luminosity distribution for the mid-infrared emission-lines, [O IV] 25.89 micron, [Ne II] 12.81 micron, [Ne III] 15.56 micron and [Ne V] 14.32/24.32 micron, and hard X-ray continuum show no differences between Seyfert 1 and Seyfert 2 populations, however six newly discovered BAT AGNs are under-luminous in [O IV], most likely the result of dust extinction in the host galaxy. The overall tightness of the mid-infrared correlations and BAT fluxes and luminosities suggests that the emission lines primarily arise in gas ionized by the AGN. We also compare the mid-infrared emission-lines in the BAT AGNs with those from published studies of ULIRGs, PG QSOs, star-forming galaxies and LINERs. We find that the BAT AGN sample fall into a distinctive region when comparing the [Ne III]/[Ne II] and the [O IV]/[Ne III] ratios. These line ratios are lower in sources that have been previously classified in the mid-infrared/optical as AGN than those found for the BAT AGN, suggesting that, in our X-ray selected sample, the AGN represents the main contribution to the observed line emission. These ratios represent a new emission line diagnostic for distinguishing between AGN and star forming galaxies.

This is the second of two papers examining Spitzer Infrared Spectrograph (IRS) observations of the ultraluminous X-ray source (ULX) in Holmberg II. Here we perform detailed photoionization modeling of the infrared lines. Our analysis suggests that the luminosity and morphology of the [O IV] 25.89 $\mu$m emission line is consistent with photoionization by the soft X-ray and far ultraviolet (FUV) radiation from the accretion disk of the binary system and inconsistent with narrow beaming. We show that the emission nebula is matter-bounded both in the line of sight direction and to the east, and probably radiation-bounded to the west. A bolometric luminosity in excess of 10$^{40}$ erg s$^{-1}$ would be needed to produce the measured [O IV] flux. We use modeling and previously published studies to conclude that shocks likely contribute very little, if at all, to the high-ionization line fluxes observed in the Holmberg II ULX. Additionally, we find that the spectral type of the companion star has a surprisingly strong effect on the predicted strength of the [O IV] emission. This finding could explain the origin of [O IV] in some starburst systems containing black hole binaries.

We present the first Spitzer Infrared Spectrograph (IRS) observations of the [O IV] 25.89 $\mu$m emission line detected from the ultraluminous X-ray source (ULX) in Holmberg II. This line is a well established signature of high ionization, usually associated with AGN. Its detection suggests that the ULX has a strong impact on the surrounding gas. A Spitzer high resolution spectral map shows that the [O IV] is coincident with the X-ray position of the ULX. The ratios of the [O IV] to lower ionization lines are similar to those observed in AGN, suggesting that a strong UV and X-ray source is responsible for the photoionization. The best XMM-Newton data is used to model the X-ray band which is then extrapolated into the UV. We perform infrared and ultraviolet photometry, and use previously published optical and radio data to construct the full SED for the ULX and its companion. The preferred model to describe the SED includes an accretion disk which dominates the soft X-rays but contributes little at UV and optical wavelengths. The optical counterpart is consistent with a B supergiant as previously suggested in other studies. The bolometric luminosity of the ULX suggests the presence of an intermediate-mass black hole with mass $>$85 M$_\odot$ for sub-Eddington accretion or, alternatively, a stellar-mass black hole that is accreting at super-Eddington rates. In a follow-up second paper we perform detailed photoionization modeling of the infrared lines in order to constrain the bolometric luminosity of the ULX.

To test the idea that ultraluminous X-ray sources (ULXs) in external galaxies represent a class of accreting intermediate-mass black holes (IMBHs), we have undertaken a program to identify ULXs and a lower luminosity X-ray comparison sample with the highest quality data in the {\it Chandra} archive. We establish as a general property of ULXs that the most X-ray-luminous objects possess the flattest X-ray spectra (in the {\it Chandra} bandpass). No prior sample studies have established the general hardening of ULX spectra with luminosity. This hardening occurs at the highest luminosities (absorbed luminosity $\geq5\times10^{39}$~erg~s$^{-1}$) and is in line with recent models arguing that ULXs are actually stellar-mass black holes. From spectral modeling, we show that the evidence originally taken to mean that ULXs are IMBHs - i.e., the "simple IMBH model" - is nowhere near as compelling when a large sample of ULXs is looked at properly. During the last couple of years, {\it XMM-Newton} spectroscopy of ULXs has to a large extent begun to negate the simple IMBH model based on fewer objects. We confirm and expand these results, which validates the {\it XMM-Newton} work in a broader sense with independent X-ray data. We find that (1) cool-disk components are present with roughly equal probability and total flux fraction for any given ULX, regardless of luminosity, and (2) cool-disk components extend below the standard ULX luminosity cutoff of 10$^{39}$~erg~s$^{-1}$, down to our sample limit of 10$^{38.3}$~erg~s$^{-1}$. The fact that cool disk components are not correlated with luminosity damages the argument that cool disks indicate IMBHs in ULXs, for which strong statistical support was never found.