Spectra of methyl cyanide were recorded to analyze interactions in low-lying vibrational states and to construct line lists for radio astronomical observations as well as for infrared spectroscopic investigations of planetary atmospheres. The rotational spectra cover large portions of the 36$-$1627 GHz region. In the infrared (IR), a spectrum was recorded for this study in the region of 2$\nu _8$ around 717 cm$^{-1}$ with assignments covering 684$-$765 cm$^{-1}$. Additional spectra in the $\nu _8$ region were used to validate the analysis. The large amount and the high accuracy of the rotational data extend to much higher $J$ and $K$ quantum numbers and allowed us to investigate for the first time in depth local interactions between these states which occur at high $K$ values. In particular, we have detected several interactions between $v_8 = 1$ and 2. Notably, there is a strong $\Delta v_8 = \pm1$, $\Delta K = 0$, $\Delta l = \pm3$ Fermi resonance between $v_8 = 1^{-1}$ and $v_8 = 2^{+2}$ at $K$ = 14. Pronounced effects in the spectrum are also caused by resonant $\Delta v_8 = \pm1$, $\Delta K = \mp2$, $\Delta l = \pm1$ interactions between $v_8 = 1$ and 2. An equivalent resonant interaction occurs between $K$ = 14 of the ground vibrational state and $K$ = 12, $l = +1$ of $v_8 = 1$ for which we present the first detailed account. A preliminary account was given in an earlier study on the ground vibrational state. From data pertaining to $v_8 = 2$, we also investigated rotational interactions with $v_4 = 1$ as well as $\Delta v_8 = \pm1$, $\Delta K = 0$, $\Delta l = \pm3$ Fermi interactions between $v_8 = 2$ and 3. We have derived N$_2$- and self-broadening coefficients for the $\nu _8$, 2$\nu _8 - \nu _8$, and 2$\nu _8$ bands from previously determined nu4 values. Subsequently, we determined transition moments and intensities for the three IR bands.

(abridged) We report the discovery of three new transiting planets: WASP-85 A b, WASP-116 b, and WASP-149 b. WASP-85 b orbits its host star every 2.66 days, and has a mass of 1.25 M_Jup and a radius of 1.25 R_Jup. The host star is of G5 spectral type, with magnitude V = 11.2, and lies 141 pc distant. The system has a K-dwarf binary companion, WASP-85 B, at a separation of ~1.5". The close proximity of this companion leads to contamination of our photometry, decreasing the apparent transit depth that we account for during our analysis. Analysis of the Ca II H+K lines shows strong emission that implies that both binary components are strongly active. WASP-116 b is a warm, mildly inflated super-Saturn, with a mass of 0.59 M_Jup and a radius of 1.43 R_Jup. It was discovered orbiting a metal-poor ([Fe/H] = -0.28 dex), cool (T_eff = 5950 K) G0 dwarf every 6.61 days. WASP-149 b is a typical hot Jupiter, orbiting a G6 dwarf with a period of 1.33 days. The planet has a mass and radius of 1.05 M_Jup and 1.29 R_Jup, respectively. The stellar host has an effective temperature of T_eff = 5750 K and has a metallicity of [Fe/H] = 0.16 dex. WASP photometry of the system is contaminated by a nearby star; we therefore corrected the depth of the WASP transits using the measured dilution. WASP-149 lies inside the 'Neptune desert' identified in the planetary mass-period plane by Mazeh, Holczer & Faigler (2016). We model the modulation visible in the K2 lightcurve of WASP-85 using a simple three-spot model consisting of two large spots on WASP-85 A, and one large spot on WASP-85 B, finding rotation periods of 13.1+/-0.1 days for WASP-85 A and 7.5+/-0.03 days for WASP-85 B. We estimate stellar inclinations of I_A = 66.8+/-0.7 degrees and I_B = 39.7+/-0.2 degrees, and constrain the obliquity of WASP-85 A b to be psi<27 degrees. We therefore conclude that WASP-85 A b is very likely to be aligned.

