BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Jeremy Goodman

Jeremy Goodman appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-phastro-ph.EPastro-ph.HEphysics.flu-dynastro-ph.COastro-ph.IMphysics.plasm-phastro-ph.SRastro-ph.GAgr-qchep-thphysics.comp-phquant-ph

Catalog footprint

What is connected

24works
13topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

24 published item(s)

preprint2021arXiv

Evolution of Retrograde Orbiters in an AGN Disk

AGN disks have been proposed as promising locations for the mergers of stellar mass black hole binaries (BBHs). Much recent work has been done on this merger channel, but the majority focuses on stellar mass black holes (BHs) orbiting in the prograde direction. Little work has been done to examine the impact of retrograde orbiters (ROs) on the formation and mergers of BBHs in AGN disks. Quantifying the retrograde contribution is important, since roughly half of all orbiters should initially be on retrograde orbits when the disk forms. We perform an analytic calculation of the evolution of ROs in an AGN disk. Because this evolution could cause the orbits of ROs to cross those of prograde BBHs, we derive the collision rate between a given RO and a given BBH orbiting in the prograde direction. In the examples given here, ROs in the inner region of the disk experience a rapid decrease in the semimajor axis of their orbits while also becoming highly eccentric in less than a million years. This rapid orbital evolution could lead to extreme mass ratio inspirals detectable by the Laser Interferometer Space Antenna. The collision rates of our example ROs with prograde BBHs in the migration trap depend strongly on the volume of the inner radiation-pressure-dominated region which depends on the mass of the supermassive black hole (SMBH). Rates are lowest for larger mass SMBHs, which dominate the AGN merger channel, suggesting that merger rates for this channel may not be significantly altered by ROs.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Numerical simulations of the Princeton magneto-rotational instability experiment with conducting axial boundaries

We investigate numerically the Princeton magneto-rotational instability (MRI) experiment and the effect of conducting axial boundaries or endcaps. MRI is identified and found to reach a much higher saturation than for insulating endcaps. This is probably due to stronger driving of the base flow by the magnetically rather than viscously coupled boundaries. Although the computations are necessarily limited to lower Reynolds numbers ($\Re$) than their experimental counterparts, it appears that the saturation level becomes independent of $\Re$ when $\Re$ is sufficiently large, whereas it has been found previously to decrease roughly as $\Re^{-1/4}$ with insulating endcaps. The much higher saturation levels will allow for the first positive detection of MRI beyond its theoretical and numerical predictions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Magneto-thermal Disk Wind from Protoplanetary Disks

Global evolution and dispersal of protoplanetary disks (PPDs) is governed by disk angular momentum transport and mass-loss processes. Recent numerical studies suggest that angular momentum transport in the inner region of PPDs is largely driven by magnetized disk wind, yet the wind mass-loss rate remains unconstrained. On the other hand, disk mass loss has conventionally been attributed to photoevaporation, where external heating on the disk surface drives a thermal wind. We unify the two scenarios by developing a 1D model of magnetized disk winds with a simple treatment of thermodynamics as a proxy for external heating. The wind properties largely depend on 1) the magnetic field strength at the wind base, characterized by the poloidal Alfvén speed $v_{Ap}$, 2) the sound speed $c_s$ near the wind base, and 3) how rapidly poloidal field lines diverge (achieve $R^{-2}$ scaling). When $v_{Ap}\gg c_s$, corotation is enforced near the wind base, resulting in centrifugal acceleration. Otherwise, the wind is accelerated mainly by the pressure of the toroidal magnetic field. In both cases, the dominant role played by magnetic forces likely yields wind outflow rates that well exceed purely hydrodynamical mechanisms. For typical PPD accretion-rate and wind-launching conditions, we expect $v_{Ap}$ to be comparable to $c_s$ at the wind base. The resulting wind is heavily loaded, with total wind mass loss rate likely reaching a considerable fraction of wind-driven accretion rate. Implications for modeling global disk evolution and planet formation are also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On obliquely magnetized and differentially rotating stars

