We present results from a data challenge posed to the radial velocity (RV) community: namely, to quantify the Bayesian "evidence" for n={0,1,2,3} planets in a set of synthetically generated RV datasets containing a range of planet signals. Participating teams were provided the same likelihood function and set of priors to use in their analysis. They applied a variety of methods to estimate Z, the marginal likelihood for each n-planet model, including cross-validation, the Laplace approximation, importance sampling, and nested sampling. We found the dispersion in Z across different methods grew with increasing n-planet models: ~3 for 0-planets, ~10 for 1-planet, ~100-1000 for 2-planets, and >10,000 for 3-planets. Most internal estimates of uncertainty in Z for individual methods significantly underestimated the observed dispersion across all methods. Methods that adopted a Monte Carlo approach by comparing estimates from multiple runs yielded plausible uncertainties. Finally, two classes of numerical algorithms (those based on importance and nested samplers) arrived at similar conclusions regarding the ratio of Zs for n and (n+1)-planet models. One analytic method (the Laplace approximation) demonstrated comparable performance. We express both optimism and caution: we demonstrate that it is practical to perform rigorous Bayesian model comparison for <=3-planet models, yet robust planet discoveries require researchers to better understand the uncertainty in Z and its connections to model selection.

ALP-photon couplings are modeled in large ensembles of string vacua and random matrix theories. In all cases, the effective coupling increases polynomially in the number of ALPs, of which hundreds or thousands are expected in the string ensembles, many of which are ultralight. The expected value of the couplings $g_{a\gamma\gamma}\simeq 10^{-12}\text{GeV}^{-1} - 10^{-10}\text{GeV}^{-1}$ provide viable targets for future x-ray telescopes and axion helioscopes, and in some cases are already in tension with existing data.

The presence of mean motion resonances (MMRs) complicates analysis and fitting of planetary systems observed through the radial velocity (RV) technique. MMR can allow planets to remain stable in regions of phase space where strong planet-planet interactions would otherwise destabilize the system. These stable orbits can occupy small phase space volumes, allowing MMRs to strongly constrain system parameters, but making searches for stable orbital parameters challenging. Furthermore, libration of the resonant angle and dynamical interaction between the planets introduces another, long period variation into the observed RV signal, complicating analysis of the periods of the planets in the system. We discuss this phenomenon using the example of HD 200964. By searching through parameter space and numerically integrating each proposed set of planetary parameters to test for long term stability, we find stable solutions in the 7:5 and 3:2 MMRs in addition to the originally identified 4:3 MMR. The 7:5 configuration provides the best match to the data, while the 3:2 configuration provides the most easily understood formation scenario. In reanalysis of the originally published shorter-baseline data, we find fits in both the 4:3 and 3:2 resonances, but not the 7:5. Because the time baseline of the data is less than the resonant libration period, the current best fit to the data may not reflect the actual resonant configuration. In the absence of a full sample of the longer libration period, we find that it is of paramount importance to incorporate long term stability when fitting for the system's orbital configuration.

We infer the number of planets-per-star as a function of orbital period and planet size using $Kepler$ archival data products with updated stellar properties from the $Gaia$ Data Release 2. Using hierarchical Bayesian modeling and Hamiltonian Monte Carlo, we incorporate planet radius uncertainties into an inhomogeneous Poisson point process model. We demonstrate that this model captures the general features of the outcome of the planet formation process around GK stars, and provides an infrastructure to use the $Kepler$ results to constrain analytic planet distribution models. We report an increased mean and variance in the marginal posterior distributions for the number of planets per $GK$ star when including planet radius measurement uncertainties. We estimate the number of planets-per-$GK$ star between 0.75 and 2.5 $R_{\oplus}$ and 50 to 300 day orbital periods to have a $68\%$ credible interval of $0.49$ to $0.77$ and a posterior mean of $0.63$. This posterior has a smaller mean and a larger variance than the occurrence rate calculated in this work and in Burke et al. (2015) for the same parameter space using the $Q1-Q16$ (previous $Kepler$ planet candidate and stellar catalog), and a larger mean and variance than when using the $DR25$ (latest $Kepler$ planet candidate and stellar catalog). We find that the accuracy and precision of our hierarchical Bayesian model posterior distributions are less sensitive to the total number of planets in the sample, and more so on the characteristics of the catalog completeness and reliability and the span of the planet parameter space.

We demonstrate that asymmetric reheating arises in a large ensemble of string compactifications with many axions and gauged dark sectors. This phenomenon may help avoid numerous cosmological problems that may arise if the sectors were reheated democratically. Distributions of couplings are presented for two classes of axion reheatons, both of which exhibit very small couplings to most of the gauge sectors. In one class, ratios of reheating couplings and also preferred gauge groups are frequently determined by local regions in the string geometry.

We apply classical machine vision and machine deep learning methods to prototype signal classifiers for the search for extraterrestrial intelligence. Our novel approach uses two-dimensional spectrograms of measured and simulated radio signals bearing the imprint of a technological origin. The studies are performed using archived narrow-band signal data captured from real-time SETI observations with the Allen Telescope Array and a set of digitally simulated signals designed to mimic real observed signals. By treating the 2D spectrogram as an image, we show that high quality parametric and non-parametric classifiers based on automated visual analysis can achieve high levels of discrimination and accuracy, as well as low false-positive rates. The (real) archived data were subjected to numerous feature-extraction algorithms based on the vertical and horizontal image moments and Huff transforms to simulate feature rotation. The most successful algorithm used a two-step process where the image was first filtered with a rotation, scale and shift-invariant affine transform followed by a simple correlation with a previously defined set of labeled prototype examples. The real data often contained multiple signals and signal ghosts, so we performed our non-parametric evaluation using a simpler and more controlled dataset produced by simulation of complex-valued voltage data with properties similar to the observed prototypes. The most successful non-parametric classifier employed a wide residual (convolutional) neural network based on pre-existing classifiers in current use for object detection in ordinary photographs. These results are relevant to a wide variety of research domains that already employ spectrogram analysis from time-domain astronomy to observations of earthquakes to animal vocalization analysis.