The averaged resonant equations of motion for the planar circular restricted three-body problem are solved on the linearization basis taking into account also non-gravitational effects. The averaged resonant equations are derived from Lagrange's planetary equations with additional Gauss's terms caused by the non-gravitational effects. The time depending solution has the standard form with exponential, quadratic, linear and constant terms. The existence of a rotational symmetry in the action of the non-gravitational effects around the star determines the order of a characteristic equation of the linearized system. In the symmetrical case (order 3) the considered non-gravitational effects are the stellar electromagnetic radiation and the radial stellar wind (stellar radiation). In the asymmetrical case (order 4) the stellar radiation and interstellar gas flow are considered. It is investigated how well the linearization solution describes real solution obtained from an equation of motion by a comparison of the resonant libration frequency found analytically and numerically. It is found that from initial values of the evolving orbital parameters (semimajor axis, eccentricity, longitude of pericenter, and resonant angular variable) in the averaged phase space the linearization frequency depends most sensitively on the initial value of the resonant angular variable. For small libration amplitudes of the resonant angular variable the best match of the real libration frequency and the linearization frequency is located approximately at the solution of the resonant condition ($da / dt$ $=$ 0). If the initial averaged conditions are chosen close to the solution of resonant condition, then the linearization frequency for practically all simple oscillatory evolutions matches the real libration frequency and the linearization solution very well approximates the real evolution.

This paper deals with a new formulation of the creep tide theory (Ferraz-Mello, Cel. Mech. Dyn. Astron. {\bf 116}, 109, 2013 $-$ Paper I) and with the tidal dissipation predicted by the theory in the case of stiff bodies whose rotation is not synchronous but is oscillating around the synchronous state with a period equal to the orbital period. We show that the tidally forced libration influences the amount of energy dissipated in the body and the average perturbation of the orbital elements. This influence depends on the libration amplitude and is generally neglected in the study of planetary satellites. However, they may be responsible for a 27 percent increase in the dissipation of Enceladus. The relaxation factor necessary to explain the observed dissipation of Enceladus ($\gamma=1.2-3.8\times 10^{-7}\ {\rm s}^{-1}$) has the expected order of magnitude for planetary satellites and corresponds to the viscosity $0.6-1.9 \times 10^{14}$ Pa s, which is in reasonable agreement with the value recently estimated by Efroimsky (2018) ($0.24 \times 10^{14}$ Pa s) and with the value adopted by Roberts and Nimmo (2008) for the viscosity of the ice shell ($10^{13}-10^{14}$ Pa s). For comparison purposes, the results are extended also to the case of Mimas and are consistent with the negligible dissipation and the absence of observed tectonic activity. The corrections of some mistakes and typos of paper II (Ferraz-Mello, Cel. Mech. Dyn. Astron. {\bf 122}, 359, 2015) are included at the end of the paper.

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

We present radial-velocity measurement of eight stars observed with the HARPS Echelle spectrograph mounted on the 3.6-m telescope in La Silla (ESO, Chile). Data span more than ten years and highlight the long-term stability of the instrument. We search for potential planets orbiting HD20003, HD20781, HD21693, HD31527, HD45184, HD51608, HD134060 and HD136352 to increase the number of known planetary systems and thus better constrain exoplanet statistics. After a preliminary phase looking for signals using generalized Lomb-Scargle periodograms, we perform a careful analysis of all signals to separate \emph{bona-fide} planets from spurious signals induced by stellar activity and instrumental systematics. We finally secure the detection of all planets using the efficient MCMC available on the Data and Analysis Center for Exoplanets (DACE web-platform), using model comparison whenever necessary. In total, we report the detection of twenty new super-Earth to Neptune-mass planets, with minimum masses ranging from 2 to 30 M$_{\rm Earth}$, and periods ranging from 3 to 1300 days. By including CORALIE and HARPS measurements of HD20782 to the already published data, we also improve the characterization of the extremely eccentric Jupiter orbiting this host.

