We study the formation and evolution of a convective layer when a stably-stratified fluid with a composition gradient is cooled from above. We perform a series of 2D simulations using the Bousinessq approximation with Prandtl number ranging from Pr = 0.1 to 7, extending previous work on salty water to low Pr. We show that the evolution of the convection zone is well-described by an entrainment prescription in which a fixed fraction of the kinetic energy of convective motions is used to mix fluid at the interface with the stable layer. We measure the entrainment efficiency and find that it grows with decreasing Prandtl number or increased applied heat flux. The kinetic energy flux that determines the entrainment rate is a small fraction of the total convective luminosity. In this time-dependent situation, the density ratio at the interface is driven to a narrow range that depends on the value of Pr, and with low enough values that advection dominates the interfacial transport. We characterize the interfacial flux ratio and how it depends on the interface stability. We present an analytic model that accounts for the growth of the convective layer with two parameters, the entrainment efficiency and the interfacial heat transport, both of which can be measure from the simulations.

The primary advantage of moderately superheated bubble chamber detectors is their simultaneous sensitivity to nuclear recoils from WIMP dark matter and insensitivity to electron recoil backgrounds. A comprehensive analysis of PICO gamma calibration data demonstrates for the first time that electron recoils in C$_3$F$_8$ scale in accordance with a new nucleation mechanism, rather than one driven by a hot-spike as previously supposed. Using this semi-empirical model, bubble chamber nucleation thresholds may be tuned to be sensitive to lower energy nuclear recoils while maintaining excellent electron recoil rejection. The PICO-40L detector will exploit this model to achieve thermodynamic thresholds as low as 2.8 keV while being dominated by single-scatter events from coherent elastic neutrino-nucleus scattering of solar neutrinos. In one year of operation, PICO-40L can improve existing leading limits from PICO on spin-dependent WIMP-proton coupling by nearly an order of magnitude for WIMP masses greater than 3 GeV c$^{-2}$ and will have the ability to surpass all existing non-xenon bounds on spin-independent WIMP-nucleon coupling for WIMP masses from 3 to 40 GeV c$^{-2}$.

Next generation surveys will be capable of determining cosmological parameters beyond percent level. To match this precision, theoretical descriptions should look beyond the linear perturbations to approximate the observables in large scale structure. A quantity of interest is the Number density of galaxies detected by our instruments. This has been focus of interest recently, and several efforts have been made to explain relativistic effects theoretically, thereby testing the full theory. However, the results at nonlinear level from previous works are in disagreement. We present a new and independent approach to computing the relativistic galaxy number counts to second order in cosmological perturbation theory. We derive analytical expressions for the full second order relativistic observed redshift, for the angular diameter distance and for the volume spanned by a survey. Finally, we compare our results with previous works which compute the general distance-redshift relation, finding that our result is in agreement at linear order.

Glitches are rare spin-up events that punctuate the smooth slow-down of the rotation of pulsars. For the Vela pulsar and PSR J0537-6910, the glitch sizes and the times between consecutive events have clear preferred scales, contrary to the handful of other pulsars with enough glitches for such a study. Moreover, PSR J0537-6910 is the only pulsar showing a strong positive correlation between the size of each glitch and the waiting time until the following one. We attempt to understand this behavior through a detailed study of the distributions and correlations of glitch properties for the eight pulsars with at least ten detected glitches. We model the distributions of glitch sizes and times between consecutive glitches for this sample and perform Monte Carlo simulations to explore two hypotheses that could explain why the correlation is so much weaker in other pulsars than in PSR J0537-6910. We confirm the above results for the Vela pulsar and PSR J0537-6910, and verify that the latter is the only pulsar with a strong correlation between glitch size and waiting time to the following glitch. For the remaining six pulsars, the waiting time distributions are best fitted by exponentials, and the size distributions either by power laws, exponentials, or log-normal functions. Some pulsars in the sample yield significant Pearson and Spearman coefficients ($r_p$ and $r_s$) for the aforementioned correlation. Moreover, for all except the Crab, both coefficients are positive. Our simulations show that the weaker correlations in pulsars other than PSR J0537$-$6910 cannot be due to missing glitches too small to be detected. We also tested the hypothesis that each pulsar may have two kinds of glitches, namely large, correlated ones and small, uncorrelated ones. We found for all pulsars that their glitches are not consistent with a pure uncorrelated population either.

We study a class of almost scale-invariant modified gravity theories, using a particular form of $f(R, G) = \alpha R^2 + \beta G \log G$ where $R$ and $G$ are the Ricci and Gauss-Bonnet scalars, respectively and $\alpha$, $\beta$ are arbitrary constants. We derive the Einstein-like field equations to first order in cosmological perturbation theory in longitudinal gauge.

We present a statistical study of the glitch population and the behaviour of the glitch activity across the known population of neutron stars. An unbiased glitch database was put together based on systematic searches of radio timing data of 898 rotation-powered pulsars obtained with the Jodrell Bank and Parkes observatories. Glitches identified in similar searches of 5 magnetars were also included. The database contains 384 glitches found in the rotation of 141 of these neutron stars. We confirm that the glitch size distribution is at least bimodal, with one sharp peak at approximately $20\, \rm{\mu\,Hz}$, which we call large glitches, and a broader distribution of smaller glitches. We also explored how the glitch activity $\dot{\nu}_{\rm{g}}$, defined as the mean frequency increment per unit of time due to glitches, correlates with the spin frequency $\nu$, spin-down rate $|\dot{\nu}|$, and various combinations of these, such as energy loss rate, magnetic field, and spin-down age. It is found that the activity is insensitive to the magnetic field and that it correlates strongly with the energy loss rate, though magnetars deviate from the trend defined by the rotation-powered pulsars. However, we find that a constant ratio $\dot\nu_{\rm{g}}/|\dot\nu| = 0.010 \pm 0.001$ is consistent with the behaviour of all rotation-powered pulsars and magnetars. This relation is dominated by large glitches, which occur at a rate directly proportional to $|\dot{\nu}|$. The only exception are the rotation-powered pulsars with the highest values of $|\dot{\nu}|$, such as the Crab pulsar and PSR B0540$-$69, which exhibit a much smaller glitch activity, intrinsically different from each other and from the rest of the population. The activity due to small glitches also shows an increasing trend with $|\dot\nu|$, but this relation is biased by selection effects.