We visually analyzed the Transit Timing Variation (TTV) data of 5930 Kepler Objects of Interest (KOIs) homogeneously. Using data from Rowe et al. 2014 and Holczer et al. 2016, we investigated TTVs for KOIs in Kepler's Data Release 24 catalog. Using TTV plots, periodograms, and folded quadratic sinusoid fits, we visually rated each KOI's TTV data in five categories. Our ratings emphasize the hundreds of planets with TTVs that are weaker than the 200 or so that have been studied in detail. Our findings are consistent with statistical methods for identifying strong TTVs (Holczer et al. 2016), though we find some additional systems worth investigation. Between about 3-50 days and 1.3-6 Earth radii, the frequency of strong TTVs increases with period and radius. As expected, strong TTVs are very common when period ratios are near a resonance, but there is not a one-to-one correspondence. The observed planet-by-planet frequency of strong TTVs is only somewhat lower in systems with 1-2 known planets (7 +/- 1 percent) than in systems with 3+ known planets (11 +/- 2 percent). We attribute TTVs to known planets in multi-transiting systems, but find about 30 cases where the perturbing planet is unknown. Our conclusions are valuable as an ensemble for learning about planetary system architectures and individually as stepping stones towards more detailed mass-radius constraints. We also discuss Data Release 25 TTVs, about 100 KOIs with Transit Duration and/or Depth Variations, and that the Transiting Exoplanet Survey Satellite (TESS) will likely find roughly 10 planets with strong TTVs.