We run a suite of dissipative N-body simulations to determine which regions of phase space for smooth disk migration are consistent with the GJ876 system, an M-dwarf hosting three planets orbiting in a chaotic 4:2:1 Laplace resonance. We adopt adaptive mesh refinement (AMR) methods which are commonly used in hydrodynamical simulations to efficiently explore the parameter space defined by the semi-major axis and eccentricity damping timescales. We find that there is a large region of phase space which produces systems in the chaotic Laplace resonance and a smaller region consistent with the observed eccentricities and libration amplitudes for the resonant angles. Under the assumptions of Type I migration for the outer planet, we translate these damping timescales into constraints on the protoplanetary disk surface density and thickness. When we strongly (weakly) damp the eccentricities of the inner two Laplace planets, these timescales correspond to disk surface densities around ten thousand (a few hundred) grams per square centimeter and disk aspect ratios between 1-10%. Additionally, smooth migration produces systems with a range of chaotic timescales, from decades and centuries to upwards of thousands of years. In agreement with previous studies, the less chaotic regions of phase space coincide with the system being in a low energy double apsidal corotation resonance. Our detailed modeling of multi-planetary systems coupled with our AMR exploration method enhances our ability to map out the parameter space of planet formation models, and is well suited to study other resonant chain systems such as Trappist-1, Kepler-60, and others.

There is a degeneracy in the radial velocity exoplanet signal between a single planet on an eccentric orbit and a two-planet system with a period ratio of 2:1. This degeneracy could lead to misunderstandings of the dynamical histories of planetary systems as well as measurements of planetary abundances if the correct architecture is not established. We constrain the rate of mischaracterization by analyzing a sample of 60 non-transiting, radial velocity systems orbiting main sequence stars from the NASA Exoplanet Archive (NASA Archive) using a new Bayesian model comparison pipeline. We find that 15 systems (25% of our sample) show compelling evidence for the two-planet case with a confidence level of 95%.

Self-organization properties of sustained magnetized plasma are applied to selected solar data to understand solar magnetic fields. Torsional oscillations are speed-up and slow-down bands of the azimuthal flow that correlate with the solar cycle, and they imply the existence of a symmetric solar dynamo with a measured polar flux of 3x10^14 Wb. It is shown that the solar dynamo is thin (~0.1 Mm gradient scale size) and powerful (~10^23 W). These properties are found from the amplitude of the torsional oscillations and the relationship of their velocity contours to solar magnetograms supports the result. The dynamo has enough power to heat the chromosphere and to power the corona and the solar wind. The dynamo also causes a rigid rotation of the heliosphere out to at least the corona and the relationship of the rotation of the corona to solar magnetograms supports this result as well. The thin solar dynamo sustains a thin stable minimum energy state that seems to be covering most of the solar surface just below the photosphere. The magnetic field lines of the minimum energy state should be parallel to the solar surface and rotate with distance from the surface with 2{\pi} radians of rotation in ~1 Mm Resistive diffusion helps to push the magnetic fields to the surface and the global magnetic structure (GMS) seems to lose {\pi} radians every 11 years, causing the observed 180 degree flipping of the solar magnetic field. The thin sheets of magnetized plasma in solar prominences may be the lost thin sheets of the GMS. For completeness, the formation of sunspots, CMEs and flares is discussed.

We report the results of a search for candidate Li-rich giants among 569,738 stars of the SDSS DR10 dataset. With small variations, our approach is based on that taken in an earlier search for EMP/CEMP stars and uses the same dataset. As part of our investigation, we demonstrate a method for separating post-main sequence and main sequence stars cooler than Teff ~ 5800 K using our feature strength measures of the Sr II 4078, Fe I 4072, and Ca I 4227 lines. By taking carefully selected cuts in a multi-dimensional phase space, we isolate a sample of potential Li-rich giant stars. From these, using detailed comparison with dwarf and giant MILES stars, and our own individual spectral classifications, we identify a set of high likelihood candidate Li-rich giant stars. We offer these for further study to promote an understanding of these enigmatic objects.

We outline the important role that ground-based, Doppler monitoring of exoplanetary systems will play in advancing our theories of planet formation and dynamical evolution. A census of planetary systems requires a well designed survey to be executed over the course of a decade or longer. A coordinated survey to monitor several thousand targets each at ~1000 epochs (~3-5 million new observations) will require roughly 40 dedicated spectrographs. We advocate for improvements in data management, data sharing, analysis techniques, and software testing, as well as possible changes to the funding structures for exoplanet science.

