Identification of habitable planets beyond our solar system is a key goal of current and future space missions. Yet habitability depends not only on the stellar irradiance, but equally on constituent parts of the planetary atmosphere. Here we show, for the first time, that radiatively active mineral dust will have a significant impact on the habitability of Earth-like exoplanets. On tidally-locked planets, dust cools the day-side and warms the night-side, significantly widening the habitable zone. Independent of orbital configuration, we suggest that airborne dust can postpone planetary water loss at the inner edge of the habitable zone, through a feedback involving decreasing ocean coverage and increased dust loading. The inclusion of dust significantly obscures key biomarker gases (e.g. ozone, methane) in simulated transmission spectra, implying an important influence on the interpretation of observations. We demonstrate that future observational and theoretical studies of terrestrial exoplanets must consider the effect of dust.

A key factor in determining the potential habitability of synchronously rotating planets is the strength of the atmospheric boundary layer inversion between the dark side surface and the free atmosphere. Here we analyse data obtained from polar night measurements at the South Pole and Alert Canada, which are the closest analogues on Earth to conditions on the dark sides of synchronously rotating exoplanets without and with a maritime influence, respectively. On Earth, such inversions rarely exceed 30 K in strength, because of the effect of turbulent mixing induced by phenomena such as so-called mesoscale slope winds, which have horizontal scales of 10s to 100s of km, suggesting a similar constraint to near-surface dark side inversions. We discuss the sensitivity of inversion strength to factors such as orography and the global-scale circulation, and compare them to a simulation of the planet Proxima Centauri b. Our results demonstrate the importance of comparisons with Earth data in exoplanet research, and highlight the need for further studies of the exoplanet atmospheric collapse problem using mesoscale and eddy-resolving models.

We report on a multi-band variability and correlation study of the TeV blazar Mrk 421 during an exceptional flaring activity observed from 2013 April 11 to 2013 April 19. The study uses, among others, data from GASP-WEBT, Swift, NuSTAR, Fermi-LAT, VERITAS, and MAGIC. The large blazar activity, and the 43 hours of simultaneous NuSTAR and MAGIC/VERITAS observations, permitted variability studies on 15 minute time bins, and over three X-ray bands (3-7 keV, 7-30 keV and 30-80 keV) and three very-high-energy (>0.1 TeV, hereafter VHE) gamma-ray bands (0.2-0.4 TeV, 0.4-0.8 TeV and >0.8 TeV). We detected substantial flux variations on multi-hour and sub-hour timescales in all the X-ray and VHE gamma-ray bands. The characteristics of the sub-hour flux variations are essentially energy-independent, while the multi-hour flux variations can have a strong dependence on the energy of the X-ray and the VHE gamma rays. The three VHE bands and the three X-ray bands are positively correlated with no time-lag, but the strength and the characteristics of the correlation changes substantially over time and across energy bands. Our findings favour multi-zone scenarios for explaining the achromatic/chromatic variability of the fast/slow components of the light curves, as well as the changes in the flux-flux correlation on day-long timescales. We interpret these results within a magnetic reconnection scenario, where the multi-hour flux variations are dominated by the combined emission from various plasmoids of different sizes and velocities, while the sub-hour flux variations are dominated by the emission from a single small plasmoid moving across the magnetic reconnection layer.

We report the results of decade-long (2008-2018) $\gamma$-ray to 1 GHz radio monitoring of the blazar 3C 279, including GASP/WEBT, $\it{Fermi}$ and $\it{Swift}$ data, as well as polarimetric and spectroscopic data. The X-ray and $\gamma$-ray light curves correlate well, with no delay > 3 hours, implying general co-spatiality of the emission regions. The $\gamma$-ray-optical flux-flux relation changes with activity state, ranging from a linear to a more complex dependence. The behaviour of the Stokes parameters at optical and radio wavelengths, including 43 GHz VLBA images, supports either a predominantly helical magnetic field or motion of the radiating plasma along a spiral path. Apparent speeds of emission knots range from 10 to 37c, with the highest values requiring bulk Lorentz factors close to those needed to explain $\gamma$-ray variability on very short time scales. The Mg II emission line flux in the `blue' and `red' wings correlates with the optical synchrotron continuum flux density, possibly providing a variable source of seed photons for inverse Compton scattering. In the radio bands we find progressive delays of the most prominent light curve maxima with decreasing frequency, as expected from the frequency dependence of the $\tau=1$ surface of synchrotron self-absorption. The global maximum in the 86 GHz light curve becomes less prominent at lower frequencies, while a local maximum, appearing in 2014, strengthens toward decreasing frequencies, becoming pronounced at $\sim5$ GHz. These tendencies suggest different Doppler boosting of stratified radio-emitting zones in the jet.

We present a multiwavelength study of the flat-spectrum radio quasar CTA 102 during 2013-2017. We use radio-to-optical data obtained by the Whole Earth Blazar Telescope, 15 GHz data from the Owens Valley Radio Observatory, 91 and 103 GHz data from the Atacama Large Millimeter Array, near-infrared data from the Rapid Eye Monitor telescope, as well as data from the Swift (optical-UV and X-rays) and Fermi ($\gamma$ rays) satellites to study flux and spectral variability and the correlation between flux changes at different wavelengths. Unprecedented $\gamma$-ray flaring activity was observed during 2016 November-2017 February, with four major outbursts. A peak flux of (2158 $\pm$ 63)$\times$10$^{-8}$ ph cm$^{-2}$ s$^{-1}$, corresponding to a luminosity of (2.2 $\pm$ 0.1)$\times$10$^{50}$ erg s$^{-1}$, was reached on 2016 December 28. These four $\gamma$-ray outbursts have corresponding events in the near-infrared, optical, and UV bands, with the peaks observed at the same time. A general agreement between X-ray and $\gamma$-ray activity is found. The $\gamma$-ray flux variations show a general, strong correlation with the optical ones with no time lag between the two bands and a comparable variability amplitude. This $\gamma$-ray/optical relationship is in agreement with the geometrical model that has successfully explained the low-energy flux and spectral behaviour, suggesting that the long-term flux variations are mainly due to changes in the Doppler factor produced by variations of the viewing angle of the emitting regions. The difference in behaviour between radio and higher energy emission would be ascribed to different viewing angles of the jet regions producing their emission.

The energy balance and climate of planets can be affected by the reflective properties of their land, ocean, and frozen surfaces. Here we investigate the effect of host star spectral energy distribution (SED) on the albedo of these surfaces using a one-dimensional (1-D) energy balance model (EBM). Incorporating spectra of M-, K-, G- and F-dwarf stars, we determined the effect of varying fractional and latitudinal distribution of land and ocean surfaces as a function of host star SED on the overall planetary albedo, climate, and ice-albedo feedback response. While noting that the spatial distribution of land masses on a given planet will have an effect on the overall planetary energy balance, we find that terrestrial planets with higher average land/ocean fractions are relatively cooler and have higher albedo regardless of star type. For Earth-like planets orbiting M-dwarf stars the increased absorption of water ice in the near-infrared (NIR), where M-dwarf stars emit much of their energy, resulted in warmer global mean surface temperatures, ice lines at higher latitudes, and increased climate stability as the ice-albedo feedback became negative at high land fractions. Conversely, planets covered largely by ocean, and especially those orbiting bright stars, had a considerably different energy balance due to the contrast between the reflective land and the absorptive ocean surface, which in turn resulted in warmer average surface temperatures than land-covered planets and a stronger potential ice-albedo feedback. While dependent on the properties of individual planetary systems, our results place so constraints on a range of climate states of terrestrial exoplanets based on albedo and incident flux.

