In this work, we present results that simultaneously constrain the dust size distribution, dust-to-gas ratio, fraction of dust re-deposition, and total mass production rates for comet 67P/Churyumov-Gerasimenko. We use a 3D Direct Simulation Monte Carlo (DSMC) gas dynamics code to simulate the inner gas coma of the comet for the duration of the Rosetta mission. The gas model is constrained by ROSINA/COPS data. Further, we simulate for different epochs the inner dust coma using a 3D dust dynamics code including gas drag and the nucleus' gravity. Using advanced dust scattering properties these results are used to produce synthetic images that can be compared to the OSIRIS data set. These simulations allow us to constrain the properties of the dust coma and the total gas and dust production rates. We determined a total volatile mass loss of $(6.1 \pm 1.5) \cdot 10^9$~kg during the 2015 apparition. Further, we found that power-laws with $q=3.7^{+0.57}_{-0.078}$ are consistent with the data. This results in a total of $5.1^{+6.0}_{-4.9}\cdot10^9$~kg of dust being ejected from the nucleus surface, of which $4.4^{+4.9}_{-4.2}\cdot10^9$~kg escape to space and $6.8^{+11}_{-6.8}\cdot10^8$~kg (or an equivalent of $14^{+22}_{-14}$~cm over the smooth regions) is re-deposited on the surface. This leads to a dust-to-gas ratio of $0.73^{+1.3}_{-0.70}$ for the escaping material and $0.84^{+1.6}_{-0.81}$ for the ejected material. We have further found that the smallest dust size must be strictly smaller than $\sim30\mu$m and nominally even smaller than $\sim12\mu$m.

Context: Remote light scattering and thermal infrared observations provide clues about the physical properties of cometary and interplanetary dust particles. Identifying these properties will lead to a better understanding of the formation and evolution of the Solar System. Aims: We present a numerical solution for the radiative and conductive heat transport in a random particulate medium enclosed by an arbitrarily shaped surface. The method will be applied to study thermal properties of cometary dust particles. Methods: The recently introduced incoherent Monte Carlo radiative transfer method developed for scattering, absorption, and propagation of electromagnetic waves in dense discrete random media is extended for radiative heat transfer and thermal emission. The solution is coupled with the conductive Fourier transport equation that is solved with the finite-element method. Results: The proposed method allows the synoptic analysis of light scattering and thermal emission by large cometary dust particles consisting of submicrometer-sized grains. In particular, we show that these particles can sustain significant temperature gradients resulting in the superheating factor phase function observed for the coma of comet 67P/Churyumov-Gerasimenko.

We show that the scattering phase functions of the coma and the nucleus of the comet 67P/Churyumov-Gerasimenko measured by the Rosetta/OSIRIS instrument can be reproduced by a particle model involving clustered densely packed submicrometer-sized grains composed of organic material and larger micrometer-sized silicate grains. The simulated and measured coma phase functions suggest that near the nucleus scattering is dominated by large particles, and the size distribution of dust particles varies with time and/or local coma environment. Further, we show that the measured nucleus phase function is consistent with the coma phase function by modelling a nucleus-sized object consisting of the same particles that explain the coma phase functions.

Weak interstellar scintillations of pulsar B0809+74 were observed at two epochs using a 30m EISCAT antenna at 933 MHz. These have been used to constrain the spectrum, the distribution and the transverse velocity of the scattering plasma with respect to the local standard of rest (LSR). The Kolmogorov power law is a satisfactory model for the electron density spectrum at scales between 20 megameters and 1 gigameter. We compare the observations with model calculations from weak scintillation theory and the known transverse velocities of the pulsar and the Earth. The simplest model is that the scattering is uniformly distributed along the 310 pc line of sight (l=140 deg, b=32 deg) and is stationary in the LSR. With the scattering measure as the only free parameter, this model fits the data within the errors and a range of about +/-10 km/s in velocity is also allowed. The integrated level of turbulence is low, being comparable to that found toward PSR B0950+08, and suggests a region of low local turbulence over as much as 90 deg. in longitude including the galactic anti-center. If, on the other hand, the scattering occurs in a compact region, the observed time scales require a specific velocity-distance relation. In particular, enhanced scattering in a shell at the edge of the local bubble, proposed by Bhat et al. (1998), near 72 pc toward the pulsar, must be moving at about ~ 17 km/s; however, the low scattering measure argues against a shell of enhanced scattering in this direction. The analysis also excludes scattering in the termination shock of the solar wind or in a nebula associated with the pulsar.