Researcher profile

Pavel Kroupa

Pavel Kroupa contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
27works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

27 published item(s)

preprint2026arXiv

Intermediate-mass black hole incubators. Gas accretion onto stellar black hole clusters in galactic central molecular zones

The stellar dynamical evolution of massive star clusters formed during starburst periods leads to the segregation of $\gtrsim10^4 M_\odot$ stellar-mass black hole sub-clusters in their centres. In gas-rich environments, such as galactic central molecular zones, these black hole clusters are likely to accrete large amounts of the gas from their surroundings, which in turn affects their internal dynamics. In this Letter we estimated the corresponding accretion rate onto the black hole cluster and its radiative feedback. We assessed whether such an accretion flow can lead to the collapse of the black hole cluster into an intermediate-mass black hole. The estimates were obtained analytically, considering the astrophysical conditions and star formation history reported for the central molecular zone of our Galaxy. We find that a stellar black hole cluster with mass $\gtrsim10^4 M_\odot$ located in the twisted ring of molecular clouds with radius $\approx100$ pc that is observed in the central molecular zone of our Galaxy can accrete about the same mass in gas on a timescale of a few million years. We suggest that this is sufficient for its subsequent collapse into an intermediate-mass black hole. Based on an estimate of the dynamical friction inspiral time, we further argue that the locations of the intermediate-mass black hole candidates recently observed in the central molecular zone are compatible with their formation therein during the last starburst period reported to have occurred $\approx1$ Gyr ago.

preprint2025arXiv

Tidal Tails of Nearby Open Clusters. II. A Review of Simulated Properties and the Reliability of Observational Catalogues

Context: Recent studies using Gaia data have reported tidal tail detections for tens to hundreds of open clusters. However, a comprehensive assessment of the reliability & completeness of these detections is lacking. Aims: This work aims to summarise the expected properties of tidal tails based on N-body simulations, review the reliability of tidal tail detections in the literature, and grade them according to a set of diagnostic tests. Methods: We used a grid of 68-20000 Msun simulated clusters & analysed the formation & evolution of the tails. We compiled 122 catalogues (58 unique clusters) from literature, within 500 pc of the Sun. We employed tests based on photometric, morphological & dynamical signatures, & comparisons with simulations, to grade their tails. Results: Based on the simulations, we analysed the complex morphology of the tails & their properties (e.g., size, span, stellar types, number density, & mass function) at various cluster masses & ages. During the first 200 Myr of evolution, the tails typically form a characteristic S shape, with an amplitude that scales with cluster mass. The tail span increases at a rate of 4 times the initial velocity dispersion, & the near-tail is predominantly populated by recent escapees. Conclusions: In the 122 published tail catalogues, we found that 15 gold & 55 silver catalogues passed the majority of the tests. The remaining 51 were graded as bronze; care should be taken before using these for further analysis. The age, metallicity, binary fraction, & mass function of stars in the tails were generally consistent with those of their parent clusters. The gold/silver catalogues (69 catalogues of 40 clusters) represent reliable samples for detailed analyses of tails. Future data will be essential for further validation & for leveraging tidal tails as tracers of cluster dissolution & the Galactic potential.

preprint2022arXiv

3D hydrodynamic simulations for the formation of the Local Group satellite planes

The existence of mutually correlated thin and rotating planes of satellite galaxies around both the Milky Way (MW) and Andromeda (M31) calls for an explanation. Previous work in Milgromian dynamics (MOND) indicated that a past MW-M31 encounter might have led to the formation of these satellite planes. We perform the first-ever hydrodynamical MOND simulation of the Local Group using Phantom of RAMSES. We show that an MW-M31 encounter at $z \approx 1$, with a perigalactic distance of about 80 kpc, can yield two disc galaxies at $z=0$ oriented similarly to the observed galactic discs and separated similarly to the observed M31 distance. Importantly, the tidal debris are distributed in phase space similarly to the observed MW and M31 satellite planes, with the correct preferred orbital pole for both. The MW-M31 orbital geometry is consistent with the presently observed M31 proper motion despite this not being considered as a constraint when exploring the parameter space. The mass of the tidal debris around the MW and M31 at $z=0$ compare well with the mass observed in their satellite systems. The remnant discs of the two galaxies have realistic radial scale lengths and velocity dispersions, and the simulation naturally produces a much hotter stellar disc in M31 than in the MW. However, reconciling this scenario with the ages of stellar populations in satellite galaxies would require that a higher fraction of stars previously formed in the outskirts of the progenitors ended up within the tidal debris, or that the MW-M31 interaction occurred at $z>1$.

preprint2022arXiv