The solar dynamo and the solar Global internal Magnetic Structure (GMS) appear to be a thin ($\sim$2 Mm thick) structure near ($\sim$1 Mm below) the solar surface. Evidence for these properties are found from the amplitude of the torsional oscillations and in their velocity contours relationship to solar magnetogram; the power to the chromosphere; power to the corona and the solar wind; the current in the helio-current-sheet measured at the radius of the orbit of Earth; the calculated size ($\sim$1 Mm) of the expanding polar flux when it enters the photosphere; and from the observation that solar magnetic activity is generated near the surface. A thin stable minimum energy state seems to be covering most of the solar surface just below the photosphere. The magnetic field lines should be parallel to the solar surface and rotate with distance from the surface for 2$\pi$ radians in $\sim$2 Mm. Resistive diffusion helps to push the magnetic fields to the surface and the GMS seems to lose $\pi$ radians every 11 years, causing the observed 180$^\circ$ flipping of the solar magnetic fields including the flipping of the polar flux. Further evidence for this GMS and its loss is that solar prominences are made of thin sheets of magnetized plasma, which are, likely, remnants of the lost thin sheet of the GMS. The loss process is consistent with the butterfly pattern of the sunspots and with the differences observed between solar maximum and solar minimum in the corona. The solar dynamo drives current parallel to the polar flux, which, in turn, sustains the GMS using cross-field current drive. For completeness, the formation of sunspots, CMEs and flares is discussed.

We present the results of a search for extremely metal-poor (EMP), carbon-enhanced metal-poor (CEMP), and cataclysmic variable (CV) stars using a new exploration tool based on linked scatter plots (LSPs). Our approach is especially designed to work with very large spectrum data sets such as the SDSS, LAMOST, RAVE, and Gaia data sets, and it can be applied to stellar, galaxy, and quasar spectra. As a demonstration, we conduct our search using the SDSS DR10 data set. We first created a 3326-dimensional phase space containing nearly 2 billion measures of the strengths of over 1600 spectral features in 569,738 SDSS stars. These measures capture essentially all the stellar atomic and molecular species visible at the resolution of SDSS spectra. We show how LSPs can be used to quickly isolate and examine interesting portions of this phase space. To illustrate, we use LSPs coupled with cuts in selected portions of phase space to extract EMP stars, CEMP stars, and CV stars. We present identifications for 59 previously unrecognized candidate EMP stars and 11 previously unrecognized candidate CEMP stars. We also call attention to 2 candidate He~II emission CV stars found by the LSP approach that have not yet been discussed in the literature.

Disk migration and high-eccentricity migration are two well-studied theories to explain the formation of hot Jupiters. The former predicts that these planets can migrate up until the planet-star Roche separation ($a_{Roche}$) and the latter predicts they will tidally circularize at a minimum distance of 2$a_{Roche}$. Considering long-running radial velocity and transit surveys have identified a couple hundred hot Jupiters to date, we can revisit the classic question of hot Jupiter formation in a data-driven manner. We approach this problem using data from several exoplanet surveys (radial velocity, Kepler, HAT, and WASP) allowing for either a single population or a mixture of populations associated with these formation channels, and applying a hierarchical Bayesian mixture model of truncated power laws of the form $x^{\gamma-1}$ to constrain the population-level parameters of interest (e.g., location of inner edges, $\gamma$, mixture fractions). Within the limitations of our chosen models, we find the current radial velocity and Kepler sample of hot Jupiters can be well explained with a single truncated power law distribution with a lower cutoff near 2$a_{Roche}$, a result that still holds after a decade, and $\gamma=-0.51\pm^{0.19}_{0.20}$. However, the HAT and WASP data show evidence for multiple populations (Bayes factor $\approx 10^{21}$). We find that $15\pm^{9}_{6}\%$ reside in a component consistent with disk migration ($\gamma=-0.04\pm^{0.53}_{1.27}$) and $85\pm^{6}_{9}\%$ in one consistent with high-eccentricity migration ($\gamma=-1.38\pm^{0.32}_{0.47}$). We find no immediately strong connections with some observed host star properties and speculate on how future exoplanet surveys could improve upon hot Jupiter population inference.

We report 1212 radial-velocity (RV) measurements obtained in the years 2009-2013 using an iodine cell for the spectroscopic binary nu Octantis (K1III/IV). This system (a_bin~2.6 au, P~1050 days) is conjectured to have a Jovian planet with a semi-major axis half that of the binary host. The extreme geometry only permits long-term stability if the planet is in a retrograde orbit. Whilst the reality of the planet (P~415 days) remains uncertain, other scenarios (stellar variability or apsidal motion caused by a yet unobserved third star) continue to appear substantially less credible based on CCF bisectors, line-depth ratios and many other independent details. If this evidence is validated but the planet is disproved, the claims of other planets using RVs will be seriously challenged. We also describe a significant revision to the previously published RVs and the full set of 1437 RVs now encompasses nearly 13 years. The sensitive orbital dynamics allow us to constrain the three-dimensional architecture with a broad prior probability distribution on the mutual inclination, which with posterior samples obtained from an N-body Markov chain Monte Carlo is found to be 158.4 +/- 1.2 deg. None of these samples are dynamically stable beyond 1 Myr. However, a grid search around the best-fitting solution finds a region that has many models stable for 10 Myr, and includes one model within 1-sigma that is stable for at least 100 Myr. The planet's exceptional nature demands robust independent verification and makes the theoretical understanding of its formation a worthy challenge.

