Source author record

Sergei Nayakshin

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

ResearcherUnclaimed source record

astro-ph.EP astro-ph.GA astro-ph.CO astro-ph.SR astro-ph.HE astro-ph astro-ph.IM

Catalog footprint

What is connected

58works

7topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2026arXiv

Disc fragmentation. II. Ejection of low mass Free Floating Planets from growing binary systems

Observations indicate that disc fragmentation due to Gravitational Instability (GI) is the likely origin of massive companions to stars, such as giant planets orbiting M-dwarf stars, Brown Dwarf (BD) companions to FGK stars, and binary stars with separations smaller than 100 au. Additionally, we have recently showed that disc fragmentation in young rapidly evolving binary systems ejects an abundant population of massive Jupiter-mass Free-Floating Planets (FFPs). In this model, a massive disc around an initially single protostar undergoes GI and hatches a number of fragments; the most massive oligarch grows by runaway accretion into the secondary star. As the system rearranges itself from a single to a binary star configuration, a dramatic "pincer movement" by the binary ejects planets through dynamical interactions with the stars. Here we propose that the same scenario applies to an even more abundant population of smaller FFPs discovered by the microlensing surveys. Although disc fragmentation is usually believed to form only massive objects, several pathways for forming small core-dominated planets at separations of tens of au exist. We present results from three complementary simulation approaches, all of which confirm planet ejection efficiency as high as 0.5 for secondaries more massive than $\sim 10$\% of the primary star mass. On the other hand, Jovian mass planets migrate through the region of tens of au too rapidly to eject planets from that region. We discuss implications of this scenario, concluding that microlensing FFPs may be the most convincing evidence yet that disc fragmentation forms planets much less massive than Jupiter.

preprint2022arXiv

ALMA constraints on assembly of Core Accretion planets

Resolved dust continuum and CO line ALMA imaging, and in some cases detection of H$α$ emission, hint that young massive planets are abundant at wide separations in protoplanetary discs. Here, we show how these observations can probe the runaway phase of planetary growth in the Core Accretion theory. Planets in this phase have the right range of masses to account for the predominantly moderate contrast gaps and rings seen in ALMA observations. However, we find that these planets gain mass and migrate inward very rapidly. As a result, the phase when they could produce gaps with properties similar to those observed is very short, i.e., $t_{\rm gap} \lesssim 0.1$~Myr, independently of the disc viscosity parameter. This would require many tens to hundreds of gas giant planets to be born per ALMA system, violating the available mass budget of solids in realistic discs. This also predicts preponderance of discs with very wide gaps or complete inner disc holes, which is not observed. We show that suppression of both planet accretion and migration by a factor of at least ten is a possible solution to these serious problems. Future population synthesis models of planet formation should aim to address both exoplanetary data of older discless planetary systems and ALMA discs with embedded planets in one framework.

preprint2022arXiv

Gap opening by planets in discs with magnetised winds

Planets open deep gaps in protoplanetary discs when their mass exceeds a gap opening mass, $M_{\rm gap}$. We use one- and two-dimensional simulations to study planet gap opening in discs with angular momentum transport powered by MHD disc winds. We parameterise the efficiency of the MHD disc wind angular momentum transport through a dimensionless parameter $α_{\rm dw}$, which is an analogue to the turbulent viscosity $α_{\rm v}$. We find that magnetised winds are much less efficient in counteracting planet tidal torques than turbulence is. For discs with astrophysically realistic values of $α_{\rm dw}$, $M_{\rm gap}$ is always determined by the residual disc turbulence, and is a factor of a few to ten smaller than usually obtained for viscous discs. We introduce a gap opening criterion applicable for any values of $α_{\rm v}$ and $α_{\rm dw}$ that may be useful for planet formation population synthesis. We show that in discs powered by magnetised winds, growing planets detach from the disc at planet masses below $\sim 0.1M_{\rm Jup}$ inside 10 AU. This promotes formation of super-Earth planets rather than gas giants in this region, in particular precluding formation of hot jupiters in situ. On larger scales, ALMA gap opening planet candidates may be less massive than currently believed. Future high-resolution observations with instruments such as the extended ALMA, ngVLA, and SKA are likely to show abundant narrow annular features at $R < 10$ AU due to ubiquitous super-Earth planets.

preprint2022arXiv

The growth and migration of massive planets under the influence of external photoevaporation

The formation of gas giant planets must occur during the first few Myr of a star's lifetime, when the protoplanetary disc still contains sufficient gas to be accreted onto the planetary core. The majority of protoplanetary discs are exposed to strong ultraviolet irradiation from nearby massive stars, which drives winds and depletes the mass budget for planet formation. It remains unclear to what degree external photoevaporation affects the formation of massive planets. In this work, we present a simple one dimensional model for the growth and migration of a massive planet under the influence of external FUV fields. We find that even moderate FUV fluxes $F_\mathrm{FUV}\gtrsim 100 \, G_0$ have a strong influence on planet mass and migration. By decreasing the local surface density and shutting off accretion onto the planet, external irradiation suppresses planet masses and halts migration early. The distribution of typical stellar birth environments can therefore produce an anti-correlation between semi-major axis and planet mass, which may explain the apparent decrease in planet occurrence rates at orbital periods $P_\mathrm{orb}\gtrsim 10^3$ days. Even moderate fluxes $F_\mathrm{FUV}$ strongly suppress giant planet formation and inward migration for any initial semi-major axis if the stellar host mass $M_*\lesssim 0.5\, M_\odot$, consistent with findings that massive planet occurrence is much lower around such stars. The outcomes of our prescription for external disc depletion show significant differences to the current approximation adopted in state-of-the-art population synthesis models, motivating future careful treatment of this important process.

preprint2020arXiv

The paradox of youth for ALMA planet candidates

Recent ALMA observations indicate that the majority of bright protoplanetary discs show signatures of young moderately massive planets. I show that this result is paradoxical. The planets should evolve away from their observed states by radial migration and gas accretion in about 1\% of the system age. These systems should then hatch tens of giant planets in their lifetime, and there should exist a very large population of bright planet-less discs; none of this is observationally supported. An alternative scenario, in which the population of bright ALMA discs is dominated by secondary discs recently rejuvenated by deposition of new gas, is proposed. The data are well explained if the gaseous mass of the discs is comparable to a Jovian planet mass, and they last a small fraction of a Million years. Self-disruptions of dusty gas giant protoplanets, previously predicted in the context of the Tidal Downsizing theory of planet formation, provide a suitable mechanism for such injections of new fuel, and yield disc and planet properties commensurate with ALMA observations. If this scenario is correct, then the secondary discs have gas-to-dust ratios considerably smaller than 100, and long look ALMA and NIR/optical observations of dimmer targets should uncover dusty, not yet disrupted, gas clumps with sizes of order an AU. Alternatively, secondary discs could originate from late external deposition of gas into the system, in which case we expect widespread signatures of warped outer discs that have not yet come into alignment with the planets.

preprint2020arXiv

TW Hya: an old protoplanetary disc revived by its planet

Dark rings with bright rims are the indirect signposts of planets embedded in protoplanetary discs. In a recent first, an azimuthally elongated AU-scale blob, possibly a planet, was resolved with ALMA in TW Hya. The blob is at the edge of a cliff-like rollover in the dust disc rather than inside a dark ring. Here we build time-dependent models of TW Hya disc. We find that the classical paradigm cannot account for the morphology of the disc and the blob. We propose that ALMA-discovered blob hides a Neptune mass planet losing gas and dust. We show that radial drift of mm-sized dust particles naturally explains why the blob is located on the edge of the dust disc. Dust particles leaving the planet perform a characteristic U-turn relative to it, producing an azimuthally elongated blob-like emission feature. This scenario also explains why a 10 Myr old disc is so bright in dust continuum. Two scenarios for the dust-losing planet are presented. In the first, a dusty pre-runaway gas envelope of about 40 Earth mass Core Accretion planet is disrupted, e.g., as a result of a catastrophic encounter. In the second, a massive dusty pre-collapse gas giant planet formed by Gravitational Instability is disrupted by the energy released in its massive core. Future modelling may discriminate between these scenarios and allow us to study planet formation in an entirely new way -- by analysing the flows of dust and gas recently belonging to planets, informing us about the structure of pre-disruption planetary envelopes.

preprint2016arXiv

A desert of gas giant planets beyond tens of au

Direct imaging observations constrain the fraction of stars orbited by gas giant planets with separations greater than 10 au to about 0.01 only. This is widely believed to indicate that massive protoplanetary discs rarely fragment on planetary mass objects. I use numerical simulations of gas clumps embedded in massive gas discs to show that these observations are consistent with $\sim 0.2 - 10$ planetary mass clumps per star being born in young gravitationally unstable discs. A trio of processes -- rapid clump migration, tidal disruption and runaway gas accretion -- destroys or transforms all of the simulated clumps into other objects, resulting in a desert of gas giants beyond separation of approximately 10 au. The cooling rate of the disc controls which of the three processes is dominant. For cooling rates faster than a few local dynamical times, clumps always grow rapidly and become massive brown dwarfs or low mass stars. For longer cooling times, post-collapse (high density) planets migrate inward to $\sim 10-20$ au where they open a gap in the disc and then continue to migrate inward much less rapidly. Pre-collapse (low density) planets are tidally disrupted and may leave massive solid cores behind. Gas giant planets observed inside the desert, such as those in HR 8799, must have followed an unusual evolutionary path, e.g., their host disc being dispersed in a catastrophic fashion.

