In about 6 Giga years our Sun will evolve into a red giant and finally end its life as a white dwarf. This stellar metamorphosis will occur to virtually all known host stars of exo-planetary systems and is therefore crucial for their final fate. It is clear that the innermost planets will be engulfed and evaporated during the giant phase and that planets located farther out will survive. However, the destiny of planets in-between, at ~ 1 - 10 au, has not yet been investigated with a multi-planet tidal treatment. We here combine for the first time multi-planet interactions, stellar evolution, and tidal effects in an N-body code to study the evolution of a Neptune-Jupiter planetary system. We report that the fate of the Neptune-mass planet, located closer to the star than the Jupiter-mass planet, can be very different from the fate of a single Neptune. The simultaneous effects of gravitational interactions, mass loss and tides can drive the planetary system towards mean motion resonances. Crossing these resonances affects particularly the eccentricity of the Neptune and thereby also its fate, which can be engulfment, collision with the Jupiter-mass planet, ejection from the system, or survival at a larger separation.

The diversity in mass and composition of planetary atmospheres stems from the different building blocks present in protoplanetary discs and from the different physical and chemical processes that these experience during the planetary assembly and evolution. This review aims to summarise, in a nutshell, the key concepts and processes operating during planet formation, with a focus on the delivery of volatiles to the inner regions of the planetary system.

Recent observations of protoplanetary disks have revealed ring-like structures that can be associated to pressure maxima. Pressure maxima are known to be dust collectors and planet migration traps. Most of planet formation works are based either on the pebble accretion model or on the planetesimal accretion model. However, recent studies proposed the possible formation of Jupiter by the hybrid accretion of pebbles and planetesimals. We aim to study the full process of planet formation consisting of dust evolution, planetesimal formation and planet growth at a pressure maximum in a protoplanetary disk. We compute, through numerical simulations, the gas and dust evolution, including dust growth, fragmentation, radial drift and particle accumulation at a pressure bump. We also consider the formation of planetesimals by streaming instability and the formation of a moon-size embryo that grows into a giant planet by the hybrid accretion of pebbles and planetesimals. We find that pressure maxima in protoplanetary disks are efficient collectors of dust drifting inwards. The condition of planetesimal formation by streaming instability is fulfilled due to the large amount of dust accumulated at the pressure bump. Then, a massive core is quickly formed (in $\sim 10^4$ yr) by the accretion of pebbles. After the pebble isolation mass is reached, the growth of the core slowly continues by the accretion of planetesimals. The energy released by planetesimal accretion delays the onset of runaway gas accretion, allowing a gas giant to form after $\sim$1 Myr of disk evolution. The pressure maximum also acts as a migration trap. Pressure maxima in protoplanetary disks are preferential locations for dust traps, planetesimal formation by streaming instability and planet migration traps. All these conditions allow the fast formation of a giant planet by the hybrid accretion of pebbles and planetesimals.

As planets grow the exchange of angular momentum with the gaseous component of the protoplanetary disc produces a net torque resulting in a variation of the semi-major axis of the planet. For low-mass planets not able to open a gap in the gaseous disc this regime is known as type I migration. Pioneer works studied this mechanism in isothermal discs finding fast inward type I migration rates that were unable to reproduce the observed properties of extrasolar planets. In the last years, several improvements have been made in order to extend the study of type I migration rates to non-isothermal discs. Moreover, it was recently shown that if the planet's luminosity due to solid accretion is taken into account, inward migration could be slowed down and even reversed. In this work, we study the planet formation process incorporating, and comparing, updated type I migration rates for non-isothermal discs and the role of planet's luminosity over such rates. We find that the latter can have important effects on planetary evolution, producing a significant outward migration for the growing planets.

Most of planet formation models that incorporate planetesimal fragmentation consider a catastrophic impact energy threshold for basalts at a constant velocity of 3 km/s during all the process of the formation of the planets. However, as planets grow the relative velocities of the surrounding planetesimals increase from velocities of the order of m/s to a few km/s. In addition, beyond the ice line where giant planets are formed, planetesimals are expected to be composed roughly by 50 percentage of ices. We aim to study the role of planetesimal fragmentation on giant planet formation considering planetesimal catastrophic impact energy threshold as a function of the planetesimal relative velocities and compositions. We improve our model of planetesimal fragmentation incorporating a functional form of the catastrophic impact energy threshold with the planetesimal relative velocities and compositions. We also improve in our model the accretion of small fragments produced by the fragmentation of planetesimals during the collisional cascade considering specific pebble accretion rates. We find that a more accurate and realistic model for the calculation of the catastrophic impact energy threshold tends to slow down the formation of massive cores. Only for reduced grain opacity values at the envelope of the planet, the cross-over mass is achieved before the disk time-scale dissipation. While planetesimal fragmentation favors the quick formation of massive cores of 5-10 Earth masses the cross-over mass could be inhibited by planetesimal fragmentation. However, grain opacity reduction or pollution by the accreted planetesimals together with planetesimal fragmentation could explain the formation of giant planets with low-mass cores.

