Researcher profile

Francesco Marzari

Francesco Marzari contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
19works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

19 published item(s)

preprint2022arXiv

Dust resurgence in protoplanetary disks due to planetesimal-planet interactions

Observational data on the dust content of circumstellar disks show that the median dust content in disks around pre-main sequence stars in nearby star forming regions seem to increase from about 1 Myr to about 2 Myr, and then decline with time. This behaviour challenges the models where the small dust grains steadily decline by accumulating into larger bodies and drifting inwards on a short timescale (less than about 1 Myr). In this Letter we explore the possibility to reconcile this discrepancy in the framework of a model where the early formation of planets dynamically stirs the nearby planetesimals and causes high energy impacts between them, resulting in the production of second-generation dust. We show that the observed dust evolution can be naturally explained by this process within a suite of representative disk-planet architectures.

preprint2020arXiv

Evolution of an asteroid family under YORP, Yarkovsky and collisions

Any population of asteroids, like asteroid families, will disperse in semi-major axis due to the Yarkovsky effect. The amount of drift is modulated by the asteroid spin state evolution which determines the balance between the diurnal and seasonal Yarkovsky force. The asteroid's spin state is, in turn, controlled in part by the YORP effect. The otherwise smooth evolution of an asteroid can be abruptly altered by collisions, which can cause impulsive changes in the spin state and can move the asteroid onto a different YORP track. In addition, collisions may also alter the YORP parameters by changing the superficial features and overall shape of the asteroid. Thus, the coupling between YORP and Yarkovsky is also strongly affected by the impact history of each body. To investigate this coupling we developed a statistical code modeling the time evolution of semi--major axis under YORP-Yarkovsky coupling. It includes the contributions of NYORP (normal YORP), TYORP (tangential YORP) and collisions whose effects are deterministically calculated and not added in a statistical way. We find that both collisions and TYORP increase the dispersion of a family in semi-major axis by making the spin axis evolution less smooth and regular. We show that the evolution of a family's structure with time is complex and collisions randomize the YORP evolution. In our test families we do not observe the formation of a 'YORP-eye' in the semi-major axis vs. diameter distribution, even after a long period of time. If present, the 'YORP-eye' might be a relic of an initial ejection velocity pattern of the collisional fragments.

preprint2020arXiv

Planets in Binaries: Formation and Dynamical Evolution

Binary systems are very common among field stars. While this relatively small number of planets in binaries is probably partly due to strong observational biases, there is, however, statistical evidence that planets are indeed less frequent in binaries with separations smaller than 100 au, strongly suggesting that the presence of a close in companion star has an adverse effect on planet formation. It is indeed possible for the gravitational pull of the second star to affect all the different stages of planet formation, from proto-planetary disk formation to dust accumulation into planetesimals, to the accretion of these planetesimals into large planetary embryos and, eventually, the final growth of these embryos into planets. For the crucial planetesimal accretion phase, the complex coupling between dynamical perturbations from the binary and friction due to gas in the protoplanetary disk suggests that planetesimal accretion might be hampered due to increased, accretion hostile impact velocities. Likewise, the interplay between the binary secular perturbations and mean motion resonances lead to unstable regions, where not only planet formation is inhibited, but where a massive body would be ejected from the system on a hyperbolic orbit. The amplitude of these two main effects is different for S and P type planets, so that a comparison between the two populations might outline the influence of the companion star on the planet formation process. Unfortunately, at present the two populations (circumstellar or circumbinary) are not known equally well and different biases and uncertainties prevent a quantitative comparison. We also highlight the long term dynamical evolution of both S and P type systems and focus on how these different evolutions influence the final architecture of planetary systems in binaries.

preprint2020arXiv

Secular evolution of close in planets: the effects of general relativity