preprint2016arXiv

Planets, debris and their host metallicity correlations

Recent observations of debris discs, believed to be made up of remnant planetesimals, brought a number of surprises. Debris disc presence does not correlate with the host star's metallicity, and may anti-correlate with the presence of gas giant planets. These observations contradict both assumptions and predictions of the highly successful Core Accretion model of planet formation. Here we explore predictions of the alternative Tidal Downsizing (TD) scenario of planet formation. In TD, small planets and planetesimal debris is made only when gas fragments, predecessors of giant planets, are tidally disrupted. We show that these disruptions are rare in discs around high metallicity stars but release more debris per disruption than their low [M/H] analogs. This predicts no simple relation between debris disc presence and host star's [M/H], as observed. A detected gas giant planet implies in TD that its predecessor fragment was not disputed, potentially explaining why DDs are less likely to be found around stars with gas giants. Less massive planets should correlate with DD presence, and sub-Saturn planets ($M_{p} \sim 50 \M_{\oplus}$) should correlate with DD presence stronger than sub-Neptunes ($M_{p} \leq 15 \M_{\oplus}$). These predicted planet-DD correlations will be diluted and weakened in observations by planetary systems' long term evolution and multi-fragment effects neglected here. Finally, although presently difficult to observe, DDs around M dwarf stars should be more prevalent than around Solar type stars.

preprint2015arXiv

A simple way to improve AGN feedback prescription in SPH simulations

AGN feedback is an important ingredient in galaxy evolution, however its treatment in numerical simulations is necessarily approximate, requiring subgrid prescriptions due to the dynamical range involved in the calculations. We present a suite of SPH simulations designed to showcase the importance of the choice of a particular subgrid prescription for AGN feedback. We concentrate on two approaches to treating wide-angle AGN outflows: thermal feedback, where thermal and kinetic energy is injected into the gas surrounding the SMBH particle, and virtual particle feedback, where energy is carried by tracer particles radially away from the AGN. We show that the latter model produces a far more complex structure around the SMBH, which we argue is a more physically correct outcome. We suggest a simple improvement to the thermal feedback model - injecting the energy into a cone, rather than spherically symmetrically - and show that this markedly improves the agreement between the two prescriptions, without requiring any noticeable increase in the computational cost of the simulation.

preprint2015arXiv

The resolution bias: low resolution feedback simulations are better at destroying galaxies

Feedback from super-massive black holes (SMBHs) is thought to play a key role in regulating the growth of host galaxies. Cosmological and galaxy formation simulations using smoothed particle hydrodynamics (SPH), which usually use a fixed mass for SPH particles, often employ the same sub-grid Active galactic nuclei (AGN) feedback prescription across a range of resolutions. It is thus important to ask how the impact of the simulated AGN feedback on a galaxy changes when only the numerical resolution (the SPH particle mass) changes. We present a suite of simulations modelling the interaction of an AGN outflow with the ambient turbulent and clumpy interstellar medium (ISM) in the inner part of the host galaxy at a range of mass resolutions. We find that, with other things being equal, degrading the resolution leads to feedback becoming more efficient at clearing out all gas in its path. For the simulations presented here, the difference in the mass of the gas ejected by AGN feedback varies by more than a factor of ten between our highest and lowest resolution simulations. This happens because feedback-resistant high density clumps are washed out at low effective resolutions. We also find that changes in numerical resolution lead to undesirable artifacts in how the AGN feedback affects the AGN immediate environment.

preprint2015arXiv

The role of feedback in accretion on Low Luminosity AGN: Sgr A* case study

We present numerical models of the gas dynamics in the inner parsec of the Galactic centre. We follow the gas from its origin as stellar winds of several observed young massive stars, until it is either captured by the central black hole, or leaves the system. Unlike our previous models, we include an outflow from the inner accretion flow. Two different kinds of outflows are modelled: (i) an instantaneous-response feedback mode, in which the outflow rate is directly proportional to the current black hole gas capture rate; and (ii) an outburst mode, which is stronger but lasts for a limited time. The latter situation may be particularly relevant to Sgr A*, since there is evidence that Sgr A* was much brighter in the recent past. We find that both types of outflow perturb the gas dynamics near the Bondi radius and the black hole capture rate significantly. The effects persist longer than the outflow itself. We also compare the effects of spherically symmetric and collimated outflows, and find that the latter are far less efficient in transferring its energy to the surrounding gas near the capture radius. Our results imply that accretion feedback is important for non-radiative accretion flows not only within but also outside the capture radius. Steady-state Bondi accretion rate estimates that do not account for feedback outflows over-predict not only the accretion rate onto the black hole but also the capture rate at the Bondi radius itself. Finally, the steady-state assumption under which non-radiative flows have been routinely studied in the literature may have to be abandoned if accretion feedback is bursty in nature.

preprint2015arXiv

Tidal Downsizing model. II. Planet-metallicity correlations

Core Accretion (CA), the de-facto accepted theory of planet formation, requires formation of massive solid cores as a prerequisite for assembly of gas giant planets. The observed metallicity correlations of exoplanets are puzzling in the context of CA. While gas giant planets are found preferentially around metal-rich host stars, planets smaller than Neptune orbit hosts with a wide range of metallicities. We propose an alternative interpretation of these observations in the framework of a recently developed planet formation hypothesis called Tidal Downsizing (TD). We perform population synthesis calculations based on TD, and find that the connection between the populations of the gas giant and the smaller solid-core dominated planets is non linear and not even monotonic. While gas giant planets formed in the simulations in the inner few AU region follow a strong positive correlation with the host star metallicity, the smaller planets do not. The simulated population of these smaller planets shows a shallow peak in their formation efficiency at around the Solar metallicity. This result is driven by the fact that at low metallicities the solid core's growth is damped by the scarcity of metals, whereas at high metallicities the fragments within which the cores grow contract too quickly, cutting the core's growth time window short. Finally, simulated giant gas planets do not show a strong host star metallicity preference at large separations, which may explain why one of the best known directly imaged gas giant planet systems, HR 8799, is metal poor.

preprint2015arXiv

Tidal Downsizing Model. III. Planets from sub-Earths to Brown Dwarfs: structure and metallicity preferences

We present population synthesis calculations of the Tidal Downsizing (TD) hypothesis for planet formation. Our models address the following observations: (i) most abundant planets being Super Earths; (ii) cores more massive than $\sim 5-15 M_\oplus$ are enveloped by massive atmospheres; (iii) the frequency of occurrence of close-in gas giant planets correlates strongly with metallicity of the host star; (iv) no such correlation is found for sub-Neptune planets; (v) presence of massive cores in giant planets; (vi) gas giant planets are over-abundant in metals compared to their host stars; (vii) this over-abundance decreases with planet's mass; (viii) a deep valley in the planet mass function between masses of $\sim 10-20 M_\oplus$ and $\sim 100 M_\oplus$. A number of observational predictions distinguish the model from Core Accretion: (a) composition of the massive cores is always dominated by rocks not ices; (b) the core mass function is smooth with no minimum at $\sim 3 M_\oplus$ and has no ice-dominated cores; (c) gas giants beyond 10 AU are insensitive to the host star metallicity; (d) objects more massive than $\sim 10 M_{\rm Jup}$ do not correlate or even anti-correlate with metallicity. The latter prediction is consistent with observations of low mass stellar companions. TD can also explain formation of planets in close binary systems. TD model is a viable alternative to the Core Accretion scenario in explaining many features of the observed population of exoplanets.

preprint2015arXiv

Tidal Downsizing Model. IV. Destructive feedback in planets

I argue that feedback is as important to formation of planets as it is to formation of stars and galaxies. Energy released by massive solid cores puffs up pre-collapse gas giant planets, making them vulnerable to tidal disruptions by their host stars. I find that feedback is the ultimate reason for some of the most robust properties of the observed exoplanet populations: the rarity of gas giants at all separations from $\sim 0.1$ to $\sim 100$~AU, the abundance of $\sim 10 M_\oplus$ cores but dearth of planets more massive than $\sim 20 M_\oplus$. Feedback effects can also explain (i) rapid assembly of massive cores at large separations as needed for Uranus, Neptune and the suspected HL Tau planets; (ii) the small core in Jupiter yet large cores in Uranus and Neptune; (iii) the existence of rare "metal monster" planets such as CoRoT-20b, a gas giant made of heavy elements by up to $\sim 50$\%.

preprint2014arXiv

Black hole feedback in a multiphase interstellar medium

