Aerosols are common in the atmospheres of exoplanets across a wide swath of temperatures, masses, and ages. These aerosols strongly impact observations of transmitted, reflected, and emitted light from exoplanets, obfuscating our understanding of exoplanet thermal structure and composition. Knowing the dominant aerosol composition would facilitate interpretations of exoplanet observations and theoretical understanding of their atmospheres. A variety of compositions have been proposed, including metal oxides and sulphides, iron, chromium, sulphur, and hydrocarbons. However, the relative contributions of these species to exoplanet aerosol opacity is unknown. Here we show that the aerosol composition of giant exoplanets observed in transmission is dominated by silicates and hydrocarbons. By constraining an aerosol microphysics model with trends in giant exoplanet transmission spectra, we find that silicates dominate aerosol opacity above planetary equilibrium temperatures of 950 K due to low nucleation energy barriers and high elemental abundances, while hydrocarbon aerosols dominate below 950 K due to an increase in methane abundance. Our results are robust to variations in planet gravity and atmospheric metallicity within the range of most giant transiting exoplanets. We predict that spectral signatures of condensed silicates in the mid-infrared are most prominent for hot (>1600 K), low-gravity (<10 m s$^{-2}$) objects.

We present a study of the stellar host galaxy, CO (1$-$0) molecular gas distribution and AGN emission on 50 to 500 pc-scales of the gravitationally lensed dust-obscured AGN MG J0751+2716 and JVAS B1938+666 at redshifts 3.200 and 2.059, respectively. By correcting for the lensing distortion using a grid-based lens modelling technique, we spatially locate the different emitting regions in the source plane for the first time. Both AGN host galaxies have 300 to 500 pc-scale size and surface brightness consistent with a bulge/pseudo-bulge, and 2 kpc-scale AGN radio jets that are embedded in extended molecular gas reservoirs that are 5 to 20 kpc in size. The CO (1$-$0) velocity fields show structures possibly associated with discs (elongated velocity gradients) and interacting objects (off-axis velocity components). There is evidence for a decrement in the CO (1$-$0) surface brightness at the location of the host galaxy, which may indicate radiative feedback from the AGN, or offset star formation.We find CO-H$_2$ conversion factors of around $\alpha_{\rm CO} = 1.5\pm0.5$ (K km s$^{-1}$ pc$^2$)$^{-1}$, molecular gas masses of $> 3\times10^{10}$ M$_{\odot}$, dynamical masses of $\sim 10^{11}$ M$_{\odot}$ and gas fractions of around 60 per cent. The intrinsic CO line luminosities are comparable to those of unobscured AGN and dusty star-forming galaxies at similar redshifts, but the infrared luminosities are lower, suggesting that the targets are less efficient at forming stars. Therefore, they may belong to the AGN feedback phase predicted by galaxy formation models, because they are not efficiently forming stars considering their large amount of molecular gas.

We present a new gravitational lens modelling technique designed to model high-resolution interferometric observations with large numbers of visibilities without the need to pre-average the data in time or frequency. We demonstrate the accuracy of the method using validation tests on mock observations. Using small data sets with $\sim 10^3$ visibilities, we first compare our approach with the more traditional direct Fourier transform (DFT) implementation and direct linear solver. Our tests indicate that our source inversion is indistinguishable from that of the DFT. Our method also infers lens parameters to within 1 to 2 per cent of both the ground truth and DFT, given sufficiently high signal-to-noise ratio (SNR). When the SNR is as low as 5, both approaches lead to errors of several tens of per cent in the lens parameters and a severely disrupted source structure, indicating that this is an issue related to the data quality rather than the modelling technique of choice. We then analyze a large data set with $\sim 10^8$ visibilities and a SNR matching real global Very Long Baseline Interferometry observations of the gravitational lens system MG J0751+2716. The size of the data is such that it cannot be modelled with traditional implementations. Using our novel technique, we find that we can infer the lens parameters and the source brightness distribution, respectively, with an RMS error of 0.25 and 0.97 per cent relative to the ground truth.