We study local and non-local Polyakov Nambu-Jona-Lasinio models and analyze their respective phase transition diagram. We construct hybrid stars using the zero temperature limit of the local and non-local versions of Nambu-Jona-Lasinio model for quark matter and the modern GM1(L) parametrization of the non-linear relativistic mean field model for hadronic matter. We compare our models with data from PSR J1614-2230 and PSR J0343+0432 and also from GW170817 and its electromagnetic counterpart GRB170817A and AT2017gfo. We study observational signatures of the appearance of a mixed phase as a result of modeling a phase transition that mimics the Gibbs formalism and compare the results with the sharp first-order phase transition obtained using the Maxwell construction. We also study in detail the $g$-mode associated with discontinuities in the equation of state, and calculate non-radial oscillation modes using relativistic Cowling approximation.

We present preliminary results of the study of intense magnetic fields effects on hybrid stars. For the description of the hadronic phase, we use the relativistic mean field approximation. For the quark matter phase, we employ the Field Correlator Method formalism. Once the the equation of state is built, we discuss the pressure anisotropy due the presence of the strong magnetic field. Finally, we calculate the structure of the compact stars using magnetized hybrid equations of state and their oscillation modes related with the emission of gravitational waves.

The first direct detection of gravitational waves has opened a new window to study the Universe and would probably start a new era: the gravitational wave Astronomy. Gravitational waves emitted by compact objects like neutron stars could provide significant information about their structure, composition and evolution. In this paper we calculate, using the relativistic Cowling approximation, the oscillations of compact stars focusing on hybrid stars, with and without a mixed phase in their cores. We study the existence of a possible hadron-quark phase transition in the central regions of neutron stars and the changes it produces on the gravitational modes frequencies emitted by these stars. We pay particular attention to the $g$-modes, which are extremely important as they could signal the existence of pure quark matter inside neutron stars. Our results show a relationship between the frequency of the $g$-modes and the constant speed of sound parametrization for the quark matter phase. We also show that the inclusion of color superconductivity produces an increase on the oscillation frequencies. We propose that observations of $g$-modes with frequencies $f_{\rm g}$ between $1$ kHz and $1.5$ kHz should be interpreted as an evidence of a sharp hadron-quark phase transition in the core of a compact object.

Within the framework of the Cowling approximation, we present the frequencies of the f (fundamental), p 1 (first pressure) and g (gravitational) of quadrupolar oscillation modes of compact objects constructed using different equations of state. Special attention will be put in hybrid stars, formed by a quark core and a hadronic envelope. In this context, g modes are only present if a sharp first order phase transition occurs between hadronic and quark phases.

Aims. We study the formation and water delivery of planets in the habitable zone (HZ) around solar-type stars. In particular, we study different dynamical environments that are defined by the most massive body in the system. Methods. First of all, a semi-analytical model was used to define the mass of the protoplanetary disks that produce each of the five dynamical scenarios of our research. Then, we made use of the same semi-analytical model to describe the evolution of embryos and planetesimals during the gaseous phase. Finally, we carried out N-body simulations of planetary accretion in order to analyze the formation and water delivery of planets in the HZ in the different dynamical environments. Results. Water worlds are efficiently formed in the HZ in different dynamical scenarios. In systems with a giant planet analog to Jupiter or Saturn around the snow line, super-Earths tend to migrate into the HZ from outside the snow line as a result of interactions with other embryos and accrete water only during the gaseous phase. In systems without giant planets, Earths and super-Earths with high water by mass contents can either be formed in situ in the HZ or migrate into it from outer regions, and water can be accreted during the gaseous phase and in collisions with water-rich embryos and planetesimals. Conclusions. The formation of planets in the HZ with very high water by mass contents seems to be a common process around Sun- like stars. Our research suggests that such planets are still very efficiently produced in different dynamical environments. Moreover, our study indicates that the formation of planets in the HZ with masses and water contents similar to those of Earth seems to be a rare process around solar-type stars in the systems under consideration.

