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.
64863479
Dynamics of Core Accretion
We perform 3-dimensional hydrodynamic simulations of gas flowing around a plan- etary core of mass Mpl=10Membedded in a near Keplerian background flow, using a modified shearing box approximation. We assume an ideal gas behavior following an equation of state with a fixed ratio of the specific heats, γ = 1.42, consistent with the conditions of a moderate temperature background disk with solar composition. No radiative heating or cooling is included in the models. We employ a nested grid hydrodynamic code implementing the 'Piecewise Parabolic Method' with as many as six fixed 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 en- velope 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. Each of these factors influence the magnitude and character of hydrodynamic feedback of the small scale flow on the background, and we conclude that accurate modeling of such feedback is critical to a complete understanding of the core accretion process. 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. No net mass accretion into the envelope is observed over the course of the simulation and none is expected, due to our neglect of cooling. Except in cases of very rapid cooling however, as defined by locally isothermal or isentropic treatments, any cooling that does affect the envelope material will have limited consequences for the dynamics, since the flow quickly carries cooled material out of the core's environment entirely. 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. The dynamical activity contrasts with the largely static behavior typically assumed within the framework of the core accretion model for Jovian planet formation. 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. Shorter durations are not observable in our work due to the limits of our resolution. Peak temperatures in these events range from T � 1000 K to as high as T � 3 4000 K, with densities ρ � 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 therefore propose that dynamical activity in the Jovian planet formation environment could be responsible for the production of chondrules and other annealed silicates in the solar nebula.
8914218
Rossby Wave Instability of Keplerian Accretion Disks
We find a linear instability of nonaxisymmetric Rossby waves in a thin nonmagnetized Keplerian disk when there is a local maximum in the radial profile of a key function (r)≡(r)S2/Γ(r), where −1=(∇×v)/Σ is the potential vorticity, S=P/ΣΓ is the entropy, Σ is the surface mass density, P is the vertically integrated pressure, and Γ 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 max(S)/Δ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.
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.
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
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
14533841
New multi-scale supersymmetric models with flavor changing neutral current suppression
Abstract We discuss the phenomenology of a class of supersymmetric models in which some of the quark and lepton superfields are an integral part of a dynamical supersymmetry breaking sector. The corresponding squarks and sleptons are much heavier than any other superpartners, and could naturally have masses as high as ∼ 40 TeV. We discuss a general set of conditions for acceptable flavor-changing neutral currents and natural electroweak symmetry breaking, and identify two particularly interesting new classes of theories. We discuss how phenomenological signatures of such multi-scale models at the CERN LEP II and Fermilab Tevatron colliders could significantly differ from previously considered scenarios. In particular, we give experimental signals which could be present if the left-handed selectron is much lighter than the right-handed one.
31227927
Flavor, compositeness, and dynamical breaking of supersymmetry
We should be taking advantage of recent gains in our nonperturbative understanding of supersymmetric gauge theories to find the “standard” model of of dynamical supersymmetry breaking, and possibly of flavor as well. As an illustration of the possibilities for understanding the flavor hierarchy, I describe a realistic, renormalizable, supersymmetric model with a compositeness scale of ∼ 1–3 TeV for the top quark, the left handed bottom quark, and the up-type Higgs. The top-Higgs Yukawa coupling is a dynamically generated strong interaction effect, and is naturally large, while the other Yukawa couplings are suppressed.
117900789
CP Violating Contribution to $\Delta \Gamma$ in the $B_s$ System from Mixing with a Hidden Pseudoscalar
Recent evidence for a CP violating asymmetry in the semileptonic decays of Bs mesons cannot be accommodated within the Standard Model. Such an asymmetry can be explained by new physics contributions to B = 2 components of either the mass matrix or the decay matrix. We show that mixing with a hidden pseudoscalar meson with a mass around 5 GeV can result in a new CP violating contribution to the mixing and can resolve several anomalies in this system including the width dierence, the average width and the charge asymmetry. We also discuss the eects of the hidden meson on other b physics observables, and present viable decay modes for the hidden meson. We make predictions for new decay channels of B hadrons, which can be tested at the Tevatron, the LHC and B-factories.
121782542
An effective up quark mass from new light particles
Abstract I consider a model where an effective current mass for the up quark is provided by the interactions of a new neutral light colored scalar and a charge 1 3 colorless fermion. The model has no strong CP problem. A new confining U (1) prevents the occurence of physical states with fractional charges. There are new light (a few GeV) states which rapidly decay into hadrons. Contributions to R and to the Z hadronic width are within experimental bounds.
17324735
Solar and atmospheric neutrino oscillations from bilinear R-parity violation
We discuss general predictions for neutrino masses and mixing angles from R parity violation in the Minimal Supersymmetric Standard Model. If the soft supersymmetry breaking terms are flavor blind at short distance, then the leptonic analogue of the CKM matrix depends on only two real parameters, which are completely determined by fits to solar and atmospheric neutrino oscillations. Either the small angle MSW, large angle MSW, or ``just-so'' solutions to the solar neutrino problem are allowed, although the large angle MSW solution requires substantial fine-tuning. The latter two cases require significant νμ→νe oscillations of atmospheric neutrinos. We present a model which could explain bilinear R parity violation as a consequence of spontaneous symmetry violation by a dynamical supersymmetry breaking sector. The decay length and branching ratios of the LSP are estimated.