Surface Density Profile Research Articles

New Insights into the Nature of Transition Disks from a Complete Disk Survey of the Lupus Star-forming Region

Abstract Transition disks with large dust cavities around young stars are promising targets for studying planet formation. Previous studies have revealed the presence of gas cavities inside the dust cavities, hinting at recently formed, giant planets. However, many of these studies are biased toward the brightest disks in the nearby star-forming regions, and it is not possible to derive reliable statistics that can be compared with exoplanet populations. We present the analysis of 11 transition disks with large cavities (≥20 au radius) from a complete disk survey of the Lupus star-forming region, using ALMA Band 7 observations at 0.″3 (22–30 au radius) resolution of the 345 GHz continuum, 13CO and C18O 3–2 observations, and the spectral energy distribution of each source. Gas and dust surface density profiles are derived using the physical–chemical modeling code DALI. This is the first study of transition disks of large cavities within a complete disk survey within a star-forming region. The dust cavity sizes range from 20 to 90 au radius, and in three cases, a gas cavity is resolved as well. The deep drops in gas density and large dust cavity sizes are consistent with clearing by giant planets. The fraction of transition disks with large cavities in Lupus is , which is inconsistent with exoplanet population studies of giant planets at wide orbits. Furthermore, we present a hypothesis of an evolutionary path for large massive disks evolving into transition disks with large cavities.

SDSS-IV MaNGA: the spatial distribution of star formation and its dependence on mass, structure, and environment