Several observational works have shown the existence of Jupiter-mass planets covering a wide range of semi-major axes around Sun-like stars. We aim to analyse the planetary formation processes around Sun-like stars that host a Jupiter-mass planet at intermediate distances ranging from $\sim$1 au to 2 au. Our study focusses on the formation and evolution of terrestrial-like planets and water delivery in the habitable zone (HZ) of the system. Our goal is also to analyse the long-term dynamical stability of the resulting systems. A semi-analytic model was used to define the properties of a protoplanetary disk that produces a Jupiter-mass planet around the snow line, which is located at $\sim$2.7 au for a solar-mass star. Then, it was used to describe the evolution of embryos and planetesimals during the gaseous phase up to the formation of the Jupiter-mass planet, and we used the results as the initial conditions to carry out N-body simulations of planetary accretion. Our simulations produce three different classes of planets in the HZ: 'water worlds', with masses between 2.75 $M_{\oplus}$ and 3.57 $M_{\oplus}$ and water contents of 58% and 75% by mass, terrestrial-like planets, with masses ranging from 0.58 $M_{\oplus}$ to 3.8 $M_{\oplus}$ and water contents less than 1.2% by mass, and 'dry worlds', simulations of which show no water. A relevant result suggests the efficient coexistence in the HZ of a Jupiter-mass planet and a terrestrial-like planet with a percentage of water by mass comparable to the Earth. Moreover, our study indicates that these planetary systems are dynamically stable for at least 1 Gyr. Systems with a Jupiter-mass planet located at 1.5 au - 2 au around solar-type stars are of astrobiological interest. These systems are likely to harbour terrestrial-like planets in the HZ with a wide diversity of water contents.

Population synthesis models of planetary systems developed during the last $\sim$15 years could reproduce several of the observables of the exoplanet population, and also allowed to constrain planetary formation models. We present our planet formation model, which calculates the evolution of a planetary system during the gaseous phase. The code incorporates relevant physical phenomena for the formation of a planetary system, like photoevaporation, planet migration, gas accretion, water delivery in embryos and planetesimals, a detailed study of the orbital evolution of the planetesimal population, and the treatment of the fusion between embryos, considering their atmospheres. The main goal of this work, unlike other works of planetary population synthesis, is to find suitable scenarios and physical parameters of the disc to form solar system analogs. We are specially interested in the final planet distributions, and in the final surface density, eccentricity and inclination profiles for the planetesimal population. These final distributions will be used as initial conditions for N-body simulations, to study the post-oligarchic formation in a second work. We then consider different formation scenarios, with different planetesimal sizes and different type I migration rates. We find that solar system analogs are favored in massive discs, with low type I migration rates, and small planetesimal sizes. Besides, those rocky planets within their habitables zones are dry when discs dissipate. At last, the final configurations of solar system analogs include information about the mass and semimajor-axis of the planets, water contents, and the properties of the planetesimal remnants.

We analyze the process of planetary scattering around M0-type stars. To do this, we carry out N-body simulations with three Jupiter-mass planets close to their instability limit together with an outer planetesimal disk. This paper focuses on the analysis of systems in which a single Jupiter-mass planet survives after the dynamical instability event. The small body reservoirs show different dynamical behaviors. In fact, our simulations produce particles on prograde and retrograde orbits, as well as particles whose orbital plane flips from prograde to retrograde and back again along their evolution. Such particles are called "Type-F particles". We find strong correlations between the inclination $i$ and the ascending node longitude $\Omega$ of such particles. First, $\Omega$ librates around 90$^{\circ}$ or/and 270$^{\circ}$. This property is very important since it represents a necessary and sufficient condition for the flipping of an orbit. Moreover, the libration periods of $i$ and $\Omega$ are equal and they are out to phase by a quarter period. We also remark that the larger the libration amplitude of $i$, the larger the libration amplitude of $\Omega$. Finally, we analyze the initial conditions of Type-F particles of all our simulations immediately after the dynamical instability event, when a single Jupiter-mass planet survives in the system. We carry out this study with the goal to determine the parameter space that lead to the flipping of an orbit. Our results suggest that the orbit of a test particle can flip for any value of its initial eccentricity, although we found only two Type-F particles with initial inclinations $i <$ 17$^{\circ}$. Moreover, our study indicates that the minimum value of the inclination of the Type-F particles in a given system decreases with an increase in the eccentricity of the giant planet.

