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

Eugene Vasiliev

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

ResearcherUnclaimed source record
astro-ph.GAastro-ph.HEgr-qcastro-ph.CO

Catalog footprint

What is connected

21works
4topics
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

21 published item(s)

preprint2022arXiv

Measuring the Milky Way mass distribution in the presence of the LMC

The ongoing interaction between the Milky Way (MW) and its largest satellite - the Large Magellanic Cloud (LMC) - creates a significant perturbation in the distribution and kinematics of distant halo stars, globular clusters and satellite galaxies, and leads to biases in MW mass estimates from these tracer populations. We present a method for compensating these perturbations for any choice of MW potential by computing the past trajectory of LMC and MW and then integrating the orbits of tracer objects back in time until the influence of the LMC is negligible, at which point the equilibrium approximation can be used with any standard dynamical modelling approach. We add this orbit-rewinding step to the mass estimation approach based on simultaneous fitting of the potential and the distribution function of tracers, and apply it to two datasets with the latest Gaia EDR3 measurements of 6d phase-space coordinates: globular clusters and satellite galaxies. We find that models with LMC mass in the range (1-2)x10^11 Msun better fit the observed distribution of tracers, and measure MW mass within 100 kpc to be (0.75+-0.1)x10^12 Msun, while neglecting the LMC perturbation increases it by ~15%.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Non-Parametric Spherical Jeans Mass Estimation with B-splines

Spherical Jeans modeling is widely used to estimate mass profiles of systems from star clusters to galactic stellar haloes to clusters of galaxies. It derives the cumulative mass profile, M(<r), from kinematics of tracers of the potential under the assumptions of spherical symmetry and dynamical equilibrium. We consider the application of Jeans modeling to mapping the dark matter distribution in the outer reaches of the Milky Way using field halo stars. We present a novel non-parametric routine for solving the spherical Jeans equation by fitting B-splines to the velocity and density profiles of halo stars. While most implementations assume parametric forms for these profiles, B-splines provide non-parametric fitting curves with analytical derivatives. Our routine recovers the mass profiles of equilibrium systems with flattened haloes or a stellar disc and bulge excellently (<~ 10% error at most radii). Tests with non-equilibrium, Milky Way-like galaxies from the Latte suite of FIRE-2 simulations perform quite well (<~ 15% error for r <~ 100 kpc). We also create observationally motivated datasets for the Latte suite by imposing selection functions and errors on phase space coordinates characteristic of Gaia and the DESI Milky Way Survey. The resulting imprecise and incomplete data require us to introduce an MCMC-based subroutine to obtain deconvolved density and velocity dispersion profiles from the tracer population. With these observational effects taken into account, the accuracy of the Jeans mass estimate remains at the level 20% or better.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Radialization of satellite orbits in galaxy mergers

We consider the orbital evolution of satellites in galaxy mergers, focusing on the evolution of eccentricity. Using a large suite of N-body simulations, we study the phenomenon of satellite orbital radialization -- a profound increase in the eccentricity of its orbit as it decays under dynamical friction. While radialization is detected in a variety of different setups, it is most efficient in the cases of high satellite mass, not very steep host density profiles, and high initial eccentricity. To understand the origin of this phenomenon, we run additional simulations with various physical factors selectively turned off: satellite mass loss, reflex motion and distortion of the host, etc. We find that all these factors are important for radialization, since it does not occur for point-mass satellites or when the host potential is replaced with an unperturbed initial profile. The analysis of forces and torques acting on both galaxies confirms the major role of self-gravity of both host and satellite in the reduction of orbital angular momentum. The classical Chandrasekhar dynamical friction formula, which accounts only for the forces between the host and the satellite, but not for internal distortions of both galaxies, does not match the evolution of eccentricity observed in N-body simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Local Group Mass in the light of Gaia