Do the majority of stars form as gravitationally unbound?

Some of the youngest stars (age $\lesssim 10$ Myr) are clustered, while many others are observed scattered throughout star forming regions or in complete isolation. It has been intensively debated whether the scattered or isolated stars originate in star clusters, or if they form truly isolated, which could help constrain the possibilities how massive stars are formed. We adopt the assumption that all stars form in gravitationally bound star clusters embedded in molecular cloud cores ($Γ$-$1\;$ model), which expel their natal gas, and compare the fraction of stars found in clusters with observational data. The star clusters are modelled by the code nbody6, which includes stellar and circumbinary evolution, gas expulsion, and the external gravitational field of their host galaxy. We find that small changes in the assumptions in the current theoretical model estimating the fraction, $Γ$, of stars forming in embedded clusters have a large influence on the results, and we present a counterexample as an illustration. This calls into question theoretical arguments about $Γ$ in embedded clusters, and it suggests that there is no firm theoretical ground for low $Γ$ in galaxies with lower star formation rates (SFRs). Instead, the assumption that all stars form in embedded clusters is in agreement with observational data for the youngest stars (age $\lesssim 10$ Myr). In the $Γ$-$1\;$ scenario, the observed fraction of the youngest stars in clusters increases with the SFR only weakly; the increase is caused by the presence of more massive clusters in galaxies with higher SFRs, which release fewer stars to the field in proportion to their mass. The $Γ$-$1\;$ model yields a higher fraction of stars in clusters for older stars (age between $10$ and $300$ Myr) than what is observed. This discrepancy can be caused by interactions with molecular clouds.

preprint2022arXiv

Far-ultraviolet investigation into the galactic globular cluster M30 (NGC 7099): I. Photometry and radial distributions

We present a far-ultraviolet (FUV) study of the globular cluster M30 (NGC 7099). The images were obtained using the Advanced Camera for Surveys (ACS/SBC, F150LP, FUV) and the Wide Field Planetary Camera 2 (WFPC2, F300W, UV) which were both on board the Hubble Space Telescope (HST). The FUV-UV colour-magnitude diagram (CMD) shows a main sequence (MS) turnoff at FUV $\approx$ 22 mag and FUV-UV $\approx$ 3 mag. The MS extends 4 mag below the turnoff, and a prominent horizontal branch (HB) and blue straggler (BS) sequence can be seen. A total of 1218 MS stars, 185 red giant branch stars, 47 BS stars and 41 HB stars are identified, along with 78 sources blueward of the MS which consist of white dwarfs (WDs) and objects in the gap between the WDs and the MS that include potential cataclysmic variable (CV) candidates. The radial distribution of the BS population is concentrated towards the cluster centre, indicating that mass segregation has occurred. The blue and red sub-populations of the double BS sequence appear mixed in the ultraviolet CMD, and no significant central concentration of CV candidates is seen in this cluster.

preprint2022arXiv

Far-ultraviolet investigation into the galactic globular cluster M30 (NGC 7099): II. Potential X-ray counterparts and variable sources

We present a far-ultraviolet (FUV) study of the globular cluster M30 (NGC 7099). The images were obtained using the Advanced Camera for Surveys (ACS/SBC, F150LP, FUV) and the Wide Field Planetary Camera 2 (WFPC2, F300W, UV) on board the Hubble Space Telescope (HST). We compare the catalogue of FUV objects to ten known X-ray sources and find six confident matches of two cataclysmic variables (CVs), one RS CVn, one red giant with strong FUV emission and two sources only detected in the FUV. We also searched for variable sources in our dataset and found a total of seven blue stragglers (BSs), four horizontal branch (HB) stars, five red giant branch stars, 28 main sequence stars and four gap objects that demonstrated variability. One BS star is a known W-UMa contact binary, one of the gap objects is a known CV identified in this work to be a dwarf nova, and the three other gap sources are weak variables. The periods and positions of two of the variable HB stars match them to two previously known RR Lyrae variables of types RRab and RRc.

preprint2022arXiv

Origin of the spectacular tidal shells of galaxy NGC474

