Contribution of dust produced by binary merger ejecta

By means of a population synthesis code and by constructing a simple toy model of dust produced by asymptotic giant branch(AGB) stars, common envelope(CE) ejecta and binary merger ejecta, we estimate the dust product rates(DPRs) of these processes in the Milky Way. The total DPR from AGB stars is～ 6.7 × 10-4M yr-1, in which about 73% of dust grains are carbon, 24% are silicates and 3% are iron. The total DPR from CE ejecta is ～ 4.2 × 10-4M yr-1, in which about 83% of dust grains are silicates, about 12% are carbon and 5% are iron.The DPR from binary merger ejecta is less than 1/3 that from AGB stars or CE ejecta,and it could even be negligible under certain circumstances. Therefore, compared with AGB stars and CE ejecta, the contribution of dust produced by binary merger ejecta to total dust grains in the Milky Way is smaller or can be negligible.     

Hubble sequence,angular momenta and time-scales of the early evolution of galaxies

The empirical evidence for a connection between type and relative angular momentum of galaxies is reviewed and some constraints for the theoretical explanation are discussed.

---

.

     

An axion-like scalar field environment effect on binary black hole merger

The environment, such as an accretion disk, could modify the signal of the gravitational wave from astrophysical black hole binaries. In this article, we model the matter field around intermediatemass binary black holes by means of an axion-like scalar field and investigate their joint evolution. In detail, we consider equal mass binary black holes surrounded by a shell of axion-like scalar field both in spherically symmetric and non-spherically symmetric cases, and with different strengths of the scalar field. Our result shows that the environmental scalar field could essentially modify the dynamics. Firstly,in the spherically symmetric case, with increase of the scalar field strength, the number of circular orbits for the binary black hole is reduced. This means that the scalar field could significantly accelerate the merger process. Secondly, once the scalar field strength exceeds a certain critical value, the scalar field could collapse into a third black hole with its mass being larger than that of the binary. Consequently,the new black hole that collapses from the environmental scalar field could accrete the binary promptly and the binary collides head-on with each other. In this process, there is almost no quadrupole signal produced, and, consequently, the gravitational wave is greatly suppressed. Thirdly, when the scalar field strength is relatively smaller than the critical value, the black hole orbit could develop eccentricity through accretion of the scalar field. Fourthly, during the initial stage of the inspiral, the gravitational attractive force from the axion-like scalar field could induce a sudden turn in the binary orbits, hence resulting in a transient wiggle in the gravitational waveform. Finally, in the non-spherical case, the scalar field could gravitationally attract the binary moving toward the center of mass for the scalar field and slow down the merger process.     

Problems of low-mass binary evolution

Some problems connected with low-mass binary evolution (from contact binaries to cataclysmic variables and origin of bursters) are considered. Most attention is given to contact W UMa-stars and to (still unclear) scenarios where the angular momentum loss by magnetic braking may, at least partly, control the contact binary evolution.     

Dynamical friction and galaxy merging time-scales

The time-scale for galaxies within merging dark matter haloes to merge with each other is an important ingredient in galaxy formation models. Accurate estimates of merging time-scales are required for predictions of astrophysical quantities such as black hole binary merger rates, the build-up of stellar mass in central galaxies and the statistical properties of satellite galaxies within dark matter haloes. In this paper, we study the merging time-scales of extended dark matter haloes using N -body simulations. We compare these results to standard estimates based on the Chandrasekhar theory of dynamical friction. We find that these standard predictions for merging time-scales, which are often used in semi-analytic galaxy formation models, are systematically shorter than those found in simulations. The discrepancy is approximately a factor of 1.7 for M _sat/ M _host≈ 0.1 and becomes larger for more disparate satellite-to-host mass ratios, reaching a factor of ∼3.3 for M _sat/ M _host≈ 0.01. Based on our simulations, we propose a new, easily implementable fitting formula that accurately predicts the time-scale for an extended satellite to sink from the virial radius of a host halo down to the halo's centre for a wide range of M _sat/ M _host and orbits. Including a central bulge in each galaxy changes the merging time-scale by ≲10 per cent. To highlight one concrete application of our results, we show that merging time-scales often used in the literature overestimate the growth of stellar mass by satellite accretion by ≈40 per cent, with the extra mass gained in low mass ratio mergers.     

Orbital evolution of small binary asteroids