The discovery of giant planets in wide orbits represents a major challenge for planet formation theory. In the standard core accretion paradigm planets are expected to form at radial distances $\lesssim 20$ au in order to form massive cores (with masses $\gtrsim 10~\textrm{M}_{\oplus}$) able to trigger the gaseous runaway growth before the dissipation of the disc. This has encouraged authors to find modifications of the standard scenario as well as alternative theories like the formation of planets by gravitational instabilities in the disc to explain the existence of giant planets in wide orbits. However, there is not yet consensus on how these systems are formed. In this letter, we present a new natural mechanism for the formation of giant planets in wide orbits within the core accretion paradigm. If photoevaporation is considered, after a few Myr of viscous evolution a gap in the gaseous disc is opened. We found that, under particular circumstances planet migration becomes synchronised with the evolution of the gap, which results in an efficient outward planet migration. This mechanism is found to allow the formation of giant planets with masses $M_p\lesssim 1 M_{\rm Jup}$ in wide stable orbits as large as $\sim$130 au from the central star.

In the classical core-accretion planet formation scenario, rapid inward migration and accretion timescales of kilometer size planetesimals may not favor the formation of massive cores of giant planets before the dissipation of protoplanetary disks. On the other hand, the existence of pressure maxima in the disk could act as migration traps and locations for solid material accumulation, favoring the formation of massive cores. We aim to study the radial drift of pebbles and planetesimals and planet migration at pressure maxima in a protoplanetary disk and their implications for the formation of massive cores as triggering a gaseous runaway accretion phase. The time evolution of a viscosity driven accretion disk is solved numerically introducing a a dead zone as a low-viscosity region in the protoplanetary disk. A population of pebbles and planetesimals evolving by radial drift and accretion by the planets is also considered. Finally, the embryos embedded in the disk grow by the simultaneous accretion of pebbles, planetesimals and the surrounding gas. Our simulations show that the pressure maxima generated at the edges of the low-viscosity region of the disk act as planet migration traps, and that the pebble and planetesimal surface densities are significantly increased due to the radial drift towards pressure maxima locations. However, our simulations also show that migration trap locations and solid material accumulation locations are not exactly at the same positions. Thus, a planet's semi-major axis oscillations around zero torque locations, predicted by MHD and HD simulations, are needed for the planet to accrete all the available material accumulated at the pressure maxima. Pressure maxima generated at the edges of a low-viscosity region of a protoplanetary disk seem to be preferential locations for the formation and trap of massive cores.

In the standard model of core accretion, the formation of giant planets occurs by two main processes: first, a massive core is formed by the accretion of solid material; then, when this core exceeds a critical value (typically greater than 10 Earth masses) a gaseous runaway growth is triggered and the planet accretes big quantities of gas in a short period of time until the planet achieves its final mass. Thus, the formation of a massive core has to occur when the nebular gas is still available in the disk. This phenomenon imposes a strong time-scale constraint in giant planet formation due to the fact that the lifetimes of the observed protoplanetary disks are in general lower than 10 Myr. The formation of massive cores before 10 Myr by accretion of big planetesimals (with radii > 10 km) in the oligarchic growth regime is only possible in massive disks. However, planetesimal accretion rates significantly increase for small bodies, especially for pebbles, particles of sizes between mm and cm, which are strongly coupled with the gas. In this work, we study the formation of giant planets incorporating pebble accretion rates in our global model of planet formation.

To study the terrestrial-type planet formation during the post oligarchic growth, the initial distributions of planetary embryos and planetesimals used in N-body simulations play an important role. Most of these studies typically use ad hoc initial distributions based on theoretical and numerical studies. We analyze the formation of planetary systems without gas giants around solar-type stars focusing on the sensitivity of the results to the particular initial distributions of planetesimals and embryos. The formation of terrestrial planets in the habitable zone (HZ) and their final water contents are topics of interest. We developed two different sets of N-body simulations from the same protoplanetary disk. The first set assumes ad hoc initial distributions for embryos and planetesimals and the second set obtains these distributions from the results of a semi-analytical model which simulates the evolution of the gaseous phase of the disk. Both sets form planets in the HZ. Ad hoc initial conditions form planets in the HZ with masses from $0.66M_{\oplus}$ to $2.27M_{\oplus}$. More realistic initial conditions obtained from a semi-analytical model, form planets with masses between $1.18M_{\oplus}$ and $2.21M_{\oplus}$. Both sets form planets in the HZ with water contents between 4.5% and 39.48% by mass. Those planets with the highest water contents respect to those with the lowest, present differences regarding the sources of water supply. We suggest that the number of planets in the HZ is not sensitive to the particular initial distribution of embryos and planetesimals and thus, the results are globally similar between both sets. However, the main differences are associated to the accretion history of the planets in the HZ. These discrepancies have a direct impact in the accretion of water-rich material and in the physical characteristics of the resulting planets.