We present a new transit timing catalog of 2599 Kepler Objects of Interest (=KOIs), using the PDC-MAP long-cadence light curves that include the full seventeen quarters of the mission (ftp://wise- ftp.tau.ac.il/pub/tauttv/TTV/ver_112). The goal is to produce an easy-to-use catalog that can stimulate further analyses of interesting systems. For 779 KOIs with high enough SNRs, we derived the timing, duration and depth of 69,914 transits. For 1820 KOIs with lower SNR, we derived only the timing of 225,273 transits. After removal of outlier timings, we derived various statistics for each KOI that were used to indicate significant variations. Including systems found by previous works, we have detected 260 KOIs which showed significant TTVs with long-term variations (>100 day), and another fourteen KOIs with periodic modulations shorter than 100 day and small amplitudes. For five of those, the periodicity is probably due to the crossing of rotating stellar spots by the transiting planets.

A few studies have reported a significant dearth of exoplanets with Neptune mass and radius with orbital periods below $2$--$4$ d. This cannot be explained by observational biases because many Neptunian planets with longer orbital periods have been detected. The existence of this desert is similar to the appearance of the so-called brown-dwarf desert that suggests different formation mechanisms of planets and stellar companions with short orbital periods. Similarly, the Neptunian desert might indicate different mechanisms of formation and evolution for hot Jupiters and short-period super-Earths. We here follow a previous study and examine the location and shape of the desert in both the period-mass and period-radius planes, using the currently available large samples of planets. The desert in the period-mass plane has a relatively sharp upper edge, with a planetary mass that is inversely proportional to the planetary orbital period, while the lower, somewhat blurred, boundary is located along masses that are apparently linearly proportional to the period. The desert in the period-radius plane of the transiting planets is less clear. It seems as if the radius along the upper boundary is inversely proportional to the period to the power of one-third, while the lower boundary shows a radius that is proportional to the period to the power of two-thirds. The combination of the two upper bounds of the desert, in the period-mass and period-radius planes, yields a planetary mass-radius relation of $R_{\rm p}/R_{\rm Jup}\simeq (1.2\pm0.3)(M_{\rm p}/M_{\rm Jup})^{0.27\pm0.11}\,$ for $ 0.1\lesssim M_{\rm p}/M_{\rm Jup}\lesssim 1$. The derived shape of the desert, which might extend up to periods of $5$--$10$ d, could shed some light on the formation and evolution of close-in planets.

We infer dynamical masses in eight multi-planet systems using transit times measured from Kepler's complete dataset, including short-cadence data where available. Of the eighteen dynamical masses that we infer, ten pass multiple tests for robustness. These are in systems; Kepler-26 (KOI-250), Kepler-29 (KOI-738), Kepler-60 (KOI-2086), Kepler-105 (KOI-115), and Kepler-307 (KOI-1576). Kepler-105 c has a density consistent with an Earth-like composition. Strong TTV signals were detected from additional planets, but their inferred masses were sensitive to outliers or consistent solutions could not be found with independently-measured transit times, including planets at; Kepler-49 (KOI-248), Kepler-57 (KOI-1270), Kepler-105 (KOI-115) and Kepler-177 (KOI-523). Strong upper limits on the mass of Kepler-177 c imply an extremely low density ~0.1 g cm$^{-3}$. In most cases, individual orbital eccentricities were poorly constrained due to degeneracies in TTV inversion. For five planet pairs in our sample, strong secular interactions imply a moderate-to-high likelihood of apsidal alignment over a wide range of possible eccentricities. We also find solutions for the three planets known to orbit Kepler-60 in a Laplace-like resonance chain. However, non-librating solutions also match the transit-timing data. For six systems, we calculate more precise stellar parameters than previously known, enabling useful constraints on planetary densities where we have robust mass measurements. Placing these exoplanets on the mass-radius diagram, we find a wide range of densities is observed among sub-Neptune mass planets and that the range in observed densities is anti-correlated with incident flux.

We model the observed X-ray spectral continuum shape, ionic column densities, and absorption measure distribution (AMD) of the warm absorber in the Seyfert galaxy NGC 3783. We assume a photo-ionized medium with a uniform total (gas+radiation) pressure. The irradiation causes the wind to be radiation pressure compressed (RPC). We compare the observational characteristics derived from the 900 ksec Chandra observation to radiative transfer computations in pressure equilibrium using the radiative transfer code TITAN. We explore different values of the ionization parameter xi of the incident flux and adjust the hydrogen-equivalent column density, N_H0 of the warm absorber to match the observed soft X-ray continuum. We derive theoretical column densities for a broad range of ionic species of iron and neon and an AMD that we compare to the observations. We find an extension of the degeneracy between xi and N_H0 for the constant pressure models previously discussed for NGC 3783. Including the ionic column densities of iron and neon in the comparison between observations and data we conclude that a range of ionization parameters between 4000 and 8000 ergs cm/s is preferred. For the first time, we present theoretical AMD for a constant pressure wind in NGC 3783 that correctly reproduces the observed level and is in approximate agreement with the observational appearance of an instability region. Using a variety of observational indicators, we confirm that the X-ray outflow of NGC 3783 can be described as an RPC medium in pressure equilibrium. The observed AMD agrees with a uniformly hot or a uniformly cold thermal state. The measured ionic column densities suggest that the wind tends to the uniformly cold thermal state. The occurrence of thermal instability in the warm absorber model may depend on the computational method and the spatial scale on which the radiative transfer is solved.

Mazeh, Holczer, and Shporer (2015) have presented an approach that can, in principle, use the derived transit timing variation (TTV) of some transiting planets observed by the $Kepler$ mission to distinguish between prograde and retrograde motion of their orbits with respect to their parent stars' rotation. The approach utilizes TTVs induced by spot-crossing events that occur when the planet moves across a spot on the stellar surface, looking for a correlation between the derived TTVs and the stellar brightness derivatives at the corresponding transits. This can work even in data that cannot temporally resolve the spot-crossing events themselves. Here we apply this approach to the $Kepler$ KOIs, identifying nine systems where the photometric spot modulation is large enough and the transit timing accurate enough to allow detection of a TTV-brightness-derivatives correlation. Of those systems five show highly significant prograde motion (Kepler-17b, Kepler-71b, KOI-883.01, KOI-895.01, and KOI-1074.01), while no system displays retrograde motion, consistent with the suggestion that planets orbiting cool stars have prograde motion. All five systems have impact parameter $0.2\lesssim b\lesssim0.5$, and all systems within that impact parameter range show significant correlation, except HAT-P-11b where the lack of a correlation follows its large stellar obliquity. Our search suffers from an observational bias against detection of high impact parameter cases, and the detected sample is extremely small. Nevertheless, our findings may suggest that stellar spots, or at least the larger ones, tend to be located at a low stellar latitude, but not along the stellar equator, similar to the Sun.

We report the spectroscopic confirmation of the Kepler object of interest KOI-183.01 (Kepler-423b), a half-Jupiter mass planet transiting an old solar-like star every 2.7 days. Our analysis is the first to combine the full Kepler photometry (quarters 1-17) with high-precision radial velocity measurements taken with the FIES spectrograph at the Nordic Optical Telescope. We simultaneously modelled the photometric and spectroscopic data-sets using Bayesian approach coupled with Markov chain Monte Carlo sampling. We found that the Kepler pre-search data conditioned (PDC) light curve of KOI-183 exhibits quarter-to-quarter systematic variations of the transit depth, with a peak-to-peak amplitude of about 4.3 % and seasonal trends reoccurring every four quarters. We attributed these systematics to an incorrect assessment of the quarterly variation of the crowding metric. The host star KOI-183 is a G4 dwarf with $M_\star=0.85\pm0.04$ M$_\rm{Sun}$, $R_\star=0.95\pm0.04$ R$_\rm{Sun}$, $T_\mathrm{eff}=5560\pm80$ K, $[M/H]=-0.10\pm0.05$ dex, and with an age of $11\pm2$ Gyr. The planet KOI-183b has a mass of $M_\mathrm{p}=0.595\pm0.081$ M$_\mathrm{Jup}$ and a radius of $R_\mathrm{p}=1.192\pm0.052$ R$_\mathrm{Jup}$, yielding a planetary bulk density of $\rho_\mathrm{p}=0.459\pm0.083$ g/cm$^{3}$. The radius of KOI-183b is consistent with both theoretical models for irradiated coreless giant planets and expectations based on empirical laws. The inclination of the stellar spin axis suggests that the system is aligned along the line of sight. We detected a tentative secondary eclipse of the planet at a 2-$\sigma$ confidence level ($\Delta F_{\mathrm{ec}}=14.2\pm6.6$ ppm) and found that the orbit might have a small non-zero eccentricity of $e=0.019^{+0.028}_{-0.014}$. With a Bond albedo of $A_\mathrm{B}=0.037\pm0.019$, KOI-183b is one of the gas-giant planets with the lowest albedo known so far.

Some transiting planets discovered by the Kepler mission display transit timing variations (TTVs) induced by stellar spots that rotate on the visible hemisphere of their parent stars. An induced TTV can be observed when a planet crosses a spot and modifies the shape of the transit light curve, even if the time resolution of the data does not allow to detect the crossing event itself. We present an approach that can, in some cases, use the derived TTVs of a planet to distinguish between a prograde and a retrograde planetary motion with respect to the stellar rotation. Assuming a single spot darker than the stellar disc, spot crossing by the planet can induce measured positive (negative) TTV, if the crossing occurs in the first (second) half of the transit. On the other hand, the motion of the spot towards (away from) the center of the stellar visible disc causes the stellar brightness to decrease (increase). Therefore, for a planet with prograde motion, the induced TTV is positive when the local slope of the stellar flux at the time of transit is negative, and vice versa. Thus, we can expect to observe a negative (positive) correlation between the TTVs and the photometric slopes for prograde (retrograde) motion. Using a simplistic analytical approximation, and also the publicly available SOAP-T tool to produce light curves of transits with spot-crossing events, we show for some cases how the induced TTVs depend on the local stellar photometric slopes at the transit timings. Detecting this correlation in Kepler transiting systems with high enough signal-to-noise ratio can allow us to distinguish between prograde and retrograde planetary motions. In coming papers we present analyses of the KOIs and Kepler eclipsing binaries, following the formalism developed here.

The winds of ionized gas driven by Active Galactic Nuclei (AGN) can be studied through absorption lines in their X-ray spectra. A recurring feature of these outflows is their broad ionization distribution, including essentially all ionization levels (e.g., Fe^0+ to Fe^25+). This characteristic feature can be quantified with the absorption measure distribution (AMD), defined as the distribution of column density with ionization parameter |dN / dlog xi|. Observed AMDs extend over 0.1 < xi < 10^4 (cgs), and are remarkably similar in different objects. Power-law fits (|dN /dlog xi| ~ N_1 xi^a) yield N_1 = 3x10^{21} cm^-2 +- 0.4 dex and a = 0 -- 0.4. What is the source of this broad ionization distribution, and what sets the small range of observed $N_1$ and $a$? A common interpretation is a multiphase outflow, with a wide range of gas densities in a uniform gas pressure medium. However, the incident radiation pressure leads to a gas pressure gradient in the photoionized gas, and therefore to a broad range of ionization states within a single slab. We show that this compression of the gas by the radiation pressure leads to an AMD with |dN / dlog xi| = 8x10^21 xi^0.03 cm^-2, remarkably similar to that observed. The calculated values of $N_1$ and $a$ depend weakly on the gas metallicity, the ionizing spectral slope, the distance from the nucleus, the ambient density, and the total absorber column. Thus, radiation pressure compression (RPC) of the photoionized gas provides a natural explanation for the observed AMD. RPC predicts that the gas pressure increases with decreasing ionization, which can be used to test the validity of RPC in ionized AGN outflows.

Following Ford et al. (2011, 2012) and Steffen et al. (2012) we derived the transit timing of 1960 Kepler KOIs using the pre-search data conditioning (PDC) light curves of the first twelve quarters of the Kepler data. For 721 KOIs with large enough SNRs, we obtained also the duration and depth of each transit. The results are presented as a catalog for the community to use. We derived a few statistics of our results that could be used to indicate significant variations. Including systems found by previous works, we have found 130 KOIs that showed highly significant TTVs, and 13 that had short-period TTV modulations with small amplitudes. We consider two effects that could cause apparent periodic TTV - the finite sampling of the observations and the interference with the stellar activity, stellar spots in particular. We briefly discuss some statistical aspects of our detected TTVs. We show that the TTV period is correlated with the orbital period of the planet and with the TTV amplitude.

By analyzing the X-ray spectra of NGC 3516 from 2001 and 2006 obtained with the HETGS spectrometer on board the Chandra observatory, we find that the kinematic structure of the outflow can be well represented by four outflow components intrinsic to NGC 3516. The outflow velocities of the different components are 350 +-100 km s-1, 1500 +-150 km s-1, 2600 +-200 km s-1 and 4000 +-400 km s-1 for components 1, 2, 3 and 4, respectively. A local component at z = 0 could be confused with intrinsic component 3. Components 1 and 2 have a broad range of ionization manifested by absorption from 23 different charge states of Fe. Component 3 and 4 are more highly ionized and show absorption from only 9 different charge states of Fe, however we were able to reconstruct the absorption measure distribution (AMD) for all four. The total column density of each component is NH = (1.8+- 0.5) X10^22 cm-2, NH = (2.5+- 0.3) X10^22 cm-2, NH = (6.9+- 4.3) X10^22 cm-2 and NH = (5.4+- 1.2) X10^22 cm-2, respectively. The fast components 3 and 4 appear only in the high state of 2006 and not in 2001, while the slower components persist during both epochs. On the other hand, there is no significant absorption variability within days during 2001 or during 2006. We find that covering factor plays a minor role for the line absorption.

Context. The study of abundances in the nucleus of active galaxies allows us to investigate the evolution of abundance by comparing local and higher redshift galaxies. However, the methods used so far have substantial drawbacks or rather large uncertainties. Some of the measurements are at odds with the initial mass function derived from the older stellar population of local elliptical galaxies. Aims. We determine accurate and reliable abundances of C, N, Ne, and Fe relative to O from the narrow absorption lines observed in the X-ray spectra of Mrk 509. Methods. We use the stacked 600 ks XMM-Newton RGS and 180 ks Chandra LETGS spectra. Thanks to simultaneous observations with INTEGRAL and the optical monitor on-board XMM-Newton for the RGS observations and HST-COS and Swift for the LETGS observations, we have an individual spectral energy distribution for each dataset. Owing to the excellent quality of the RGS spectrum, the ionisation structure of the absorbing gas is well constrained, allowing for a reliable abundance determination using ions over the whole observed range of ionisation parameters. Results. We find that the relative abundances are consistent with the proto-solar abundance ratios: C/O = 1.19$\pm$0.08, N/O = 0.98$\pm$0.08, Ne/O = 1.11$\pm$0.10, Mg/O = 0.68$\pm$0.16, Si/O = 1.3$\pm$0.6, Ca/O = 0.89$\pm$0.25, and Fe/O = 0.85$\pm$0.06, with the exception of S, which is slightly under-abundant, S/O = 0.57$\pm$0.14. Our results, and their implications, are discussed and compared to the results obtained using other techniques to derive abundances in galaxies.

We present the results of our 600 ks RGS observation as part of the multiwavelength campaign on Mrk 509. The very high quality of the spectrum allows us to investigate the ionized outflow with an unprecedented accuracy due to the long exposure and the use of the RGS multipointing mode. We detect multiple absorption lines from the interstellar medium and from the ionized absorber in Mrk 509. A number of emission components are also detected, including broad emission lines consistent with an origin in the broad line region, the narrow OVII forbidden emission line and also (narrow) radiative recombination continua. The ionized absorber consists of two velocity components (v = -13 \pm 11 km/s and v = -319 \pm 14 km/s), which both are consistent with earlier results, including UV data. There is another tentative component outflowing at high velocity, -770 \pm 109 km/s, which is only seen in a few highly ionized absorption lines. The outflow shows discrete ionization components, spanning four orders of magnitude in ionization parameter. Due to the excellent statistics of our spectrum, we demonstrate for the first time that the outflow in Mrk 509 in the important range of log xi between 1-3 cannot be described by a smooth, continuous absorption measure distribution, but instead shows two strong, discrete peaks. At the highest and lowest ionization parameters we cannot differentiate smooth and discrete components.

By analyzing the X-ray spectrum of MCG-6-30-15 obtained with the HETG spectrometer on board the Chandra Observatory, we identify three kinematically distinct absorption systems; two outflow components intrinsic to MCG-6-30-15, and one local at z = 0. The slow outflow at -100 +/- 50 km s^-1 has a large range of ionization manifested by absorption from 24 different charge states of Fe, which enables a detailed reconstruction of the absorption measure distribution (AMD). This AMD spans five orders of magnitude in ionization parameter: -1.5 < log xi < 3.5 (cgs units), with a total column density of N_H = (5.3 +/- 0.7) x 10^21 cm^-2. The fast outflow at -1900 +/- 150 km s^-1 has a well defined ionization parameter with log xi = 3.82 +/- 0.03 (cgs units) and column density N_H = 8.1 +/- 0.7 x 10^22 cm^-2. Assuming this component is a thin, uniform, spherical shell, it can be estimated to lie within 11 light days of the \agn center. The third component, most clearly detected in the lower oxygen charge states O^+1 - O^+6, has been confused in the past with the fast outflow, but is identified here with local gas z = 0 and a total column density N_H of a few 10^20 cm^-2. Finally, we exploit the excellent spectral resolution of the HETGS and use the present spectrum to determine the rest-frame wavelengths of oxygen inner-shell lines that were previously uncertain.

We analyze the Chandra X-ray spectrum obtained with the HETGS grating spectrometer of IRAS 13349+2438, which has one of the richest absorption spectra of a quasar outflow. Absorption from almost all charge states of Fe is detected. This allows for a detailed reconstruction of the absorption measure distribution (AMD), which we define as the continuous distribution of column density as a function of ionization parameter. We find a double peaked AMD for IRAS 13349+2438 with a total (ionized) column density of N_H = (1.2 +- 0.3)* 10^{22} cm^-2} assuming solar iron abundance. For comparison, we perform a similar analysis on the well studied HETGS spectrum of NGC 3783. Both sources feature a deep minimum in column density consistent with no absorption from gas at temperatures of 4.5 < log T < 5 (K). We interpret the minima as observational evidence for thermal instability in this temperature regime.

