sort results by

Use logical operators AND, OR, NOT and round brackets to construct complex queries. Whitespace-separated words are treated as ANDed.

Show articles per page in mode

Malik, Mohit

Normalized to: Malik, M.

26 article(s) in total. 144 co-authors, from 1 to 11 common article(s). Median position in authors list is 4,0.

[1]  oai:arXiv.org:1910.01567  [pdf] - 2027894
Evidence for H$_{2}$ Dissociation and Recombination Heat Transport in the Atmosphere of KELT-9b
Comments:
Submitted: 2019-10-03, last modified: 2020-01-07
Phase curve observations provide an opportunity to study the full energy budgets of exoplanets by quantifying the amount of heat redistributed from their daysides to their nightsides. Theories explaining the properties of phase curves for hot Jupiters have focused on the balance between radiation and dynamics as the primary parameter controlling heat redistribution. However, recent phase curves have shown deviations from the trends that emerge from this theory, which has led to work on additional processes that may affect hot Jupiter energy budgets. One such process, molecular hydrogen dissociation and recombination, can enhance energy redistribution on ultra-hot Jupiters with temperatures above $\sim2000$ K. In order to study the impact of H$_{2}$ dissociation on ultra-hot Jupiters, we present a phase curve of KELT-9b observed with the Spitzer Space Telescope at 4.5 $\mu$m. KELT-9b is the hottest known transiting planet, with a 4.5-$\mu$m dayside brightness temperature of $4566^{+140}_{-136}$ K and a nightside temperature of $2556^{+101}_{-97}$ K. We observe a phase curve amplitude of $0.609 \pm 0.020$ and a hot spot offset of $18.7^{+2.1}_{-2.3}$ degrees. The observed amplitude is too small to be explained by a simple balance between radiation and advection. General circulation models (GCMs) and an energy balance model that include the effects of H$_{2}$ dissociation and recombination provide a better match to the data. The GCMs, however, predict a maximum hot spot offset of $5$ degrees, which disagrees with our observations at $>5\sigma$ confidence. This discrepancy may be due to magnetic effects in the planet's highly ionized atmosphere.
[2]  oai:arXiv.org:1910.11795  [pdf] - 2026138
Supervised Machine Learning for Intercomparison of Model Grids of Brown Dwarfs: Application to GJ 570D and the Epsilon Indi B Binary System
Comments: Accepted for publication in The Astronomical Journal
Submitted: 2019-10-25, last modified: 2019-12-18
Self-consistent model grids of brown dwarfs involve complex physics and chemistry, and are often computed using proprietary computer codes, making it challenging to identify the reasons for discrepancies between model and data as well as between the models produced by different research groups. In the current study, we demonstrate a novel method for analyzing brown dwarf spectra, which combines the use of the Sonora, AMES-Cond and HELIOS model grids with the supervised machine learning method of the random forest. Besides performing atmospheric retrieval, the random forest enables information content analysis of the three model grids as a natural outcome of the method, both individually on each grid and by comparing the grids against one another, via computing large suites of mock retrievals. Our analysis reveals that the different choices made in modelling the alkali line shapes hinder the use of the alkali lines as gravity indicators. Nevertheless, the spectrum longward of 1.2 micron encodes enough information on the surface gravity to allow its inference from retrieval. Temperature may be accurately and precisely inferred independent of the choice of model grid, but not the surface gravity. We apply random forest retrieval to three objects: the benchmark T7.5 brown dwarf GJ 570D; and Epsilon Indi Ba (T1.5 brown dwarf) and Bb (T6 brown dwarf), which are part of a binary system and have measured dynamical masses. For GJ 570D, the inferred effective temperature and surface gravity are consistent with previous studies. For Epsilon Indi Ba and Bb, the inferred surface gravities are broadly consistent with the values informed by the dynamical masses.
[3]  oai:arXiv.org:1907.13135  [pdf] - 2010011
Analyzing Atmospheric Temperature Profiles and Spectra of M dwarf Rocky Planets
Comments: published in ApJ. Also see these three companion papers: 1. Mansfield et al. 2019 (ApJ, 886, 2), "Identifying Atmospheres on Rocky Exoplanets Through Inferred High Albedo", 2. Koll et al. 2019 (ApJ, 886, 2), "Identifying candidate atmospheres on rocky M dwarf planets via eclipse photometry", 3. Koll (submitted) "A Scaling Theory for Atmospheric Heat Redistribution on Rocky Exoplanets"
Submitted: 2019-07-30, last modified: 2019-12-05
The James Webb Space Telescope (JWST) will open up the possibility of comprehensively measuring the emission spectra of rocky exoplanets orbiting M dwarfs to detect and characterize their atmospheres. In preparation for this opportunity, we present model atmospheres for three M-dwarf planets particularly amenable to secondary eclipse spectroscopy -- TRAPPIST-1b, GJ 1132b, and LHS 3844b. Using three limiting cases of candidate atmospheric compositions (pure H2O, pure CO2 and solar abundances) we calculate temperature-pressure profiles and emission and reflection spectra in radiative-convective equilibrium, including the effects of a solid surface. We find that the atmospheric radiative transfer is significantly influenced by the cool M-star irradiation; H2O and CO2 absorption bands in the near-infrared are strong enough to absorb a sizeable fraction of the incoming stellar light at low pressures, which leads to temperature inversions in the upper atmosphere. The non-gray band structure of gaseous opacities in the infrared is hereby an important factor. Opacity windows are muted at higher atmospheric temperatures, so we expect temperature inversions to be common only for sufficiently cool planets. We also find that pure CO2 atmospheres exhibit lower overall temperatures and stronger reflection spectra compared to models of the other two compositions. We estimate that for GJ 1132b and LHS 3844b we should be able to distinguish between different atmospheric compositions with JWST. The emission lines from the predicted temperature inversions are currently hard to measure, but high resolution spectroscopy with future Extremely Large Telescopes may be able to detect them.
[4]  oai:arXiv.org:1907.13150  [pdf] - 2025658
Identifying Atmospheres on Rocky Exoplanets Through Inferred High Albedo
Comments: Accepted to ApJ. Also see these three companion papers: 1. Koll et al (accepted to ApJ), "Identifying Candidate Atmospheres on Rocky M Dwarf Planets Via Eclipse Photometry", 2. Malik et al (accepted to ApJ), "Analyzing Atmospheric Temperature Profiles and Spectra of M dwarf Rocky Planets", and 3. Koll (accepted to ApJ) "A Scaling Theory for Atmospheric Heat Redistribution on Rocky Exoplanets"
Submitted: 2019-07-30, last modified: 2019-11-22
The upcoming launch of the James Webb Space Telescope (JWST) means that we will soon have the capability to characterize the atmospheres of rocky exoplanets. However, it is still unknown whether such planets orbiting close to M dwarf stars can retain their atmospheres, or whether high-energy irradiation from the star will strip the gaseous envelopes from these objects. We present a new method to detect an atmosphere on a synchronously rotating rocky exoplanet around a K/M dwarf, by using thermal emission during secondary eclipse to infer a high dayside albedo that could only be explained by bright clouds. Based on calculations for plausible surface conditions, we conclude that a high albedo could be unambiguously interpreted as a signal of an atmosphere for planets with substellar temperatures of $T_{sub}=$ 410-1250 K. This range corresponds to equilibrium temperatures of $T_{eq}=$ 300-880 K. We compare the inferred albedos of eight possible planet surface compositions to cloud albedo calculations. We determine that a layer of clouds with optical depths greater than $\tau=0.5$-$7$, would have high enough albedos to be distinguishable from a bare rock surface. This method of detecting an atmosphere on a rocky planet is complementary to existing methods for detecting atmospheres, because it provides a way to detect atmospheres with pressures below 1 bar (e.g. Mars), which are too tenuous to transport significant heat but thick enough to host high-albedo clouds.
[5]  oai:arXiv.org:1907.13138  [pdf] - 2025657
Identifying Candidate Atmospheres on Rocky M dwarf Planets via Eclipse Photometry
Comments: Submitted to ApJ. Also see these three companion papers: 1. Mansfield et al (submitted), "Identifying Atmospheres on Rocky Exoplanets Through Inferred High Albedo", 2. Malik et al (submitted), "Analyzing Atmospheric Temperature Profiles and Spectra of M dwarf Rocky Planets", and 3. Koll (submitted) "A Scaling Theory for Atmospheric Heat Redistribution on Rocky Exoplanets"
Submitted: 2019-07-30
Most rocky planets in the galaxy orbit a cool host star, and there is large uncertainty among theoretical models whether these planets can retain an atmosphere. The James Webb Space Telescope (JWST) might be able to settle this question empirically, but most proposals for doing so require large observational effort because they are based on spectroscopy. Here we show that infrared photometry of secondary eclipses could quickly identify "candidate" atmospheres, by searching for rocky planets with atmospheres thick enough that atmospheric heat transport noticeably reduces their dayside thermal emission compared to that of a bare rock. For a planet amenable to atmospheric follow-up, we find that JWST should be able to confidently detect the heat redistribution signal of an O(1) bar atmosphere with one to two eclipses. One to two eclipses is generally much less than the effort needed to infer an atmosphere via transmission or emission spectroscopy. Candidate atmospheres can be further validated via follow-up spectroscopy or phase curves. In addition, because this technique is fast it could enable a first atmospheric survey of rocky exoplanets with JWST. We estimate that the TESS mission will find ~100 planets that are too hot to be habitable but that can be quickly probed via eclipse photometry. Knowing whether hot, rocky planets around M dwarfs have atmospheres is important not only for understanding the evolution of uninhabitable worlds: if atmospheres are common on hot planets, then cooler, potentially habitable planets around M dwarfs are also likely to have atmospheres.
[6]  oai:arXiv.org:1812.02027  [pdf] - 1859942
Evolutionary models of cold and low-mass planets: Cooling curves, magnitudes, and detectability
Comments: Language corrected version and improved arrangements of figures, online data at: http://www.space.unibe.ch/research/research_groups/planets_in_time/numerical_data/index_eng.html
Submitted: 2018-12-05, last modified: 2019-03-27
Future instruments like NIRCam and MIRI on JWST or METIS at the ELT will be able to image exoplanets that are too faint for current direct imaging instruments. Evolutionary models predicting the planetary intrinsic luminosity as a function of time have traditionally concentrated on gas-dominated giant planets. We extend these cooling curves to Saturnian and Neptunian planets. We simulate the cooling of isolated core-dominated and gas giant planets with masses of 5 Earthmasses to 2 Jupitermasses. The luminosity includes the contribution from the cooling and contraction of the core and of the H/He envelope, as well as radiogenic decay. For the atmosphere we use grey, AMES-Cond, petitCODE, and HELIOS models. We consider solar and non-solar metallicities as well as cloud-free and cloudy atmospheres. The most important initial conditions, namely the core-to-envelope ratio and the initial luminosity are taken from planet formation simulations based on the core accretion paradigm. We first compare our cooling curves for Uranus, Neptune, Jupiter, Saturn, GJ 436b, and a 5 Earthmass-planet with a 1% H/He envelope with other evolutionary models. We then present the temporal evolution of planets with masses between 5 Earthmasses and 2 Jupitermasses in terms of their luminosity, effective temperature, radius, and entropy. We discuss the impact of different post formation entropies. For the different atmosphere types and initial conditions magnitudes in various filter bands between 0.9 and 30 micrometer wavelength are provided. Using black body fluxes and non-grey spectra, we estimate the detectability of such planets with JWST. It is found that a 20 (100) Earthmass-planet can be detected with JWST in the background limit up to an age of about 10 (100) Myr with NIRCam and MIRI, respectively.
[7]  oai:arXiv.org:1903.06794  [pdf] - 1868170
Self-luminous and irradiated exoplanetary atmospheres explored with HELIOS
Comments: 33 pages, 18 figures, accepted for publication in AJ
Submitted: 2019-03-15
We present new methodological features and physical ingredients included in the 1D radiative transfer code HELIOS, improving the hemispheric two-stream formalism. We conduct a thorough intercomparison survey with several established forward models, including COOLTLUSTY, PHOENIX, and find satisfactory consistency with their results. Then, we explore the impact of (i) different groups of opacity sources, (ii) a stellar path length adjustment, and (iii) a scattering correction on self-consistently calculated atmospheric temperatures and planetary emission spectra. First, we observe that temperature-pressure (T-P) profiles are very sensitive to the opacities included, with metal oxides, hydrides, the alkali atoms (and ionized hydrogen) playing an important role for the absorption of shortwave radiation (in very hot surroundings). Moreover, if these species are sufficiently abundant, they are likely to induce non-monotonic T-P profiles. Second, without the stellar path length adjustment, the incoming stellar flux is significantly underestimated for zenith angles above 80{\deg}, which somewhat affects the upper atmospheric temperatures and the planetary emission. Third, the scattering correction improves the accuracy of the computation of the reflected stellar light by ~10%. We use HELIOS to calculate a grid of cloud-free atmospheres in radiative-convective equilibrium for self-luminous planets for a range of effective temperatures, surface gravities, metallicities, and C/O ratios, to be used by planetary evolution studies. Furthermore, we calculate dayside temperatures and secondary eclipse spectra for a sample of exoplanets for varying chemistry and heat redistribution. These results may be used to make predictions on the feasibility of atmospheric characterizations with future observations.
[8]  oai:arXiv.org:1808.00501  [pdf] - 1806051
Three-Dimensional Circulation Driving Chemical Disequilibrium in WASP-43b
Comments: Accepted for publication in ApJ
Submitted: 2018-08-01, last modified: 2018-10-30
Spectral features in the observed spectra of exoplanets depend on the composition of their atmospheres. A good knowledge of the main atmospheric processes that drive the chemical distribution is therefore essential to interpret exoplanetary spectra. An atmosphere reaches chemical equilibrium if the rates of the forward and backward chemical reactions converge to the same value. However, there are atmospheric processes, such as atmospheric transport, that destabilize this equilibrium. In this work we study the changes in composition driven by a 3D wind field in WASP-43b using our Global Circulation Model, THOR. Our model uses validated temperature- and pressure-dependent chemical timescales that allow us to explore the disequilibrium chemistry of CO, CO$_2$, H$_2$O and CH$_4$. In WASP-43b the formation of the equatorial jet has an important impact in the chemical distribution of the different species across the atmosphere. At low latitudes the chemistry is longitudinally quenched, except for CO$_2$ at solar abundances. The polar vortexes have a distinct chemical distribution since these are regions with lower temperature and atmospheric mixing. Vertical and latitudinal mixing have a secondary impact in the chemical transport. We determine graphically the effect of disequilibrium on observed emission spectra. Our results do not show any significant differences in the emission spectra between the equilibrium and disequilibrium solutions for C/O = 0.5. However, if C/O is increased to 2.0, differences in the spectra due to the disequilibrium chemistry of CH$_4$ become non-negligible. In some spectral ranges the emission spectra can have more than 15$\%$ departures from the equilibrium solution.
[9]  oai:arXiv.org:1803.07087  [pdf] - 1761472
Confronting phantom inflation with Planck data
Comments: 9 pages, 4 figures, 2 tables, Accepted in Astrophysics and Space Science
Submitted: 2018-03-19, last modified: 2018-09-25
The latest Planck results are in excellent agreement with the theoretical expectations predicted from standard normal inflation based on slow-roll approximation which assumes equation-of-state $\omega \geq-1$. In this work, we study the phantom inflation ($\omega<-1$) as an alternative cosmological model within the slow-climb approximation using two hybrid inflationary fields. We perform Chain Monte Carlo analysis to determine the posterior distribution and best fit values for the cosmological parameters using Planck data and show that current CMB data does not discriminate between normal and phantom inflation. Interestingly, unlike in normal inflation, $\omega$ in phantom induced inflation evolves very slowly away from $-1$ during the inflation. Furthermore, in contrast to the standard normal inflation for which only upper bound on tensor-to-scalar ratio $r$ are possible, we obtain both upper and lower bounds for the two hybrid fields in the phantom scenario. Finally, we discuss prospects of future high precision polarization measurements and show that it may be possible to establish the dominance of one model over the other.
[10]  oai:arXiv.org:1807.10769  [pdf] - 1775614
Diagnosing aerosols in extrasolar giant planets with cross-correlation function of water bands
Comments: Accepted for publication in A&A; 19 pages, 15 figures, 5 tables
Submitted: 2018-07-27
Transmission spectroscopy with ground-based, high-resolution instruments provides key insight into the composition of exoplanetary atmospheres. Molecules such as water and carbon monoxide have been unambiguously identified in hot gas giants through cross- correlation techniques. A maximum in the cross-correlation function (CCF) is found when the molecular absorption lines in a binary mask or model template match those contained in the planet. Here, we demonstrate how the CCF method can be used to diagnose broad-band spectroscopic features such as scattering by aerosols in high-resolution transit spectra. The idea consists in exploiting the presence of multiple water bands from the optical to the near-infrared. We produced a set of models of a typical hot Jupiter spanning various conditions of temperature and aerosol coverage. We demonstrate that comparing the CCFs of individual water bands for the models constrains the presence and the properties of the aerosol layers. The contrast difference between the CCFs of two bands can reach about 100 ppm, which could be readily detectable with current or upcoming high-resolution stabilized spectrographs spanning a wide spectral range, such as ESPRESSO, CARMENES, HARPS-N+GIANO, HARPS+NIRPS, SPIRou or CRIRES+.
[11]  oai:arXiv.org:1804.07137  [pdf] - 1739927
The Peculiar Atmospheric Chemistry of KELT-9b
Comments: Accepted by ApJ. 9 pages, 6 figures. Corrected minor errors in Figures 1a and 1b (some line styles were switched by accident), text and conclusions unchanged, these minor changes will be updated in final ApJ proof
Submitted: 2018-04-19, last modified: 2018-07-12
The atmospheric temperatures of the ultra-hot Jupiter KELT-9b straddle the transition between gas giants and stars, and therefore between two traditionally distinct regimes of atmospheric chemistry. Previous theoretical studies assume the atmosphere of KELT-9b to be in chemical equilibrium. Despite the high ultraviolet flux from KELT-9, we show using photochemical kinetics calculations that the observable atmosphere of KELT-9b is predicted to be close to chemical equilibrium, which greatly simplifies any theoretical interpretation of its spectra. It also makes the atmosphere of KELT-9b, which is expected to be cloudfree, a tightly constrained chemical system that lends itself to a clean set of theoretical predictions. Due to the lower pressures probed in transmission (compared to emission) spectroscopy, we predict the abundance of water to vary by several orders of magnitude across the atmospheric limb depending on temperature, which makes water a sensitive thermometer. Carbon monoxide is predicted to be the dominant molecule under a wide range of scenarios, rendering it a robust diagnostic of the metallicity when analyzed in tandem with water. All of the other usual suspects (acetylene, ammonia, carbon dioxide, hydrogen cyanide, methane) are predicted to be subdominant at solar metallicity, while atomic oxygen, iron and magnesium are predicted to have relative abundances as high as 1 part in 10,000. Neutral atomic iron is predicted to be seen through a forest of optical and near-infrared lines, which makes KELT-9b suitable for high-resolution ground-based spectroscopy with HARPS-N or CARMENES. We summarize future observational prospects of characterizing the atmosphere of KELT-9b.
[12]  oai:arXiv.org:1804.04961  [pdf] - 1728934
Analytical Models of Exoplanetary Atmospheres. VI. Full Solutions for Improved Two-stream Radiative Transfer Including Direct Stellar Beam
Comments: Accepted by ApJS. 8 pages, 1 figure, 1 table, 40 equations
Submitted: 2018-04-13, last modified: 2018-07-03
Two-stream radiative transfer is a workhorse in the Earth, planetary and exoplanetary sciences due to its simplicity and ease of implementation. However, a longstanding limitation of the two-stream approximation is its inaccuracy in the presence of medium-sized or large aerosols. This limitation was lifted by the discovery of the improved two-stream technique, where the accuracy of the scattering greenhouse effect is matched to that of multi-stream calculations by construction. In the present study, we derive the full solutions for improved two-stream radiative transfer, following its introduction by Heng & Kitzmann (2017), and include contributions from the direct stellar beam. The generalization of the original two-stream flux solutions comes in the form of a correction factor, traditionally set to unity, which is the ratio of a pair of first Eddington coefficients. We derive an analytical expression for this correction factor and also provide a simple fitting function for its ease of use by other workers. We prove that the direct stellar beam is associated with a second Eddington coefficient that is on the order of unity. Setting this second Eddington coefficient to 2/3 and $1/\sqrt{3}$ reproduces the Eddington and quadrature closures, respectively, associated with the direct beam. We use our improved two-stream solutions for the fluxes to derive two-stream source function solutions for both the intensity and fluxes.
[13]  oai:arXiv.org:1709.09678  [pdf] - 1678459
Combining low- to high-resolution transit spectroscopy of HD 189733b. Linking the troposphere and the thermosphere of a hot gas giant
Comments: A&A 612, A53 (2018); 13 pages, 8 figures, 2 tables
Submitted: 2017-09-27, last modified: 2018-05-08
Space-borne low-to medium-resolution (R~10^2-10^3) transmission spectroscopy of atmospheres detect the broadest spectral features (alkali doublets, molecular bands, scattering), while high-resolution (R~10^5), ground-based observations probe the sharpest features (cores of the alkali lines, molecular lines).The two techniques differ by:(1) The LSF of ground-based observations is 10^3 times narrower than for space-borne observations;(2)Space-borne transmission spectra probe up to the base of thermosphere, while ground-based observations can reach pressures down to 10^(-11);(3)Space-borne observations directly yield the transit depth of the planet, while ground-based observations measure differences in the radius of the planet at different wavelengths.It is challenging to combine both techniques.We develop a method to compare theoretical models with observations at different resolutions.We introduce PyETA, a line-by-line 1D radiative transfer code to compute transmission spectra at R~10^6 (0.01 A) over a broad wavelength range.An hybrid forward modeling/retrieval optimization scheme is devised to deal with the large computational resources required by modeling a broad wavelength range (0.3-2 $\mu$m) at high resolution.We apply our technique to HD189733b.Here, HST observations reveal a flattened spectrum due to scattering by aerosols, while high-resolution ground-based HARPS observations reveal the sharp cores of sodium lines.We reconcile these results by building models that reproduce simultaneously both data sets, from the troposphere to the thermosphere. We confirm:(1)the presence of scattering by tropospheric aerosols;(2)that the sodium core feature is of thermospheric origin.Accounting for aerosols, the sodium cores indicate T up to 10000K in the thermosphere.The precise value of the thermospheric temperature is degenerate with the abundance of sodium and altitude of the aerosol deck.
[14]  oai:arXiv.org:1802.03047  [pdf] - 1656310
Revisiting the Phase Curves of WASP-43b: Confronting Reanalyzed Spitzer Data with Cloudy Atmospheres
Comments: Accepted for publication in ApJ
Submitted: 2018-02-08
Recently acquired Hubble and Spitzer phase curves of the short-period hot Jupiter WASP-43b make it an ideal target for confronting theory with data. On the observational front, we re-analyze the 3.6 and 4.5 $\mu$m Spitzer phase curves and demonstrate that our improved analysis better removes residual red noise due to intra-pixel sensitivity, which leads to greater fluxes emanating from the nightside of WASP-43b, thus reducing the tension between theory and data. On the theoretical front, we construct cloudfree and cloudy atmospheres of WASP-43b using our Global Circulation Model (GCM), THOR, which solves the non-hydrostatic Euler equations (compared to GCMs that typically solve the hydrostatic primitive equations). The cloudfree atmosphere produces a reasonable fit to the dayside emission spectrum. The multi-phase emission spectra constrain the cloud deck to be confined to the nightside and have a finite cloud-top pressure. The multi-wavelength phase curves are naturally consistent with our cloudy atmospheres, except for the 4.5 $\mu$m phase curve, which requires the presence of enhanced carbon dioxide in the atmosphere of WASP-43b. Multi-phase emission spectra at higher spectral resolution, as may be obtained using the James Webb Space Telescope, and a reflected-light phase curve at visible wavelengths would further constrain the properties of clouds in WASP-43b.
[15]  oai:arXiv.org:1705.04444  [pdf] - 1604213
AGN feedback with the Square Kilometer Array (SKA) and implications for cluster physics and cosmology
Comments: 14 pages, 10 figures, Review article accepted in Journal of Astrophysics and Astronomy (JOAA)
Submitted: 2017-05-12, last modified: 2017-10-17
AGN feedback is regarded as an important non-gravitational process in galaxy clusters, providing useful constraints on large-scale structure formation. It modifies the structure and energetics of the intra-cluster medium (ICM) and hence its understanding is crucially needed in order to use clusters as high precision cosmological probes. In this context, particularly keeping in mind the upcoming high quality radio data expected from radio surveys like SKA with its higher sensitivity, high spatial and spectral resolutions, we review our current understanding of AGN feedback, its cosmological implications and the impact that SKA can have in revolutionizing our understanding of AGN feedback in large-scale structures. Recent developments regarding the AGN outbursts and its possible contribution to excess entropy in the hot atmospheres of groups and clusters, its correlation with the feedback energy in ICM, quenching of cooling flows and the possible connection between cool core clusters and radio mini-halos, are discussed. We describe current major issues regarding modeling of AGN feedback and its impact on the surrounding medium. With regard to the future of AGN feedback studies, we examine the possible breakthroughs that can be expected from SKA observations. In the context of cluster cosmology, for example, we point out the importance of SKA observations for cluster mass calibration by noting that most of $z>1$ clusters discovered by eROSITA X-ray mission can be expected to be followed up through a 1000 hour SKA-1 mid programme. Moreover, approximately $1000$ radio mini halos and $\sim 2500$ radio halos at $z<0.6$ can be potentially detected by SKA1 and SKA2 and used as tracers of galaxy clusters and determination of cluster selection function.
[16]  oai:arXiv.org:1710.02427  [pdf] - 1589389
Aerosol Constraints on the Atmosphere of the Hot Saturn-mass planet WASP-49b
Comments: Accepted for publication in ApJ
Submitted: 2017-10-06
The strong, nearly wavelength-independent absorption cross section of aerosols produces featureless exoplanet transmission spectra, limiting our ability to characterize their atmospheres. Here we show that even in the presence of featureless spectra, we can still characterize certain atmospheric properties. Specifically, we constrain the upper and lower pressure boundaries of aerosol layers, and present plausible composition candidates. We study the case of the bloated Saturn-mass planet WASP-49b, where near-infrared observations reveal a flat transmission spectrum between 0.7 and 1.0 {\microns}. First, we use a hydrodynamic upper-atmosphere code to estimate the pressure reached by the ionizing stellar high-energy photons at $10^{-8}$ bar, setting the upper pressure boundary where aerosols could exist. Then, we combine HELIOS and Pyrat Bay radiative-transfer models to constrain the temperature and photospheric pressure of atmospheric aerosols, in a Bayesian framework. For WASP-49b, we constrain the transmission photosphere (hence, the aerosol deck boundaries) to pressures above $10^{-5}$ bar (100$\times$ solar metallicity), $10^{-4}$ bar (solar), and $10^{-3}$ bar (0.1$\times$ solar) as lower boundary, and below $10^{-7}$ bar as upper boundary. Lastly, we compare condensation curves of aerosol compounds with the planet's pressure-temperature profile to identify plausible condensates responsible for the absorption. Under these circumstances, we find as candidates: Na$_{2}$S (at 100$\times$ solar metallicity); Cr and MnS (at solar and 0.1$\times$ solar); and forsterite, enstatite, and alabandite (at 0.1$\times$ solar).
[17]  oai:arXiv.org:1709.00338  [pdf] - 1587807
Retrieval Analysis of the Emission Spectrum of WASP-12b: Sensitivity of Outcomes to Prior Assumptions and Implications for Formation History
Comments: Accepted by ApJL. 7 pages, 5 figures
Submitted: 2017-09-01
We analyze the emission spectrum of the hot Jupiter WASP-12b using our HELIOS-R retrieval code and HELIOS-K opacity calculator. When interpreting Hubble and Spitzer data, the retrieval outcomes are found to be prior-dominated. When the prior distributions of the molecular abundances are assumed to be log-uniform, the volume mixing ratio of HCN is found to be implausibly high. A VULCAN chemical kinetics model of WASP-12b suggests that chemical equilibrium is a reasonable assumption even when atmospheric mixing is implausibly rigorous. Guided by (exo)planet formation theory, we set Gaussian priors on the elemental abundances of carbon, oxygen and nitrogen with the Gaussian peaks being centered on the measured C/H, O/H and N/H values of the star. By enforcing chemical equilibrium, we find substellar O/H and stellar to slightly superstellar C/H for the dayside atmosphere of WASP-12b. The superstellar carbon-to-oxygen ratio is just above unity, regardless of whether clouds are included in the retrieval analysis, consistent with Madhusudhan et al. (2011). Furthermore, whether a temperature inversion exists in the atmosphere depends on one's assumption for the Gaussian width of the priors. Our retrieved posterior distributions are consistent with the formation of WASP-12b in a solar-composition protoplanetary disk, beyond the water iceline, via gravitational instability or pebble accretion (without core erosion) and migration inwards to its present orbital location via a disk-free mechanism, and are inconsistent with both in-situ formation and core accretion with disk migration, as predicted by Madhusudhan et al. (2017). We predict that the interpretation of James Webb Space Telescope WASP-12b data will not be prior-dominated.
[18]  oai:arXiv.org:1610.03216  [pdf] - 1580376
HELIOS-Retrieval: An Open-source, Nested Sampling Atmospheric Retrieval Code, Application to the HR 8799 Exoplanets and Inferred Constraints for Planet Formation
Comments: 27 pages, 21 figures, 3 tables, published in AJ
Submitted: 2016-10-11, last modified: 2017-08-21
We present an open-source retrieval code named HELIOS-Retrieval (hereafter HELIOS-R), designed to obtain chemical abundances and temperature-pressure profiles from inverting the measured spectra of exoplanetary atmospheres. In the current implementation, we use an exact solution of the radiative transfer equation, in the pure absorption limit, in our forward model, which allows us to analytically integrate over all of the outgoing rays (instead of performing Gaussian quadrature). Two chemistry models are considered: unconstrained chemistry (where the mixing ratios are treated as free parameters) and equilibrium chemistry (enforced via analytical formulae, where only the elemental abundances are free parameters). The nested sampling algorithm allows us to formally implement Occam's Razor based on a comparison of the Bayesian evidence between models. We perform a retrieval analysis on the measured spectra of the HR 8799b, c, d and e directly imaged exoplanets. Chemical equilibrium is disfavored by the Bayesian evidence for HR 8799b, c and d. We find supersolar C/O, C/H and O/H values for the outer HR 8799b and c exoplanets, while the inner HR 8799d and e exoplanets have substellar C/O, substellar C/H and superstellar O/H values. If these retrieved properties are representative of the bulk compositions of the exoplanets, then they are inconsistent with formation via gravitational instability (without late-time accretion) and consistent with a core accretion scenario in which late-time accretion of ices occurred differently for the inner and outer exoplanets. For HR 8799e, we find that spectroscopy in the K band is crucial for constraining C/O and C/H. HELIOS-R is publicly available as part of the Exoclimes Simulation Platform (ESP; www.exoclime.org).
[19]  oai:arXiv.org:1703.00028  [pdf] - 1581687
Excess entropy and energy feedback from within cluster cores up to r$_{200}$
Comments: 17 pages, 15 figures, 5 tables, Accepted in MNRAS
Submitted: 2017-02-28, last modified: 2017-08-03
We estimate the "non-gravitational" entropy-injection profiles, $\Delta K$, and the resultant energy feedback profiles, $\Delta E$, of the intracluster medium for 17 clusters using their Planck SZ and ROSAT X-Ray observations, spanning a large radial range from $0.2r_{500}$ up to $r_{200}$. The feedback profiles are estimated by comparing the observed entropy, at fixed gas mass shells, with theoretical entropy profiles predicted from non-radiative hydrodynamic simulations. We include non-thermal pressure and gas clumping in our analysis. The inclusion of non-thermal pressure and clumping results in changing the estimates for $r_{500}$ and $r_{200}$ by 10\%-20\%. When clumpiness is not considered it leads to an under-estimation of $\Delta K\approx300$ keV cm$^2$ at $r_{500}$ and $\Delta K\approx1100$ keV cm$^2$ at $r_{200}$. On the other hand, neglecting non-thermal pressure results in an over-estimation of $\Delta K\approx 100$ keV cm$^2$ at $r_{500}$ and under-estimation of $\Delta K\approx450$ keV cm$^2$ at $r_{200}$. For the estimated feedback energy, we find that ignoring clumping leads to an under-estimation of energy per particle $\Delta E\approx1$ keV at $r_{500}$ and $\Delta E\approx1.5$ keV at $r_{200}$. Similarly, neglect of the non-thermal pressure results in an over-estimation of $\Delta E\approx0.5$ keV at $r_{500}$ and under-estimation of $\Delta E\approx0.25$ keV at $r_{200}$. We find entropy floor of $\Delta K\approx300$ keV cm$^2$ is ruled out at $\approx3\sigma$ throughout the entire radial range and $\Delta E\approx1$ keV at more than 3$\sigma$ beyond $r_{500}$, strongly constraining ICM pre-heating scenarios. We also demonstrate robustness of results w.r.t sample selection, X-Ray analysis procedures, entropy modeling etc.
[20]  oai:arXiv.org:1610.05776  [pdf] - 1560179
Low-$\ell$ power suppression in punctuated inflation
Comments: 18 pages, 6 figures, 3 tables, Accepted in JCAP
Submitted: 2016-10-18, last modified: 2017-04-03
Motivated by Planck confirmation of an anomalously low value of the CMB temperature fluctuations up to multipole $\ell<40$, we in this paper try to explain such feature by investigating case of punctuated inflation scenario. This form of inflation potential is inspired by Minimal Super-symmetric Standard Model (MSSM) wherein suppression of curvature perturbation power at large scales is produced by introducing period of fast-roll phase of the inflation sandwiched between two stages of slow-roll phase. We apply Markov Chain Monte Carlo analysis to determine posterior distribution and the best fit values of the model parameters using recent WMAP9 and Planck data. We show that WMAP9 and Planck results are consistent with each other and that with Planck data we obtain tighter constraints for punctuated inflation parameters. We find that punctuated inflation leads to better fit in CMB data compared to simple power law model. The improvement in the fit to the WMAP9 data is $\Delta \chi^2\sim 3.6$ and for Planck the improvement is $\Delta \chi^2\sim 5.4$. We find that $AIC$ does not discriminate between punctuated inflation and simple power law model for WMAP9 data. However, for Planck data we find that punctuated inflation is moderately preferred over a simple power law model.
[21]  oai:arXiv.org:1606.05474  [pdf] - 1530855
HELIOS: An Open-source, GPU-accelerated Radiative Transfer Code For Self-consistent Exoplanetary Atmospheres
Comments: 24 pages, 14 figures, accepted for publication in ApJ
Submitted: 2016-06-17, last modified: 2016-11-19
We present the open-source radiative transfer code named HELIOS, which is constructed for studying exoplanetary atmospheres. In its initial version, the model atmospheres of HELIOS are one-dimensional and plane-parallel, and the equation of radiative transfer is solved in the two-stream approximation with non-isotropic scattering. A small set of the main infrared absorbers is employed, computed with the opacity calculator HELIOS-K and combined using a correlated-$k$ approximation. The molecular abundances originate from validated analytical formulae for equilibrium chemistry. We compare HELIOS with the work of Miller-Ricci & Fortney using a model of GJ 1214b, and perform several tests, where we find: model atmospheres with single-temperature layers struggle to converge to radiative equilibrium; $k$-distribution tables constructed with $\gtrsim 0.01$ cm$^{-1}$ resolution in the opacity function ($ \lesssim 10^3$ points per wavenumber bin) may result in errors $\gtrsim 1$-10 % in the synthetic spectra; and a diffusivity factor of 2 approximates well the exact radiative transfer solution in the limit of pure absorption. We construct "null-hypothesis" models (chemical equilibrium, radiative equilibrium and solar element abundances) for 6 hot Jupiters. We find that the dayside emission spectra of HD 189733b and WASP-43b are consistent with the null hypothesis, while it consistently under-predicts the observed fluxes of WASP-8b, WASP-12b, WASP-14b and WASP-33b. We demonstrate that our results are somewhat insensitive to the choice of stellar models (blackbody, Kurucz or PHOENIX) and metallicity, but are strongly affected by higher carbon-to-oxygen ratios. The code is publicly available as part of the Exoclimes Simulation Platform (ESP; exoclime.net).
[22]  oai:arXiv.org:1606.00014  [pdf] - 1522507
Little evidence for entropy and energy excess beyond $r_{500}$ - An end to ICM preheating?
Comments: 6 pages, 4 figures, 1 table, Accepted in MNRAS Letters
Submitted: 2016-05-31, last modified: 2016-10-27
Non-gravitational feedback affects the nature of the intra-cluster medium (ICM). X-ray cooling of the ICM and in situ energy feedback from AGN's and SNe as well as {\it preheating} of the gas at epochs preceding the formation of clusters are proposed mechanisms for such feedback. While cooling and AGN feedbacks are dominant in cluster cores, the signatures of a preheated ICM are expected to be present even at large radii. To estimate the degree of preheating, with minimum confusion from AGN feedback/cooling, we study the excess entropy and non-gravitational energy profiles upto $r_{200}$ for a sample of 17 galaxy clusters using joint data sets of {\it Planck} SZ pressure and {\it ROSAT/PSPC} gas density profiles. The canonical value of preheating entropy floor of $\gtrsim 300$ keV cm$^2$, needed in order to match cluster scalings, is ruled out at $\approx 3\sigma$. We also show that the feedback energy of 1 keV/particle is ruled out at 5.2$\sigma$ beyond $r_{500}$. Our analysis takes both non-thermal pressure and clumping into account which can be important in outer regions. Our results based on the direct probe of the ICM in the outermost regions do not support any significant preheating.
[23]  oai:arXiv.org:1501.02647  [pdf] - 1126203
Joint Planck and WMAP Assessment of Low CMB Multipoles
Comments: 27 pages, 10 figures, 3 tables, matches with the published version, abstract is shortened to keep it within arXiv's limit (1920 characters)
Submitted: 2015-01-12, last modified: 2015-06-08
The remarkable progress in cosmic microwave background (CMB) studies over past decade has led to the era of precision cosmology in striking agreement with the $\Lambda$CDM model. However, the lack of power in the CMB temperature anisotropies at large angular scales (low-$\ell$), as has been confirmed by the recent Planck data also (up to $\ell=40$), although statistically not very strong (less than $3\sigma$), is still an open problem. One can avoid to seek an explanation for this problem by attributing the lack of power to cosmic variance orcan look for explanations i.e., different inflationary potentials or initial conditions for infl ation to begin with, non-trivial topology, ISW effect etc. Features in the primordial power spectrum (PPS) motivated by the early universe physics has been the most common solution to address this problem. In the present work we also follow this approach and consider a set of PPS which have features and constrain the parameters of those using WMAP 9 year and Planck data employing Markov-Chain Monte Carlo (MCMC) analysis. The prominent feature of all the models of PPS that we consider is an infra-red cut off which leads to suppression of power at large angular scales. We consider models of PPS with maximum three extra parameters and use Akaike information criterion ($AIC$) and Bayesian information criterion ($BIC$) of model and Bayesian information criterion ($BIC$) of model selection to compare the models. For most models, we find good constraints for the cut off scale $k_c$, however, for other parameters our constraints are not that good. We find that sharp cut off model gives best likelihood value for the WMAP 9 year data, but is as good as power law model according to $AIC$. For the joint WMAP 9+Planck data set, Starobinsky model is slightly preferred by $AIC$ which is also able to produce CMB power suppression up to $\ell\leq30$ to some extent.
[24]  oai:arXiv.org:1502.06597  [pdf] - 1231836
On the gap-opening criterion of migrating planets in protoplanetary disks
Comments: Accepted by ApJ, 22 pages, 13 figures, 6 tables
Submitted: 2015-02-23
We perform two-dimensional hydrodynamical simulations to quantitatively explore the torque balance criterion for gap-opening (as formulated by Crida et al. 2006) in a variety of disks when considering a migrating planet. We find that even when the criterion is satisfied, there are instances when planets still do not open gaps. We stress that gap-opening is not only dependent on whether a planet has the ability to open a gap, but whether it can do so quickly enough. This can be expressed as an additional condition on the gap-opening timescale versus the crossing time, i.e. the time it takes the planet to cross the region which it is carving out. While this point has been briefly made in the previous literature, our results quantify it for a range of protoplanetary disk properties and planetary masses, demonstrating how crucial it is for gap-opening. This additional condition has important implications for the survival of planets formed by core accretion in low mass disks as well as giant planets or brown dwarfs formed by gravitational instability in massive disks. It is particularly important for planets with intermediate masses susceptible to Type III-like migration. For some observed transition disks or disks with gaps, we expect that estimates on the potential planet masses based on the torque balance gap-opening criterion alone may not be sufficient. With consideration of this additional timescale criterion theoretical studies may find a reduced planet survivability or that planets may migrate further inwards before opening a gap.
[25]  oai:arXiv.org:1304.3125  [pdf] - 1165885
Debris discs around M stars: non-existence versus non-detection
Comments: Accepted by MNRAS. 11 pages, 9 figures
Submitted: 2013-04-10
Motivated by the reported dearth of debris discs around M stars, we use survival models to study the occurrence of planetesimal discs around them. These survival models describe a planetesimal disc with a small number of parameters, determine if it may survive a series of dynamical processes and compute the associated infrared excess. For the WISE satellite, we demonstrate that the dearth of debris discs around M stars may be attributed to the small semi-major axes generally probed if either: 1. the dust grains behave like blackbodies emitting at a peak wavelength coincident with the observed one; 2. or the grains are hotter than predicted by their blackbody temperatures and emit at peak wavelengths that are shorter than the observed one. At these small distances from the M star, planetesimals are unlikely to survive or persist for time scales of 300 Myr or longer if the disc is too massive. Conversely, our survival models allow for the existence of a large population of low-mass debris discs that are too faint to be detected with current instruments. However, our interpretation becomes less clear and large infrared excesses are allowed if only one of these scenarios holds: 3. the dust grains are hotter than blackbody and predominantly emit at the observed wavelength; 4. or are blackbody in nature and emit at peak wavelengths longer than the observed one. Both scenarios imply that the parent planetesimals reside at larger distances from the star than inferred if the dust grains behaved like blackbodies. In all scenarios, we show that the infrared excesses detected at 22 and 70 microns from AU Mic are easily reconciled with its young age. We elucidate the conditions under which stellar wind drag may be neglected when considering dust populations around M stars. The WISE satellite should be capable of detecting debris discs around young M stars with ages on the order of 10 Myr.
[26]  oai:arXiv.org:1011.5170  [pdf] - 1042156
MARVELS-1b: A Short-Period, Brown Dwarf Desert Candidate from the SDSS-III MARVELS Planet Search
Comments: 47 pages, 11 figures. Accepted for publication in ApJ
Submitted: 2010-11-23
We present a new short-period brown dwarf candidate around the star TYC 1240-00945-1. This candidate was discovered in the first year of the Multi-object APO Radial Velocity Exoplanets Large-area Survey (MARVELS), which is part of the third phase of the Sloan Digital Sky Survey (SDSS-III), and we designate the brown dwarf as MARVELS-1b. MARVELS uses the technique of dispersed fixed-delay interferometery to simultaneously obtain radial velocity measurements for 60 objects per field using a single, custom-built instrument that is fiber fed from the SDSS 2.5-m telescope. From our 20 radial velocity measurements spread over a ~370 d time baseline, we derive a Keplerian orbital fit with semi-amplitude K=2.533+/-0.025 km/s, period P=5.8953+/-0.0004 d, and eccentricity consistent with circular. Independent follow-up radial velocity data confirm the orbit. Adopting a mass of 1.37+/-0.11 M_Sun for the slightly evolved F9 host star, we infer that the companion has a minimum mass of 28.0+/-1.5 M_Jup, a semimajor axis 0.071+/-0.002 AU assuming an edge-on orbit, and is probably tidally synchronized. We find no evidence for coherent instrinsic variability of the host star at the period of the companion at levels greater than a few millimagnitudes. The companion has an a priori transit probability of ~14%. Although we find no evidence for transits, we cannot definitively rule them out for companion radii ~<1 R_Jup.