The lenticular galaxy NGC474 hosts a rich system of tidal shells and streams, some of which are exceptionally bright. Two teams recently presented spectroscopic observations of the brightest shells. These were the first shell spectra ever observed in integrated starlight. The authors studied the stellar populations of the shell, of the center of the galaxy and of its globular clusters. The precise formation scenario for the tidal features of this prominent galaxy however still remained unclear. Here, we add further clues on their formation from the radii of the shells, and we present a scenario for the formation of the tidal features that seems to be unique and explaining all available data. Shell radii are analyzed with the shell identification method, and we run self-consistent simulations of the formation of the tidal features. We consider Newtonian as well as MOND gravity. Observations suggest that the tidal features originate from the accretion of a spiral galaxy. The shell identification method yields that the merging galaxies collided first 1.3Gyr ago and then again 0.9Gyr ago, thereby forming the shells in two generations. This would also explain the young ages of stellar populations in the center of the galaxy and the young age of the globular clusters. The analytic models of shell propagation, that underlie the shell identification method, are verified by a simulation. The simulations reproduce well the observed morphology of the tidal features. The accreted spiral likely reached NGC474 nearly radially, in the plane of the sky, from the south, its rotation axis pointing toward us. It should have had a stellar mass of around 1/6 of NGC474, i.e. $10^{9.8}\,M_\odot$. It seems that all tidal features in the galaxy originate from one merger.

preprint2022arXiv

Overestimated inclinations of Milgromian disc galaxies: the case of the ultradiffuse galaxy AGC 114905

We present hydrodynamical star-forming simulations in the Milgromian dynamics (MOND) framework of a gas-rich disc galaxy with properties similar to AGC 114905, which has recently been argued to have a rotation curve (RC) that is inconsistent with the MOND prediction. Our first model considers the galaxy in isolation, while our second model includes an external field of $0.05 \, a_{_0}$, the estimated gravitational field from large-scale structure. We show that isophotes in the face-on view can differ from circular at the 50% level. This could mislead observers into overestimating the inclination $i$ between disc and sky planes. Because RCs require a correction factor of $1/\sin i$, the actual RC could be much higher than reported by observers. This plausibly reconciles AGC 114905 with MOND expectations.

preprint2022arXiv

The formation of early-type galaxies through monolithic collapse of gas clouds in Milgromian gravity

Studies of stellar populations in early-type galaxies (ETGs) show that the more massive galaxies form earlier and have a shorter star formation history (SFH). In this study, we investigate the initial conditions of ETG formation. The study begins with the collapse of non-rotating post-Big-Bang gas clouds in Milgromian (MOND) gravitation. These produce ETGs with star-forming timescales (SFT) comparable to those observed in the real Universe. Comparing these collapse models with observations, we set constraints on the initial size and density of the post-Big-Bang gas clouds in order to form ETGs. The effective-radius-mass relation of the model galaxies falls short of the observed relation. Possible mechanisms for later radius expansion are discussed. Using hydrodynamic MOND simulations this work thus for the first time shows that the SFTs observed for ETGs may be a natural occurrence in the MOND paradigm. We show that different feedback algorithms change the evolution of the galaxies only to a very minor degree in MOND. The first stars have, however, formed more rapidly in the real Universe than possible just from the here studied gravitational collapse mechanism. Dark-matter-based cosmological structure formation simulations disagree with the observed SFTs at more than 5 sigma confidence.

preprint2022arXiv

The global stability of M33 in MOND

The dynamical stability of disk galaxies is sensitive to whether their anomalous rotation curves are caused by dark matter halos or Milgromian Dynamics (MOND). We investigate this by setting up a MOND model of M33. We first simulate it in isolation for 6 Gyr, starting from an initial good match to the rotation curve (RC). Too large a bar and bulge form when the gas is too hot, but this is avoided by reducing the gas temperature. A strong bar still forms in 1 Gyr, but rapidly weakens and becomes consistent with the observed weak bar. Previous work showed this to be challenging in Newtonian models with a live dark matter halo, which developed strong bars. The bar pattern speed implies a realistic corotation radius of 3 kpc. However, the RC still rises too steeply, and the central line of sight velocity dispersion (LOSVD) is too high. We then add a constant external acceleration field of $8.4 \times 10^{-12}$ m/s$^2$ at $30^\circ$ to the disk as a first order estimate for the gravity exerted by M31. This suppresses buildup of material at the centre, causing the RC to rise more slowly and reducing the central LOSVD. Overall, this simulation bears good resemblance to several global properties of M33, and highlights the importance of including even a weak external field on the stability and evolution of disk galaxies. Further simulations with a time-varying external field, modeling the full orbit of M33, will be needed to confirm its resemblance to observations.

preprint2022arXiv

The High Fraction of Thin Disk Galaxies Continues to Challenge ΛCDM Cosmology