Over the course of 3.1 years we determined the position of Mars on 75 dates using a handheld cross staff and two to five bright reference stars of known right ascension and declination on each occasion. On average the observed positions of Mars are within +/- 11 arc minutes of the true ecliptic latitudes and ecliptic longitudes. After deriving the two dates of opposition to the Sun, we were able to carry out a two stage analytical experiment on Mars' orbit. From the 2015/2016 data we obtain a value of the eccentricity of 0.093 +/- 0.012. The 2015/2016 data can be fit reasonably well by adopting a circular orbit for the Earth, but the 2017/2018 data must be fit with an ellipse for Mars and an ellipse for the Earth. Applying the two ellipse model to the 2015/2016 data, we obtain an RMS error of the ecliptic longitudes of only +/- 7.5 arc minutes. While Kepler was able to derive the shape of Mars' orbit while relying on data of Tycho Brahe accurate to +/- 2 arc minutes, today we may assume an ellptical orbit, and we can show that much less accurate data are consistent with an ellipse having an eccentricity equal to the modern accepted value of 0.0934, within the errors.

The disruption of an icy moon in a collision with an interloping comet a few hundred million years ago is a simple way to create Saturn's rings. A ring parent moon with a mass comparable to Mimas could be trapped in mean motion resonance with Enceladus and Dione in an orbit near the current outer edge of the rings just beyond the Roche zone. I present collisional N-body simulations of cometary impacts that lead to the partial disruption of a differentiated moon with a rocky core and icy mantle. The core can survive largely intact while the debris from the mantle settles into a ring of predominantly ice particles straddling the orbital radius of the parent moon. The nascent ring spreads radially due to collisional viscosity while mass re-accretes onto the remnant rocky core to form a new moon that can be identified as Mimas. The icy debris that migrates into the Roche zone evolves into Saturn's ring system. Torques from tidal interaction with Saturn and resonant interactions with the rings push the recently formed Mimas outward to its current position on the same timescale of a few hundred million years. This scenario accounts for the high ice fraction observed in Saturn's rings and explains why the ring mass is comparable to the mass of Mimas. The prior existence of a ring parent moon in mean motion resonance results in a tidal heating rate for Enceladus in the recent past that is significantly larger than the current rate.

Lakes existed on Mars later than 3.6 billion years ago, according to sedimentary evidence for deltaic deposition. The observed fluvio-lacustrine deposits suggest that individual lake-forming climates persisted for at least several thousand years (assuming dilute flow). But the lake watersheds' little weathered soils indicate a largely dry climate history, with intermittent runoff events. Here we show that these observational constraints, while inconsistent with many previously-proposed triggers for lake-forming climates, are consistent with a methane burst scenario. In this scenario, chaotic transitions in mean obliquity drive latitudinal shifts in temperature and ice loading that destabilize methane clathrate. Using numerical simulations, we find that outgassed methane can build up to atmospheric levels sufficient for lake forming climates, for past clathrate hydrate stability zone occupancy fractions >0.04. Such occupancy fractions are consistent with methane production by water-rock reactions due to hydrothermal circulation on early Mars. We further estimate that photochemical destruction of atmospheric methane curtails the duration of individual lake-forming climates to less than a million years, consistent with observations. We conclude that methane bursts represent a potential pathway for intermittent excursions to a warm, wet climate state on early Mars.

Large alluvial fan deposits on Mars record relatively recent habitable surface conditions ($\lesssim$3.5 Ga, Late Hesperian - Amazonian). We find net sedimentation rate <(4-8) {\mu}m/yr in the alluvial-fan deposits, using the frequency of craters that are interbedded with alluvial-fan deposits as a fluvial-process chronometer. Considering only the observed interbedded craters sets a lower bound of >20 Myr on the total time interval spanned by alluvial-fan aggradation, >10^3-fold longer than previous lower limits. A more realistic approach that corrects for craters fully entombed in the fan deposits raises the lower bound to >(100-300) Myr. Several factors not included in our calculations would further increase the lower bound. The lower bound rules out fan-formation by a brief climate anomaly. Therefore, during the Late Hesperian - Amazonian on Mars, persistent or repeated processes permitted habitable surface conditions.

