Slow modes in stellar systems with nearly harmonic potentials: II. Gravitational loss-cone instability

Using a consistent perturbation theory for collisionless disk-like and spherical star clusters, we construct a theory of slow modes for systems having an extended central region with a nearly harmonic potential due to the presence of a fairly homogeneous (on the scales of the stellar system) heavy, dynamically passive halo. In such systems, the stellar orbits are slowly precessing, centrally symmetric ellipses (2: 1 orbits). We consider star clusters with monoenergetic distribution functions that monotonically increase with angular momentum in the entire range of angular momenta (from purely radial orbits to circular ones) or have a growing region only at low angular momenta. In these cases, there are orbits with a retrograde precession, i.e., in a direction opposite to the orbital rotation of the star. The presence of a gravitational loss-cone instability, which is also observed in systems of 1: 1 orbits in near-Keplerian potentials, is associated with such orbits. In contrast to 1: 1 systems, the loss-cone instability takes place even for distribution functions monotonically increasing with angular momentum, including those for systems with circular orbits. The regions of phase space with retrograde orbits do not disappear when the distribution function is smeared in energy. We investigate the influence of a weak inhomogeneity of a heavy halo with a density that decreases with distance from the center.     

Gravitational loss-cone instability in stellar systems with retrograde orbit precession

We study spherical and disc clusters in a near-Keplerian potential of galactic centres or massive black holes. In such a potential orbit precession is commonly retrograde, that is, the direction of the orbit precession is opposite to the orbital motion. It is assumed that stellar systems consist of nearly-radial orbits. We show that if there is a loss-cone at low angular momentum (e.g. due to consumption of stars by a black hole), an instability similar to loss-cone instability in plasma may occur. The gravitational loss-cone instability is expected to enhance black hole feeding rates. For spherical systems, the instability is possible for the number of spherical harmonics   l ≥ 3  . If there is some amount of counter-rotating stars in flattened systems, they generally exhibit the instability independent of azimuthal number m . The results are compared with those obtained recently by Tremaine for distribution functions monotonically increasing with angular momentum.
The analysis is based on simple characteristic equations describing small perturbations in a disc or a sphere of stellar orbits highly elongated in radius. These characteristic equations are derived from the linearized Vlasov equations (combining the collisionless Boltzmann kinetic equation and the Poisson equation), using the action-angle variables. We use two techniques for analysing the characteristic equations: the first one is based on preliminary finding of neutral modes, and the second one employs a counterpart of the plasma Penrose–Nyquist criterion for disc and spherical gravitational systems.     

Semiclassical approximation for low-degree stellar p modes – I. The classical eigenfrequency equation

A new eigenfrequency equation for low-degree solar-like oscillations in stars is developed, based on the assumption of purely classical propagation in the stellar interior of acoustic waves modified by buoyancy and gravity . Compared with high-frequency asymptotic analysis, the eigenfrequency equation has a new functional form, with expansion in powers of ℓ(ℓ+1) instead of 1/ ω . Basic observable quantities, the 'large' and 'small' frequency separations , are interpreted as the dependence on frequency and refraction angle of a classical action integral for waves propagating close to the stellar diameter. The new eigenfrequency equation gives a significant improvement in accuracy over previous analyses when tested with solar p modes, suggesting this as an alternative and more powerful tool for applications in stellar seismology.     

Dissipative models of colliding stellar winds — I. Effects of thermal conduction in wide binary systems

The influence of electron thermal conduction on the 2D gas dynamics of colliding stellar winds is investigated. It is shown that, as a result of the non-linear dependence of the electron thermal flux on the temperature, the pre-heating zones (in which the hot gas in the interaction region heats the cool winds in front of the shocks) have finite sizes. The dependence of the problem of the structure of the flow in the interaction region on the dimensionless parameters is studied, and a simple expression is derived for the size of the pre-heating zones at the axis of symmetry. It is shown that small values of the thermal conductivity do not suppress the Kelvin–Helmholtz instability if the adiabatic flow is subject to it. Further studies, both numerical and analytical, in this direction will be of great interest. The influence of thermal conduction on the X-ray emission from the interaction region is also estimated.     

Realistic simulations of stellar surface convection with ANTARES: I. Boundary conditions and model relaxation

We have implemented open boundary conditions into the ANTARES code to increase the realism of our simulations of stellar surface convection. Even though we greatly benefit from the high accuracy of our fifth order numerical scheme (WENO5) the broader stencils needed for the numerical scheme complicate the implementation of boundary conditions. We show that the effective temperature of a numerical simulation cannot be changed by corrections at the lower boundary since the thermal stratification does only change on the Kelvin–Helmholtz time scale. Except for very shallow models this time scale cannot be covered by multidimensional simulations due to the enormous computational requirements. We demonstrate to what extent numerical simulations of stellar surface convection are sensitive to the initial conditions and the boundary conditions. An ill-conceived choice of parameters for the boundary conditions can have a severe impact. Numerical simulations of stellar surface convection will only be (physically) meaningful and realistic if the initial model, the extent and position of the simulation box and the parameters from the boundary conditions are chosen adequately.     