We report constraints on the three-dimensional orbital architecture for all four planets known to orbit the nearby M dwarf Gliese 876 based solely on Doppler measurements and demanding long-term orbital stability. Our dataset incorporates publicly available radial velocities taken with the ELODIE and CORALIE spectrographs, HARPS, and Keck HIRES as well as previously unpublished HIRES velocities. We first quantitatively assess the validity of the planets thought to orbit GJ 876 by computing the Bayes factors for a variety of different coplanar models using an importance sampling algorithm. We find that a four-planet model is preferred over a three-planet model. Next, we apply a Newtonian MCMC algorithm to perform a Bayesian analysis of the planet masses and orbits using an n-body model in three-dimensional space. Based on the radial velocities alone, we find that a 99% credible interval provides upper limits on the mutual inclinations for the three resonant planets ($\Phi_{cb}<6.20^\circ$ for the "c" and "b" pair and $\Phi_{be}<28.5^\circ$ for the "b" and "e" pair). Subsequent dynamical integrations of our posterior sample find that the GJ 876 planets must be roughly coplanar ($\Phi_{cb}<2.60^\circ$ and $\Phi_{be}<7.87^\circ$), suggesting the amount of planet-planet scattering in the system has been low. We investigate the distribution of the respective resonant arguments of each planet pair and find that at least one argument for each planet pair and the Laplace argument librate. The libration amplitudes in our three-dimensional orbital model supports the idea of the outer-three planets having undergone significant past disk migration.

We perform a comprehensive study of models of dark matter (DM) in a Universe with a non-thermal cosmological history, i.e with a phase of pressure-less matter domination before the onset of big-bang nucleosynethesis (BBN). Such cosmological histories are generically predicted by UV completions that contain gravitationally coupled scalar fields (moduli). We classify the different production mechanisms for DM in this framework, generalizing previous works by considering a wide range of DM masses/couplings and allowing for DM to be in equilibrium with a "dark" sector. We identify four distinct parametric regimes for the production of relic DM, and derive accurate semi-analytic approximations for the DM relic abundance. Our results are particularly relevant for supersymmetric theories, in which the standard non-thermally produced DM candidates are disfavored by indirect detection constraints. We also comment on experimental signals in this framework, focusing on novel effects involving the power spectrum of DM density perturbations. In particular, we identify a class of models where the spectrum of DM density perturbations is sensitive to the pressure-less matter dominated era before BBN, giving rise to interesting astrophysical signatures to be looked for in the future. A worthwhile future direction would be to study well-motivated theoretical models within this framework and carry out detailed studies of the pattern of expected experimental signals.

We present an update to seven stars with long-period planets or planetary candidates using new and archival radial velocities from Keck-HIRES and literature velocities from other telescopes. Our updated analysis better constrains orbital parameters for these planets, four of which are known multi-planet systems. HD 24040 b and HD 183263 c are super-Jupiters with circular orbits and periods longer than 8 yr. We present a previously unseen linear trend in the residuals of HD 66428 indicative on an additional planetary companion. We confirm that GJ 849 is a multi-planet system and find a good orbital solution for the c component: it is a $1 M_{\rm Jup}$ planet in a 15 yr orbit (the longest known for a planet orbiting an M dwarf). We update the HD 74156 double-planet system. We also announce the detection of HD 145934 b, a $2 M_{\rm Jup}$ planet in a 7.5 yr orbit around a giant star. Two of our stars, HD 187123 and HD 217107, at present host the only known examples of systems comprising a hot Jupiter and a planet with a well constrained period $> 5$ yr, and with no evidence of giant planets in between. Our enlargement and improvement of long-period planet parameters will aid future analysis of origins, diversity, and evolution of planetary systems.

We present an updated study of the planets known to orbit 55 Cancri A using 1,418 high-precision radial velocity observations from four observatories (Lick, Keck, Hobby-Eberly Telescope, Harlan J. Smith Telescope) and transit time/durations for the inner-most planet, 55 Cancri "e" (Winn et al. 2011). We provide the first posterior sample for the masses and orbital parameters based on self-consistent n-body orbital solutions for the 55 Cancri planets, all of which are dynamically stable (for at least $10^8$ years). We apply a GPU version of Radial velocity Using N-body Differential evolution Markov Chain Monte Carlo (RUN DMC; B. Nelson et al. 2014) to perform a Bayesian analysis of the radial velocity and transit observations. Each of the planets in this remarkable system has unique characteristics. Our investigation of high-cadence radial velocities and priors based on space-based photometry yields an updated mass estimate for planet "e" ($8.09\pm0.26$ M$_\oplus$), which affects its density ($5.51\pm^{1.32}_{1.00}$ g cm$^{-3}$) and inferred bulk composition. Dynamical stability dictates that the orbital plane of planet "e" must be aligned to within $60^o$ of the orbital plane of the outer planets (which we assume to be coplanar). The mutual interactions between the planets "b" and "c" may develop an apsidal lock about $180^o$. We find 36-45% of all our model systems librate about the anti-aligned configuration with an amplitude of $51^o\pm^{6^o}_{10^o}$. Other cases showed short-term perturbations in the libration of $\varpi_b-\varpi_c$, circulation, and nodding, but we find the planets are not in a 3:1 mean-motion resonance. A revised orbital period and eccentricity for planet "d" pushes it further toward the closest known Jupiter analog in the exoplanet population.

