SMBH in Galactic Nuclei with Tidal Disruption of Stars

Tidal Disruption of stars by super massive central black holes from dense star clusters is modeled by high-accuracy direct $N$-body simulation. The time evolution of the stellar tidal disruption rate, the effect of tidal disruption on the stellar density profile and for the first time the detailed origin of tidally disrupted stars are carefully examined and compared with classic papers in the field. Up to 128k particles are used in simulation to model the star cluster around the super massive black hole, we use the particle number and the tidal radius of black hole as free parameters for a scaling analysis. The transition from full to empty loss-cone is analyzed in our data, the tidal disruption rate scales with the particle number $N$ in the expected way for both cases. For the first time in numerical simulations (under certain conditions) we can support the concept of a critical radius of Frank & Rees (1976), which claims that most stars are tidally accreted on highly eccentric orbits originating from regions far outside the tidal radius. Due to the consumption of stars moving on radial orbits, a velocity anisotropy is founded inside the cluster. Finally we make an estimation for the real galactic center based on our simulation results and the scaling analysis.Comment: 15 pages, 16 figures, accepted by Ap

Three-coloring graphs with no induced seven-vertex path II : using a triangle

In this paper, we give a polynomial time algorithm which determines if a given graph containing a triangle and no induced seven-vertex path is 3-colorable, and gives an explicit coloring if one exists. In previous work, we gave a polynomial time algorithm for three-coloring triangle-free graphs with no induced seven-vertex path. Combined, our work shows that three-coloring a graph with no induced seven-vertex path can be done in polynomial time.Comment: 26 page

Supermassive Black Holes in Galactic Nuclei with Tidal Disruption of Stars: Paper II - Axisymmetric Nuclei

Tidal Disruption of stars by supermassive central black holes from dense rotating star clusters is modelled by high-accuracy direct N-body simulation. As in a previous paper on spherical star clusters we study the time evolution of the stellar tidal disruption rate and the origin of tidally disrupted stars, now according to several classes of orbits which only occur in axisymmetric systems (short axis tube and saucer). Compared with that in spherical systems, we found a higher TD rate in axisymmetric systems. The enhancement can be explained by an enlarged loss-cone in phase space which is raised from the fact that total angular momentum $\bf J$ is not conserved. As in the case of spherical systems, the distribution of the last apocenter distance of tidally accreted stars peaks at the classical critical radius. However, the angular distribution of the origin of the accreted stars reveals interesting features. Inside the influence radius of the supermassive black hole the angular distribution of disrupted stars has a conspicuous bimodal structure with a local minimum near the equatorial plane. Outside the influence radius this dependence is weak. We show that the bimodal structure of orbital parameters can be explained by the presence of two families of regular orbits, namely short axis tube and saucer orbits. Also the consequences of our results for the loss cone in axisymmetric galactic nuclei are presented.Comment: 14 pages, 16 figures, accepted by Ap

The limits of Hamiltonian structures in three-dimensional elasticity, shells, and rods

This paper uses Hamiltonian structures to study the problem of the limit of three-dimensional (3D) elastic models to shell and rod models. In the case of shells, we show that the Hamiltonian structure for a three-dimensional elastic body converges, in a sense made precise, to that for a shell model described by a one-director Cosserat surface as the thickness goes to zero. We study limiting procedures that give rise to unconstrained as well as constrained Cosserat director models. The case of a rod is also considered and similar convergence results are established, with the limiting model being a geometrically exact director rod model (in the framework developed by Antman, Simo, and coworkers). The resulting model may or may not have constraints, depending on the nature of the constitutive relations and their behavior under the limiting procedure. The closeness of Hamiltonian structures is measured by the closeness of Poisson brackets on certain classes of functions, as well as the Hamiltonians. This provides one way of justifying the dynamic one-director model for shells. Another way of stating the convergence result is that there is an almost-Poisson embedding from the phase space of the shell to the phase space of the 3D elastic body, which implies that, in the sense of Hamiltonian structures, the dynamics of the elastic body is close to that of the shell. The constitutive equations of the 3D model and their behavior as the thickness tends to zero dictates whether the limiting 2D model is a constrained or an unconstrained director model. We apply our theory in the specific case of a 3D Saint Venant-Kirchhoff material andderive the corresponding limiting shell and rod theories. The limiting shell model is an interesting Kirchhoff-like shell model in which the stored energy function is explicitly derived in terms of the shell curvature. For rods, one gets (with an additional inextensibility constraint) a one-director Kirchhoff elastic rod model, which reduces to the well-known Euler elastica if one adds an additional single constraint that the director lines up with the Frenet frame