Monte Carlo simulations of star clusters – II. Tidally limited, multimass systems with stellar evolution

A revision of Stodółkiewicz's Monte Carlo code is used to simulate evolution of large star clusters. The new method treats each superstar as a single star and follows the evolution and motion of all individual stellar objects. A survey of the evolution of N -body systems influenced by the tidal field of a parent galaxy and by stellar evolution is presented. The process of energy generation is realized by means of appropriately modified versions of Spitzer's and Mikkola's formulae for the interaction cross-section between binaries and field stars and binaries themselves. The results presented are in good agreement with theoretical expectations and the results of other methods (Fokker–Planck, Monte Carlo and N -body). The initial rapid mass loss, resulting from stellar evolution of the most massive stars, causes expansion of the whole cluster and eventually leads to the disruption of less bound systems ( W ₀=3). Models with larger W ₀ survive this phase of evolution and then undergo core collapse and subsequent post-collapse expansion, like isolated models. The expansion phase is eventually reversed when tidal limitation becomes important. The results presented are the first major step in the direction of simulating evolution of real globular clusters by means of the Monte Carlo method.     

The structure of clusters with bimodal distributions of galaxy radial velocities,I. A1035

The structure of the A1035 cluster of galaxies (10^h32^m + 40°13′, cz ～ 22000 kms^?1), which exhibits a bimodal distribution of galaxy radial velocities (ΔV ≈ 3000 kms^?1), is analyzed using three methods of determining the relative distances to clusters from early-type galaxies: the Kormendy relation corrected for the dependence of residuals on galaxy magnitude, the photometric plane, and the fundamental plane. We use the data obtained with the l-m telescope of the Special Astrophysical Observatory of the Russian Academy of Sciences and SDSS (DR5) data to show that A1035 consists of two gravitationally unbound independent clusters. These clusters with the velocity dispersions of 566 km s^?1 and 610 km s^?1 and masses within R ₂₀₀ equal to 2.7 · 10¹⁴ and 3.5 · 10¹⁴ M _⊙, respectively, obey the Hubble law.     

Searching for exoplanets with HEPS:I.detection probability of Earth-like planets in multiple systems

The astrometry method has great advantages in searching for exoplanets in the habitable zone around solar-like stars. However, the presence of multiple planets may cause a problem with degeneracy when trying to compute accurate planet parameters from observation data and reduce detectability. The degeneracy problem is extremely critical, especially in a space mission which has limited observation time and cadence. In this series of papers, we study the detectability of habitable Earth-mass planets in different types of multi-planet systems, aiming to find the most favorable targets for the potential space mission–Habitable ExoPlanet Survey(HEPS). In the first paper, we present an algorithm to find planets in the habitable zone around solar-like stars using astrometry. We find the detectability can be well described by planets' signal-to-noise ratio(SNR) and a defined parameter S = M2/(T1-T2)2, where M2 and T2are the mass and period of the second planet, respectively. T1 is the period of the planet in the habitable zone. The parameter S represents the influence of planetary architectures. We fit the detectability as a function of both the SNR of the planet in the habitable zone and the parameter S. An Earth-like planet in a habitable zone is harder to detect(with detectability PHP 80%) in a system with a hot Jupiter or warm Jupiter(within2 AU), in which the parameter S is large. These results can be used in target selections and to determine the priority of target stars for HEPS, especially when we select and rank nearby planet hosts with a single planet.     

A 2.5-dimensional viscous,resistive, advective magnetized accretion-outflow coupling in black hole systems: a higher order polynomial approximation

The correlated and coupled dynamics of accretion and outflow around black holes(BHs) are essentially governed by the fundamental laws of conservation as outflow extracts matter, momentum and energy from the accretion region. Here we analyze a robust form of 2.5-dimensional viscous, resistive,advective magnetized accretion-outflow coupling in BH systems. We solve the complete set of coupled MHD conservation equations self-consistently, through invoking a generalized polynomial expansion in two dimensions. We perform a critical analysis of the accretion-outflow region and provide a complete quasi-analytical family of solutions for advective flows. We obtain the physically plausible outflow solutions at high turbulent viscosity parameter α(0.3), and at a reduced scale-height, as magnetic stresses compress or squeeze the flow region. We found that the value of the large-scale poloidal magnetic field BPis enhanced with the increase of the geometrical thickness of the accretion flow. On the other hand,differential magnetic torque(-r~2 B_φB_z) increases with the increase in M.B_P,-r~2 B_φB_z as well as the plasma beta β_P get strongly augmented with the increase in the value of α, enhancing the transport of vertical flux outwards. Our solutions indicate that magnetocentrifugal acceleration plausibly plays a dominant role in effusing out plasma from the radial accretion flow in a moderately advective paradigm which is more centrifugally dominated. However in a strongly advective paradigm it is likely that the thermal pressure gradient would play a more contributory role in the vertical transport of plasma.     