In the 20+ years of Doppler observations of stars, scientists have uncovered a diverse population of extrasolar multi-planet systems. A common technique for characterizing the orbital elements of these planets is Markov chain Monte Carlo (MCMC), using a Keplerian model with random walk proposals and paired with the Metropolis-Hastings algorithm. For approximately a couple of dozen planetary systems with Doppler observations, there are strong planet-planet interactions due to the system being in or near a mean-motion resonance (MMR). An N-body model is often required to accurately describe these systems. Further computational difficulties arise from exploring a high-dimensional parameter space ($\sim$7 x number of planets) that can have complex parameter correlations. To surmount these challenges, we introduce a differential evolution MCMC (DEMCMC) applied to radial velocity data while incorporating self-consistent N-body integrations. Our Radial velocity Using N-body DEMCMC (RUN DMC) algorithm improves upon the random walk proposal distribution of the traditional MCMC by using an ensemble of Markov chains to adaptively improve the proposal distribution. We describe the methodology behind the algorithm, along with results of tests for accuracy and performance. We find that most algorithm parameters have a modest effect on the rate of convergence. However, the size of the ensemble can have a strong effect on performance. We show that the optimal choice depends on the number of planets in a system, as well as the computer architecture used and the resulting extent of parallelization. While the exact choices of optimal algorithm parameters will inevitably vary due to the details of individual planetary systems, we offer recommendations for choosing the DEMCMC algorithm's algorithmic parameters that result in excellent performance for a wide variety of planetary systems.

We present Swarm-NG, a C++ library for the efficient direct integration of many n-body systems using highly-parallel Graphics Processing Unit (GPU), such as NVIDIA's Tesla T10 and M2070 GPUs. While previous studies have demonstrated the benefit of GPUs for n-body simulations with thousands to millions of bodies, Swarm-NG focuses on many few-body systems, e.g., thousands of systems with 3...15 bodies each, as is typical for the study of planetary systems. Swarm-NG parallelizes the simulation, including both the numerical integration of the equations of motion and the evaluation of forces using NVIDIA's "Compute Unified Device Architecture" (CUDA) on the GPU. Swarm-NG includes optimized implementations of 4th order time-symmetrized Hermite integration and mixed variable symplectic integration, as well as several sample codes for other algorithms to illustrate how non-CUDA-savvy users may themselves introduce customized integrators into the Swarm-NG framework. To optimize performance, we analyze the effect of GPU-specific parameters on performance under double precision. Applications of Swarm-NG include studying the late stages of planet formation, testing the stability of planetary systems and evaluating the goodness-of-fit between many planetary system models and observations of extrasolar planet host stars (e.g., radial velocity, astrometry, transit timing). While Swarm-NG focuses on the parallel integration of many planetary systems,the underlying integrators could be applied to a wide variety of problems that require repeatedly integrating a set of ordinary differential equations many times using different initial conditions and/or parameter values.

We report the discovery via radial velocity of a short-period (P = 2.430420 \pm 0.000006 days) companion to the F-type main sequence star TYC 2930-00872-1. A long-term trend in the radial velocities indicates the presence of a tertiary stellar companion with $P > 2000$ days. High-resolution spectroscopy of the host star yields T_eff = 6427 +/- 33 K, log(g) = 4.52 +/- 0.14, and [Fe/H]=-0.04 +/- 0.05. These parameters, combined with the broad-band spectral energy distribution and parallax, allow us to infer a mass and radius of the host star of M_1=1.21 +/- 0.08 M_\odot and R_1=1.09_{-0.13}^{+0.15} R_\odot. We are able to exclude transits of the inner companion with high confidence. The host star's spectrum exhibits clear Ca H and K core emission indicating stellar activity, but a lack of photometric variability and small v*sin(I) suggest the primary's spin axis is oriented in a pole-on configuration. The rotational period of the primary from an activity-rotation relation matches the orbital period of the inner companion to within 1.5 \sigma, suggesting they are tidally locked. If the inner companion's orbital angular momentum vector is aligned with the stellar spin axis, as expected through tidal evolution, then it has a stellar mass of M_2 ~ 0.3-0.4 M_\odot. Direct imaging limits the existence of stellar companions to projected separations < 30 AU. No set of spectral lines and no significant flux contribution to the spectral energy distribution from either companion are detected, which places individual upper mass limits of M < 1.0 M_\odot, provided they are not stellar remnants. If the tertiary is not a stellar remnant, then it likely has a mass of ~0.5-0.6 M_\odot, and its orbit is likely significantly inclined from that of the secondary, suggesting that the Kozai-Lidov mechanism may have driven the dynamical evolution of this system.

We present a focused study of a predictive unified model whose measurable consequences are immediately relevant to early discovery prospects of supersymmetry at the LHC. ATLAS and CMS have released their analysis with 35~pb$^{-1}$ of data and the model class we discuss is consistent with this data. It is shown that with an increase in luminosity the LSP dark matter mass and the gluino mass can be inferred from simple observables such as kinematic edges in leptonic channels and peak values in effective mass distributions. Specifically, we consider cases in which the neutralino is of low mass and where the relic density consistent with WMAP observations arises via the exchange of Higgs bosons in unified supergravity models. The magnitudes of the gaugino masses are sharply limited to focused regions of the parameter space, and in particular the dark matter mass lies in the range $\sim (50-65) ~\rm GeV$ with an upper bound on the gluino mass of $575~{\rm GeV}$, with a typical mass of $450~{\rm GeV}$. We find that all model points in this paradigm are discoverable at the LHC at $\sqrt s = 7 \rm ~TeV$. We determine lower bounds on the entire sparticle spectrum in this model based on existing experimental constraints. In addition, we find the spin-independent cross section for neutralino scattering on nucleons to be generally in the range of $\sigma^{\rm SI}_{\na p} = 10^{-46 \pm 1}~\rm cm^2$ with much higher cross sections also possible. Thus direct detection experiments such as CDMS and XENON already constrain some of the allowed parameter space of the low mass gaugino models and further data will provide important cross-checks of the model assumptions in the near future.

