We compute inclusive electron-nucleus cross sections using ab initio spectral functions of $^4$He and $^{16}$O obtained within the Self Consistent Green's Function approach. The formalism adopted is based on the factorization of the spectral function and the nuclear transition matrix elements. This allows to provide an accurate description of nuclear dynamics and to account for relativistic effects in the interaction vertex. Our calculations use a saturating chiral Hamiltonian in order reproduce the correct nuclear sizes. When final state interactions for the struck particle are accounted for, we find nice agreement between the data and the theory for the inclusive electron-$^{16}$O cross section. The results lay the foundations for future applications of the Self Consistent Green's Function method, in both closed and open shell nuclei, to neutrino data analysis. This work also presents results for the point-proton, charge and single-nucleon momentum distribution of the same two nuclei. The center of mass can affect these quantities for light nuclei and cannot be separated cleanly in most ab initio post-Hartree-Fock methods. In order to address this, we developed a Metropolis Monte Carlo calculation in which the center of mass coordinate can be subtracted exactly from the trial wave function and the expectation values. We gauged this effect for $^4$He by removing the center of mass effect from the Optimal Reference State wave function that is generated during the Self Consistent Green's Function calculations. Our findings clearly indicate that the residual center of mass contribution strongly modifies calculated matter distributions with respect to those obtained in the intrinsic frame. Hence, its subtraction is crucial for a correct description of light nuclei.

We perform ab initio self-consistent Green's function calculations of the closed shell nuclei $^{\rm 4}$He, $^{\rm 16}$O and $^{\rm 40}$Ca, based on two-nucleon potentials derived from Lattice QCD simulations, in the flavor SU(3) limit and at the pseudo-scalar meson mass of 469~MeV/c$^{\rm 2}$. The nucleon-nucleon interaction is obtained using the HAL QCD method and its short-distance repulsion is treated by means of ladder resummations outside the model space. Our results show that this approach diagonalises ultraviolet degrees of freedom correctly. Therefore, ground state energies can be obtained from infrared extrapolations even for the relatively hard potentials of HAL QCD. Comparing to previous Brueckner Hartree-Fock calculations, the total binding energies are sensibly improved by the full account of many-body correlations. The results suggest an interesting possible behaviour in which nuclei are unbound at very large pion masses and islands of stability appear at first around the traditional doubly-magic numbers when the pion mass is lowered toward its physical value. The calculated one-nucleon spectral distributions are qualitatively close to those of real nuclei even for the pseudo-scalar meson mass considered here.

A precision mass investigation of the neutron-rich titanium isotopes $^{51-55}$Ti was performed at TRIUMF's Ion Trap for Atomic and Nuclear science (TITAN). The range of the measurements covers the $N=32$ shell closure and the overall uncertainties of the $^{52-55}$Ti mass values were significantly reduced. Our results confirm the existence of a weak shell effect at $N=32$, establishing the abrupt onset of this shell closure. Our data were compared with state-of-the-art \textit{ab-initio} shell model calculations which, despite very successfully describing where the $N=32$ shell gap is strong, overpredict its strength and extent in titanium and heavier isotones. These measurements also represent the first scientific results of TITAN using the newly commissioned Multiple-Reflection Time-of-Flight Mass Spectrometer (MR-TOF-MS), substantiated by independent measurements from TITAN's Penning trap mass spectrometer.

We directly measure twenty overhanging cliffs on the surface of comet 67P/Churyumov-Gerasimenko extracted from the latest shape model and estimate the minimum tensile strengths needed to support them against collapse under the comet's gravity. We find extremely low strengths of around one Pa or less (one to five Pa, when scaled to a metre length). The presence of eroded material at the base of most overhangs, as well as the observed collapse of two features and implied previous collapse of another, suggests that they are prone to failure and that true material strengths are close to these lower limits (although we only consider static stresses and not dynamic stress from, for example, cometary activity). Thus, a tensile strength of a few pascals is a good approximation for the tensile strength of 67P's nucleus material, which is in agreement with previous work. We find no particular trends in overhang properties with size, over the $\sim10-100$ m range studied here, or location on the nucleus. There are no obvious differences, in terms of strength, height or evidence of collapse, between the populations of overhangs on the two cometary lobes, suggesting that 67P is relatively homogenous in terms of tensile strength. Low material strengths are supportive of cometary formation as a primordial rubble pile or by collisional fragmentation of a small (tens of km) body.