We investigate the interaction of differential rotation and a misaligned magnetic field. The incompressible magnetohydrodynamic equations are solved numerically for a free-decay problem. In the kinematic limit, differential rotation annihilates the non-axisymmetric field on a timescale proportional to the cube root of magnetic Reynolds number ($Rm$), as predicted by Rädler. Nonlinearly, the outcome depends upon the initial energy in the non-axisymmetric part of the field. Sufficiently weak fields approach axisymmetry as in the kinematic limit; some differential rotation survives across magnetic surfaces, at least on intermediate timescales. Stronger fields enforce uniform rotation and remain non-axisymmetric. The initial field strength that divides these two regimes does not follow the scaling $Rm^{-1/3}$ predicted by quasi-kinematic arguments, perhaps because our $Rm$ is never sufficiently large or because of reconnection. We discuss the possible relevance of these results to tidal synchronization and tidal heating of close binary stars, particularly double white dwarfs.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A method for the microlensed flux variance of QSOs

A fast and practical method is described for calculating the microlensed flux variance of an arbitrary source by uncorrelated stars. The required inputs are the mean convergence and shear due to the smoothed potential of the lensing galaxy, the stellar mass function, and the absolute square of the Fourier transform of the surface brightness in the source plane. The mathematical approach follows previous authors but has been generalized, streamlined, and implemented in publicly available code. Examples of its application are given for Dexter and Agol's inhomogeneous-disk models as well as the usual gaussian sources. Since the quantity calculated is a second moment of the magnification, it is only logarithmically sensitive to the sizes of very compact sources. However, for the inferred sizes of actual QSOs, it has some discriminatory power and may lend itself to simple statistical tests. At the very least, it should be useful for testing the convergence of microlensing simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A New Way to Conserve Total Energy for Eulerian Hydrodynamic Simulations with Self-Gravity

We propose a new method to conserve the total energy to round-off error in grid-based codes for hydrodynamic simulations with self-gravity. A formula for the energy flux due to the work done by the the self-gravitational force is given, so the change in total energy can be written in conservative form. Numerical experiments with the code Athena show that the total energy is indeed conserved with our new algorithm and the new algorithm is second order accurate. We have performed a set of tests that show the numerical errors in the traditional, non-conservative algorithm can affect the dynamics of the system. The new algorithm only requires one extra solution of the Poisson equation, as compared to the traditional algorithm which includes self-gravity as a source term. If the Poisson solver takes a negligible fraction of the total simulation time, such as when FFTs are used, the new algorithm is almost as efficient as the original method. This new algorithm is useful in Eulerian hydrodynamic simulations with self-gravity, especially when results are sensitive to small energy errors, as for radiation pressure dominated flow.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Observation of a Free-Shercliff-Layer Instability in Cylindrical Geometry

We report on observations of a free-Shercliff-layer instability in a Taylor-Couette experiment using a liquid metal over a wide range of Reynolds numbers, $Re\sim 10^3-10^6$. The free Shercliff layer is formed by imposing a sufficiently strong axial magnetic field across a pair of differentially rotating axial endcap rings. This layer is destabilized by a hydrodynamic Kelvin-Helmholtz-type instability, characterized by velocity fluctuations in the $r-θ$ plane. The instability appears with an Elsasser number above unity, and saturates with an azimuthal mode number $m$ which increases with the Elsasser number. Measurements of the structure agree well with 2D global linear mode analyses and 3D global nonlinear simulations. These observations have implications for a range of rotating MHD systems in which similar shear layers may be produced.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Ruling Out Bosonic Repulsive Dark Matter in Thermal Equilibrium

Self-interacting dark matter (SIDM), especially bosonic, has been considered a promising candidate to replace cold dark matter (CDM) as it resolves some of the problems associated with CDM. Here, we rule out the possibility that dark matter is a repulsive boson in thermal equilibrium. We develop the model first proposed by Goodman (2000) and derive the equation of state at finite temperature. Isothermal spherical halo models indicate a Bose-Einstein condensed core surrounded by a non-degenerate envelope, with an abrupt density drop marking the boundary between the two phases. Comparing this feature with observed rotation curves constrains the interaction strength of our model's DM particle, and Bullet Cluster measurements constrain the scattering cross section. Both ultimately can be cast as constraints on the particle's mass. We find these two constraints cannot be satisfied simultaneously in any realistic halo model---and hence dark matter cannot be a repulsive boson in thermal equilibrium. It is still left open that DM may be a repulsive boson provided it is not in thermal equilibrium; this requires that the mass of the particle be significantly less than a millivolt.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Stability of quasi-Keplerian shear flow in a laboratory experiment

