We present a new power spectrum emulator named EuclidEmulator that estimates the nonlinear correction to the linear dark matter power spectrum. It is based on a spectral decomposition method called polynomial chaos expansion. All steps in the construction of the emulator have been tested and optimized: the large high-resolution N-body simulations carried out with PKDGRAV3 were validated using a simulation from the Euclid Flagship campaign and demonstrated to have converged up to wavenumbers $k\approx 5\,h\,{\rm Mpc}^{-1}$ for redshifts $z\leq 5$. The emulator is constructed using the uncertainty quantification software UQLab and it has been optimized first by creating mock emulators based on Takahashi's HALOFIT. We show that it is possible to successfully predict the performance of the final emulator in this way prior to performing any N-body simulations. We provide a C-code to calculate the nonlinear correction at a relative accuracy of $\sim0.3\%$ with respect to N-body simulations within 50 ms. The absolute accuracy of the final nonlinear power spectrum is comparable to one obtained with N-body simulations, i.e. $\sim 1\%$ for $k\lesssim 1\,h\,{\rm Mpc}^{-1}$ and $z\lesssim 3.5$. This enables efficient forward modeling in the nonlinear regime allowing for maximum likelihood estimation of cosmological parameters. EuclidEmulator has been compared to HALOFIT and CosmicEmu, an alternative emulator based on the Mira-Titan Universe, and shown to be more accurate than these other approaches. This work paves a new way for optimal construction of future emulators that also consider other cosmological observables, use higher resolution input simulations and investigate higher dimensional cosmological parameter spaces.

Theoretical models for the expected merger rates of intermediate-mass black holes (IMBHs) are vital for planned gravitational-wave detection experiments such as the Laser Interferometer Space Antenna (LISA). Using collisionless $N$-body simulations of dwarf galaxy (DG) mergers, we examine how the orbital decay of IMBHs and the efficiency of IMBH binary formation depend on the central dark matter (DM) density profile of the merging DGs. Specifically, we explore various asymptotic inner slopes $\gamma$ of the DG's DM density distribution, ranging from steep cusps ($\gamma=1$) to shallower density profiles ($\gamma<1$), motivated by well-known baryonic-feedback effects as well as by DM models that differ from cold DM at the scales of DGs. We find that the inner DM slope is crucial for the formation (or lack thereof) of an IMBH binary; only mergers between DGs with cuspy DM profiles ($\gamma=1$) are favourable to forming a hard IMBH binary, whereas when $\gamma<1$ the IMBHs stall at a separation of 50-100 pc. Consequently, the rate of LISA signals from IMBH coalescence will be determined by the fraction of DGs with a cuspy DM profile. Conversely, the LISA event rates at IMBH mass scales offer in principle a novel way to place constraints on the inner structure of DM halos in DGs and address the core-cusp controversy. We also show that, with spatial resolutions of $\sim$0.1 kpc, as often adopted in cosmological simulations, all IMBHs stall, independent of $\gamma$. This suggests caution in employing cosmological simulations of galaxy formation to study BH dynamics in DGs.

We report on the successful completion of a 2 trillion particle cosmological simulation to z=0 run on the Piz Daint supercomputer (CSCS, Switzerland), using 4000+ GPU nodes for a little less than 80h of wall-clock time or 350,000 node hours. Using multiple benchmarks and performance measurements on the US Oak Ridge National Laboratory Titan supercomputer, we demonstrate that our code PKDGRAV3, delivers, to our knowledge, the fastest time-to-solution for large-scale cosmological N-body simulations. This was made possible by using the Fast Multipole Method in conjunction with individual and adaptive particle time steps, both deployed efficiently (and for the first time) on supercomputers with GPU-accelerated nodes. The very low memory footprint of PKDGRAV3 allowed us to run the first ever benchmark with 8 trillion particles on Titan, and to achieve perfect scaling up to 18000 nodes and a peak performance of 10 Pflops.