We present ground-based optical observations of the September 2009 and January 2010 transits of HD 80606b. Based on 3 partial light curves of the September 2009 event, we derive a midtransit time of T_c [HJD] = 2455099.196 +- 0.026, which is about 1 sigma away from the previously predicted time. We observed the January 2010 event from 9 different locations, with most phases of the transit being observed by at least 3 different teams. We determine a midtransit time of T_c [HJD] = 2455210.6502 +- 0.0064, which is within 1.3 sigma of the time derived from a Spitzer observation of the same event.

The Two Micron All-Sky Survey (2MASS) has provided a uniform photometric catalog to search for previously unknown red AGN and QSOs. We have extended the search to the southern equatorial sky by obtaining spectra for 1182 AGN candidates using the 6dF multifibre spectrograph on the UK Schmidt Telescope. These were scheduled as auxiliary targets for the 6dF Galaxy Redshift Survey. The candidates were selected using a single color cut of J - Ks > 2 to Ks ~ 15.5 and a galactic latitude of |b|>30 deg. 432 spectra were of sufficient quality to enable a reliable classification. 116 sources (or ~27%) were securely classified as type 1 AGN, 20 as probable type 1s, and 57 as probable type 2 AGN. Most of them span the redshift range 0.05<z<0.5 and only 8 (or ~6%) were previously identified as AGN or QSOs. Our selection leads to a significantly higher AGN identification rate amongst local galaxies (>20%) than in any previous galaxy survey. A small fraction of the type 1 AGN could have their optical colors reddened by optically thin dust with A_V<2 mag relative to optically selected QSOs. A handful show evidence for excess far-IR emission. The equivalent width (EW) and color distributions of the type 1 and 2 AGN are consistent with AGN unified models. In particular, the EW of the [OIII] emission line weakly correlates with optical--near-IR color in each class of AGN, suggesting anisotropic obscuration of the AGN continuum. Overall, the optical properties of the 2MASS red AGN are not dramatically different from those of optically-selected QSOs. Our near-IR selection appears to detect the most near-IR luminous QSOs in the local universe to z~0.6 and provides incentive to extend the search to deeper near-IR surveys.

The PAMELA positron excess is analyzed within the framework of nonuniversal SUGRA models with an extended $U(1)^n$ gauge symmetry in the hidden sector leading to neutralino dark matter with either a mixed Higgsino-wino LSP or an essentially pure wino dominated LSP. The Higgsino-wino LSP can produce the observed PAMELA positron excess and satisfy relic density constraints in the extended class of models due to a near degeneracy of the mass spectrum of the extended neutralino sector with the LSP mass. The simultaneous satisfaction of the WMAP relic density data and the PAMELA data is accomplished through a co-annihilation mechanism ($B_{\rm Co}-mechanism$), and leads to predictions of a neutralino and a chargino in the mass range (180-200) GeV as well as low lying sparticles accessible at colliders. We show that the models are consistent with the antiproton constraints from PAMELA as well as photon flux data from EGRET and FERMI-LAT. Predictions for the scalar neutralino proton cross section relevant for the direct detection of dark matter are also discussed and signatures at the LHC for these PAMELA inspired models are analyzed. It is shown that the mixed Higgsino-wino LSP model will be discoverable with as little as 1 fb$^{-1}$ of data and is thus a prime candidate for discovery in the low luminosity runs at the LHC.

We present infrared (IR) to X-ray spectral energy distributions (SEDs) for 44 red AGN selected from the 2MASS survey on the basis of their red J-K$_S$ color (>2 mag) and later observed by Chandra. In comparison with optically-, radio-, and X-ray selected AGN, their median SEDs are red in the optical and near-IR with little/no blue bump. It thus seems that near-IR color selection isolates the reddest subset of AGN that can be classified optically. The shape of the SEDs is generally consistent with modest absorption by gas (in the X-ray) and dust (in the optical-IR). The levels of obscuration, estimated from X-rays, far-IR and our detailed optical/near-IR color modeling are all consistent implying N_H < few*10^{22} cm^{-2}. We present SED models that show how the AGN optical/near-IR colors change due to differing amounts of reddening, AGN to host galaxy ratio, redshift and scattered light emission and apply them to the sources in the sample. We find that the 2MASS AGN optical color, B-R, and to a lesser extent the near-IR color, J-K$_S$, are strongly affected by reddening, host galaxy emission, redshift, and in few, highly polarized objects, also by scattered AGN light. The obscuration/inclination of the AGN allows us to see weaker emission components which are generally swamped by the AGN.

We have obtained optical imaging with the Canada France Hawaii Telescope (CFHT) of 21 2MASS-selected QSOs of redshift greater than 0.3. This paper complements the sample of lower redshift 2MASS QSOs previously published. The QSOs have higher overall and nuclear luminosity, bluer colours, and higher ratio of nuclear to host flux than the lower redshift sample. From these and other properties, we argue that the sample is consistent with the emergence of the AGN from dusty starbursts following major tidal interactions between galaxies.

XMM-Newton spectra of five red, 2MASS AGN, selected from a sample observed by Chandra to be relatively X-ray bright and to cover a range of hardness ratios, confirm the presence of substantial absorbing material in three sources with optical classifications ranging from Type 1 to Type 2. A flat (hard), power law continuum is observed in the other two. The combination of X-ray absorption and broad optical emission lines suggests either a small (nuclear) absorber or a favored viewing angle so as to cover the X-ray source but not the broad emission line region (BELR). A soft excess is detected in all three Type 1 sources. We speculate that this may arise in an extended region of ionised gas, perhaps linked with the polarised (scattered) optical light present in these sources. The spectral complexity revealed by XMM-Newton emphasizes the limitations of the low S/N \chandra data. The new results strengthen our earlier conclusions that the observed X-ray continua of red AGN are unusually hard at energies >2 keV. Their observed spectra are consistent with contributing significantly to the missing hard/absorbed population of the Cosmic X-ray Background (CXRB) although their intrinsic power law slopes are typical of broad-line (Type 1) AGN (Gamma ~1.7-1.9). This suggests that the missing X-ray-absorbed CXRB population may include Type 1 AGN/QSOs in addition to the Type 2 AGN generally assumed.