We present results of a search for variable stars in a region of the globular cluster NGC 4147 based on photometric observations with 4k x 4k CCD imager mounted at the axial port of the recently installed 3.6 m Devasthal optical telescope at Aryabhatta Research Institute of Observational Sciences, Nainital, India. We performed time series photometry of NGC 4147 in V and R bands, and identified 42 periodic variables in the region of NGC 4147, 28 of which have been detected for the first time. Seventeen variable stars are located within the half light radius $\lesssim$ 0.48 arcmin, of which 10 stars are newly identified variables. Two of 10 variables are located within the core radius $\lesssim$ 0.09 arcmin. Based on the location in the $V/(V-R)$ colour magnitude diagram and variability characteristics, 7, 8, 5 and 1 newly identified probable member variables are classified as RRc, EA/E, EW and SX Phe, respectively. The metallicity of NGC 4147 estimated from light curves of RRab and RRc stars with the help of Fourier decomposition is found to be characteristics of Oosterhoff II. The distance derived using light curves of RRab stars is consistent with that obtained from the observed $V/(V-R)$ colour-magnitude diagram.

We present the results of extensive modeling of the spectral energy distributions (SEDs) of the closest blazar (z=0.031) Markarian 421 (Mrk 421) during a giant outburst in February 2010. The source underwent rapid flux variations in both X-rays and very high energy (VHE) gamma-rays as it evolved from a low-flux state on 2010 February 13-15 to a high-flux state on 2010 February 17. During this period, the source exhibited significant spectral hardening from X-rays to VHE gamma-rays while exhibiting a "harder when brighter" behavior in these energy bands. We reproduce the broadband SED using a time-dependent multi-zone leptonic jet model with radiation feedback. We find that an injection of the leptonic particle population with a single power-law energy distribution at shock fronts followed by energy losses in an inhomogeneous emission region is suitable for explaining the evolution of Mrk 421 from low- to high-flux state in February 2010. The spectral states are successfully reproduced by a combination of a few key physical parameters, such as the maximum $\&$ minimum cutoffs and power-law slope of the electron injection energies, magnetic field strength, and bulk Lorentz factor of the emission region. The simulated light curves and spectral evolution of Mrk 421 during this period imply an almost linear correlation between X-ray flux at 1-10 keV energies and VHE gamma-ray flux above 200 GeV, as has been previously exhibited by this source. Through this study, a general trend that has emerged for the role of physical parameters is that, as the flare evolves from a low- to a high-flux state, higher bulk kinetic energy is injected into the system with a harder particle population and a lower magnetic field strength.

We compare the rise and decay timescales of $\sim$200 long-term ($\sim$weeks-months) GeV and R-band outbursts and $\sim$25 short-term ($\sim$hr-day) GeV flares in a sample of 10 blazars using light curves from the Fermi-LAT and the Yale/SMARTS monitoring project. We find that most of the long-term outbursts are symmetric, indicating that the observed variability is dominated by the crossing timescale of a disturbance, e.g., a shock. A larger fraction of short-term flares are asymmetric with an approximately equal fraction of longer and shorter decay than rise timescale. We employ the MUlti-ZOne Radiation Feedback (MUZORF) model to interpret the above results. We find that the outbursts with slow rise times indicate a gradual acceleration of the particles to GeV energy. A change in the bulk Lorentz factor of the plasma or the width of the shocked region can lead to an increase of the cooling time causing a faster rise than decay time. Parameters such as the luminosity or the distance of the broad line region (BLR) affects the cooling time strongly if a single emission mechanism, e.g., external Compton scattering of BLR photons is considered but may not if other mechanisms, e.g., synchrotron self-compton and external Compton scattering of the torus photon are included. This work carries out a systematic study of the symmetry of flares, which can be used to estimate relevant geometric and physical parameters of blazar jets in the context of the MUZORF model.

We analyze the linear polarization of the relativistic jet in BL Lacertae object OJ~287 as revealed by multi-epoch Very Long Baseline Array (VLBA) images at 43 GHz and monitoring observations at optical bands. The electric-vector position angle (EVPA) of the optical polarization matches that at 43 GHz at locations that are often in the compact millimeter-wave "core" or, at other epochs, coincident with a bright, quasi-stationary emission feature $\sim0.2$~milliarcsec ($\sim$0.9~pc projected on the sky) downstream from the core. This implies that electrons with high enough energies to emit optical synchrotron and $\gamma$-ray inverse Compton radiation are accelerated both in the core and at the downstream feature, the latter of which lies $\geq10$~pc from the central engine. The polarization vector in the stationary feature is nearly parallel to the jet axis, as expected for a conical standing shock capable of accelerating electrons to GeV energies.

We present the results of a study of the time-averaged spectral energy distributions (SEDs) of eight flat spectrum radio quasars (FSRQs) present in the second catalog of high energy sources detected beyond 50 GeV by the {\it Fermi}-Large Area Telescope (2FHL). Both leptonic and hadronic scenarios are adopted to explain the multi-wavelength SEDs and we find them to be marginally consistent with the 2FHL spectra above 50 GeV. We derive the expected degree of X-ray and $\gamma$-ray polarizations both for the average and elevated activity states and note that: (i) a hadronic radiative model consistently predicts a higher degree of high energy polarization compared to leptonic ones, and (ii) the X-ray polarization degree is higher than the $\gamma$-ray polarization in the leptonic scenario, but similar to the $\gamma$-ray polarization if the observed radiation is powered by hadronic processes. From the leptonic modeling, the location of the $\gamma$-ray emitting region is found to be at the outer edge of the broad line region (BLR) and is consistent with the $\gamma\gamma$ opacity estimates for the $\gamma$-ray absorption by the BLR. We conclude that a majority of the FSRQs could be detected by the upcoming Cherenkov Telescope Array, though future high energy polarimeters will be able to detect them only during elevated activity states, which could provide supportive evidence for the hadronic origin of the X-ray and $\gamma$-ray emission.