Although AGN feedback through ionised winds is of great importance in models of AGN/galaxy coevolution, the mass and energy output via these winds, even in the nearby universe, is poorly understood. The issue is complicated by the wide range of ionisation in the winds, which means that multiwavelength observational campaigns are required to obtain the complete picture. In this paper, we use a ~ 160 ks XMM-Newton RGS spectrum to get the most accurate view yet of the ionised outflow (warm absorber) in NGC 7469 as seen in X-rays, finding that there is a wide range of ionisation, with log xi in the range ~ 0.5-3.5 erg cm s^-1, and two main velocity regimes, at 580-720 and 2300 km s^-1, with the highest velocity gas being the least ionised. The total absorbing column density in the X-rays is of order 3 x 10^21 cm^-2. We find that the lowest ionisation phase of the absorber is probably identical with one of the phases of the UV absorber discovered in previous studies. We show that both X-ray and UV absorbers are consistent with an origin near the base of a torus wind, where matter is being launched and accelerated. Calculating the mass outflow rate and kinetic luminosity of all the absorber phases, we demonstrate that the X-ray absorbing gas carries respectively ~ 90% and 95% of the mass and kinetic energy output of the ionised outflow.

We provide improved atomic calculation of wavelengths, oscillator strengths, and autoionization rates relevant to the $2\to 3$ inner-shell transitions of Fe~VI--XVI, the so-called Fe~M-shell unresolved transition array (UTA). A second order many-body perturbation theory is employed to obtain accurate transition wavelengths, which are systematically larger than previous theoretical results by 15--45~m{\AA}. For a few transitions of Fe~XVI and Fe~XV where laboratory measurements exist, our new wavelengths are accurate to within a few m{\AA}. Using these new calculations, the apparent discrepancy in the velocities between the Fe~M-shell UTA and other highly ionized absorption lines in the outflow of NGC 3783 disappears. The oscillator strengths in our new calculation agree well with the previous theoretical data, while the new autoionization rates are significantly larger, especially for lower charge states. We attribute this discrepancy to the missing autoionization channels in the previous calculation. The increased autoionization rates may slightly affect the column density analysis of the Fe~M-shell UTA for sources with high column density and very low turbulent broadening. The complete set of atomic data is provided as an electronic table.

The X-ray emission of many active galactic nuclei (AGNs) is absorbed between 15 and 17 Angstrom by means of unresolved (inner-shell) transition arrays (UTAs) of Fe M-shell ions. The outflow velocities implied by the Doppler shifts of the individual UTAs in the spectrum have never before been measured. Thus, the Fe-M absorber has been commonly assumed to be part of the ionized AGN outflow, whose velocities are readily obtained from more easily measured spectral lines. The best spectrum of Fe-M absorption is available from the integrated 900 ks Chandra HETGS observations of NGC 3783, in which some Fe-M ions are clearly resolved. We measure the velocities of the individual Fe-M ions in NGC 3783 for the first time. Surprisingly, we find that the Fe-M absorber, most noticeably Fe$^{+8}$, Fe$^{+9}$, and Fe$^{+10}$, is not outflowing at the same velocity as the previously known wind. In fact, it appears to be stationary and therefore not part of the outflow at all. It could, alternatively, be ascribed to the skin of the dusty torus. This reduces appreciably the mass loss rate estimated for the NGC 3783 outflow and perhaps for other similar sources as well, in which the various Fe-M ions are not resolved.