Models of planet formation are mainly focused on the accretion and dynamical processes of the planets, neglecting their chemical composition. In this work, we calculate the condensation sequence of the different chemical elements for a low-mass protoplanetary disk around a solar-type star. We incorporate this sequence of chemical elements (refractory and volatile elements) in our semi-analytical model of planet formation which calculates the formation of a planetary system during its gaseous phase. The results of the semi-analytical model (final distributions of embryos and planetesimals) are used as initial conditions to develope N-body simulations that compute the post-oligarchic formation of terrestrial-type planets. The results of our simulations show that the chemical composition of the planets that remain in the habitable zone has similar characteristics to the chemical composition of the Earth. However, exist differences that can be associated to the dynamical environment in which they were formed.

In general, most of the studies of terrestrial-type planet formation typically use ad hoc initial conditions. In this work we improved the initial conditions described in Ronco & de El\'ia (2014) starting with a semi-analytical model wich simulates the evolution of the protoplanetary disk during the gas phase. The results of the semi-analytical model are then used as initial conditions for the N-body simulations. We show that the planetary systems considered are not sensitive to the particular initial distribution of embryos and planetesimals and thus, the results are globally similar to those found in the previous work.

In the standard model of core accretion, the cores of the giant planets form by the accretion of planetesimals. In this scenario, the evolution of the planetesimal population plays an important role in the formation of massive cores. Recently, we studied the role of planetesimal fragmentation in the in situ formation of a giant planet. However, the exchange of angular momentum between the planet and the gaseous disk causes the migration of the planet in the disk. In this new work, we incorporate the migration of the planet and globally study the role of planet migration in the formation of a massive core when the population of planetesimals evolves by planet accretion, migration due to the nebular drag, and fragmentation due to planetesimal collisions.

In the standard scenario of planet formation, terrestrial planets and the cores of the giant planets are formed by accretion of planetesimals. As planetary embryos grow the planetesimal velocity dispersion increases due to gravitational excitations produced by embryos. The increase of planetesimal relative velocities causes the fragmentation of them due to mutual collisions. We study the role of planetesimal fragmentation on giant planet formation. We analyze how planetesimal fragmentation modifies the growth of giant planet's cores for a wide range of planetesimal sizes and disk masses. We incorporate a model of planetesimal fragmentation into our model of in situ giant planet formation. We calculate the evolution of the solid surface density (planetesimals plus fragments) due to the accretion by the planet, migration and fragmentation. The incorporation of planetesimal fragmentation significantly modifies the process of planetary formation. If most of the mass loss in planetesimal collisions is distributed in the smaller fragments, planetesimal fragmentation inhibits the growth of the embryo for initial planetesimals of radii lower than 10 km. Only for initial planetesimals of 100 km of radius, and disks greater than 0.06 Msun, embryos achieve masses greater than the mass of the Earth. However, even for such big planetesimals and massive disks, planetesimal fragmentation induces the quickly formation of massive cores only if most of the mass loss in planetesimal collisions is distributed in the bigger fragments. Planetesimal fragmentation seems to play an important role in giant planet formation. The way in which the mass loss in planetesimal collisions is distributed leads to different results, inhibiting or favoring the formation of massive cores.

We develop a semi-analytical model for computing planetary system formation with the aim of explaining the observed diversity of planetary systems architectures and relate this primordial diversity with the initial properties of the disc where they were born. We adopt different initial conditions based on recent results in protoplanetary discs observations, to generate a variety of planetary systems and analyze them statistically. We explore the relevance of the mass and size of the disc, its metallicity, the mass of the central star and the time-scale of gaseous disc dissipation, in defining the architecture of the planetary system. We also test different values of some key parameters of our model, to find out which factors best reproduce the diverse sample of observed planetary systems. According to this, we predict which systems are the most common in the solar neighbourhood. Our results show that planetary systems with only terrestrial planets are the most common, being the only planetary systems formed when considering low metallicity discs and which also represent the best environment for the developing of rocky, potentially habitable planets. We also found that planetary systems like our own are not rare in the solar neighbourhood, being its formation favoured in massive discs where there is not a large accumulation of solids in the inner region of the disc. Regarding the planetary systems that harbor hot and warm Jupiter planets, we found that this systems are born in very massive, metal-rich discs. Also a fast migration rate is required in order to form these systems. According to our results, most of the hot and warm Jupiter systems are composed by only one giant planet, which is also a tendency of the current observational data.