Context: Subcritical transition to turbulence has been proposed as a source of turbulent viscosity required for the associated angular momentum transport for fast accretion in Keplerian disks. Previously cited laboratory experiments in supporting this hypothesis were performed either in a different type of flow than Keplerian or without quantitative measurements of angular momentum transport and mean flow profile, and all of them appear to suffer from Ekman effects, secondary flows induced by nonoptimal axial boundary conditions. Such Ekman effects are expected to be absent from astronomical disks, which probably have stress-free vertical boundaries unless strongly magnetized. Aims: To quantify angular momentum transport due to subcritical hydrodynamic turbulence, if exists, in a quasi-Keplerian flow with minimized Ekman effects. Methods: We perform a local measurement of the azimuthal--radial component of the Reynolds stress tensor in a novel laboratory apparatus where Ekman effects are minimized by flexible control of axial boundary conditions. Results: We find significant Ekman effects on angular momentum transport due to nonoptimal axial boundary conditions in quasi-Keplerian flows. With the optimal control of Ekman effects, no statistically meaningful angular momentum transport is detected in such flows at Reynolds number up to two millions. Conclusions: Either a subcritical transition does not occur, or, if a subcritical transition does occur, the associated radial transport of angular momentum in optimized quasi-Keplerian laboratory flows is too small to directly support the hypothesis that subcritical hydrodynamic turbulence is responsible for accretion in astrophysical disks. Possible limitations in applying laboratory results to astrophysical disks due to experimental geometry are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The role of boundaries in the MagnetoRotational Instability

In this paper, we investigate numerically the flow of an electrically conducting fluid in a cylindrical Taylor-Couette flow when an axial magnetic field is applied. To minimize Ekman recirculation due to vertical no-slip boundaries, two independently rotating rings are used at the vertical endcaps. This configuration reproduces setup used in laboratory experiments aiming to observe the MagnetoRotational Instability (MRI). Our 3D global simulations show that the nature of the bifurcation, the non-linear saturation, and the structure of axisymmetric MRI modes are significantly affected by the presence of boundaries. In addition, large scale non-axisymmetric modes are obtained when the applied field is sufficiently strong. We show that these modes are related to Kelvin-Helmoltz destabilization of a free Shercliff shear layer created by the combined action of the applied field and the rotating rings at the endcaps. Finally, we compare our numerical simulations to recent experimental results obtained in the Princeton MRI experiment.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Chance and Chandra (and repulsive dark matter)

A few examples are given of Chandra's work on statistical and stochastic problems that relate to open questions in astrophysics, in particular his theory of dynamical relaxation in systems with inverse-square interparticle forces. The roles of chaos and integrability in this theory require clarification, especially for systems having a dominant central mass. After this prelude, a hypothetical form of repulsive bosonic dark matter is discussed. The repulsion leads to nontrivial thermodynamic behavior, including superfluidity, and would tend to suppress dynamical friction, greatly reducing the drag exerted on rotating galactic bars. However, this form of dark matter can probably be ruled out, at least for parameters that allow halos to reach thermal equilibria within a Hubble time. One combination of the particle mass and interparticle repulsion determines the minimum core radius of dark halos. Bounds on dark-matter collisionality inferred from the Bullet Cluster constrain a second combination. It is possible to satisfy both constraints only for parameters that predict unacceptable rotation curves outside the halo core.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Instabilities in magnetized spherical Couette flow

We report 3D numerical simulations of the flow of an electrically conducting fluid in a spherical shell when a magnetic field is applied. Different spherical Couette configurations are investigated, by varying the rotation ratio between the inner and the outer sphere, the geometry of the imposed field, and the magnetic boundary conditions on the inner sphere. Either a Stewartson layer or a Shercliff layer, accompanied by a radial jet, can be generated depending on the rotation speeds and the magnetic field strength, and various non-axisymmetric destabilizations of the flow are observed. We show that instabilities arising from the presence of boundaries present striking similarities with the magnetorotational instability (MRI). To this end, we compare our numerical results to experimental observations of the Maryland experiment, who claimed to observe MRI in a similar setup.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Circulation and Dissipation on Hot Jupiters

