2655191
A proposed origin for chondrule-forming shocks in the solar nebula
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.
5533361
The Range of Masses and Periods Explored by Radial Velocity Searches for Planetary Companions
Radial velocity measurements have proven a powerful tool for finding planets in short-period orbits around other stars. 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 and the sensitivity to periods longer than the total length of observations are explored. We show that current observations are not yet long or accurate enough to make unambiguous detection of planets with the same mass and period as Jupiter. However, if measurements are continued at the current levels of accuracy (5 m s-1) for a decade, then planets of Jovian mass and brown dwarfs will either be detected or ruled out for orbits with periods less than ~15 yr. 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 12 yr period by Walker et al. Had any of these stars shown motion like that caused by the exoplanets recently detected, it would have been easily detected. The data set interesting limits on the presence of brown dwarfs at orbital radii of 5-10 AU. The most significant features in the Walker et al. data are apparent long-term velocity trends in 36 UMa and β Vir, consistent with super planets of mass of 2 MJ in a 10 yr period, or 20-30 MJ in a 50 yr period. If the data are free of long-term systematic errors, the probability of just one of the 14 stars showing this signal by chance is about 15%. Finally, we suggest an observing strategy for future large radial velocity surveys that, 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 2 yr of the survey, then 2-3 measurements yr-1 thereafter, provided no (or slow) variation is observed. More frequent observations would of course be indicated if such variations were present.
17175349
Numerical requirements for simulations of self-gravitating and non-self-gravitating discs
We define three requirements for accurate simulations that attempt to model circumstellar discs 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 disc, where a Jeans analysis does not apply. 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=δx/λT, where the wavelength, λT, is the wavelength for neutral stability for waves in discs. The requirements are analogues of the conditions for cloud collapse simulations as discussed in Bate & Burkert and Truelove et al., where the required minimum resolution was shown to be twice the number of neighbours per Jeans mass or four–five times the local Jeans wavelength, λJ, for particle or grid simulations, respectively. 
 
 
 
We apply our criterion to particle simulations of disc evolution and find that in order to prevent numerically induced fragmentation of the disc, the Toomre mass must be resolved by a minimum of six times the average number of neighbour particles used. We investigate the origin of the apparent discrepancy between the number of particles required by the cloud and disc fragmentation criteria and find that it is due largely to ambiguities in the definition of the Jeans mass, as used by different authors. We reconcile the various definitions, and when an identical definition of the Jeans mass is used, the condition that J≲ 1/4 in the Truelove condition is equivalent to requiring about 10–12 times the average number of neighbour particles per Jeans mass in a smoothed particle hydrodynamics (SPH) simulation, reducing the difference between simulations of discs and clouds to about two. While the numbers of particles per critical mass are similar for both the Jeans and Toomre formalisms, the Toomre requirement is more restrictive than the Jeans requirement when the local value of the Toomre stability parameter Q falls below about one half. 
 
 
 
Second, we require that particle-based simulations with self-gravity use a variable gravitational softening, in order to avoid inducing fragmentation by an inappropriate choice of softening length. We show that using a fixed gravitational softening length for all particles can lead either to artificial suppression or enhancement of structure (including fragmentation) in a given disc, or both in different locations of the same disc, depending on the value chosen for the softening length. Unphysical behaviour can occur whether or not the system is properly resolved by the new Toomre criterion. 
 
 
 
Third, we require that three-dimensional SPH simulations resolve the vertical structure with at least ∼4 particle smoothing lengths per scaleheight at the disc mid-plane, a value which implies a substantially larger number per vertical column because the disc itself extends over many scaleheights. We suggest that a similar criterion applies to grid-based simulations. We demonstrate that failure to meet this criterion leads to underestimates in the mid-plane density of up to 30–50 per cent at resolutions common in the literature. As a direct consequence, gas pressures will be dramatically underestimated and simulations of self-gravitating systems may artificially and erroneously inflate the likelihood of fragmentation. We outline an additional condition on the vertical resolution in simulations that include radiative transfer in order to ensure a correct description of the cooling, specifically that the temperature structure near the disc photosphere must be well resolved. As an example, we demonstrate that for an isentropic vertical structure, the criterion translates into resolution comparable to H/20 near the disc photosphere, to avoid serious errors in transfer rates of thermal energy in and out of the disc. 
 
 
 