Ultrafast outflows (UFOs) from supermassive black holes (SMBHs) are thought to regulate the growth of SMBHs and host galaxies, resulting in a number of observational correlations. We present high-resolution numerical simulations of the impact of a thermalized UFO on the ambient gas in the inner part of the host galaxy. Our results depend strongly on whether the gas is homogeneous or clumpy. In the former case all of the ambient gas is driven outward rapidly as expected based on commonly used energy budget arguments, while in the latter the flows of mass and energy decouple. Carrying most of the energy, the shocked UFO escapes from the bulge via paths of least resistance, taking with it only the low-density phase of the host. Most of the mass is however in the high-density phase, and is affected by the UFO much less strongly, and may even continue to flow inwards. We suggest that the UFO energy leakage through the pores in the multiphase interstellar medium (ISM) may explain why observed SMBHs are so massive despite their overwhelmingly large energy production rates. The multiphase ISM effects reported here are probably under-resolved in cosmological simulations but may be included in prescriptions for active galactic nuclei feedback in future simulations and in semi-analytical models.

preprint2014arXiv

Core-assisted gas capture instability: a new mode of giant planet formation by gravitationally unstable discs

Giant planet formation in the core accretion (CA) paradigm is predicated by the formation of a core, assembled by the coagulation of grains and later by planetesimals within a protoplanetary disc. In contrast, in the disc instability paradigm, giant planet formation is believed to be independent of core formation: massive self-gravitating gas fragments cool radiatively and collapse as a whole. We show that giant planet formation in the disc instability model may be also enhanced by core formation for reasons physically very similar to the CA paradigm. In the model explored here, efficient grain sedimentation within an initial fragment (rather than the disc) leads to the formation of a core composed of heavy elements. We find that massive atmospheres form around cores and undergo collapse as a critical core mass is exceeded, analogous to CA theory. The critical mass of the core to initiate such a collapse depends on the fragment mass and metallicity, as well as core luminosity, but ranges from less than 1 to as much as $\sim80$ Earth masses. We therefore suggest that there are two channels for the collapse of a gaseous fragment to planetary scales within the disc instability model: (i) H$_2$ dissociative collapse of the entire gaseous clump, and (ii) core-assisted gas capture, as presented here. We suggest that the first of these two is favoured in metal-poor environments and for fragments more massive than $\sim 5-10$ Jupiter masses, whereas the second is favored in metal-rich environments and fragments of lower mass. [Abridged]

preprint2014arXiv

Energy- and momentum-conserving AGN feedback outflows

It is usually assumed that outflows from luminous AGN are either in the energy-conserving (non-radiative) or in the momentum-conserving (radiative) regime. We show that in a non-spherical geometry the effects of both regimes may manifest at the same time, and that it is the momentum of the outflow that sets the $M_{\rm BH}-σ$ relation. Considering an initially elliptical distribution of gas in the host galaxy, we show that a non-radiative outflow opens up a wide ``escape route'' over the paths of least resistance. Most of the outflow energy escapes in that direction. At the same time, in the directions of higher resistance, the ambient gas is affected mainly by the incident momentum from the outflow. Quenching SMBH growth requires quenching gas delivery along the paths of highest resistance, and therefore, it is the momentum of the outflow that limits the black hole growth. We present an analytical argument showing that such energy-conserving feedback bubbles driving leaky ambient shells will terminate SMBH growth once its mass reaches roughly the $M_σ$ mass derived earlier by King (2003) for momentum-conserving AGN outflows. Our simulations also have potentially important implications for observations of AGN jet feedback and starburst galaxy feedback. The collimation of the wide angle AGN outflow away from the symmetry plane, as found in our simulations, indicates that credit for work done by such outflows may sometimes be mistakenly given to AGN jets or star formation feedback since wide angle $v \sim 0.1 c$ outflows are harder to observe and the phase when they are present may be short.

preprint2014arXiv

Metal loading of giant gas planets

One of many challenges in forming giant gas planets via Gravitational disc Instability model (GI) is an inefficient radiative cooling of the pre-collapse fragments. Since fragment contraction times are as long at $10^5 -10^7$ years, the fragments may be tidally destroyed sooner than they contract into gas giant planets. Here we explore the role of "pebble accretion" onto the pre-collapse giant planets and find an unexpected result. Despite larger dust opacity at higher metallicities, addition of metals actually accelerates -- rather than slows down -- collapse of high opacity, relatively low mass giant gas planets ($M_p$ below a few Jupiter masses). A simple analytical theory that explains this result exactly in idealised simplified cases is presented. The theory shows that planets with the central temperature in the range between $\sim$ 1000 to 2000K are especially sensitive to pebble accretion: addition of just $\sim 5$ to 10 % of metals by weight is sufficient to cause their collapse. These results show that dust grain physics and dynamics is essential for an accurate modelling of self-gravitating disc fragments and their near environments in the outer massive and cold protoplanetary discs.

preprint2014arXiv

Positive metallicity correlation for coreless giant planets

Frequency of detected giant planets is observed to increase rapidly with metallicity of the host star. This is usually interpreted as evidence in support of the Core Accretion (CA) theory, which assembles giant planets as a result of formation of a massive solid core. A strong positive planet-metallicity correlation for giant planets formed in the framework of Gravitational disc Instability (GI) model is found here. The key novelty of this work is "pebble accretion" onto GI fragments which has been recently demonstrated to accelerate contraction of GI fragments. Driven closer to the star by the inward migration, only the fragments that accrete metals rapidly enough collapse and survive the otherwise imminent tidal disruption. The survival fraction of simulated planets correlates strongly with the metallicity of the host star, as observed.

preprint2014arXiv

Tidal Downsizing model. I. Numerical methods: saving giant planets from tidal disruptions

Tidal Downsizing (TD) is a recently developed planet formation theory that supplements the classical Gravitational disc Instability (GI) model with planet migration inward and tidal disruptions of GI fragments in the inner regions of the disc. Numerical methods for a detailed population synthesis of TD planets are presented here. As an example application, the conditions under which GI fragments collapse faster than they migrate into the inner $a\sim$ few AU disc are considered. It is found that most gas fragments are tidally or thermally disrupted unless (a) their opacity is $\sim 3$ orders of magnitude less than the interstellar dust opacity at metallicities typical of the observed giant planets, or (b) the opacity is high but the fragments accrete large dust grains (pebbles) from the disc. Case (a) models produce very low mass solid cores ($M_{\rm core} < 0.1$ Earth masses) and follow a negative correlation of giant planet frequency with host star metallicity. In contrast, case (b) models produce massive solid cores, correlate positively with host metallicity and explain naturally while giant gas planets are over-abundant in metals.

preprint2013arXiv

A link between feedback outflows and satellite galaxy suppression

We suggest a direct link between the two "missing" baryon problems of contemporary galaxy formation theory: (1) that large galaxies are known to contain too little gas and stars and (2) that too few dwarf satellite galaxies are observed around large galaxies compared with cosmological simulations. The former can be explained by invoking some energetic process -- most likely AGN or star formation feedback -- which expels to infinity a significant fraction of the gas initially present in the proto-galaxy, while the latter problem is usually explained by star formation feedback inside the dwarf or tidal and ram pressure stripping of the gas from the satellite galaxy by its parent. Here we point out that the host galaxy "missing" baryons, if indeed ejected at velocities of hundreds to a thousand km s$^{-1}$, must also affect smaller satellite galaxies by stripping or shocking the gas there. We estimate the fraction of gas removed from the satellites as a function of the satellite galaxy's properties. Applying these results to a Milky Way like dark matter halo, we find that this singular shock ram pressure stripping event may be quite efficient in removing the gas from the satellites provided that they are close enough. We also use the orbital and mass modelling data for eight Galactic dwarf spheroidal (dSph) satellites, and find that it is likely that many of them have been affected by the Galactic outflow. We point out that galactic outflows of the host may also trigger a starburst in the satellite galaxies by over-pressuring their gas discs. This process may be responsible for the formation of the globular clusters observed in some of the Milky Way's dSphs (e.g. the Fornax and Sagittarius dSphs).

preprint2013arXiv

AGN outflows trigger starbursts in gas-rich galaxies

Recent well resolved numerical simulations of AGN feedback have shown that its effects on the host galaxy may be not only negative but also positive. In the late gas poor phase, AGN feedback blows the gas away and terminates star formation. However, in the gas-rich phase(s), AGN outflows trigger star formation by over-compressing cold dense gas and thus provide positive feedback on their hosts. In this paper we study this AGN-triggered starburst effect. We show that star formation rate in the burst increases until the star formation feedback counteracts locally the AGN outflow compression. Globally, this predicts a strong nearly linear statistical correlation between the AGN and starburst bolometric luminosities in disc galaxies, L_* \propto L_{AGN}^{5/6}. The correlation is statistical only because AGN activity may fluctuate on short time scales (as short as tens of years), and because AGN may turn off but its effects on the host may continue to last until the AGN-driven outflow leaves the host, which may be up to 10 times longer than the duration of the AGN activity. The coefficient in front of this relation depends on the clumpiness and morphology of the cold gas in the galaxy. A "maximum starburst" takes place in am azimuthally uniform gas disc, for which we derive an upper limit of L_* \sim 50 times larger than L_{AGN} for typical quasars. For more clumpy and/or compact cold gas distributions, the starburst luminosity decreases. We also suggest that similar AGN-triggerred starbursts are possible in hosts of all geometries, including during galaxy mergers, provided the AGN is activated. Finally, we note that due to the short duration of the AGN activity phase the accelerating influence of AGN on starbursts may be much more common than observations of simultaneous AGN and starbursts would suggest.

