Heavy metal rules. I. Exoplanet incidence and metallicity

Discovery of only handful of exoplanets required to establish a correlation between giant planet occurrence and metallicity of their host stars. More than 20 years have already passed from that discovery, however, many questions are still under lively debate: What is the origin of that relation? what is the exact functional form of the giant planet -- metallicity relation (in the metal-poor regime)?, does such a relation exist for terrestrial planets? All these question are very important for our understanding of the formation and evolution of (exo)planets of different types around different types of stars and are subject of the present manuscript. Besides making a comprehensive literature review about the role of metallicity on the formation of exoplanets, I also revisited most of the planet -- metallicity related correlations reported in the literature using a large and homogeneous data provided by the SWEET-Cat catalog. This study lead to several new results and conclusions, two of which I believe deserve to be highlighted in the abstract: i) The hosts of sub-Jupiter mass planets ($\sim$0.6 -- 0.9~M$_{\jupiter}$) are systematically less metallic than the hosts of Jupiter-mass planets. This result might be related to the longer disk lifetime and higher amount of planet building materials available at high metallicities, which allow a formation of more massive Jupiter-like planets. ii) Contrary to the previous claims, our data and results do not support the existence of a breakpoint planetary mass at 4~M$_{\jupiter}$ above and below which planet formation channels are different. However, the results also suggest that planets of the same (high) mass can be formed through different channels depending on the (disk) stellar mass i.e. environmental conditions.


Introduction
Understanding the origin of Earth and life on it was one of the most important and daring questions since ancient times. These questions still did not loose their actuality. Predictions about the existence of planets orbiting other stars, currently called extrasolar planets or exoplanets, have been expressed already for several centuries [e.g. 1]. Moreover, the RV and transiting methods, the two methods that are currently the most successful ones in terms of planet detection counts and are responsible for about 95% of all the known exoplanets, were proposed to detect these planets around other stars already 80 years ago [2,3].
Starting from 1980s, thanks to the developments of high-precision spectrographs and CCD cameras, the search for exoplanets received a new impulse. Given the predictions that planets outside of our Solar System should exist numerously and the intensive search for these extrasolar planets, the detection of a planet around other stars should not have came as a surprise. However, it came: the first exoplanets 1 were discovered around a pulsar PSR1257 + 12 [4]. The existence of planets around neutron stars made a very difficult question to answer: i) did these planets form around a massive precursor star and survived the supernova explosion, ii) or they formed in a protoplanetary disk which was left over immediately after the neutron star was formed, iii) or these planets were formed from a disk consisting of an already disrupted binary companion [see 6, for different formation scenarios]. The first confident detection of an exoplanet around a solar-type stars was made in 1995: Michel Mayor and Didier Queloz discovered a giant planet orbiting its host star with a period of only about 4 days [7]. The discovery of such a planet, which does not have an analog in our Solar System, rose more questions about formation of planets in and out of our system.