Recently, Squire & Hopkins (2017) showed any coupled dust-gas mixture is subject to a class of linear 'resonant drag instabilities' (RDI). These can drive large dust-to-gas ratio fluctuations even at arbitrarily small dust-to-gas mass ratios. Here, we identify and study both resonant and new non-resonant instabilities, in the simple case where the gas satisfies neutral hydrodynamics and supports acoustic waves ($\omega^{2}=c_{s}^{2}\,k^{2}$). The gas and dust are coupled via an arbitrary drag law and subject to external accelerations (e.g. gravity, radiation pressure). If there is any dust drift velocity, the system is unstable. The instabilities exist for {\em all} dust-to-gas ratios $\mu$ and their growth rates depend only weakly on $\mu$ around resonance, as $\sim\mu^{1/3}$ or $\sim \mu^{1/2}$ (depending on wavenumber). The behavior changes depending on whether the drift velocity is larger or smaller than the sound speed $c_{s}$. In the supersonic regime a 'resonant' instability appears with growth rate increasing without limit with wavenumber, even for vanishingly small $\mu$ and values of the coupling strength ('stopping time'). In the subsonic regime non-resonant instabilities always exist, but their growth rates no longer increase indefinitely towards small wavelengths. The dimensional scalings and qualitative behavior of the instability do not depend sensitively on the drag law or equation-of-state of the gas. The instabilities directly drive exponentially growing dust-to-gas-ratio fluctuations, which can be large even when the modes are otherwise weak. We discuss physical implications for cool-star winds, AGN-driven winds and torii, and starburst winds: the instabilities alter the character of these outflows and could drive clumping and/or turbulence in the dust and gas.

Direct imaging has paved the way for atmospheric characterization of young and self-luminous gas giants. Scattering in a horizontally-inhomogeneous atmosphere causes the disk-integrated polarization of the thermal radiation to be linearly polarized, possibly detectable with the newest generation of high-contrast imaging instruments. We aim to investigate the effect of latitudinal and longitudinal cloud variations, circumplanetary disks, atmospheric oblateness, and cloud particle properties on the integrated degree and direction of polarization in the near-infrared. We have developed a three-dimensional Monte Carlo radiative transfer code (ARTES) for scattered light simulations in (exo)planetary atmospheres. The code is applicable to calculations of reflected light and thermal radiation in a spherical grid with a parameterized distribution of gas, clouds, hazes, and circumplanetary material. The disk-integrated degree of polarization of a horizontally-inhomogeneous atmosphere is maximal when the planet is flattened, the optical thickness of the equatorial clouds is large compared to the polar clouds, and the clouds are located at high altitude. For a flattened planet, the integrated polarization can both increase or decrease with respect to a spherical planet which depends on the horizontal distribution and optical thickness of the clouds. The direction of polarization can be either parallel or perpendicular to the projected direction of the rotation axis when clouds are zonally distributed. Rayleigh scattering by submicron-sized cloud particles will maximize the polarimetric signal whereas the integrated degree of polarization is significantly reduced with micron-sized cloud particles as a result of forward scattering. The presence of a cold or hot circumplanetary disk may also produce a detectable degree of polarization ($\lesssim$1%) even with a uniform cloud layer in the atmosphere.