Any viable cosmological framework has to match the observed proportion of early- and late-type galaxies. In this contribution, we focus on the distribution of galaxy morphological types in the standard model of cosmology (Lambda cold dark matter, $Λ$CDM). Using the latest state-of-the-art cosmological $Λ$CDM simulations known as Illustris, IllustrisTNG, and EAGLE, we calculate the intrinsic and sky-projected aspect ratio distribution of the stars in subhalos with stellar mass $M_* > 10^{10}\,M_\odot$ at redshift $z=0$. There is a significant deficit of intrinsically thin disk galaxies, which however comprise most of the locally observed galaxy population. Consequently, the sky-projected aspect ratio distribution produced by these $Λ$CDM simulations disagrees with the Galaxy And Mass Assembly (GAMA) survey and Sloan Digital Sky Survey at $\geq 12.52σ$ (TNG50-1) and $\geq 14.82σ$ (EAGLE50) confidence. The deficit of intrinsically thin galaxies could be due to a much less hierarchical merger-driven build-up of observed galaxies than is given by the $Λ$CDM framework. It might also arise from the implemented sub-grid models, or from the limited resolution of the above-mentioned hydrodynamical simulations. We estimate that an $8^5$ times better mass resolution realization than TNG50-1 would reduce the tension with GAMA to the $5.58σ$ level. Finally, we show that galaxies with fewer major mergers have a somewhat thinner aspect ratio distribution. Given also the high expected frequency of minor mergers in $Λ$CDM, the problem may be due to minor mergers. In this case, the angular momentum problem could be alleviated in Milgromian dynamics (MOND) because of a reduced merger frequency arising from the absence of dynamical friction between extended dark matter halos.

preprint2021arXiv

Estimating the ages of open star clusters from properties of their extended tidal tails

The most accurate current methods for determining the ages of open star clusters, stellar associations and stellar streams are based on isochrone fitting or the lithium depletion boundary. We propose another method for dating these objects based on the morphology of their extended tidal tails, which have been recently discovered around several open star clusters. Assuming that the early-appearing tidal tails, the so called tidal tails I, originate from the stars released from the cluster during early gas expulsion, or that they form in the same star forming region as the cluster (i.e. being coeval with the cluster), we derive the analytical formula for the tilt angle $β$ between the long axis of the tidal tail and the orbital direction for clusters or streams on circular trajectories. Since at a given Galactocentric radius, $β$ is only a function of age $t$ regardless of the initial properties of the cluster, we estimate the cluster age by inverting the analytical formula $β= β(t)$. We illustrate the method on a sample of $12$ objects, which we compiled from the literature, and we find a reasonable agreement with previous dating methods in $\approx 70$% of the cases. This can probably be improved by taking into account the eccentricity of the orbits and by revisiting the dating methods based on stellar evolution. The proposed morphological method is suitable for relatively young clusters (age $\lesssim 300$ Myr), where it provides a relative age error of the order of $10$ to $20$% for an error in the observed tilt angle of $5 ^\circ$.

preprint2021arXiv

HI mapping of the Leo Triplet: Morphologies and kinematics of tails and bridges