(Abridged) We present an unbiassed near-IR selected AGN sample, covering 12.56 square degrees down to K ~ 15.5, selected from the Two Micron All Sky Survey (2MASS). Our only selection effect is a moderate color cut (J-K>1.2) designed to reduce contamination from galactic stars. We observed both point-like and extended sources. Using the brute-force capabilities of the 2dF multi-fiber spectrograph on the Anglo-Australian Telescope, we obtained spectra of 65% of the target list: an unbiassed sub-sample of 1526 sources. 80% of the 2MASS sources in our fields are galaxies, with a median redshift of 0.15. The remainder are K- and M-dwarf stars. Seyfert-2 Galaxies are roughly three times more common in this sample than in optically selected galaxy samples (once corrections have been made for the equivalent width limit and for different aperture sizes). We find 14 broad-line (Type-1) AGNs, giving a surface density down to K<15 comparable to that of optical samples down to B=18.5. Half of our Type-1 AGNs could not have been found by normal color selection techniques. In all cases this was due host galaxy light contamination rather than intrinsically red colors. We conclude that the Type-1 AGN population found in the near-IR is not dramatically different from that found in optical samples. There is no evidence for a large population of AGNs that could not be found at optical wavelengths, though we can only place very weak constraints on any population of dusty high-redshift QSOs.

We present and discuss optical imaging of 76 AGN which represent the 2MASS-selected sample for z<0.3, from a full list of 243. They are found to have dust-obscured nuclei, residing in host galaxies that show a high fraction (>70%) of tidal interactions. The derived luminosities of the AGN and host galaxies are similar to traditionally-selected AGN, and they may comprise some 40% of the total AGN population at low redshift. We have measured a number of host galaxy properties, and discuss their distributions and interrelations. We compare the 2MASS AGN with optically selected samples and the IRAS-selected galaxy samples, and discuss the differences in terms of merger processes and initial conditions.

Using simple infrared color selection, 2MASS has found a large number of red, previously unidentified, radio-quiet QSOs. Although missed by UV/optical surveys, the 2MASS QSOs have K-band luminosities that are comparable to "classical" QSOs. This suggests the possible discovery of a previously predicted large population of dust-obscured radio-quiet QSOs. We present the results of an imaging survey of 29 2MASS QSOs observed with WFPC2 onboard the Hubble Space Telescope. I-band images, which benefit from the relative faintness of the nuclei at optical wavelengths, are used to characterize the host galaxies, measure the nuclear contribution to the total observed I-band emission, and to survey the surrounding environments. The 2MASS QSOs are found to lie in galaxies with a variety of morphologies, luminosities, and dynamical states, not unlike those hosting radio-quiet PG QSOs. Our analysis suggests that the extraordinary red colors of the 2MASS QSOs are caused by extinction of an otherwise typical QSO spectrum due to dust near the nucleus.

We extend previous analyses on the interplay between nonuniversalities in the gaugino mass sector and the thermal relic densities of LSP neutralinos, in particular to the case of moderate to large tan beta. We introduce a set of parameters that generalizes the standard unified scenario to cover the complete allowed parameter space in the gaugino mass sector. We discuss the physical significance of the cosmologically preferred degree of degeneracy between charginos and the LSP and study the effect this degree of degeneracy has on the prospects for direct detection of relic neutralinos in the next round of dark matter detection experiments. Lastly, we compare the fine tuning required to achieve a satisfactory relic density with the case of universal gaugino masses, as in minimal supergravity, and find it to be of a similar magnitude. The sensitivity of quantifiable measures of fine-tuning on such factors as the gluino mass and top and bottom masses is also examined.

The Two Micron All Sky Survey (2MASS) is finding previously unidentified, luminous red active galactic nuclei (AGN). This new sample has a space density similar to, or greater than, previously known AGN, suggesting that a large fraction of the overall population has been missed. Chandra observations of a well-defined subset of these objects reveal that all are X-ray faint, with the reddest sources being the faintest in X-rays. The X-ray hardness ratios cover a wide range, generally indicating log NH ~21-23 /cm*2, but the softest sources show no spectral evidence for intrinsic absorption. These characteristics suggest that a mix of absorbed, direct emission and unabsorbed, scattered and/or extended emission contributes to the X-ray flux, although we cannot rule out the possibility that they are intrinsically X-ray weak. This population of X-ray faint, predominantly broad-line objects could provide the missing population of X-ray absorbed AGN required by current models of the cosmic X-ray background. The existence of AGN which display both broad emission lines and absorbed X-rays has important implications for unification schemes and emphasizes the need for care in assigning classifications to individual AGN.

