Low-scale baryogenesis and dark matter generation can occur via the production of neutral $B$ mesons at MeV temperatures in the early Universe, which undergo CP-violating oscillations and subsequently decay into a dark sector. In this work, we discuss the consequences of realizing this mechanism in a supersymmetric model with an unbroken $U(1)_R$ symmetry which is identified with baryon number. $B$ mesons decay into a dark sector through a baryon number conserving operator mediated by TeV scale squarks and a GeV scale Dirac bino. The dark sector particles can be identified with sterile neutrinos and their superpartners in a type-I seesaw framework for neutrino masses. The sterile sneutrinos are sufficiently long lived and constitute the dark matter. The produced matter-antimatter asymmetry is directly related to observables measurable at $B$ factories and hadron colliders, the most relevant of which are the semileptonic-leptonic asymmetries in neutral $B$ meson systems and the inclusive branching fraction of $B$ mesons into hadrons and missing energy. We discuss model independent constraints on these experimental observables before quoting predictions made in the supersymmetric context. Constraints from astrophysics, neutrino physics and flavor observables are studied, as are potential LHC signals with a focus on novel long lived particle searches which are directly linked to properties of the dark sector.

emcee is a Python library implementing a class of affine-invariant ensemble samplers for Markov chain Monte Carlo (MCMC). This package has been widely applied to probabilistic modeling problems in astrophysics where it was originally published, with some applications in other fields. When it was first released in 2012, the interface implemented in emcee was fundamentally different from the MCMC libraries that were popular at the time, such as PyMC, because it was specifically designed to work with "black box" models instead of structured graphical models. This has been a popular interface for applications in astrophysics because it is often non-trivial to implement realistic physics within the modeling frameworks required by other libraries. Since emcee's release, other libraries have been developed with similar interfaces, such as dynesty (Speagle 2019). The version 3.0 release of emcee is the first major release of the library in about 6 years and it includes a full re-write of the computational backend, several commonly requested features, and a set of new "move" implementations.

We present a new mechanism of Baryogenesis and dark matter production in which both the dark matter relic abundance and the baryon asymmetry arise from neutral $B$ meson oscillations and subsequent decays. This set-up is testable at hadron colliders and $B$-factories. In the early Universe, decays of a long lived particle produce $B$ mesons and anti-mesons out of thermal equilibrium. These mesons/anti-mesons then undergo CP violating oscillations before quickly decaying into visible and dark sector particles. Dark matter will be charged under Baryon number so that the visible sector baryon asymmetry is produced without violating the total baryon number of the Universe. The produced baryon asymmetry will be directly related to the leptonic charge asymmetry in neutral $B$ decays; an experimental observable. Dark matter is stabilized by an unbroken discrete symmetry, and proton decay is simply evaded by kinematics. We will illustrate this mechanism with a model that is unconstrained by di-nucleon decay, does not require a high reheat temperature, and would have unique experimental signals -- a positive leptonic asymmetry in $B$ meson decays, a new decay of $B$ mesons into a baryon and missing energy, and a new decay of $b$-flavored baryons into mesons and missing energy. These three observables are testable at current and upcoming collider experiments, allowing for a distinct probe of this mechanism.

The default assumption of early universe cosmology is that the postinflationary universe was radiation dominated until it was about 47000 years old. Direct evidence for the radiation dominated epoch extends back until nucleosynthesis, which began during the first second. However there are theoretical reasons to prefer a period of earlier matter domination, prior to nucleosynthesis, e.g. due to late decaying massive particles needed to explain baryogenesis. Axion cosmology is quantitatively affected by an early period of matter domination, with a different axion mass range preferred and greater inhomogeneity produced on small scales. In this work we show that such increased inhomogeneity can lead to the formation of axion miniclusters in axion parameter ranges that are different from those usually assumed. If the reheating temperature is below $58$ MeV, axion miniclusters can form even if the axion field is present during inflation and has been previously homogenized. The upper bound on the typical initial axion minicluster mass is raised from $10^{-10} M_{\odot}$ to $10^{-7} M_{\odot}$, where $M_{\odot}$ is a solar mass. These results may have consequences for indirect detection of axion miniclusters, and could conceivably probe the thermal history of the universe before nucleosynthesis.

Exotic new particles carrying baryon number and with mass of order the nucleon mass have been proposed for various reasons including baryogenesis, dark matter, mirror worlds, and the neutron lifetime puzzle. We show that the existence of neutron stars with mass greater than 0.7 $M_\odot$ places severe constraints on such particles, requiring them to be heavier than 1.2 GeV or to have strongly repulsive self-interactions.

The relaxion mechanism provides a potentially elegant solution to the hierarchy problem without resorting to anthropic or other fine-tuning arguments. This mechanism introduces an axion-like field, dubbed the relaxion, whose expectation value determines the electroweak hierarchy as well as the QCD strong CP violating $\bar{\theta}$ parameter. During an inflationary period, the Higgs mass squared is selected to be negative and hierarchically small in a theory which is consistent with 't Hooft's technical naturalness criteria. However, in the original model proposed by Graham, Kaplan and Rajendran (2015), the relaxion does not solve the strong CP problem, and in fact contributes to it, as the coupling of the relaxion to the Higgs field and the introduction of a linear potential for the relaxion produces large strong CP violation. We resolve this tension by considering inflation with a Hubble scale which is above the QCD scale but below the weak scale, and estimating the Hubble temperature dependence of the axion mass. The relaxion potential is thus very different during inflation than it is today. We find that provided the inflationary Hubble scale is between the weak scale and about 3 GeV, the relaxion resolves the hierarchy, strong CP, and dark matter problems in a way that is technically natural.

Neutrino masses and the number of light neutrino species can be tested in a variety of laboratory experiments and also can be constrained by particle astrophysics and precision cosmology. A conflict between these various results could be an indication of new physics in the neutrino sector. In this paper we explore the possibility for reconciliation of otherwise discrepant results in a simple model containing a light scalar field which produces Mass Varying Neutrinos (MaVaNs). We extend previous work on MaVaNs to consider issues of neutrino clumping, the effects of additional contributions to neutrino mass, and reconciliation of eV mass sterile neutrinos with cosmology.