In the framework of the planar and circular restricted three-body problem, we consider an asteroid that orbits the Sun in quasi-satellite motion with a planet. A quasi-satellite trajectory is a heliocentric orbit in co-orbital resonance with the planet, characterized by a non zero eccentricity and a resonant angle that librates around zero. Likewise, in the rotating frame with the planet it describes the same trajectory as the one of a retrograde satellite even though the planet acts as a perturbator. In the last few years, the discoveries of asteroids in this type of motion made the term "quasi-satellite" more and more present in the literature. However, some authors rather use the term "retrograde satellite" when referring to this kind of motion in the studies of the restricted problem in the rotating frame. In this paper we intend to clarify the terminology to use, in order to bridge the gap between the perturbative co-orbital point of view and the more general approach in the rotating frame. Through a numerical exploration of the co-orbital phase space, we describe the quasi-satellite domain and highlight that it is not reachable by low eccentricities by averaging process. We will show that the quasi-satellite domain is effectively included in the domain of the retrograde satellites and neatly defined in terms of frequencies. Eventually, we highlight a remarkable high eccentric quasi-satellite orbit corresponding to a frozen ellipse in the heliocentric frame. We extend this result to the eccentric case (planet on an eccentric motion) and show that two families of frozen ellipses originate from this remarkable orbit.

Exoplanet transit spectroscopy enables the characterization of distant worlds, and will yield key results for NASA's James Webb Space Telescope. However, transit spectra models are often simplified, omitting potentially important processes like refraction and multiple scattering. While the former process has seen recent development, the effects of light multiple scattering on exoplanet transit spectra has received little attention. Here, we develop a detailed theory of exoplanet transit spectroscopy that extends to the full refracting and multiple scattering case. We explore the importance of scattering for planet-wide cloud layers, where the relevant parameters are the slant scattering optical depth, the scattering asymmetry parameter, and the angular size of the host star. The latter determines the size of the "target" for a photon that is back-mapped from an observer. We provide results that straightforwardly indicate the potential importance of multiple scattering for transit spectra. When the orbital distance is smaller than 10-20 times the stellar radius, multiple scattering effects for aerosols with asymmetry parameters larger than 0.8-0.9 can become significant. We provide examples of the impacts of cloud/haze multiple scattering on transit spectra of a hot Jupiter-like exoplanet. For cases with a forward and conservatively scattering cloud/haze, differences due to multiple scattering effects can exceed 200 ppm, but shrink to zero at wavelength ranges corresponding to strong gas absorption or when the slant optical depth of the cloud exceeds several tens. We conclude with a discussion of types of aerosols for which multiple scattering in transit spectra may be important.

One of the biggest successes of the Cassini mission is the detection of small moons (moonlets) embedded in Saturn's rings which cause S-shaped density structures in their close vicinity, called propellers (Spahn and Sremcevic 2000; Tiscareno et al. 2006; Sremcevic et al. 2007). Here, we present isothermal hydrodynamic simulations of moonlet-induced propellers in Saturn's A ring which denote a further development of the original model (Spahn and Sremcevic 2000). We find excellent agreement between these new hydrodynamic and corresponding N-body simulations. Furthermore, the hydrodynamic simulations confirm the predicted scaling laws (Spahn and Sremcevic 2000) and the analytical solution for the density in the propeller gaps (Sremcevic et al. 2002). Finally, this mean field approach allows us to simulate the pattern of the giant propeller Bleriot, which is too large to be modeled by direct N-body simulations. Our results are compared to two stellar occultation observations by the Cassini Ultraviolet Imaging Spectrometer (UVIS), that intersect the propeller Bleriot. Best fits to the UVIS optical depth profiles are achieved for a Hill radius of 590 m, which implies a moonlet diameter of about 860 m. Furthermore, the model favours a kinematic shear viscosity of the surrounding ring material of $\nu_0 = 340$ cm^2/s, a dispersion velocity in the range of 0.3 cm/s $< c_0 <$ 1.5 cm/s, and a fairly high bulk viscosity $7 < \xi_0/\nu_0 < 17$. These large transport values might be overestimated by our isothermal ring model and should be reviewed by an extended model including thermal fluctuations.