The James Webb Space Telescope (JWST) is expected to revolutionize the field of exoplanets. The broad wavelength coverage and the high sensitivity of its instruments will allow characterization of exoplanetary atmospheres with unprecedented precision. Following the Call for the Cycle 1 Early Release Science Program, the Transiting Exoplanet Community was awarded time to observe several targets, including WASP-43b. The atmosphere of this hot Jupiter has been intensively observed but still harbors some mysteries, especially concerning the day-night temperature gradient, the efficiency of the atmospheric circulation, and the presence of nightside clouds. We will constrain these properties by observing a full orbit of the planet and extracting its spectroscopic phase curve in the 5--12 $\mu$m range with JWST/MIRI. To prepare for these observations, we performed an extensive modeling work with various codes: radiative transfer, chemical kinetics, cloud microphysics, global circulation models, JWST simulators, and spectral retrieval. Our JWST simulations show that we should achieve a precision of 210 ppm per 0.1 $\mu$m spectral bin on average, which will allow us to measure the variations of the spectrum in longitude and measure the night-side emission spectrum for the first time. If the atmosphere of WASP-43b is clear, our observations will permit us to determine if its atmosphere has an equilibrium or disequilibrium chemical composition, providing eventually the first conclusive evidence of chemical quenching in a hot Jupiter atmosphere. If the atmosphere is cloudy, a careful retrieval analysis will allow us to identify the cloud composition.

We determine the observability in transmission of inhomogeneous cloud cover on the limbs of hot Jupiters through post processing a general circulation model to include cloud distributions computed using a cloud microphysics model. We find that both the east and west limb often form clouds, but that the different properties of these clouds enhances the limb to limb differences compared to the clear case. Using JWST it should be possible to detect the presence of cloud inhomogeneities by comparing the shape of the transit lightcurve at multiple wavelengths because inhomogeneous clouds impart a characteristic, wavelength dependent signature. This method is statistically robust even with limited wavelength coverage, uncertainty on limb darkening coefficients, and imprecise transit times. We predict that the short wavelength slope varies strongly with temperature. The hot limb of the hottest planets form higher altitude clouds composed of smaller particles leading to a strong rayleigh slope. The near infrared spectral features of clouds are almost always detectable, even when no spectral slope is visible in the optical. In some of our models a spectral window between 5 and 9 microns can be used to probe through the clouds and detect chemical spectral features. Our cloud particle size distributions are not log-normal and differ from species to species. Using the area or mass weighted particle size significantly alters the relative strength of the cloud spectral features compared to using the predicted size distribution. Finally, the cloud content of a given planet is sensitive to a species' desorption energy and contact angle, two parameters that could be constrained experimentally in the future.

We present a novel method for determining the surface density of protoplanetary disks through consideration of disk 'dust lines' which indicate the observed disk radial scale at different observational wavelengths. This method relies on the assumption that the processes of particle growth and drift control the radial scale of the disk at late stages of disk evolution such that the lifetime of the disk is equal to both the drift timescale and growth timescale of the maximum particle size at a given dust line. We provide an initial proof of concept of our model through an application to the disk TW Hya and are able to estimate the disk dust-to-gas ratio, CO abundance, and accretion rate in addition to the total disk surface density. We find that our derived surface density profile and dust-to-gas ratio are consistent with the lower limits found through measurements of HD gas. The CO ice line also depends on surface density through grain adsorption rates and drift and we find that our theoretical CO ice line estimates have clear observational analogues. We further apply our model to a large parameter space of theoretical disks and find three observational diagnostics that may be used to test its validity. First we predict that the dust lines of disks other than TW Hya will be consistent with the normalized CO surface density profile shape for those disks. Second, surface density profiles that we derive from disk ice lines should match those derived from disk dust lines. Finally, we predict that disk dust and ice lines will scale oppositely, as a function of surface density, across a large sample of disks.

We present new determinations of disk surface density, independent of an assumed dust opacity, for a sample of 7 bright, diverse protoplanetary disks using measurements of disk dust lines. We develop a robust method for determining the location of dust lines by modeling disk interferometric visibilities at multiple wavelengths. The disks in our sample have newly derived masses that are 9-27% of their host stellar mass, substantially larger than the minimum mass solar nebula. All are stable to gravitational collapse except for one which approaches the limit of Toomre-Q stability. Our mass estimates are 2-15 times larger than estimates from integrated optically thin dust emission. We derive depleted dust-to-gas ratios with typical values of ~$10^{-3}$ in the outer disk. Using coagulation models we derive dust surface density profiles that are consistent with millimeter dust observations. In these models, the disks formed with an initial dust mass that is a factor of ~10 greater than is presently observed. Of the three disks in our sample with resolved CO line emission, the masses of HD 163296, AS 209, and TW Hya are roughly 3, 115, and 40 times more massive than estimates from CO respectively. This range indicates that CO depletion is not uniform across different disks and that dust is a more robust tracer of total disk mass. Our method of determining surface density using dust lines is robust even if particles form as aggregates and is useful even in the presence of dust substructure caused by pressure traps. The low Toomre-Q values observed in this sample indicate that at least some disks do not accrete efficiently.