Many global circulation models predict supersonic zonal winds and large vertical shears in the atmospheres of short-period jovian exoplanets. Using linear analysis and nonlinear local simulations, we investigate hydrodynamic dissipation mechanisms to balance the thermal acceleration of these winds. The adiabatic Richardson criterion remains a good guide to linear stability, although thermal diffusion allows some modes to violate it at very long wavelengths and very low growth rates. Nonlinearly, wind speeds saturate at Mach numbers $\approx 2$ and Richardson numbers $\lesssim 1/4$ for a broad range of plausible diffusivities and forcing strengths. Turbulence and vertical mixing, though accompanied by weak shocks, dominate the dissipation, which appears to be the outcome of a recurrent Kelvin-Helmholtz instability. An explicit shear viscosity, as well as thermal diffusivity, is added to ZEUS to capture dissipation outside of shocks. The wind speed is not monotonic nor single valued for shear viscosities larger than about $10^{-3}$ of the sound speed times the pressure scale height. Coarsening the numerical resolution can also increase the speed. Hence global simulations that are incapable of representing vertical turbulence and shocks, either because of reduced physics or because of limited resolution, may overestimate wind speeds. We recommend that such simulations include artificial dissipation terms to control the Mach and Richardson numbers and to capture mechanical dissipation as heat.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Star Formation in Quasar Disk

Using a version of the ZEUS code, we carry out two-dimensional simulations of self-gravitating shearing sheets, with application to QSO accretion disks at a few thousand Schwarzschild radii, corresponding to a few hundredths of a parsec for a 10^8 solar-mass black hole. Radiation pressure and optically thick radiative cooling are implemented via vertical averages. We determine dimensionless versions of the maximum surface density, accretion rate, and effective viscosity that can be sustained by density-wave turbulence without fragmentation. Where fragments do form, we study the final masses that result. The maximum Shakura-Sunyaev viscosity parameter is approximately 0.4. Fragmentation occurs when the cooling time is less than about twice the shearing time, as found by Gammie and others, but can also occur at very long cooling times in sheets that are strongly radiation-pressure dominated. For accretion at the Eddington rate onto a 10^8 solar-mass black hole, fragmentation occurs beyond four thousand Schwarzschild radii, r_s. Near this radius, initial fragment masses are several hundred suns, consistent with estimates from linear stability; final masses after merging increase with the size of the sheet, reaching several thousand suns in our largest simulations. With increasing black-hole mass at a fixed Eddington ratio, self-gravity prevails to smaller multiples of r_s, where radiation pressure is more important and the cooling time is longer compared to the dynamical time; nevertheless, fragmentation can occur and produces larger initial fragment masses. We also find energy conservation is likely to be a challenge for all eulerian codes in self-gravitating regimes where radiation pressure dominates.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

On the Saturation of the Magnetorotational Instability via Parasitic Modes

We investigate the stability of incompressible, exact, non-ideal magnetorotational (MRI) modes against parasitic instabilities. Both Kelvin-Helmholtz and tearing-mode parasitic instabilities may occur in the dissipative regimes accessible to current numerical simulations. We suppose that a primary MRI mode saturates at an amplitude such that its fastest parasite has a growth rate comparable to its own. The predicted alpha parameter then depends critically on whether the fastest primary and parasitic modes fit within the computational domain and whether non-axisymmetric parasitic modes are allowed. Hence even simulations that resolve viscous and resistive scales may not saturate properly unless the numerical domain is large enough to allow the free evolution of both MRI and parasitic modes. To minimally satisfy these requirements in simulations with vertical background fields, the vertical extent of the domain should accommodate the fastest growing MRI mode while the radial and azimuthal extents must be twice as large. The fastest parasites have horizontal wavelengths roughly twice as long as the vertical wavelength of the primary.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Dynamical Tides in Rotating Planets and Stars