We study the spatially resolved star formation of 1494 galaxies in the SDSSIV-MaNGA Survey. SFRs are calculated using a two-step process, using $H_\alpha$ in star forming regions and $D_n4000$ in regions identified as AGN/LI(N)ER or lineless. The roles of secular and environmental quenching processes are investigated by studying the dependence of the radial profiles of specific star formation rate on stellar mass, galaxy structure and environment. We report on the existence of `Centrally Suppressed' galaxies, which have suppressed SSFR in their cores compared to their disks. The profiles of centrally suppressed and unsuppressed galaxies are distibuted in a bimodal way. Galaxies with high stellar mass and core velocity dispersion are found to be much more likely to be centrally suppressed than low mass galaxies, and we show that this is related to morphology and the presence of AGN/LI(N)ER like emission. Centrally suppressed galaxies also display lower star formation at all radii compared to unsuppressed galaxies. The profiles of central and satellite galaxies are also compared, and we find that satellite galaxies experience lower specific star formation rates at all radii than central galaxies. This uniform suppression could be a signal of the stripping of hot halo gas in the process known as strangulation. We find that satellites are not more likely to be suppressed in their cores than centrals, indicating that the core suppression is an entirely internal process. We find no correlation between the local environment density and the profiles of star formation rate surface density.

On the diversity and statistical properties of protostellar discs

AbstractWe present results from the first population synthesis study of protostellar discs. We analyse the evolution and properties of a large sample of protostellar discs formed in a radiation hydrodynamical simulation of star cluster formation. Due to the chaotic nature of the star formation process, we find an enormous diversity of young protostellar discs, including misaligned discs, and discs whose orientations vary with time. Star–disc interactions truncate discs and produce multiple systems. Discs may be destroyed in dynamical encounters and/or through ram-pressure stripping, but reform by later gas accretion. We quantify the distributions of disc mass and radii for protostellar ages up to ≈105 yr. For low-mass protostars, disc masses tend to increase with both age and protostellar mass. Disc radii range from of order 10 to a few hundred au, grow in size on time-scales ≲ 104 yr, and are smaller around lower mass protostars. The radial surface density profiles of isolated protostellar discs are flatter than the minimum mass solar nebula model, typically scaling as Σ ∝ r−1. Disc to protostar mass ratios rarely exceed two, with a typical range of Md/M* = 0.1–1 to ages ≲ 104 yr and decreasing thereafter. We quantify the relative orientation angles of circumstellar discs and the orbit of bound pairs of protostars, finding a preference for alignment that strengths with decreasing separation. We also investigate how the orientations of the outer parts of discs differ from the protostellar and inner disc spins for isolated protostars and pairs.

Schrödinger evolution of self-gravitating discs

An understanding of the long-term evolution of self-gravitating disks ranks among the classic outstanding problems of astrophysics. In this work, we show that the secular inclination dynamics of a geometrically thin quasi-Keplerian disk, with a surface density profile that scales as the inverse square-root of the orbital radius, are described by the time-dependent Schr\"odinger equation. Within the context of this formalism, nodal bending waves correspond to the eigenmodes of a quasiparticle's wavefunction, confined in an infinite square well with boundaries given by the radial extent of the disk. We further show that external secular perturbations upon self-gravitating disks exhibit a mathematical similarity to quantum scattering theory. Employing this framework, we derive an analytic criterion for the gravitational rigidity of a nearly-Keplerian disk under external perturbations. Applications of the theory to circumstellar disks and Galactic nuclei are discussed.

Probing the distant galaxy cluster JKCS 041 on the L − T − M scaling relations

The detailed X-ray analysis of the distant galaxy cluster JKCS 041 is presented. We use deep (~75 ks) archived data of X-ray Chandra Observatory to extract the main physical characteristic for one of the most distant galaxy cluster known to date. We investigate the imaging and spectral properties of JKCS 041. We explore its surface brightness, density, entropy, cooling time, and mass profiles. The temperature of JKCS 041 is equal to 7.4 ± 2.9 keV while the total virial mass is M200 = (4.6 ± 2.9) × 1014MSun. The gas fraction is ~10% while the dark matter is ~90% at R200. We use the obtained physical parameters of JKCS 041 to build numerous X-ray scaling relations. By adding JKCS 041 parameters we increase the redshift of our previous cluster’s sample from 1.4 to 1.8. We study the three classical relations between temperature, luminosity and total mass, and two additional. We find the concentration parameter of JKCS 041, build c − M relation and compare them with current hydrodynamic simulations. In addition, we explore, for the first time in the case of distant objects, the M − Y = T · Mg relation which is one of the most robust mass estimators. We conclude that concentration parameter, c, of JKCS 041 is in a good agreement with theoretical predictions. The obtained X-ray scaling relations were used to probe their evolution. We find that our results show inconsistent with self-similar evolution models.

The distribution of stars around the Milky Way’s central black hole

Context. The existence of dynamically relaxed stellar density cusps in dense clusters around massive black holes is a long-standing prediction of stellar dynamics, but it has so far escaped unambiguous observational confirmation. Aims. In this paper we aim to revisit the problem of inferring the innermost structure of the Milky Way’s nuclear star cluster via star counts, to clarify whether it displays a core or a cusp around the central black hole. Methods. We used judiciously selected adaptive optics assisted high angular resolution images obtained with the NACO instrument at the ESO VLT. Through image stacking and improved point spread function fitting we pushed the completeness limit about one magnitude deeper than in previous, comparable work. Crowding and extinction corrections were derived and applied to the surface density estimates. Known young, and therefore dynamically not relaxed stars, are excluded from the analysis. Contrary to previous work, we analyse the stellar density in well-defined magnitude ranges in order to be able to constrain stellar masses and ages. Results. We focus on giant stars, with observed magnitudes K = 12.5−16, and on stars with observed magnitudes K ≈ 18, which may have similar mean ages and masses than the former. The giants display a core-like surface density profile within a projected radius R ≤ 0.3 pc of the central black hole, in agreement with previous studies, but their 3D density distribution is not inconsistent with a shallow cusp if we take into account the extent of the entire cluster, beyond the radius of influence of the central black hole. The surface density of the fainter stars can be described well by a single power-law at R < 2 pc. The cusp-like profile of the faint stars persists even if we take into account the possible contamination of stars in this brightness range by young pre-main sequence stars. The data are inconsistent with a core-profile for the faint stars. Finally, we show that a 3D Nuker law provides a good description of the cluster structure. Conclusions. We conclude that the observed density of the faintest stars detectable with reasonable completeness at the Galactic centre, is consistent with the existence of a stellar cusp around the Milky Way’s central black hole, Sagittarius A*. This cusp is well developed inside the influence radius of Sagittarius A* and can be described by a single three-dimensional power-law with an exponent γ = 1.43 ± 0.02 ± 0.1sys. This corroborates existing conclusions from Nbody simulations performed in a companion paper. An important caveat is that the faint stars analysed here may be contaminated significantly by dynamically unrelaxed stars that formed about 100 Myr ago. The apparent lack of giants at projected distances of R ≲ 0.3 pc (R ≲ 8′′) of the massive black hole may indicate that some mechanism may have altered their distribution or intrinsic luminosity. We roughly estimate the number of possibly missing giants to about 100.

First Results on the Cluster Galaxy Population from the Subaru Hyper Suprime-Cam Survey. III. Brightest Cluster Galaxies, Stellar Mass Distribution, and Active Galaxies

Abstract The unprecedented depth and area surveyed by the Subaru Strategic Program with the Hyper Suprime-Cam (HSC-SSP) have enabled us to construct and publish the largest distant cluster sample out to to date. In this exploratory study of cluster galaxy evolution from z = 1 to z = 0.3, we investigate the stellar mass assembly history of brightest cluster galaxies (BCGs), the evolution of stellar mass and luminosity distributions, the stellar mass surface density profile, as well as the population of radio galaxies. Our analysis is the first high-redshift application of the top N richest cluster selection, which is shown to allow us to trace the cluster galaxy evolution faithfully. Over the 230 deg2 area of the current HSC-SSP footprint, selecting the top 100 clusters in each of the four redshift bins allows us to observe the buildup of galaxy population in descendants of clusters whose mass is about . Our stellar mass is derived from a machine-learning algorithm, which is found to be unbiased and accurate with respect to the COSMOS data. We find very mild stellar mass growth in BCGs (about 35% between z = 1 and 0.3), and no evidence for evolution in both the total stellar mass–cluster mass correlation and the shape of the stellar mass surface density profile. We also present the first measurement of the radio luminosity distribution in clusters out to , and show hints of changes in the dominant accretion mode powering the cluster radio galaxies at .

Measurement of Circumstellar Disk Sizes in the Upper Scorpius OB Association with ALMA

Abstract We present detailed modeling of the spatial distributions of gas and dust in 57 circumstellar disks in the Upper Scorpius OB Association observed with ALMA at submillimeter wavelengths. We fit power-law models to the dust surface density and CO J = 3–2 surface brightness to measure the radial extent of dust and gas in these disks. We found that these disks are extremely compact: the 25 highest signal-to-noise disks have a median dust outer radius of 21 au, assuming an dust surface density profile. Our lack of CO detections in the majority of our sample is consistent with these small disk sizes assuming the dust and CO share the same spatial distribution. Of seven disks in our sample with well-constrained dust and CO radii, four appear to be more extended in CO, although this may simply be due to the higher optical depth of the CO. Comparison of the Upper Sco results with recent analyses of disks in Taurus, Ophiuchus, and Lupus suggests that the dust disks in Upper Sco may be approximately three times smaller in size than their younger counterparts, although we caution that a more uniform analysis of the data across all regions is needed. We discuss the implications of these results for disk evolution.

ALMA Observations of Elias 2–24: A Protoplanetary Disk with Multiple Gaps in the Ophiuchus Molecular Cloud

Abstract We present ALMA 1.3 mm continuum observations at (25 au) resolution of Elias 2–24, one of the largest and brightest protoplanetary disks in the Ophiuchus Molecular Cloud, and we report the presence of three partially resolved concentric gaps located at ∼20, 52, and 87 au from the star. We perform radiative transfer modeling of the disk to constrain its surface density and temperature radial profile and place the disk structure in the context of mechanisms capable of forming narrow gaps such as condensation fronts and dynamical clearing by actively forming planets. In particular, we estimate the disk temperature at the locations of the gaps to be 23, 15, and 12 K (at 20, 52, and 87 au, respectively), very close to the expected snowlines of CO (23–28 K) and N2 (12–15 K). Similarly, by assuming that the widths of the gaps correspond to 4–8× the Hill radii of forming planets (as suggested by numerical simulations), we estimate planet masses in the range of , , and for the inner, middle, and outer gap, respectively. Given the surface density profile of the disk, the amount of “missing mass” at the location of each one of these gaps (between 4 and 20 ) is more than sufficient to account for the formation of such planets.

Deep Subaru Hyper Suprime-Cam Observations of Milky Way Satellites Columba I and Triangulum II*

Abstract We present deep, wide-field Subaru Hyper Suprime-Cam photometry of two recently discovered satellites of the Milky Way (MW): Columba I (Col I) and Triangulum II (Tri II). The color–magnitude diagrams of both objects point to exclusively old and metal-poor stellar populations. We re-derive structural parameters and luminosities of these satellites, and find for Col I and for Tri II, with corresponding half-light radii of pc and pc. The properties of both systems are consistent with observed scaling relations for MW dwarf galaxies. Based on archival data, we derive upper limits on the neutral gas content of these dwarfs, and find that they lack H i, as do the majority of observed satellites within the MW virial radius. Neither satellite shows evidence of tidal stripping in the form of extensions or distortions in matched-filter stellar density maps or surface-density profiles. However, the smaller Tri II system is relatively metal-rich for its luminosity (compared to other MW satellites), possibly because it has been tidally stripped. Through a suite of orbit simulations, we show that Tri II is approaching pericenter of its eccentric orbit, a stage at which tidal debris is unlikely to be seen. In addition, we find that Tri II may be on its first infall into the MW, which helps explain its unique properties among MW dwarfs. Further evidence that Tri II is likely an ultra-faint dwarf comes from its stellar mass function, which is similar to those of other MW dwarfs.

Tracing the assembly history of NGC 1395 through its Globular Cluster System

We used deep Gemini-South/GMOS g'r'i'z' images to study the globular cluster (GC) system of the massive elliptical galaxy NGC 1395, located in the Eridanus supergroup. The photometric analysis of the GC candidates reveals a clear colour bimodality distribution, indicating the presence of "blue" and "red" GC subpopulations. While a negative radial colour gradient is detected in the projected spatial distribution of the red GCs, the blue GCs display a shallow colour gradient. The blue GCs also display a remarkable shallow and extended surface density profile, suggesting a significant accretion of low-mass satellites in the outer halo of the galaxy. In addition, the slope of the projected spatial distribution of the blue GCs in the outer regions of the galaxy, is similar to that of the X-ray halo emission. Integrating up to 165 kpc the profile of the projected spatial distribution of the GCs, we estimated a total GC population and specific frequency of 6000$\pm$1100 and $S_N$=7.4$\pm$1.4, respectively. Regarding NGC 1395 itself, the analysis of the deep Gemini/GMOS images shows a low surface brightness umbrella-like structure indicating, at least, one recent merger event. Through relations recently published in the literature, we obtained global parameters, such as $M_\mathrm{stellar}=9.32\times10^{11}$ M$\odot$ and $M_h=6.46\times10^{13}$ M$\odot$. Using public spectroscopic data, we derive stellar population parameters of the central region of the galaxy by the full spectral fitting technique. We have found that, this region, seems to be dominated for an old stellar population, in contrast to findings of young stellar populations from the literature.

The properties of the inner disk around HL Tau: Multi-wavelength modeling of the dust emission

We conducted a detailed radiative transfer modeling of the dust emission from the circumstellar disk around HL Tau. The goal of our study is to derive the surface density profile of the inner disk and its structure. In addition to the Atacama Large Millimeter/submillimeter Array images at Band 3 (2.9mm), Band 6 (1.3mm), and Band 7 (0.87mm), the most recent Karl G. Jansky Very Large Array (VLA) observations at 7mm were included in the analysis. A simulated annealing algorithm was invoked to search for the optimum model. The radiative transfer analysis demonstrates that most radial components (i.e., >6AU) of the disk become optically thin at a wavelength of 7mm, which allows us to constrain, for the first time, the dust density distribution in the inner region of the disk. We found that a homogeneous grain size distribution is not sufficient to explain the observed images at different wavelengths simultaneously, while models with a shallower grain size distribution in the inner disk work well. We found clear evidence that larger grains are trapped in the first bright ring. Our results imply that dust evolution has already taken place in the disk at a relatively young (i.e., ~1Myr) age. We compared the midplane temperature distribution, optical depth, and properties of various dust rings with those reported previously. Using the Toomre parameter, we briefly discussed the gravitational instability as a potential mechanism for the origin of the dust clump detected in the first bright ring via the VLA observations.

Dwarf galaxy mass estimators versus cosmological simulations.

We use a suite of high-resolution cosmological dwarf galaxy simulations to test the accuracy of commonly used mass estimators from Walker et al. (2009) and Wolf et al. (2010), both of which depend on the observed line-of-sight velocity dispersion and the 2D half-light radius of the galaxy, R<inf>e</inf> . The simulations are part of the Feedback in Realistic Environments (fire) project and include 12 systems with stellar masses spanning 105–107M<inf>⊙</inf> that have structural and kinematic properties similar to those of observed dispersion-supported dwarfs. Both estimators are found to be quite accurate: and , with errors reflecting the 68 per cent range over all simulations. The excellent performance of these estimators is remarkable given that they each assume spherical symmetry, a supposition that is broken in our simulated galaxies. Though our dwarfs have negligible rotation support, their 3D stellar distributions are flattened, with short-to-long axis ratios c/a ≃ 0.4–0.7. The median accuracy of the estimators shows no trend with asphericity. Our simulated galaxies have sphericalized stellar profiles in 3D that follow a nearly universal form, one that transitions from a core at small radius to a steep fall-off ∝r−42 at large r; they are well fit by Sérsic profiles in projection. We find that the most important empirical quantity affecting mass estimator accuracy is R<inf>e</inf> . Determining R<inf>e</inf> by an analytic fit to the surface density profile produces a better estimated mass than if the half-light radius is determined via direct summation.

Radial distribution of dust, stars, gas, and star-formation rate in DustPedia face-on galaxies

The purpose of this work is the characterization of the radial distribution of dust, stars, gas, and star-formation rate (SFR) in a sub-sample of 18 face-on spiral galaxies extracted from the DustPedia sample. This study is performed by exploiting the multi-wavelength, from UV to sub-mm bands, DustPedia database, in addition to molecular (12CO) and atomic (HI) gas maps and metallicity abundance information available in the literature. We fitted the surface brightness profiles of the tracers of dust and stars, the mass surface density profiles of dust, stars, molecular gas, and total gas, and the SFR surface density profiles with an exponential curve and derived their scale-lengths. We also developed a method to solve for the CO-to-H2 conversion factor (alpha_CO) per galaxy by using dust and gas mass profiles. Although each galaxy has its own peculiar behaviour, we identified a common trend of the exponential scale-lengths vs. wavelength. On average, the scale-lengths normalized to the B-band 25 mag/arcsec^2 radius decrease from UV to 70 micron, from 0.4 to 0.2, and then increase back up to 0.3 at 500 microns. The main result is that, on average, the dust mass surface density scale-length is about 1.8 times the stellar one derived from IRAC data and the 3.6 micron surface brightness, and close to that in the UV. We found a mild dependence of the scale-lengths on the Hubble stage T: the scale-lengths of the Herschel bands and the 3.6 micron scale-length tend to increase from earlier to later types, the scale-length at 70 micron tends to be smaller than that at longer sub-mm wavelength with ratios between longer sub-mm wavelengths and 70 micron that decrease with increasing T. The scale-length ratio of SFR and stars shows a weak increasing trend towards later types.

Modelling hard and soft states of Cygnus X-1 with propagating mass accretion rate fluctuations

We present a timing analysis of threeRossi X-ray Timing Explorer observations of the black hole binary Cygnus X-1 with the propagating mass accretion rate fluctuations model PROPFLUC. The model simultaneously predicts power spectra, time lags, and coherence of the variability as a function of energy. The observations cover the soft and the hard state of the source, and the transition between the two. We find good agreement between model pre-dictions and data in the hard and in the soft state. Our analysis suggests that in the soft state the fluctuations propagate in an optically thin hot flow extending up to large radii above and below a stable optically thick disc. In the hard state, our results are consistent with a truncated disc geometry, where the hot flow extends radially inside the inner radius of the disc. In the transition from soft to hard state, the characteristics of the rapid variability are too complex to be successfully described with PROPFLUC. The surface density profile of the hot flow predicted by our model and the lack of QPOs in the soft and hard state, suggest that the spin of the black hole is aligned with the inner accretion disc and therefore probably with the rotational axis of the binary system.

The structure of young embedded protostellar discs

Young protostellar discs provide the initial conditions for planet formation. The properties of these discs may be different from those of late-phase (T Tauri) discs due to continuing infall from the envelope and protostellar variability resulting from irregular gas accretion. We use a set of hydrodynamic simulations to determine the structure of discs forming in collapsing molecular clouds. We examine how radiative feedback from the host protostar affects the disc properties by examining three regimes: without radiative feedback, with continuous radiative feedback and with episodic feedback, similar to FU Ori-type outbursts. We find that the radial surface density and temperature profiles vary significantly as the disc accretes gas from the infalling envelope. These profiles are sensitive to the presence of spiral structure, induced by gravitational instabilities, and the radiative feedback provided by the protostar, especially in the case when the feedback is episodic. We also investigate whether mass estimates from position-velocity (PV) diagrams are accurate for early-phase discs. We find that the protostellar system mass (i.e. the mass of the protostar and its disc) is underestimated by up to 20%, due to the impact of an enhanced radial pressure gradient on the gas. The mass of early-phase discs is a significant fraction of the mass of the protostar, so position-velocity diagrams cannot accurately provide the mass of the protostar alone. The enhanced radial pressure gradient expected in young discs may lead to an increased rate of dust depletion due to gas drag, and therefore to a reduced dust-to-gas ratio.

Red Clump Stars from LAMOST II: the outer disc of the Milky Way

We present stellar density maps of the Galactic outer disc with red clump stars from LAMOST data. These samples are separated into younger (mean age ∼ 2.7 Gyr) and older (mean age ∼ 4.6 Gyr) populations so that they can trace the variation of structures with ages in the range of Galactocentric radius R from 9 to 13.5 kpc. We show that scale heights for both of the two populations increase with R and display radial gradients of 48±6 and 40±4 pc kpc for the older and younger populations, respectively. It is evident that flaring occurs in the thin disc populations with a wide range of ages. Moreover, the intensity of flaring does not seem to be significantly related to the age of the thin disc populations. On the other hand, scale lengths of the radial surface density profiles are 4.7±0.5 kpc for the younger population and 3.4±0.2 kpc for the older population, meaning that the younger disc population is more radially extended than the older one. Although the fraction of the younger population mildly increases from 28% at R ∼ 9 to about 35% at R ∼ 13 kpc, the older population is prominent with a fraction of no less than 65% in the outer disc.

Dynamics of Dwarf Galaxies Disfavor Stellar-Mass Black Holes as Dark Matter.

We study the effects of black hole dark matter on the dynamical evolution of stars in dwarf galaxies. We find that mass segregation leads to a depletion of stars in the center of dwarf galaxies and the appearance of a ring in the projected stellar surface density profile. Using Segue 1 as an example we show that current observations of the projected surface stellar density rule out at the 99.9%confidence level the possibility that more than 6% of the dark matter is composed of black holes with a mass of few tens of solar masses.

Galaxy–galaxy lensing in EAGLE: comparison with data from 180 deg2 of the KiDS and GAMA surveys

We present predictions for the galaxy-galaxy lensing profile from the EAGLE hydrodynamical cosmological simulation at redshift z=0.18, in the spatial range 0.02 < R/(Mpc/h) < 2, and for five logarithmically equi-spaced stellar mass bins in the range 10.3 < $\log_{10}$(Mstar/ $M_{\odot}$) < 11.8. We compare these excess surface density profiles to the observed signal from background galaxies imaged by the Kilo Degree Survey around spectroscopically confirmed foreground galaxies from the GAMA survey. Exploiting the GAMA galaxy group catalogue, the profiles of central and satellite galaxies are computed separately for groups with at least five members to minimise contamination. EAGLE predictions are in broad agreement with the observed profiles for both central and satellite galaxies, although the signal is underestimated at R$\approx$0.5-2 Mpc/h for the highest stellar mass bins. When central and satellite galaxies are considered simultaneously, agreement is found only when the selection function of lens galaxies is taken into account in detail. Specifically, in the case of GAMA galaxies, it is crucial to account for the variation of the fraction of satellite galaxies in bins of stellar mass induced by the flux-limited nature of the survey. We report the inferred stellar-to-halo mass relation and we find good agreement with recent published results. We note how the precision of the galaxy-galaxy lensing profiles in the simulation holds the potential to constrain fine-grained aspects of the galaxy-dark matter connection.

Formation of solar system analogues – I. Looking for initial conditions through a population synthesis analysis

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.