We present the first application of a bin-scheme microphysical and vertical transport model to determine the size distribution of titanium and silicate cloud particles in the atmospheres of hot Jupiters. We predict particle size distributions from first principles for a grid of planets at four representative equatorial longitudes, and investigate how observed cloud properties depend on the atmospheric thermal structure and vertical mixing. The predicted size distributions are frequently bimodal and irregular in shape. There is a negative correlation between total cloud mass and equilibrium temperature as well as a positive correlation between total cloud mass and atmospheric mixing. The cloud properties on the east and west limbs show distinct differences that increase with increasing equilibrium temperature. Cloud opacities are roughly constant across a broad wavelength range with the exception of features in the mid-infrared. Forward scattering is found to be important across the same wavelength range. Using the fully resolved size distribution of cloud particles as opposed to a mean particle size has a distinct impact on the resultant cloud opacities. The particle size that contributes the most to the cloud opacity depends strongly on the cloud particle size distribution. We predict that it is unlikely that silicate or titanium clouds are responsible for the optical Rayleigh scattering slope seen in many hot Jupiters. We suggest that cloud opacities in emission may serve as sensitive tracers of the thermal state of a planet's deep interior through the existence or lack of a cold trap in the deep atmosphere.

We present a new method for generating initial conditions for numerical cosmological simulations in which massive neutrinos are treated as an extra set of N-body (collisionless) particles. It allows us to accurately follow the density field for both Cold Dark Matter (CDM) and neutrinos at both high and low redshifts. At high redshifts, the new method is able to reduce the shot noise in the neutrino power spectrum by a factor of more than $10^7$ compared to previous methods, where the power spectrum was dominated by shot noise at all scales. We find that our new approach also helps to reduce the noise on the total matter power spectrum on large scales, whereas on small scales the results agree with previous simulations. Our new method also allows for a systematic study of clustering of the low velocity tail of the distribution function of neutrinos. This method also allows for the study of the evolution of the overall velocity distribution as a function of the environment determined by the CDM field.

Indirect detection of dark matter is a major avenue for discovery. However, baryonic backgrounds are diverse enough to mimic many possible signatures of dark matter. In this work, we study the newly proposed technique of dark matter velocity spectroscopy\,\cite{speckhard2016}. The non-rotating dark matter halo and the Solar motion produce a distinct longitudinal dependence of the signal which is opposite in direction to that produced by baryons. Using collisionless dark matter only simulations of Milky Way like halos, we show that this new signature is robust and holds great promise. We develop mock observations by high energy resolution X-ray spectrometer on a sounding rocket, the Micro-X experiment, to our test case, the 3.5 keV line. We show that by using six different pointings, Micro-X can exclude a constant line energy over various longitudes at $\geq$ 3$\sigma$. The halo triaxiality is an important effect and it will typically reduce the significance of this signal. We emphasize that this new {\it smoking gun in motion} signature of dark matter is general, and is applicable to any dark matter candidate which produces a sharp photon feature in annihilation or decay.

We study evolution of voids in cosmological simulations using a new method for tracing voids over cosmic time. The method is based on tracking watershed basins (contiguous regions around density minima) of well developed voids at low redshift, on a regular grid of density field. It enables us to construct a robust and continuous mapping between voids at different redshifts, from initial conditions to the present time. We discuss how the new approach eliminates strong spurious effects of numerical origin when voids evolution is traced by matching voids between successive snapshots (by analogy to halo merger trees). We apply the new method to a cosmological simulation of a standard LambdaCDM cosmological model and study evolution of basic properties of typical voids (with effective radii between 6Mpc/h and 20Mpc/h at redshift z=0) such as volumes, shapes, matter density distributions and relative alignments. The final voids at low redshifts appear to retain a significant part of the configuration acquired in initial conditions. Shapes of voids evolve in a collective way which barely modifies the overall distribution of the axial ratios. The evolution appears to have a weak impact on mutual alignments of voids implying that the present state is in large part set up by the primordial density field. We present evolution of dark matter density profiles computed on iso-density surfaces which comply with the actual shapes of voids. Unlike spherical density profiles, this approach enables us to demonstrate development of theoretically predicted bucket-like shape of the final density profiles indicating a wide flat core and a sharp transition to high-density void walls.