preprint2013arXiv

Hydrogen-loosing planets in transition discs around young protostars

We point out that protoplanets created in the framework of the Tidal Downsizing (TD) theory for planet formation play a very important role for the evolution of accretion discs hosting them. Since all TD protoplanets are initially as massive as $\sim 10$ Jupiter masses, they are able to open very deep gaps in their discs, and even completely isolate the inner disc flows from the outer ones. Furthermore, in contrast to other planet formation theories, TD protoplanets are mass donors for their protostars. One potentially observable signature of planets being devoured by their protostars are FU Ori like outbursts, and episodic protostar accretion more generally, as discussed by a number of authors recently. Here we explore another observational implication of TD hypothesis: dust poor inner accretion flows, which we believe may be relevant to some of the observed mm-bright transitional discs around protostars. In our model, a massive protoplanet interrupts the flow of the outer dust-rich disc on its protostar, and at the same time looses a part of its dust-poor envelope into the inner disc. This then powers the observed gas-but-no-dust accretion onto the star. Upon a more detailed investigation, we find that this scenario is quite natural for young massive discs but is less so for older discs, e.g., those whose self-gravitating phase has terminated a fraction of a Million year or more ago. This stems from the fact that TD protoplanets of such an age should have contracted significantly, and so are unlikely to loose much mass. Therefore, we conclude that either (i) the population of "transition discs" with large holes and dust-poor accretion is much younger than generally believed; or (ii) there is a poorly understood stage for late removal of dust-poor envelopes from TD planets; (iii) another explanation for the observations is correct.

preprint2013arXiv

Inverse Compton X-ray signature of AGN feedback

Bright AGN frequently show ultra-fast outflows (UFOs) with outflow velocities vout ~0.1c. These outflows may be the source of AGN feedback on their host galaxies sought by galaxy formation modellers. The exact effect of the outflows on the ambient galaxy gas strongly depends on whether the shocked UFOs cool rapidly or not. This in turn depends on whether the shocked electrons share the same temperature as ions (one temperature regime; 1T) or decouple (2T), as has been recently suggested. Here we calculate the Inverse Compton spectrum emitted by such shocks, finding a broad feature potentially detectable either in mid-to-high energy X-rays (1T case) or only in the soft X-rays (2T). We argue that current observations of AGN do not seem to show evidence for the 1T component. The limits on the 2T emission are far weaker, and in fact it is possible that the observed soft X-ray excess of AGN is partially or fully due to the 2T shock emission. This suggests that UFOs are in the energy-driven regime outside the central few pc, and must pump considerable amounts of not only momentum but also energy into the ambient gas. We encourage X-ray observers to look for the Inverse Compton components calculated here in order to constrain AGN feedback models further.

preprint2013arXiv

Outflows of stars due to quasar feedback

Quasar feedback outflows are commonly invoked to drive gas out of galaxies in the early gas-rich epoch to terminate growth of galaxies. Here we present simulations that show that AGN feedback may drive not only gas but also stars out of their host galaxies under certain conditions. The mechanics of this process is as following: (1) AGN-driven outflows accelerate and compress gas filling the host galaxy; (2) the accelerated dense shells become gravitationally unstable and form stars on radial trajectories. For the spherically symmetric initial conditions explored here, the black hole needs to exceed the host's M_sigma mass by a factor of a few to accelerate the shells and the new stars to escape velocities. We discuss potential implications of these effects for the host galaxies: (i) radial mixing of bulge stars with the rest of the host; (ii) contribution of quasar outflows to galactic fountains as sources of high-velocity clouds; (iii) wholesale ejection of hyper velocity stars out of their hosts, giving rise to type II supernovae on galactic outskirts, and contributing to reionization and metal enrichment of the Universe; (iv) bulge erosion and even complete destruction in extreme cases resulting in overweight or bulgeless SMBHs.

preprint2013arXiv

Radiative feedback from protoplanets in self-gravitating protoplanetary discs

It is well known that massive protoplanetary disc are gravitationally unstable beyond tens of AU from their parent star. The eventual fate of the self-gravitating gas clumps born in the disc is currently not understood, although the range of uncertainty is well known. If clumps migrate inward rapidly, they are tidally disrupted, which may leave behind giant or terrestrial like planets. On the other hand, if clumps migrate less rapidly, they tend to accrete gas, becoming proto brown dwarfs or low mass companions to the parent star. Here we argue that radiative feedback of contracting clumps (protoplanets) on their discs is an important effect that has been overlooked in previous calculations. We show analytically that temperature in clump's vicinity may be high enough to support a quasi-static atmosphere if the clump mass is below a critical value, $M_{\rm cr} \sim 6$ Jupiter masses ($M_J$). This may arrest further gas accretion onto the clump and thus promote formation of planets rather than low mass companions. We use numerical simulations to evaluate these analytical conclusions by studying migration and accretion of gas clumps as a function of their initial mass, $M_i$. Simulations neglecting the radiative preheating effect show that gas clumps with mass less than $\sim 2 M_J$ migrate inward rapidly; more massive clumps result in low mass companions. In contrast, simulations that include radiative preheating from the clump show that clumps as massive as $8 M_J$ migrate inward rapidly and end up tidally disrupted in the inner disc. We conclude that, with other parameters being equal, previous simulations neglecting radiative feedback from self-gravitating clumps over-estimated the population of brown dwarfs and low mass stellar companions and under-estimated the population of planets.

preprint2013arXiv

Two-phase model for Black Hole feeding and feedback

We study effects of AGN feedback outflows on multi-phase inter stellar medium (ISM) of the host galaxy. We argue that SMBH growth is dominated by accretion of dense cold clumps and filaments. AGN feedback outflows overtake the cold medium, compress it, and trigger a powerful starburst -- a positive AGN feedback. This predicts a statistical correlation between AGN luminosity and star formation rate at high luminosities. Most of the outflow's kinetic energy escapes from the bulge via low density voids. The cold phase is pushed outward only by the ram pressure (momentum) of the outflow. The combination of the negative and positive forms of AGN feedback leads to an $M-σ$ relation similar to the result of King (2003). Due to porosity of cold ISM in the bulge, SMBH influence on the low density medium of the host galaxy is significant even for SMBH well below the $M-σ$ mass. The role of SMBH feedback in our model evolves in space and time with the ISM structure. In the early gas rich phase, SMBH accelerates star formation in the bulge. During later gas poor (red-and-dead) phases, SMBH feedback is mostly negative everywhere due to scarcity of the cold ISM.

preprint2012arXiv

An alternative origin for debris rings of planetesimals

Core Accretion, the most widely accepted scenario for planet formation, postulates existence of km-sized solid bodies, called planetesimals, arranged in a razor-thin disc in the earliest phases of planet formation. In the Tidal Downsizing hypothesis, an alternative scenario for formation of planets, grain growth, sedimentation and formation of planetary cores occur inside dense and massive gas clumps formed in the outer cold disc by gravitational instability. As a clump migrates inward, tidal forces of the star remove all or most of the gas from the clump, downsizing it to a planetary mass body. Here we argue that such a clump may form not only the planetary core but also numerous smaller bodies. As an example, we consider the simplest case of bodies on circular orbits around the planetary core in the centre of the gas clump. Bodies smaller than 1 km suffer a strong enough aerodynamic drag, spiral in and accrete onto the solid core rapidly; bodies in the planetesimal size range lose their centrifugal support very slowly. We find that planetesimals orbiting the protoplanetary core closely remain gravitationally bound to it; these may be relevant to formation of satellites of giant planets. Planetesimals on more distant orbits within the host clump are unbound from the protoplanet and are set on mildly eccentric heliocentric orbits, generically forming wide rings. These may correspond to debris discs around main sequence stars and the Kuiper belt in the Solar System. For the latter in particular, our hypothesis naturally explains the observed sharp outer edge and the "mass deficit" of the Kuiper belt.

preprint2012arXiv

Fermi Bubbles in the Milky Way: the closest AGN feedback laboratory courtesy of Sgr A*?

Deposition of a massive ($10^4$ to $10^5 \msun$) giant molecular cloud (GMC) into the inner parsec of the Galaxy is widely believed to explain the origin of over a hundred unusually massive young stars born there $\sim 6$ Myr ago. An unknown fraction of that gas could have been accreted by Sgr A*, the supermassive black hole (SMBH) of the Milky Way. It has been recently suggested that two observed $γ$-ray-emitting bubbles emanating from the very center of our Galaxy were inflated by this putative activity of Sgr A*. We run a suite of numerical simulations to test whether the observed morphology of the bubbles could be due to the collimation of a wide angle outflow from Sgr A* by the disc-like Central Molecular Zone (CMZ), a well known massive repository of molecular gas in the central $\sim 200$ pc. We find that an Eddington-limited outburst of Sgr A* lasting $\simeq 1$ Myr is required to reproduce the morphology of the {\it Fermi} bubbles, suggesting that the GMC mass was $\sim 10^5 \msun$ and it was mainly accreted by Sgr A* rather than used to make stars. We also find that the outflow from Sgr A* enforces strong angular momentum mixing in the CMZ disc, robustly sculpting it into a much narrower structure -- a ring -- perhaps synonymous with the recently reported "Herschel ring". In addition, we find that Sgr A* outflow is likely to have induced formation of massive star-forming GMCs in the CMZ. In this scenario, the Arches and Quintuplet clusters, the two observed young star clusters in the central tens of parsecs of the Galaxy, and also GMCs such as Sgr B2, owe their existence to the recent Sgr A* activity.