We investigate the atmospheric dynamics of terrestrial planets in synchronous rotation within the habitable zone of low-mass stars using the Community Atmosphere Model (CAM). The surface temperature contrast between day and night hemispheres decreases with an increase in incident stellar flux, which is opposite the trend seen on gas giants. We define three dynamical regimes in terms of the equatorial Rossby deformation radius and the Rhines length. The slow rotation regime has a mean zonal circulation that spans from day to night side, with both the Rossby deformation radius and the Rhines length exceeding planetary radius, which occurs for planets around stars with effective temperatures of 3300 K to 4500 K (rotation period > 20 days). Rapid rotators have a mean zonal circulation that partially spans a hemisphere and with banded cloud formation beneath the substellar point, with the Rossby deformation radius is less than planetary radius, which occurs for planets orbiting stars with effective temperatures of less than 3000 K (rotation period < 5 days). In between is the Rhines rotation regime, which retains a thermally-direct circulation from day to night side but also features midlatitude turbulence-driven zonal jets. Rhines rotators occur for planets around stars in the range of 3000 K to 3300 K (rotation period ~ 5 to 20 days), where the Rhines length is greater than planetary radius but the Rossby deformation radius is less than planetary radius. The dynamical state can be observationally inferred from comparing the morphology of the thermal emission phase curves of synchronously rotating planets.

We present the results of a multi-frequency, time-averaged analysis of blazars included in the Candidate Gamma-ray Blazar Survey catalog. Our sample consists of 324 $\gamma$-ray detected ($\gamma$-ray loud) and 191 non $\gamma$-ray detected ($\gamma$-ray quiet) blazars, and we consider all the data up to 2016 April 1. We find that both the $\gamma$-ray loud and the $\gamma$-ray quiet blazar populations occupy similar regions in the WISE color-color diagram, and in the radio and X-ray bands $\gamma$-ray loud sources are brighter. A simple one-zone synchrotron inverse-Compton emission model is applied to derive the physical properties of both populations. We find that the central black hole mass and the accretion disk luminosity ($L_{\rm disk}$) computed from the modeling of the optical-UV emission with a Shakura-Sunyaev disk reasonably matches with that estimated from the optical spectroscopic emission-line information. A significantly larger Doppler boosting in the $\gamma$-ray loud blazars is noted, and their jets are more radiatively efficient. On the other hand, the $\gamma$-ray quiet objects are more MeV-peaked, thus could be potential targets for next-generation MeV missions. Our results confirm the earlier findings about the accretion-jet connection in blazars; however, many of the $\gamma$-ray quiet blazars tend to deviate from the recent claim that the jet power exceeds $L_{\rm disk}$ in blazars. A broadband study, considering a larger set of $\gamma$-ray quiet objects and also including BL Lacs, will be needed to confirm/reject this hypothesis and also to verify the evolution of the powerful high-redshift blazars into their low-power nearby counterparts.

The knowledge of the structure of the magnetic field inside a blazar jet, as deduced from polarization observations at radio to optical wavelengths, is closely related to the formation and propagation of relativistic jets that result from accretion onto supermassive black holes. However, a largely unexplored aspect of the theoretical understanding of radiation transfer physics in blazar jets has been the magnetic field geometry as revealed by the polarized emission and the connection between the variability in polarization and flux across the spectrum. Here, we explore the effects of various magnetic geometries that can exist inside a blazar jet: parallel, oblique, toroidal, and tangled. We investigate the effects of changing the orientation of the magnetic field, according to the above-mentioned geometries, on the resulting high-energy spectral energy distributions (SEDs) and spectral variability patterns (SVPs) of a typical blazar. We use the MUlti-ZOne Radiation Feedback (MUZORF) model of Joshi et al. (2014) to carry out this study and to relate the geometry of the field to the observed SEDs at X-ray and gamma-ray energies. One of the goals of the study is to understand the relationship between synchrotron and inverse Compton peaks in blazar SEDs and the reason for the appearance of gamma-ray "orphan flares" observed in some blazars. This can be associated with the directionality of the magnetic field, which creates a difference in the radiation field as seen by an observer versus that seen by the electrons in the emission region.

Over the past few years, several occasions of large, continuous rotations of the electric vector position angle (EVPA) of linearly polarized optical emission from blazars have been reported. These events are often coincident with high energy gamma-ray flares and they have attracted considerable attention, as they could allow one to probe the magnetic field structure in the gamma-ray emitting region of the jet. The flat-spectrum radio quasar 3C279 is one of the most prominent examples showing this behaviour. Our goal is to study the observed EVPA rotations and to distinguish between a stochastic and a deterministic origin of the polarization variability. We have combined multiple data sets of R-band photometry and optical polarimetry measurements of 3C279, yielding exceptionally well-sampled flux density and polarization curves that cover a period of 2008-2012. Several large EVPA rotations are identified in the data. We introduce a quantitative measure for the EVPA curve smoothness, which is then used to test a set of simple random walk polarization variability models against the data. 3C279 shows different polarization variation characteristics during an optical low-flux state and a flaring state. The polarization variation during the flaring state, especially the smooth approx. 360 degrees rotation of the EVPA in mid-2011, is not consistent with the tested stochastic processes. We conclude that during the two different optical flux states, two different processes govern the polarization variation, possibly a stochastic process during the low-brightness state and a deterministic process during the flaring activity.

We present coordinated multiwavelength observations of the bright, nearby BL Lac object Mrk 421 taken in 2013 January-March, involving GASP-WEBT, Swift, NuSTAR, Fermi-LAT, MAGIC, VERITAS, and other collaborations and instruments, providing data from radio to very-high-energy (VHE) gamma-ray bands. NuSTAR yielded previously unattainable sensitivity in the 3-79 keV range, revealing that the spectrum softens when the source is dimmer until the X-ray spectral shape saturates into a steep power law with a photon index of approximately 3, with no evidence for an exponential cutoff or additional hard components up to about 80 keV. For the first time, we observed both the synchrotron and the inverse-Compton peaks of the spectral energy distribution (SED) simultaneously shifted to frequencies below the typical quiescent state by an order of magnitude. The fractional variability as a function of photon energy shows a double-bump structure which relates to the two bumps of the broadband SED. In each bump, the variability increases with energy which, in the framework of the synchrotron self-Compton model, implies that the electrons with higher energies are more variable. The measured multi-band variability, the significant X-ray-to-VHE correlation down to some of the lowest fluxes ever observed in both bands, the lack of correlation between optical/UV and X-ray flux, the low degree of polarization and its significant (random) variations, the short estimated electron cooling time, and the significantly longer variability timescale observed in the NuSTAR light curves point toward in-situ electron acceleration, and suggest that there are multiple compact regions contributing to the broadband emission of Mrk 421 during low-activity states.

We report on simultaneous broadband observations of the TeV-emitting blazar Markarian 501 between 1 April and 10 August 2013, including the first detailed characterization of the synchrotron peak with Swift and NuSTAR. During the campaign, the nearby BL Lac object was observed in both a quiescent and an elevated state. The broadband campaign includes observations with NuSTAR, MAGIC, VERITAS, the Fermi Large Area Telescope (LAT), Swift X-ray Telescope and UV Optical Telescope, various ground-based optical instruments, including the GASP-WEBT program, as well as radio observations by OVRO, Mets\"ahovi and the F-Gamma consortium. Some of the MAGIC observations were affected by a sand layer from the Saharan desert, and had to be corrected using event-by-event corrections derived with a LIDAR (LIght Detection And Ranging) facility. This is the first time that LIDAR information is used to produce a physics result with Cherenkov Telescope data taken during adverse atmospheric conditions, and hence sets a precedent for the current and future ground-based gamma-ray instruments. The NuSTAR instrument provides unprecedented sensitivity in hard X-rays, showing the source to display a spectral energy distribution between 3 and 79 keV consistent with a log-parabolic spectrum and hard X-ray variability on hour timescales. None (of the four extended NuSTAR observations) shows evidence of the onset of inverse-Compton emission at hard X-ray energies. We apply a single-zone equilibrium synchrotron self-Compton model to five simultaneous broadband spectral energy distributions. We find that the synchrotron self-Compton model can reproduce the observed broadband states through a decrease in the magnetic field strength coinciding with an increase in the luminosity and hardness of the relativistic leptons responsible for the high-energy emission.