High accuracy proper motions (PMs) of M31 and other Local Group satellites have now been provided by the {\it Gaia} satellite. We revisit the Timing Argument to compute the total mass $M$ of the Local Group from the orbit of the Milky Way and M31, allowing for the Cosmological Constant. We rectify for a systematic effect caused by the presence of the Large Magellanic Cloud (LMC). The interaction of the LMC with the Milky Way induces a motion towards the LMC. This contribution to the measured velocity of approach of the Milky Way and M31 must be removed. We allow for cosmic bias and scatter by extracting correction factors tailored to the accretion history of the Local Group. The distribution of correction factors is centered around $0.63$ with a scatter $\pm 0.2$, indicating that the Timing Argument significantly overestimates the true mass. Adjusting for all these effects, the estimated mass of the Local Group is $ M = 3.4^{+1.4}_{-1.1} \times 10^{12} M_{\odot}$ (68 % CL) when using the M31 tangential velocity $ 82^{+38}_{-35}$ km/s. Lower tangential velocity models with $59^{+42}_{-38}$ km/s (derived from the same PM data with a flat prior on the tangential velocity) lead to an estimated mass of $ M = 3.1^{+1.3}_{-1.0} \times 10^{12} M_{\odot}$ (68 % CL). By making an inventory of the total mass associated with the 4 most substantial LG members (the Milky Way, M31, M33 and the LMC), we estimate the known mass is in the range $3.7^{+0.5}_{-0.5} \times 10^{12} \, M_{\odot}$.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The stellar mass distribution of the Milky Way's bar: an analytic model

We present an analytic model of the stellar mass distribution of the Milky Way bar. The model is obtained by fitting a multi-component parametric density distribution to a made-to-measure N-body model of Portail et al., constructed to match a variety of density and kinematics observational data. The analytic model reproduces in detail the 3D density distribution of the N-body bar including the X-shape. The model and the gravitational potential it generates are available as part of the software package AGAMA for galactic dynamics, and can be readily used for orbit integrations, hydrodynamical simulations or other applications.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Tango for three: Sagittarius, LMC, and the Milky Way

We assemble a catalogue of candidate Sagittarius stream members with 5d and 6d phase-space information, using astrometric data from Gaia DR2, distances estimated from RR Lyrae stars, and line-of-sight velocities from various spectroscopic surveys. We find a clear misalignment between the stream track and the direction of the reflex-corrected proper motions in the leading arm of the stream, which we interpret as a signature of a time-dependent perturbation of the gravitational potential. A likely cause of this perturbation is the recent passage of the most massive Milky Way satellite - the Large Magellanic Cloud (LMC). We develop novel methods for simulating the Sagittarius stream in the presence of the LMC, using specially tailored N-body simulations and a flexible parametrization of the Milky Way halo density profile. We find that while models without the LMC can fit most stream features rather well, they fail to reproduce the misalignment and overestimate the distance to the leading arm apocentre. On the other hand, models with an LMC mass in the range (1.3+-0.3)x10^11 Msun rectify these deficiencies. We demonstrate that the stream can not be modelled adequately in a static Milky Way. Instead, our Galaxy is required to lurch toward the massive in-falling Cloud, giving the Sgr stream its peculiar shape and kinematics. By exploring the parameter space of Milky Way potentials, we determine the enclosed mass within 100 kpc to be (5.6+-0.4)x10^11 Msun, and the virial mass to be (9.0+-1.3)x10^11 Msun, and find tentative evidence for a radially-varying shape and orientation of the Galactic halo.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Rates of Stellar Tidal Disruption

Tidal disruption events occur rarely in any individual galaxy. Over the last decade, however, time-domain surveys have begun to accumulate statistical samples of these flares. What dynamical processes are responsible for feeding stars to supermassive black holes? At what rate are stars tidally disrupted in realistic galactic nuclei? What may we learn about supermassive black holes and broader astrophysical questions by estimating tidal disruption event rates from observational samples of flares? These are the questions we aim to address in this Chapter, which summarizes current theoretical knowledge about rates of stellar tidal disruption, and compares theoretical predictions to the current state of observations.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A new implementation of the Schwarzschild method for constructing observationally-driven dynamical models of galaxies of all morphological types

We present Forstand, a new code for constructing dynamical models of galaxies with the Schwarzschild orbit-superposition method. These models are constrained by line-of-sight kinematic observations and applicable to galaxies of all morphological types, including disks and triaxial rotating bars. Our implementation has several novel and improved features, is computationally efficient, and made publicly available. Using mock datasets taken from N-body simulations, we demonstrate that the pattern speed of a bar can be recovered with an accuracy of 10-20%, regardless of orientation, if the 3D shape of the galaxy is known or inferred correctly.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Schwarzschild modeling of barred galaxies

We review the Schwarzschild orbit-superposition approach and present a new implementation of this method, which can deal with a large class of systems, including rotating barred disk galaxies. We discuss two conceptuals problems in this field: the intrinsic degeneracy of determining the potential from line-of-sight kinematics, and the non-uniqueness of deprojection and related biases in potential inference, especially acute for triaxial bars. When applied to mock datasets with known 3d shape, our method correctly recovers the pattern speed and other potential parameters. However, more work is needed to systematically address these two problems for real observational datasets.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Using Gaia for studying Milky Way star clusters