About 15% of both near-Earth and main-belt asteroids with diameters below 10 km are now known to be binary. These small asteroid binaries are relatively uniform and typically contain a fast-spinning, flattened primary and a synchronously rotating, elongated secondary that is 20-40% as large (in diameter) as the primary. The principal formation mechanism for these binaries is now thought to be YORP (Yarkovsky-O’Keefe-Radzievskii-Paddack) effect induced spin-up of the primary followed by mass loss and accretion of the secondary from the released material. It has previously been suggested (?uk, M. [2007]. Astrophys. J. 659, L57-L60) that the present population of small binary asteroids is in a steady state between production through YORP and destruction through binary YORP (BYORP), which should increase or decrease secondary’s orbit, depending on the satellite’s shape. However, BYORP-driven evolution has not been directly modeled until now. Here we construct a simple numerical model of the binary’s orbital as well the secondary’s rotational dynamics which includes BYORP and selected terms representing main solar perturbations. We find that many secondaries should be vulnerable to chaotic rotation even for relatively low-eccentricity mutual orbits. We also find that the precession of the mutual orbit for typical small binary asteroids might be dominated by the perturbations from the prolate and librating secondary, rather than the oblate primary. When we evolve the mutual orbit by BYORP we find that the indirect effects on the binary’s eccentricity (through the coupling between the orbit and the secondary’s spin) dominate over direct ones caused by the BYORP acceleration. In particular, outward evolution causes eccentricity to increase and eventually triggers chaotic rotation of the secondary. We conclude that the most likely outcome will be reestablishing of the synchronous lock with a “flipped” secondary which would then evolve back in. For inward evolution we find an initial decrease of eccentricity and secondary’s librations, to be followed by later increase. We think that it is likely that various forms of dissipation we did not model may damp the secondary’s librations close to the primary, allowing for further inward evolution and a possible merger. We conclude that a merger or a tidal disruption of the secondary are the most likely outcomes of the BYORP evolution. Dissociation into heliocentric pairs by BYORP alone should be very difficult, and satellite loss might be restricted to the minority of systems containing more than one satellite at the time.     

Tidal evolution in close binary systems

The aim of the present paper will be to give a mathematical outline of the theory of tidal evolution in close binary systems of secularly constant total momentum — an evolution activated by viscous friction of dynamical tides raised by the two components on each other. The first section contains a general outline of the problem; and in Section 2 we shall establish the basic expressions for the energy and momenta of close binaries consisting of components of arbitrary internal structure. In Section 3 we shall investigate the maximum and minimum values of the energy (kinetic and potential) which such systems can attain for given amount of total momentum; while in Section 4 we shall compare these results with the actual facts encountered in binaries with components whose internal structure (and, therefore, rotational momenta) are known to us from evidence furnished by the observed rates of apsidal advance.The results show that all such systems — be these of detached or semi-detached type — disclose that more than 99% of their total momenta are stored in the orbital momentum. The sum of the rotational momenta of the constituent components amounts to less than a percent of the total — a situation characteristic of a state close to the minimum energy for given total momentum. This appears, moreover, to be true not only of the systems with both components on the Main Sequence, but also of those possessing evolved components in contact with their Roche limits.Under such conditions, a synchronism between rotation and revolution (characteristic of both extreme states of maximum and minimum energy) is not only possible, but appears to have been actually approached — if not attained — in the majority of cases. In other words, it would appear that — in at least a large majority of known cases — the existing close binaries have already attained orbits of maximum distension consistent with their momenta; and tidal evolution alone can no longer increase the present separations of the components to any appreciable extent.The virtual absence, in the sky, of binary systems intermediate between the stages of maximum and minimum energy for given momentum leads us to conjecture that the process of dynamical evolution activated by viscous tides may enroll on a time-scale which is relatively short in comparison with their total age — even for systems like Y Cygni or AG Persei, whose total age can scarcely exceed 10⁷ yr. A secular increase of the semi-major axes of relative orbits is dynamically coupled with a corresponding variation in the velocity of axial rotation of both components through the tidal lag arising from the viscosity of stellar material. The differential equations of so coupled a system are given in Section 5; but their solution still constitutes a task for the future.The Lunar Science Institute Contribution No. 90. The Lunar Science Institute is operated by the Universities Space Research Association under Contract No. NSR 09-051-001 with the National Aeronautics and Space Administration.     

Theoretical models for wide binary evolution

Summary Wide binary evolution under the influence of stars, giant molecular clouds, and the galactic tidal field is reviewed. Recent results show that binary lifetimes ata≳0.1 pc are determined by an interplay between relatively gentle evolution due to stellar encounters and catastrophic collisions with cloud subclumps. It is emphasized that binary semimajor axis distributions may not be understood from lifetimes alone. To determine the theoretical distribution binary ‘birthrate functions’ must be convolved with probability distributions for the evolution of an ensemble of binaries from given initial separations. Simple models show no sharp breaks below the tidal cut-off imposed by the Galaxy.     