Pairs of planets in a system may end up close to their host star on eccentric orbits as a consequence of planet planet scattering, Kozai or secular migration. In this scenario, general relativity and secular perturbations have comparable timescales and may interfere with each other with relevant effects on the eccentricity and pericenter evolution of the two planets. We explore, both analytically and via numerical integration, how the secular evolution is changed by general relativity for a wide range of different initial conditions. We find that when the faster secular frequency approaches the general relativity precession rate, which tipically occurs when the outer planet moves away from the inner one, it relaxes to it and a significant damping of the proper eccentricity of the inner planet occurs. The proper eccentricity of the outer planet is reduced as well due to the changes in the secular interaction of the bodies. The lowering of the peak eccentricities of the two planets during their secular evolution has important implications on their stability. A significant number of two planet systems, otherwise chaotic because of the mutual secular perturbations, are found stable when general relativity is included.

preprint2020arXiv

The GAPS Programme at TNG XXVIII -- A pair of hot-Neptunes orbiting the young star TOI-942

Both young stars and multi-planet systems are primary objects that allow us to study, understand and constrain planetary formation and evolution theories. We validate the physical nature of two Neptune-type planets transiting TOI-942 (TYC 5909-319-1), a previously unacknowledged young star (50+30-20 Myr) observed by the TESS space mission in Sector 5. Thanks to a comprehensive stellar characterization, TESS light curve modelling and precise radial-velocity measurements, we validated the planetary nature of the TESS candidate and detect an additional transiting planet in the system on a larger orbit. From photometric and spectroscopic observations we performed an exhaustive stellar characterization and derived the main stellar parameters. TOI-942 is a relatively active K2.5V star (logR'hk = -4.17+-0.01) with rotation period Prot = 3.39+-0.01 days, a projected rotation velocity vsini=13.8+-0.5 km/s and a radius of ~0.9 Rsun. We found that the inner planet, TOI-942b, has an orbital period Pb=4.3263+-0.0011 days, a radius Rb=4.242-0.313+0.376 Rearth and a mass upper limit of 16 Mearth at 1-sigma confidence level. The outer planet, TOI-942c, has an orbital period Pc=10.1605-0.0053+0.0056 days, a radius Rc=4.793-0.351+0.410 Rearth and a mass upper limit of 37 Mearth at 1-sigma confidence level.

preprint2016arXiv

Matching asteroid population characteristics with a model constructed from the YORP-induced rotational fission hypothesis

From the results of a comprehensive asteroid population evolution model, we conclude that the YORP-induced rotational fission hypothesis can be consistent with the observed population statistics of small asteroids in the main belt including binaries and contact binaries. The foundation of this model is the asteroid rotation model of Marzari et al. (2011), which incorporates both the YORP effect and collisional evolution. This work adds to that model the rotational fission hypothesis and the binary evolution model of Jacobson & Scheeres (2011). The asteroid population evolution model is highly constrained by these and other previous works, and therefore it has only two significant free parameters: the ratio of low to high mass ratio binaries formed after rotational fission events and the mean strength of the binary YORP (BYORP) effect. We successfully reproduce characteristic statistics of the small asteroid population: the binary fraction, the fast binary fraction, steady-state mass ratio fraction and the contact binary fraction. We find that in order for the model to best match observations, rotational fission produces high mass ratio ($> 0.2$) binary components with four to eight times the frequency as low mass ratio ($< 0.2$) components, where the mass ratio is the mass of the secondary component divided by the mass of the primary component. This is consistent with post-rotational fission binary system mass ratio being drawn from either a flat or a positive and shallow distribution, since the high mass ratio bin is four times the size of the low mass ratio bin; this is in contrast to the observed steady-state binary mass ratio, which has a negative and steep distribution. This can be understood in the context of the BYORP-tidal equilibrium hypothesis, which predicts that low mass ratio binaries survive for a significantly longer period of time than high mass ratio systems.

preprint2015arXiv

Decoupling of a giant planet from its disk in an inclined binary system

We explore the dynamical evolution of a planet embedded in a disk surrounding a star part of a binary system where the orbital plane of the binary is significantly tilted respect to the initial disk plane. Our aim is to test whether the planet remains within the disk and continues to migrate towards the star in a Type I/II mode in spite of the secular perturbations of the companion star. This would explain observed exoplanets with significant inclination respect to the equatorial plane of their host star. We have used two different SPH codes, vine and phantom, to model the evolution of a system star+disk+planet and companion star with time. After an initial coupled evolution, the inclination of the disk and that of the planet begin to differ significantly. The period of oscillation of the disk inclination, respect to the initial plane, is shorter than that of the planet which evolves independently after about 10^4 yr following a perturbed N-body behavior. However, the planet keeps migrating towards the star because during its orbital motion it crosses the disk plane and the friction with the gas causes angular momentum loss. Disk and planet in a significantly inclined binary system are not dynamically coupled for small binary separations but evolve almost independently. The planet abandons the disk and, due to the onset of a significant mutual inclination, it interacts with the gas only when its orbit intersects the disk plane. The drift of the planet towards the star is not due to type I/II with the planet embedded in the disk but to the friction with the gas during the disk crossing.

preprint2015arXiv

Detection of Sharp Symmetric Features in the Circumbinary Disk Around AK Sco

The Search for Planets Orbiting Two Stars (SPOTS) survey aims to study the formation and distribution of planets in binary systems by detecting and characterizing circumbinary planets and their formation environments through direct imaging. With the SPHERE Extreme Adaptive Optics instrument, a good contrast can be achieved even at small (<300 mas) separations from bright stars, which enables studies of planets and disks in a separation range that was previously inaccessible. Here, we report the discovery of resolved scattered light emission from the circumbinary disk around the well-studied young double star AK Sco, at projected separations in the ~13--40 AU range. The sharp morphology of the imaged feature is surprising, given the smooth appearance of the disk in its spectral energy distribution. We show that the observed morphology can be represented either as a highly eccentric ring around AK Sco, or as two separate spiral arms in the disk, wound in opposite directions. The relative merits of these interpretations are discussed, as well as whether these features may have been caused by one or several circumbinary planets interacting with the disk.

preprint2015arXiv

The Complex History of Trojan Asteroids

The Trojan asteroids provide a unique perspective on the history of Solar System. As a large population of small bodies, they record important gravitational interactions and dynamical evolution of the Solar System. In the past decade, significant advances have been made in understanding physical properties, and there has been a revolution in thinking about the origin of Trojans. The ice and organics generally presumed to be a significant part of Trojan compositions have yet to be detected directly, though low density of the binary system Patroclus (and possibly low density of the binary/moonlet system Hektor) is consistent with an interior ice component. By contrast, fine-grained silicates that appear to be similar to cometary silicates in composition have been detected, and a color bimodality may indicate distinct compositional groups among the Trojans. Whereas Trojans had traditionally been thought to have formed near 5 AU, a new paradigm has developed in which the Trojans formed in the proto-Kuiper Belt, and they were scattered inward and captured in the Trojan swarms as a result of resonant interactions of the giant planets. Whereas the orbital and population distributions of current Trojans are consistent with this origin scenario, there are significant differences between current physical properties of Trojans and those of Kuiper Belt objects. These differences may be indicative of surface modification due to the inward migration of objects that became the Trojans, but understanding of appropriate modification mechanisms is poor and would benefit from additional laboratory studies. Many open questions remain, and the future promises significant strides in our understanding of Trojans. The time is ripe for a spacecraft mission to the Trojans, to turn these objects into geologic worlds that can be studied in detail to unravel their complex history.

preprint2014arXiv

Effect of rotational disruption on the size-frequency distribution of the Main Belt asteroid population

The size distribution of small asteroids in the Main Belt is assumed to be determined by an equilibrium between the creation of new bodies out of the impact debris of larger asteroids and the destruction of small asteroids by collisions with smaller projectiles. However, for a diameter less than 6 km we find that YORP-induced rotational disruption significantly contributes to the erosion even exceeding the effects of collisional fragmentation. Including this additional grinding mechanism in a collision evolution model for the asteroid belt, we generate size-frequency distributions from either an accretional (Weidenschilling, 2011) or an &#34;Asteroids were born big&#34; (Morbidelli, 2009) initial size-frequency distribution that are consistent with observations reported in Gladman et al. (2009). Rotational disruption is a new mechanism that must be included in all future collisional evolution models of asteroids.

preprint2014arXiv

Effects of stellar flybys on planetary systems: 3D modeling of the circumstellar disks damping effects

Stellar flybys in star clusters are suspected to affect the orbital architecture of planetary systems causing eccentricity excitation and orbital misalignment between the planet orbit and the equatorial plane of the star. We explore whether the impulsive changes in the orbital elements of planets, caused by an hyperbolic stellar flyby, can be fully damped by the circumstellar disk surrounding the star. The time required to disperse stellar clusters is in fact comparable to circumstellar disk&#39;s lifetime. We have modelled in 3D a system made of a solar type star surrounded by a low density disk with a giant planet embedded in it approached on a hyperbolic encounter trajectory by a second star, of similar mass and with its own disk. We focus on extreme configurations where a very deep stellar flyby perturbs a Jovian planet on an external orbit. This allows to test in full the ability of the disk to erase the effects of the stellar encounter. We find that the amount of mass lost by the disk during the stellar flyby is less than in 2D models where a single disk was considered due to the mass exchange between the two disks at the encounter. The damping in eccentricity is slightly faster than in 2D models and it occurs on timescales of the order of a few kyr. The only trace of the flyby left in the planet system, after about 10^4 yr, is a small misalignment, lower than 9 degrees, between the star equatorial plane and the planet orbit. In a realistic model based on 3D simulations of star--planet--disk interactions, we find that stellar flybys cannot excite significant eccentricities and inclinations of planets in stellar clusters. The circumstellar disks hosting the planets damp on a short timescale all the step changes in the two orbital parameters produced during any stellar encounter. All records of past encounters are erased.

preprint2014arXiv

SPOTS: The Search for Planets Orbiting Two Stars. I. Survey description and first observations

Direct imaging surveys for exoplanets commonly exclude binary stars from their target lists, leaving a large part of the overall planet demography unexplored. To address this gap in our understanding of planet formation and evolution, we have launched the first direct imaging survey dedicated to circumbinary planets: SPOTS, the Search for Planets Orbiting Two Stars. In this paper, we discuss the theoretical context, scientific merit, and technical feasibility of such observations, describe the target sample and observational strategy of our survey, and report on the first results from our pilot survey of 26 targets with the VLT NaCo facility. While we have not found any confirmed substellar companions to date, a number of promising candidate companions remain to be tested for common proper motion in upcoming follow-up observations. We also report on the astrometry of the three resolved binaries in our target sample. This pilot survey constitutes a successful proof of concept for our survey strategy and paves the way for a second stage of exploratory observations with VLT SPHERE.

preprint2014arXiv

TRADES: A new software to derive orbital parameters from observed transit times and radial velocities. Revisiting Kepler-11 and Kepler-9

Aims. With the purpose of determining the orbital parameters of exoplanetary systems from observational data, we have developed a software, named TRADES (TRAnsits and Dynamics of Exoplanetary Systems), to simultaneously fit observed radial velocities and transit times data. Methods. We implemented a dynamical simulator for N-body systems, which also fits the available data during the orbital integration and determines the best combination of the orbital parameters using grid search, $χ^2$ minimization, genetic algorithms, particle swarm optimization, and bootstrap analysis. Results. To validate TRADES, we tested the code on a synthetic three-body system and on two real systems discovered by the Kepler mission: Kepler-9 and Kepler-11. These systems are good benchmarks to test multiple exoplanet systems showing transit time variations (TTVs) due to the gravitational interaction among planets. We have found that orbital parameters of Kepler-11 planets agree well with the values proposed in the discovery paper and with a a recent work from the same authors. We analyzed the first three quarters of Kepler-9 system and found parameters in partial agreement with discovery paper. Analyzing transit times (T0s) covering 12 quarters of Kepler data, that we have found a new best-fit solution. This solution outputs masses that are about 55% of the values proposed in the discovery paper; this leads to a reduced semi-amplitude of the radial velocities of about 12.80 m/s.

preprint2013arXiv

Influence of the circumbinary disk gravity on planetesimal accumulation in the Kepler 16 system

Recent observations from NASA&#39;s Kepler mission detected the first planets in circumbinary orbits. The question we try to answer is where these planets formed in the circumbinary disk and how far inside they migrated to reach their present location. We investigate the first and more delicate phase of planet formation when planetesimals accumulate to form planetary embryos. We use the hydrodynamical code FARGO to study the evolution of the disk and of a test population of planetesimals embedded in it. With this hybrid hydrodynamical--N--body code we can properly account for the gas drag force on the planetesimals and for the gravitational force of the disk on them. The numerical simulations show that the gravity of the eccentric disk on the planetesimal swarm excites their eccentricities to values much larger than those induced by the binary perturbations only within 10 AU from the stars. Moreover, the disk gravity prevents a full alignment of the planetesimal pericenters. Both these effects lead to large impact velocities, beyond the critical value for erosion. Planetesimals accumulation in circumbinary disks appears to be prevented close to the stellar pair by the gravitational perturbations of the circumbinary disk. The observed planets possibly formed in the outer regions of the disk and then migrated inside by tidal interaction with the disk.

preprint2013arXiv

Three-dimensional modeling of radiative disks in binaries

Circumstellar disks in binaries are perturbed by the companion gravity causing significant alterations of the disk morphology. Spiral waves due to the companion tidal force also develop in the vertical direction and affect the disk temperature profile. These effects may significantly influence the process of planet formation. We perform 3D numerical simulations of disks in binaries with different initial dynamical configurations and physical parameters. Our goal is to investigate their evolution and their propensity to grow planets. We use an improved version of the SPH code VINE modified to better account for momentum and energy conservation. The energy equation includes a flux--limited radiative transfer algorithm and the disk cooling is obtained via &#34;boundary particles&#34;. We model a system made of star/disk + star/disk where the secondary star (and relative disk) is less massive than the primary. The numerical simulations performed for different values of binary separation and disk density show that the disk morphology is substantially affected by the companion perturbations. Trailing spiral shock waves develop when the stars approach their pericenter. Strong hydraulic jumps occur at the shock front creating breaking waves and a consistent mass stream between the two disks, significantly heating them. The high gas temperature may prevent the ice condensation by moving outward the &#34;snow line&#34;. The hydraulic jumps may slow down or even halt the dust coagulation process. At apocenter these perturbations are reduced and the disks are cooled down and less eccentric. The strength of the hydraulic jumps, disk heating, and mass exchange depends on the binary separation, and for larger semi-major axes, the tidal spiral pattern is substantially reduced.

preprint2012arXiv

Outward Migration of Jupiter and Saturn in Evolved Gaseous Disks

The outward migration of a pair of resonant-orbit planets, driven by tidal interactions with a gas-dominated disk, is studied in the context of evolved solar nebula models. The planets&#39; masses, M1 and M2, correspond to those of Jupiter and Saturn. Hydrodynamical calculations in two and three dimensions are used to quantify the migration rates and analyze the conditions under which the outward migration mechanism may operate. The planets are taken to be fully formed after 1e6 and before 3e6 years. The orbital evolution of the planets in an evolving disk is then calculated until the disk&#39;s gas is completely dissipated. Orbital locking in the 3:2 mean motion resonance may lead to outward migration under appropriate conditions of disk viscosity and temperature. However, resonance locking does not necessarily result in outward migration. This is the case, for example, if convergent migration leads to locking in the 2:1 mean motion resonance, as post-formation disk conditions seem to suggest. Accretion of gas on the planets may deactivate the outward migration mechanism by raising the mass ratio M2/M1 and/or by reducing the accretion rate toward the star, hence depleting the inner disk. For migrating planets locked in the 3:2 mean motion resonance, there are stalling radii that depend on disk viscosity and on stellar irradiation, when it determines the disk&#39;s thermal balance. Planets locked in the 3:2 orbital resonance that start moving outward from within 1-2 AU may reach beyond 5 AU only under favorable conditions. However, within the explored space of disk parameters, only a small fraction - less than a few per cent - of the models predict that the interior planet reaches beyond 4 AU.

preprint2011arXiv

On how the optical depth tunes the effects of ISM neutral atom flow on debris disks

The flux of ISM neutral atoms surrounding stars and their environment affects the motion of dust particles in debris disks, causing a significant dynamical evolution. Large values of eccentricity and inclination can be excited and strong correlations settle in among the orbital angles. This dynamical behaviour, in particular for bound dust grains, can potentially cause significant asymmetries in dusty disks around solar type stars which might be detected by observations. However, the amount of orbital changes due to this non--gravitational perturbation is strongly limited by the collisional lifetime of dust particles. We show that for large values of the disk&#39;s optical depth the influence of ISM flow on the disk shape is almost negligible because the grains are collisionally destroyed before they can accumulate enough orbital changes due to the ISM perturbations. On the other hand, for values smaller than $10^{-3}$, peculiar asymmetric patterns appear in the density profile of the disk when we consider 1-10 mum grains, just above the blow-out threshold. The extent and relevance of these asymmetries grow for lower values of the optical depth. An additional sink mechanism, which may prevent the formation of large clumps and warping in the disks is related to the fast inward migration due to the drag component of the forces. When a significant eccentricity is pumped up by the ISM perturbations, the drag forces (Poynting-Robertson and in particular ISM drag) drive the disk particles on fast migrating tracks leading them into the star on a short timescale. It is then expected that disks with small optical depth expand inside the parent body ring all the way towards the star while disks with large optical depth would not significantly extend inside.

preprint2010arXiv

Debris discs in binaries: a numerical study

Debris disc analysis and modelling provide crucial information about the structure and the processes at play in extrasolar planetary systems. In binary systems, this issue is more complex because the disc should in addition respond to the companion star&#39;s perturbations. We explore the dynamical evolution of a collisionally active debris disc for different initial parent body populations, diverse binary configurations and optical depths. We focus on the radial extent and size distribution of the disc at a stationary state. We numerically follow the evolution of $10^{5}$ massless small grains, initially produced from a circumprimary disc of parent bodies following a size distribution in $dN \propto s^{-3.5}$ds . Grains are submitted to both stars&#39; gravity as well as radiation pressure. In addition, particles are assigned an empirically derived collisional lifetime. For all the binary configurations the disc extends far beyond the critical semimajor axis $a_crit$ for orbital stability. This is due to the steady production of small grains, placed on eccentric orbits reaching beyond $a_crit$ by radiation pressure. The amount of matter beyond acrit depends on the balance between collisional production and dynamical removal rates: it increases for more massive discs as well as for eccentric binaries. Another important effect is that, in the dynamically stable region, the disc is depleted from its smallest grains. Both results could lead to observable signatures. We have shown that a companion star can never fully truncate a collisionally active disc. For eccentric companions, grains in the unstable regions can significantly contribute to the thermal emission in the mid-IR. Discs with sharp outer edges, especially bright ones such as HR4796A, are probably shaped by other mechanisms.

preprint2010arXiv

Effects of Interplanetary Dust on the LISA drag-free Constellation

The analysis of non-radiative sources of static or time-dependent gravitational fields in the Solar System is crucial to accurately estimate the free-fall orbits of the LISA space mission. In particular, we take into account the gravitational effects of Interplanetary Dust (ID) on the spacecraft trajectories. The perturbing gravitational field has been calculated for some ID density distributions that fit the observed zodiacal light. Then we integrated the Gauss planetary equations to get the deviations from the LISA keplerian orbits around the Sun. This analysis can be eventually extended to Local Dark Matter (LDM), as gravitational fields are expected to be similar for ID and LDM distributions. Under some strong assumptions on the displacement noise at very low frequency, the Doppler data collected during the whole LISA mission could provide upper limits on ID and LDM densities.