We review the implications of the Gaia Data Release 2 catalogue for studying the dynamics of Milky Way globular clusters, focusing on two separate topics. The first one is the analysis of the full 6-dimensional phase-space distribution of the entire population of Milky Way globular clusters: their mean proper motions (PM) can be measured with an exquisite precision (down to 0.05 mas/yr, including systematic errors). Using these data, and a suitable ansatz for the steady-state distribution function (DF) of the cluster population, we then determine simultaneously the best-fit parameters of this DF and the total Milky Way potential. We also discuss possible correlated structures in the space of integrals of motion. The second topic addresses the internal dynamics of a few dozen of the closest and richest globular clusters, again using the Gaia PM to measure the velocity dispersion and internal rotation, with a proper treatment of spatially correlated systematic errors. Clear rotation signatures are detected in 10 clusters, and a few more show weaker signatures at a level >=0.05$ mas/yr. PM dispersion profiles can be reliably measured down to 0.1 mas/yr, and agree well with the line-of-sight velocity dispersion profiles from the literature.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Tidal disruption rates in non-spherical galactic nuclei formed by galaxy mergers

We explore the rates of tidal disruption events (TDEs) of stars by supermassive black holes (SBHs) in galactic nuclei formed in mergers followed by a formation and coalescence of a binary SBH. Such systems initially have a deficit of stars on low-angular-momentum orbits caused by the slingshot process during the binary SBH stage, which tends to reduce the flux of stars into the SBH compared to the steady-state value. On the other hand, a newly formed galactic nucleus has a non-spherical shape which enhances the mixing of stars in angular momentum and thus the TDE rate. In galaxies with relatively low SBH masses (<10^7 Msun), relaxation times are short enough to wash out the anisotropy in initial conditions, and for more massive SBH the enhancement of flux due to non-sphericity turns out to be more important than its suppression due to initial anisotropy. Therefore, the present-day TDE rates generally exceed conventional steady-state estimates based on a spherical isotropic approximation. We thus conjecture that the lower observationally inferred TDE rates compared to theoretical predictions cannot be attributed to the depletion of low-angular-momentum stars by SBH binaries.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Applying Schwarzschild's orbit superposition method to barred or non-barred disc galaxies

We present an implementation of the Schwarzschild orbit superposition method which can be used for constructing self-consistent equilibrium models of barred or non-barred disc galaxies, or of elliptical galaxies with figure rotation. This is a further development of the publicly available code SMILE; its main improvements include a new efficient representation of an arbitrary gravitational potential using two-dimensional spline interpolation of Fourier coefficients in the meridional plane, as well as the ability to deal with rotation of the density profile and with multicomponent mass models. We compare several published methods for constructing composite axisymmetric disc--bulge--halo models and demonstrate that our code produces the models that are closest to equilibrium. We also apply it to create models of triaxial elliptical galaxies with cuspy density profiles and figure rotation, and find that such models can be found and are stable over many dynamical times in a wide range of pattern speeds and angular momenta, covering both slow- and fast-rotator classes. We then attempt to create models of strongly barred disc galaxies, using an analytic three-component potential, and find that it is not possible to make a stable dynamically self-consistent model for this density profile. Finally, we take snapshots of two N-body simulations of barred disc galaxies embedded in nearly-spherical haloes, and construct equilibrium models using only information on the density profile of the snapshots. We demonstrate that such reconstructed models are in near-stationary state, in contrast with the original N-body simulations, one of which displayed significant secular evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Suppression of stellar tidal disruption rates by anisotropic initial conditions

We compute the rates of capture of stars by supermassive black holes, using time-dependent Fokker--Planck equation with initial conditions that have a deficit of stars on low-angular-momentum orbits. One class of initial conditions has a gap in phase space created by a binary black hole, and the other has a globally tangentially-anisotropic velocity distribution. We find that for galactic nuclei that are younger than ~0.1 relaxation times, the flux of stars into the black hole is suppressed with respect to the steady-state value. This effect may substantially reduce the number of observable tidal disruption flares in galaxies with black hole masses M>10^7 Msun.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The final-parsec problem in the collisionless limit