We present a statistical analysis of the distribution of large scale topographic features on comet 67P/Churyumov-Gerasimenko. We observe that the cumulative cliff height distribution across the surface follows a power law with a slope equal to -1.69 +- 0.02. When this distribution is studied independently for each region, we find a good correlation between the slope of the power law and the orbital erosion rate of the surface. For instance, the northern hemisphere topography is dominated by structures on the 100~m scale while the southern hemisphere topography, illuminated at perihelion, is dominated by 10~m scale terrain features. Our study suggest that the current size of a cliff is controlled not only by material cohesion but by the dominant erosional process in each region. This observation can be generalized to other comets, where we argue that primitive nuclei are characterized by the presence of large cliffs with a cumulative height power index equal to or above -1.5, while older, eroded cometary surfaces have a power index equal to or below -2.3. In effect, our model shows that a measure of the topography provides a quantitative assessment of a comet's erosional history, i.e. its evolutionary age.

We collect observational evidence that supports the scheme of mass transfer on the nucleus of comet 67P/Churyumov-Gerasimenko. The obliquity of the rotation axis of 67P causes strong seasonal variations. During perihelion the southern hemisphere is four times more active than the north. Northern territories are widely covered by granular material that indicates back fall originating from the active south. Decimetre sized chunks contain water ice and their trajectories are influenced by an anti-solar force instigated by sublimation. OSIRIS observations suggest that up to 20 % of the particles directly return to the nucleus surface taking several hours of travel time. The back fall covered northern areas are active if illuminated but produce mainly water vapour. The decimetre chunks from the nucleus surface are too small to contain more volatile compounds such as CO 2 or CO. This causes a north-south dichotomy of the composition measurements in the coma. Active particles are trapped in the gravitational minimum of Hapi during northern winter. They are "shock frozen" and only reactivated when the comet approaches the sun after its aphelion passage. The insolation of the big cavity is enhanced by self-heating, i. e. reflection and IR radiation from the walls. This, together with the pristinity of the active back fall, explains the early observed activity of the Hapi region. Sobek may be a role model for the consolidated bottom of Hapi. Mass transfer in the case of 67P strongly influences the evolution of the nucleus and the interpretation of coma measurements.

The Southern hemisphere of the 67P/Churyumov-Gerasimenko comet has become visible from Rosetta only since March 2015. It was illuminated during the perihelion passage and therefore it contains the regions that experienced the strongest heating and erosion rate, thus exposing the subsurface most pristine material. In this work we investigate, thanks to the OSIRIS images, the geomorphology, the spectrophotometry and some transient events of two Southern hemisphere regions: Anhur and part of Bes. Bes is dominated by outcropping consolidated terrain covered with fine particle deposits, while Anhur appears strongly eroded with elongated canyon-like structures, scarp retreats, different kinds of deposits, and degraded sequences of strata indicating a pervasive layering. We discovered a new 140 m long and 10 m high scarp formed in the Anhur/Bes boundary during/after the perihelion passage, close to the area where exposed CO$_2$ and H$_2$O ices were previously detected. Several jets have been observed originating from these regions, including the strong perihelion outburst, an active pit, and a faint optically thick dust plume. We identify several areas with a relatively bluer slope (i.e. a lower spectral slope value) than their surroundings, indicating a surface composition enriched with some water ice. These spectrally bluer areas are observed especially in talus and gravitational accumulation deposits where freshly exposed material had fallen from nearby scarps and cliffs. The investigated regions become spectrally redder beyond 2 au outbound when the dust mantle became thicker, masking the underlying ice-rich layers.