A low-mass-ratio and deep contact binary as the progenitor of the merger V1309 Sco

Nova Sco 2008(=V1309 Sco) is an example of a V838 Mon type eruption rather than a typical classical nova. This enigmatic object was recently shown to have resulted from the merger of two stars in a contact binary. It is the first stellar merger that was identified to be undergoing a common envelope transient. To understand the properties of its binary progenitor, the pre-outburst light curves were analyzed by using the W-D method. The photometric solution of the 2002 light curve shows that it is a deep contact binary(f = 89.5(±40.5)%) with a mass ratio of 0.094. The asymmetry of the light curve is explained by the presence of a dark spot on the more massive component. The extremely high fill-out factor suggests that the merging of the contact binary is driven by dynamical mass loss from the outer Lagrange point. However,the analysis of the 2004 light curve indicates that no solutions were obtained even at an extremely low mass ratio of q = 0.03. This suggests that the common convective envelope of the binary system disappeared and the secondary component spiraled into the envelope of the primary in 2004. Finally, the ejection of the envelope of the primary produced the outburst.     

Tidal fluxes and evolution of binary star systems

Explicit analytical expressions for the tidal velocities in the stars which are components of the binary system are obtained, taking into account the viscosity of the stellar matter, the star's rotation, the eccentricity of the orbit and tilt of the axis of rotation to the orbital plane.Equations which determine the tidal evolution of the star's axial and orbital motions are obtained and analyzed. Numerical solutions of these equations are analyzed.     

Local subgiants and time-scales of disc formation

This paper concentrates on the relationship between the rate of gas emission from the nucleus of Comet 9P/Tempel 1, the fraction f of the nucleus that is active, and the crater damage inflicted by the recent 2005 July 4 Deep Impact space mission. The cometary nucleus has a surface area of about  1.7 × 10⁸ m²  and a mean radius of about 3700 m. Before the impact it is estimated that only a fraction f = 0.0056 of the nucleus surface was actively producing gas and dust. The active area was about  9.4 × 10⁵ m²  . Absolute magnitudes obtained at recent perihelion passages of this comet indicate that variations in the 0.0074 > f > 0.0039 range can occur from apparition to apparition. Because of the low size of the original active area, the production of a new impact crater in the diameter range 40 to 300 m would lead to a long-term change in the cometary visual magnitude in the range 0.0018 to 0.098 respectively. This is below the limit of detectability. It has been suggested that the cometary dust is in the form of 'talcum powder' not 'beach sand'. We suggest that the dust ejected from the impact site has been broken up by the energetic impact process and thus has a different size distribution from dust locked in the snowy matrix of the nucleus and normally lifted off the nucleus by gentle sublimation processes.     

A binary merger model for the formation of the Supernova 1987A triple-ring nebula

We examine a binary merger model for the formation of the mysterious triple-ring nebula surrounding Supernova 1987A, which still has not been convincingly explained in detailed hydrodynamical calculations. During the merger of 15 and  5 M_⊙  binary systems, mass is ejected primarily at mid-latitudes for a sufficiently evolved primary, as demonstrated by Morris & Podsiadlowski. This material is swept up by the fast wind of the central star during its post-merger blue supergiant phase, leading to a density contrast of ∼150 in the outer rings at the time of the supernova. The equatorial ring probably formed later when the star contracted to become a blue supergiant. The asymmetry between the northern and southern outer rings can be explained by a 10 per cent asymmetry during the merger, perhaps due to a pulsational instability in the common envelope.
We present a parameter study from which we determine a mass-loss rate in the blue supergiant wind in the range  1.5–3 × 10⁻⁷ M_⊙ yr⁻¹  in agreement with previous estimates. The morphology of the best model is consistent with the well-known Hubble Space Telescope image at better than 5 per cent and is also in broad agreement with light-echo observations. The circumstellar environment on larger scales (up to 3 pc) is also investigated. We conclude with a brief discussion of the bipolar nebulae surrounding the Galactic stars, Sheridan 25, HD 168625 and η Carinae.     

CaseA evolution of massive close binary systems

The caseA evolution of close binary systems with total mass of 20.4M is investigated by following the evolution of both components simultaneously. The evolution is followed up to the stage at which a system overflows the outer critical surface or evolves into the phase of modeBr mass-transfer. It is found that the evolution of the systems can be classified into six types. The ranges of initial parameters which lead systems to each type of evolution are shown on the initial-parameter plane. The evolutionary features of each evolution type are described in detail.     

Tidal evolution of close binary asteroid systems