A binary supermassive black hole loses energy via ejection of stars in a galactic nucleus, until emission of gravitational waves becomes strong enough to induce rapid coalescence. Evolution via the gravitational slingshot requires that stars be continuously supplied to the binary, and it is known that in spherical galaxies the reservoir of such stars is quickly depleted, leading to stalling of the binary at parsec-scale separations. Recent N-body simulations of galaxy mergers and isolated nonspherical galaxies suggest that this stalling may not occur in less idealized systems. However, it remains unclear to what degree these conclusions are affected by collisional relaxation, which is much stronger in the numerical simulations than in real galaxies. In this study, we present a novel Monte Carlo method that can efficiently deal with both collisional and collisionless dynamics, and with galaxy models having arbitrary shapes. We show that without relaxation, the final-parsec problem may be overcome only in triaxial galaxies. Axisymmetry is not enough, but even a moderate departure from axisymmetry is sufficient to keep the binary shrinking. We find that the binary hardening rate is always substantially lower than the maximum possible, "full-loss-cone" rate, and that it decreases with time, but that stellar-dynamical interactions are nevertheless able to drive the binary to coalescence on a timescale <=1 Gyr in any triaxial galaxy.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A new Monte Carlo method for dynamical evolution of non-spherical stellar systems

We have developed a novel Monte Carlo method for simulating the dynamical evolution of stellar systems in arbitrary geometry. The orbits of stars are followed in a smooth potential represented by a basis-set expansion and perturbed after each timestep using local velocity diffusion coefficients from the standard two-body relaxation theory. The potential and diffusion coefficients are updated after an interval of time that is a small fraction of the relaxation time, but may be longer than the dynamical time. Thus our approach is a bridge between the Spitzer's formulation of the Monte Carlo method and the temporally smoothed self-consistent field method. The primary advantages are the ability to follow the secular evolution of shape of the stellar system, and the possibility of scaling the amount of two-body relaxation to the necessary value, unrelated to the actual number of particles in the simulation. Possible future applications of this approach in galaxy dynamics include the problem of consumption of stars by a massive black hole in a non-spherical galactic nucleus, evolution of binary supermassive black holes, and the influence of chaos on the shape of galaxies, while for globular clusters it may be used for studying the influence of rotation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Rates of capture of stars by supermassive black holes in non-spherical galactic nuclei

We consider the problem of star consumption by supermassive black holes in non-spherical (axisymmetric, triaxial) galactic nuclei. We review the previous studies of the loss-cone problem and present a novel simulation method which allows to separate out the collisional (relaxation-related) and collisionless (related to non-conservation of angular momentum) processes and determine their relative importance for the capture rates in different geometries. We show that for black holes more massive than 10^7 Msun, the enhancement of capture rate in non-spherical galaxies is substantial, with even modest triaxiality being capable of keeping the capture rate at the level of few percent of black hole mass per Hubble time.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The final-parsec problem in non-spherical galaxies revisited

We consider the evolution of supermassive black hole binaries at the center of spherical, axisymmetric, and triaxial galaxies, using direct N-body integrations as well as analytic estimates. We find that the rates of binary hardening exhibit a significant N-dependence in all the models, at least for N in the investigated range of 10^5<=N<=10^6. Binary hardening rates are also substantially lower than would be expected if the binary loss cone remained full, as it would be if the orbits supplying stars to the binary were being efficiently replenished. The difference in binary hardening rates between the spherical and nonspherical models is less than a factor of two even in the simulations with the largest N. By studying the orbital populations of our models, we conclude that the rate of supply of stars to the binary via draining of centrophilic orbits is indeed expected to be much lower than the full-loss-cone rate, consistent with our simulations. We argue that the binary's evolution in the simulations is driven in roughly equal amounts by collisional and collisionless effects, even at the highest N-values currently accessible. While binary hardening rates would probably reach a limiting value for large N, our results suggest that we cannot approach that rate with currently available algorithms and computing hardware. The extrapolation of results from N-body simulations to real galaxies is therefore not straightforward, casting doubt on recent claims that triaxiality or axisymmetry alone are capable of solving the final-parsec problem in gas-free galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A new code for orbit analysis and Schwarzschild modelling of triaxial stellar systems

We review the methods used to study the orbital structure and chaotic properties of various galactic models and to construct self-consistent equilibrium solutions by Schwarzschild's orbit superposition technique. These methods are implemented in a new publicly available software tool, SMILE, which is intended to be a convenient and interactive instrument for studying a variety of 2D and 3D models, including arbitrary potentials represented by a basis-set expansion, a spherical-harmonic expansion with coefficients being smooth functions of radius (splines), or a set of fixed point masses. We also propose two new variants of Schwarzschild modelling, in which the density of each orbit is represented by the coefficients of the basis-set or spline spherical-harmonic expansion, and the orbit weights are assigned in such a way as to reproduce the coefficients of the underlying density model. We explore the accuracy of these general-purpose potential expansions and show that they may be efficiently used to approximate a wide range of analytic density models and serve as smooth representations of discrete particle sets (e.g. snapshots from an N-body simulation), for instance, for the purpose of orbit analysis of the snapshot. For the variants of Schwarzschild modelling, we use two test cases - a triaxial Dehnen model containing a central black hole, and a model re-created from an N-body snapshot obtained by a cold collapse. These tests demonstrate that all modelling approaches are capable of creating equilibrium models.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The loss cone problem in axisymmetric nuclei