The range of validity of the two-body approximation in models of terrestrial planet accumulation: I. Gravitational perturbations

Gravitational perturbations in semimajor axis, eccentricity, and inclination resulting from close planetesimal encounters (near 1 AU) out to 10 Tisserand sphere of influence radii were calculated by two- and three-dimensional numerical integration. These are compared with the results of treating the encounter as a two-body problem, as is customary in Monte Carlo calculations of orbital evolution and in numerical and analytical studies of planetary accumulation. It is found that for values of $\text{(}\text{V}\text{V}{}_{\mathrm{e}}\text{) ? 0.35 (V =}\text{relative velocity}\text{, V}{}_{\mathrm{e}}\text{=}\text{escape velocity of largest body}\text{)}$, the two-body body approximation fails to describe the outcome of individual encounters. In this low-velocity region, the two-body “gravitational focusing” cross section is no longer valid; “anomalous gravitational focusing” often leads to bodies on distant unperturbed trajectories becoming close encounters and vice versa. In spite of these differences, average perturbations given by the two-body approximation are valid within a factor of 2 when $\text{V}\text{V}{}_{\mathrm{e}}\text{> 0.07}$. In this same velocity range the “Arnold extrapolation,” whereby a few very close encounters are used to estimate the effect of many more distant encounters, is found to be a useful approximation.     

A new titanium-bearing calcium aluminosilicate phase: I. Meteoritic occurrences and formation in synthetic systems

Abstract A new titanium-bearing calcium aluminosilicate mineral has been identified in coarse-grained calcium-aluminum-rich inclusions (CAIs) from carbonaceous chondrites. The formula for this phase, which we have temporarily termed “UNK,” is Ca₃Ti(Al, Ti)₂(Si, Al)₃O¹⁴, and it is present in at least 8 of the 20 coarse-grained CAIs from the Allende CV3 chondrite examined as part of this project. The phase occurs in Types A and B1 inclusions as small tabular crystals oriented along two mutually perpendicular planes in melilite. UNK crystallizes from melts in dynamic crystallization experiments conducted in air from four bulk compositions modeled after Types A, B1, B2 and C inclusions. Cooling rates resulting in crystallization of UNK ranged from 0.5 to 200 °C/h from maximum (initial) temperatures of 1375 to 1580 °C. Only below 1190 °C does UNK itself begin to crystallize. To first order, the presence or absence of UNK from individual experiments can be understood in terms of the compositions of residual melts and nucleation probabilities. Compositions of synthetic and meteoritic UNK are very similar in terms of major oxides, differing only in the small amounts of trivalent Ti (7–13% of total Ti) in meteoritic samples. UNK crystallized from the Type A analog is similar texturally to that found in CAIs, although glass, which is typically associated with synthetic UNK, is not observed in meteoritic occurrences. A low Ti end-member of UNK (“Si-UNK”) with a composition near that of Ca₃Al₂Si₄O₁₄ was produced in a few samples from the Type B1 analog. This phase has not been found in the meteoritic inclusions.     

Space Plasma Phenomena : Laboratory Modeling I. Fields,Streams, Matter in Plasma

Earth, Moon, and Planets -     

Close binary systems in multiple stars: III. Eclipsing variables with LITE effect

The list of 240 eclipsing variables with the light equation (LITE effect) has been compiled. Some characteristics of multiple systems have been determined assuming that those stars may contain extra components. Quantity distributions have been obtained for Close Binary Systems (CBSs) with LITE effect and for those without it. The observed excess of short period stars (<0.3 days) is shown to be realistic. The eccentricities of Long Periodic Orbits (LPOs) have a parabolic distribution with a weak maximum at e = 0.28. The ratios of LPO to CBS periods cannot be considered to be random. The dependence of those values on LPO periods has been revealed. The distributions of periods and semimajor axes of LPOs are shown to have maxima at P ₃ = 32 yrs and a = 16 AU, respectively. The latter distribution can be described by the Gaussian function. Approximately 85% of tertiary companions in stellar systems have a mass which is half as much as the total CBS mass. Except for one CBS (RR Lyn), all multiple systems of the obtained list are dynamically stable configurations. The possibilities to detect tertiary companions using speckle interferometry observations and proper motions of CBSs, as well as by extra light contribution when analyzing light curves and radial velocity changes of the CBS center of mass are discussed. Other causes are considered which may be a reason for the LITE effect, such as the mass exchange in CBSs and the influence of the magnetic field of components. Certain difficulties may arise when interpreting the LITE effect caused by the magnetic field of components. The most promising research directions in detecting tertiary companions in CBSs with the LITE effect are suggested.