We perform a multi-wavelength polarimetric study of the quasar CTA 102 during an extraordinarily bright $\gamma$-ray outburst detected by the {\it Fermi} Large Area Telescope in September-October 2012 when the source reached a flux of F$_{>100~\mathrm{MeV}} =5.2\pm0.4\times10^{-6}$ photons cm$^{-2}$ s$^{-1}$. At the same time the source displayed an unprecedented optical and NIR outburst. We study the evolution of the parsec scale jet with ultra-high angular resolution through a sequence of 80 total and polarized intensity Very Long Baseline Array images at 43 GHz, covering the observing period from June 2007 to June 2014. We find that the $\gamma$-ray outburst is coincident with flares at all the other frequencies and is related to the passage of a new superluminal knot through the radio core. The powerful $\gamma$-ray emission is associated with a change in direction of the jet, which became oriented more closely to our line of sight ($\theta\sim$1.2$^{\circ}$) during the ejection of the knot and the $\gamma$-ray outburst. During the flare, the optical polarized emission displays intra-day variability and a clear clockwise rotation of EVPAs, which we associate with the path followed by the knot as it moves along helical magnetic field lines, although a random walk of the EVPA caused by a turbulent magnetic field cannot be ruled out. We locate the $\gamma$-ray outburst a short distance downstream of the radio core, parsecs from the black hole. This suggests that synchrotron self-Compton scattering of near-infrared to ultraviolet photons is the probable mechanism for the $\gamma$-ray production.

Blazars exhibit flares across the electromagnetic spectrum. Many $\gamma$-ray flares are highly correlated with flares detected at optical wavelengths; however, a small subset appears to occur in isolation, with little or no variability detected at longer wavelengths. These "orphan" $\gamma$-ray flares challenge current models of blazar variability, most of which are unable to reproduce this type of behavior. We present numerical calculations of the time-variable emission of a blazar based on a proposal by Marscher et al. (2010) to explain such events. In this model, a plasmoid ("blob") propagates relativistically along the spine of a blazar jet and passes through a synchrotron-emitting ring of electrons representing a shocked portion of the jet sheath. This ring supplies a source of seed photons that are inverse-Compton scattered by the electrons in the moving blob. The model includes the effects of radiative cooling, a spatially varying magnetic field, and acceleration of the blob's bulk velocity. Synthetic light curves produced by our model are compared to the observed light curves from an orphan flare that was coincident with the passage of a superluminal knot through the inner jet of the blazar PKS 1510$-$089. In addition, we present Very Long Baseline Array polarimetric observations that point to the existence of a jet sheath in PKS 1510$-$089, thus providing further observational support for the plausibility of our model. An estimate of the bolometric luminosity of the sheath within PKS 1510$-$089 is made, yielding $L_{\rm sh} \sim 3 \times 10^{45} ~ \rm erg ~ \rm s^{-1}$. This indicates that the sheath within PKS 1510$-$089 is potentially a very important source of seed photons.

The blazar 1156+295 was active at gamma-ray energies, exhibiting three prominent flares during the year 2010. Here, we present results using the combination of broadband (X-ray through mm single dish) monitoring data and radio band imaging data at 43 GHz on the connection of gamma-ray events to the ejections of superluminal components and other changes in the jet of 1156+295. The kinematics of the jet over the interval 2007.0-2012.5 using 43 GHz Very Long Baseline Array observations, reveal the presence of four moving and one stationary component in the inner region of the blazar jet. The propagation of the third and fourth components in the jet corresponds closely in time to the active phase of the source in gamma rays. We briefly discuss the implications of the structural changes in the jet for the mechanism of gamma-ray production during bright flares. To localise the gamma-ray emission site in the blazar, we performed the correlation analysis between the 43 GHz radio core and the gamma-ray light curve. The time lag obtained from the correlation constrains the gamma-ray emitting region in the parsec-scale jet.

We present $\gamma$-ray, X-ray, ultraviolet, optical, and near-infrared light curves of 33 $\gamma$-ray bright blazars over four years that we have been monitoring since 2008 August with multiple optical, ground-based telescopes and the Swift satellite, and augmented by data from the Fermi Gamma-ray Space Telescope and other publicly available data from Swift. The sample consists of 21 flat-spectrum radio quasars (FSRQs) and 12 BL Lac objects (BL Lacs). We identify quiescent and active states of the sources based on their $\gamma$-ray behavior. We derive $\gamma$-ray, X-ray, and optical spectral indices, $\alpha_\gamma$, $\alpha_X$, and $\alpha_o$, respectively ($F_\nu\propto\nu^\alpha$), and construct spectral energy distributions (SEDs) during quiescent and active states. We analyze the relationships between different spectral indices, blazar classes, and activity states. We find (i) significantly steeper $\gamma$-ray spectra of FSRQs than for BL Lacs during quiescent states, but a flattening of the spectra for FSRQs during active states while the BL Lacs show no significant change; (ii) a small difference of $\alpha_X$ within each class between states, with BL Lac X-ray spectra significantly steeper than in FSRQs; (iii) a highly peaked distribution of X-ray spectral slopes of FSRQs at $\sim-$0.60, but a very broad distribution of $\alpha_X$ of BL Lacs during active states; (iv) flattening of the optical spectra of FSRQs during quiescent states, but no statistically significant change of $\alpha_o$ of BL Lacs between states; and (v) a positive correlation between optical and $\gamma$-ray spectral slopes of BL Lacs, with similar values of the slopes. We discuss the findings with respect to the relative prominence of different components of high-energy and optical emission as the flux state changes.

We present the results of optical (R band) photometric and polarimetric monitoring and Very Long Baseline Array (VLBA) imaging of the blazar S4 0954+658, along with Fermi and gamma;-ray data during a multi-waveband outburst in 2011 March-April. After a faint state with a brightness level R ~17.6 mag registered in the first half of January 2011, the optical brightness of the source started to rise and reached ~14.8 mag during the middle of March, showing flare-like behavior. The most spectacular case of intranight variability was observed during the night of 2011 March 9, when the blazar brightened by ~0.7 mag within ~7 hours. During the rise of the flux the position angle of optical polarization rotated smoothly over more than 300$\deg$. At the same time, within 1$\sigma$ uncertainty a new superluminal knot appeared with an apparent speed of 19.0$\pm$0.3 c. We have very strong evidence for association of this knot with the multi-waveband outburst in 2011 March-April. We also analyze the multi-frequency behavior of S4 0954+658 during a number of minor outbursts from August 2008 to April 2012. We find some evidence of connections between at least two more superluminal ejecta and near-simultaneous optical flares.