Diversity of exoplanet properties and their formation scenarios
Since the discoveries of first exoplanets via radial velocity [51 Pegasi b 7], transiting method [HD209458 b 8,9], and by direct/indirect imaging [2M1207 b/Fomalhaut b 10,11] 2 almost 4000 planets have been detected 3 and more than 4000 candidates awaiting for validation [12]. Right panel of Fig. 1 shows the masses (for planets detected by RV method only the minimum mass is available) of the so far discovered planets as a function of discovery time labeled according to detection technique. The data is extracted from Extrasolar Planets Encyclopaedia [13]. Thanks to the increasing number of ground based and space based planet search surveys the number of exoplanets has increased significantly during the last years. The plot also shows that the detection of planets with masses as low as the one of Earth became possible. 10 1 10 0 10 1 10 2 10 3 10 4 10 5 10 6 Period (day)  Mass of the known planets as a function of the discovery year (right). The planets detected by different detection techniques are shown by different symbols. The planets of our Solar System are shown as a reference. The two low-mass and short-period planets detected by direct imaging, Kepler-70b and Kepler-70c, are found to orbit a post-red-giant star [14].
Left panel of Fig. 1 the distribution of exoplanets in the period -mass diagram is plotted. The planets of our Solar System are also shown for a visual comparison. This figure suffers from selection effects and observational biases of different planet detection techniques. However, it still contains 1 Two planets were first discovered around this neutron star in 1992 [4] and one additional planet was detected later in 1994 [5]. 2 The "List of exoplanet firsts" can be found from https://en.wikipedia.org/wiki/List_of_exoplanet_firsts 3 http://exoplanet.eu/ very valuable information about the processes of formation and evolution of the detected planets. In particular one can easily see that most of the exoplanets are clustered in three groups: i) hot-Jupiters with M p ∼1-2 M , ii) gas and ice giants with M p ∼1-2 M and P ∼1000 days, and iii) hot super-Earths with M p ∼10 M ⊕ and P < 100 days. The observed lack of giant planets with periods between 10-100 days, so called "period-valley" [15,16] is probably related to the disk migration [see 17, for a discussion] and the relative paucity of intermediate mass planets with masses between 10 to 100 M ⊕ (so called "planetary desert" [e.g. 18]) 4 is explained as a consequence of the rapid increase in mass through runaway accretion when the planet core and envelope mass reach to about 30M ⊕ [e.g. 18,21]. Note, that the recent CA population synthesis of Mordasini [22] suggest not as strong "planetary desert" as was previously proposed [e.g . 18]. The reason that the planets of our SS occupy the empty regions of the diagram is partly due to detection biases. At the same time, in our SS there are no short period (shorter than about 100 days) super-Earths 5 , which are very common around other stars. This property makes our SS somehow special [24] 6 . For more details about the distribution of planets on the periodmass diagram and on the architecture of exoplanetary systems I refer the reader to the recent excellent review by Winn and Fabrycky [26]. Radius (R ) Figure 2. The mass-radius relation for the discovered exoplanets (exoplanet.eu) with mass and radius determination better than 30%. The planets of our Solar System are shown for comparison. 4 The recent surveys of microlensing planets suggest about 10 times more intermediate-mass (20 -80 M ⊕ for a median host star mass of 0.6 M ) planets than the CA models of Ida and Lin [18] and Mordasini et al. [19] predict [20]. 5 I note that using the term 'super-Earth' can be misleading or confusing for a reader (especially for a reader outside of the field or for public audience) [23]. That is not our intention. By super-Earth I mean (as almost everywhere in the literature) planets more massive or larger than the Earth. I do not imply that more information about the properties of these planets are known that makes them similar to the Earth. 6 It was also proposed that Jupiter's periastron is atypical making the SS an outlier [e.g. 25]. However, the Jupiter's 'atypical' periastron is probably due to the observational biases [24]. Fig. 2 shows the mass-radius diagram of the extrasolar planets with precise measurements of mass and radius better than 30%. Similar to Fig. 1 I see large diversity of the physical characteristics of these exoplanets. In particular it is apparent the large number of Jupiter mass planets with radii much larger than that of Jupiter. Most of these bloated planets have very short orbits receiving incident flux of 2-3 magnitudes higher than the Earth receives [27,28]. Several theoretical explanations have been proposed [29][30][31] that can efficiently explain the radii of individual planets, but they have difficulties to explain the properties of the entire population of these puffed-up planets. The recent work by Sestovic et al. [32] suggests that the inflation extent depends on the mass of the planets. In the figure it is also clear that there is a transition mass at which the mass-radius relation changes its functional form. Statistical studies of the mass-radius relation can provide valuable information about the bulk composition and structure of these planets, and even information about their formation and evolution [33][34][35].   [36] and TD/GI [37] planet formation models against the metallicity of their host stars. The running mean of [Fe/H] as a function of mass of CA planets (in red), TD/GI planets (in blue), and observed planets (in black) is shown in the top panel. The running means were calculated using windows of 50, 500, and 1000 for observed, TD/GI, and CA planets, respectively. These numbers reflect the sizes of each sample.
The radius-mass ( Fig. 1) and period-mass ( Fig. 2) diagrams (among other parameter spaces) contain enormous information about the physical processes acting during the formation and evolution of (exo)planets. Thus the reproduction of these diagrams is the ultimate goal of any planet formation theory. The models developed to explain the formation of SS planets did not predict such an enormous diversity of exoplanets that we observe in Figs. 1 and 2 [e.g. 38,39]. The two most commonly used exoplanet formation scenarios are the core accretion [CA; e.g. 18,19,[40][41][42][43] and gravitational instability [GI; e.g. [44][45][46][47][48]. In the CA scenario the low-mass planets and cores of massive planets form from the coagulation of small grains in the protoplanetary disc. When these cores reach to a critical mass of about 5-10M ⊕ [49][50][51] they undergo runaway accretion of gas. Giant planets then will be formed if the critical core mass is reached before the dissipation of the protoplanetary disk [e.g. 52], which has a typical lifetime of several Myr [53,54]. Since the phase of runaway accretion is very fast the planet formation time is essentially defined by the time of core-formation [e.g . 42]. In the GI scenario, the massive and cold disks fragment into a few Jupiter mass clump which then contracts to form giant planets [e.g. 55]. These massive planets are formed quickly, much before the gas in the disk depletes [see 56, and references therein]. If the GI is followed by tidal stripping during the migration of the clumps, then even a rocky planet (including objects in the SS [57,58]), can be formed [see 59, for a recent review on "Tidal Downsizing"].
Both the classical CA and GI theories have experienced substantial development and modifications during the last decades [e.g. [60][61][62][63]. They include important processes such as pebble accretion [e.g. [64][65][66] and/or migration in the disk [e.g. 67]. Both theories have their strengths and weaknesses [68] and might operate efficiently under different physical conditions and different parameter space [69]. The planetary population synthesis calculations [18,[70][71][72]  In the bottom panel of Fig. 3 I plot the mass of synthetic planets formed through CA of Mordasini et al. [36] 7 and through TD/GI of Nayakshin [37] 8 against the metallicity of their host stars. The plot shows that the two conceptually different models cover very similar areas in this diagram. On the top panel of the figure I plot running mean of the host metallicity as a function of planetary mass for CA synthetic planets, TD/GI synthetic planet, and confirmed planets with hosts metallicities derived in a homogeneous way in SWEET-Cat (see Sect. 2.1). This panel clearly shows that when one looks at the populations of planets as a whole, the two models show significant deviation from each other and from the observed planet populations in some metallicity regions. Although in the two considered models the absolute values of the mean [Fe/H] may change with population synthesis parameters, the trends should still remain since they are based on solid physical principles. Thus testing the predictions of the planet formation models with observations is the way to identify the weaknesses in the models and to provide insightful information for their further development.

Motivation and the outline of this review
As discussed in the previous section, the modern planet formation theories need to consider many different physical phenomena (e.g. migration of planets, evaporation, pebble accretion) in order to reproduce the main properties of the detected exoplanets. Besides these physical processes that directly affect the characteristics of the planets and their orbital architecture, one should also consider the environmental conditions where these planets are formed [e.g. 75].
Linking the properties of exoplanets with the properties of protoplanetary disks where they have been formed will provide important insights about planet formation in different environments. Studying the dependencies of exoplanets properties on environmental conditions are usually performed indirectly by looking at the physical properties of the host stars that are linked to the characteristics of the proto-stellar/planetary disk. Since the planets and their host stars are formed from the same molecular clouds, some global environmental properties such us chemical composition of some important mineralogical abundance ratios [e.g. Mg/Si, Fe/Si 76-78] 9 , can be determined from the spectroscopic studies of their host stars. Moreover, precise characterization of the planet host stars is crucial for the characterization of the planets themselves [81].
In this manuscript I will review the main observational correlations between the properties of exoplanets and the metallicity of their host stars. I will discuss the most-studied correlation between 7 The data is downloaded from http://www.mpia.de/~mordasini/Site7.html. 8 The data is kindly provided by Sergei Nayakshin. 9 Note that the C/O ratio, that controls the amount of the carbides and silicates in planet building blocks, varies in the protoplanetary disk and might be different from the value of the parent star [79,80]. occurrence rate of different types of planets with metallicity (iron content) of their host stars and its link to different planet formation processes and models. Sine some of the correlations reported in the literature are based on heterogeneous data (for example compilation of properties of hosts stars from different literature sources), whenever possible, I will revisit them using the largest homogeneous catalog of exoplanet host stars called SWEET-Cat [see Sect. 2.1, 82].

SWEET-Cat
There are four exoplanet catalogs that exoplanetologists commonly use [83]. Despite small biases and discrepancies between these catalogs, on average they provide a good compilation of properties of exoplanets and their host stars [83]. All these catalogs, however, suffer from heterogeneous literature compilations of stellar parameters. This heterogeneous compilation makes significant discrepancies when comparing with parameters derived in a homogeneous way [e.g. 84]. Moreover, homogeneous and uniform spectral analysis is essential to minimize the uncertainties in stellar atmospheric properties [e.g. 85].
Santos et al. [82], presented a catalog of stellar parameters for stars with planets (SWEET-Cat 10 ) listed in the Extrasolar Planets Encyclopaedia [exoplanet. eu 13]. This catalog provides a compilation of atmospheric stellar parameters (see Fig. 4) from literature and, whenever possible, derived using the same uniform methodology [see e.g. [86][87][88]. The catalog is regularly updated [84,89] and after the last major update [90] it contains uniformly derived stellar parameters of ∼80% of bright stars (V <12 mag) hosting RV detected planets. Besides the main stellar parameters, SWEET-Cat will soon provide abundance of different chemical abundances of the host stars again derived in a homogeneous and uniform way [Sousa et al. in prep.]. It is worth to note that different groups also intensively work on homogeneous derivation of stellar parameters of Kepler stars hosting planets and planet candidates e.g. The Kepler Follow-up Observation Program [91] and The California-Kepler Survey [92]. Moreover, during the past few years dedicated communities 11 and web interfaces 12 have been created to optimize the resources in exoplanet follow-up studies and characterization of their host stars.

Nomenclature: the term 'metallicity' in astronomy
The term 'metallicity' is very widely used in the literature. Although it is always used to quantify the amount of 'metals', its nomenclature and scales are not the same in different fields of astronomy. In this subsection I will very briefly present the definitions of metallicity, its scales, and the quantities that are derived in practice as a proxy of metallicity in different fields of stellar astronomy.

Mass fraction
In astronomy, all the chemical elements heavier than hydrogen and helium are called metals 13 . The composition of the stars can be characterized by the mass fraction of hydrogen (denoted as X), helium (denoted as Y), and metals (denoted as Z). The sum of these three abundances is normalized to unity i.e. X + Y + Z = 1. This definition is convenient when dealing with the models of stars. However, direct derivations of these quantities from observations is impossible. Note, that in the literature Z is sometimes called "metal mass fraction" [e.g. 93] and sometimes "metallicity" [e.g. 94].
The solar (log ε(H) ≡ 12) scale Atomic abundances of different metals (elements with atomic number greater than 2) in stars vary over several magnitudes, and thus, it is convenient to express them in a logarithmic scale. In stellar astrophysics, the atomic abundances of elements are usually presented relative to H 14 . Traditionally the logarithmic abundance (number of atoms) of H is set to 12 i.e. log ε(H) ≡ 12. Atomic abundance of any other element El, then will be where N El and N H are the number of atoms of element X and H (log(N H ) = 12). The unit used for the abundances is the "dex" which is contraction of "decimal exponent".
Sometimes it is more convenient to derive the atomic abundances of stars relative to Sun using the log ε(H) ≡ 12 scale. In this case square braces are used to express the abundance of an element El The conversion from solar scaled to absolute abundances requires the knowledge of the reference solar abundances [e.g. [95][96][97][98].
Metallicity of a star [m/H] is related to the metal mass fraction Z in the following form Unfortunately, it is very difficult to derive atomic abundances of all the metals for a given star. Thus, in stellar astrophysics, the iron content ([Fe/H]) is commonly used as a tracer of overall stellar metallicity. Although iron is not the most abundant metal in the Universe, the visible spectra of solar-type stars contain many strong iron lines which are easy the measure 15 . The usage of iron abundance as a proxy of overall meallicity assumes that the abundances of all the metals change proportionally to iron content. Thus under the assumption that the distribution of heavy metals in a given star is the same as in the  The assumption about the universality of solar-scaled heavy element distribution is not always valid, especially for metal-poor stars that are enhanced in α elements [e.g. [101][102][103][104][105]. To account for the α-enrichment the following correction is proposed [106] [ For a recent review on the derivation of stellar metallicities and other spectroscopic parameters of exoplanet hosting stars I refer the reader to Adibekyan et al. [81].

Metallicity of planet host stars
The first observed correlation that linked the presence of planets and a property of their host stars was the giant planet -stellar metallicity correlation [e.g. 107,108]. The discovery of this correlation played a crucial role for the advancement of planet formation theories [e.g. 18,36,59]. Until now, studies of the possible correlations between different types of planets (e.g. low-mass or massive planets at short or long period orbits) and metallicities of stars of different spectral types and evolutionary stages intensively continue, and provide new constraints for the models of planet formation and evolution. Planet-metallicity correlation is the subject of the discussion of this section.
It is very important to note, that in this manuscript the iron content is used as a proxy of overall metallicity 16 . This approximation is very commonly used in the literature and is well justified for solar metallicity stars. However, as already mentioned, at low metallicities the iron abundance does not necessarily equal to the overall metal content in stellar atmospheres. There is another limitation that is typically ignored when studying the metallicity distribution of planet hosting stars. There are some astrophysical processes, such as atomic diffusion due to concentration and thermal gradients [e.g. 109], that can affect the stellar metallicity. Moreover, engulfment of planetisimals under specific conditions [e.g. depending on the number of accretion episodes 110] and/or non-accretion of metal-rich material as a consequence of planet formation [again, depending on the time of planet formation and the composition of the formed planets 111] may also change the original metallicity of the planet hosting stars. Thus, the present-day stellar metallicity can be slightly different from the metallicity of the stars at the time of their formation [e.g. 96].. Finally, it is important to note, that the key parameter in the core-accretion models of planet formation is the solid surface density [e.g. 18,19]. This parameter is proportional to the dust-to-gas ratio in the disk and is very difficult to constrain by observations. In turn, average dust-to-gas ratio in the protoplanetary disks is proportional to the metal abundance of the ISM [e.g. 112] and is usually assumed to be proportional to the atmospheric metallicity of the host stars [e.g. 19,113,114]. The validity of this assumption was very recently revisited by Liu et al. [115]. The authors concluded that although the metallicity of the disk can change significantly as the disks evolve, the dust-to-gas ratio (within an order of magnitude) is comparable with the heavy element content in the atmosphere of the host stars. However, for particular cases and in particular locations the dust-to-gas ratio can be significantly different from the host star metallicity [e.g. 116] because of particles being decoupled from gas and radially drifting inwards [e.g. 117,118].

Giant planets and metallicity
Detection of only four giant planets required from Gonzalez [107] to notice that these planets tend to appear around metal-rich stars. This hint was very soon confirmed using relatively larger samples [e.g. 108,119,120] and become a well established correlation [e.g. 86,121-124] known as giant planetmetallicity correlation.
At the end of 1990s and beginning of 2000s, when only handful of hot-Jupiters have been detected, two hypotheses have been proposed to explain the metal excess of massive exoplanet host stars: self enrichment (aka pollution) mechanism and primordial origin [e.g. 107]. In the self enrichment scenario the high metallicity of planet hosts was explained by the pollution of outer convective envelope of the hosts due to accretion of gas-depleted material. This pollution could be a result of inward migration of massive planets that sweep metal-rich material with them. In the alternative, primordial scenario, the over-abundance of metals in massive planet hosts was considered to represent the high metallicity of the primordial cloud where the stars have been formed and was assumed that giant planets form more readily in high metallicity environments.
The two aforementioned scenarios provide different observationally testable signatures for the host stars. In particular, in case of pollution by rocky material, the increase of atmospheric metallicity of the host stars should depend on the mass of convective zone, thus on the mass and temperature of the stars [107,113,125,126]. Going one step further, Gonzalez [107] proposed that if the accreting material is fractionated then a trend between chemical abundances and condensation temperature of the elements is expected. Gonzalez [107] also noted that the high-metallicity accretion signature should be lost for evolved stars which have larger convective zones where the accreted heavy metals will be mixed and diluted. When the number of detected planets increased, several systematic abundance studies on stars with and without planets have been performed focused on the aforementioned predictions. The results of these studies favored the primordial cloud as the most likely origin for the metal-rich nature of giant planet host stars [e.g. 86,121,122,127]. This hypothesis is supported by models of planet formation and evolution based on CA [e.g. 18,36] and tidal downsizing [63].
In Fig. 5 I plot the metallicity distribution of stars hosting high mass planets (hereafter, HMPH) and stars known to host no detected planet (hereafter, SnoP). Only FGK dwarf hosts (log g > 4.0 dex, and 0.6 M < M < 1.5 M ) stars with with stellar parameters derived in a homogeneous way in SWEET-Cat are considered. The comparison sample of stars without detected planets consists of 954 FGK dwarf stars observed within the HARPS GTO program and are taken from Adibekyan et al. [103]. In this figure only planets with masses between 50 M ⊕ and 13 M 17 are considered. In case of multiplanetary systems the mass of the most massive planet is considered. I note that there is no clear definition for what should be the low mass limit of giant planets, thus the choice for lower limit of mass M = 50 M ⊕ is somewhat semi-arbitrary. Different authors chose different definitions 18 for giant planets and the choice of the low-mass limit varies (typically between 30 to 90 M ⊕ ) depending on the aim of the study [e.g. 35,90,[131][132][133][134][135]. In any case I choose this limit since it is close to the mass 'gap' observed in Fig. 1 and 2. At the same time the planet core-accretion based population synthesis models by Mordasini et al. [19] predicted such a minimum in the mass-distribution at about this mass. It is important to note that presented results are not sensitive to the choice of these limits. Here I 17 13M approximately corresponds to the deuterium-burning mass limit for BDs at solar metallicity [e.g. 128] and is usually used as a border line between giant planets and BDs. The exact deuterium-burning mass limit depends on several parameters [see e.g. 129,130, for more details]. 18 I refer the readers to Chabrier et al. [129] and Schneider et al. [13] for an interesting discussion about the classifications of planets and BDs based on physical observable properties and formation mechanisms.
should also comment that the classifications of different types of planets based on their radii is not straightforward as well, and varies from work to work [e.g. 136-140]. Although some theoretical works based on CA models propose that the critical metallicity for giant planet formation is in the range of −1.8 to −1.5 dex 19 [50,142], the lowest metallicity of the giant planet host in Fig. 5 is -0.65 dex. In fact, the lowest metallicity of a confirmed sub-stellar companion host currently listed in exoplanet.eu is −1.00±0.07 dex derived for BD+202457 20 [143]. However, the more recent spectroscopic analysis of this star suggests a higher metallicity of [Fe/H] = −0.79 [144]. Besides BD+202457, the lowest (spectroscopically derived) metallicity star hosting a giant planet in exoplanet.eu is the 24 Boötis with a metallicity of −0.77±0.03 dex [145]. Observationally determining the metallicity limit below which giant planets do not form can provide a very important insights for the planet formation theories. Indeed, some programs are focused on the metal-poor regime trying to tackle this issue [e.g. [146][147][148].
From Fig. 5 it is clear that the SnoP and HMPH samples have different metallicity distributions, the latter ones being more metallic. In particular the mean metallicity of SnoP is −0.159±0.009 dex, while the average metallicity of the stars hosting high mass transiting and RV planets are 0.117±0.013 dex and 0.112±0.013 dex, respectively. Here the errors represent the standard error of the mean i.e. standard deviation divided by the square root of the sample size. To evaluate the significance of this difference more quantitatively I applied a two-sample KS test to the samples. The results presented in Table 1 show that stars hosting giant planets have significantly different metallicity distribution when compared to the stars without planets. Interestingly, the two samples of HMPH of transiting and RV planets have metallicity distributions that come from the same parent distribution. This is somehow surprising since the distributions of the orbital periods of the transiting and RV planets are statistically different. The transiting planets orbit their stars on average in 11 days (std = 66 days), while the RV detected planets have an average orbital period of 1202 days (std = 1775 days), and the KS tests suggests a p-value of 3.5×10 −61 for the two orbital period distributions being similar.
As already mentioned the giant planet -metallicity relation provides a general but very important information for planet formation models. The knowledge of the exact functional form of this relation can provide additional important constraints for the existing planet formation theories. However, a large, unbiased, and volume-limited sample of stars surveyed for planets is required to study the details of this relation. Usually, the dependence of planet fraction ( f ) on stellar metallicity is described by the following functional form: f = α10 β[Fe/H] , where α and β are the coefficients to be derived. To include the dependence of planet occurrence rate on stellar metallicity usually On the left panel of Fig. 6 I plot the metallicity distribution of SnoP and HMPH (again only planets 50 M ⊕ < M < 13 M are considered) taken from the HARPS volume-limited (up to about 60 pc from the Sun) survey [149]. This sample consist of 582 FGK stars, for which the stellar parameters, including the stellar metallicities, were homogeneously derived in Sousa et al. [123]. The figure shows that the mean metallicity of the HMPH sample (0.046±0.035) is again higher than the mean metallicity of the SnoP (−0.112±0.010). The performed KS test also suggest a significantly different metallicity distributions for the two samples (see Table 1). The right panel of Fig. 6 shows the dependence of giant planet frequency (number of stars with planets relative to number of all surveyed stars) on stellar metallicity. The plot clearly shows that the frequency sharply increase with metallicty. If for sub-solar metallcities about 5% of the stars have giant planets, at super-solar metallicities the relative frequency reaches to about 20%. It is important to note that this giant-planet -metallicity relation is practically unaffected by the possible dependence of RV precision on [Fe/H]. In principle, metallicity affects the strength of spectral lines, and thus RV precision. However, for the stars with [Fe/H] −0.8 dex, the RV precision is practically independent of metallicity 21 [127]. In Fig. 6 (right panel) I also show the functional form of the planet occurrence -metallicity relation derived in several works using different samples [121,122,152]. One can see that there is, at least, qualitative agreement between the fitted curves obtained by different authors for stars with metallicities above the solar value. These curves also describe well the increase of planet frequency for solar and super-solar metallicity stars observed for the HARPS sample in the figure. The main difficulty and disagreement between the different works come for sub-solar metallicities. Udry and Santos [152], based on the CORALIE sample from [86] and the sample from Fischer and Valenti [121], proposed that at sub-solar metallicities the giant planet frequency might be constant rather than exponential, which is clearly the case for solar and super-solar metallicities. This result has been contested by Johnson et al. [122] who used a larger sample of 1266 stars with and without planets from the California Planet Survey. Based on a Bayesian analysis the authors found that giant planet occurrence rate is a strong exponential function of metallicity that continues to sub-solar metallicities. More recently, Mortier et al. [124] performed a statistical Bayesian analysis on large volume-limited samples of CORALIE, HARPS 22 , and CORALIE+HARPS (the combined sample) to test between different functional forms of the planet-metallicity correlation. While confirming the general correlation between giant planet incident and metallicity, they concluded that the samples are simply too small to be able to distinguish 21 To my knowledge there is no study on the dependence of photometric precision on stellar metallicity for transiting planet surveys. However, the detectability of transiting planets slightly depends on metallicity through the relationships between transit depth, stellar radius, and stellar metallicity: transiting planets are slightly easier to detect around metal-poor dwarf stars than around metal-rich ones [150]. 22 This HARPS sample is the same as the one shown in Fig. 6. The only difference is that a few planets have been discovered around these stars after the study of Mortier et al. [124].
between a constant or an exponential form in the low-metallicity region. The authors pointed out that a sample of about 5000 stars (about a factor of 3 larger than the sample they used) is required to distinguish between different models (functional forms). Hopefully, this question will get its answer soon with the arrival of the ongoing and upcoming missions such as Gaia 23 [153], TESS [154], and PLATO [155].
Using well defined samples one can determine the occurrence rate of given type of planets around given type of stars. Occurrence rate calculations require corrections for the detectability of a given type of planets. In case of transit surveys in addition one should take into account the geometric probability for seeing a transit of a given type of planet. For a recent review on planet occurrence in RV and transit surveys I refer the reader to Winn [156]. The main properties of giant planets detected by RV and transit methods was recently reviewed by Santerne [157]. The latter author compiled a list of occurrence rate estimations of giant planet at different orbital periods: hot-Jupiters (P < 10 days), period-valley giants (10 < P < 85 days), and temperate giants (85 < P < 400 days). Their results show that the average occurrence rate increases from about 1% for the shortest period giants to about 3% for giant planets with orbital periods of several hundred days [see 157, and references therein]. It is interesting to note that the occurrence rate calculations for hot-Jupiters detected by transit method are about factor of two lower than the one calculated for RV planets [e.g. 16]. The exact reason for this discrepancy is still unclear, but can be related to properties of the stars (e.g. binary rate) monitored in RV and transit surveys and/or (in)ability of accurate characterization of these stars [e.g. 158]. The integrated occurrence rate of giant planets (M > 50 M ⊕ ) with orbital periods shorter than 10 years is estimated to be 13.9±1.7% [159].

Metallicity of sub-Jupiters
In Fig. 3, where I qualitatively compare the metallicities of different types of observed and synthetic planets, one can clearly see that there are two dips in the metallicity of observed planets at masses of ∼ 0.8 M and 2 M . In this subsection I study these dips more in detail to understand whether they have an astrophysical origin or are due to different observational biases and selection effects.
In the leftmost panel of Fig. 7 I show the dependence of planet host metallicity on the mass of planets with masses between 50 M ⊕ and 4 M . These planets are selected to be orbited around FGK dwarf stars (log g > 4.0 dex, 0.6 M < M < 1.5 M ) with stellar parameters derived in a homogeneous way in SWEET-Cat. The running means of hot-(P < 10 days) and cold-Jupiters (P > 100 days) with masses above about 1.5 M have very different behaviors (see the middle panel). This is most probably due to small number of such planets in the sample. The small number of planets with masses above ∼ 1.5 M and the change of the relative number of hot-and cold-Jupiters at 2 M are probably responsible for the dip in metallicity observed at 2 M . The number of hot-and cold-Jupiters with masses between 0.3 and 1.5 M (see the middle panel) is relatively large and their relative fractions are quite balanced. This means that the metallicity dip observed at ∼ 0.8 M might have a physical origin and deserves further exploration.
In the middle column of Fig. 7 one can see that both hot-and cold-Jupiters show a dip in metallicity at about ∼ 0.8 M , although the dip of cold Jupiters is more apparent. The main reason of separating hot and cold Jupiters is that orbital periods of giant planets correlate with the host star metallicity: long-period Jovians tend to orbit more metal-poor stars than hot-Jupiters [e.g. 138, [160][161][162]. This correlation can be seen in Fig. 7 as well, where the red squares and red curve lie above the blue circles and the blue curve. To make the description of planets inside and outside of the metallicity dip easier and shorter, the planets inside the dip will be called 'sub-Jupiters', and the planets with masses higher and lower than the sub-Jupiters will be called 'Jupiters' and 'Saturns', respectively. 23 The name "Gaia" was initially derived as an acronym for Global Astrometric Interferometer for Astrophysics. I tried to identify the lower and upper mass limits of sub-Jupiters by looking at a mass range for which the metallicity of their hosts stars is significantly lower than the metallicity of the'Jupiters' and 'Saturns'. Performed KS tests showed that the most significant difference between sub-Jupiters and non-sub-Jupiters is observed when sub-Jupiters are defined as planets with masses between about 200 M ⊕ (∼0.6 M ) and 290 M ⊕ (∼0.9 M ). However, when comparing the metallicities of the hosts of sub-Jupiters with those of the Saturns no statistically significant difference can be seen, while the difference from the hosts of Jupiters is always significant (see Table 2).  The next step was to investigate whether the sub-Jupiters (and/or their hosts) are different from the Jupiters and Saturns (and/or their hosts) in parameters that are linked to planet formation processes (e.g. orbital period, eccentricity, multiplicity, host star mass). To test whether the obtained results are sensitive to the multiplicity of the planets, I also considered the most massive planet for each planetary system (see rightmost panel of Fig. 7 and Table 2). The results did not differ significantly from the case where all the planets in each planetary system were considered. The most significant difference was observed for orbital periods of hot giant planets (Table 2). In the left panel of Fig. 8 the dependence of orbital periods on the mass of the sample of hot giant planets is plotted. The observed clear correlation is the cause of the different orbital periods observed for planets with different masses. The lack of very short period (P < 2 days) sub-Jupiters probably cannot be explained by photoevaportation as these planets loose very small fraction of their mass [e.g. [163][164][165]. Worth mentioning that this process is invoked to explain the lack of super-Earths and Neptunes at very short orbits [e.g. [166][167][168]. Most of the mechanisms proposed to explain the orbital distribution and the observed lower boundary of hot-Jupiters include disk migration [e.g. 17,169] and high-eccentricity migration [e.g. [169][170][171], both followed by tidal dissipation. Very recently, Bailey and Batygin [172] proposed that the lower boundary of the periods of hot-Jupiters can be a result of their in-situ formation followed by tidal decay for the most massive planets. Whatever is the mechanism shaping the orbital architecture of the hot-Jupiters, it is not straightforward to link it with the significant metallicitiy differences observed between the sub-Jupiter and Jupiter host stars.  179], which makes the formation of massive cores of gas giant planets difficult. At the same time, the migration of giant planets formed in these short-lived disks is supposed to be very limited [180,181]. To make the picture more complex one should remember that the disk mass also correlates with the disk accretion rate [182], a parameter which might be important for the fast formation of the cores of giant planets [115].
Interestingly, no significant difference of stellar mass was found for the subsamples of giant planets suggesting that for a given stellar mass Jupiters tend to form around more metallic stars than their lower mass counterparts. This is somehow surprising since the stellar mass correlates with the stellar metallicity as shown in the right panel of Fig. 8. This correlation is a consequence of stellar evolution, stellar age -metallicity relation, and the typical color and magnitude cuts used in the selection of target in planet search programs [121,122,183,184]. With such results at hand, one might argue that the higher metallicity of Jupiter hosts would mean longer disk lifetime [114] and higher amount of planet building material [e.g. 36]. These two parameters can directly affect the formation efficiency of giant planets and their migration rate. For cold giant planets no significant difference in orbital periods between sub-Jupiters and Jupiters is observed. This might be due to difference in starting positions of planet formation in disks with different metallicities [e.g. 36]. While in this section a very qualitative and somewhat speculative discussion was developed to explain the low metallicity of sub-Jupiter hosts when compared to the stars hosting Jupiter-mass planets, a more quantitative analysis is necessary to understand the origin of this metallicity difference. Perhaps, even before such an analysis it is necessary to study the influence of different selection and observational biases on this finding.

Very massive giants and metallicity
In this subsection I will discuss the dependence of very massive planet formation on stellar metallicity. However, such a discussion requires a knowledge of the upper mass boundary of planets, which is still an open question. As mentioned in Sect. 3.1, depending on the goal of the study, different authors use different lower and upper mass limits when defining giant planets. 13 M is usually used as an upper mass limit for giant planets to separate from sub-stellar mass objects with significant central deuterium burning [e.g. 128]. However, as already mentioned in the manuscript, the exact deuterium burning mass-limit depends on several parameters of the object. Another possible limitation of mass-limit based definition is that it does not involve formation mechanism of the object. Finally, the mass distribution of sub-stellar objects show no any special feature around this mass limit [e.g. 185]. Schlaufman [186] proposed that the separation between giant planets and BDs should be somewhere between 4 M to 10 M , assuming that the sub-stellar objects formed via CA and GI can be statistically well separated. In particular, Schlaufman [186] assumed that the giant planets preference for metal-rich host stars is an indication of their formation through CA [but see 59]. Schneider et al. [13], in exoplanet.eu set an upper limit of 25 M for giant planets. Although appreciating the importance of having a formation-based definition for planets, the authors correctly noted about the difficulty (or at this point even probably impossibility) of knowledge about the formation mechanism for the planets. The choice of 25 M the authors justified by the start of a dip at around this mass in the observed distribution in M sin i of substellar mass objects [e.g. 185,[187][188][189].
Exoplanets are sometimes categorized into different populations defined by distinct metallicity distributions [e.g. 136]. Ribas and Miralda-Escudé [190] found that the average metallicity of stars with planets of mass < 4 M is lower than that of higher mass planet hosts by about 0.15 dex. However, in their comparison the authors did not separate the super-Earth and Neptune mass planets from giant planets. Moreover, these low mass planets in their sample mostly had super-solar metallicities, while the majority of later-on detected low-mass planets hosts have subsolar or solar metallicities [e.g. 159,191] if not orbiting their stars at very short periods [161,192]. Using a significantly larger sample of giant planets with masses greater than 0.1 M and with homogeneously derived host metallicities derived in SWEET-Cat, Adibekyan et al. [161] revisited the metallicity distribution of very massive giants. The authors found that the stars hosting massive Jupiters (mass > 4 M ) have an average metallicity of 0.083 ± 0.032 dex, while the stars hosting planets with masses between 1 and 4 M had an average metallicity of 0.149 ± 0.016 dex. Despite the lower average metallicity (lower by 0.066 ± 0.036 dex) of massive Jupiter hosts when compared to the relatively lower mass Jupiter (1 < M p < 4 M ) hosts, the performed KS statistics suggested a probability of about 15% that the two subsamples have the same underlying metallicity distribution. Adibekyan et al. [161] concluded that a simple separation in mass at 4 M does not reveal two different populations in metallicity.
Recently, Santos et al. [193] using the SWEET-Cat database also addressed the question whether very massive Jupiters tend to form at metallicities lower than their less massive counterparts. The authors found a convincing observational evidence that giant planets with masses above and below ∼4 M constitute two distinct populations 24 . They showed that the hosts of giant planets with masses < 4 M have metallicities on average higher than that of the fields stars without planets. At the same time, the hosts of their most massive planets were more massive stars with relatively low metallicities, similar to that observed in field stars of similar mass. Compensation of the lack of metallicity by higher stellar mass (thus high disk mass and availability of larger amount of planet building material) was not considered as a convincing explanation. Based on the aforementioned observational results Santos et al. [193] proposed that the planets in the two mass regimes may form in different ways: low-mass giant planets are formed in metal-rich disks through CA [e.g. 36] and the more massive planets are formed in massive and metal-poor disks via GI [e.g. 194,195].
It is interesting to see that the transiting planets with a radii larger than about 10 R ⊕ (note that the radius of Jupiter is 11.2 R ⊕ ) also seem to have an average metallicity slightly (perhaps statistically not significant) lower than the giant planets with smaller radii [see Fig.3, Fig. 7, and Fig. 6 of 138,139,196, respectively]. However, the number of transiting planets with very large radii and precise host metallicities is small and it is difficult to conclude whether there is a break-point radius above which planet hosts show low metallicities. Perhaps, besides the sample size and precision in metallicity, the main limitation of testing the aforementioned hypothesis is the very weak correlation between mass and radius for these H/He dominated giant planets [e.g. 35,197,198].
Very recently, Schlaufman [186] also studied the upper mass limit of exoplanets in connection with metallicity distributions of planets and sub-stellar mass objects. The author used a sample of 119 systems with planetary-mass companion having both RV and transit signals and true masses 25 > 0.1 M . These planets had orbits mostly within 10 days. His sample of subs-stellar objects consisted of 27 objects with masses < 300 M . The constrain of having RV and transit signal restricted the orbital periods of the subs-stellar mass objects to be mostly within 50 days. First, Schlaufman [186] applied a clustering algorithm to separate the giant planets and non-planets based on the metallicities of their host/primary stars. The objects with M < 4 M and M > 10 M the algorithm unanimously classified as belonging to two different groups: giant planets and BD/low-mass stars, respectively. The clustering analysis of Schlaufman [186] using the SWEET-Cat sample suggested a separation at a similar mass of 4 M (note the similarity in the separation mass proposed by Santos et al. [193]). In the second part of the analysis, "Moving Median Analysis", the author estimated the mass that separates the giant planets that are formed exclusively through CA and objects that do not show preference towards the high metallicity of their hosts and thus, cannot be formed through CA. Besides the assumption that the giant planets form only through CA, Schlaufman [186] also build the analysis on the theoretical prediction of Mordasini et al. [36] that "massive objects formed by core accretion should only occur around the most metal-rich primaries". He found that the point at which secondaries do not preferentially orbit metal-rich primary stars occurs at about 10 M .
Inclusion of stellar mass binary systems in the sample of Schlaufman [186] might decrease the average metallicity of the massive secondaries sample, which in turn, would impact on the estimation of the breakpoint mass separating the giant planets from BDs. Recent observations suggest that the close (period < 100 -10,000 days, depending on the work) binary fraction of solar type stars strongly anti-correlates with metallicity [199][200][201]. In particular, Moe et al. [201] showed that the binary fraction decreases from 40% to 10% in the mtallicity range of −1.0 dex and 0.5 dex. It is interesting to note that the wide orbit binary fraction of solar-type stars and close binary fraction of O-and B-type stars does not significantly depend on the metallicity across −1.5 < [Fe/H] < 0.5 dex range.
In the bottom panel of Fig. 9 I show the dependence of mass of giant planets (yellow) and sub-stellar objects (green) on the metallicity of the host star from the sample of Schlaufman [186]. In the same figure I also plot the giant planets (masses between 50 M ⊕ and 4 M ) and super-giant planets (masses greater than 4 M ) orbiting around stars with metallicities homogeneously derived in SWEET-Cat. I note that while the masses of the objects from the sample of Schlaufman [186] are true masses, the masses for the SWEET-Cat objects are mostly the minimum masses (M sin i).  [36]: at least not with the parameters (e.g. disk lifetime, disk mass) that were adopted in that specific model. It is interesting to see that almost all the sub-stellar objects with masses below about 60-70 M of the sample of Schlaufman [186] are observed at super-solar metallicities. At the same time, almost all the low-mass stars 26 with masses above about 150 M have sub-solar metallicities. This latter trend is expected as the stellar binary fraction decreases with metallicity [e.g . 201]. It is also interesting to see (compare the yellow and black/blue curves on the top panel of Fig. 9) that the giant planets of Schlaufman [186] are on average more metallic (0.119 ± 0.015 dex) when compared with the giant planets from SWEET-Cat (0.091 ± 0.009 dex). This difference is probably due to the different orbital period distributions of the two samples. As discussed in Sect. 3.2, the short period Jupiters tend to orbit around more metallic stars than their longer period counterparts [e.g. 161,162].   Table 1) than the hosts of more massive planets. The average metallicity of GPH is 0.091 ± 0.009 dex and the average metallicity of SGPH is 0.021 ± 0.022 dex. Despite the significant differences in the average metallicity distributions of giant planets and super-massive planets and/or BDs one probably should not make a bold conclusion about their different formation channels. These planets might have formed via the same mechanisms in different environments. As discussed in Santos et al. [193], the hosts of super-giant planets on average are more massive than the hosts of Jupiter-like planets which means that the super-giant planets have been formed in more massive disks. In order to study the impact of stellar mass on the obtained results I divide the samples into dwarf host stars (M < 1.5 M ) and giant host stars (M > 1.5 M ). The stellar masses were derived using the calibration presented in Torres et al. [205]. This calibration is based on stellar temperature, surface gravity, and metallicity of the stars. Fig. 10 and Table 1 show that when considering only dwarf stars the metallicity distribution of GPH and SGPH are almost identical. The average metallicity of SGP dwarf hosts is 0.106 ± 0.024 dex and the average metallicity of GP dwarf hosts is 0.100 ± 0.009 dex. Almost all the dwarf stars hosting super-giant planets have solar and super-solar metallicities (see Fig. 11). Such a metallicity distribution is expected for planets formed through CA [e.g. 36] rather than GI/TD [59,206]. Only four dwarf stars (HD111232, HD114762, HD181720, and HD22781) hosting SGP have metallicities bellow −0.3 dex. It is worth noting that most of the dwarf stars hosting planets with [Fe/H] < −0.3 dex are enhanced in α-elements, such as Mg and Si [207][208][209]. α-enhanced iron-poor stars are also enhanced in oxygen [e.g. 210] and carbon [e.g . 211]. Enhancement in such abundant elements as O, C, Mg, and Si makes these stars not very metal-poor, but iron-poor. Interestingly, the picture is drastically different when comparing the metallicity distributions of SGP and GP giant hosts. The average metallicity of SGP giant hosts is −0.122 ± 0.031 dex and the average metallicity of GP giant hosts is −0.005 ± 0.027 dex. The performed KS test suggest a probability of 0.06 that the metallicies of the two samples come from the same parent distribution. In general, the data shows that giant stars hosting planets (giant and super-giant) are less metallic than the dwarf hosts and that practically there are no giant stars hosting planets with super solar metallicities [Fe/H] > 0.2 dex (see Fig. 11).
To understand better the reason of dwarf and giant stars hosting planets having different metallicity distributions, in Fig. 11 I show the distribution of stars with and without planets in the stellar mass -metallicity diagram. It is important to note that the stars without planets have been searched for planets. The FGK field dwarf stars without detected planets have been observed within the HARPS GTO planet search program [103] and GK field giant stars without planets have been observed as part of the CORALIE planet search program [212] and the Okayama Planet search program [213]. The masses of these field stars were derived using the same calibration formulae [205] as for the exoplanet host stars. Fig. 11 clearly shows that all the giant stars with and without planets have a limiting maximum metallicity of about 0.2 dex and the fraction of very low metallicity (e.g. [Fe/H] < -0.4 dex) giant stars is much less when compared with the lower mass dwarf stars in the same metallicity region. Several studies have already observed this tendency of evolved, giant stars lacking the metal-rich and very metal-poor stars [e.g. [213][214][215][216][217]. Due to their shorter main sequence lifetimes, most of these giant and evolved (the stars with M > 1.5 M have surface gravities between 1.5 and 3.0 dex) stars are younger than their dwarf counterparts. The younger age together with the age-metallicity dispersion relation [e.g. 218,219] might explain the narrower metallicity distribution of the giant stars. These young stars are mostly local since they do not have time to migrate within the Galaxy [e.g. 220,221]. If the lifetime of the stars is long, then the radial migration in the disc would bring them from the inner, metal-rich Galaxy [e.g. 220,221]. In addition, it is important to note that most large programs to search for planets around giant stars make a sample selection based on a cut-off in the B -V colour. This B -V cut-off result in the lack of low-gravity and massive stars with high-metallicities [222]. Besides comparing the metallicitis of the hosts of massive and very massive planets, it is very informative to compare their metallicity distributions with the distribution of stars with not detected planets. The cumulative distribution of [Fe/H] for FGK field dwarf stars of Adibekyan et al. [103], and GK field giant stars of Alves et al. [212] and Takeda et al. [213] is shown in Fig. 10. The figure and the performed KS statistics (see Table 1) suggest that the dwarf hosts of both GP and SGP are significantly more metallic than the single field dwarf stars. This is a very important result suggesting that metallicity plays a very important role in the formation of massive planets independent of mass.
When comparing the metallicity distributions of GP and SGP hosting giant stars with the metallicity distribution of filed giant stars without planets, the KS statistics suggests probabilities of 0.007 (for GPH) and 0.105 (for SGPH) that the samples come from the same parent distribution. However, it is very important to note that while the GPH are more metallic than the field giant stars (the mean metallicity of the [103] + [212] sample is −0.094 ± 0.007 dex), the giant hosts of SGP are slightly less metallic. This result suggest that giant planet (at least up to 4 M ) formation around massive stars is more efficient if the disk is more metallic. Contrary, it seems that the most massive giant planets and/or BDs do not need metal-rich disks to be formed around massive stars.
Since our sample of planet hosting and single stars do not consist a single, well defined sample, it is not possible to determine the occurrence rate of very massive planets and test whether it depends on metallicity and stellar mass. Perhaps, a precise knowledge about a dependence (or absence of a dependence) of very massive planets and BD occurrence rate on stellar mass would help to better understand planet formation in massive disks. Different independent studies suggest that giant planet occurrence rate increases with the mass of the host star [e.g. [223][224][225] for up to masses of about 2 M and rapidly drops for masses larger than about 2.5 M [e.g. 223] 27 . This dependence is well explained by planet formation models based on CA [e.g. 36,178,229]. GI models also predict an increase of giant planet formation efficiency with stellar mass 28 [230][231][232]. However, the short lifetime of massive disks (for stars with intermediate masses) affecting the formation efficiency of giant planets in the CA models should not be a problem for the giant planets formed through GI. Unfortunately, the limited number of BD detections does not permit to firmly conclude whether BD occurrence rate is a function of stellar mass or not. The occurrence rate of relatively short-period (within a few years) BDs around solar-mass stars is estimated to be between 0.5% and 1% [e.g. 188,233,234]. Borgniet et al. [235] did not find any BD mass companion within about 1000 days around a sample of 225 AF-type main sequence stars observed with SOPHIE and/or HARPS spectrographs. These authors estimated the upper limit of the occurrence rate of BDs around these stars to be below 4%. An occurrence rate of 1.6% for BDs around intermediate-mass giant and main-sequence stars was estimated by Jones et al. [236], which is slightly higher than the rate around lower mass stars. A similar occurrence rate of about 2% for BDs orbiting white dwarfs (these are typically progeny of AF-type main sequence stars) was determined by Girven et al. [237] 29 , again giving a tentative evidence that BD formation might be more efficient around massive stars. Interestingly, the parent stars of the BD candidates detected in Jones et al. [236] are metal-rich. According to the CA model of Mordasini et al. [19] the formation of super-massive planets and BDs are possible in massive and metal-rich disks at large distances.
In summary, our data and results do not support the previous hints and claims [e.g. 138,186,193] about the existence of a breakpoint planetary mass at 4 M above and bellow which planet formation channels are different. iant planet formation (independent of their mass) around solar-like stars preferentially occurs in metal-rich disks. These planets thus can be result of CA process or be formed by TD model of Nayakshin [59]. In contrast, it seems that the most massive planets and BDs orbiting massive stars can be formed regardless of the disk metallicity. Such planets are not typical outcome of CA models and it is natural to suggest that they might have been formed via GI which efficiently produce very massive planets in massive and metal-poor disks [59,194]. Interestingly, several studies suggest that BDs with masses above and bellow ∼ 42 M might have formed by different processes [e.g. 189,[238][239][240]. In particular, low-mass BDs can be formed by disk instability and the high-mass BDs via cloud fragmentation as stars [238].

Low-mass planets and metallicity
Ever since the first giant exoplanet was discovered orbiting a Sun-like star, the search has been ongoing for small, rocky planets around other stars, evocative of Earth and other terrestrial planets in the Solar System. Several ongoing and upcoming missions (e.g. TESS: Ricker et al. [154]; CHEOPS Fortier et al. [241]; PLATO Rauer et al. [155]) and instruments (e.g. HARPS Mayor et al. [242]; HARPS-N Cosentino et al. [243]; SOPHIE Perruchot et al. [244]; CARMENES Quirrenbach et al. [245]; ESPRESSO Pepe et al. [246]; SPIRou Artigau et al. [247]) are developed to search and characterize low-mass exoplanets, which should ultimately help us to understand the formation and evolution of these planets. However, the detection of low-mass/small-sized planets is not an easy task not only because of the very small photoelectric and RV signals, but also because of the activity signals coming from the host stars that often have the same magnitude as the planetary signal and can strongly perturb the detection of these planets and/or mimic a planetary signal [e.g. [248][249][250][251][252][253][254]. These difficulties not only limits the number of detected low-mass planets, but also make very hard to construct a control sample of stars without low-mass planets for RV surveys and make practically impossible 30 for transit surveys [255]. Obviously, this makes the comparison of the properties of stars with and without low-mass/small-sized planets very difficult. Such a comparison is crucial to constrain the models of low-mass planet formation and for understanding environmental conditions required for their formation. In fact no strong correlation between low-mass planet frequency and metallicity is predicted by most of the CA based models [e.g. 36,50]. However, there are subtle differences in the relation between the frequency of these planets and metallicity predicted by different CA models. For example, recent N-body simulations of planet formation via pebble accretion by Matsumura et al. [256] suggest that the formation efficiency of low-mass planets subtly depends on the stellar metallicitiy. Their simulations show that at high metallicities large number of low-mass planets can be formed, but majority of them leave the systems because of dynamical instabilities. In the CA model of Hasegawa and Pudritz [50] formation of low-mass planets is almost independent and their frequency slightly decrease towards a minimum (at [Fe/H] of about −0.2 dex) and then increases again with increasing metallicity. In their model low-mass planets are considered as failed cores of giant planets and the observed minimum in their frequency is related with the rise in the frequency of giant planets at that metallicity. Unlike most of the GI models [e.g. 55,257], GI/TD model of Nayakshin [63] can also explain the formation of low-mass planets. The model of Nayakshin [63], as the CA models, also predicts high occurrence rate of these planets at low metallicities. However, in their model the frequency of low-mass planets decreases at super-solar metallicities.
Despite the aforementioned difficulties many groups tried to study the metallicity distribution of low-mass stars. Based only on a sample of seven short-period Neptune-mass (M sin i < 21 M ⊕ ) planets, Udry et al. [258] suggested that their formation efficiency may not depend on the host metallicity. This hint was later supported by the results of Sousa et al. [87] who used a sample of 17 low-mass planets (M sin i < 25 M ⊕ ) orbiting FGK and M-type dwarf stars. With the increasing number of RV detected low-mass planets, several independent studies confirmed that these planets (with masses bellow 30-40 M ⊕ ) do not show any preference towards metal-rich FGK [e.g. 90,123,159,184,259,260] and M-type stars [e.g. [261][262][263][264]. However, it seems that the maximum mass of super-Earth/Neptune-class planets (M sin i 40 M ⊕ ) shows a dependence on metallicity [265,266]. As for the massive planets, when studying the metallicity distributions of low-mass planet hosts, the chosen upper planetary mass limit varies from work to work, mostly depending on the sample size. The mass below which planet formation does not depend on metallicity (if there is such a limit) is still not well known. It is worth to note, that both observations [159] and CA based models [e.g. 36,267] suggest a strong minimum starting at about 30 M ⊕ in the mass distribution of low-mass planets.
The left and right panels of Fig. 12 show the metallicity distributions of FGK type stars hosting planets or planetary systems with only low-mass detected planets (< 30 M ⊕ ) and systems where at least one of the planets have a mass of less than 30 M ⊕ . The comparison of these distributions with that of FGK type single dwarfs from the HARPS sample [103] shows that the systems with only low-mass planets are less metallic (−0.101 ± 0.031 dex) than the systems of low-mass planets accompanied by more massive planets (−0.021 ± 0.026 dex), but are slightly more metallic than the single stars (−0.159±0.009 dex). KS tests show that the difference in the metallicity distributions between stars with only low-mass planets and single stars is not statistically significant, while the metal excess of systems including low-mass planets when compared to the single stars is statistically significant ( Table 1). As it is well known, and already discussed in this manuscript, giant planets tend to orbit stars with high metallicities. Therefore, exclusion of low-mass planet hosting systems that also contain giant planets will artificially reduce the number of metallic planetary systems with low-mass planets. However, the occurrence rate of giant planets is several times less than the occurrence rate of low-mass planets [e.g. 159]. Thus, when using a well defined volume limited sample of stars with and without planets, where the fraction of low-and high-mass planet hosts reflect their occurrence rates, exclusion of stars hosting high-mass planets probably should not significantly affect the study of low-mass planet -metallicity relation, at least for metallicities (e.g. [Fe/H] 0.1 dex) at which the giant planet frequency is not extremely high. When excluding super-metallic stars ([Fe/H] > 0.1 dex) and comparing the metallicity distributions of single stars and low-mass planet hosts accompanied by massive planets, the KS tests provides a p-value of 0.17 that the two samples come from the same parent distribution. This test shows that the statistically significant difference in metallicity obtained for the single stars and stars hosting simultaneously low-and high-mass planets was due to the presence of super-metallic stars with giant planets. In Fig. 12 one can see that in our SWEET-Cat sample, the number of stars hosting only low-mass planets is 49 and the number of stars hosting low-and high-mass planets is 75-49 = 26. This relatively large number of multiple planetary systems with low-mass and giant planets is due to selection effects towards giant planets and the fact that about 90% of cold giant planets accompanied by a low-mass planet [268]. Majority of the planets in the Earth-like to Neptune-like regime have been discovered via transit method (see Fig. 1). The large number of such planets detected by Kepler provides a good base for a statistical studies of these small worlds. However, an important limitation of Kepler planet studies is the faintness of their host stars which makes difficult to characterize them very accurately. Most of the early studies based on transiting planets reached to a conclusion that stars hosting exoplanets with radii smaller than about 4 R ⊕ show wide range of metallicities, indistinguishable from the distribution of field stars without known planets [e.g. 191,255,269]. Here I should remind again that a star without known planets does not necessarily mean that the star does not have planets.
Mass determination of the small-sized planets (R < 4 R ⊕ ) revealed the existence of two sub-populations: rocky and gaseous. A division between rocky planets and small-sized planets with gaseous envelopes was suggested to occur at about 1.6 R ⊕ [270,271]. Recently, Fulton et al. [137] found a bimodal distribution for radii of Kepler short period (< 100 d) small-sized planets with a clear lack of planets with radii between 1.5 and 2.0 R ⊕ . Later, it was shown that the position of the gap (sometimes called "Fulton gap") depends on orbital distance of the planets [272] and the mass/type of the hosts [272,273]. In fact, Jin and Mordasini [274] showed that the location of the gap should depend on the composition of these planets. They argued that the observed gap is best explained by their cores being rocky (rather than icy) and that most close-in low-mass planets may be formed within the snow line. The presence of the gap is mostly explained by photoevaporaton of the atmospheres of these low-mass planets [166,272,[275][276][277], but see also [278,279] for alternative explanations], which may depend on metallicity [280]. In addition, CA based in-situ planet formation model of Dawson et al. [116] suggest that the presence or absence of gaseous atmosphere depends on the solid surface density of the protoplanetary discs where they have been formed. Similarly, Owen and Murray-Clay [280] found that the core masses of low-mass planets are larger at high metallicities allowing to accreate larger H/He envelopes. If the photoevaporaton is the mechanism producing the "Fulton gap" then the aforementioned results indicates that the metallicity distribution of the rocky planet hosts should be contaminated by the metallicities of the hosts of the photoevaporated cores that before the photoevaporaton had larger sizes.
Using a sample of about 400 Kepler exoplanet candidate host stars with homogeneously derived metallicities Buchhave et al. [136] found that there are three regimes of exoplanets -terrestrial-like planets (R < 1.7 R ⊕ ), gas dwarf planets (1.7 R ⊕ < R < 3.9 R ⊕ ), and ice or gas giant planets (R > 3.9 R ⊕ ) -that show different mean metallicities. In particular, Buchhave et al. [136] showed that the terrestrial-like planets have a mean metallicity close to the solar value, while the hosts of gas dwarfs and ice/gas giants are more metallic than stars without known planets [255]. Schlaufman [281] using the sample of Buchhave et al. [136] did not find a convincing evidence of the boundary at 1.7 R ⊕ , and argued that the data are better fitted by a continuous relation between radius and metallicity.
Most of the studies based on a relatively large number of transiting exoplanets with spectroscopically derived metallicities indicate that small-sized rocky planets (smaller than about 1.7 R ⊕ ) show no preference towards high metallicities [e.g. 136,138,139,192], while the occurrence rate of larger transiting planets (1.7-3.9 R ⊕ ) show a correlation with metallicity [e.g. 136,139,192]. By comparing the metallicity distribution of long-period (> 25 days) small-sized (1-1.8 R ⊕ ) planet hosts with the metallicity distribution of stellar hosts of planets with radii between 1 and 6 R ⊕ (independent of the orbital period) Owen and Murray-Clay [280] found that the terrestrial-like planets are more common around low-metallicity stars and might have been formed after the gas disc was dispersed. However, to confirm or infirm whether rocky planet formation is more prevalent around lower metallicity stars one needs to study the dependence of occurrence rate of this planets on metallicity. Petigura et al. [e.g. 139] showed that the occurrence rate of warm super-Earths (P = 10-100 days and R = 1.0-1.7 R ⊕ ) almost independent of metallicity, showing slight (insignificant) decrease with metallicity. However, it is important to note that for transiting planets it is very hard to estimate the fraction of stars with planets [282] and Petigura et al. [e.g. 139], as most of the works, used the average number of planets per star when studying the planet occurrence -metalliciy dependence. Obviously, these two estimates of the occurrence rates are different (they would be the same if there was no star with more than one planet), and show different dependencies on stellar metallicity [283]. Contrary to aforementioned results, Wang and Fischer [284] suggested a universal planet-metallicity correlation for planet of all sizes. It is important to note that, while the stellar metallicites of their planet candidate hosts were derived spectroscopically, the metallicities of their control sample stars were based on photometric stellar parameters with further conversion to more "representative" stellar parameters [284]. Zhu et al. [285] suggested that the high occurrence rate and low detection efficiency of low-mass planets can explain the discrepant results regarding the small-sized planet -metallicity correlations obtained by different authors.
When studying the metallicity dependence of low-mass/small size planets it is very important to take into account the relation between orbital periods of planets and their host stars metallicities [116,139,161,192,[286][287][288][289]. In particular, Adibekyan et al. [161] found that the super-Earth-like planets (M sin i < 10 M ⊕ ) orbiting metal-rich stars have orbital periods shorter than about 20 days, whereas planets orbiting metal-poor stars span a wide range of orbital periods [see also 116,287]. However, these results were contested by Mulders et al. [192] who argued that the observed trends might be due to selection effects. A systematic excess of short period ( 10 days) rocky planets (< 1.7 R ⊕ ) around metal-rich stars was reported in several works [139,192,288]. The metallicity preference of hot rocky planets is explained by the possible dependence of efficient inward migration of solids on metallicity [139], higher rates of planet-planet scattering in metal-rich disks [139], and/or dependence of planet trap at the inner age on metallicity [192,288]. An alternative possibility is that hot rocky planets and hot Jupiters might share the same formation mechanism [192].
Despite aforementioned results suggesting that low-mass/small planet formation can be efficiently occur around stars with wide range of metallicities there seems to be an observational evidence that at a given metallicity these planets prefer to orbit stars with high [Mg/Si] abundance ratio [290,291] and at low-iron regime preferably orbit stars enhanced in α-elements such as Mg and Si [208]. These results suggest that low-mass planet formation is not completely independent of the composition of the protoplanetary disk.

Summary and conclusion
In this manuscript I reviewed and discussed the dependence of planet occurrence on stellar metallicity (iron abundance). In the upcoming manuscripts I will review i) the dependencies observed between exoplanet properties (e.g. orbital properties, multiplicity) and the metallicity (iron abundance) of their host stars (Paper II: 'Heavy metal rules II. Exoplanet properties and stellar metallicity') and ii) the role of individual heavy and light elements in the formation of exoplanets and the link between composition of exoplanets and their hosts stars (Paper III: 'Heavy metal rules III. Exoplanets and elements other than iron').

Questions to answer or to think about
Before starting to write this review, I was hoping that by gathering/collecting all the available information and evidences coming from observations and theoretical predictions about formation of exoplanets of different types, I and/or the reader will get closer in answering to some important questions that I had in my mind. I hope the reader will be more successful in answering to this questions than the writer.
The origin of giant planet -metallicity correlation. As discussed in the manuscript, formation of giant planets depends on several characteristics of the protoplanetary disks that are inter-correlated and all can produce a correlation between giant planet occurrence and metallicity. The most direct explanation of the correlation would go as, large amount of metals in the disk means large amount of material to build cores of giants which would translate into higher probability to form gas giants. However, at high metallicities the disk lifetime is longer [114] which makes easier (enough time) to build cores of giants and hence means higher probability to form gas giants as well. Besides the disk lifetime and metallicity, the mass of the disk also plays an important role for the formation of giant planets [e.g. 178]. As discussed in the manuscript, the disk mass, in turn correlates with the metallicity (due to selection effects) and disk lifetime [e.g . 179]. The frequency of massive planets is lower than the frequency of systems with super-Earth-like planets. The latter systems, in principle, should have had enough material to build a core for giant planets. This may suggest that not the availability of core building blocks in the disk, but the conditions for planetary embryos to merge and acquire the critical core mass for giant planet formation is what might determine the formation probability of gas giants. Liu et al. [115] suggested that the parameter determining these conditions is the critical disk accretion rate, which correlates with metallicity and disk mass. Finally, as a very alternative explanation, it was proposed that the giant planet -metallicity correlation can have a secondary origin Haywood [292]. Because of the radial metallicity gradient in the Galaxy [e.g. 293,294], the giant planet -metallicity correlation can be easily recovered if the planet formation is more efficient in the inner Galaxy than in the solar neighborhood. Haywood [292] suggested that giant planet formation might be related to the surface density of molecular hydrogen H 2 in the Galaxy.
Low-mass/small-sized planets and metallicity. Until the recent advancement of the TD model [e.g. 59] terrestrial planet formation was explained exclusively by CA. Most of the CA models and the TD model of Nayakshin [63] 31 suggest practically no dependence of low-mass planet formation on disk metallicity, although still small differences exist in the predicted exact form of this dependence. The current observations seem to support the absence of correlation between the lowest-mass/smallest-sized planets and metallicity, at least for orbital periods longer than about 10 days [e.g. 139,192]. However, because of the very high frequency and low detectability of these planets it is still very difficult to firmly conclude whether terrestrial planet formation depends on metallicity or not [e.g. 285]. Very large samples of stars with and without planets will help to tackle this issue and help to select the model that explains the observations best.Hopefully such samples can be constructed in the near future with the ongoing and upcoming surveys and missions.

The main new results and conclusions
Besides making questions and extensive literature review about the role of metallicity on the formation of planets of different types, in this manuscript I also revisited most of the correlations related to metallicity of planet host stars reported in the literature. I used the SWEET-Cat database [82] that provides homogeneously derived stellar parameters (including metallicities) for a very large sample of planet host stars, especially planets detected through RV method. Among the many findings discussed in the manuscript, I highlight two main new results and findings.
• The study of the dependence of planet masses on the mean metallicity of their hosts revealed (see Sect. 3.2) that the hosts of sub-Jupiter mass planets (∼0.6 -0.9 M ) are systematically less metallic than the hosts of the Jupiter-like planets. The results seem to suggest that at high metallicities, the longer disk lifetime [114] and higher amount of planet building material [e.g. 36], allows a formation of more massive Jupiter-like planets than at lower metallicities. • Recently, several authors suggested that planets more massive than about 4 M tend to orbit around low-metallicity stars and might have formed in a different way than Jupiter-mass planets with < 4 M [e.g. 186,193]. The results of this study show that giant planets with masses above and bellow 4 M orbiting solar-like stars are preferentially metal-rich, which does not support the previous hints about the different formation mechanisms of these two groups of planets. Formation of these planets and the observed metallicity trend can be explained by the CA models [e.g. 36] more easily than by the GI and TD models [e.g. 55,59]. The results also show statistically significant difference in the metallicity distributions of giant stars (> 1.5 M ) hosting planets with masses greater or less than 4 M . Perhaps, GI based models should be able to explain the formation of the most massive giant planets (> 4 M ) in low-meallicity environment more easily than CA based models. It is thus suggested that planets of the same mass can be formed through different channels depending on the disk mass i.e. environmental conditions.
Summarizing all the results and trends discussed in the manuscript one can conclude that there is no yet a single CA nor GI/TD model that can effectively explain the formation (and the observed trends with metallicity) of planets of all types. CA models can qualitatively explain the formation and observed trends with metallicity for all planets, but the super-massive planets (> 4 M ) orbiting metal-poor giant stars. GI can be responsible for the formation of these later planets, but will not explain the formation of low-mass planets and the giant planet -metallicity correlation. The TD model of Nayakshin [59] can qualitatively explain most of the metallicity trends discussed in this manuscript, but will fail in explaining the super-giant planet -metallicity correlation observed for solar-type stars.
The planet formation models constrained by the observations undergo significant improvements on how different physical processes are implemented. Observational exoplanetology does not stand still as well. Several ongoing and up-coming ground-based surveys and space-based missions (sometimes driven by theoretical predictions) are planned that will provide the exoplanet community with unprecedented amount of high quality data. Everything speaks about the bright future of exoplanetology. My very general conclusion, thus is that heavy metal rules! I also acknowledge the support from FCT through Investigador FCT contract nr. IF/00650/2015/CP1273/CT0001. Acknowledgments: I would like to thank the EXO-Earth team members for working together on the subject of this manuscript for many years. It is my special pleasure to thank Pedro Figueira, João Faria, Nuno Santos, Elisa Delgado Mean, Sérgio Sousa, Olivier Demangeon, Susana Barros, and Vardan Elbakyan for very interesting conversations and discussions. My double thanks to Mahmoud Oshagh for the very interesting discussions and for inviting me to review this topic. I also thank Sergei Nayakshin for the provided data and for very interesting discussion about the formation of planets through GI. I would also like to thank David Adibekyan for his critical comments that were sometimes difficult to figure out and to implement. Finally, I would like to thank the referees for very constructive comments that helped to improve the quality of the manuscript.

Conflicts of Interest:
The author declares no conflict of interest.

Abbreviations
The following abbreviations are used in this manuscript: