  We study lithium depletion in low-mass and solar-like stars as a function of
time, using a new diffusion coefficient describing extra-mixing taking place at
the bottom of a convective envelope. This new form is motivated by
multi-dimensional fully compressible, time implicit hydrodynamic simulations
performed with the MUSIC code. Intermittent convective mixing at the convective
boundary in a star can be modeled using extreme value theory, a statistical
analysis frequently used for finance, meteorology, and environmental science.
In this letter, we implement this statistical diffusion coefficient in a
one-dimensional stellar evolution code, using parameters calibrated from
multi-dimensional hydrodynamic simulations of a young low-mass star. We propose
a new scenario that can explain observations of the surface abundance of
lithium in the Sun and in clusters covering a wide range of ages, from $\sim$
50 Myr to $\sim$ 4 Gyr. Because it relies on our physical model of convective
penetration, this scenario has a limited number of assumptions. It can explain
the observed trend between rotation and depletion, based on a single additional
assumption, namely that rotation affects the mixing efficiency at the
convective boundary. We suggest the existence of a threshold in stellar
rotation rate above which rotation strongly prevents the vertical penetration
of plumes and below which rotation has small effects. In addition to providing
a possible explanation for the long standing problem of lithium depletion in
pre-main sequence and main sequence stars, the strength of our scenario is that
its basic assumptions can be tested by future hydrodynamic simulations.
