Elemental abundances, particularly the C/O ratio, are seen as a way to connect the composition of planetary atmospheres with planet formation scenario and the disc chemical environment. We model the chemical composition of gas and ices in a self-gravitating disc on timescales of 0.5 Myr since its formation to study the evolution of C/O ratio due to dust dynamics and growth and phase transitions of the volatile species. We use the thin-disc hydrodynamic code FEOSAD, which includes disc self-gravity, thermal balance, dust evolution, and turbulent diffusion, and treats dust as a dynamically different and evolving component interacting with the gas. It also describes freeze-out, sublimation, and advection of four most abundant volatile species: H$_2$O, CO$_2$, CH$_4$, and CO. We demonstrate the effect of gas and dust substructures such as spirals and rings on the distribution of volatiles and C/O ratios, including the formation of multiple snowlines of one species, and point out the anticorrelation between dust-to-gas ratio and total C/O ratio emerging due to the contribution of oxygen-rich ice mantles. We identify time and spatial locations where two distinct trigger mechanisms for planet formation are operating and differentiate them by C/O ratio range: wide range of the C/O ratios of $0-1.4$ for streaming instability, and a much narrower range $0.3-0.6$ for gravitational instability (with the initial value of 0.34). This conclusion is corroborated by observations, showing that transiting exoplanets, which possibly experienced migration through a variety of disc conditions, have significantly larger spread of C/O in comparison with directly imaged exoplanets likely formed in gravitationally unstable outer disk regions. We show that the ice-phase $\textrm{C/O}\approx$0.2–0.3 between the CO, CO$_2$, and CH$_4$ snowlines corresponds to the composition of the Solar system comets, that represent primordial planetesimals.

Planets form in protoplanetary accretion discs around young protostars. These discs are driven by internal turbulence and the gas flow is not laminar but has stochastic components. For weakly ionised discs the turbulence can be generated purely hydrodynamically through the vertical shear instability (VSI). Embedded particles (dust/pebbles) experience a hydrodynamic drag and drift inward radially and are stirred up vertically by the turbulent motion of the disc. We study the accretion of particles onto a forming planet embedded in a VSI turbulent protoplanetary disc through a series of 3D hydrodynamical simulations for locally isothermal discs with embedded planets in the mass range from 5 to 100 Earth masses (M2295).