preprint2012arXiv

Modelling supermassive black hole growth: towards an improved sub-grid prescription

Accretion onto supermassive black holes (SMBHs) in galaxy formation simulations is frequently modelled by the Bondi-Hoyle formalism. Here we examine the validity of this approach analytically and numerically. We argue that the character of the flow where one evaluates the gas properties is unlikely to satisfy the simple Bondi-Hoyle model. Only in the specific case of hot virialised gas with zero angular momentum and negligible radiative cooling is the Bondi-Hoyle solution relevant. In the opposite extreme, where the gas is in a state of free-fall at the evaluation radius due to efficient cooling and the dominant gravity of the surrounding halo, the Bondi-Hoyle formalism can be erroneous by orders of magnitude in either direction. This may impose artificial trends with halo mass in cosmological simulations by being wrong by different factors for different halo masses. We propose an expression for the sub-grid accretion rate which interpolates between the free-fall regime and the Bondi-Hoyle regime, therefore taking account of the contribution of the halo to the gas dynamics.

preprint2012arXiv

Quasar feedback: accelerated star formation and chaotic accretion

Growing Supermassive Black Holes (SMBH) are believed to influence their parent galaxies in a negative way, terminating their growth by ejecting gas out before it could turn into stars. Here we present some of the most sophisticated SMBH feedback simulations to date showing that quasar's effects on galaxies are not always negative. We find that when the ambient shocked gas cools rapidly, the shocked gas is compressed into thin cold dense shells, filaments and clumps. Driving these high density features out is much more difficult than analytical models predict since dense filaments are resilient to the feedback. However, in this regime quasars have another way of affecting the host -- by triggering a massive star formation burst in the cold gas by over-pressurising it. Under these conditions SMBHs actually accelerate star formation in the host, having a positive rather than negative effect on their host galaxies. The relationship between SMBH and galaxies is thus even more complex and symbiotic than currently believed. We also suggest that the instabilities found here may encourage the chaotic AGN feeding mode.

preprint2012arXiv

The observed $M - σ$ relations imply that SMBHs grow by cold chaotic accretion

We argue that current observations of $M - σ$ relations for galaxies can be used to constrain theories of super-massive black holes (SMBH) feeding. In particular, assuming that SMBH mass is limited only by the feedback on the gas that feeds it, we show that SMBHs fed via a planar galaxy scale gas flow, such as a disc or a bar, should be much more massive than their counterparts fed by quasi-spherical inflows. This follows from the relative inefficiency of AGN feedback on a flattened inflow. We find that even under the most optimistic conditions for SMBH feedback on flattened inflows, the mass at which the SMBH expels the gas disc and terminates its own growth is a factor of several higher than the one established for quasi-spherical inflows. Any beaming of feedback away from the disc and any disc self-shadowing strengthens this result further. Contrary to this theoretical expectation, recent observations have shown that SMBH in pseudobulge galaxies (which are associated with barred galaxies) are typically under- rather than over-massive when compared with their classical bulge counterparts at a fixed value of $σ$. We conclude from this that SMBHs are not fed by large (100 pc to many kpc) scale gas discs or bars, most likely because such planar flows are turned into stars too efficiently to allow any SMBH growth. Based on this and other related observational evidence, we argue that most SMBHs grow by chaotic accretion of gas clouds with a small and nearly randomly distributed direction of angular momentum.

preprint2012arXiv

What's in a Fermi Bubble: a quasar episode in the Galactic centre

Fermi bubbles, the recently observed giant (~10 kpc high) gamma-ray emitting lobes on either side of our Galaxy (Su et al. 2010), appear morphologically connected to the Galactic center, and thus offer a chance to test several models of supermassive black hole (SMBH) evolution, feedback and relation with their host galaxies. We use a physical feedback model (King 2003, 2010) and novel numerical techniques (Nayakshin et al. 2009) to simulate a short burst of activity in Sgr A*, the central SMBH of the Milky Way, ~6 Myr ago, temporally coincident with a star formation event in the central parsec. We are able to reproduce the bubble morphology and energetics both analytically (Zubovas et al. 2011) and numerically (Zubovas & Nayakshin, in prep). These results provide strong support to the model, which was also used to simulate more extreme environments (Nayakshin & Power 2010).

preprint2011arXiv

A numerical simulation of a "super-Earth" core delivery from ~ 100 AU to ~ 8 AU

We use SPH simulations with an approximate radiative cooling prescription to model evolution of a massive and large ($\sim 100$ AU) very young protoplanetary disc. We also model dust growth and gas-grain dynamics with a second fluid approach. It is found that the disc fragments onto a large number of $\sim 10$ Jupiter mass clumps that cool and contract slowly. Some of the clumps evolve onto eccentric orbits delivering them into the inner tens of AU, where they are disrupted by tidal forces from the star. Dust grows and sediments inside the clumps, displaying a very strong segregation, with the largest particles forming dense cores in the centres. The density of the dust cores may exceed that of the gas and is limited only by the numerical constraints, indicating that these cores should collapse into rocky planetary cores. One particular giant planet embryo migrates inward close enough to be disrupted at about 10 AU, leaving a self-bound solid core of about 7.5 $\mearth$ mass on a low eccentricity orbit at a radius of $\sim$ 8 AU. These simulations support the recent suggestions that terrestrial and giant planets may be the remnants of tidally disrupted giant planet embryos.

preprint2011arXiv

A relativistic jetted outburst from a massive black hole fed by a tidally disrupted star

While gas accretion onto some massive black holes (MBHs) at the centers of galaxies actively powers luminous emission, the vast majority of MBHs are considered dormant. Occasionally, a star passing too near a MBH is torn apart by gravitational forces, leading to a bright panchromatic tidal disruption flare (TDF). While the high-energy transient Swift J164449.3+573451 ("Sw 1644+57") initially displayed none of the theoretically anticipated (nor previously observed) TDF characteristics, we show that the observations (Levan et al. 2011) suggest a sudden accretion event onto a central MBH of mass ~10^6-10^7 solar masses. We find evidence for a mildly relativistic outflow, jet collimation, and a spectrum characterized by synchrotron and inverse Compton processes; this leads to a natural analogy of Sw 1644+57 with a smaller-scale blazar. The phenomenologically novel Sw 1644+57 thus connects the study of TDFs and active galaxies, opening a new vista on disk-jet interactions in BHs and magnetic field generation and transport in accretion systems.

preprint2011arXiv

Are SMBHs shrouded by "super-Oort" clouds of comets and asteroids?

The last decade has seen a dramatic confirmation that an in situ star formation is possible inside the inner parsec of the Milky Way. Here we suggest that giant planets, solid terrestrial-like planets, comets and asteroids may also form in these environments, and that this may have observational implications for Active Galactic Nuclei (AGN). Like in debris discs around main sequence stars, collisions of large solid objects should initiate strong fragmentation cascades. The smallest particles in such a cascade - the microscopic dust - may provide a significant opacity. We put a number of observational and physical constraints on AGN obscuring torii resulting from such fragmentation cascades. We find that torii fed by fragmenting asteroids disappear at both low and high AGN luminosities. At high luminosities, $L \sim L_{\rm Edd}$, where $L_{\rm Edd}$ is the Eddington limit, the AGN radiation pressure blows out the microscopic dust too rapidly. At low luminosities, on the other hand, the AGN discs may avoid gravitational fragmentation into stars and solids. We also note that these fragmentation cascades may be responsible for astrophysically "large" dust particles of approximately micrometer sizes that were postulated by some authors to explain unusual absorption properties of the AGN torii.

preprint2011arXiv

Feeding SMBHs through supersonic turbulence and ballistic accretion

It has long been recognised that the main obstacle to accretion of gas onto supermassive black holes (SMBHs) is large specific angular momentum. It is feared that the gas settles in a large scale disc, and that accretion would then proceed too inefficiently to explain the masses of the observed SMBHs. Here we point out that, while the mean angular momentum in the bulge is very likely to be large, the deviations from the mean can also be significant. Indeed, cosmological simulations show that velocity and angular momentum fields of gas flows onto galaxies are very complex. Furthermore, inside bulges the gas velocity distribution can be further randomised by the velocity kicks due to feedback from star formation. We perform hydrodynamical simulations of gaseous rotating shells infalling onto an SMBH, attempting to quantify the importance of velocity dispersion in the gas at relatively large distances from the black hole. We implement this dispersion by means of a supersonic turbulent velocity spectrum. We find that, while in the purely rotating case the circularisation process leads to efficient mixing of gas with different angular momentum, resulting in a low accretion rate, the inclusion of turbulence increases this accretion rate by up to several orders of magnitude. We show that this can be understood based on the notion of "ballistic" accretion, whereby dense filaments, created by convergent turbulent flows, travel through the ambient gas largely unaffected by hydrodynamical drag. We derive a simple analytical formula that captures the numerical results to within a factor of a few. Rescaling our results to astrophysical bulges, we argue that this "ballistic" mode of accretion could provide the SMBHs with a sufficient supply of fuel without the need to channel the gas via large-scale discs or bars, and therefore that star formation in bulges can be a strong catalyst for SMBH accretion.