Distant planets in globally ice-covered, "snowball", states may depend on increases in their host stars' luminosity to become hospitable for surface life. Using a General Circulation Model (GCM), we simulated the equilibrium climate response of a planet to a range of instellations from an F-, G-, or M- dwarf star. The range of instellation that permits both complete ice cover and at least partially ice-free climate states is a measure of the climate hysteresis that a planet can exhibit. An ice-covered planet with high climate hysteresis would show a higher resistance to the initial loss of surface ice coverage with increases in instellation, and abrupt, extreme ice loss once deglaciation begins. Our simulations indicate that the climate hysteresis depends sensitively on the host star spectral energy distribution. Under fixed CO2 conditions, a planet orbiting an M-dwarf star exhibits a smaller climate hysteresis, requiring a smaller instellation to initiate deglaciation than planets orbiting hotter, brighter stars. This is due to the higher absorption of near-IR radiation by ice on the surfaces and greenhouse gases and clouds in the atmosphere of an M-dwarf planet. Increases in atmospheric CO2 further lower the climate hysteresis, as M-dwarf snowball planets exhibit a larger radiative response than G-dwarf snowball planets for the same increase in CO2. For a smaller hysteresis, planets near the outer edge of the habitable zone will thaw earlier in their evolutionary history, and will experience a less abrupt transition out of global ice cover.

We extend our approach of modeling spectral energy distribution (SED) and lightcurves of blazars to include external Compton (EC) emission due to inverse Compton scattering of an external anisotropic target radiation field. We describe the time-dependent impact of such seed photon fields on the evolution of multifrequency emission and spectral variability of blazars using a multi-zone time-dependent leptonic jet model, with radiation feedback, in the internal shock model scenario. We calculate accurate EC-scattered high-energy spectra produced by relativistic electrons throughout the Thomson and Klein-Nishina regimes. We explore the effects of varying the contribution of (1) a thermal Shakura-Sunyaev accretion disk, (2) a spherically symmetric shell of broad-line clouds, the broad line region (BLR), and (3) a hot infrared emitting dusty torus (DT), on the resultant seed photon fields. We let the system evolve to beyond the BLR and within the DT and study the manifestation of the varying target photon fields on the simulated SED and lightcurves of a typical blazar. The calculations of broadband spectra include effects of gamma-gamma absorption as gamma-rays propagate through the photon pool present inside the jet due to synchrotron and inverse Compton processes, but neglect gamma-gamma absorption by the BLR and DT photon fields outside the jet. Thus, our account of gamma-gamma absorption is a lower limit to this effect. Here, we focus on studying the impact of parameters relevant for EC processes on high-energy (HE) emission of blazars.

Quasar 3C 279 is known to exhibit episodes of optical polarization angle rotation. We present new, well-sampled optical polarization data for 3C 279 and introduce a method to distinguish between random and deterministic electric vector position angle (EVPA) variations. We observe EVPA rotations in both directions with different amplitudes and find that the EVPA variation shows characteristics of both random and deterministic cases. Our analysis indicates that the EVPA variation is likely dominated by a random process in the low brightness state of the jet and by a deterministic process in the flaring state.

We analyze the multifrequency behavior of the quasar 3C 454.3 during three prominent \gamma-ray outbursts: 2009 Autumn, 2010 Spring, and 2010 Autumn. The data reveal a repeating pattern, including a triple flare structure, in the properties of each \gamma-ray outburst, which implies similar mechanism(s) and location for all three events. The multi-frequency behavior indicates that the lower frequency events are co-spatial with the \gamma-ray outbursts, although the \gamma-ray emission varies on the shortest timescales. We determine that the variability from UV to IR wavelengths during an outburst results from a single synchrotron component whose properties do not change significantly over the different outbursts. Despite a general increase in the degree of optical linear polarization during an outburst, the polarization drops significantly at the peak of the \gamma-ray event, which suggests that both shocks and turbulent processes are involved. We detect two disturbances (knots) with superluminal apparent speeds in the parsec-scale jet associated with the outbursts in 2009 Autumn and 2010 Autumn. The kinematic properties of the knots can explain the difference in amplitudes of the \gamma-ray events, while their millimeter-wave polarization is related to the optical polarization during the outbursts. We interpret the multi-frequency behavior within models involving either a system of standing conical shocks or magnetic reconnection events located in the parsec-scale millimeter-wave core of the jet. We argue that \gamma-ray outbursts with variability timescales as short as ~ 3 hr can occur on parsec scales if flares take place in localized regions such as turbulent cells.

Planetary climate can be affected by the interaction of the host star spectral energy distribution with the wavelength-dependent reflectivity of ice and snow. Here we explore this effect using a one dimensional (1-D), line-by-line, radiative-transfer model to calculate broadband planetary albedos as input to a seasonally varying, 1-D energy-balance climate model. A three-dimensional general circulation model is also used to explore the atmosphere's response to changes in incoming stellar radiation, or instellation, and surface albedo. Using this hierarchy of models we simulate planets covered by ocean, land, and water ice of varying grain size, with incident radiation from stars of different spectral types. Terrestrial planets orbiting stars with higher near-UV radiation exhibit a stronger ice-albedo feedback. We find that ice-covered conditions occur on an F-dwarf planet with only a 2% reduction in instellation relative to the present instellation on Earth, assuming fixed CO2 (present atmospheric level on Earth). A similar planet orbiting the Sun at an equivalent flux distance requires an 8% reduction in instellation, while a planet orbiting an M-dwarf star requires an additional 19% reduction in instellation to become ice-covered, equivalent to 73% of the modern solar constant. The surface ice-albedo feedback effect becomes less important at the outer edge of the habitable zone, where atmospheric CO2 can be expected to be high in order to maintain clement conditions for surface liquid water. We show that 3-10 bars of CO2 will entirely mask the climatic effect of ice and snow, leaving the outer limits of the habitable zone unaffected by the spectral dependence of water ice and snow albedo. However, less CO2 is needed to maintain open water for a planet orbiting an M-dwarf star, than would be the case for hotter main-sequence stars.

We present multiwavelength data of the blazar 3C 454.3 obtained during an extremely bright outburst from November 2010 through January 2011. These include flux density measurements with the Herschel Space Observatory at five submillimeter-wave and far-infrared bands, the Fermi Large Area Telescope at gamma-ray energies, Swift at X-ray, ultraviolet (UV), and optical frequencies, and the Submillimeter Array at 1.3 mm. From this dataset, we form a series of 52 spectral energy distributions (SEDs) spanning nearly two months that are unprecedented in time coverage and breadth of frequency. Discrete correlation anlaysis of the millimeter, far-infrared, and gamma-ray light curves show that the variations were essentially simultaneous, indicative of co-spatiality of the emission, at these wavebands. In contrast, differences in short-term fluctuations at various wavelengths imply the presence of inhomegeneities in physical conditions across the source. We locate the site of the outburst in the parsec-scale core, whose flux density as measured on 7 mm Very Long Baseline Array images increased by 70 percent during the first five weeks of the outburst. Based on these considerations and guided by the SEDs, we propose a model in which turbulent plasma crosses a conical standing shock in the parsec-scale region of the jet. Here, the high-energy emission in the model is produced by inverse Compton scattering of seed photons supplied by either nonthermal radiation from a Mach disk, thermal emission from hot dust, or (for X-rays) synchrotron radiation from plasma that crosses the standing shock. For the two dates on which we fitted the model SED to the data, the model corresponds very well to the observations at all bands except at X-ray energies, where the spectrum is flatter than observed.