We present initial results from the "Ponos" zoom-in numerical simulations of dark matter substructures in massive ellipticals. Two very highly resolved dark matter halos with $M_{\rm vir}=1.2\times 10^{13}$ $M_{\odot}$ and $M_{\rm vir}=6.5\times 10^{12}$ $M_{\odot}$ and different ("violent" vs. "quiescent") assembly histories have been simulated down to $z=0$ in a $\Lambda$CDM cosmology with a total of 921,651,914 and 408,377,544 particles, respectively. Within the virial radius, the total mass fraction in self-bound $M_{\rm sub}>10^6$ $M_{\odot}$ subhalos at the present epoch is 15% for the violent host and 16.5% for the quiescent one. At $z=0.7$, these fractions increase to 19 and 33%, respectively, as more recently accreted satellites are less prone to tidal destruction. In projection, the average fraction of surface mass density in substructure at a distance of $R/R_{\rm vir}=0.02$ ($\sim 5-10$ kpc) from the two halo centers ranges from 0.6% to $\gtrsim 2$%, significantly higher than measured in simulations of Milky Way-sized halos. The contribution of subhalos with $M_{\rm sub} < 10^9$ $M_{\odot}$ to the projected mass fraction is between one fifth and one third of the total, with the smallest share found in the quiescent host. We assess the impact of baryonic effects via twin, lower-resolution hydrodynamical simulations that include metallicity-dependent gas cooling, star formation, and a delayed-radiative-cooling scheme for supernova feedback. Baryonic contraction produces a super-isothermal total density profile and increases the number of massive subhalos in the inner regions of the main host. The host density profiles and projected subhalo mass fractions appear to be broadly consistent with observations of gravitational lenses.

Future galaxy surveys require one percent precision in the theoretical knowledge of the power spectrum over a large range including very nonlinear scales. While this level of accuracy is easily obtained in the linear regime with perturbation theory, it represents a serious challenge for small scales where numerical simulations are required. In this paper we quantify the precision of present-day $N$-body methods, identifying main potential error sources from the set-up of initial conditions to the measurement of the final power spectrum. We directly compare three widely used $N$-body codes, Ramses, Pkdgrav3, and Gadget3 which represent three main discretisation techniques: the particle-mesh method, the tree method, and a hybrid combination of the two. For standard run parameters, the codes agree to within one percent at $k\leq1$ $h\,\rm Mpc^{-1}$ and to within three percent at $k\leq10$ $h\,\rm Mpc^{-1}$. We also consider the bispectrum and show that the reduced bispectra agree at the sub-percent level for $k\leq 2$ $h\,\rm Mpc^{-1}$. In a second step, we quantify potential errors due to initial conditions, box size, and resolution using an extended suite of simulations performed with our fastest code Pkdgrav3. We demonstrate that the simulation box size should not be smaller than $L=0.5$ $h^{-1}\rm Gpc$ to avoid systematic finite-volume effects (while much larger boxes are required to beat down the statistical sample variance). Furthermore, a maximum particle mass of $M_{\rm p}=10^{9}$ $h^{-1}\rm M_{\odot}$ is required to conservatively obtain one percent precision of the matter power spectrum. As a consequence, numerical simulations covering large survey volumes of upcoming missions such as DES, LSST, and Euclid will need more than a trillion particles to reproduce clustering properties at the targeted accuracy.

We explore fundamental properties of the distribution of low mass dark matter halos within the cosmic web using warm dark matter (WDM) and cold dark matter (CDM) cosmological simulations. Using self abundance-matched mock galaxy catalogs, we show that the distribution of dwarf galaxies in a WDM universe, wherein low mass halo formation is heavily suppressed, is nearly indistinguishable to that of a CDM universe whose low mass halos are not seen because galaxy formation is suppressed below some threshold halo mass. However, if the scatter between dwarf galaxy luminosity and halo properties is large enough, low mass CDM halos would sometimes host relatively bright galaxies thereby populating CDM voids with the occasional isolated galaxy and reducing the numbers of completely empty voids. Otherwise, without high mass to light scatter, all mock galaxy clustering statistics that we consider--the auto-correlation function, the numbers and radial profiles of satellites, the numbers of isolated galaxies, and the PDF of small voids--are nearly identical in CDM and WDM. WDM voids are neither larger nor emptier than CDM voids, when constructed from abundance-matched halo catalogs. It is thus a challenge to determine whether the CDM problem of the over-abundance of small halos with respect to the number density of observed dwarf galaxies has a cosmological solution or an astrophysical solution. However, some clues about the dark matter particle and the scatter between the properties of dwarf galaxies and their dark matter halo hosts might be found in the cosmic web of galaxies in future surveys of the local volume.