Images of the nucleus and the coma (gas and dust) of comet 67P/Churyumov- Gerasimenko have been acquired by the OSIRIS (Optical, Spectroscopic, and Infrared Remote Imaging System) cameras since March 2014 using both the Wide Angle Camera (WAC) and the Narrow Angle Camera (NAC). We use images from the NAC camera to study a bright outburst observed in the southern hemisphere on July 29, 2015. The high spatial resolution of the NAC is needed to localize the source point of the outburst on the surface of the nucleus. The heliocentric distance is 1.25 au and the spacecraft-comet distance is 186 km. Aiming to better understand the physics that led to the outgassing, we used the Direct Simulation Monte Carlo (DSMC) method to study the gas flow close to the nucleus and the dust trajectories. The goal is to understand the mechanisms producing the outburst. We reproduce the opening angle of the outburst in the model and constrain the outgassing ratio between the outburst source and the local region. The outburst is in fact a combination of both gas and dust, in which the active surface is approximately 10 times more active than the average rate found in the surrounding areas. We need a number of dust particles 7.83 $\times$ 10$^{11}$ - 6.90 $\times$ 10$^{15}$ (radius 1.97 - 185 {\mu}m), which corresponds to a mass of dust 220 - 21 $\times$ 10$^{3}$kg.

The ESA Rosetta spacecraft has been tracking its target, the Jupiter-family comet 67P/Churyumov-Gerasimenko, in close vicinity for over two years. It hosts the OSIRIS instruments: the Optical, Spectroscopic, and Infrared Remote Imaging System composed of two cameras, see e.g. Keller et al. (2007). In some imaging sequences dedicated to observe dust particles in the comet's coma, the two cameras took images at the same time. The aim of this work is to use these simultaneous double camera observations to calculate the dust particles' distance to the spacecraft. As the two cameras are mounted on the spacecraft with an offset of 70 cm, the distance of particles observed by both cameras can be determined by a shift of the particles' apparent trails on the images. This paper presents first results of the ongoing work, introducing the distance determination method for the OSIRIS instrument and the analysis of an example particle. We note that this method works for particles in the range of about 500 m - 6000 m from the spacecraft.

The inclusion of the three-nucleon forces (3NFs) in \textit{ab initio} many-body approaches is a formidable task, due to the computational load implied by the treatment of their matrix elements. For this reason, practical applications have mostly been limited to contributions where 3NFs enter as effective two-nucleon interactions. In this contribution, we derive the algebraic diagrammatic construction (ADC) working equations for a specific Feynman diagram of the self-energy that contains a fully irreducible three-nucleon force. This diagram is expected to be the most important among those previously neglected, because it connects dominant excited intermediate state configurations.

We show the first results for the elastic scattering of neutrons off oxygen and calcium isotopes obtained from ab initio optical potentials. The potential is derived using self consistent Green's function theory (SCGF) with the saturating chiral interaction NNLO$_{\textrm{sat}}$. Our calculations are compared to available scattering data and show that it is possible to reproduce low energy scattering observables in medium mass nuclei from first principles.

The possibility that an unconventional depletion in the center of the charge density distribution of certain nuclei occurs due to a purely quantum mechanical effect has attracted theoretical and experimental attention in recent years. We report on ab initio self-consistent Green's function calculations of one of such candidates, $^{34}$Si, together with its Z+2 neighbour $^{36}$S. Binding energies, rms radii and density distributions of the two nuclei as well as low-lying spectroscopy of $^{35}$Si, $^{37}$S, $^{33}$Al and $^{35}$P are discussed. The interpretation of one-nucleon removal and addition spectra in terms of the evolution of the underlying shell structure is also provided. The study is repeated using several chiral effective field theory Hamiltonians as a way to test the robustness of the results with respect to input inter-nucleon interactions. The prediction regarding the (non-)existence of the bubble structure in $^{34}$Si varies significantly with the nuclear Hamiltonian used. However, demanding that the experimental charge density distribution and the root mean square radius of $^{36}$S are well reproduced, along with $^{34}$Si and $^{36}$S binding energies, only leaves the NNLO$_{\text{sat}}$ Hamiltonian as a serious candidate to perform this prediction. In this context, a bubble structure, whose fingerprint should be visible in an electron scattering experiment of $^{34}$Si, is predicted. Furthermore, a clear correlation is established between the occurrence of the bubble structure and the weakening of the 1/2$^-$-3/2$^-$ splitting in the spectrum of $^{35}$Si as compared to $^{37}$S.