The quasar-type blazar 4C 38.41 (B3 1633+382) experienced a large outburst in 2011, which was detected throughout the entire electromagnetic spectrum. We present the results of low-energy multifrequency monitoring by the GASP project of the WEBT consortium and collaborators, as well as those of spectropolarimetric/spectrophotometric monitoring at the Steward Observatory. We also analyse high-energy observations of the Swift and Fermi satellites. In the optical-UV band, several results indicate that there is a contribution from a QSO-like emission component, in addition to both variable and polarised jet emission. The unpolarised emission component is likely thermal radiation from the accretion disc that dilutes the jet polarisation. We estimate its brightness to be R(QSO) ~ 17.85 - 18 and derive the intrinsic jet polarisation degree. We find no clear correlation between the optical and radio light curves, while the correlation between the optical and \gamma-ray flux apparently fades in time, likely because of an increasing optical to \gamma-ray flux ratio. As suggested for other blazars, the long-term variability of 4C 38.41 can be interpreted in terms of an inhomogeneous bent jet, where different emitting regions can change their alignment with respect to the line of sight, leading to variations in the Doppler factor \delta. Under the hypothesis that in the period 2008-2011 all the \gamma-ray and optical variability on a one-week timescale were due to changes in \delta, this would range between ~ 7 and ~ 21. If the variability were caused by changes in the viewing angle \theta\ only, then \theta\ would go from ~ 2.6 degr to ~ 5 degr. Variations in the viewing angle would also account for the dependence of the polarisation degree on the source brightness in the framework of a shock-in-jet model.

We have observed 3C~279 in a gamma-ray flaring state in November 2008. We construct quasi-simultaneous spectral energy distributions (SEDs) of the source for the flaring period of 2008 and during a quiescent period in May 2010. Data have been compiled from observations with Fermi, Swift, RXTE, the VLBA, and various ground-based optical and radio telescopes. The objective is to comprehend the correspondence between the flux and polarization variations observed during these two time periods by carrying out a detailed spectral analyses of 3C~279 in the internal shock scenario, and gain insights into the role of intrinsic parameters and interplay of synchrotron and inverse Compton radiation processes responsible for the two states. As a first step, we have used a multi-slice time-dependent leptonic jet model, in the framework of the internal shock scenario, with radiation feedback to simulate the SED of 3C~279 observed in an optical high state in early 2006. We have used physical jet parameters obtained from the VLBA monitoring to guide our modeling efforts. We briefly discuss the effects of intrinsic parameters and various radiation processes in producing the resultant SED.