Tidal dissipation may be important for the internal evolution as well as the orbits of short-period massive planets--hot Jupiters. We revisit a mechanism proposed by Ogilvie and Lin for tidal forcing of inertial waves, which are short-wavelength, low-frequency disturbances restored primarily by Coriolis rather than buoyancy forces. This mechanism is of particular interest for hot Jupiters because it relies upon a rocky core, and because these bodies are otherwise largely convective. Compared to waves excited at the base of the stratified, externally heated atmosphere, waves excited at the core are more likely to deposit heat in the convective region and thereby affect the planetary radius. However, Ogilvie and Lin's results were numerical, and the manner of the wave excitation was not clear. Using WKB methods, we demonstrate the production of short waves by scattering of the equilibrium tide off the core at critical latitudes. The tidal dissipation rate associated with these waves scales as the fifth power of the core radius, and the implied tidal $Q$ is of order ten million for nominal values of the planet's mass, radius, orbital period, and core size. We comment upon an alternative proposal by Wu for exciting inertial waves in an unstratified fluid body by means of compressibility rather than a core. We also find that even a core of rock is unlikely to be rigid. But Ogilvie and Lin's mechanism should still operate if the core is substantially denser than its immediate surroundings.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Statistical Description of a Magnetized Corona above a Turbulent Accretion Disk

We present a physics-based statistical theory of a force-free magnetic field in the corona above a turbulent accretion disk. The field is represented by a statistical ensemble of loops tied to the disk. Each loop evolves under several physical processes: Keplerian shear, turbulent random walk of the disk footpoints, and reconnection with other loops. To build a statistical description, we introduce the distribution function of loops over their sizes and construct a kinetic equation that governs its evolution. This loop kinetic equation is formally analogous to Boltzmann's kinetic equation, with loop-loop reconnection described by a binary collision integral. A dimensionless parameter is introduced to scale the (unknown) overall rate of reconnection relative to Keplerian shear. After solving for the loop distribution function numerically, we calculate self-consistently the distribution of the mean magnetic pressure and dissipation rate with height, and the equilibrium shapes of loops of different sizes. We also compute the energy and torque associated with a given loop, as well as the total magnetic energy and torque in the corona. We explore the dependence of these quantities on the reconnection parameter and find that they can be greatly enhanced if reconnection between loops is suppressed.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Cassini States with Dissipation: Why Obliquity Tides Cannot Inflate Hot Jupiters

Some short-period exoplanets (hot Jupiters) are observed by their transits to have anomalously large radii. It has been suggested that these planets are in a resonance involving persistent misalignment and synchronous precession of their spin and orbital angular momenta, a Cassini state, and that the attendant tidal heating inflates the planet. We argue against this. Using explicit tidal integrations, we show that although an oblique Cassini state can dissipate many times the rotational energy of the planet, the rate of dissipation must be much less than hypothesized. Dissipation causes the planetary spin to lie at an angle to the plane containing the orbital and total angular momenta. If dissipation is too rapid, this angle becomes so large that Cassini equilibrium is lost. A separate consideration limits the total energy that can be extracted from the orbit. The source of the torque on the orbit, either an oblique parent star or an inclined third body, aligns with the orbit as it absorbs the angular momentum shed by the planet. Alignment removes the orbital precession required by the Cassini state. In combination with observational bounds on the mass and semimajor axis of a possible second planet and with bounds on the stellar rotation and obliquity, these constraints make it very unlikely that obliquity tides can be the explanation for inflated hot Jupiters, especially HD 209458b.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Ultra-relativistic geometrical shock dynamics and vorticity

Geometrical shock dynamics, also called CCW theory, yields approximate equations for shock propagation in which only the conditions at the shock appear explicitly; the post-shock flow is presumed approximately uniform and enters implicitly via a Riemann invariant. The nonrelativistic theory, formulated by G. B. Whitham and others, matches many experimental results surprisingly well. Motivated by astrophysical applications, we adapt the theory to ultra-relativistic shocks advancing into an ideal fluid whose pressure is negligible ahead of the shock, but one third of its proper energy density behind the shock. Exact results are recovered for some self-similar cylindrical and spherical shocks with power-law pre-shock density profiles. Comparison is made with numerical solutions of the full hydrodynamic equations. We review relativistic vorticity and circulation. In an ultrarelativistic ideal fluid, circulation can be defined so that it changes only at shocks, notwithstanding entropy gradients in smooth parts of the flow.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

Stochastic Coagulation and the Timescale for Runaway Growth

We study the stochastic coagulation equation using simplified models and efficient Monte Carlo simulations. It is known that (i) runaway growth occurs if the two-body coalescence kernel rises faster than linearly in the mass of the heavier particle; and (ii) for such kernels, runaway is instantaneous in the limit that the number of particles tends to infinity at fixed collision time per particle. Superlinear kernels arise in astrophysical systems where gravitational focusing is important, such as the coalescence of planetesimals to form planets or of stars to form supermassive black holes. We find that the time required for runaway decreases as a power of the logarithm of the the initial number of particles. Astrophysical implications are briefly discussed.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Radio scintillation of gamma-ray-burst afterglows