A fully-sampled and hitherto highest resolution and sensitivity observation of neutral hydrogen (HI) in the Leo Triplet (NGC 3628, M 65/NGC 3623, and M 66/NGC 3627) reveals six HI structures beyond the three galaxies. We present detailed results of the morphologies and kinematics of these structures, which can be used for future simulations. In particular, we detect a two-arm structure in the plume of NGC 3628 for the first time, which can be explained by a tidal interaction model. The optical counterpart of the plume is mainly associated with the southern arm. The connecting part (base) of the plume (directed eastwards) with NGC 3628 is located at the blueshifted (western) side of NGC 3628. Two bases appear to be associated with the two arms of the plume. A clump with reversed velocity gradient (relative to the velocity gradient of M 66) and a newly detected tail, i.e. M 66SE, is found in the southeast of M 66. We suspect that M 66SE represents gas from NGC 3628 which was captured by M 66 in the recent interaction between the two galaxies. Meanwhile gas is falling toward M 66, resulting in features already previously observed in the southeastern part of M 66, e.g. large line widths and double peaks. An upside-down `Y'-shaped HI gas component (M 65S) is detected in the south of M 65 which suggests that M 65 may also have been involved in the interaction. We strongly encourage modern hydrodynamical simulations of this interacting group of galaxies to reveal the origin of the gaseous debris surrounding all three galaxies.

preprint2021arXiv

On the dynamical evolution of Cepheids in star clusters

We investigate the occurrence of classical (type-I) Cepheid variable stars (henceforth: Cepheids) in dynamically evolving star clusters from birth to an age of approximately 300 Myr. The clusters are modelled by the Aarseth code nbody6, and they feature a realistic stellar initial mass function and initial binary star population, single star and binary star evolution, expulsion of the primordial gas, and the tidal field of the galaxy. Our simulations provide the first detailed dynamical picture of how frequently Cepheids remain gravitationally bound to their birth clusters versus how frequently they occur in the field. They allow us to quantify the relevance of various cluster escape mechanisms and how they depend on stellar mass. Overall, the simulations agree with the empirical picture that a small fraction ($\approx 10\%$) of Cepheids reside in clusters, that cluster halo membership is relatively common, and that the majority of Cepheid hosting clusters have only a single Cepheid member. Additionally, the simulations predict that a) Cepheid progenitors are much more likely to escape from low-mass than higher-mass clusters; b) higher-mass (long-period) Cepheids are $\approx 30\%$ more likely to be found in clusters than low-mass (short-period) Cepheids; c) the clustered Cepheid fraction increases with galactocentric radius since cluster dispersal is less efficient at greater radii; d) a lower metallicity reduces the overall clustered Cepheid fraction; e) high-mass clusters are much more likely to have more than one Cepheid member at any given time, in particular at a lower metallicity. We interpret the results as outcomes of various aspects of star cluster dynamics. The comparison of predicted and observed clustered Cepheid fractions, $f_{\rm CC}$, highlights the need for additional cluster disruption mechanisms, most likely encounters with giant molecular clouds.

preprint2020arXiv

A correlation between the number of satellites and the baryonic bulge-to-total ratio extending beyond the Local Group

Recent observations of the fields surrounding a few Milky-Way-like galaxies in the local Universe have become deep enough to enable investigations of the predictions of the standard LCDM cosmological model down to small scales outside the Local Group. Motivated by an observed correlation between the number of dwarf satellites (N_sat) and the bulge-to-total baryonic mass ratios (B/T) of the three main galaxies in the Local Group, i.e. the Milky Way, Andromeda, and Triangulum (M33), we use published data of three well-studied galaxies outside the Local Group, namely M81, Centaurus A, and M101, and their confirmed satellites, and we find a strong and significant correlation between N_sat and B/T. This presents itself in contradiction with the hitherto published results from cosmological simulations reporting an absence of a correlation between N_sat and B/T in the LCDM model. We conclude that, based on the current data, the N_sat vs. B/T correlation is no longer a property confined to only the Local Group.

preprint2020arXiv

Chemical evolution of ultra-faint dwarf galaxies in the self-consistently calculated IGIMF theory

The galaxy-wide stellar initial mass function (gwIMF) of a galaxy in dependence of its metallicity and star formation rate (SFR) can be calculated by the integrated galactic IMF (IGIMF) theory. Lacchin et al. (2019) apply the IGIMF theory for the first time to study the chemical evolution of the ultra-faint dwarf (UFD) satellite galaxies and failed to reproduce the data. Here, we find that the IGIMF theory is naturally consistent with the data. We apply the time-evolving gwIMF calculated at each timestep. The number of type Ia supernova explosions per unit stellar mass formed is renormalised according to the gwIMF. The chemical evolution of Boötes I, one of the best observed UFD, is calculated. Our calculation suggests a mildly bottom-light and top-light gwIMF for Boötes I, and that this UFD has the same gas-consumption timescale as other dwarfs but was quenched about 0.1 Gyr after formation, being consistent with independent estimations and similar to Dragonfly 44. The recovered best fitting input parameters in this work are not covered in the work of Lacchin et al. (2019), creating the discrepancy between our conclusions. In addition, a detailed discussion of uncertainties is presented addressing how the results of chemical evolution models depend on applied assumptions. This study demonstrates the power of the IGIMF theory in understanding the star-formation in extreme environments and shows that UDFs are a promising pathway to constrain the variation of the low-mass stellar IMF.

preprint2020arXiv

Constraints on the star formation histories of galaxies in the Local Cosmological Volume

The majority of galaxies with current star-formation rates (SFRs), SFRo >= 10^-3 Msun/yr, in the Local Cosmological Volume where observations should be reliable, have the property that their observed SFRo is larger than their average star formation rate. This is in tension with the evolution of galaxies described by delayed-tau models, according to which the opposite would be expected. The tension is apparent in that local galaxies imply the star formation timescale tau approx 6.7 Gyr, much longer than the 3.5-4.5 Gyr obtained using an empirically determined main sequence at several redshifts. Using models where the SFR is a power law in time of the form propto (t - t1)^eta for t1 = 1.8 Gyr (with no stars forming prior to t1) implies that eta = 0.18 +- 0.03. This suggested near-constancy of a galaxy's SFR over time raises non-trivial problems for the evolution and formation time of galaxies, but is broadly consistent with the observed decreasing main sequence with increasing age of the Universe.

preprint2020arXiv

Scale-invariant dynamics in the Solar System

The covariant scale-invariant dynamics (SID) theory has recently been proposed as a possible explanation for the observed dynamical discrepancies in galaxies (Maeder & Gueorguiev 2020). SID implies that these discrepancies $-$ commonly attributed to dark matter $-$ arise instead from a non-standard velocity-dependent force that causes two-body near-Keplerian orbits to expand. We show that the predicted expansion of the Earth-Moon orbit is incompatible with lunar laser ranging data at $>200σ$. Moreover, SID predicts that the gravitating mass of any object was much smaller in the past. If true, a low-mass red giant star must be significantly older than in standard theory. This would make it much older than the conventional age of the Universe, which however is expected to be similarly old in SID. Moreover, it is not completely clear whether SID truly contains new physics beyond General Relativity, with several previous works arguing that the extra degree of freedom is purely mathematical. We conclude that the SID model is falsified at high significance by observations across a range of scales, even if it is theoretically well formulated.

preprint2020arXiv

Solar System limits on gravitational dipoles

The gravitational dipole theory of Hadjukovic (2010) is based on the hypothesis that antimatter has a negative gravitational mass and thus falls upwards on Earth. Astrophysically, the model is similar to but more fundamental than Modified Newtonian Dynamics (MOND), with the Newtonian gravity $g_{_N}$ towards an isolated point mass boosted by the factor $ν= 1 + \left( α/x \right) \tanh \left( \sqrt{x}/α\right)$, where $x \equiv g_{_N}/a_{_0}$ and $a_{_0} = 1.2 \times 10^{-10}$ m/s$^2$ is the MOND acceleration constant. We show that $α$ must lie in the range ${0.4-1}$ to acceptably fit galaxy rotation curves. In the Solar System, this interpolating function implies an extra Sunwards acceleration of ${αa_{_0}}$. This would cause Saturn to deviate from Newtonian expectations by ${7000 \left( α/0.4 \right)}$ km over 15 years, starting from known initial position and velocity on a near-circular orbit. We demonstrate that this prediction should not be significantly altered by the postulated dipole haloes of other planets due to the rather small region in which each planet's gravity dominates over that of the Sun. The orbit of Saturn should similarly be little affected by a possible ninth planet in the outer Solar System and by the Galactic gravity causing a non-spherical distribution of gravitational dipoles several kAU from the Sun. Radio tracking of the Cassini spacecraft orbiting Saturn yields a ${5σ}$ upper limit of 160 metres on deviations from its conventionally calculated trajectory. These measurements imply a much more stringent upper limit on $α$ than the minimum required for consistency with rotation curve data. Therefore, no value of $α$ can simultaneously match all available constraints, falsifying the gravitational dipole theory in its current form at extremely high significance.

preprint2020arXiv

The failure of testing for cosmic opacity via the distance-duality relation

The distance-duality relation (DDR) between the luminosity distance $D_L$ and the angular diameter distance $D_A$ is viewed as a powerful tool for testing for the opacity of the Universe, being independent of any cosmological model. It was applied by many authors, who mostly confirm its validity and report a negligible opacity of the Universe. Nevertheless, a thorough analysis reveals that applying the DDR in cosmic opacity tests is tricky. Its applicability is strongly limited because of a non-unique interpretation of the $D_L$ data in terms of cosmic opacity and a rather low accuracy and deficient extent of currently available $D_A$ data. Moreover, authors usually assume that cosmic opacity is frequency independent and parametrize it in their tests by a prescribed phenomenological function. In this way, they only prove that cosmic opacity does not follow their assumptions. As a consequence, no convincing evidence of transparency of the universe using the DDR has so far been presented.

preprint2020arXiv

The formation of exponential disk galaxies in MOND

The formation and evolution of galaxies is highly dependent on the dynamics of stars and gas, which is governed by the underlying law of gravity. To investigate how the formation and evolution of galaxies takes place in Milgromian gravity (MOND), we present full hydrodynamical simulations with the Phantom of Ramses (POR) code. These are the first-ever galaxy formation simulations done in MOND with detailed hydrodynamics, including star formation, stellar feedback, radiative transfer and supernovae. These models start from simplified initial conditions, in the form of isolated, rotating gas spheres in the early Universe. These collapse and form late-type galaxies obeying several scaling relations, which was not a priori expected. The formed galaxies have a compact bulge and a disk with exponentially decreasing surface mass density profiles and scale lengths consistent with observed galaxies, and vertical stellar mass distributions with distinct exponential profiles (thin and thick disk). This work thus shows for the first time that disk galaxies with exponential profiles in both gas and stars are a generic outcome of collapsing gas clouds in MOND. These models have a slight lack of stellar angular momentum because of their somewhat compact stellar bulge, which is connected to the simple initial conditions and the negligible later gas accretion. We also analyse how the addition of more complex baryonic physics changes the main resulting properties of the models and find this to be negligibly so in the Milgromian framework.

preprint2020arXiv

Tidal tails of open star clusters as probes to early gas expulsion I: A semi-analytic model

Star clusters form out of the densest parts of infrared dark clouds. The emergence of massive stars expels the residual gas, which has not formed stars yet. Gas expulsion lowers the gravitational potential of the embedded cluster, unbinding many of the cluster stars. These stars then move on their own trajectories in the external gravitational field of the Galaxy, forming a tidal tail. We investigate the formation and evolution of the tidal tail forming due to expulsion of primordial gas under various scenarios of gas expulsion to provide predictions for tidal tails around dynamically evolved (age > 100 Myr) galactic star clusters, which can be possibly detected by the Gaia mission. We provide a semi-analytical model for the tail evolution. We find that tidal tails released during gas expulsion have different kinematic properties than the tails gradually forming due to evaporation. The gas expulsion tidal tail shows non-monotonic expansion with time, where longer epochs of expansion are interspersed with shorter epochs of contraction. The tail thickness and velocity dispersions strongly, but not exactly periodically, vary with time. The times of minima of tail thickness and velocity dispersions are given only by the properties of the galactic potential, and not by the properties of the cluster. The estimates provided by the (semi-)analytical model for the extent of the tail, the minima of tail thickness, and velocity dispersions are in a very good agreement with the nbody6 simulations. This implies that the semi-analytic model can be used for estimating the properties of the gas expulsion tidal tail for a cluster of a given age and orbital parameters without the necessity of performing numerical simulations. A study with a more extended parameter space of the initial conditions is performed in the follow up paper.

preprint2020arXiv

Tidal tails of open star clusters as probes to early gas expulsion II: Predictions for Gaia

We study the formation and evolution of the tidal tail released from a young Pleiades-like star cluster, due to expulsion of primordial gas in a realistic gravitational field of the Galaxy. The tidal tails (as well as clusters) are integrated from their embedded phase for 300 Myr. We vary star formation efficiencies (SFEs) from 33% to 100% and the timescales of gas expulsion as free parameters, and provide predictions for the morphology and kinematics of the evolved tail for each of the models. The resulting tail properties are intended for comparison with anticipated Gaia observations in order to constrain the poorly understood early conditions during the gas phase and gas expulsion. The simulations are performed with the code Nbody6 including a realistic external gravitational potential of the Galaxy, and an analytical approximation for the natal gaseous potential. Assuming that the Pleiades formed with rapid gas expulsion and an SFE of $\approx 30$%, the current Pleiades are surrounded by a rich tail extending from $\approx 150$ to $\approx 350$ pc from the cluster and containing 0.7 to 2.7 times the number of stars in the present-day cluster. If the Pleiades formed with an SFE close to 100%, then the tail is shorter ($\lesssim 90$ pc) and substantially poorer, containing only $\approx 2$% of the number of present-day cluster stars. If the Pleiades formed with an SFE of $\approx 30$%, but the gas expulsion was adiabatic, the tail signatures are indistinguishable from the case of the model with 100% SFE. The model takes into account the estimated contamination due to the field stars and the Hyades-Pleiades stream, which constitutes a more limiting factor than the accuracy of the Gaia measurements.

preprint2020arXiv

Very high redshift quasars and the rapid emergence of super-massive black holes

The observation of quasars at very high redshift such as Poniuaena is a challenge for models of super-massive black hole (SMBH) formation. This work presents a study of SMBH formation via known physical processes in star-burst clusters formed at the onset of the formation of their hosting galaxy. While at the early stages hyper-massive star-burst clusters reach the luminosities of quasars, once their massive stars die, the ensuing gas accretion from the still forming host galaxy compresses its stellar black hole (BH) component to a compact state overcoming heating from the BH--BH binaries such that the cluster collapses, forming a massive SMBH-seed within about a hundred Myr. Within this scenario the SMBH--spheroid correlation emerges near-to-exactly. The highest-redshift quasars may thus be hyper-massive star-burst clusters or young ultra-compact dwarf galaxies (UCDs), being the precursors of the SMBHs that form therein within about 200 Myr of the first stars. For spheroid masses <10^9.6 Msun a SMBH cannot form and instead only the accumulated nuclear cluster remains. The number evolution of the quasar phases with redshift is calculated and the possible problem of missing quasars at very high redshift is raised. SMBH-bearing UCDs and the formation of spheroids are discussed critically in view of the high redshift observations. A possible tension is found between the high star-formation rates (SFRs) implied by downsizing and the observed SFRs, which may be alleviated within the IGIMF theory and if the downsizing times are somewhat longer.

preprint2019arXiv

The Milky Way&#39;s Disk of Classical Satellite Galaxies in Light of Gaia DR2

We study the correlation of orbital poles of the 11 classical satellite galaxies of the Milky Way, comparing results from previous proper motions with the independent data by Gaia DR2. Previous results on the degree of correlation and its significance are confirmed by the new data. A majority of the satellites co-orbit along the Vast Polar Structure, the plane (or disk) of satellite galaxies defined by their positions. The orbital planes of eight satellites align to $<20^\circ$ with a common direction, seven even orbit in the same sense. Most also share similar specific angular momenta, though their wide distribution on the sky does not support a recent group infall or satellites-of-satellites origin. The orbital pole concentration has continuously increased as more precise proper motions were measured, as expected if the underlying distribution shows true correlation that is washed out by observational uncertainties. The orbital poles of the up to seven most correlated satellites are in fact almost as concentrated as expected for the best-possible orbital alignment achievable given the satellite positions. Combining the best-available proper motions substantially increases the tension with $Λ$CDM cosmological expectations: <0.1 per cent of simulated satellite systems in IllustrisTNG contain seven orbital poles as closely aligned as observed. Simulated systems that simultaneously reproduce the concentration of orbital poles and the flattening of the satellite distribution have a frequency of <0.1 per cent for any number of k > 3 combined orbital poles, indicating that these results are not affected by a look-elsewhere effect. This compounds the Planes of Satellite Galaxies Problem.

preprint2019arXiv

The systematically varying stellar IMF

Some ultra-compact dwarf galaxies have large dynamical mass to light (M/L) ratios and also appear to contain an overabundance of LMXB sources, and some Milky Way globular clusters have a low concentration and appear to have a deficit of low-mass stars. These observations can be explained if the stellar IMF becomes increasingly top-heavy with decreasing metallicity and increasing gas density of the forming object. The thus constrained stellar IMF then accounts for the observed trend of metallicity and M/L ratio found amongst M31 globular star clusters. It also accounts for the overall shift of the observationally deduced galaxy-wide IMF from top-light to top-heavy with increasing star formation rate amongst galaxies. If the IMF varies similarly to deduced here, then extremely young very massive star-burst clusters observed at a high redshift would appear quasar-like (Jerabkova et al. 2017) .

preprint2019arXiv

The ultra-diffuse dwarf galaxies NGC 1052-DF2 and 1052-DF4 are in conflict with standard cosmology

Recently van Dokkum et al. (2018b) reported that the galaxy NGC 1052-DF2 (DF2) lacks dark matter if located at $20$ Mpc from Earth. In contrast, DF2 is a dark-matter-dominated dwarf galaxy with a normal globular cluster population if it has a much shorter distance near $10$ Mpc. However, DF2 then has a high peculiar velocity wrt. the cosmic microwave background of $886$ $\rm{km\,s^{-1}}$, which differs from that of the Local Group (LG) velocity vector by $1298$ $\rm{km\,s^{-1}}$ with an angle of $117 \, ^{\circ}$. Taking into account the dynamical $M/L$ ratio, the stellar mass, half-light radius, peculiar velocity, motion relative to the LG, and the luminosities of the globular clusters, we show that the probability of finding DF2-like galaxies in the lambda cold dark matter ($Λ$CDM) TNG100-1 simulation is at most $1.0\times10^{-4}$ at $11.5$ Mpc and is $4.8\times10^{-7}$ at $20.0$ Mpc. At $11.5$ Mpc, the peculiar velocity is in significant tension with the TNG100-1, TNG300-1, and Millennium simulations, but occurs naturally in a Milgromian cosmology. At $20.0$ Mpc, the unusual globular cluster population would challenge any cosmological model. Estimating that precise measurements of the internal velocity dispersion, stellar mass, and distance exist for $100$ galaxies, DF2 is in $2.6σ$ ($11.5$ Mpc) and $4.1σ$ ($20.0$ Mpc) tension with standard cosmology. Adopting the former distance for DF2 and assuming that NGC 1052-DF4 is at $20.0$ Mpc, the existence of both is in tension at $\geq4.8σ$ with the $Λ$CDM model. If both galaxies are at $20.0$ Mpc the $Λ$CDM cosmology has to be rejected by $\geq5.8σ$.