Merging haloes with similar masses (i.e., major mergers) pose significant challenges for halo finders. We compare five halo finding algorithms' (AHF, HBT, Rockstar, SubFind, and VELOCIraptor) recovery of halo properties for both isolated and cosmological major mergers. We find that halo positions and velocities are often robust, but mass biases exist for every technique. The algorithms also show strong disagreement in the prevalence and duration of major mergers, especially at high redshifts (z>1). This raises significant uncertainties for theoretical models that require major mergers for, e.g., galaxy morphology changes, size changes, or black hole growth, as well as for finding Bullet Cluster analogues. All finders not using temporal information also show host halo and subhalo relationship swaps over successive timesteps, requiring careful merger tree construction to avoid problematic mass accretion histories. We suggest that future algorithms should combine phase-space and temporal information to avoid the issues presented.

Future large-scale structure surveys of the Universe will aim to constrain the cosmological model and the true nature of dark energy with unprecedented accuracy. In order for these surveys to achieve their designed goals, they will require predictions for the nonlinear matter power spectrum to sub-percent accuracy. Through the use of a large ensemble of cosmological N-body simulations, we demonstrate that if we do not understand the uncertainties associated with simulating structure formation, i.e. knowledge of the `true' simulation parameters, and simply seek to marginalize over them, then the constraining power of such future surveys can be significantly reduced. However, for the parameters {n_s, h, Om_b, Om_m}, this effect can be largely mitigated by adding the information from a CMB experiment, like Planck. In contrast, for the amplitude of fluctuations sigma8 and the time-evolving equation of state of dark energy {w_0, w_a}, the mitigation is mild. On marginalizing over the simulation parameters, we find that the dark-energy figure of merit can be degraded by ~2. This is likely an optimistic assessment, since we do not take into account other important simulation parameters. A caveat is our assumption that the Hessian of the likelihood function does not vary significantly when moving from our adopted to the 'true' simulation parameter set. This paper therefore provides strong motivation for rigorous convergence testing of N-body codes to meet the future challenges of precision cosmology.

We present a study of the substructure finder dependence of subhalo clustering in the Aquarius Simulation. We run 11 different subhalo finders on the haloes of the Aquarius Simulation and we study their differences in the density profile, mass fraction and 2-point correlation function of subhaloes in haloes. We also study the mass and vmax dependence of subhalo clustering. As the Aquarius Simulation has been run at different resolutions, we study the convergence with higher resolutions. We find that the agreement between finders is at around the 10% level inside R200 and at intermediate resolutions when a mass threshold is applied, and better than 5% when vmax is restricted instead of mass. However, some discrepancies appear in the highest resolution, underlined by an observed resolution dependence of subhalo clustering. This dependence is stronger for the smallest subhaloes, which are more clustered in the highest resolution, due to the detection of subhaloes within subhaloes (the sub-subhalo term). This effect modifies the mass dependence of clustering in the highest resolutions. We discuss implications of our results for models of subhalo clustering and their relation with galaxy clustering.

The ever increasing size and complexity of data coming from simulations of cosmic structure formation demands equally sophisticated tools for their analysis. During the past decade, the art of object finding in these simulations has hence developed into an important discipline itself. A multitude of codes based upon a huge variety of methods and techniques have been spawned yet the question remained as to whether or not they will provide the same (physical) information about the structures of interest. Here we summarize and extent previous work of the "halo finder comparison project": we investigate in detail the (possible) origin of any deviations across finders. To this extent we decipher and discuss differences in halo finding methods, clearly separating them from the disparity in definitions of halo properties. We observe that different codes not only find different numbers of objects leading to a scatter of up to 20 per cent in the halo mass and Vmax function, but also that the particulars of those objects that are identified by all finders differ. The strength of the variation, however, depends on the property studied, e.g. the scatter in position, bulk velocity, mass, and the peak value of the rotation curve is practically below a few per cent, whereas derived quantities such as spin and shape show larger deviations. Our study indicates that the prime contribution to differences in halo properties across codes stems from the distinct particle collection methods and -- to a minor extent -- the particular aspects of how the procedure for removing unbound particles is implemented. We close with a discussion of the relevance and implications of the scatter across different codes for other fields such as semi-analytical galaxy formation models, gravitational lensing, and observables in general.

Cosmological surveys aim to use the evolution of the abundance of galaxy clusters to accurately constrain the cosmological model. In the context of LCDM, we show that it is possible to achieve the required percent level accuracy in the halo mass function with gravity-only cosmological simulations, and we provide simulation start and run parameter guidelines for doing so. Some previous works have had sufficient statistical precision, but lacked robust verification of absolute accuracy. Convergence tests of the mass function with, for example, simulation start redshift can exhibit false convergence of the mass function due to counteracting errors, potentially misleading one to infer overly optimistic estimations of simulation accuracy. Percent level accuracy is possible if initial condition particle mapping uses second order Lagrangian Perturbation Theory, and if the start epoch is between 10 and 50 expansion factors before the epoch of halo formation of interest. The mass function for halos with fewer than ~1000 particles is highly sensitive to simulation parameters and start redshift, implying a practical minimum mass resolution limit due to mass discreteness. The narrow range in converged start redshift suggests that it is not presently possible for a single simulation to capture accurately the cluster mass function while also starting early enough to model accurately the numbers of reionisation era galaxies, whose baryon feedback processes may affect later cluster properties. Ultimately, to fully exploit current and future cosmological surveys will require accurate modeling of baryon physics and observable properties, a formidable challenge for which accurate gravity-only simulations are just an initial step.

We present a study of a comparison of spin distributions of subhaloes found associated with a host halo. The subhaloes are found within two cosmological simulation families of Milky Way-like galaxies, namely the Aquarius and GHALO simulations. These two simulations use different gravity codes and cosmologies. We employ ten different substructure finders, which span a wide range of methodologies from simple overdensity in configuration space to full 6-d phase space analysis of particles.We subject the results to a common post-processing pipeline to analyse the results in a consistent manner, recovering the dimensionless spin parameter. We find that spin distribution is an excellent indicator of how well the removal of background particles (unbinding) has been carried out. We also find that the spin distribution decreases for substructure the nearer they are to the host halo's, and that the value of the spin parameter rises with enclosed mass towards the edge of the substructure. Finally subhaloes are less rotationally supported than field haloes, with the peak of the spin distribution having a lower spin parameter.

[abridged] We present a detailed comparison of fundamental dark matter halo properties retrieved by a substantial number of different halo finders. These codes span a wide range of techniques including friends-of-friends (FOF), spherical-overdensity (SO) and phase-space based algorithms. We further introduce a robust (and publicly available) suite of test scenarios that allows halo finder developers to compare the performance of their codes against those presented here. This set includes mock haloes containing various levels and distributions of substructure at a range of resolutions as well as a cosmological simulation of the large-scale structure of the universe. All the halo finding codes tested could successfully recover the spatial location of our mock haloes. They further returned lists of particles (potentially) belonging to the object that led to coinciding values for the maximum of the circular velocity profile and the radius where it is reached. All the finders based in configuration space struggled to recover substructure that was located close to the centre of the host halo and the radial dependence of the mass recovered varies from finder to finder. Those finders based in phase space could resolve central substructure although they found difficulties in accurately recovering its properties. Via a resolution study we found that most of the finders could not reliably recover substructure containing fewer than 30-40 particles. However, also here the phase space finders excelled by resolving substructure down to 10-20 particles. By comparing the halo finders using a high resolution cosmological volume we found that they agree remarkably well on fundamental properties of astrophysical significance (e.g. mass, position, velocity, and peak of the rotation curve).

The velocity distribution function of dark matter particles is expected to show significant departures from a Maxwell-Boltzmann distribution. This can have profound effects on the predicted dark matter - nucleon scattering rates in direct detection experiments, especially for dark matter models in which the scattering is sensitive to the high velocity tail of the distribution, such as inelastic dark matter (iDM) or light (few GeV) dark matter (LDM), and for experiments that require high energy recoil events, such as many directionally sensitive experiments. Here we determine the velocity distribution functions from two of the highest resolution numerical simulations of Galactic dark matter structure (Via Lactea II and GHALO), and study the effects for these scenarios. For directional detection, we find that the observed departures from Maxwell-Boltzmann increase the contrast of the signal and change the typical direction of incoming DM particles. For iDM, the expected signals at direct detection experiments are changed dramatically: the annual modulation can be enhanced by more than a factor two, and the relative rates of DAMA compared to CDMS can change by an order of magnitude, while those compared to CRESST can change by a factor of two. The spectrum of the signal can also change dramatically, with many features arising due to substructure. For LDM the spectral effects are smaller, but changes do arise that improve the compatibility with existing experiments. We find that the phase of the modulation can depend upon energy, which would help discriminate against background should it be found.

We use the recently completed one billion particle Via Lactea II LambdaCDM simulation to investigate local properties like density, mean velocity, velocity dispersion, anisotropy, orientation and shape of the velocity dispersion ellipsoid, as well as structure in velocity space of dark matter haloes. We show that at the same radial distance from the halo centre, these properties can deviate by orders of magnitude from the canonical, spherically averaged values, a variation that can only be partly explained by triaxiality and the presence of subhaloes. The mass density appears smooth in the central relaxed regions but spans four orders of magnitude in the outskirts, both because of the presence of subhaloes as well as of underdense regions and holes in the matter distribution. In the inner regions the local velocity dispersion ellipsoid is aligned with the shape ellipsoid of the halo. This is not true in the outer parts where the orientation becomes more isotropic. The clumpy structure in local velocity space of the outer halo can not be well described by a smooth multivariate normal distribution. Via Lactea II also shows the presence of cold streams made visible by their high 6D phase space density. Generally, the structure of dark matter haloes shows a high degree of graininess in phase space that cannot be described by a smooth distribution function.

Our recently completed one billion particle Via Lactea II simulation of a Milky Way-sized dark matter halo resolves over 50,000 gravitationally bound clumps orbiting today within the virialized region of the main host. About 2,300 of these subhalos have one or more "progenitors" above 1e6 Msun at redshift 11, i.e. massive enough for their gas to have cooled via excitation of H2 and fragmented prior to the epoch of cosmic reionization. We count 4,500 such progenitors: if these were able to convert a fraction of their gas content into very metal-poor stars with a Salpeter initial mass function (IMF), they would be shining today with a visual magnitude M_V=6.7 per solar mass in stars. Assuming a universal baryon fraction, we show that mean star formation efficiencies as low as 0.1% in progenitors below 1e8 Msun would overproduce the abundance of the faint Galatic dwarf spheroidals observed by the Sloan Digital Sky Survey. Star formation at first light must have occurred either with an IMF lacking stars below 0.9 Msun, or was intrinsically very inefficient in small dark matter halos. If the latter, our results may be viewed as another hint of a minimum scale in galaxy formation.

We perform a series of simulations of a Galactic mass dark matter halo at different resolutions, our largest uses over three billion particles and has a mass resolution of 1000 M_sun. We quantify the structural properties of the inner dark matter distribution and study how they depend on numerical resolution. We can measure the density profile to a distance of 120 pc (0.05% of R_vir) where the logarithmic slope is -0.8 and -1.4 at (0.5% of R_vir). We propose a new two parameter fitting function that has a linearly varying logarithmic density gradient which fits the GHALO and VL2 density profiles extremely well. Convergence in the density profile and the halo shape scales as N^(-1/3), but the shape converges at a radius three times larger at which point the halo becomes more spherical due to numerical resolution. The six dimensional phase-space profile is dominated by the presence of the substructures and does not follow a power law, except in the smooth under-resolved inner few kpc.

In cold dark matter cosmological models, structures form and grow by merging of smaller units. Numerical simulations have shown that such merging is incomplete; the inner cores of halos survive and orbit as "subhalos" within their hosts. Here we report a simulation that resolves such substructure even in the very inner regions of the Galactic halo. We find hundreds of very concentrated dark matter clumps surviving near the solar circle, as well as numerous cold streams. The simulation reveals the fractal nature of dark matter clustering: Isolated halos and subhalos contain the same relative amount of substructure and both have cuspy inner density profiles. The inner mass and phase-space densities of subhalos match those of recently discovered faint, dark matter-dominated dwarf satellite galaxies and the overall amount of substructure can explain the anomalous flux ratios seen in strong gravitational lenses. Subhalos boost gamma-ray production from dark matter annihilation, by factors of 4-15, relative to smooth galactic models. Local cosmic ray production is also enhanced, typically by a factor 1.4, but by more than a factor of ten in one percent of locations lying sufficiently close to a large subhalo. These estimates assume that gravitational effects of baryons on dark matter substructure are small.

We use a series of cosmological N-body simulations for a flat LCDM cosmology to investigate the properties of dark matter haloes in the mass range 3.0e9-3.0e13 Msun. These properties include the concentration parameter (c), the spin parameter (lambda) and the mean axis ratio (q). For the concentration-mass relation we find c~M^(-0.11) in agreement with the model proposed by Bullock et al. even if we find a lower normalization (15%). The results for lambda and q are in good agreement with previous studies, while c and lambda are anti-correlated. In an attempt to remove unrelaxed haloes, we use the offset parameter (xoff), defined as the distance between the most bound particle and the center of mass. Removing haloes with large xoff increases the c by ~10%, lowers the lambda by ~15%, and removes the most prolate haloes. In addition, it largely removes the anti-correlation between c and lambda though not entirely. We also investigate the effects of the large-scale environment. We find that more concentrated haloes live in denser environments. Note, however, that the trend is weak compared to the scatter. For the spin parameters we find no environment dependence, while there is a weak indication that the most spherical haloes reside in denser region. Finally, using a simple model for disk galaxy formation we show that haloes that host low surface brightness galaxies are expected to be hosted by a biased sub-set of haloes. Not only do these haloes have spin parameters that are larger than average, they also have c that are 15% percent lower than the average at a given halo mass. We discuss the implications of all these findings for the claimed disagreement between halo concentrations inferred from LSB rotation curves, and those expected for a LCDM cosmology. (abridged)

N-body simulations sample their initial conditions on an initial particle distribution, which for cosmological simulations is usually a glass or grid, whilst a Poisson distribution is used for galaxy models, spherical collapse etc. These pre-initial conditions have inherent correlations, noise due to discreteness and preferential alignments, whilst the glass distribution is poorly defined and computationally expensive to construct. We present a novel particle distribution which can be useful as a pre-initial condition for N-body simulations, using a simple construction based on a ``quaquaversal'' tiling of space. This distribution has little preferred orientation (i.e. is statistically isotropic), has a rapidly vanishing large scale power-spectrum (P(k) ~ k^4), and is trivial to create. It should be particularly useful for warm dark matter and cold collapse simulations.

We have carried out a hydrodynamical code comparison study of interacting multiphase fluids. The two commonly used techniques of grid and smoothed particle hydrodynamics (SPH) show striking differences in their ability to model processes that are fundamentally important across many areas of astrophysics. Whilst Eulerian grid based methods are able to resolve and treat important dynamical instabilities, such as Kelvin-Helmholtz or Rayleigh-Taylor, these processes are poorly or not at all resolved by existing SPH techniques. We show that the reason for this is that SPH, at least in its standard implementation, introduces spurious pressure forces on particles in regions where there are steep density gradients. This results in a boundary gap of the size of the SPH smoothing kernel over which information is not transferred.

A possible explanation of the present discrepancy between the abundance of galactic subhaloes predicted by N-Body simulations with those observed in the Local Group is presented. We study the impact of an early reionization on the baryonic component of the Universe using SPH simulations on group and galactic scales. We use a simplified model for reionization described as an instantaneous increment of the IGM temperature (with 1e5<T_{IGM}<1e6 K) at a given redshift (z_r > 8). We find that a key role is played by compton cooling (interaction between hot electrons and the CMB photons); at high redshift (z_r>10) this cooling is very efficient and it is able to counteract any heating of the gas within few Myrs. This means that a late reionization is needed (z_r<9) to sufficiently reduce the number of luminous dwarf satellites around our Galaxy. For a reionization redshift z_r=8 and a reionization temperature of T_{IGM}~1e5 K we are able to reproduce the observed number of Local Group dwarf galaxies in our simulations.

We assess the constraints on the evolutionary models of young low-mass objects that are provided by the measurements of the companion AB Dor C by Close and coworkers and by a new comparison of model-derived IMFs of star-forming regions to the well-calibrated IMF of the solar neighborhood. After performing an independent analysis of Close's imaging and spectroscopic data for AB Dor C, we find that AB Dor C is not detected at a significant level (SN 1.2) in the SDI images when one narrow-band image is subtracted from another, but that it does appear in the individual SDI frames as well as the images at JHK. Using the age of 75-150 Myr for AB Dor from Luhman, Stauffer, & Mamajek, the luminosity predicted by the models of Chabrier & Baraffe is consistent with the value that we estimate. We measure a spectral type of M6+/-1 from the K-band spectrum of AB Dor C, which is earlier than the value of M8+/-1 from Close and is consistent with the model predictions when a dwarf temperature scale is adopted. In a test of these models at much younger ages, we show that the low-mass IMFs that they produce for star-forming regions are similar to the IMF of the solar neighborhood. If the masses of the low-mass stars and brown dwarfs in these IMFs of star-forming regions were underestimated by a factor of two as suggested by Close, then the IMF characterizing the current generation of Galactic star formation would have to be radically different from the IMF of the solar neighborhood.

We present high resolution (~0.1"), very high Strehl ratio (0.97+-0.03) mid-infrared (IR) adaptive optics (AO) images of the AGB star RV Boo utilizing the MMT adaptive secondary AO system. RV Boo was observed at a number of wavelengths over two epochs (9.8 um in May 2003, 8.8, 9.8 and 11.7 um in February 2004) and appeared slightly extended at all wavelengths. While the extension is very slight at 8.8 and 11.7 um data, the extension is somewhat more pronounced at 9.8 um. With such high Strehls we can achieve super-resolutions of 0.1" by deconvolving RV Boo with a point-spread function (PSF) derived from an unresolved star. We tentatively resolve RV Boo into a 0.16" FWHM extension at a position angle of 120 degrees. At a distance of 390(+250)(-100) pc, this corresponds to a FWHM of 60(+40)(-15) AU. We measure a total flux at 9.8 um of 145+-24 Jy for the disk and star. Based on a dust thermal emission model for the observed IR spectral energy distribution and the 9.8 um AO image, we derive a disk dust mass of 1.6x10^-6 Msun and an inclination of 30 to 45 degrees from edge-on. We discuss whether the dust disk observed around RV Boo is an example of the early stages in the formation of asymmetric structure in planetary nebula.

We present high-contrast imaging data on the disk of the classical T Tauri star TW Hya. The images were obtained through the polarimetric differential imaging technique with the adaptive optics system NACO. Our commissioning data show the presence of polarized disk emission between 0.1" and 1.4" from the star. We derive the first Ks-band radial polarized intensity distribution. We show that the polarized intensity compares well to shorter wavelengths surface brightness observations and confirm the previously reported gradual slope change around 0.8". These results show the potential of the new polarimetric differential imaging technique at 8m-class telescopes to map the inner regions of protoplanetary disks.

We report the discovery of two faint ultracool companions to the nearby (d~17.9 pc) young G2V star HD 130948 (HR 5534, HIP 72567) using the Hokupa'a adaptive optics instrument mounted on the Gemini North 8-meter telescope. Both objects have the same common proper motion as the primary star as seen over a 7 month baseline and have near-IR photometric colors that are consistent with an early-L classification. Near-IR spectra taken with the NIRSPEC AO instrument on the Keck II telescope reveal K I lines, FeH, and water bandheads. Based on these spectra, we determine that both objects have spectral type dL2 with an uncertainty of 2 spectral subclasses. The position of the new companions on the H-R diagram in comparison with theoretical models is consistent with the young age of the primary star (<0.8 Gyr) estimated on the basis of X-ray activity, lithium abundance and fast rotation. HD 130948 B and C likely constitute a pair of young contracting brown dwarfs with an orbital period of about 10 years, and will yield dynamical masses for L dwarfs in the near future.

Use of the highly sensitive Hokupa'a/Gemini curvature wavefront sensor has allowed for the first time direct adaptive optics (AO) guiding on M8-M9 very low mass (VLM) stars. An initial survey of 20 such objects (SpT=M8-M9) discovered 4 binaries. Three of the systems have separations of less than 4.2 AU and similar mass ratios (Delta K<0.8 mag; 0.85<q<1.0). One system, however, did have the largest Delta K=2.38 mag and sep=14.4 AU yet observed for a VLM star with a brown dwarf companion. Based on our initial flux limited (Ks<12 mag) survey of 20 M8-M9 stars over 14:26<RA<4:30 hours from the sample of Gizis et al. (2000) we find a binary fraction in the range 14-24% for M8-M9 binaries with sep>3 AU. This is likely consistent with the 23+/-5% measured for more massive (M0-M6) stars over the same separation range. It appears M8-M9 binaries have a much smaller semi-major axis distribution peak (~4 AU; with no systems wider than 15 AU) compared to M and G stars which have a broad peak at larger \~30 AU separations.

Use of the highly sensitive Hokupa'a curvature wavefront sensor has allowed for the first time direct adaptive optics (AO) guiding on brown dwarfs and VLM stars (SpT=M7-L2). An initial survey of 9 such objects discovered one 0.15" binary (2MASSJ 1426316+155701). The companion is about half as bright as the primary (Delta K = 0.61+/-0.05$, Delta H = 0.70+/-0.05) and has even redder colors H-K=0.59+/-0.14 than the primary. The blended spectrum of the binary has been previously determined to be M9.0. We modeled a blend of an M8.5 template and a L1-L3 template reproducing a M9.0 spectrum in the case of Delta K = 0.61+/-0.05,Delta H = 0.70\pm0.05$. These spectral types also match the observed H-K colors of each star. Based the previously observed low space motion and $H_{\alpha}$ activity we assign an age of $0.8^{+6.7}_{-0.3} Gyr$. Utilizing this age range and the latest DUSTY models of the Lyon group we assign a photometric distance of $18.8^{+1.44}_{-1.02} pc$ and masses of $M_{A}=0.074^{+0.005}_{-0.011} M_\odot$ and $M_{B}=0.066^{+0.006}_{-0.015} M_\odot$. We therefore estimate a system separation of $2.92_{+0.22}^{-0.16}AU$ and a period of $13.3{+3.18}^{-1.51} yr$ respectively. Hence, 2M1426 is among the smallest separation brown dwarf binaries resolved to date.

There is strong evidence that the planets in the solar system evolved from a disk-shaped solar nebula ~4.56 Gyr ago. By studying young stars in various evolutionary stages, one aims at tracing back the early history of the solar system, in particular the timescales for disk dissipation and for the formation of planetary systems. We used the VLT & ISAAC, and ESA's Infrared Space Observatory & ISOCAM to study the circumstellar environment of young low-mass stars.

Many extra-solar planets and a few planetary systems have been found indirectly by small periodic radial velocity variations around old nearby stars. The orbital characteristics of most of them are different from the planets in our solar system. Hence, planet formation theories have to be revised. Therefore, observational constraints regarding young planets would be very valuable. We have started a ground-based direct imaging search for giant planets in orbit around young nearby stars. Here, we will motivate the sample selection and will present our direct imaging observation of the very low-mass (15 to 40 Jupiter masses) brown dwarf companion TWA-5 B in orbit around the nearby young star TWA-5 A, recently obtained with the 36-element curvature-sensing AO instrument Hokupa'a of the University of Hawai'i at the 8.3m Gemini-North telescope on Mauna Kea. We could achieve a FWHM of 64 mas and 25 % Strehl. We find significance evidence for orbital motion of B around A.

Few circumstellar disks have been directly observed. Here we use sensitive differential polarimetric techniques to overcome atmospheric speckle noise in order to image the circumstellar material around HD 169142. The detected envelope or disk is considerably smaller than expectations based on the measured strength of the far-IR excess from this system.

The bright stellar content near the center of the Local Group elliptical galaxy M32 is investigated with 0.12 arcsec FWHM H and K images obtained with the Gemini Mauna Kea telescope. Stars with K = 15.5, which are likely evolving near the tip of the asymptotic giant branch (AGB), are resolved to within 2 arcsec of the nucleus, and it is concluded that the peak stellar brightness near the center of M32 is similar to that in the outer regions of the galaxy. Moreover, the projected density of bright AGB stars follows the visible light profile to within 2 arcsec of the nucleus, indicating that the brightest stars are well mixed throughout the galaxy. Thus, there is no evidence for an age gradient, and the radial variations in spectroscopic indices and ultraviolet colors that have been detected previously must be due to metallicity and/or some other parameter. We suggest that either the bright AGB stars formed as part of a highly uniform and coherent galaxy-wide episode of star formation, or they originated in a separate system that merged with M32.