Stars twinkle to the eye through atmospheric turbulence, but planets, because of their larger angular size, do not. Similarly, scintillation due to the local interstellar medium will modulate the radio flux of gamma-ray-burst afterglows and may permit indirect measurements of their angular sizes. The amplitude of refractive scintillation is of order ten percent at ten gigahertz unless the source size is much larger than the expected size, of order ten microarcseconds. Diffractive scintillation is marginally possible, depending sensitively on the source size, observing frequency, and scattering measure of the interstellar medium.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Two-Component Fokker-Planck Models for the Evolution of Isolated Globular Clusters

Two-component (normal and degenerate stars) models are the simplest realization of clusters with a mass spectrum because high mass stars evolve quickly into degenerates, while low mass stars remain on the main-sequence for the age of the universe. Here we examine the evolution of isolated globular clusters using two-component Fokker-Planck (FP) models that include heating by binaries formed in tidal capture and in three-body encounters. Three-body binary heating dominates and the postcollapse expansion is self-similar, at least in models with total mass M <= 3 x 10^5 M_\odot, initial half-mass radius r_{h,i} >= 5 pc, component mass ratio m_2/m_1 <= 2, and number ratio N_1/N_2 <= 300 when m_2=1.4 M_\odot. We derive scaling laws for ρ_c, v_c, r_c, and r_h as functions of m_1/m_2, N, M, and time t from simple energy-balance arguments, and these agree well with the FP simulations. We have studied the conditions under which gravothermal oscillations (GTOs) occur. If E_{tot} and E_c are the energies of the cluster and of the core, respectively, and t_{rh} and t_c are their relaxation times, then ε\equiv (E_{tot}/t_{rh})/(E_c/t_{rc}) is a good predictor of GTOs: all models with ε>0.01 are stable, and all but one with ε< 0.01 oscillate. We derive a scaling law for εagainst N and m_1/m_2 and compared with our numerical results. Clusters with larger m_2/m_1 or smaller N are stabler.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

MEASURING SPIRAL ARM TORQUES: RESULTS FOR M100

Spiral arms, if they are massive, exert gravitational torques that transport angular momentum radially within galactic disks. These torques depend not on the pattern speed or permanence of the arms but only on the nonaxisymmetric mass distribution. Hence the torques can be measured from photometry. We demonstrate this using $gri$ CCD data for M100 (NGC 4321). Since we find consistency among the three bands, we believe that dust and young stars in the arms do not seriously bias our results. If the present epoch is representative, the timescale for redistribution of angular momentum in M100 is $5-10$ Gyr, the main uncertainty being the mass-to-light ratio of the disk.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

Nonlinear Damping of Oscillations in Tidal-Capture Binaries

We calculate the damping of quadrupole f and low order g modes (primary modes) by nonlinear coupling to other modes of the star. This damping is orders of magnitude more rapid than direct radiative damping when the primary amplitude is large, as in tidal capture. Primary modes destabilize high degree g-modes of half their frequency (daughter modes) by 3-mode coupling in radiative zones. In sunlike stars, the growth time $\equivη^{-1}\approx 4 E_{0,42}^{-1/2}$ days, where $E_{0,42}$ is the initial energy of the primary mode in units of $10^{42}~$erg, and of order $10^{10}E_{0,42}^{5/4}$ daughters are unstable. The growth rate is approximately equal to the angular frequency of the primary mode times its dimensionless radial amplitude, $δR/R_*\approx 0.002E_{0,42}^{1/2}$. Although the daughter modes are limited by their own nonlinearities, collectively they absorb most of the primary mode's energy after a time $\sim 10η^{-1}$ provided $E_{0}> 10^{40}~\mbox{erg}$. In fact nonlinear mode interaction may be the dominant damping process if $E_0\gtrsim 10^{37}~\mbox{erg}$. Our results have application to tidally captured main sequence globular cluster stars of mass $\ge 0.5 M_{\sun}$; the tidal energy is dissipated in the radiative core of the star in about a month, which is less than the initial orbital period.

0 citations0 reviewsNo code linkedOpen access