We derive a self-similar description for the 2D streamline topology and flow structure of an axi-symmetric, thermally driven wind originating from a disc in which the density is a power law function of radius. Our scale-free solution is strictly only valid in the absence of gravity or centrifugal support; comparison with 2D hydrodynamic simulations of winds from Keplerian discs however demonstrates that the scale-free solution is a good approximation also in the outer regions of such discs, and can provide a reasonable description even for launch radii well within the gravitational radius of the flow. Although other authors have considered the flow properties along streamlines whose geometry has been specified in advance, this is the first isothermal calculation in which the flow geometry and variation of flow variables along streamlines is determined self-consistently. It is found that the flow trajectory is very sensitive to the power-law index of radial density variation in the disc: the steeper the density gradient, the stronger is the curvature of streamlines close to the flow base that is required in order to maintain momentum balance perpendicular to the flow. Steeper disc density profiles are also associated with more rapid acceleration, and a faster fall-off of density, with height above the disc plane. The derivation of a set of simple governing equations for the flow structure of thermal winds from the outer regions of power law discs offers the possibility of deriving flow observables without having to resort to hydrodynamical simulation.

Using spatio-temporal spectra we show direct evidence of excitation of magnetosonic and Alfv\'en waves in three-dimensional compressible magnetohydrodynamic turbulence at small Mach numbers. For the plasma pressure dominated regime, or high $\beta$ regime (with $\beta$ the ratio between fluid and magnetic pressure), and for the magnetic pressure dominated regime, or low $\beta$ regime, we study magnetic field fluctuations parallel and perpendicular to a guide magnetic field $\textbf{B}_0$. In the low $\beta$ case we find excitation of compressible and incompressible fluctuations, with a transfer of energy towards Alfv\'enic modes and to a lesser extent towards magnetosonic modes. In particular, we find signatures of the presence of fast magnetosonic waves in a scenario compatible with that of weak turbulence. In the high $\beta$ case, fast and slow magnetosonic waves are present, with no clear trace of Alfv\'en waves, and a significant part of the energy is carried by two-dimensional turbulent eddies.

GravityCam is a new concept of ground-based imaging instrument capable of delivering significantly sharper images from the ground than is normally possible without adaptive optics. Advances in optical and near infrared imaging technologies allow images to be acquired at high speed without significant noise penalty. Aligning these images before they are combined can yield a 2.5 to 3 fold improvement in image resolution. By using arrays of such detectors, survey fields may be as wide as the telescope optics allows. Consequently, GravityCam enables both wide-field high-resolution imaging and high-speed photometry. We describe the instrument and detail its application to provide demographics of planets and satellites down to Lunar mass (or even below) across the Milky Way. GravityCam is also suited to improve the quality of weak shear studies of dark matter distribution in distant clusters of galaxies and multiwavelength follow-ups of background sources that are strongly lensed by galaxy clusters. The photometric data arising from an extensive microlensing survey will also be useful for asteroseismology studies, while GravityCam can be used to monitor fast multiwavelength flaring in accreting compact objects, and promises to generate a unique data set on the population of the Kuiper belt and possibly the Oort cloud.

Distinct gap features in the nearest protoplanetary disk, TW Hya (distance of 59.5$\pm$0.9 pc), may be signposts of ongoing planet formation. We performed long-exposure thermal infrared coronagraphic imaging observations to search for accreting planets especially within dust gaps previously detected in scattered light and submm-wave thermal emission. Three nights of observations with the Keck/NIRC2 vortex coronagraph in $L^\prime$ (3.4-4.1$\mu$m) did not reveal any statistically significant point sources. We thereby set strict upper limits on the masses of non-accreting planets. In the four most prominent disk gaps at 24, 41, 47, and 88 au, we obtain upper mass limits of 1.6-2.3, 1.1-1.6, 1.1-1.5, and 1.0-1.2 Jupiter masses ($M_J$) assuming an age range of 7-10 Myr for TW Hya. These limits correspond to the contrast at 95\% completeness (true positive fraction of 0.95) with a 1\% chance of a false positive within $1^{\prime\prime}$ of the star. We also approximate an upper limit on the product of planet mass and planetary accretion rate of $M_p\dot{M}\lesssim10^{-8} M_J^2/yr$ implying that any putative $\sim0.1 M_J$ planet, which could be responsible for opening the 24 au gap, is presently accreting at rates insufficient to build up a Jupiter mass within TW Hya's pre-main sequence lifetime.