preprint2011arXiv

Fu Ori outbursts and the planet-disc mass exchange

It has been recently proposed that giant protoplanets migrating inward through the disc more rapidly than they contract could be tidally disrupted when they fill their Roche lobes ~ 0.1 AU away from their parent protostars. Here we consider the process of mass and angular momentum exchange between the tidally disrupted planet and the surrounding disc in detail. We find that the planet's adiabatic mass-radius relation and its ability to open a deep gap in the disc determine whether the disruption proceeds as a sudden runaway or a balanced quasi-static process. In the latter case the planet feeds the inner disc through its Lagrangian L1 point like a secondary star in a stellar binary system. As the planet loses mass it gains specific angular momentum and normally migrates in the outward direction until the gap closes. Numerical experiments show that planet disruption outbursts are preceded by long "quiescent" periods during which the disc inward of the planet is empty. The hole in the disc is created when the planet opens a deep gap, letting the inner disc to drain onto the star while keeping the outer one stalled behind the planet. We find that the mass-losing planet embedded in a realistic protoplanetary disc spawns an extremely rich set of variability patterns. In a subset of parameter space, there is a limit cycle behaviour caused by non-linear interaction between the planet mass loss and the disc hydrogen ionisation instability. We suggest that tidal disruptions of young massive planets near their stars may be responsible for the observed variability of young accreting protostars such as FU Ori, EXor and T Tauri stars in general.

preprint2011arXiv

Hot Super Earths: disrupted young jupiters?

Recent {\em Kepler} observations revealed an unexpected abundance of "hot" Earth-size to Neptune-size planets in the inner $0.02-0.2$ AU from their parent stars. We propose that these smaller planets are the remnants of massive giant planets that migrated inward quicker than they could contract. We show that such disruptions naturally occur in the framework of the Tidal Downsizing hypothesis for planet formation. We find that the characteristic planet-star separation at which such "hot disruptions" occur is $R \approx 0.03-0.2$ AU. This result is independent of the planet's embryo mass but is dependent on the accretion rate in the disc. At high accretion rates, $\dot M \simgt 10^{-6}\msun$ yr$^{-1}$, the embryo is unable to contract quickly enough and is disrupted. At late times, when the accretion rate drops to $\dot M \simlt 10^{-8} \msun$ yr$^{-1}$, the embryos migrate sufficiently slow to not be disrupted. These "late arrivals" may explain the well known population of hot jupiters. If type I migration regime is inefficient, then our model predicts a pile-up of planets at $R\sim 0.1$ AU as the migration rate suddenly switches from the type II to type I in that region.

preprint2011arXiv

Hot Super Earths: Tidally Downsized hot jupiters? The hot jupiter around the metal poor star explained

Recent Kepler observations revealed an unexpected abundance of "hot" Earth-size to Neptune-size planets in the inner ~0.02-0.2 AU from their parent stars. Another "oddity" of recent observations is the discovery of a hot jupiter planet around a very metal poor star of extragalactic origin. We point out that these observations have natural explanations in the framework of the Tidal Downsizing hypothesis, a new and an alternative model for planet formation. We show that disassociation of molecular hydrogen leads to the "second collapse" of the giant planet embryos. The collapse forms much denser fluid configurations, closely resembling isolated proto-giant planets on the Hayashi track before they enter the degenerate "main sequence" phase. We find that, in the early gas-rich phase of proto-planetary discs, these configurations migrate inward more rapidly than they contract. They are therefore tidally disrupted at characteristic distance of ~ 0.05 -0.1 AU from their parent stars. The disruption is likely to leave solid cores behind, which we identify with the hot Super-Earths observed by Kepler. In addition, we show that single hot jupiters (also found by the Kepler telescope) are not surprising in the context of TD hypothesis. Finally, formation of coreless giant planets by the second collapse of the giant embryos occurs for either very massive embryos (more than ~ 20 Jupiter masses) or very metal poor ones. We suggest that this process explains the extra-galactic hot jupiter.

preprint2011arXiv

Sgr A* flares: tidal disruption of asteroids and planets?

It is theoretically expected that a supermassive black hole (SMBH) in the centre of a typical nearby galaxy disrupts a Solar-type star every ~ 10^5 years, resulting in a bright flare lasting for months. Sgr A*, the resident SMBH of the Milky Way, produces (by comparison) tiny flares that last only hours but occur daily. Here we explore the possibility that these flares could be produced by disruption of smaller bodies - asteroids. We show that asteroids passing within an AU of Sgr A* could be split into smaller fragments which then vaporise by bodily friction with the tenuous quiescent gas accretion flow onto Sgr A*. The ensuing shocks and plasma instabilities may create a transient population of very hot electrons invoked in several currently popular models for Sgr A* flares, thus producing the required spectra. We estimate that asteroids larger than ~ 10 km in size are needed to power the observed flares, with the maximum possible luminosity of the order 10^39 erg s^-1. Assuming that the asteroid population per parent star in the central parsec of the Milky Way is not too dissimilar from that around stars in the Solar neighbourhood, we estimate the asteroid disruption rates, and the distribution of the expected luminosities, finding a reasonable agreement with the observations. We also note that planets may be tidally disrupted by Sgr A* as well, also very infrequently. We speculate that one such disruption may explain the putative increase in Sgr A* luminosity ~ 300 yr ago.

preprint2011arXiv

The Tidal Downsizing hypothesis for planet formation and the composition of Solar System comets

Comets are believed to be born in the outer Solar System where the temperature is assumed to have never exceeded T ~ 100 K. Surprisingly, observations and samples of cometary dust particles returned to Earth showed that they are in fact made of a mix of ices, as expected, but also of materials forged at high-temperatures (T ~ 1500 K). We propose a radically new view regarding the origin of the high-temperature processed materials in comets, based on the recent "Tidal Downsizing" (TD) hypothesis for planet formation. In the latter, the outer proto-planetary disc is gravitationally unstable and forms massive giant planet embryos (GEs). These hot (T ~ hundreds to 2000 K) and dense regions, immersed in the background cold and low density disc, are eventually disrupted. We propose that both planets and the high-T materials in comets are synthesised inside the GEs. Disruption of GEs separates planets and small solids as the latter are "frozen-in" into gas and are peeled off together with it. These small solids are then mixed with the ambient cold disc containing ices before being incorporated into comets. Several predictions of this picture may be testable with future observations of the Solar System and exoplanets.

preprint2010arXiv

A new view on planet formation

The standard picture of planet formation posits that giant gas planets are over-grown rocky planets massive enough to attract enormous gas atmospheres. It has been shown recently that the opposite point of view is physically plausible: the rocky terrestrial planets are former giant planet embryos dried of their gas "to the bone" by the influences of the parent star. Here we provide a brief overview of this "Tidal Downsizing" hypothesis in the context of the Solar System structure.

preprint2010arXiv

Feedback from Nuclear Star Clusters and SMBHs

The observed super-massive black hole (SMBH) mass -- galaxy velocity dispersion ($M_{\rm cmo} - σ$) correlation, and the similar correlation for nuclear star clusters, may be established when winds/outflows from the CMO ("central massive object") drive gas out of the potential wells of classical bulges. Timescales of growth for these objects may explain why smaller bulges appear to host preferentially NCs while larger ones contain SMBHs only. Despite much recent progress, feedback processes in bulge/galaxy formation are far from being understood. Our numerical simulations show that that understanding how the CMO feeds is as important a piece of the puzzle as understanding how its feedback affects its host galaxy.

preprint2010arXiv

Formation of planets by tidal downsizing of giant planet embryos

We hypothesise that planets are made by tidal downsizing of migrating giant planet embryos. The proposed scheme for planet formation consists of these steps: (i) a massive young protoplanetary disc fragments at R ~ several tens to hundreds of AU on gaseous clumps with masses of a few Jupiter masses; (ii) the clumps cool and contract, and simultaneously migrate closer in to the parent star; (iii) as earlier suggested by Boss (1998), dust sediments inside the gas clumps to form terrestrial mass solid cores; (iv) if the solid core becomes more massive than ~ 10 Earth masses, a massive gas atmosphere collapses onto the solid core; (v) when the gas clumps reach the inner few AU from the star, tidal shear and evaporation due to stellar irradiation peel off the outer metal-poor envelope of the clump. If tidal disruption occurs quickly, while the system is still in stage (iii), a terrestrial planet core is left. If it happens later, in stage (iv), a metal rich gas giant planet with a solid core emerges from the envelope.

