This is the third paper of a series where we study the stellar population gradients (SP; ages, metallicities, $\alpha$-element abundance ratios and stellar initial mass functions) of early type galaxies (ETGs) at $z\le 0.08$ from the MaNGA-DR15 survey. In this work we focus on the S0 population and quantify how the SP varies across the population as well as with galactocentric distance. We do this by measuring Lick indices and comparing them to stellar population synthesis models. This requires spectra with high signal-to-noise which we achieve by stacking in bins of luminosity (L$_r$) and central velocity dispersion ($\sigma_0$). We find that: 1) There is a bimodality in the S0 population: S0s more massive than $3\times 10^{10}M_\odot$ show stronger velocity dispersion and age gradients (age and $\sigma_r$ decrease outwards) but little or no metallicity gradient, while the less massive ones present relatively flat age and velocity dispersion profiles, but a significant metallicity gradient (i.e. [M/H] decreases outwards). Above $2\times10^{11}M_\odot$ the number of S0s drops sharply. These two mass scales are also where global scaling relations of ETGs change slope. 2) S0s have steeper velocity dispersion profiles than fast rotating elliptical galaxies (E-FRs) of the same luminosity and velocity dispersion. The kinematic profiles and stellar population gradients of E-FRs are both more similar to those of slow rotating ellipticals (E-SRs) than to S0s, suggesting that E-FRs are not simply S0s viewed face-on. 3) At fixed $\sigma_0$, more luminous S0s and E-FRs are younger, more metal rich and less $\alpha$-enhanced. Evidently for these galaxies, the usual statement that 'massive galaxies are older' is not true if $\sigma_0$ is held fixed.

Early-type galaxies -- slow and fast rotating ellipticals (E-SRs and E-FRs) and S0s/lenticulars -- define a Fundamental Plane (FP) in the space of half-light radius $R_e$, enclosed surface brightness $I_e$ and velocity dispersion $\sigma_e$. Since $I_e$ and $\sigma_e$ are distance-independent measurements, the thickness of the FP is often expressed in terms of the accuracy with which $I_e$ and $\sigma_e$ can be used to estimate sizes $R_e$. We show that: 1) The thickness of the FP depends strongly on morphology. If the sample only includes E-SRs, then the observed scatter in $R_e$ is $\sim 16\%$, of which only $\sim 9\%$ is intrinsic. Removing galaxies with $M_*<10^{11}M_\odot$ further reduces the observed scatter to $\sim 13\%$ ($\sim 4\%$ intrinsic). The observed scatter increases to the $\sim 25\%$ usually quoted in the literature if E-FRs and S0s are added. If the FP is defined using the eigenvectors of the covariance matrix of the observables, then the E-SRs again define an exceptionally thin FP, with intrinsic scatter of only $5\%$ orthogonal to the plane. 2) The structure within the FP is most easily understood as arising from the fact that $I_e$ and $\sigma_e$ are nearly independent, whereas the $R_e-I_e$ and $R_e-\sigma_e$ correlations are nearly equal and opposite. 3) If the coefficients of the FP differ from those associated with the virial theorem the plane is said to be `tilted'. If we multiply $I_e$ by the global stellar mass-to-light ratio $M_*/L$ and we account for non-homology across the population by using S\'ersic photometry, then the resulting stellar mass FP is less tilted. Accounting self-consistently for $M_*/L$ gradients will change the tilt. The tilt we currently see suggests that the efficiency of turning baryons into stars increases and/or the dark matter fraction decreases as stellar surface brightness increases.

We propose an analytical approach to study non-Markov random walks by employing an exact enumeration method. Using the method, we derive an exact expansion for the first-passage time (FPT) distribution for any continuous, differentiable non-Markov random walk with Gaussian or non-Gaussian multivariate distribution. As an example, we study the FPT distribution of a fractional Brownian motion with a Hurst exponent $H\in(1/2,1)$ that describes numerous non-Markov stochastic phenomena in physics, biology and geology, and for which the limit $H=1/2$ represents a Markov process.

The excursion set approach is a framework for estimating how the number density of nonlinear structures in the cosmic web depends on the expansion history of the universe and the nature of gravity. A key part of the approach is the estimation of the first crossing distribution of a suitably chosen barrier by random walks having correlated steps: The shape of the barrier is determined by the physics of nonlinear collapse, and the correlations between steps by the nature of the initial density fluctuation field. We describe analytic and numerical methods for calculating such first up-crossing distributions. While the exact solution can be written formally as an infinite series, we show how to approximate it efficiently using the Stratonovich approximation. We demonstrate its accuracy using Monte-Carlo realizations of the walks, which we generate using a novel Cholesky-decomposition based algorithm, which is significantly faster than the algorithm that is currently in the literature.

Excursion set theory (EST) is an analytical framework to study the large-scale structure of the Universe. EST introduces a procedure to calculate the number density of structures by relating the cosmological linear perturbation theory to the nonlinear structures in late time. In this work, we introduce a novel approach to reformulate the EST in matrix formalism. We propose that the matrix representation of EST will facilitate the calculations in this framework. The method is to discretize the two-dimensional plane of variance and density contrast of EST, where the trajectories for each point in the Universe lived there. The probability of having a density contrast in a chosen variance is represented by a probability ket. Naturally, the concept of the transition matrix pops up to define the trajectories. We also define the probability transition rate which is used to obtain the first up-crossing of trajectories and the number count of the structures. In this work we show that the discretization let us study the non-Markov processes by forcing them to look like a Wiener process. Also we discuss that the zero drift processes with Gaussian and also non-Gaussian initial conditions can be studied by this formalism. The continuous limit of the formalism is discussed, and the known Fokker-Planck dispersion equation is recovered. Finally we show that the probability of the most massive progenitors can be extracted in this framework.