There have been many efforts to correct systematic effects in astronomical light curves to improve the detection and characterization of planetary transits and astrophysical variability. Algorithms like the Trend Filtering Algorithm (TFA) use simultaneously-observed stars to remove systematic effects, and binning is used to reduce high-frequency random noise. We present TFAW, a wavelet-based modified version of TFA. TFAW aims to increase the periodic signal detection and to return a detrended and denoised signal without modifying its intrinsic characteristics. We modify TFA's frequency analysis step adding a Stationary Wavelet Transform filter to perform an initial noise and outlier removal and increase the detection of variable signals. A wavelet filter is added to TFA's signal reconstruction to perform an adaptive characterization of the noise- and trend-free signal and the noise contribution at each iteration while preserving astrophysical signals. We carried out tests over simulated sinusoidal and transit-like signals to assess the effectiveness of the method and applied TFAW to real light curves from TFRM. We also studied TFAW's application to simulated multiperiodic signals, improving their characterization. TFAW improves the signal detection rate by increasing the signal detection efficiency (SDE) up to a factor ~2.5x for low SNR light curves. For simulated transits, the transit detection rate improves by a factor ~2-5x in the low-SNR regime compared to TFA. TFAW signal approximation performs up to a factor ~2x better than bin averaging for planetary transits. The standard deviations of simulated and real TFAW light curves are ~40x better than TFA. TFAW yields better MCMC posterior distributions and returns lower uncertainties, less biased transit parameters and narrower (~10x) credibility intervals for simulated transits. We present a newly-discovered variable star from TFRM.

Previous estimates of the annual mean surface temperature of Jupiter's moon, Europa, neglected the effect of the eccentricity of Jupiter's orbit around the Sun, the effect of the emissivity and heat capacity of Europa's ice, the effect of the eclipse of Europa (i.e., the relative time that Europa is within the shadow of Jupiter), the effect of Jupiter's radiation, and the effect of Europa's internal heating. Other studies concentrated on the diurnal cycle but neglected some of the above factors. In addition, to our knowledge, the seasonal cycle of the surface temperature of Europa was not estimated. Here we systematically estimate the diurnal, seasonal and annual mean surface temperature of Europa, when Europa's obliquity, emissivity, heat capacity, and eclipse, as well as Jupiter's radiation, internal heating, and eccentricity, are all taken into account. For a typical internal heating rate of 0.05 W m$^{-2}$, the equator, pole, and the global and mean annual mean surface temperatures are 96K, 46K, and 90K, respectively. We found that the temperature at the high latitudes is significantly affected by the internal heating, especially during the winter solstice, suggesting that measurements of high latitude surface temperatures can be used to constrain the internal heating. We also estimate the incoming solar radiation to Enceladus, the moon of Saturn.

The Transiting Exoplanet Survey Satellite (TESS) will be conducting a nearly all-sky photometric survey over two years, with a core mission goal to discover small transiting exoplanets orbiting nearby bright stars. It will obtain 30-minute cadence observations of all objects in the TESS fields of view, along with 2-minute cadence observations of 200,000 to 400,000 selected stars. The choice of which stars to observe at the 2-min cadence is driven by the need to detect small transiting planets, which leads to the selection of primarily bright, cool dwarfs. We describe the catalogs assembled and the algorithms used to populate the TESS Input Catalog (TIC). We also describe a ranking system for prioritizing stars according to the smallest transiting planet detectable, and assemble a Candidate Target List (CTL) using that ranking. We discuss additional factors that affect the ability to photometrically detect and dynamically confirm small planets, and we note additional stellar populations of interest that may be added to the final target list. The TIC is available on the STScI MAST server, and an enhanced CTL is available through the Filtergraph data visualization portal system at the URL https://filtergraph.vanderbilt.edu/tess_ctl .