We provide a generalized discussion of tidal evolution to arbitrary order in the expansion of the gravitational potential between two spherical bodies of any mass ratio. To accurately reproduce the tidal evolution of a system at separations less than 5 times the radius of the larger primary component, the tidal potential due to the presence of a smaller secondary component is expanded in terms of Legendre polynomials to arbitrary order rather than truncated at leading order as is typically done in studies of well-separated system like the Earth and Moon. The equations of tidal evolution including tidal torques, the changes in spin rates of the components, and the change in semimajor axis (orbital separation) are then derived for binary asteroid systems with circular and equatorial mutual orbits. Accounting for higher-order terms in the tidal potential serves to speed up the tidal evolution of the system leading to underestimates in the time rates of change of the spin rates, semimajor axis, and mean motion in the mutual orbit if such corrections are ignored. Special attention is given to the effect of close orbits on the calculation of material properties of the components, in terms of the rigidity and tidal dissipation function, based on the tidal evolution of the system. It is found that accurate determinations of the physical parameters of the system, e.g., densities, sizes, and current separation, are typically more important than accounting for higher-order terms in the potential when calculating material properties. In the scope of the long-term tidal evolution of the semimajor axis and the component spin rates, correcting for close orbits is a small effect, but for an instantaneous rate of change in spin rate, semimajor axis, or mean motion, the close-orbit correction can be on the order of tens of percent. This work has possible implications for the determination of the Roche limit and for spin-state alteration during close flybys.     

Particle-acceleration time-scales in TeV blazar flares

Observations of minute-scale flares in TeV Blazars place constraints on particle-acceleration mechanisms in those objects. The implications for a variety of radiation mechanisms have been addressed in the literature; in this paper, we compare four different acceleration mechanisms: diffusive shock acceleration, second-order Fermi, shear acceleration and the converter mechanism. When the acceleration time-scales and radiative losses are taken into account, we can exclude shear acceleration and the neutron-based converted mechanism as possible acceleration processes in these systems. The first-order Fermi process and the converter mechanism working via synchrotron self-Compton (SSC) photons are still practically instantaneous, however, provided sufficient turbulence is generated on the time-scale of seconds. We propose stochastic acceleration as a promising candidate for the energy-dependent time delays in recent gamma-ray flares of Markarian 501.     

Spin evolution of neutron stars in binary systems

I review our understanding of the evolution of the spin periods of neutron stars in binary stellar systems, from their birth as fast, spin-powered pulsars, through their middle life as accretion-powered pulsars, upto their recycling or “rebirth” as spin-powered pulsars with relatively low magnetic fields and fast rotation. I discuss how the new-born neutron star is spun down by electromagnetic and “propeller” torques, until accretion of matter from the companion star begins, and the neutron star becomes an accretion-powered X-ray pulsar. Detailed observations of massive radio pulsar binaries like PSR 1259-63 will yield valuable information about this phase of initial spindown. I indicate how the spin of the neutron star then evolves under accretion torques during the subsequent phase as an accretion-powered pulsar. Finally, I describe how the neutron star is spun up to short periods again during the subsequent phase of recycling, with the accompanying reduction in the stellar magnetic field, the origins of which are still not completely understood.     

Tidal evolution in close binary systems,II

The aim of the present investigation will be to determine the explicit forms of differential equations which govern secular perturbations of the orbital elements of close binary systems in the plane of the orbit (i.e., of the semi-major axisA, eccentricitye, and longitude of the periastron ), arising from the lag of dynamical tides due to viscosity of stellar material. The results obtained are exact for any value of orbital eccentricity comprised between 0e<1; and include the effects produced by the second, third and fourth-harmonic dynamical tides, as well as by axial rotation with arbitrary inclination of the equator to the orbital plane.In Section 2 following brief introductory remarks the variational equations of the problem of plane motion will be set up in terms of the rectangular componentsR, S, W of disturbing accelerations with respect to a revolving system of coordinates. The explicit form of these coefficients will be established in Section 3 to the degree of accuracy to which squares and higher powers of quantities of the order of superficial distortion can be ignored. Section 4 will be devoted to a derivation of the explicit form of the variational equations for the case of a perturbing function arising from axial rotation; and in Section 5 we shall derive variational equations which govern the perturbation of orbital elements caused by lagging dynamical tides.Numerical integrations of these equations, which govern the tidal evolution of close binary systems prompted by viscous friction at constant mass, are being postponed for subsequent investigations.Prepared at the Lunar Science Institute, Houston, Texas, under the joint support of the Universities Space Research Association, Charlottesville, Virginia, and the National Aeronautics and Space Administration Manned Spacecraft Center, Houston, Texas, under Contract No. NSR 09-051-001. This paper constitutes Lunar Science Institute Contribution no. 100.Normally at the Department of Astronomy, University of Manchester, England.