preprint2010arXiv

Formation of terrestrial planet cores inside giant planet embryos

Giant planet embryos are believed to be spawned by gravitational instability in massive extended (R ~ 100 AU) protostellar discs. In a recent paper we have shown that dust can sediment inside the embryos, as argued earlier by Boss (1998) in a slightly different model. Here we study numerically the next stage of this process -- the formation of a solid core. If conditions are conducive to solid core formation, the centre of the gas cloud goes through the following sequence of phases: (i) becomes grain (and metal) rich; (ii) forms a terrestrial mass solid core via a rapid collapse driven by self-gravity of the grains; (iii) starts to accrete a gaseous atmosphere when the solid core reaches mass of a few to 10 Earth masses. This sequence of events may build either terrestrial planet cores or metal rich giant planets inside the larger gas reservoir of the giant planet embryo. In a companion Letter we argue that tidal and irradiation effects from the parent star should disrupt the outer metal poor layers of the embryo, releasing nearly "ready to use" planets. We propose this as an alternative way to build planets.

preprint2010arXiv

Grain sedimentation inside giant planet embryos

In the context of massive fragmenting protoplanetary discs, Boss (1998) suggested that grains can grow and sediment inside giant planet embryos formed at R ~ 5 AU away from the star. Several authors since then criticised the suggestion. Convection may prevent grain sedimentation, and the embryos cannot even form so close to the parent star as cooling is too inefficient at these distances. Here we reconsider the grain sedimentation process suggested by Boss (1998) but inside an embryo formed, as expected in the light of the cooling constraints, at R ~ 100 AU. Such embryos are much less dense and are also cooler. We make analytical estimates of the process and also perform simple spherically symmetric radiation hydrodynamics simulations to test these ideas. We find that convection in our models does not become important before a somewhat massive (~ an Earth mass, this is clarified in a followup paper) solid core is built. Turbulent mixing slows down dust sedimentation but is overwhelmed by grain sedimentation when the latter grow to a centimetres size. The minimum time required for dust sedimentation to occur is a few thousand years, and is a strong function of the embryo's mass, dust content and opacity. An approximate analytical criterion is given to delineate conditions in which a giant embryo contracts and heats up faster than dust can sediment. As Boss et al (2002), we argue that core formation through grain sedimentation inside the giant planet embryos may yield an unexplored route to form giant gas and giant ice planets. The present model also stands at the basis of paper III, where we study the possibility of forming terrestrial planet cores by tidal disruption and photoevaporation of the planetary envelope.

preprint2010arXiv

Rotation of the Solar System planets and the origin of the Moon in the context of the tidal downsizing hypothesis

It has been proposed recently that the first step in the formation of both rocky and gas giant planets is dust sedimentation into a solid core inside a gas clump (giant planet embryo). The clumps are then assumed to migrate closer to the star where their metal poor envelopes are sheared away by the tidal forces or by an irradiation-driven mass loss. We consider the implications of this hypothesis for natal rotation rates of both terrestrial and gas giant planets. It is found that both types of planets may rotate near their break up angular frequencies at birth. The direction of the spin should coincide with that of the parent disc and the star, except in cases of embryos that had close interactions or mergers with other embryos in the past. Furthermore, the large repository of specific angular momentum at birth also allows formation of close binary rocky planets inside the same embryos. We compare these predictions with rotation rates of planets in the Solar System and also question whether the Earth-Moon pair could have been formed within the same giant planet embryo.

preprint2010arXiv

The Accretion Disc Particle Method for Simulations of Black Hole Feeding and Feedback

Black holes grow by accreting matter from their surroundings. However, angular momentum provides an efficient natural barrier to accretion and so only the lowest angular momentum material will be available to feed the black holes. The standard sub-grid model for black hole accretion in galaxy formation simulations - based on the Bondi-Hoyle method - does not account for the angular momentum of accreting material, and so it is unclear how representative the black hole accretion rate estimated in this way is likely to be. In this paper we introduce a new sub-grid model for black hole accretion that naturally accounts for the angular momentum of accreting material. Both the black hole and its accretion disc are modelled as a composite accretion disc particle. Gas particles are captured by the accretion disc particle if and only if their orbits bring them within its accretion radius R_acc, at which point their mass is added to the accretion disc and feeds the black hole on a viscous timescale t_visc. The resulting black hole accretion rate (dM/dt)_BH powers the accretion luminosity L_acc ~ (dM/dt)_BH, which drives black hole feedback. Using a series of controlled numerical experiments, we demonstrate that our new accretion disc particle method is more physically self-consistent than the Bondi-Hoyle method. We also discuss the physical implications of the accretion disc particle method for systems with a high degree of rotational support, and we argue that the M_BH-sigma relation in these systems should be offset from the relation for classical bulges and ellipticals, as appears to be observed.

preprint2009arXiv

Competitive feedback in galaxy formation

It is now well established that many galaxies have nuclear star clusters (NCs) whose total masses correlate with the velocity dispersion (sigma) of the galaxy spheroid in a very similar way to the well--known supermassive black hole (SMBH) M - sigma relation. Previous theoretical work suggested that both correlations can be explained by a momentum feedback argument. Observations further show that most known NCs have masses < 10^8 Msun, while SMBHs frequently have masses > 10^8 Msun, which remained unexplained in previous work. We suggest here that this changeover reflects a competition between the SMBH and nuclear clusters in the feedback they produce. When one of the massive objects reaches its limiting M-sigma value, it drives the gas away and hence cuts off its own mass and also the mass of the ``competitor''. The latter is then underweight with respect to the expected M-sigma mass (abridged).

preprint2009arXiv

Dynamic Monte Carlo radiation transfer in SPH. Radiation pressure force implementation

We present a new framework for radiation hydrodynamics simulations. Gas dynamics is modelled by the Smoothed Particle Hydrodynamics (SPH) method, whereas radiation transfer is simulated via a time-dependent Monte-Carlo approach that traces photon packets. As a first step in the development of the method, in this paper we consider the momentum transfer between radiation field and gas, which is important for systems where radiation pressure is high. There is no fundamental limitations on the number of radiation sources, geometry or the optical depth of the problems that can be studied with the method. However, as expected for any Monte-Carlo transfer scheme, stochastic noise presents a serious limitation. We present a number of tests that show that the errors of the method can be estimated accurately by considering Poisson noise fluctuations in the number of photon packets that SPH particles interact with per dynamical time. It is found that for a reasonable accuracy the momentum carried by photon packets must be much smaller than a typical momentum of SPH particles. We discuss numerical limitations of the code, and future steps that can be taken to improve performance and applicability of the method.

preprint2009arXiv

Kelvin-Helmholtz instabilities with Godunov SPH

Numerical simulations for the non-linear development of Kelvin-Helmholtz instability in two different density layers have been performed with the particle-based method (Godunov SPH) developed by Inutsuka (2002). The Godunov SPH can describe the Kelvin-Helmholtz instability even with a high density contrast, while the standard SPH shows the absence of the instability across a density gradient (Agertz et al. 2007). The interaction of a dense blob with a hot ambient medium has been performed also. The Godunov SPH describes the formation and evolution of the fingers due to the combinations of Rayleigh-Taylor, Richtmyer-Meshkov, and Kelvin-Helmholtz instabilities. The blob test result coincides well with the results of the grid-based codes. An inaccurate handling of a density gradient in the standard SPH has been pointed out as the direct reason of the absence of the instabilities. An unphysical force happens at the density gradient even in a pressure equilibrium, and repulses particles from the initial density discontinuity. Therefore, the initial perturbation damps, and a gap forms at the discontinuity. The unphysical force has been studied in terms of the consistency of a numerical scheme. Contrary to the standard SPH, the momentum equation of the Godunov SPH doesnt use the particle approximation, and has been derived from the kernel convolution or a new Lagrangian function. The new Lagrangian function used in the Godunov SPH is more analogous to the real Lagrangian function for continuum. The momentum equation of the Godunov SPH has much better linear consistency, so the unphysical force is greatly reduced compared to the standard SPH in a high density contrast.

preprint2009arXiv

Simulations of momentum feedback by black hole winds

The observed super-massive black hole (SMBH) mass -- galaxy velocity dispersion ($M_{\rm bh} - σ$) correlation may be established when winds/outflows from the SMBH drive gas out of the potential wells of classical bulges. Here we present numerical simulations of this process in a static isothermal potential. Simple spherically symmetric models of SMBH feedback at the Eddington luminosity can successfully explain the $M_{\rm bh} - σ$ and nuclear cluster mass $M_{\rm NC}-σ$ correlations, as well as why larger bulges host SMBHs while smaller ones host nuclear star clusters. However these models do not specify how SMBHs feed on infalling gas whilst simultaneously producing feedback that drives gas out of the galaxy. More complex models with rotation and/or anisotropic feedback allow SMBHs to feed via a disc or regions not exposed to SMBH winds, but in these more realistic cases it is not clear why a robust $M_{\rm bh} - σ$ relation should be established. In fact, some of the model predictions contradict observations. For example, an isotropic SMBH wind impacting on a disc (rather than a shell) of aspect ratio $H/R \ll 1$ requires the SMBH mass to be larger by a factor $\sim R/H$, which is opposite to what is observed. We conclude that understanding how a SMBH feeds is as important a piece of the puzzle as understanding how its feedback affects its host galaxy. Finally, we note that in aspherical cases the SMBH outflows induce differential motions in the bulge. This may pump turbulence that is known to hinder star formation in star forming regions. SMBH feedback thus may not only drive gas out of the bulge but also reduce the fraction of gas turned into stars.