Optical broad-band polarimetry is presented for near-infrared color-selected AGN classified QSOs based on their K-band luminosity. More than 10% of a sample of 70 QSOs discovered in the Two Micron All Sky Survey (2MASS) with J-K > 2 and M_K < -23 show high broad-band linear polarization (P > 3%), and values range to a maximum of P = 11%. High polarization tends to be associated with the most luminous objects at K, and with QSOs having the highest near-IR-to-optical flux ratios. The 2MASS QSO sample includes objects possessing a wide range of optical spectral types. High polarization is seen in two of 22 broad emission-line (Type 1) objects, but 1/4 of the QSOs of intermediate spectral type (Type 1.5-1.9) are highly polarized. None of the nine QSOs classified as Type 2 exhibit P > 3%. It is likely that the unavoidable inclusion of unpolarized starlight from the host galaxy within the observation aperture results in reduced polarization for the narrow emission-line objects. The high polarization of 2MASS-discovered QSOs supports the conclusion inferred from their near-IR and optical colors, that the nuclei of many of these objects are obscured to some degree by dust. The 2MASS QSO sample is compared to other, predominantly radio-quiet, QSOs and is found to be consistent with the idea that the orientation of AGN to the line of sight plays a major role in determining their observed properties.

We investigate the dark matter prospects of supersymmetric models with nonuniversal gaugino masses. We find that for very particular values of the ratio of soft supersymmetry-breaking gaugino masses, M2/M1, an enhanced coannihilation efficiency between the lightest chargino and the lightest neutralino occurs, allowing for scalars with masses well above the normally accepted limits for viable dark matter in the universal case. As a specific example, we investigate models of hidden sector gaugino condensation. These models exhibit high scalar masses, previously thought dangerous, and the requisite freedom in the ratio of gaugino masses. The cosmologically viable regions of parameter space are investigated, allowing very specific statements to be made about the content of the supersymmetry-breaking hidden sector.

We present optical observations of two remarkable new AGN discovered by the Two-Micron All Sky Survey. Both are classified as QSOs based on their optical spectra, near-IR colors, and near-IR and [O III]5007 luminosities, but their optical polarizations are among the highest seen for non-blazar AGN; approaching 15% for 2MASSI J151653.2+190048. The polarized light spectrum for each object exhibits broad Balmer emission lines, but lacks the narrow lines that are evident in the total light spectrum. This is most pronounced for the Type-1.5 object 2MASSI J165939.7+183436, where broad lines dominate only in polarized light. The polarization properties of these AGN suggest that dust near the nucleus at least partially obscures the AGN and that material probably intermixed with the narrow line-emitting gas scatters nuclear light into our line of sight. These QSOs illustrate the variety of highly polarized AGN that have been missed by traditional optical search techniques, and demonstrate that such objects are exposed by surveys in the near-IR.

We present JHKs photometry, far red spectra, and spectral classifications for an additional 67 L dwarfs discovered by the Two Micron All Sky Survey. One of the goals of this new search was to locate more examples of the latest L dwarfs. Of the 67 new discoveries, 17 have types of L6 or later. Analysis of these new discoveries shows that H-alpha emission has yet to be convincingly detected in any L dwarf later than type L4.5, indicating a decline or absence of chromospheric activity in the latest L dwarfs. Further analysis shows that 16 (and possibly 4 more) of the new L dwarfs are lithium brown dwarfs and that the average line strength for those L dwarfs showing lithium increases until roughly type L6.5 V then declines for later types. This disappearance may be the first sign of depletion of atomic lithium as it begins to form into lithium-bearing molecules. Another goal of the search was to locate nearer, brighter L dwarfs of all subtypes. Using absolute magnitudes for 17 L dwarf systems with trigonometric parallax measurements, we develop spectrophotometric relations to estimate distances to the other L dwarfs. Of the 67 new discoveries, 21 have photometric distances placing them within 25 parsecs of the Sun. A table of all known L and T dwarfs believed to lie within 25 parsecs - 53 in total - is also presented. Using the distance measurement of the coolest L dwarf known, we calculate that the gap in temperature between L8 and the warmest known T dwarfs is less than 350K and probably much less. If the transition region between the two classes spans a very small temperature interval, this would explain why no transition objects have yet been uncovered. This evidence, combined with model fits to low-resolution spectra of late-M and early-L dwarfs, indicates that L dwarfs span the range 1300K < Teff < 2000K.

We present 263 optical and near-infrared (NIR) images for 42 Seyfert 1s and 48 Seyfert 2s, selected from the Extended 12-micron Galaxy Sample. Elliptically-averaged profiles are derived from the images, and isophotal radii and magnitudes are calculated from these. We also report virtual aperture photometry, that judging from comparison with previous work, is accurate to roughly 0.05mag in the optical, and 0.07mag in the NIR. Our B-band isophotal magnitude and radii, obtained from ellipse fitting, are in good agreement with those of RC3. When compared with the B band, V, I, J, and K isophotal diameters show that the colors in the outer regions of Seyferts are consistent with the colors of normal spirals. Differences in the integrated isophotal colors and comparison with a simple model show that the active nucleus+bulge is stronger and redder in the NIR than in the optical. Finally, roughly estimated Seyfert disk surface brightnesses are significantly brighter in B and K than those in normal spirals of similar morphological type.

We present the results of three years of ground-based observations of the Seyfert 1 galaxy NGC 5548, which combined with previously reported data, yield optical continuum and broad-line H-beta light curves for a total of eight years. The light curves consist of over 800 points, with a typical spacing of a few days between observations. During this eight-year period, the nuclear continuum has varied by more than a factor of seven, and the H-beta emission line has varied by a factor of nearly six. The H-beta emission line responds to continuum variations with a time delay or lag of 10-20 days, the precise value varying somewhat from year to year. We find some indications that the lag varies with continuum flux in the sense that the lag is larger when the source is brighter.