We consider the problem of consumption of stars by a supermassive black hole (SBH) at the center of an axisymmetric galaxy. Inside the SBH sphere of influence, motion of stars in the mean field is regular and can be described analytically in terms of three integrals of motion: the energy E, the z-component of angular momentum L_z, and the secular Hamiltonian H. There exist two classes of orbits, tubes and saucers; saucers occupy the low-angular-momentum parts of phase space and their fraction is proportional to the degree of flattening of the nucleus. Perturbations due to gravitational encounters lead to diffusion of stars in integral space, which can be described using the Fokker-Planck equation. We calculate the diffusion coefficients and solve this equation in the two-dimensional phase space (L_z, H), for various values of the capture radius and the degree of flattening. Capture rates are found to be modestly higher than in the spherical case, up to a factor of a few, and most captures take place from saucer orbits. We also carry out a set of collisional N-body simulations to confirm the predictions of the Fokker-Planck models. We discuss the implications of our results for rates of tidal disruption and capture in the Milky Way and external galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Spin Evolution of Supermassive Black Holes and Galactic Nuclei

The spin angular momentum S of a supermassive black hole (SBH) precesses due to torques from orbiting stars, and the stellar orbits precess due to dragging of inertial frames by the spinning hole. We solve the coupled post-Newtonian equations describing the joint evolution of S and the stellar angular momenta Lj, j = 1...N in spherical, rotating nuclear star clusters. In the absence of gravitational interactions between the stars, two evolutionary modes are found: (1) nearly uniform precession of S about the total angular momentum vector of the system; (2) damped precession, leading, in less than one precessional period, to alignment of S with the angular momentum of the rotating cluster. Beyond a certain distance from the SBH, the time scale for angular momentum changes due to gravitational encounters between the stars is shorter than spin-orbit precession times. We present a model, based on the Ornstein-Uhlenbeck equation, for the stochastic evolution of star clusters due to gravitational encounters and use it to evaluate the evolution of S in nuclei where changes in the Lj are due to frame dragging close to the SBH and to encounters farther out. Long-term evolution in this case is well described as uniform precession of the SBH about the cluster's rotational axis, with an increasingly important stochastic contribution when SBH masses are small. Spin precessional periods are predicted to be strongly dependent on nuclear properties, but typical values are 10-100 Myr for low-mass SBHs in dense nuclei, 100 Myr - 10 Gyr for intermediate mass SBHs, and > 10 Gyr for the most massive SBHs. We compare the evolution of SBH spins in stellar nuclei to the case of torquing by an inclined, gaseous accretion disk.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Orbits Around Black Holes in Triaxial Nuclei

We discuss the properties of orbits within the influence sphere of a supermassive black hole (BH), in the case that the surrounding star cluster is nonaxisymmetric. There are four major orbit families; one of these, the pyramid orbits, have the interesting property that they can approach arbitrarily closely to the BH. We derive the orbit-averaged equations of motion and show that in the limit of weak triaxiality, the pyramid orbits are integrable: the motion consists of a two-dimensional libration of the major axis of the orbit about the short axis of the triaxial figure, with eccentricity varying as a function of the two orientation angles, and reaching unity at the corners. Because pyramid orbits occupy the lowest angular momentum regions of phase space, they compete with collisional loss cone repopulation and with resonant relaxation in supplying matter to BHs. General relativistic advance of the periapse dominates the precession for sufficiently eccentric orbits, and we show that relativity imposes an upper limit to the eccentricity: roughly the value at which the relativistic precession time is equal to the time for torques to change the angular momentum. We argue that this upper limit to the eccentricity should apply also to evolution driven by resonant relaxation, with potentially important consequences for the rate of extreme-mass-ratio inspirals in low-luminosity galaxies. In giant galaxies, we show that capture of stars on pyramid orbits can dominate the feeding of BHs, at least until such a time as the pyramid orbits are depleted; however this time can be of order a Hubble time.

0 citations0 reviewsNo code linkedOpen access