preprint2008arXiv

Simulations of the formation of a gaseous disc and young stars near Sgr A* via cloud-cloud collisions

Young massive stars in the central parsec of our Galaxy are best explained by star formation within at least one, and possibly two, massive self-gravitating gaseous discs. With help of numerical simulations, we here consider whether the observed population of young stars could have originated from a large angle collision of two massive gaseous clouds at R ~ 1 pc from Sgr A*. In all the simulations performed, the post-collision gas flow forms an inner nearly circular gaseous disc and one or two eccentric outer filaments, consistent with the observations. Furthermore, the radial stellar mass distribution is always very steep, sigma proportional to R^-2, again consistent with the data. The 3D velocity structure of the stellar distribution is however sensitive to initial conditions (e.g., the impact parameter of the clouds) and gas cooling details. In all the cases the amount of gas accreted by our inner boundary condition is large, enough to allow Sgr A* to radiate near its Eddington limit during ~ 10^5 years. This suggests that a refined model would have physically larger clouds (or a cloud and a disc such as the CND) colliding at a distance of a few pc rather than 1 pc as in our simulations.

preprint2006arXiv

The M-sigma Relation for Nucleated Galaxies

Momentum feedback from super-Eddington accretion offers a simple explanation for the observed M-sigma and M-M_{spheroid} relations between supermassive black holes (SMBH) and the spheroids of their host galaxies. Recently Ferrarese et al. and Wehner & Harris observed analogous relations between the masses of central star clusters and their hosts. We show that stellar winds and supernovae from such nuclear clusters (NC) give similar feedback explanations for this case also, and we discuss the connection to the Faber-Jackson relation for the spheroids themselves.

preprint1999arXiv

On the X-ray heated skin of Accretion Disks

We present a simple analytical formula for the Thomson depth of the X-ray heated skin of accretion disks valid at any radius and for a broad range of spectral indices of the incident X-rays, accretion rates and black hole masses. We expect that this formula may find useful applications in studies of geometry of the inner part of accretion flows around compact objects, and in several other astrophysically important problems, such as the recently observed X-ray ``Baldwin'' effect (i.e., monotonic decrease of Fe line's equivalent width with the X-ray luminosity of AGN), the problem of missing Lyman edge in AGN, and line and continuum variability studies in accretion disks around compact objects. We compute the reflected X-ray spectra for several representative cases and show that for hard X-ray spectra and large ionizing fluxes the skin represents a perfect mirror that does not produce any Fe lines or absorption features. At the same time, for soft X-ray spectra or small ionizing fluxes, the skin produces very strong ionized absorption edge and highly ionized Fe lines that should be observable in the reflected spectra.

preprint1999arXiv

Time-Dependent Disk Models for the Microquasar GRS1915+105

We investigate different accretion disk models and viscosity prescriptions in order to provide a basic explanation for the exotic temporal behavior in GRS 1915+105. Based on the fact that the overall cycle times are very much longer than the rise/fall time scales in GRS 1915, we rule out the geometry of ADAF or a hot quasi-spherical region plus a cold outer disk for this source. We thus concentrate on geometrically thin Shakura-Sunyaev type disks (Shakura & Sunyaev 1973; hereafter SS73). We have devised a modified viscosity law that has a quasi-stable upper branch. Via numerical simulations, we show that the model does account for several gross observational features of GRS 1915+105. On the other hand, the rise/fall time scales are not short enough, and no rapid oscillations on time scales $\simlt$ 10 s emerge naturally from the model. We then consider and numerically test a more elaborate model that includes the cold disk, a corona, and plasma ejections from the inner disk region and show that this model allows us to reproduce several additional observed features of GRS 1915+105. We conclude that the most likely structure of the accretion flow in this source is that of a cold disk with a modified viscosity prescription, plus a corona that accounts for much of the X-ray emission, and unsteady plasma ejections that occur when the luminosity of the source is high.

preprint1997arXiv

Self-Consistent Fokker-Planck Treatment Of Particle Distributions in Astrophysical Plasmas

High-energy, multi-component plasmas in which pair creation and annihilation, lepton-lepton scattering, lepton-proton scattering, and Comptonization all contribute to establishing the particle and photon distributions, are present in a broad range of compact astrophysical objects. Earlier work has included much of the microphysics needed to account for electron-photon and electron-proton interactions, but little has been done to handle the redistribution of the particles as a result of their Coulomb interaction with themselves in an arbitrary case. Our goal here is to use a Fokker-Planck approach in order to develop a fully self-consistent theory for the interaction of arbitrarily distributed particles and radiation to arrive at an accurate representation of the high-energy plasma in these sources. We conduct several tests representative of two dominant segments of parameter space and discuss physical implications of the non-Maxwellian distribution function. Approximate analytical forms for the electron distribution in the case of strongly non-Maxwellian plasmas are given.

Sergei Nayakshin

What is connected

Connect this record

See the researcher in context

Building this map preview

58 published item(s)

Disc fragmentation. II. Ejection of low mass Free Floating Planets from growing binary systems

ALMA constraints on assembly of Core Accretion planets

Gap opening by planets in discs with magnetised winds

The growth and migration of massive planets under the influence of external photoevaporation

The paradox of youth for ALMA planet candidates

TW Hya: an old protoplanetary disc revived by its planet

A desert of gas giant planets beyond tens of au

Planets, debris and their host metallicity correlations

A simple way to improve AGN feedback prescription in SPH simulations

The resolution bias: low resolution feedback simulations are better at destroying galaxies

The role of feedback in accretion on Low Luminosity AGN: Sgr A* case study

Tidal Downsizing model. II. Planet-metallicity correlations

Tidal Downsizing Model. III. Planets from sub-Earths to Brown Dwarfs: structure and metallicity preferences

Tidal Downsizing Model. IV. Destructive feedback in planets

Black hole feedback in a multiphase interstellar medium

Core-assisted gas capture instability: a new mode of giant planet formation by gravitationally unstable discs

Energy- and momentum-conserving AGN feedback outflows

Metal loading of giant gas planets

Positive metallicity correlation for coreless giant planets

Tidal Downsizing model. I. Numerical methods: saving giant planets from tidal disruptions

A link between feedback outflows and satellite galaxy suppression

AGN outflows trigger starbursts in gas-rich galaxies

Hydrogen-loosing planets in transition discs around young protostars

Inverse Compton X-ray signature of AGN feedback

Outflows of stars due to quasar feedback

Radiative feedback from protoplanets in self-gravitating protoplanetary discs

Two-phase model for Black Hole feeding and feedback

An alternative origin for debris rings of planetesimals

Fermi Bubbles in the Milky Way: the closest AGN feedback laboratory courtesy of Sgr A*?

Modelling supermassive black hole growth: towards an improved sub-grid prescription

Quasar feedback: accelerated star formation and chaotic accretion

The observed $M - σ$ relations imply that SMBHs grow by cold chaotic accretion

What's in a Fermi Bubble: a quasar episode in the Galactic centre

A numerical simulation of a "super-Earth" core delivery from ~ 100 AU to ~ 8 AU

A relativistic jetted outburst from a massive black hole fed by a tidally disrupted star

Are SMBHs shrouded by "super-Oort" clouds of comets and asteroids?

Feeding SMBHs through supersonic turbulence and ballistic accretion

Fu Ori outbursts and the planet-disc mass exchange

Hot Super Earths: disrupted young jupiters?

Hot Super Earths: Tidally Downsized hot jupiters? The hot jupiter around the metal poor star explained

Sgr A* flares: tidal disruption of asteroids and planets?

The Tidal Downsizing hypothesis for planet formation and the composition of Solar System comets

A new view on planet formation

Feedback from Nuclear Star Clusters and SMBHs

Formation of planets by tidal downsizing of giant planet embryos

Formation of terrestrial planet cores inside giant planet embryos

Grain sedimentation inside giant planet embryos

Rotation of the Solar System planets and the origin of the Moon in the context of the tidal downsizing hypothesis

The Accretion Disc Particle Method for Simulations of Black Hole Feeding and Feedback

Competitive feedback in galaxy formation

Dynamic Monte Carlo radiation transfer in SPH. Radiation pressure force implementation

Kelvin-Helmholtz instabilities with Godunov SPH

Simulations of momentum feedback by black hole winds

Simulations of the formation of a gaseous disc and young stars near Sgr A* via cloud-cloud collisions

The M-sigma Relation for Nucleated Galaxies

On the X-ray heated skin of Accretion Disks

Time-Dependent Disk Models for the Microquasar GRS1915+105

Self-Consistent Fokker-Planck Treatment Of Particle Distributions in Astrophysical Plasmas