During its two years mission around comet 67P/Churyumov-Gerasimenko, ESA's Rosetta spacecraft had the unique opportunity to follow closely a comet in the most active part of its orbit. Many studies have presented the typical features associated to the activity of the nucleus, such as localized dust and gas jets. Here we report on series of more energetic transient events observed during the three months surrounding the comet's perihelion passage in August 2015. We detected and characterized 34 outbursts with the Rosetta cameras, one every 2.4 nucleus rotation. We identified 3 main dust plume morphologies associated to these events: a narrow jet, a broad fan, and more complex plumes featuring both previous types together. These plumes are comparable in scale and temporal variation to what has been observed on other comets. We present a map of the outbursts source locations, and discuss the associated topography. We find that the spatial distribution sources on the nucleus correlates well with morphological region boundaries, especially in areas marked by steep scarps or cliffs. Outbursts occur either in the early morning or shortly after the local noon, indicating two potential processes: Morning outbursts may be triggered by thermal stresses linked to the rapid change of temperature, afternoon events are most likely related to the diurnal or seasonal heat wave reaching volatiles buried under the first surface layer. In addition, we propose that some events can be the result of a completely different mechanism, in which most of the dust is released upon the collapse of a cliff.

Since a number of years our group is engaged in the design, construction and operation of instruments with very high time resolution in the optical band for applications to Quantum Astronomy and more conventional Astrophysics. Two instruments were built to perform photon counting with sub-nanosecond temporal accuracy. The first of the two, Aqueye+, is regularly mounted at the 1.8 m Copernicus telescope in Asiago, while the second one, Iqueye, was mounted at the ESO New Technology Telescope in Chile, and at the William Herschel Telescope and Telescopio Nazionale Galileo on the Roque (La Palma, Canary Islands). Both instruments deliver extraordinarily accurate results in optical pulsar timing. Recently, Iqueye was moved to Asiago to be mounted at the 1.2 m Galileo telescope to attempt, for the first time ever, experiments of optical intensity interferometry (a` la Hanbury Brown and Twiss) on a baseline of a few kilometers, together with the Copernicus telescope. This application was one of the original goals for the development of our instrumentation. To carry out these measurements, we are experimenting a new way of coupling the instruments to the telescopes, by means of moderate-aperture, low-optical-attenuation multi-mode optical fibers with a double-clad design. Fibers are housed in dedicated optical interfaces attached to the focus of another instrument of the 1.8 m telescope (Aqueye+) or to the Nasmyth focus of the 1.2 m telescope (Iqueye). This soft-mount solution has the advantage to facilitate the mounting of the photon counters, to keep them under controlled temperature and humidity conditions (reducing potential systematics related to varying ambient conditions), and to mitigate scheduling requirements. Here we will describe the first successful implementation of the Asiago intensity interferometer and future plans for improving it.

Since the orbital insertion of the Rosetta spacecraft, comet 67P/Churyumov-Gerasimenko (67P/C-G) has been mapped by OSIRIS camera and VIRTIS spectro-imager, producing a huge quantity of images and spectra of the comet's nucleus. The aim of this work is to search for the presence of H$_2$O on the nucleus which, in general, appears very dark and rich in dehydrated organic material. After selecting images of the bright spots which could be good candidates to search for H$_2$O ice, taken at high resolution by OSIRIS, we check for spectral cubes of the selected coordinates to identify these spots observed by VIRTIS. The selected OSIRIS images were processed with the OSIRIS standard pipeline and corrected for the illumination conditions for each pixel using the Lommel-Seeliger disk law. The spots with higher I/F were selected and then analysed spectrophotometrically and compared with the surrounding area. We selected 13 spots as good targets to be analysed by VIRTIS to search for the 2 micron absorption band of water ice in the VIRTIS spectral cubes. Out of the 13 selected bright spots, eight of them present positive H$_2$O ice detection on the VIRTIS data. A spectral analysis was performed and the approximate temperature of each spot was computed. The H$_2$O ice content was confirmed by modeling the spectra with mixing (areal and intimate) of H$_2$O ice and dark terrain, using Hapke's radiative transfer modeling. We also present a detailed analysis of the detected spots.

Beginning in March 2014, the OSIRIS (Optical, Spectroscopic, and Infrared Remote Imaging System) cameras began capturing images of the nucleus and coma (gas and dust) of comet 67P/Churyumov-Gerasimenko using both the wide angle camera (WAC) and the narrow angle camera (NAC). The many observations taken since July of 2014 have been used to study the morphology, location, and temporal variation of the comet's dust jets. We analyzed the dust monitoring observations shortly after the southern vernal equinox on May 30 and 31, 2015 with the WAC at the heliocentric distance Rh = 1.53 AU, where it is possible to observe that the jet rotates with the nucleus. We found that the decline of brightness as a function of the distance of the jet is much steeper than the background coma, which is a first indication of sublimation. We adapted a model of sublimation of icy aggregates and studied the effect as a function of the physical properties of the aggregates (composition and size). The major finding of this article was that through the sublimation of the aggregates of dirty grains (radius a between 5 microm and 50 microm) we were able to completely reproduce the radial brightness profile of a jet beyond 4 km from the nucleus. To reproduce the data we needed to inject a number of aggregates between 8.5 x $10^{13}$ and 8.5 x $10^{10}$ for a = 5 microm and 50 microm respectively, or an initial mass of $H_2O$ ice around 22kg.

We present OSIRIS/NAC observations of decimetre-sized, likely ice-containing aggregates ejected from a confined region on the surface of comet 67P/Churyumov-Gerasimenko. The images were obtained in January 2016 when the comet was at 2 AU from the Sun out-bound from perihelion. We measure the acceleration of individual aggregates through a two-hour image series. Approximately 50% of the aggregates are accelerated away from the nucleus, and 50% towards it, and likewise towards either horizontal direction. The accelerations are up to one order of magnitude stronger than local gravity, and are most simply explained by the combined effect of gas drag accelerating all aggregates upwards, and the recoil force from asymmetric outgassing, either from rotating aggregates with randomly oriented spin axes and sufficient thermal inertia to shift the temperature maximum away from an aggregate's subsolar region, or from aggregates with variable ice content. At least 10% of the aggregates will escape the gravity field of the nucleus and feed the comet's debris trail, while others may fall back to the surface and contribute to the deposits covering parts of the northern hemisphere. The rocket force plays a crucial role in pushing these aggregates back towards the surface. Our observations show the future back fall material in the process of ejection, and provide the first direct measurement of the acceleration of aggregates in the innermost coma (<2km) of a comet, where gas drag is still significant.

We present a systematic study of both nuclear radii and binding energies in (even) oxygen isotopes from the valley of stability to the neutron drip line. Both charge and matter radii are compared to state-of-the-art {\it ab initio} calculations along with binding energy systematics. Experimental matter radii are obtained through a complete evaluation of the available elastic proton scattering data of oxygen isotopes. We show that, in spite of a good reproduction of binding energies, {\it ab initio} calculations with conventional nuclear interactions derived within chiral effective field theory fail to provide a realistic description of charge and matter radii. A novel version of two- and three-nucleon forces leads to considerable improvement of the simultaneous description of the three observables for stable isotopes, but shows deficiencies for the most neutron-rich systems. Thus, crucial challenges related to the development of nuclear interactions remain.

The Wide Angle Camera of the OSIRIS instrument on board the Rosetta spacecraft is equipped with several narrowband filters that are centered on the emission lines and bands of various fragment species. These are used to determine the evolution of the production and spatial distribution of the gas in the inner coma of comet 67P with time and heliocentric distance, here between 2.6 - 1.3 AU pre-perihelion. Our observations indicate that the emission observed in the OH, OI, CN, NH, and NH2 filters is mostly produced by dissociative electron impact excitation of different parent species. We conclude that CO2 rather than H2O is a significant source of the [OI] 630 nm emission. A strong plume-like feature observed in the in CN and [OI] filters is present throughout our observations. This plume is not present in OH emission and indicates a local enhancement of the CO2/H2O ratio by as much as a factor of 3. We observed a sudden decrease in intensity levels after March 2015, which we attribute to decreased electron temperatures in the first kilometers above the nucleus surface.

The ESA Rosetta spacecraft, currently orbiting around comet 67P, has already provided in situ measurements of the dust grain properties from several instruments, particularly OSIRIS and GIADA. We propose adding value to those measurements by combining them with ground-based observations of the dust tail to monitor the overall, time-dependent dust-production rate and size distribution. To constrain the dust grain properties, we take Rosetta OSIRIS and GIADA results into account, and combine OSIRIS data during the approach phase (from late April to early June 2014) with a large data set of ground-based images that were acquired with the ESO Very Large Telescope (VLT) from February to November 2014. A Monte Carlo dust tail code has been applied to retrieve the dust parameters. Key properties of the grains (density, velocity, and size distribution) were obtained from Rosetta observations: these parameters were used as input of the code to considerably reduce the number of free parameters. In this way, the overall dust mass-loss rate and its dependence on the heliocentric distance could be obtained accurately. The dust parameters derived from the inner coma measurements by OSIRIS and GIADA and from distant imaging using VLT data are consistent, except for the power index of the size-distribution function, which is $\alpha$=--3, instead of $\alpha$=--2, for grains smaller than 1 mm. This is possibly linked to the presence of fluffy aggregates in the coma. The onset of cometary activity occurs at approximately 4.3 au, with a dust production rate of 0.5 kg/s, increasing up to 15 kg/s at 2.9 au. This implies a dust-to-gas mass ratio varying between 3.8 and 6.5 for the best-fit model when combined with water-production rates from the MIRO experiment.

We analyze the physical properties and dynamical origin of a curved jet of comet 67P/Churyumov-Gerasimenko that was observed repeatedly in several nucleus rotations starting on May 30 and persisting until early August, 2015. We simulated the motion of dust grains ejected from the nucleus surface under the influence of the gravity and viscous drag effect of the expanding gas flow from the rotating nucleus. The formation of the curved jet is a combination of the size of the dust particles (~0.1-1 mm) and the location of the source region near the nucleus equator. This enhances the spiral feature of the collimated dust stream after the dust is accelerated to a terminal speed on the order of m/s.

Optical observations provide convincing evidence that the optical phase of the Crab pulsar follows the radio one closely. Since optical data do not depend on dispersion measure variations, they provide a robust and independent confirmation of the radio timing solution. The aim of this paper is to find a global mathematical description of Crab pulsar's phase as a function of time for the complete set of published Jodrell Bank radio ephemerides (JBE) in the period 1988-2014. We apply the mathematical techniques developed for analyzing optical observations to the analysis of JBE. We break the whole period into a series of episodes and express the phase of the pulsar in each episode as the sum of two analytical functions. The first function is the best-fitting local braking index law, and the second function represents small residuals from this law with an amplitude of only a few turns, which rapidly relaxes to the local braking index law. From our analysis, we demonstrate that the power law index undergoes "instantaneous" changes at the time of observed jumps in rotational frequency (glitches). We find that the phase evolution of the Crab pulsar is dominated by a series of constant braking law episodes, with the braking index changing abruptly after each episode in the range of values between 2.1 and 2.6. Deviations from such a regular phase description behave as oscillations triggered by glitches and amount to fewer than 40 turns during the above period, in which the pulsar has made more than 2.0e10 turns. Our analysis does not favor the explanation that glitches are connected to phenomena occurring in the interior of the pulsar. On the contrary, timing irregularities and changes in slow down rate seem to point to electromagnetic interaction of the pulsar with the surrounding environment.

Dust jets, i.e. fuzzy collimated streams of cometary material arising from the nucleus, have been observed in-situ on all comets since the Giotto mission flew by comet 1P/Halley in 1986. Yet their formation mechanism remains unknown. Several solutions have been proposed, from localized physical mechanisms on the surface/sub-surface (see review in Belton (2010)) to purely dynamical processes involving the focusing of gas flows by the local topography (Crifo et al. 2002). While the latter seems to be responsible for the larger features, high resolution imagery has shown that broad streams are composed of many smaller features (a few meters wide) that connect directly to the nucleus surface. We monitored these jets at high resolution and over several months to understand what are the physical processes driving their formation, and how this affects the surface. Using many images of the same areas with different viewing angles, we performed a 3-dimensional reconstruction of collimated jets, and linked them precisely to their sources on the nucleus. Results.We show here observational evidence that the Northern hemisphere jets of comet 67P arise from areas with sharp topographic changes and describe the physical processes involved. We propose a model in which active cliffs are the main source of jet-like features, and therefore the regions eroding the fastest on comets. We suggest that this is a common mechanism taking place on all comets.

We study the link between gravitational slopes and the surface morphology on the nucleus of comet 67P/Churyumov-Gerasimenko and provide constraints on the mechanical properties of the cometary material. We computed the gravitational slopes for five regions on the nucleus that are representative of the different morphologies observed on the surface, using two shape models computed from OSIRIS images by the stereo-photoclinometry (SPC) and stereo-photogrammetry (SPG) techniques. We estimated the tensile, shear, and compressive strengths using different surface morphologies and mechanical considerations. The different regions show a similar general pattern in terms of the relation between gravitational slopes and terrain morphology: i) low-slope terrains (0-20 deg) are covered by a fine material and contain a few large ($>$10 m) and isolated boulders, ii) intermediate-slope terrains (20-45 deg) are mainly fallen consolidated materials and debris fields, with numerous intermediate-size boulders from $<$1 m to 10 m for the majority of them, and iii) high-slope terrains (45-90 deg) are cliffs that expose a consolidated material and do not show boulders or fine materials. The best range for the tensile strength of overhangs is 3-15 Pa (upper limit of 150 Pa), 4-30 Pa for the shear strength of fine surface materials and boulders, and 30-150 Pa for the compressive strength of overhangs (upper limit of 1500 Pa). The strength-to-gravity ratio is similar for 67P and weak rocks on Earth. As a result of the low compressive strength, the interior of the nucleus may have been compressed sufficiently to initiate diagenesis, which could have contributed to the formation of layers. Our value for the tensile strength is comparable to that of dust aggregates formed by gravitational instability and tends to favor a formation of comets by the accrection of pebbles at low velocities.

We report on the first major temporal morphological changes observed on the surface of the nucleus of comet 67P/Churyumov-Gerasimenko, in the smooth terrains of the Imhotep region. We use images of the OSIRIS cameras onboard Rosetta to follow the temporal changes from 24 May 2015 to 11 July 2015. The morphological changes observed on the surface are visible in the form of roundish features, which are growing in size from a given location in a preferential direction, at a rate of 5.6 - 8.1$\times$10$^{-5}$ m s$^{-1}$ during the observational period. The location where changes started and the contours of the expanding features are bluer than the surroundings, suggesting the presence of ices (H$_2$O and/or CO$_2$) exposed on the surface. However, sublimation of ices alone is not sufficient to explain the observed expanding features. No significant variations in the dust activity pattern are observed during the period of changes.