We present observations of WASP-63b by the Hubble Space Telescope (HST) as part of "A Preparatory Program to Identify the Single Best Transiting Exoplanet for JWST Early Release Science". WASP-63b is one of the community targets under consideration for the James Webb Space Telescope (JWST) Early Release Science (ERS) program. We present a spectrum derived from a single observation by HST Wide Field Camera 3 in the near infrared. We engaged groups across the transiting exoplanet community to participate in the analysis of the data and present results from each. There is general agreement amongst all results that we find an H2O absorption feature with 3.5-4.0 sigma significance. However, the feature is muted in comparison to a clear atmosphere at solar composition. Although the detection of the water feature is robust, the reasons for the muting of this feature are ambiguous due to a degeneracy between clouds and composition. The data does not yield robust detections of any molecular species other than H2O. The group was motivated to perform an additional set of retrieval exercises to investigate an apparent bump in the spectrum at ~ 1.55 um. We explore possible disequilibrium chemistry and find this feature is consistent with super-solar HCN abundance but it is questionable if the required mixing ratio of HCN is chemically and physically plausible. The ultimate goal of this study is to vet WASP-63b as a potential community target to best demonstrate the capabilities and systematics of JWST instruments for transiting exoplanet science. In the case of WASP-63b, the presence of a detectable water feature indicates that WASP-63b remains a plausible target for ERS observations.

This work is a theoretical exploration for facilitating the interpretation of polarimetric observations in terms of cloudiness, rotational velocities and effective temperatures of brown dwarfs (BDs). An envelope of scatterers like free-electrons, atoms/molecules, or haze/clouds affects the Stokes-vector of radiation emitted by oblate bodies. Due to high rotation rates, BDs can be considerably oblate. We present a conics-based radiative transfer (RT) scheme for computing the disc-resolved and disc-integrated polarized emission of an oblate BD or extrasolar giant planet (EGP) bearing homogenous or patchy clouds. Assuming a uniform grey atmosphere, we theoretically examine the photopolarimetric sensitivity to its scattering properties like cloud optical thickness and grain-size, concurrently with BD properties, like oblateness, inclination and effective temperature, which are all treated as free parameters. Additionally, we examine the potential effects of gravitational darkening (GD), revealing that it could significantly amplify disc-integrated polarization. GD imparts a non-linear inverse temperature-dependence to the resulting polarization. Photopolarimetric observations are sensitive to oblateness and inclination. The degree-of-polarization (DoP) increases in response to both, making it potentially useful for assessing the spatial orientation of the BD. Under our model assumptions, increasing droplet size in optically thick clouds causes a blue-ward shift in near-infrared (NIR) colors of BDs --- interesting in view of the observed J-K brightening in L/T transition. For large cloud grains, polarization decreases sharply, while transmitted intensity shows a steady increase. BD polarization is thus a potential indicator not only of the presence of clouds but also provides information on cloud grain size.

We present results of simulations of the climate of the newly discovered planet Proxima Centauri B, performed using the Met Office Unified Model (UM). We examine the responses of both an `Earth-like' atmosphere and simplified nitrogen and trace carbon dioxide atmosphere to the radiation likely received by Proxima Centauri B. Additionally, we explore the effects of orbital eccentricity on the planetary conditions using a range of eccentricities guided by the observational constraints. Overall, our results are in agreement with previous studies in suggesting Proxima Centauri B may well have surface temperatures conducive to the presence of liquid water. Moreover, we have expanded the parameter regime over which the planet may support liquid water to higher values of eccentricity (>= 0.1) and lower incident fluxes (881.7 Wm-2) than previous work. This increased parameter space arises because of the low sensitivity of the planet to changes in stellar flux, a consequence of the stellar spectrum and orbital configuration. However, we also find interesting differences from previous simulations, such as cooler mean surface temperatures for the tidally-locked case. Finally, we have produced high resolution planetary emission and reflectance spectra, and highlight signatures of gases vital to the evolution of complex life on Earth (oxygen, ozone and carbon dioxide).