(abridged) We present 2-dimensional hydrodynamic simulations using the Smoothed Particle Hydrodynamic (SPH) code, VINE, to model a self-gravitating binary system similar to the GG Tau A system. We simulate systems configured with semi-major axes of either $a=62$~AU (`wide') or $a=32$~AU (`close'), and with eccentricity of either $e=0$ or $e=0.3$. Strong spiral structures are generated with large material streams extending inwards. A small fraction accretes onto the circumstellar disks, with most returning to the torus. Structures also propagate outwards, generating net outwards mass flow and eventually losing coherence at large distances. The torus becomes significantly eccentric in shape. Accretion onto the stars occurs at a rate of a few $\times10^{-8}$\msun/yr implying disk lifetimes shorter than $\sim10^4$~yr, without replenishment. Only wide configurations retain disks by virtue of robust accretion. In eccentric configurations, accretion is episodic, occurs preferentially onto the secondary at wrates peaked near binary periapse. We conclude that the \ggtaua\ torus is strongly self gravitating and that a major contribution to its thermal energy is shock dissipation. We interpret its observed features as manifestations of spiral structures and the low density material surrounding it as an excretion disk created by outward mass flux. We interpret GG Tau A as a coplanar system with an eccentric torus, and account for its supposed mutual inclination as due to degeneracy between the interpretation of inclination and eccentricity. Although the disks persist for long enough to permit planet formation, the environment remains unfavorable due to high temperatures. We conclude that the GG Tau A system is in an eccentric, $a\sim62$~AU orbit.

Liquid argon is a bright scintillator with potent particle identification properties, making it an attractive target for direct-detection dark matter searches. The DarkSide-50 dark matter search here reports the first WIMP search results obtained using a target of low-radioactivity argon. DarkSide-50 is a dark matter detector, using two-phase liquid argon time projection chamber, located at the Laboratori Nazionali del Gran Sasso. The underground argon is shown to contain Ar-39 at a level reduced by a factor (1.4 +- 0.2) x 10^3 relative to atmospheric argon. We report a background-free null result from (2616 +- 43) kg d of data, accumulated over 70.9 live-days. When combined with our previous search using an atmospheric argon, the 90 % C.L. upper limit on the WIMP-nucleon spin-independent cross section based on zero events found in the WIMP search regions, is 2.0 x 10^-44 cm^2 (8.6 x 10^-44 cm^2, 8.0 x 10^-43 cm^2) for a WIMP mass of 100 GeV/c^2 (1 TeV/c^2 , 10 TeV/c^2).

We propose a new mechanism for baryogenesis at the 1-200 MeV scale. Enhancement of CP violation takes place via interference between oscillations and decays of mesinos--bound states of a scalar quark and antiquark and their CP conjugates. We present the mechanism in a simplified model with four new fundamental particles, with masses between 300 GeV and 10 TeV, and show that some of the experimentally allowed parameter space can give the observed baryon-to-entropy ratio.

We report the first results of DarkSide-50, a direct search for dark matter operating in the underground Laboratori Nazionali del Gran Sasso (LNGS) and searching for the rare nuclear recoils possibly induced by weakly interacting massive particles (WIMPs). The dark matter detector is a Liquid Argon Time Projection Chamber with a (46.4+-0.7) kg active mass, operated inside a 30 t organic liquid scintillator neutron veto, which is in turn installed at the center of a 1 kt water Cherenkov veto for the residual flux of cosmic rays. We report here the null results of a dark matter search for a (1422+-67) kg d exposure with an atmospheric argon fill. This is the most sensitive dark matter search performed with an argon target, corresponding to a 90% CL upper limit on the WIMP-nucleon spin-independent cross section of 6.1x10^-44 cm^2 for a WIMP mass of 100 GeV/c^2.

It is generally inferred from astronomical measurements that Dark Matter (DM) comprises approximately 27\% of the energy-density of the universe. If DM is a subatomic particle, a possible candidate is a Weakly Interacting Massive Particle (WIMP), and the DarkSide-50 (DS) experiment is a direct search for evidence of WIMP-nuclear collisions. DS is located underground at the Laboratori Nazionali del Gran Sasso (LNGS) in Italy, and consists of three active, embedded components; an outer water veto (CTF), a liquid scintillator veto (LSV), and a liquid argon (LAr) time projection chamber (TPC). This paper describes the data acquisition and electronic systems of the DS detectors, designed to detect the residual ionization from such collisions.

Cold dark matter explains a wide range of data on cosmological scales. However, there has been a steady accumulation of evidence for discrepancies between simulations and observations at scales smaller than galaxy clusters. Solutions to these small scale structure problems may indicate that simulations need to improve how they include feedback from baryonic matter, or may imply that dark matter properties differ from the standard cold, noninteracting scenario. One promising way to affect structure formation on small scales is a relatively strong coupling of dark matter to neutrinos. We construct an experimentally viable, simple, renormalizable, model with new interactions between neutrinos and dark matter. We show that addressing the small scale structure problems requires dark matter with a mass that is tens of MeV, and a present-day density determined by an initial particle-antiparticle asymmetry in the dark sector. Generating a sufficiently large dark matter-neutrino coupling requires a new heavy neutrino with a mass around 100 MeV. The heavy neutrino is mostly sterile but has a substantial $\tau$ neutrino component, while the three nearly massless neutrinos are partly sterile. We provide the first discussion of how such dark matter-neutrino interactions affect neutrino (especially $\tau$ neutrino) phenomenology. This model can be tested by future astrophysical, particle physics, and neutrino oscillation data. A feature in the neutrino energy spectrum and flavor content from a future nearby supernova would provide strong evidence of neutrino-dark matter interactions. Promising signatures include anomalous matter effects in neutrino oscillations due to nonstandard interactions and a component of the $\tau$ neutrino with mass around 100 MeV.

Evidence for an excess of gamma rays with O(GeV) energy coming from the center of our galaxy has been steadily accumulating over the past several years. Recent studies of the excess in data from the Fermi telescope have cast doubt on an explanation for the excess arising from unknown astrophysical sources. A potential source of the excess is the annihilation of dark matter into standard model final states, giving rise to gamma ray production. The spectrum of the excess is well fit by 30 GeV dark matter annihilating into a pair of b quarks with a cross section of the same order of magnitude as expected for a thermal relic. Simple models that can lead to this annihilation channel for dark matter are in strong tension with null results from direct detection experiments. We construct a renormalizable model where dark matter-standard model interactions are mediated by a pseudoscalar that mixes with the CP-odd component of a pair of Higgs doublets, allowing for the gamma ray excess to be explained while suppressing the direct detection signal. We consider implications for this scenario from Higgs decays, rare B meson decays and monojet searches and also comment on some difficulties that any dark matter model explaining the gamma ray excess via direct annihilation into quarks will encounter.

In models with a light scalar field (the `acceleron') coupled to neutrinos, neutrino masses depend on neutrino density. The resulting coupled system of mass varying neutrinos (MaVaNs) and the acceleron can act as a negative pressure fluid and is a candidate for dark energy \cite{Fardon:2003eh} . MaVaNs also allow for higher $\Sigma$m$_\nu$ than terrestrial bounds, giving late forming warm dark matter. In this paper we study the effect of increasing neutrino mass on the CMB spectrum, implementing MaVaNs cosmology using CMBEASY. We find that the CMB spectrum is only affected at very low multipoles. Cosmic variance allows for significant warm dark matter at late times. This implies that in MaVaNs cosmology $\sigma_8$ as determined by the CMB may not be a good determinant of structure evolution at late times, potentially reducing the tension between $\sigma$$_{8}$ as reported by Planck Collaboration \cite{PlanckSZ:2013} without increasing the tension in the Planck determined value of Hubble's constant. In addition, in MaVaNs theories, CMB data do not necessarily constrain possible neutrino mass results in terrestrial experiments.

These reports present the results of the 2013 Community Summer Study of the APS Division of Particles and Fields ("Snowmass 2013") on the future program of particle physics in the U.S. Chapter 4, on the Cosmic Frontier, discusses the program of research relevant to cosmology and the early universe. This area includes the study of dark matter and the search for its particle nature, the study of dark energy and inflation, and cosmic probes of fundamental symmetries.

We study how many-body effects alter the dark matter (DM) thermalization time inside neutron stars. We find that Pauli blocking, kinematic constraints, and superfluidity and superconductivity in the neutron star significantly affect the DM thermalization time, in general lengthening it. This could change the final DM mass and DM-nucleon cross section constraints by considering black hole formation in neutron stars due to DM accretion. We consider the class of models in which DM is an asymmetric, complex scalar particle with a mass between 1 keV and 5 GeV which couples to regular matter via some heavy vector boson. Interestingly, we find that the discovery of asymmetric, bosonic DM could motivate the existence of exotic neutron star cores. We apply our results to the case of mixed sneutrino DM.

Report of the CF6 Working Group at Snowmass 2013. Topics addressed include ultra-high energy cosmic rays, neutrinos, gamma rays, baryogenesis, and experiments probing the fundamental nature of spacetime.

As part of the DarkSide program of direct dark matter searches using liquid argon TPCs, a prototype detector with an active volume containing 10 kg of liquid argon, DarkSide-10, was built and operated underground in the Gran Sasso National Laboratory in Italy. A critically important parameter for such devices is the scintillation light yield, as photon statistics limits the rejection of electron-recoil backgrounds by pulse shape discrimination. We have measured the light yield of DarkSide-10 using the readily-identifiable full-absorption peaks from gamma ray sources combined with single-photoelectron calibrations using low-occupancy laser pulses. For gamma lines of energies in the range 122-1275 keV, we get consistent light yields averaging 8.887+-0.003(stat)+-0.444(sys) p.e./keVee. With additional purification, the light yield measured at 511 keV increased to 9.142+-0.006(stat) p.e./keVee.

(shortened) We perform 3D hydrodynamic simulations of gas flowing around a planetary core of mass \mplan=10\me embedded in a near Keplerian background flow, using a modified shearing box approximation. We employ a nested grid hydrodynamic code with as many as six nested grids, providing spatial resolution on the finest grid comparable to the present day diameters of Neptune and Uranus. We find that a strongly dynamically active flow develops such that no static envelope can form. The activity is not sensitive to plausible variations in the rotation curve of the underlying disk. It is sensitive to the thermodynamic treatment of the gas, as modeled by prescribed equations of state (either `locally isothermal' or `locally isentropic') and the temperature of the background disk material. The activity is also sensitive to the shape and depth of the core's gravitational potential, through its mass and gravitational softening coefficient. The varying flow pattern gives rise to large, irregular eruptions of matter from the region around the core which return matter to the background flow: mass in the envelope at one time may not be found in the envelope at any later time. The angular momentum of material in the envelope, relative to the core, varies both in magnitude and in sign on time scales of days to months near the core and on time scales a few years at distances comparable to the Hill radius. We show that material entering the dynamically active environment may suffer intense heating and cooling events the durations of which are as short as a few hours to a few days. Peak temperatures in these events range from $T \sim 1000$ K to as high as $T \sim 3-4000$ K, with densities $\rho\sim 10^{-9}-10^{-8}$ g/cm$^3$. These time scales, densities and temperatures span a range consistent with those required for chondrule formation in the nebular shock model.

We report on the first large-scale production of low radioactivity argon from underground gas wells. Low radioactivity argon is of general interest, in particular for the construction of large scale WIMP dark matter searches and detectors of reactor neutrinos for non-proliferation efforts. Atmospheric argon has an activity of about 1 Bq/kg from the decays of 39Ar; the concentration of 39Ar in the underground argon we are collecting is at least a factor of 100 lower than this value. The argon is collected from a stream of gas from a CO2 well in southwestern Colorado with a Vacuum Pressure Swing Adsorption (VPSA) plant. The gas from the well contains argon at a concentration of 400-600 ppm, and the VPSA plant produces an output stream with an argon concentration at the level of 30,000-50,000 ppm (3-5%) in a single pass. This gas is sent for further processing to Fermilab where it is purified by cryogenic distillation. The argon production rate is presently 0.5 kg/day.

This white paper addresses the hypothesis of light sterile neutrinos based on recent anomalies observed in neutrino experiments and the latest astrophysical data.

We explore the possibility that the dark matter is a condensate of a very light vector boson. Such a condensate could be produced during inflation, provided the vector mass arises via the Steuckelberg mechanism. We derive bounds on the kinetic mixing of the dark matter boson with the photon, and point out several potential signatures of this model.

We present a simple model in which dark matter couples to the standard model through a light scalar intermediary that is itself unstable. We find this model has several notable features, and allows a natural explanation for a surplus of positrons, but no surplus of anti-protons, as has been suggested by early data from PAMELA and ATIC. Moreover, this model yields a very small nucleon coupling, well below the direct detection limits. In this paper we explore the effect of this model in both the early universe and in the galaxy.

A new harassment model is presented that models the complex, and dynamical tidal field of a Virgo like galaxy cluster. The model is applied to small, late-type dwarf disc galaxies (of substantially lower mass than in previous harassment simulations) as they infall into the cluster from the outskirts. These dwarf galaxies are only mildly affected by high speed tidal encounters with little or no observable consequences; typical stellar losses are $<10\%$, producing very low surface brightness streams ($\mu_B > 31$ mag arcsec$^{-2}$), and a factor of two drop in dynamical mass-to-light ratio. Final stellar discs remain disc-like, and dominated by rotation although often with tidally induced spiral structure. By means of Monte-Carlo simulations, the statistically likely influences of harassment on infalling dwarf galaxies are determined. The effects of harassment are found to be highly dependent on the orbit of the galaxy within the cluster, such that newly accreted dwarf galaxies typically suffer only mild harassment. Strong tidal encounters, that can morphologically transform discs into spheroidals, are rare occurring in $<15 \%$ of dwarf galaxy infalls for typical orbits of sub-structure within $\Lambda$CDM cluster mass halos. For orbits with small apocentric distances ($<$250 kpc), harassment is significantly stronger resulting in complete disruption or heavy mass loss ($>90 \%$ dark matter and $> 50 \%$ stellar), however, such orbits are expected to be highly improbable for newly infalling galaxies due to the deep potential well of the cluster.

We investigate the dynamical evolution of a Jovian--mass planet injected into an orbit highly inclined with respect to its nesting gaseous disk. Planet--planet scattering induced by convergent planetary migration and mean motion resonances may push a planet into such an out of plane configuration with inclinations as large as $20^\circ-30^\circ$. In this scenario the tidal interaction of the planet with the disk is more complex and, in addition to the usual Lindblad and corotation resonances, it involves also inclination resonances responsible of bending waves. We have performed three--dimensional hydrodynamic simulations of the disk and of its interactions with the planet with a Smoothed Particle Hydrodynamics (SPH) code. A main result is that the initial large eccentricity and inclination of the planetary orbit are rapidly damped on a timescale of the order of $10^3$ yrs, almost independently of the initial semimajor axis and eccentricity of the planet. The disk is warped in response to the planet perturbations and it precesses. Inward migration occurs also when the planet is inclined and it has a drift rate which is intermediate between type I and type II migration. The planet is not able to open a gap until its inclination becomes lower than $\sim 10^\circ$ when it also begins to accrete a significant amount of mass from the disk.

We continue our presentation of VINE. We begin with a description of relevant architectural properties of the serial and shared memory parallel computers on which VINE is intended to run, and describe their influences on the design of the code itself. We continue with a detailed description of a number of optimizations made to the layout of the particle data in memory and to our implementation of a binary tree used to access that data for use in gravitational force calculations and searches for SPH neighbor particles. We describe modifications to the code necessary to obtain forces efficiently from special purpose `GRAPE' hardware. We conclude with an extensive series of performance tests, which demonstrate that the code can be run efficiently and without modification in serial on small workstations or in parallel using OpenMP compiler directives on large scale, shared memory parallel machines. We analyze the effects of the code optimizations and estimate that they improve its overall performance by more than an order of magnitude over that obtained by many other tree codes. Scaled parallel performance of the gravity and SPH calculations, together the most costly components of most simulations, is nearly linear up to maximum machine sizes available to us (120 processors on an Origin~3000). At similar accuracy, performance of VINE, used in GRAPE-tree mode, is approximately a factor two slower than that of VINE, used in host-only mode. Optimizations of the GRAPE/host communications could improve the speed by as much as a factor of three, but have not yet been implemented in VINE.

We present a Fortran 95 code for simulating the evolution of astrophysical systems using particles to represent the underlying fluid flow. The code is designed to be versatile, flexible and extensible, with modular options that can be selected either at compile time or at run time. We include a number of general purpose modules describing a variety of physical processes commonly required in the astrophysical community. The code can be used as an N-body code to evolve a set of particles in two or three dimensions using either a Leapfrog or Runge-Kutta-Fehlberg integrator, with or without individual timesteps for each particle. Particles may interact gravitationally as $N$-body particles, and all or any subset may also interact hydrodynamically, using the Smoothed Particle Hydrodynamic (SPH) method. Massive point particles (`stars') which may accrete nearby SPH or $N$-body particles may also be included. The default free boundary conditions can be replaced by a module to include periodic boundaries. Cosmological expansion may also be included. An interface with special purpose `GRAPE' hardware may also be selected. If available, forces obtained from the GRAPE coprocessors may be transparently substituted for those obtained from the default tree based calculation. The code may be run without modification on single processors or in parallel using OpenMP compiler directives on large scale, shared memory parallel machines. In comparison to the Gadget-2 code of Springel 2005, the gravitational force calculation is $\approx 3.5 - 4.8$ times faster with VINE when run on 8 Itanium~2 processors in an SGI Altix, while producing nearly identical outcomes in our test problems. We present simulations of several test problems, including a merger simulation of two elliptical galaxies with 800000 particles.

In theories of axion dark matter with large axion decay constant, temperature variations in the CMB are extremely sensitive to perturbations in the initial axion field, allowing one to place a lower bound on the total amount of inflation. The most stringent bound comes from axion strings, which for axion decay constant f=10^17 GeV would currently be observable at a distance of 6 x 10^16 light-years, nearly ten million times as far away as our horizon.

We show that for a force mediated by a vector particle coupled to a conserved U(1) charge, the apparent range and strength can depend on the size and density of the source, and the proximity to other sources. This "chameleon" effect is due to screening from a light charged scalar. Such screening can weaken astrophysical constraints on new gauge bosons. As an example we consider the constraints on chameleonic gauged B-L. We show that although Casimir measurements greatly constrain any B-L force much stronger than gravity with range longer than 0.1 microns, there remains an experimental window for a long range chameleonic B-L force. Such a force could be much stronger than gravity, and long or infinite range in vacuum, but have an effective range near the surface of the earth which is less than a micron.

We review the models and results of simulations of self-gravitating, gaseous protoplanetary disks in binary star systems. These models have been calculated by three different groups with three different computational methods, two particle-based and one grid-based. We show that interactions with the companion star can affect the temperature distribution and structural evolution of disks, and discuss the implications for giant planet formation by gravitational instability as well as by core accretion. Complete consensus has not been reached yet on whether the formation of giant planets is promoted or suppressed by tidal interactions with a companion star. While systems with binary separations of order 100 AU or larger appear to behave more or less as in isolation, systems with smaller separations exhibit an increased or decreased susceptibility to fragmentation, depending on the details of thermodynamics, in particular on the inclusion or absence of artificial viscosity, and on the initial conditions. While code comparisons on identical problems need to be carried out (some of which are already in progress) to decide which computer models are more realistic, it is already clear that relatively close binary systems, with separations of order 60 AU or less, should provide strong constraints on how giant planets form in these systems.

(highly abridged) We define three requirements for accurate simulations that attempt to model circumstellar disks and the formation of collapsed objects (e.g. planets) within them. First, we define a resolution requirement based on the wavelength for neutral stability of self gravitating waves in the disk. For particle based or grid based simulations, this criterion takes the form, respectively, of a minimum number of particles per critical (`Toomre') mass or maximum value of a `Toomre number', $T= \delta x/\lambda_T$, where the wavelength, $\lambda_T$, is the wavelength for neutral stability for waves in disks. We apply our criterion to particle simulations and find that in order to prevent numerically induced fragmentation of the disk, the Toomre mass must be resolved by a minimum of six times the average number of neighbor particles used. Second, we require that particle based simulations with self gravity use a variable gravitational softening. We show that using a fixed gravitational softening length can lead either to artificial suppression or enhancement of structure (including fragmentation) in a given disk, or both in different locations of the same disk. Third, we require that 3D SPH simulations resolve the disk's vertical structure with at least $\sim4$ particle smoothing lengths per scale height at the disk midplane and suggest that a similar criterion applies to grid based simulations. Failure to meet this criterion leads to underestimates in the midplane density of up to 30--50% at resolutions common in the literature.

We show that inclusion of an extremely small quartic coupling constant in the potential for a nearly massless scalar field greatly increases the experimentally allowed region for the mass term and the coupling of the field to matter.

Protoplanetary gas disks are likely to experience gravitational instabilites (GI's) during some phase of their evolution. Density perturbations in an unstable disk grow on a dynamic time scale into spiral arms that produce efficient outward transfer of angular momentum and inward transfer of mass through gravitational torques. In a cool disk with rapid enough cooling, the spiral arms in an unstable disk form self-gravitating clumps. Whether gas giant protoplanets can form by such a disk instability process is the primary question addressed by this review. We discuss the wide range of calculations undertaken by ourselves and others using various numerical techniques, and we report preliminary results from a large multi-code collaboration. Additional topics include -- triggering mechanisms for GI's, disk heating and cooling, orbital survival of dense clumps, interactions of solids with GI-driven waves and shocks, and hybrid scenarios where GI's facilitate core accretion. The review ends with a discussion of how well disk instability and core accretion fare in meeting observational constraints.

We present a supersymmetric model of dark energy from Mass Varying Neutrinos which is stable against radiative corrections to masses and couplings, and free of dynamical instabilities. This is the only such model of dark energy involving fields with significant couplings to any standard model particle. We briefly discuss consequences for neutrino oscillations and solar neutrinos.

We propose that the nebular shocks currently favored as a model to form chondrules and other annealed silicates in the solar nebula originate in the dynamical activity present in the envelope of forming Jovian planets. In contrast to the classic `core accretion model', our 3D hydrodynamic simulations show that this envelope is not a 1D hydrostatic structure but is instead vigorously active and contains densities and temperatures that appear similar in magnitude and spatial extent to those thought to be responsible for the production of chondrules.

We consider a class of theories in which neutrino masses depend significantly on environment, as a result of interactions with the dark sector. Such theories of mass varying neutrinos (MaVaNs) were recently introduced to explain the origin of the cosmological dark energy density and why its magnitude is apparently coincidental with that of neutrino mass splittings. In this Letter we argue that in such theories neutrinos can exhibit different masses in matter and in vacuum, dramatically affecting neutrino oscillations. Both long and short baseline experiments are essential to test for these interactions. As an example of modifications to the standard picture, we consider simple models which may simultaneously account for the LSND anomaly, KamLAND, K2K and studies of solar and atmospheric neutrinos, while providing motivation to continue to search for neutrino oscillations in short baseline experiments such as BooNE.

We show that mass varying neutrinos (MaVaNs) can behave as a negative pressure fluid which could be the origin of the cosmic acceleration. We derive a model independent relation between the neutrino mass and the equation of state parameter of the neutrino dark energy, which is applicable for general theories of mass varying particles. The neutrino mass depends on the local neutrino density and the observed neutrino mass can exceed the cosmological bound on a constant neutrino mass. We discuss microscopic realizations of the MaVaN acceleration scenario, which involve a sterile neutrino. We consider naturalness constraints for mass varying particles, and find that both ev cutoffs and ev mass particles are needed to avoid fine-tuning. These considerations give a (current) mass of order an eV for the sterile neutrino in microscopic realizations, which could be detectable at MiniBooNE. Because the sterile neutrino was much heavier at earlier times, constraints from big bang nucleosynthesis on additional states are not problematic. We consider regions of high neutrino density and find that the most likely place today to find neutrino masses which are significantly different from the neutrino masses in our solar system is in a supernova. The possibility of different neutrino mass in different regions of the galaxy and the local group could be significant for Z-burst models of ultra-high energy cosmic rays. We also consider the cosmology of and the constraints on the ``acceleron'', the scalar field which is responsible for the varying neutrino mass, and briefly discuss neutrino density dependent variations in other constants, such as the fine structure constant.

We use the results of a set of three-dimensional SPH-Treecode simulations which model the formation and early evolution of disk galaxies, including the generation of heavy elements by star formation, to investigate the effects of dust absorption in quasar absorption line systems. Using a simple prescription for the production of dust, we have compared the column density, zinc abundance and optical depth properties of our models to the known properties of Damped Lyman alpha systems. We find that a significant fraction of our model galaxy disks have a higher column density than any observed DLA system. We are also able to show that such parts of the disk tend to be optically thick, implying that any background quasar would be obscured through much of the disk. This would produce the selection effect against the denser absorption systems thought to be present in observations.

We review recent experimental tests of the gravitational inverse-square law and the wide variety of theoretical considerations that suggest the law may break down in experimentally accessible regions.

(abridged) We find that a disk can supply a forming planet with mass at an essentially infinite rate ($\sim1$\mj/25 yr) so that a gap could form very quickly. We show that mass accretion rates faster than $\sim10^{-4}$\mj/yr are not physically reasonable in the limit of either a thin, circumplanetary disk or of a spherical envelope. Planet growth and ultimately survival are therefore limited to the planet's ability to accept additional matter, not by the disk in which it resides. We find that common analytic torque approximations predict values that are a factor $\sim10$ larger than those obtained from the simulations. Accounting for the disk's vertical structure (crudely modeled through a gravitational softening parameter), small shifts in resonance positions due to pressure gradients, to disk self gravity and to inclusion of non-WKB terms in the analysis (Artymowicz 1993) reduce the difference to a factor $\sim3-6$. Torques from the corotation resonances that are positive in sign contribute 20-30% or more of the net torque on the planet. The assumption of linearity underlying theoretical analyses is recovered in the simulations with planets with masses below 0.5\mj, but the assumption that interactions occur only at the resonances is more difficult to support. The detailed shape of the disk's response varies from pattern to pattern, making its true position less clear. We speculate that the finite width allows for overlap and mixing between resonances and may be responsible for the remainder of the differences between torques from theory and simulation, but whether accounting for such overlap in a theory will improve the agreement with the simulations is not clear.

(abridged) We analyze the formation and migration of a proto-Jovian companion in a circumstellar disk in 2d, during the period in which the companion makes its transition from `Type I' to `Type II' migration, using a PPM code. Spiral waves are generated by the gravitational torque of the planet on the disk. Their effects are to cause the planet to migrate inward and the disk to form a deep (low surface density) gap. Until a transition to slower Type II migration, the migration rate of the planet is of order 1 AU/10$^3$ yr, and varies by less than a factor of two with a factor twenty change in planet mass, but depends near linearly on the disk mass. Although the disk is stable to self gravitating perturbations (Toomre $Q>5$ everywhere), migration is faster by a factor of two or more when self gravity is suppressed. Migration is equally sensitive to the disk's mass distribution within 1--2 Hill radii of the planet, as demonstrated by our simulations' sensitivity to the planet's assumed gravitational softening parameter. Rapid migration can continue after gap formation. Gaps are typically several AU in width and display the \mplan$^{2/3}$ proportionality predicted by theory. Beginning from an initially unperturbed 0.05\msun disk, planets of mass $M_{\rm pl}> 0.3$\mj can open a gap deep and wide enough to complete the transition to slower \ttwo migration. Lower mass objects continue to migrate rapidly, eventually impacting the inner boundary of our grid. This transition mass is much larger than that predicted as the `Shiva mass' discussed in Ward and Hahn (2000), making the survival of forming planets even more precarious than they would predict.

The accurate simulation of collapsed objects requires that a huge range of spatial scales be well resolved if the result is not to be contaminated by numerically induced fragmentation. In this context, `insufficient resolution' means comparable to the local instability wavelength in the system. I define a minimum resolution criterion for simulations of disk systems, using the critical wavelength appropriate for disks.

We measure the luminosity profiles of 16 brightest cluster galaxies (BCGs) at $0.4 < z < 0.8$ using high resolution F160W NICMOS and F814W WFPC2 HST imaging. The heterogeneous sample is drawn from a variety of surveys: seven from clusters in the Einstein Medium Sensitivity Survey, five from the Las Campanas Distant Cluster Survey and its northern hemisphere precursor, and the remaining four from traditional optical surveys. We find that the surface brightness profiles of all but three of these BCGs are well described by a standard de Vaucouleurs ($r^{1/4}$) profile out to at least $\sim2r_{e}$ and that the biweight-estimated NICMOS effective radius of our high redshift BCGs ($r_{e} = 8.3\pm 1.4$ kpc for $H_{0} = 80$ km s$^{-1}$ Mpc$^{-1}$, $\Omega_{m} = 0.2, \Omega_\Lambda = 0.0$) is $\sim 2$ times smaller than that measured for a local BCG sample. If high redshift BCGs are in dynamical equilibrium and satisfy the same scaling relations as low redshift ones, this change in size would correspond to a mass growth of a factor of 2 since $z \sim 0.5$. However, the biweight-estimated WFPC2 effective radius of our sample is 18 $\pm $ 5.1 kpc, which is fully consistent with the local sample. While we can rule out mass accretion rates higher than a factor of 2 in our sample, the discrepancy between our NICMOS and WFPC2 results, which after various tests we describe appears to be physical, does not yet allow us to place strong constraints on accretion rates below that level.

The Las Campanas Distant Cluster Survey (LCDCS), which contains over 1000 cluster candidates at z>0.3, is a unique sample with which to probe the evolution of both cluster galaxies and the properties of the cluster population. Programs are now underway to utilize the LCDCS for both purposes. We provide a brief overview of these programs, and also present new results for the LCDCS cluster correlation function. Utilizing well-defined, approximately dispersion-limited subsamples, we measure the angular correlation function for clusters at z~0.5. Spatial correlation lengths are then derived via Limber inversion. We find that the correlation length depends upon mass, as parameterized by the mean cluster separation, in a manner that is consistent with both local observations and CDM predictions for the clustering strength at z=0.5.

We investigate the influence of environment on brightest cluster galaxy (BCG) evolution using a sample of 63 clusters at 0.3 < z < 0.9 drawn primarily from the Las Campanas Distant Cluster Survey and follow-up V, I, and K' photometry. The luminosity evolution of the entire BCG sample is not adequately described by a single evolutionary model. Using the integrated light from the cluster detection as a proxy for cluster Lx and the suggestion by Burke, Collins, & Mann, we set Lx = 2 x 10^{44} ergs/s to be the division between high and low luminosity clusters. At high redshift (z>0.6) BCGs from low-Lx clusters are fainter, on average, than those from high-Lx clusters and are best modeled as having constant luminosity with redshift. The BCGs from high-Lx cluster are best modeled as having a stellar population that formed at large redshift (z_form > 5) and is passively evolving. However, for the entire BCG population, the observed V-I and I-K' colors are well described by a single evolutionary model in which the stellar populations have z_form > 5 and subsequently passively evolve. We conclude that accretion is proportionally more significant for BCGs in lower mass clusters at these redshifts (factor of 2 to 4 increase in mass since z ~ 1 for the low Lx systems) and that the accreted matter is in the form of systems with evolved stellar populations.

Recently, interesting 4-D Lorentz violating models have been proposed, in which all particles have a common maximum velocity $c$, but gravity propagates (in the preferred frame) with a different maximum velocity $c_g \neq c$. We show that the case $c_g < c$ is very tightly constrained by the observation of the highest energy cosmic rays. Assuming a galactic origin for the cosmic rays gives a conservative bound of $c-c_g < 2 \times 10^{-15} c$; if the cosmic rays have an extragalactic origin the bound is orders of magnitude tighter, of order $c-c_g < 2 \times 10^{-19} c$.

(Abridged) We analyze photometry and spectroscopy of a sample of 63 clusters at 0.3<z<0.9 drawn from the Las Campanas Distant Cluster Survey to empirically constrain models of cluster galaxy evolution. Specifically, by combining data on our clusters with those from the literature we parametrize the redshift dependence of 1) M*_I in the observed frame; 2) the V-I color of the E/S0 red sequence in the observed frames; and 3) the I-K' color of the E/S0 red sequence in the observed frame. Using the peak surface brightness of the cluster detection, S, as a proxy for cluster mass, we find no correlation between S and M* or the location of the red envelope in V-I. We suggest that these observations can be explained with a model in which luminous early type galaxies (or more precisely, the progenitors of current day luminous early type galaxies) form the bulk of their stellar populations at high redshift (>~ 5) and in which many of these galaxies, if not all, accrete mass either in the form of evolved stellar populations or gas that causes only a short term episode of star formation at lower redshifts (1.5 < z < 2). Our data are too crude to reach conclusions regarding the evolutionary state of any particular cluster or to investigate whether the morphological evolution of galaxies matches the simple scenario we discuss, but the statistical nature of this study suggests that the observed evolutionary trends are universal in massive clusters.

We present our study of the evolution of brightest cluster galaxies based on data from the Las Campanas Distant Cluster Survey (Gonzalez et al. 2001) and HST optical and IR imaging. We briefly discuss the technique that enabled us to use short (~ 3 min) exposures and a small (1m) telescope to efficiently survey about 130 sq. deg. of sky for high-z galaxy clusters. Follow-up imaging is used to construct a sample of several tens of BCGs at 0.3 < z < 0.9 with which to explore their evolution. In particular, we examine the luminosity and color evolution of these galaxies. We confirm the previous results that 1) BCGs in different environments evolve differently (eg. Burke et al 2000), 2) BCGs, particularly those in low-mass clusters, must be accreting significantly since z ~ 1 (Aragon-Salamanca et al. 1998). From measurments of the colors vs. redshift, we conclude that this accretion must consist of old stellar populations. Using HST NICMOS and WFPC2 images, we find preliminary evidence for evolution in the sizes of BCGs that is consistent with the accretion model presented in standard hierarchical models (Aragon-Salamanca et al. 1998), but puzzling differences between NICMOS and WFPC2 scale lengths precludes a definitive conclusion.

We present the Las Campanas Distant Cluster Survey, which has produced over a thousand galaxy cluster candidates at 0.35 < z < 1.1 (see Gonzalez et al. 2001 for the full catalog). We discuss the technique that enabled us to use short (~ 3 min) exposures and a small (1m) telescope to efficiently survey ~ 130 sq. deg. of sky. Follow-up imaging and spectroscopy using a wide array of telescopes including the Keck and VLT suggest that the bona-fide cluster fraction is ~70%. We construct methods to estimate both the redshift and cluster mass from the survey data themselves and discuss our first result on large-scale structure, the dependence of the cluster correlation length with mean cluster separation at z ~ 0.5 (Gonzalez, Zaritsky, & Wechsler 2001).

We present an optically-selected catalog of 1073 galaxy cluster and group candidates at 0.3<z<1. These candidates are drawn from the Las Campanas Distant Clusters Survey (LCDCS), a drift-scan imaging survey of a 130 square degree strip of the southern sky. To construct this catalog we utilize a novel detection process in which clusters are detected as positive surface brightness fluctuations in the background sky. This approach permits us to find clusters with significantly shallower data than other matched-filter methods that are based upon number counts of resolved galaxies. Selection criteria for the survey are fully automated so that this sample constitutes a well-defined, homogeneous sample that can be used to address issues of cluster evolution and cosmology. Estimated redshifts are derived for the entire sample, and an observed correlation between surface brightness and velocity dispersion (sigma) is used to estimate the limiting velocity dispersion of the survey as a function of redshift. We find a net surface density of 15.5 candidates per square degree at z_{est}>0.3, with a false-detection rate of 30%. At z=0.3 we probe down to the level of poor groups while by z=0.8 we detect only the most massive systems (sigma>1000 km/s). We also present a supplemental catalog of 112 candidates that fail one or more of the automated selection criteria, but appear from visual inspection to be bona fide clusters.

We show that planet formation via both gravitational collapse and core accretion is unlikely to occur in equal mass binary systems with moderate (~ 50 AU) semi-major axes. Internal thermal energy generation in the disks is sufficient to heat the gas everywhere so that spiral structures quickly decay rather than grow or fragment. This same heating will inhibit dust coagulation because the temperatures rise above the vaporization temperatures of many volatile materials. We consider other processes not included in the model and conclude that our temperatures are conservatively estimated (low), i.e. planet formation is less likely in real systems than in the model.

I construct solutions to Einstein's equations in 6 dimensions with bulk cosmological constant and intersecting 4-branes. Solutions exist for a continuous range of 4-brane tension, with long distance gravity localized to a 3+1 dimensional Minkowski intersection, provided that the additional tension of the intersection satisfies one condition.

In my thesis I present a study of the dynamics and observational characteristics of massive circumstellar disks in two dimensions ($r$, $\phi$) using two complimentary hydrodynamic codes: a `Smoothed Particle Hydrodynamic' (SPH) code and a `Piecewise Parabolic Method' (PPM) code. I also study the detection limits available to radial velocity searches for low mass companions to main sequence stars. The thesis is organized as a series of published or submitted papers, connected by introductory and concluding material and is archived `for posterity'. I recommend that readers of this abstract obtain the published versions of each of these papers. I first outline the progress which has been made in the modeling of the structure and origins of the solar system, then in chapter 2 (ApJ 502,p342) I proceed with numerical simulations of circumstellar disks using both hydrodynamic codes assuming a `locally isothermal' equation of state. In chapter 4 (ApJ submitted) I relax the locally isothermal evolution assumption and instead include simple heating and cooling prescriptions for the system. In chapter 6 (ApJ 500,p940) I examine the limits which may be placed upon the detection of planets, brown dwarfs and low mass stellar companions using radial velocity measurements. In chapter 7, I outline several problems which may be profitably addressed by building on this work.

We present a series of 2-d ($r,\phi$) hydrodynamic simulations of marginally self gravitating disks around protostars using an SPH code. We implement simple dynamical heating and we cool each location as a black body, using a photosphere temperature obtained from the local vertical structure. We synthesize SEDs from our simulations and compare them to fiducial SEDs derived from observed systems. These simulations produce less distinct spiral structure than isothermally evolved systems, especially in the inner third of the disk. Pattern are similar further from the star but do not collapse into condensed objects. The photosphere temperature is well fit to a power law in radius with index $q\sim1.1$, which is very steep. Far from the star, internal heating ($PdV$ work and shocks) are not responsible for generating a large fraction of the thermal energy contained in the disk matter. Gravitational torques responsible for such shocks cannot transport mass and angular momentum efficiently in the outer disk. Within $\sim$5--10 AU of the star, rapid break up and reformation of spiral structure causes shocks, which provide sufficient dissipation to power a larger fraction of the near IR energy output. The spatial and size distribution of grains can have marked consequences on the observed near IR SED and can lead to increased emission and variability on $\lesssim 10$ year time scales. When grains are vaporized they do not reform into a size distribution similar to that from which most opacity calculations are based. With rapid grain reformation into the original size distribution, the disk does not emit near infrared photons. With a plausible modification to the opacity, it contributes much more.

We present 2D hydrodynamic simulations of circumstellar disks around protostars using a `Piecewise Parabolic Method' (PPM) code. We include a point mass embedded within the disk and follow the migration of that point mass through the disk. Companions with masses $M_c\ga 0.5M_J$ can open a gap in the disk sufficient to halt rapid migration through the disk. Lower mass companions open gaps, but migration continues because sufficient disk mass remains close to the disk to exert large tidal torques. We find that the torques which dominate the migration of low mass planets originate within a radial region within 1-2 Hill radii of the planet's orbit radius, a distance smaller than the thickness of the disk. We conclude that a very high resolution 3D treatment will be required to adequately describe the planet's migration.

We discuss the first results from two successful simulations of galaxy formation within a cosmological volume. With over 2000 large objects forming in each we have sufficient numbers to reliably produce both galaxy correlation and luminosity functions. We find that the observed galaxy counts are well fitted by these models and that the galaxies display an almost un-evolving correlation function back to a redshift of 3 which closely resembles the featureless observed form and amplitude.

We have simulated the growth of structure in two 100 Mpc boxes for LCDM and SCDM universes. These N-body/SPH simulations include a gaseous component which is able to cool radiatively. A fraction of the gas cools into cold dense objects which we identify as galaxies. In this article we give a preliminary analysis of the clustering behaviour of these galaxies concentrating on the LCDM model. We find a galaxy correlation function which is very close to a power law and which evolves relatively little with redshift. The pairwise dispersions of the galaxies are significantly lower than the dark matter. The LCDM model gives a surprisingly good match to the observational determinations of the galaxy correlation function and pairwise dispersions.

We find a linear instability of non-axisymmetric Rossby waves in a thin non-magnetized Keplerian disk when there is a local maximum in the radial profile of a key function ${\cal L}(r) \equiv {\cal F}(r) S^{2/\Gamma}(r)$, where ${\cal F}^{-1} = \hat {\bf z}\cdot ({\bf \nabla}\times {\bf v}) /\Sigma$ is the potential vorticity, $S = P/\Sigma^\Gamma$ is the entropy, $\Sigma$ is the surface mass density, $P$ is the vertically integrated pressure, and $\Gamma$ is the adiabatic index. We consider in detail the special case where there is a local maximum in the disk entropy profile $S(r)$. This maximum acts to trap the waves in its vicinity if its height to width ratio ${\rm max}(S)/\Delta r$ is larger than a threshold value. The pressure gradient derived from this entropy variation provides the restoring force for the wave growth. We show that the trapped waves act to transport angular momentum outward. A plausible way to produce an entropy variation is when an accretion disk is starting from negligible mass and temperature, therefore negligible entropy. As mass accumulates by either tidal torquing, magnetic torquing, or Roche-lobe overflow, confinement of heat will lead to an entropy maximum at the outer boundary of the disk. Possible nonlinear developments from this instability include the formation of Rossby vortices and the formation of spiral shocks. What remains to be determined from hydrodynamic simulations is whether or not Rossby wave packets (or vortices) ``hold together'' as they propagate radially inward.

We measure the extent of 100 micron galactic emission in two independent galaxy samples using the IRAS 100 micron Sky Survey images and constrain the distribution of dust at large ($\ltsim$30 kpc) radii. The first sample consists of 90 nearby (v < 6000 km/s) galaxies from the RC3 catalog with similar angular sizes and absolute luminosities (5 arcmin $\leq $D$_{25} \leq$ 10 arcmin and $-22.5 \leq $M$_{B} \leq -18$) that are isolated in the 100 micron images. The second sample consists of 24 local galaxies (v < 1500 km/s, 10 arcmin $\leq $D$_{25} \leq$ 30 arcmin). We rescale the 100 micron images of these galaxies using their optical diameters, D$_{25}$, rotate the images using their optical major axis position angle, construct the mean and median image, and rebin the final images into polar coordinates to study the 100 micron emission as a function of radius and azimuthal angle. We find that the 100 micron emission extends at least to radii of 27 kpc (2$\sigma$ detection) for the typical galaxy in the 5 arcmin - 10 arcmin sample and to 21 kpc (2$\sigma$ detection) in the 10 arcmin - 30 arcmin sample (H$_{0} = $75 km/s/Mpc). In both samples, the emission is preferentially elongated along the optical major axis. We fit an exponential to the 100 micron emission along the major axis and measure a scale length of $2.5 \pm$ 0.8 kpc (90% confidence interval). Using a simple model that relates the far-IR emission to the stellar distribution, we examine the range of acceptable dust mass distributions allowed by our data and conclude that the dust is more extended than the starlight.

In this paper we develop an analytical expression relating the sensitivity to a periodic signal to the duration and accuracy of a given set of data. The effects of windowing are explored, and also the sensitivity to periods longer than the total length of observations. As specific examples, we outline the performance of our technique on large amplitude and large eccentricity radial velocity signals recently discussed in the literature and we delineate the region explored by the measurements of 14 single stars made over a twelve year period by Walker et al. (1995). Had any of these stars shown motion like that caused by the exo-planets recently detected, it would have been easily detected. Finally, we suggest an observing strategy for future large radial velocity surveys which, if implemented, will allow coverage of the largest range of parameter space with the smallest amount of observing time per star. We suggest that about 10-15 measurements be made of each star in the first two years of the survey, then 2--3 measurements per year thereafter, provided no (or slow) variation is observed. More frequent observations would of course be indicated if such variations were present.

We present a series of 2-dimensional hydrodynamic simulations of massive disks around protostars. We simulate the same physical problem using both a `Piecewise Parabolic Method' (PPM) code and a `Smoothed Particle Hydrodynamic' (SPH) code, and analyze their differences. The disks studied here range in mass from $0.05 M_*$ to $1.0 M_*$ and in initial minimum Toomre $Q$ value from 1.1 to 3.0. For this problem, the strengths of the codes overlap only in a limited fashion, but similarities exist in their predictions, including spiral arm pattern speeds and morphological features. Our results represent limiting cases (i.e. systems evolved isothermally) rather than true physical systems. Disks become active from the inner regions outward. From the earliest times, their evolution is a strongly dynamic process rather than a smooth progression toward eventual nonlinear behavior. We calculate approximate growth rates for the spiral patterns; the one-armed ($m=1$) spiral arm is not the fastest growing pattern of most disks. In our SPH simulations, disks with initial minimum $Q=1.5$ or lower break up into proto-binary or proto-planetary clumps. However, these simulations cannot follow the physics important for the flow and must be terminated before the system has completely evolved. At their termination, PPM simulations with similar initial conditions show uneven mass distributions within spiral arms, suggesting that clumping behavior might result if they were carried further. Concern that the point-like nature of SPH exaggerates clumping, that our representation of the gravitational potential in PPM is too coarse, and that our physics assumptions are too simple, suggest caution in interpretation of the clumping in both the disk and torus simulations.

We present an analysis of the clustering evolution of dark matter in four cold dark matter (CDM) cosmologies. We use a suite of high resolution, 17-million particle, N-body simulations which sample volumes large enough to give clustering statistics with unprecedented accuracy. We investigate both a flat and an open model with Omega_0=0.3, and two models with Omega=1, one with the standard CDM power spectrum and the other with the same power spectrum as the Omega_0=0.3 models. The amplitude of primordial fluctuations is set so that the models reproduce the observed abundance of rich galaxy clusters by the present day. The mass 2-point correlation function and power spectrum of all the simulations differ significantly from those of the observed galaxy distribution, in both shape and amplitude. Thus, for any of these models to provide an acceptable representation of reality, the distribution of galaxies must be biased relative to the mass in a non-trivial, scale-dependent, fashion. In the Omega=1 models the required bias is always greater than unity, but in the Omega_0=0.3 models an "antibias" is required on scales smaller than \sim 5\hmpc. The mass correlation functions in the simulations are well fit by recently published analytic models. The velocity fields are remarkably similar in all the models, whether they be characterised as bulk flows, 1-particle or pairwise velocity dispersions. This similarity is a direct consequence of our adopted normalisation. The small-scale pairwise velocity dispersion of the dark matter is somewhat larger than recent determinations from galaxy redshift surveys, but the bulk-flows predicted by our models are broadly in agreement with most available data.

We study the topology of the Virgo N-body simulations and compare it to the 1.2-Jy redshift survey of IRAS galaxies by means of the genus statistic. Four high-resolution simulations of variants of the CDM cosmology are considered: a flat standard model (SCDM), a variant of it with more large-scale power (tCDM), and two low density universes, one open (OCDM) and one flat (LCDM). The fully sampled N-body simulations are examined down to strongly nonlinear scales, both with spatially fixed smoothing, and with an adaptive smoothing technique. While the tCDM, LCDM, and OCDM simulations have very similar genus statistics in the regime accessible to fixed smoothing, they can be separated with adaptive smoothing at small mass scales. In order to compare the N-body models with the 1.2-Jy survey, we extract large ensembles of mock catalogues from the simulations. These mock surveys are used to test for systematic effects in the genus analysis and to establish the distribution of errors of the genus curve. We find that a simple multivariate analysis of the genus measurements is compromised both by non-Gaussian distributed errors and by noise that dominates the covariance matrix. We therefore introduce a principal components analysis of the genus curve. With a likelihood ratio test we find that the 1.2-Jy data favours the LCDM, tCDM and OCDM models compared to SCDM. When genus measurements for different smoothing scales are combined, the SCDM model can be excluded at a 99% confidence level, while the other three models fit the 1.2-Jy data well. (abridged)

We investigate the internal structure of clusters of galaxies in high-resolution N-body simulations of 4 different cosmologies. There is a higher proportion of disordered clusters in critical-density than in low-density universes, although the structure of relaxed clusters is very similar in each. Crude measures of substructure, such as the shift in the position of the centre-of-mass as the density threshold is varied, can distinguish the two in a sample of just 20 or so clusters; it is harder to differentiate between clusters in open and flat models with the same density parameter. Most clusters are in a quasi-steady state within the virial radius and are well-described by the density profile of Navarro, Frenk & White (1995).

We report on work done by the Virgo consortium, an international collaboration set up in order to study the formation and evolution of Large Scale Structure using N-body simulations on the latest generation of parallel supercomputers. We show results of 256^3 particle simulations of the formation of clusters in four Dark Matter models with different cosmological parameters. Normalizing the models such that one obtains the correct abundance of rich clusters yields an interesting result: The peculiar velocities of the clusters are almost independent of \Omega, and depend only weakly on \Gamma, the shape parameter of the power spectrum. Thus, it is nearly impossible to distinguish between high and low \Omega models on the basis of the peculiar velocities.

We present the first high redshift (0.3 < z < 1.1) galaxy clusters found by systematically identifying optical low surface brightness fluctuations in the background sky. Using spectra obtained with the Keck telescope and I-band images from the Palomar 1.5m telescope, we conclude that at least eight of the ten candidates examined are high redshift galaxy clusters. The identification of such clusters from low surface brightness fluctuations provides a complementary alternative to classic selection methods based on overdensities of resolved galaxies, and enables us to search efficiently for rich high redshift clusters over large areas of the sky. The detections described here are the first in a survey that covers a total of nearly 140 sq. degrees of the sky and should yield, if these preliminary results are representative, over 300 such clusters.

We report on work in progress by the Virgo consortium, a collaboration set up to carry out large simulations of the formation of galaxies and large-scale structure exploiting the latest generation of parallel supercomputers. We show results of $256^3$ particle N-body simulations of the clustering evolution of dark matter in four cold dark matter models with different cosmological parameters. The high resolution and large volume of these simulations allows us to determine reliably the mass autocorrelation function for pair separations in the range $40\hkpc$ to $20\hmpc$. Comparison of these with the observed galaxy correlation function shows that for any of these models to be viable, the distribution of galaxies must be biased relative to the distribution of mass in a non-trivial, scale-dependent fashion. In particular, low $\Omega_0$ models require the galaxies to be more ``weakly'' clustered than the mass at small and intermediate pair separations. Simulations which include the evolution of gas show that cold gas knots form with approximately the abundance expected on theoretical grounds, although a few excessively massive objects grow near the centres of rich clusters. The locations where these cold gas knots form are, in general, biased relative to the distribution of mass in a scale-dependent way. Some of these biases have the required sign but they are, for the most part, weaker than is necessary for agreement with observations. The antibias present in our low $\Omega_0$ N-body/SPH simulation appears to be related to the merging and disruption of galaxies in rich clusters.