Resolved submillimeter imaging of transitional disks is increasingly revealing the complexity of disk structure. Here we present the first high-resolution submillimeter image of a recently identified transitional disk around IRAS 04125+2902 in the Taurus star-forming region. We measure an inner disk hole of ~20 AU around IRAS 04125+2902 by simultaneously modeling new 880 micron Submillimeter Array (SMA) data along with an existing spectral energy distribution supplemented by new Discovery Channel Telescope (DCT) photometry. We also constrain the outer radius of the dust disk in IRAS~04125+2902 to ~50-60 AU. Such a small dust disk could be attributed to initial formation conditions, outward truncation by an unseen companion, or dust evolution in the disk. Notably, the dust distribution of IRAS 04125+2902 resembles a narrow ring (delta R ~ 35 AU) composed of large dust grains at the location of the disk wall. Such narrow dust rings are also seen in other transitional disks and may be evidence of dust trapping in pressure bumps, possibly produced by planetary companions. More sensitive submillimeter observations of the gas are necessary to further probe the physical mechanisms at work in shaping the spatial distribution of large dust in this disk. Interestingly, the IRAS 04125+2902 disk is significantly fainter than other transitional disks that have been resolved at submillimeter wavelengths, hinting that more objects with large disk holes may exist at the faint end of the submillimeter luminosity distribution that await detection with more sensitive imaging telescopes.

We present an exact general remeshing scheme to compute analytic integrals of polynomial functions over the intersections between convex polyhedral cells of old and new meshes. In physics applications this allows one to ensure global mass, momentum, and energy conservation while applying higher-order polynomial interpolation. We elaborate on applications of our algorithm arising in the analysis of cosmological N-body data, computer graphics, and continuum mechanics problems. We focus on the particular case of remeshing tetrahedral cells onto a Cartesian grid such that the volume integral of the polynomial density function given on the input mesh is guaranteed to equal the corresponding integral over the output mesh. We refer to this as "physically conservative voxelization". At the core of our method is an algorithm for intersecting two convex polyhedra by successively clipping one against the faces of the other. This algorithm is an implementation of the ideas presented abstractly by Sugihara (1994), who suggests using the planar graph representations of convex polyhedra to ensure topological consistency of the output. This makes our implementation robust to geometric degeneracy in the input. We employ a simplicial decomposition to calculate moment integrals up to quadratic order over the resulting intersection domain. We also address practical issues arising in a software implementation, including numerical stability in geometric calculations, management of cancellation errors, and extension to two dimensions. In a comparison to recent work, we show substantial performance gains. We provide a C implementation intended to be a fast, accurate, and robust tool for geometric calculations on polyhedral mesh elements.

We present radio and X-ray observations of the nearby Type IIb Supernova 2013df in NGC4414 from 10 to 250 days after the explosion. The radio emission showed a peculiar soft-to-hard spectral evolution. We present a model in which inverse Compton cooling of synchrotron emitting electrons can account for the observed spectral and light curve evolution. A significant mass loss rate, $\dot{M} \approx 8 \times 10^{-5}\,\rm M_{\odot}/yr$ for a wind velocity of 10 km/s, is estimated from the detailed modeling of radio and X-ray emission, which are primarily due to synchrotron and bremsstrahlung, respectively. We show that SN 2013df is similar to SN 1993J in various ways. The shock wave speed of SN 2013df was found to be average among the radio supernovae; $v_{sh}/c \sim 0.07$. We did not find any significant deviation from smooth decline in the light curve of SN 2013df. One of the main results of our self-consistent multiband modeling is the significant deviation from energy equipartition between magnetic fields and relativistic electrons behind the shock. We estimate $\epsilon_{e} = 200 \epsilon_{B}$. In general for Type IIb SNe, we find that the presence of bright optical cooling envelope emission is linked with free-free radio absorption and bright thermal X-ray emission. This finding suggests that more extended progenitors, similar to that of SN 2013df, suffer from substantial mass loss in the years before the supernova.