Finally, we discuss results in the literature that purport to form collapsed objects and conclude that many are likely to have violated one or more of our criteria, and have therefore made incorrect conclusions regarding the likelihood for fragmentation and planet formation.
16491613
VINE—A NUMERICAL CODE FOR SIMULATING ASTROPHYSICAL SYSTEMS USING PARTICLES. II. IMPLEMENTATION AND PERFORMANCE CHARACTERISTICS
We continue our presentation of VINE. In this paper, 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 the modifications to the code necessary to obtain forces efficiently from special purpose ‘GRAPE’ hardware, the interfaces required to allow transparent substitution of those forces in the code instead of those obtained from the tree, and the modifications necessary to use both tree and GRAPE together as a fused GRAPE/tree combination. 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 the 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 at least 120 processors on moderate sized test problems using the Origin 3000 architecture, and to the maximum machine sizes available to us on several other architectures. 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. Further 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. Finally, we find that although parallel performance on small problems may reach a plateau beyond which more processors bring no additional speedup, performance never decreases, a factor important for running large simulations on many processors with individual time steps, where only a small fraction of the total particles require updates at any given moment. Subject headings: methods: numerical — methods: N-body simulations
282379
Planet Formation is Unlikely in Equal-Mass Binary Systems with a ~ 50 AU
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) semimajor 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.
12208225
VINE—A NUMERICAL CODE FOR SIMULATING ASTROPHYSICAL SYSTEMS USING PARTICLES. I. DESCRIPTION OF THE PHYSICS AND THE NUMERICAL METHODS
We present a numerical code for simulating the evolution of astrophysical systems using particles to represent the underlying fluid flow. The code is written in Fortran 95 and is designed to be versatile, flexible and extensible, with modular options that can be selected either at the time the code is compiled or at run time through a text input file. We include a number of general purpose modules describing a variety of physical processes commonly required in the astrophysical community and we expect that the effort required to integrate additional or alternate modules into the code will small. In its simplest form the code can evolve the dynamical trajectories of a set of particles in two or three dimensions using a module which implements either a Leapfrog or Runge-Kutta-Fehlberg integrator, selected by the user at compile time. The user may choose to allow the integrator to evolve the system using individual timesteps for each particle or with a single, global time step for all. Particles may interact gravitationally as N-body particles, and all or any subset may also interact hydrodynamically, using the Smoothed Particle Hydrodynamic (SPH) method by selecting the SPH module. A third particle species can be included with a module to model massive point particles which may accrete nearby SPH or N-body particles. Such particles may be used to model, e.g., stars in a molecular cloud. Free boundary conditions are implemented by default, and a module may be selected to include periodic boundary conditions. We use a binary ‘Press’ tree to organize particles for rapid access in gravity and SPH calculations. Modules implementing an interface with special purpose ‘GRAPE’ hardware may also be selected to accelerate the gravity calculations. If available, forces obtained from the GRAPE coprocessors may be transparently substituted for those obtained from the tree, or both tree and GRAPE may be used as a combination GRAPE/tree code. The code may be run without modification on single processors or in parallel using OpenMP compiler directives on large scale, shared memory parallel machines. We present simulations of several test problems, including a merger simulation of two elliptical galaxies with 800000 particles. In comparison to the Gadget-2 code of Springel (2005), the gravitational force calculation, which is the most costly part of any simulation including self-gravity, is ∼ 4.6 − 4.9 times faster with VINE when tested on different snapshots of the elliptical galaxy merger simulation when run on an Itanium 2 processor in an SGI Altix. A full simulation of the same setup with 8 processors is a factor of 2.91 faster with VINE. The code is available to the public under the terms of the Gnu General Public License. Subject headings: methods: numerical — methods: N-body simulations — galaxies: interactions
15538906
Dynamical Electroweak Superconductivity from a Composite Little Higgs
I describe, from the bottom up, a sequence of natural effective field theories. Below a TeV we have the minimal standard model with a light Higgs, and an extra neutral scalar. In the 1-10 TeV region these scalars are part of a multiplet of pseudo Nambu-Goldstone Bosons. Interactions with additional TeV mass scalars, gauge bosons, and vector-like charge 2/3 quarks stabilize the Higgs mass squared parameter without finetuning. Electroweak superconductivity may be determined in this effective theory as a UV insensitive vacuum alignment problem. Above the 10 TeV scale we have strongly coupled new gauge interactions.
13891269
Gaugino mediated supersymmetry breaking
We consider supersymmetric theories where the standard-model quark and lepton fields are localized on a `3-brane' in extra dimensions, while the gauge and Higgs fields propagate in the bulk. If supersymmetry is broken on another 3-brane, supersymmetry breaking is communicated to gauge and Higgs fields by direct higher-dimension interactions, and to quark and lepton fields via standard-model loops. We show that this gives rise to a realistic and predictive model for supersymmetry breaking. The size of the extra dimensions is required to be of order 10-100 times larger than fundamental scale (e.g. the string scale). The spectrum is similar to (but distinguishable from) the predictions of `no-scale' models. Flavor-changing neutral currents are naturally suppressed. The μ term can be generated by the Giudice-Masiero mechanism. The supersymmetric CP problem is naturally solved if CP violation occurs only on the observable sector 3-brane. These are the simplest models in the literature that solve all supersymmetric naturalness problems.
45469605
Supersymmetric theories of neutrino dark energy
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.
119283388
CP Violation in Pseudo-Dirac Fermion Oscillations
Supersymmetric theories with a $U(1)_R$ symmetry have Dirac gauginos, solve the supersymmetric flavor and CP problems, and have distinctive collider signatures. However when supergravity is included, the $U(1)_R$ must be broken, adding small Majorana mass terms which split the mass of the two components of the Dirac gaugino and lead to oscillations between $U(1)_R$ charge eigenstates. We present a general study of fermion-antifermion oscillations in this system, including the effects of decays and CP violation. We consider the effects of such oscillations in the case where the two $U(1)_R$ charge eigenstates can decay into the same final state, and show that $O(1)$ CP violation is allowed. In the case of decays into final states containing leptons such CP violation can be observed as a same sign dilepton asymmetry.
18977934
A viable model of dynamical supersymmetry breaking in the hidden sector
Abstract We propose a simple and natural model of dynamical supersymmetry breaking, which could be used as a mechanism for spontaneous supersymmetry breaking in a gravitationally coupled hidden sector. The gaugino masses in the visible sector are naturally of the same size as the squark, slepton, and gravitino masses.