Simple toy models are often not sufficient to cover the complexity of the dust coagulation process, and a number of numerical approaches are therefore used, among which integration of the Smoluchowski equation and various versions of Monte Carlo algorithm are the most popular. In this paper, we directly compare the Smoluchowski and Monte Carlo approaches and we find a general agreement for most of the coagulation problems. However, for the sweep-up growth driven by the "lucky" breakthrough mechanism, the methods exhibit very different resolution dependencies. With too few mass bins, the Smoluchowski algorithm tends to overestimate the growth rate and the probability of breakthrough. The Monte Carlo method is less resolution dependent in the growth timescale aspect but it tends to underestimate the breakthrough chance due to its limited dynamic mass range. We discuss the features and drawbacks of both the approaches, which may limit their astrophysical applications.

Context. The growth process of dust particles in protoplanetary disks can be modeled via numerical dust coagulation codes. In this approach, physical effects that dominate the dust growth process often must be implemented in a parameterized form. Due to a lack of these parameterizations, existing studies of dust coagulation have ignored the effects a hydrodynamical gas flow can have on grain growth, even though it is often argued that the flow could significantly contribute either positively or negatively to the growth process. Aims. We intend to provide a quantification of hydrodynamical effects on the growth of dust particles, such that these effects can be parameterized and implemented in a dust coagulation code. Methods. We numerically integrate the trajectories of small dust particles in the flow of disk gas around a proto-planetesimal, sampling a large parameter space in proto-planetesimal radii, headwind velocities, and dust stopping times. Results. The gas flow deflects most particles away from the proto-planetesimal, such that its effective collisional cross section, and therefore the mass accretion rate, is reduced. The gas flow however also reduces the impact velocity of small dust particles onto a proto-planetesimal. This can be beneficial for its growth, since large impact velocities are known to lead to erosion. We also demonstrate why such a gas flow does not return collisional debris to the surface of a proto-planetesimal. Conclusions. We predict that a laminar hydrodynamical flow around a proto-planetesimal will have a significant effect on its growth. However, we cannot easily predict which result, the reduction of the impact velocity or the sweep-up cross section, will be more important. Therefore, we provide parameterizations ready for implementation into a dust coagulation code.

The early stages of planet formation are still not well understood. Coagulation models have revealed numerous obstacles to the dust growth, such as the bouncing, fragmentation and radial drift barriers. We study the interplay between dust coagulation and drift in order to determine the conditions in protoplanetary disk that support the formation of planetesimals. We focus on planetesimal formation via sweep-up and investigate whether it can take place in a realistic protoplanetary disk. We have developed a new numerical model that resolves spatial distribution of dust in the radial and vertical dimension. The model uses representative particles approach to follow the dust evolution in protoplanetary disk. The coagulation and fragmentation of solids is taken into account using Monte Carlo method. A collision model adopting the mass transfer effect, that can occur for different-sized dust aggregate collisions, is implemented. We focus on a protoplanetary disk including a pressure bump caused by a steep decline of turbulent viscosity around the snow line. Our results show that sufficient resolution of the vertical disk structure in dust coagulation codes is necessary to obtain adequately short growth timescales, especially in the case of a low turbulence region. We find that a sharp radial variation of the turbulence strength at the inner edge of dead zone promotes planetesimal formation in several ways. It provides a pressure bump that efficiently prevents the dust from drifting inwards. It also causes a radial variation in the size of aggregates at which growth barriers occur, favoring the growth of large aggregates via sweeping up of small particles. In our model, by employing an ad hoc alpha viscosity change near the snow line, it is possible to grow planetesimals by incremental growth on timescales of approximately 10^5 years.

It is unknown how far dust growth can proceed by coagulation. Obstacles to collisional growth are the fragmentation and bouncing barriers. However, in all previous simulations of the dust-size evolution in protoplanetary disks, only the mean collision velocity has been considered, neglecting that a small but possibly important fraction of the collisions will occur at both much lower and higher velocities. We study the effect of the probability distribution of impact velocities on the collisional dust growth barriers. Assuming a Maxwellian velocity distribution for colliding particles to determine the fraction of sticking, bouncing, and fragmentation, we implement this in a dust-size evolution code. We also calculate the probability of growing through the barriers and the growth timescale in these regimes. We find that the collisional growth barriers are not as sharp as previously thought. With the existence of low-velocity collisions, a small fraction of the particles manage to grow to masses orders of magnitude above the main population. A particle velocity distribution softens the fragmentation barrier and removes the bouncing barrier. It broadens the size distribution in a natural way, allowing the largest particles to become the first seeds that initiate sweep-up growth towards planetesimal sizes.

The formation of planetesimals is often accredited to collisional sticking of dust grains. The exact process is unknown, as collisions between larger aggregates tend to lead to fragmentation or bouncing rather than sticking. Recent laboratory experiments have however made great progress in the understanding and mapping of the complex physics involved in dust collisions. We want to study the possibility of planetesimal formation using the results from the latest laboratory experiments, particularly by including the fragmentation with mass transfer effect, which might lead to growth even at high impact velocities. We present a new experimentally and physically motivated dust collision model capable of predicting the outcome of a collision between two particles of arbitrary masses and velocities. It is used together with a continuum dust-size evolution code that is both fast in terms of execution time and able to resolve the dust well at all sizes, allowing for all types of interactions to be studied without biases. We find that for the general dust population, bouncing collisions prevent the growth above millimeter-sizes. However, if a small number of cm-sized particles are introduced, for example due to vertical mixing or radial drift, they can act as a catalyst and start to sweep up the smaller particles. At a distance of 3 AU, 100-meter-sized bodies are formed on a timescale of 1 Myr. We conclude that direct growth of planetesimals might be a possibility thanks to a combination of the existence of a bouncing barrier and the fragmentation with mass transfer effect. The bouncing barrier is here even beneficial, as it prevents the growth of too many large particles that would otherwise only fragment among each other, and creates a reservoir of small particles that can be swept up by larger bodies. However, for this process to work, a few seeds of cm in size or larger have to be introduced.

Context. The Gaia satellite will measure highly accurate absolute parallaxes of hundreds of millions of stars by comparing the parallactic displacements in the two fields of view of the optical instrument. The requirements on the stability of the 'basic angle' between the two fields are correspondingly strict, and possible variations (on the microarcsec level) are therefore monitored by an on-board metrology system. Nevertheless, since even very small periodic variations of the basic angle might cause a global offset of the measured parallaxes, it is important to find independent verification methods. Aims. We investigate the potential use of Galactic Cepheids as standard candles for verifying the Gaia parallax zero point. Methods. We simulate the complete population of Galactic Cepheids and their observations by Gaia. Using the simulated data, simultaneous fits are made of the parameters of the period-luminosity relation and a global parallax zero point. Results. The total number of Galactic Cepheids is estimated at about 20 000, of which nearly half could be observed by Gaia. In the most favourable circumstances, including negligible intrinsic scatter and extinction errors, the determined parallax zero point has an uncertainty of 0.2 microarcsec. With more realistic assumptions the uncertainty is several times larger, and the result is very sensitive to errors in the applied extinction corrections. Conclusions. The use of Galactic Cepheids alone will not be sufficient to determine a possible parallax zero-point error to the full potential systematic accuracy of Gaia. The global verification of Gaia parallaxes will most likely depend on a combination of many different methods, including this one.