We present time-resolved broad-band observations of the quasar 3C 279 obtained from multi-wavelength campaigns conducted during the first two years of the Fermi Gamma-ray Space Telescope mission. While investigating the previously reported gamma-ray/optical flare accompanied by a change in optical polarization, we found that the optical emission appears delayed with respect to the gamma-ray emission by about 10 days. X-ray observations reveal a pair of `isolated' flares separated by ~90 days, with only weak gamma-ray/optical counterparts. The spectral structure measured by Spitzer reveals a synchrotron component peaking in the mid-infrared band with a sharp break at the far-infrared band during the gamma-ray flare, while the peak appears in the mm/sub-mm band in the low state. Selected spectral energy distributions are fitted with leptonic models including Comptonization of external radiation produced in a dusty torus or the broad-line region. Adopting the interpretation of the polarization swing involving propagation of the emitting region along a curved trajectory, we can explain the evolution of the broad-band spectra during the gamma-ray flaring event by a shift of its location from ~ 1 pc to ~ 4 pc from the central black hole. On the other hand, if the gamma-ray flare is generated instead at sub-pc distance from the central black hole, the far-infrared break can be explained by synchrotron self-absorption. We also model the low spectral state, dominated by the mm/sub-mm peaking synchrotron component, and suggest that the corresponding inverse-Compton component explains the steady X-ray emission.

We analyze total and polarized intensity images of the quasar 3C 454.3 obtained monthly with the VLBA at 43 GHz within the ongoing Boston U. monitoring program of gamma-ray blazars started in June 2007. The data are supplemented by VLBA observations performed during intense campaigns of 2 week duration when the quasar was observed 3 times per campaign. We find a strong increase of activity in the parsec-scale jet of the quasar during high gamma-ray states in December 2009, April 2010, and November 2010. We detect new superluminal knots, K09 and K10, associated with the autumn 2009 and 2010 outbursts, respectively, and compare their kinematic parameters. We analyze optical polarimetric behavior along with polarization parameters of the parsec-scale jet and outline similarities and differences in polarization properties across wavelengths. The results of the analysis support the conclusions that the optical polarized emission is produced in a region located in the vicinity of the mm-wave core of the jet of the quasar, and that the gamma-ray outbursts occur when a superluminal disturbance passes through the core.

M-stars comprise 80% of main-sequence stars, and so their planetary systems provide the best chance for finding habitable planets, i.e.: those with surface liquid water. We have modelled the broadband albedo or reflectivity of water ice and snow for simulated planetary surfaces orbiting two observed red dwarf stars (or M-stars) using spectrally resolved data of the Earth's cryosphere. The gradual reduction of the albedos of snow and ice at wavelengths greater than 1 ?m, combined with M-stars emitting a significant fraction of their radiation at these same longer wavelengths, mean that the albedos of ice and snow on planets orbiting M-stars are much lower than their values on Earth. Our results imply that the ice/snow albedo climate feedback is significantly weaker for planets orbiting M-stars than for planets orbiting G-type stars such as the Sun. In addition, planets with significant ice and snow cover will have significantly higher surface temperatures for a given stellar flux if the spectral variation of cryospheric albedo is considered, which in turn implies that the outer edge of the habitable zone around M-stars may be 10-30% further away from the parent star than previously thought.

We present multiwavelength spectral analyses of two Fermi-LAT blazars, OJ 287 and 3C 279, that are part of the Boston University multiwaveband polarization program. The data have been compiled from observations with Fermi, RXTE, the VLBA, and various ground-based optical and radio telescopes. We simulate the dynamic spectral energy distributions (SEDs) within the framework of a multi-slice, time-dependent leptonic jet model for blazars, with radiation feedback, in the internal shock scenario. We use the physical jet parameters obtained from the VLBA monitoring to guide our modeling efforts. We discuss the role of intrinsic parameters and the interplay between synchrotron and inverse Compton radiation processes responsible for producing the resultant SEDs.

Since 2005, the blazar 3C 454.3 has shown remarkable flaring activity at all frequencies, and during the last four years it has exhibited more than one gamma-ray flare per year, becoming the most active gamma-ray blazar in the sky. We present for the first time the multi-wavelength AGILE, SWIFT, INTEGRAL, and GASP-WEBT data collected in order to explain the extraordinary gamma-ray flare of 3C 454.3 which occurred in November 2010. On 2010 November 20 (MJD 55520), 3C 454.3 reached a peak flux (E>100 MeV) of F_gamma(p) = (6.8+-1.0)E-5 ph/cm2/s on a time scale of about 12 hours, more than a factor of 6 higher than the flux of the brightest steady gamma-ray source, the Vela pulsar, and more than a factor of 3 brighter than its previous super-flare on 2009 December 2-3. The multi-wavelength data make a thorough study of the present event possible: the comparison with the previous outbursts indicates a close similarity to the one that occurred in 2009. By comparing the broadband emission before, during, and after the gamma-ray flare, we find that the radio, optical and X-ray emission varies within a factor 2-3, whereas the gamma-ray flux by a factor of 10. This remarkable behavior is modeled by an external Compton component driven by a substantial local enhancement of soft seed photons.

We present observations of a major outburst at centimeter, millimeter, optical, X-ray, and gamma-ray wavelengths of the BL Lacertae object AO 0235+164. We analyze the timing of multi-waveband variations in the flux and linear polarization, as well as changes in Very Long Baseline Array (VLBA) images at 7mm with 0.15 milliarcsecond resolution. The association of the events at different wavebands is confirmed at high statistical significance by probability arguments and Monte-Carlo simulations. A series of sharp peaks in optical linear polarization, as well as a pronounced maximum in the 7 mm polarization of a superluminal jet knot, indicate rapid fluctuations in the degree of ordering of the magnetic field. These results lead us to conclude that the outburst occurred in the jet both in the quasi-stationary "core" and in the superluminal knot, both parsecs downstream of the supermassive black hole. We interpret the outburst as a consequence of the propagation of a disturbance, elongated along the line of sight by light-travel time delays, that passes through a standing recollimation shock in the core and propagates down the jet to create the superluminal knot. The multi-wavelength light curves vary together on long time-scales (months/years), but the correspondence is poorer on shorter time-scales. This, as well as the variability of the polarization and the dual location of the outburst, agrees with the expectations of a multi-zone emission model in which turbulence plays a major role in modulating the synchrotron and inverse Compton fluxes.

We describe the time-dependent radiation transfer in blazar jets, within the internal shock model. We assume that the central engine, which consists of a black hole and an accretion disk, spews out relativistic shells of plasma with different velocity, mass, and energy. We consider a single inelastic collision between a faster (inner) and a slower (outer) moving shell. We study the dynamics of the collision and evaluate the subsequent emission of radiation via the synchrotron and synchrotron self Compton (SSC) processes after the interaction between the two shells has begun. The collision results in the formation of a forward shock (FS) and a reverse shock (RS) that convert the ordered bulk kinetic energy of the shells into magnetic field energy and accelerate the particles, which then radiate. We assume a cylindrical geometry for the emission region of the jet. We treat the self-consistent radiative transfer by taking into account the inhomogeneity in the photon density throughout the region. In this paper, we focus on understanding the effects of varying relevant input parameters on the simulated spectral energy distribution (SED) and spectral variability patterns.

We present results of five years of optical (UBVRI) observations of the very-high-energy gamma-ray blazar 1ES 1011+496 at the MDM Observatory. We calibrated UBVRI magnitudes of five comparison stars in the field of the object. Most of our observations were done during moderately faint states of 1ES 1011+496 with R > 15.0. The light curves exhibit moderate, closely correlated variability in all optical wavebands on time scales of a few days. A cross-correlation analysis between optical bands does not show significant evidence for time lags. We find a positive correlation (Pearson's r = 0.57; probability of non-correlation P(>r) ~ 4e-8) between the R-band magnitude and the B - R color index, indicating a bluer-when-brighter trend. Snap-shot optical spectral energy distributions (SEDs) exhibit a peak within the optical regime, typically between the V and B bands. We find a strong (r = 0.78; probability of non-correlation P (>r) ~ 1e-15) positive correlation between the peak flux and the peak frequency, best fit by a relation $\nu F_{\nu}^{\rm pk} \propto \nu_{\rm pk}^k$ with k = 2.05 +/- 0.17. Such a correlation is consistent with the optical (synchrotron) variability of 1ES 1011+496 being primarily driven by changes in the magnetic field.

The ZEPLIN-III experiment in the Palmer Underground Laboratory at Boulby uses a 12kg two-phase xenon time projection chamber to search for the weakly interacting massive particles (WIMPs) that may account for the dark matter of our Galaxy. The detector measures both scintillation and ionisation produced by radiation interacting in the liquid to differentiate between the nuclear recoils expected from WIMPs and the electron recoil background signals down to ~10keV nuclear recoil energy. An analysis of 847kg.days of data acquired between February 27th 2008 and May 20th 2008 has excluded a WIMP-nucleon elastic scattering spin-independent cross-section above 8.1x10(-8)pb at 55GeV/c2 with a 90% confidence limit. It has also demonstrated that the two-phase xenon technique is capable of better discrimination between electron and nuclear recoils at low-energy than previously achieved by other xenon-based experiments.

The BL Lac object 3C 66A was observed in an extensive multiwavelength monitoring campaign from July 2003 till April 2004. The spectral energy distribution (SED) was measured over the entire electromagnetic spectrum, with flux measurements from radio to X-ray frequencies and upper limits in the very high energy (VHE) gamma-ray regime. Here, we use a time-dependent leptonic jet model to reproduce the SED and optical spectral variability observed during our multiwavelength campaign. Our model simulations could successfully reproduce the observed SED and optical light curves and predict an intrinsic cutoff value for the VHE gamma-ray emission at ~ 4 GeV. The effect of the optical depth due to the intergalactic infrared background radiation (IIBR) on the peak of the high-energy component of 3C 66A was found to be negligible. Also, the presence of a broad line region (BLR) in the case of 3C 66A may play an important role in the emission of gamma-ray photons when the emission region is very close to the central engine, but further out, the production mechanism of hard X-ray and gamma-ray photons becomes rapidly dominated by synchrotron self-Compton emission. We further discuss the possibility of an observable X-ray spectral variability pattern. The simulated results do not predict observable hysteresis patterns in the optical or soft X-ray regimes for major flares on multi-day time scales.

Results are presented from the first underground data run of ZEPLIN-II, a 31 kg two phase xenon detector developed to observe nuclear recoils from hypothetical weakly interacting massive dark matter particles. Discrimination between nuclear recoils and background electron recoils is afforded by recording both the scintillation and ionisation signals generated within the liquid xenon, with the ratio of these signals being different for the two classes of event. This ratio is calibrated for different incident species using an AmBe neutron source and Co-60 gamma-ray sources. From our first 31 live days of running ZEPLIN-II, the total exposure following the application of fiducial and stability cuts was 225 kgxdays. A background population of radon progeny events was observed in this run, arising from radon emission in the gas purification getters, due to radon daughter ion decays on the surfaces of the walls of the chamber. An acceptance window, defined by the neutron calibration data, of 50% nuclear recoil acceptance between 5 keVee and 20 keVee, had an observed count of 29 events, with a summed expectation of 28.6+/-4.3 gamma-ray and radon progeny induced background events. These figures provide a 90% c.l. upper limit to the number of nuclear recoils of 10.4 events in this acceptance window, which converts to a WIMP-nucleon spin-independent cross-section with a minimum of 6.6x10^-7 pb following the inclusion of an energy dependent, calibrated, efficiency. A second run is currently underway in which the radon progeny will be eliminated, thereby removing the background population, with a projected sensitivity of 2x10^-7 pb for similar exposures as the first run.

Stable, hydrogen-burning, M dwarf stars comprise about 75% of all stars in the Galaxy. They are extremely long-lived and because they are much smaller in mass than the Sun (between 0.5 and 0.08 MSun), their temperature and stellar luminosity are low and peaked in the red. We have re-examined what is known at present about the potential for a terrestrial planet forming within, or migrating into, the classic liquid-surface-water habitable zone close to an M dwarf star. Observations of protoplanetary disks suggest that planet-building materials are common around M dwarfs, but N-body simulations differ in their estimations of the likelihood of potentially-habitable, wet planets residing within their habitable zones, which are only ~ 1/5 to 1/50 of the width of that for a G star. Particularly in light of the claimed detection of the planets with masses as small as 5.5 and 7.5 MEarth orbiting M stars, there seems no reason to exclude the possibility of terrestrial planets. Tidally locked synchronous rotation within the narrow habitable zone doesn't necessarily lead to atmospheric collapse, and active stellar flaring may not be as much of an evolutionarily disadvantageous factor as has previously been supposed. We conclude that M dwarf stars may indeed be viable hosts for planets on which the origin and evolution of life can occur. A number of planetary processes such as cessation of geothermal activity, or thermal and non-thermal atmospheric loss processes may limit the duration of planetary habitability to periods far shorter than the extreme lifetime of the M dwarf star. Nevertheless, it makes sense to include M dwarf stars in programs that seek to find habitable worlds and evidence of life.

We present details of the technical design and manufacture of the ZEPLIN-III dark matter experiment. ZEPLIN-III is a two-phase xenon detector which measures both the scintillation light and the ionisation charge generated in the liquid by interacting particles and radiation. The instrument design is driven by both the physics requirements and by the technology requirements surrounding the use of liquid xenon. These include considerations of key performance parameters, such as the efficiency of scintillation light collection, restrictions placed on the use of materials to control the inherent radioactivity levels, attainment of high vacuum levels and chemical contamination control. The successful solution has involved a number of novel design and manufacturing features which will be of specific use to future generations of direct dark matter search experiments as they struggle with similar and progressively more demanding requirements.

We present results from a GEANT4-based Monte Carlo tool for end-to-end simulations of the ZEPLIN-III dark matter experiment. ZEPLIN-III is a two-phase detector which measures both the scintillation light and the ionisation charge generated in liquid xenon by interacting particles and radiation. The software models the instrument response to radioactive backgrounds and calibration sources, including the generation, ray-tracing and detection of the primary and secondary scintillations in liquid and gaseous xenon, and subsequent processing by data acquisition electronics. A flexible user interface allows easy modification of detector parameters at run time. Realistic datasets can be produced to help with data analysis, an example of which is the position reconstruction algorithm developed from simulated data. We present a range of simulation results confirming the original design sensitivity of a few times $10^{-8}$ pb to the WIMP-nucleon cross-section.

We present lightcurves of the afterglow of GRB050502A, including very early data at t-t_{GRB} < 60s. The lightcurve is composed of unfiltered ROTSE-IIIb optical observations from 44s to 6h post-burst, R-band MDM observations from 1.6 to 8.4h post-burst, and PAIRITEL J H K_s observations from 0.6 to 2.6h post-burst. The optical lightcurve is fit by a broken power law, where t^{alpha} steepens from alpha = -1.13 +- 0.02 to alpha = -1.44 +- 0.02 at \~5700s. This steepening is consistent with the evolution expected for the passage of the cooling frequency nu_c through the optical band. Even in our earliest observation at 44s post-burst, there is no evidence that the optical flux is brighter than a backward extrapolation of the later power law would suggest. The observed decay indices and spectral index are consistent with either an ISM or a Wind fireball model, but slightly favor the ISM interpretation. The expected spectral index in the ISM interpretation is consistent within 1 sigma with the observed spectral index beta = -0.8 +- 0.1; the Wind interpretation would imply a slightly (~2 sigma) shallower spectral index than observed. A small amount of dust extinction at the source redshift could steepen an intrinsic spectrum sufficiently to account for the observed value of beta. In this picture, the early optical decay, with the peak at or below 4.7e14 Hz at 44s, requires very small electron and magnetic energy partitions from the fireball.

The BL Lac object 3C 66A was the target of an extensive multiwavelength monitoring campaign from July 2003 through April 2004, involving observations throughout the electromagnetic spectrum. Radio, infrared, and optical observations were carried out by the WEBT-ENIGMA collaboration. At higher energies, 3C 66A was observed in X-rays (RXTE), and at very-high-energy (VHE) gamma-rays (STACEE, VERITAS). In addition, the source has been observed with the VLBA at 9 epochs throughout the period September 2003 - December 2004, including 3 epochs contemporaneous with the core campaign.The WEBT campaign revealed microvariability with flux changes of ~ 5 % on time scales as short as \~ 2 hr. The source was in a relatively bright state, with several bright flares on time scales of several days. The SED indicates a $\nu F_{\nu}$ peak in the optical regime. A weak trend of optical spectral hysteresis with a trend of spectral softening throughout both the rising and decaying phases, has been found. The 3 - 10 keV X-ray flux of 3C 66A during the core campaign was historically high and its spectrum very soft, indicating that the low-frequency component of the broadband SED extends beyond ~ 10 keV. No significant X-ray flux and/or spectral variability was detected. STACEE and Whipple observations provided upper flux limits at > 150 GeV and > 390 GeV, respectively. The 22 and 43 GHz data from the 3 VLBA epochs between September 2003 and January 2004 indicate a smooth jet with only moderate internal structure. Evidence for superluminal motion (8.5 +/- 5.6 h^{-1} c) was found in only one out of 6 components. The radial radio brightness profile suggests a magnetic field decay $\propto r^{-1}$ and, thus, a predominantly perpendicular magnetic field.

The status of dark matter searches with inorganic scintillator detectors at Boulby mine is reviewed and the results of tests with a CsI(Tl) crystal are presented. The objectives of the latter experiment were to study anomalous fast events previously observed and to identify ways to remove this background. Clear indications were found that these events were due to surface contamination of crystals by alphas, probably from radon decay. A new array of unencapsulated NaI(Tl) crystals immersed either in liquid paraffin or pure nitrogen atmosphere is under construction at Boulby. Such an approach allows complete control of the surface of the crystals and the ability to remove any surface contamination. First data from the unencapsulated NaI(Tl) do not show the presence of anomalous fast events.