In the last few years, the so-called "Nice model" has got a significant importance in the study of the formation and evolution of the solar system. According to this model, the initial orbital configuration of the giant planets was much more compact than the one we observe today. We study the formation of the giant planets in connection with some parameters that describe the protoplanetary disk. The aim of this study is to establish the conditions that favor their simultaneous formation in line with the initial configuration proposed by the Nice model. We focus in the conditions that lead to the simultaneous formation of two massive cores, corresponding to Jupiter and Saturn, able to achieve the cross-over mass (where the mass of the envelope of the giant planet equals the mass of the core, and gaseous runway starts) while Uranus and Neptune have to be able to grow to their current masses. We compute the in situ planetary formation, employing the numerical code introduced in our previous work, for different density profiles of the protoplanetary disk. Planetesimal migration is taken into account and planetesimals are considered to follow a size distribution between $r_p^{min}$ (free parameter) and $r_p^{max}= 100$ km. The core's growth is computed according to the oligarchic growth regime. The simultaneous formation of the giant planets was successfully completed for several initial conditions of the disk. We find that for protoplanetary disks characterized by a power law ($\Sigma \propto r^{-p}$), smooth surface density profiles ($p \leq 1.5$) favor the simultaneous formation. However, for steep slopes ($p\sim 2$, as previously proposed by other authors) the simultaneous formation of the solar system giant planets is unlikely ...

In order to explain the main characteristics of the observed population of extrasolar planets and the giant planets in the Solar System, we need to get a clear understanding of which are the initial conditions that allowed their formation. To this end we develop a semi-analytical model for computing planetary systems formation based on the core instability model for the gas accretion of the embryos and the oligarchic growth regime for the accretion of the solid cores. With this model we explore not only different initial discs profiles motivated by similarity solutions for viscous accretion discs, but we also consider different initial conditions to generate a variety of planetary systems assuming a large range of discs masses and sizes according to the last results in protoplanetary discs observations. We form a large population of planetary systems in order to explore the effects in the formation of assuming different discs and also the effects of type I and II regimes of planetary migration, which were found to play fundamental role in reproducing the distribution of observed exoplanets. Our results show that the observed population of exoplanets and the giant planets in the Solar System are well represented when considering a surface density profile with a power law in the inner part characterized by an exponent of -1, which represents a softer profile when compared with the case most similar to the MMSN model case.

The core accretion mechanism is presently the most widely accepted cause of the formation of giant planets. For simplicity, most models presently assume that the growth of planetary embryos occurs in isolation. We explore how the simultaneous growth of two embryos at the present locations of Jupiter and Saturn affects the outcome of planetary formation. We model planet formation on the basis of the core accretion scenario and include several key physical ingredients. We consider a protoplanetary gas disk that exponentially decays with time. For planetesimals, we allow for a distribution of sizes from 100~m to 100~km with most of the mass in the smaller objects. We include planetesimal migration as well as different profiles for the surface density $\Sigma$ of the disk. The core growth is computed in the framework of the oligarchic growth regime and includes the viscous enhancement of the planetesimal capture cross-section. Planet migration is ignored. By comparing calculations assuming formation of embryos in isolation to calculations with simultaneous embryo growth, we find that the growth of one embryo generally significantly affects the other. This occurs in spite of the feeding zones of each planet never overlapping. The results may be classified as a function of the gas surface density profile $\Sigma$: if $\Sigma \propto r^{-3/2}$ and the protoplanetary disk is rather massive, Jupiter's formation inhibits the growth of Saturn. If $\Sigma \propto r^{-1}$ isolated and simultaneous formation lead to very similar outcomes; in the the case of $\Sigma \propto r^{-1/2}$ Saturn grows faster and induces a density wave that later acclerates the formation of Jupiter. Our results indicate that the simultaneous growth of several embryos impacts the final outcome and should be taken into account by planet formation models.