We present a public code to generate random fields with an arbitrary probability distribution function (PDF) and an arbitrary correlation function. The algorithm is cosmology-independent, applicable to any stationary stochastic process over a three dimensional grid. We implement it in the case of the matter density field, showing its benefits over the lognormal approximation, which is often used in cosmology for generation of mock catalogues. We find that the covariance of the power spectrum from the new fast realizations is more accurate than that from a lognormal model. As a proof of concept, we also apply the new simulation scheme to the divergence of the Lagrangian displacement field. We find that information from the correlation function and the PDF of the displacement-divergence provides modest improvement over other standard analytical techniques to describe the particle field in the simulation. This suggests that further progress in this direction should come from multi-scale or non-local properties of the initial matter distribution.

The quest for understanding the late-time acceleration is haunted by an immense freedom in the analysis of dynamical models for dark energy in extended parameter spaces. Often-times having no prior knowledge at our disposal, arbitrary choices are implemented to reduce the degeneracies between parameters. We also encounter this issue in the case of quintessence fields, where a scalar degree of freedom drives the late-time acceleration. In this study, we implement a more physical prescription, the \textit{flow} condition, to fine-tune the quintessence evolution for several field potentials. We find that this prescription agrees well with the most recent catalogue of data, namely supernovae type Ia, baryon acoustic oscillations, cosmic clocks and distance to last scattering surface, and it enables us to infer the initial conditions for the field, both potential and cosmological parameters. At $2\sigma$ we find stricter bounds on the potential parameters $f/m_{pl}>0.26$ and $n<0.15$ for the PNGB and IPL potentials, respectively, while constraints on cosmological parameters remain extremely consistent across all assumed potentials. By implementing information criteria to assess their ability to fit the data, we do not find any evidence against thawing models, which in fact are statistically equivalent to $\Lambda$CDM, and the freezing ones are moderately disfavoured. Through our analysis we place upper bounds on the slope of quintessence potentials, consequently revealing a strong tension with the recently proposed swampland criterion, finding the 2$\sigma$ upper bound of $\lambda \sim 0.31$ for the exponential potential.

Solar UV variability is extremely relevant for the stratospheric ozone. It has an impact on Earth's atmospheric structure and dynamics through radiative heating and ozone photochemistry. Our goal is to study the slope of the solar UV spectrum in two UV bands important for the stratospheric ozone production. In order to investigate the solar spectral variability, we use SOLSTICE (the Solar Stellar Irradiance Comparison Experiment) data onboard Solar Radiation and Climate Experiment (SORCE) satellite. Data sets used are far UV (115-180nm) and middle UV (180-310nm), as well as the Mg II index (the Bremen composite). We introduce the SOLSTICE [FUV - MUV] colour to study the solar spectral characteristics, as well as analysis of the colour versus Mg II index. To isolate the 11-year scale variation, we used the Empirical Mode decomposition (EMD) on the data sets. The [FUV - MUV] colour strongly correlates with the Mg II index. More in detail, the [FUV - MUV] colour shows a time dependent behavior when plotted versus Mg II index. To explain this dependence we hypothesize an efficiency reduction of SOLSTICE FUV irradiance using an exponential aging law.