# Posts Tagged Monte Carlo

## Recent Postings from Monte Carlo

### Second large-scale Monte Carlo study for the Cherenkov Telescope Array

The Cherenkov Telescope Array (CTA) represents the next generation of ground based instruments for Very High Energy gamma-ray astronomy. It is expected to improve on the sensitivity of current instruments by an order of magnitude and provide energy coverage from 20 GeV to more than 200 TeV. In order to achieve these ambitious goals Monte Carlo (MC) simulations play a crucial role, guiding the design of CTA. Here, results of the second large-scale MC production are reported, providing a realistic estimation of feasible array candidates for both Northern and Sourthern Hemisphere sites performance, placing CTA capabilities into the context of the current generation of High Energy $\gamma$-ray detectors.

### Self-consistent modelling of line-driven hot-star winds with Monte Carlo radiation hydrodynamics

Radiative pressure exerted by line interactions is a prominent driver of outflows in astrophysical systems, being at work in the outflows emerging from hot stars or from the accretion discs of cataclysmic variables, massive young stars and active galactic nuclei. In this work, a new radiation hydrodynamical approach to model line-driven hot-star winds is presented. By coupling a Monte Carlo radiative transfer scheme with a finite-volume fluid dynamical method, line-driven mass outflows may be modelled self-consistently, benefiting from the advantages of Monte Carlo techniques in treating multi-line effects, such as multiple scatterings, and in dealing with arbitrary multidimensional configurations. In this work, we introduce our approach in detail by highlighting the key numerical techniques and verifying their operation in a number of simplified applications, specifically in a series of self-consistent, one-dimensional, Sobolev-type, hot-star wind calculations. The utility and accuracy of our approach is demonstrated by comparing the obtained results with the predictions of various formulations of the so-called CAK theory and by confronting the calculations with modern sophisticated techniques of predicting the wind structure. Using these calculations, we also point out some useful diagnostic capabilities our approach provides. Finally we discuss some of the current limitations of our method, some possible extensions and potential future applications.

### Phenomenology of Large Extra Dimensions Models at Hadrons Colliders using Monte Carlo Techniques (Spin-2 Graviton)

Large Extra Dimensions Models have been proposed to remove the hierarchy problem and give an explanation why the gravity is so much weaker than the other three forces. In this work, we present an analysis of Monte Carlo data events for new physics signatures of spin-2 Graviton in context of ADD model with total dimensions $D=4+\delta,$ $\delta = 1,2,3,4,5,6$ where $\delta$ is the extra special dimension, this model involves missing momentum $P_{T}^{miss}$ in association with jet in the final state via the process $pp(\bar{p}) \rightarrow G+jet$, Also, we present an analysis in context of the RS model with 5-dimensions via the process $pp(\bar{p}) \rightarrow G+jet$, $G \rightarrow e^{+}e^{-}$ with final state $e^{+}e^{-}+jet$. We used Monte Carlo event generator Pythia8 to produce efficient signal selection rules at the Large Hadron Collider with $\sqrt{s}$=14TeV and at the Tevatron $\sqrt{s}$=1.96TeV .

### Analytic Boosted Boson Discrimination [Replacement]

Observables which discriminate boosted topologies from massive QCD jets are of great importance for the success of the jet substructure program at the Large Hadron Collider. Such observables, while both widely and successfully used, have been studied almost exclusively with Monte Carlo simulations. In this paper we present the first all-orders factorization theorem for a two-prong discriminant based on a jet shape variable, $D_2$, valid for both signal and background jets. Our factorization theorem simultaneously describes the production of both collinear and soft subjets, and we introduce a novel zero-bin procedure to correctly describe the transition region between these limits. By proving an all orders factorization theorem, we enable a systematically improvable description, and allow for precision comparisons between data, Monte Carlo, and first principles QCD calculations for jet substructure observables. Using our factorization theorem, we present numerical results for the discrimination of a boosted $Z$ boson from massive QCD background jets. We compare our results with Monte Carlo predictions which allows for a detailed understanding of the extent to which these generators accurately describe the formation of two-prong QCD jets, and informs their usage in substructure analyses. Our calculation also provides considerable insight into the discrimination power and calculability of jet substructure observables in general.

### Analytic Boosted Boson Discrimination

Observables which discriminate boosted topologies from massive QCD jets are of great importance for the success of the jet substructure program at the Large Hadron Collider. Such observables, while both widely and successfully used, have been studied almost exclusively with Monte Carlo simulations. In this paper we present the first all-orders factorization theorem for a two-prong discriminant based on a jet shape variable, $D_2$, valid for both signal and background jets. Our factorization theorem simultaneously describes the production of both collinear and soft subjets, and we introduce a novel zero-bin procedure to correctly describe the transition region between these limits. By proving an all orders factorization theorem, we enable a systematically improvable description, and allow for precision comparisons between data, Monte Carlo, and first principles QCD calculations for jet substructure observables. Using our factorization theorem, we present numerical results for the discrimination of a boosted $Z$ boson from massive QCD background jets. We compare our results with Monte Carlo predictions which allows for a detailed understanding of the extent to which these generators accurately describe the formation of two-prong QCD jets, and informs their usage in substructure analyses. Our calculation also provides considerable insight into the discrimination power and calculability of jet substructure observables in general.

### Analytic Boosted Boson Discrimination [Cross-Listing]

Observables which discriminate boosted topologies from massive QCD jets are of great importance for the success of the jet substructure program at the Large Hadron Collider. Such observables, while both widely and successfully used, have been studied almost exclusively with Monte Carlo simulations. In this paper we present the first all-orders factorization theorem for a two-prong discriminant based on a jet shape variable, $D_2$, valid for both signal and background jets. Our factorization theorem simultaneously describes the production of both collinear and soft subjets, and we introduce a novel zero-bin procedure to correctly describe the transition region between these limits. By proving an all orders factorization theorem, we enable a systematically improvable description, and allow for precision comparisons between data, Monte Carlo, and first principles QCD calculations for jet substructure observables. Using our factorization theorem, we present numerical results for the discrimination of a boosted $Z$ boson from massive QCD background jets. We compare our results with Monte Carlo predictions which allows for a detailed understanding of the extent to which these generators accurately describe the formation of two-prong QCD jets, and informs their usage in substructure analyses. Our calculation also provides considerable insight into the discrimination power and calculability of jet substructure observables in general.

### Analytic Boosted Boson Discrimination [Replacement]

Observables which discriminate boosted topologies from massive QCD jets are of great importance for the success of the jet substructure program at the Large Hadron Collider. Such observables, while both widely and successfully used, have been studied almost exclusively with Monte Carlo simulations. In this paper we present the first all-orders factorization theorem for a two-prong discriminant based on a jet shape variable, $D_2$, valid for both signal and background jets. Our factorization theorem simultaneously describes the production of both collinear and soft subjets, and we introduce a novel zero-bin procedure to correctly describe the transition region between these limits. By proving an all orders factorization theorem, we enable a systematically improvable description, and allow for precision comparisons between data, Monte Carlo, and first principles QCD calculations for jet substructure observables. Using our factorization theorem, we present numerical results for the discrimination of a boosted $Z$ boson from massive QCD background jets. We compare our results with Monte Carlo predictions which allows for a detailed understanding of the extent to which these generators accurately describe the formation of two-prong QCD jets, and informs their usage in substructure analyses. Our calculation also provides considerable insight into the discrimination power and calculability of jet substructure observables in general.

### Evolution kernels for parton shower Monte Carlo

We report on re-calculation of the next-to-leading order DGLAP evolution kernels performed in a scheme suited for Monte Carlo simulations of parton cascades (parton showers).

### A User-Friendly Dark Energy Model Generator [Replacement]

We provide software with a graphical user interface to calculate the phenomenology of a wide class of dark energy models featuring multiple scalar fields. The user chooses a subclass of models and, if desired, initial conditions, or else a range of initial parameters for Monte Carlo. The code calculates the energy density of components in the universe, the equation of state of dark energy, and the linear growth of density perturbations, all as a function of redshift and scale factor. The output also includes an approximate conversion into the average equation of state, as well as the common $(w_0, w_a)$ parametrization. The code is available here: http://github.com/kahinton/Dark-Energy-UI-and-MC

### A User-Friendly Dark Energy Model Generator

We provide software with a graphical user interface to calculate the phenomenology of a wide class of dark energy models featuring multiple scalar fields. The user chooses a subclass of models and, if desired, initial conditions, or else a range of initial parameters for Monte Carlo. The code calculates the energy density of components in the universe, the equation of state of dark energy, and the linear growth of density perturbations, all as a function of redshift and scale factor. The output also includes an approximate conversion into the average equation of state, as well as the common $(w_0, w_a)$ parametrization. The code is available here: http://github.com/kahinton/Dark-Energy-UI-and-MC

### SKIRT: the design of a suite of input models for Monte Carlo radiative transfer simulations

The Monte Carlo method is the most popular technique to perform radiative transfer simulations in a general 3D geometry. The algorithms behind and acceleration techniques for Monte Carlo radiative transfer are discussed extensively in the literature, and many different Monte Carlo codes are publicly available. On the contrary, the design of a suite of components that can be used for the distribution of sources and sinks in radiative transfer codes has received very little attention. The availability of such models, with different degrees of complexity, has many benefits. For example, they can serve as toy models to test new physical ingredients, or as parameterised models for inverse radiative transfer fitting. For 3D Monte Carlo codes, this requires algorithms to efficiently generate random positions from 3D density distributions. We describe the design of a flexible suite of components for the Monte Carlo radiative transfer code SKIRT. The design is based on a combination of basic building blocks (which can be either analytical toy models or numerical models defined on grids or a set of particles) and the extensive use of decorators that combine and alter these building blocks to more complex structures. For a number of decorators, e.g. those that add spiral structure or clumpiness, we provide a detailed description of the algorithms that can be used to generate random positions. Advantages of this decorator-based design include code transparency, the avoidance of code duplication, and an increase in code maintainability. Moreover, since decorators can be chained without problems, very complex models can easily be constructed out of simple building blocks. Finally, based on a number of test simulations, we demonstrate that our design using customised random position generators is superior to a simpler design based on a generic black-box random position generator.

### An Unbiased Hessian Representation for Monte Carlo PDFs [Replacement]

We develop a methodology for the construction of a Hessian representation of Monte Carlo sets of parton distributions, based on the use of a subset of the Monte Carlo PDF replicas as an unbiased linear basis, and of a genetic algorithm for the determination of the optimal basis. We validate the methodology by first showing that it faithfully reproduces a native Monte Carlo PDF set (NNPDF3.0), and then, that if applied to Hessian PDF set (MMHT14) which was transformed into a Monte Carlo set, it gives back the starting PDFs with minimal information loss. We then show that, when applied to a large Monte Carlo PDF set obtained as combination of several underlying sets, the methodology leads to a Hessian representation in terms of a rather smaller set of parameters (CMC-H PDFs), thereby providing an alternative implementation of the recently suggested Meta-PDF idea and a Hessian version of the recently suggested PDF compression algorithm (CMC-PDFs). The mc2hessian conversion code is made publicly available together with (through LHAPDF6) a Hessian representations of the NNPDF3.0 set, and the CMC-H PDF set.

### An Unbiased Hessian Representation for Monte Carlo PDFs

We develop a methodology for the construction of a Hessian representation of Monte Carlo sets of parton distributions, based on the use of a subset of the Monte Carlo PDF replicas as an unbiased linear basis, and of a genetic algorithm for the determination of the optimal basis. We validate the methodology by first showing that it faithfully reproduces a native Monte Carlo PDF set (NNPDF3.0), and then, that if applied to Hessian PDF set (MMHT14) which was transformed into a Monte Carlo set, it gives back the starting PDFs with minimal information loss. We then show that, when applied to a large Monte Carlo PDF set obtained as combination of several underlying sets, the methodology leads to a Hessian representation in terms of a rather smaller set of parameters (CMC-H PDFs), thereby providing an alternative implementation of the recently suggested Meta-PDF idea and a Hessian version of the recently suggested PDF compression algorithm (CMC-PDFs). The mc2hessian conversion code is made publicly available together with (through LHAPDF6) a Hessian representations of the NNPDF3.0 set, and the CMC-H PDF set.

### An Unbiased Hessian Representation for Monte Carlo PDFs [Replacement]

We develop a methodology for the construction of a Hessian representation of Monte Carlo sets of parton distributions, based on the use of a subset of the Monte Carlo PDF replicas as an unbiased linear basis, and of a genetic algorithm for the determination of the optimal basis. We validate the methodology by first showing that it faithfully reproduces a native Monte Carlo PDF set (NNPDF3.0), and then, that if applied to Hessian PDF set (MMHT14) which was transformed into a Monte Carlo set, it gives back the starting PDFs with minimal information loss. We then show that, when applied to a large Monte Carlo PDF set obtained as combination of several underlying sets, the methodology leads to a Hessian representation in terms of a rather smaller set of parameters (CMC-H PDFs), thereby providing an alternative implementation of the recently suggested Meta-PDF idea and a Hessian version of the recently suggested PDF compression algorithm (CMC-PDFs). The mc2hessian conversion code is made publicly available together with (through LHAPDF6) a Hessian representations of the NNPDF3.0 set, and the CMC-H PDF set.

### Monte Carlo Method for Calculating Oxygen Abundances and Their Uncertainties from Strong-Line Flux Measurements

We present the open-source Python code pyMCZ that determines oxygen abundance and its distribution from strong emission lines in the standard metallicity scales, based on the original IDL code of Kewley & Dopita (2002) with updates from Kewley & Ellison (2008), and expanded to include more recently developed scales. The standard strong-line diagnostics have been used to estimate the oxygen abundance in the interstellar medium through various emission line ratios in many areas of astrophysics, including galaxy evolution and supernova host galaxy studies. We introduce a Python implementation of these methods that, through Monte Carlo (MC) sampling, better characterizes the statistical reddening-corrected oxygen abundance confidence region. Given line flux measurements and their uncertainties, our code produces synthetic distributions for the oxygen abundance in up to 13 metallicity scales simultaneously, as well as for E(B-V), and estimates their median values and their 66% confidence regions. In addition, we provide the option of outputting the full MC distributions, and their kernel density estimates. We test our code on emission line measurements from a sample of nearby supernova host galaxies ($z<0.15$) and compare our metallicity results with those from previous methods. We show that our metallicity estimates are consistent with previous methods but yield smaller uncertainties. We also offer visualization tools to assess the spread of the oxygen abundance in the different scales, as well as the shape of the estimated oxygen abundance distribution in each scale, and develop robust metrics for determining the appropriate MC sample size. The code is open access and open source and can be found at https://github.com/nyusngroup/pyMCZ

### Monte-Carlo approach to particle-field interactions and the kinetics of the chiral phase transition

The kinetics of the chiral phase transition is studied within a linear quark-meson-$\sigma$ model, using a Monte-Carlo approach to semiclassical particle-field dynamics. The meson fields are described on the mean-field level and quarks and antiquarks as ensembles of test particles. Collisions between quarks and antiquarks as well as the $q\overline{q}$ annihilation to $\sigma$ mesons and the decay of $\sigma$ mesons is treated, using the corresponding transition-matrix elements from the underlying quantum field theory, obeying strictly the rule of detailed balance and energy-momentum conservation. The approach allows to study fluctuations without making ad hoc assumptions concerning the statistical nature of the random process as necessary in Langevin-Fokker-Planck frameworks.

### Polarisation spectral synthesis for Type Ia supernova explosion models

We present a Monte Carlo radiative transfer technique for calculating synthetic spectropolarimetry for multi-dimensional supernova explosion models. The approach utilises "virtual-packets" that are generated during the propagation of the Monte Carlo quanta and used to compute synthetic observables for specific observer orientations. Compared to extracting synthetic observables by direct binning of emergent Monte Carlo quanta, this virtual-packet approach leads to a substantial reduction in the Monte Carlo noise. This is vital for calculating synthetic spectropolarimetry (since the degree of polarisation is typically very small) but also useful for calculations of light curves and spectra. We first validate our approach via application of an idealised test code to simple geometries. We then describe its implementation in the Monte Carlo radiative transfer code ARTIS and present test calculations for simple models for Type Ia supernovae. Specifically, we use the well-known one-dimensional W7 model to verify that our scheme can accurately recover zero polarisation from a spherical model, and to demonstrate the reduction in Monte Carlo noise compared to a simple packet-binning approach. To investigate the impact of aspherical ejecta on the polarisation spectra, we then use ARTIS to calculate synthetic observables for prolate and oblate ellipsoidal models with Type Ia supernova compositions.

### Transition redshift in $f(T)$ cosmology and observational constraints

We extract constraints on the transition redshift $z_{tr}$, determining the onset of cosmic acceleration, predicted by an effective cosmographic construction, in the framework of $f(T)$ gravity. In particular, employing cosmography we obtain bounds on the viable $f(T)$ forms and their derivatives. Since this procedure is model independent, as long as the scalar curvature is fixed, we are able to determine intervals for $z_{tr}$. In this way we guarantee that the Solar-System constraints are preserved and moreover we extract bounds on the transition time and the free parameters of the scenario. We find that the transition redshifts predicted by $f(T)$ cosmology, although compatible with the standard $\Lambda$CDM predictions, are slightly smaller. Finally, in order to obtain observational constraints on $f(T)$ cosmology, we perform a Monte Carlo fitting using supernova data, involving the most recent union 2.1 data set.

### Transition redshift in $f(T)$ cosmology and observational constraints [Replacement]

We extract constraints on the transition redshift $z_{tr}$, determining the onset of cosmic acceleration, predicted by an effective cosmographic construction, in the framework of $f(T)$ gravity. In particular, employing cosmography we obtain bounds on the viable $f(T)$ forms and their derivatives. Since this procedure is model independent, as long as the scalar curvature is fixed, we are able to determine intervals for $z_{tr}$. In this way we guarantee that the Solar-System constraints are preserved and moreover we extract bounds on the transition time and the free parameters of the scenario. We find that the transition redshifts predicted by $f(T)$ cosmology, although compatible with the standard $\Lambda$CDM predictions, are slightly smaller. Finally, in order to obtain observational constraints on $f(T)$ cosmology, we perform a Monte Carlo fitting using supernova data, involving the most recent union 2.1 data set.

### Transition redshift in $f(T)$ cosmology and observational constraints [Cross-Listing]

We extract constraints on the transition redshift $z_{tr}$, determining the onset of cosmic acceleration, predicted by an effective cosmographic construction, in the framework of $f(T)$ gravity. In particular, employing cosmography we obtain bounds on the viable $f(T)$ forms and their derivatives. Since this procedure is model independent, as long as the scalar curvature is fixed, we are able to determine intervals for $z_{tr}$. In this way we guarantee that the Solar-System constraints are preserved and moreover we extract bounds on the transition time and the free parameters of the scenario. We find that the transition redshifts predicted by $f(T)$ cosmology, although compatible with the standard $\Lambda$CDM predictions, are slightly smaller. Finally, in order to obtain observational constraints on $f(T)$ cosmology, we perform a Monte Carlo fitting using supernova data, involving the most recent union 2.1 data set.

### Transition redshift in $f(T)$ cosmology and observational constraints [Replacement]

We extract constraints on the transition redshift $z_{tr}$, determining the onset of cosmic acceleration, predicted by an effective cosmographic construction, in the framework of $f(T)$ gravity. In particular, employing cosmography we obtain bounds on the viable $f(T)$ forms and their derivatives. Since this procedure is model independent, as long as the scalar curvature is fixed, we are able to determine intervals for $z_{tr}$. In this way we guarantee that the Solar-System constraints are preserved and moreover we extract bounds on the transition time and the free parameters of the scenario. We find that the transition redshifts predicted by $f(T)$ cosmology, although compatible with the standard $\Lambda$CDM predictions, are slightly smaller. Finally, in order to obtain observational constraints on $f(T)$ cosmology, we perform a Monte Carlo fitting using supernova data, involving the most recent union 2.1 data set.

### Transition redshift in $f(T)$ cosmology and observational constraints [Replacement]

We extract constraints on the transition redshift $z_{tr}$, determining the onset of cosmic acceleration, predicted by an effective cosmographic construction, in the framework of $f(T)$ gravity. In particular, employing cosmography we obtain bounds on the viable $f(T)$ forms and their derivatives. Since this procedure is model independent, as long as the scalar curvature is fixed, we are able to determine intervals for $z_{tr}$. In this way we guarantee that the Solar-System constraints are preserved and moreover we extract bounds on the transition time and the free parameters of the scenario. We find that the transition redshifts predicted by $f(T)$ cosmology, although compatible with the standard $\Lambda$CDM predictions, are slightly smaller. Finally, in order to obtain observational constraints on $f(T)$ cosmology, we perform a Monte Carlo fitting using supernova data, involving the most recent union 2.1 data set.

### Transition redshift in $f(T)$ cosmology and observational constraints [Cross-Listing]

We extract constraints on the transition redshift $z_{tr}$, determining the onset of cosmic acceleration, predicted by an effective cosmographic construction, in the framework of $f(T)$ gravity. In particular, employing cosmography we obtain bounds on the viable $f(T)$ forms and their derivatives. Since this procedure is model independent, as long as the scalar curvature is fixed, we are able to determine intervals for $z_{tr}$. In this way we guarantee that the Solar-System constraints are preserved and moreover we extract bounds on the transition time and the free parameters of the scenario. We find that the transition redshifts predicted by $f(T)$ cosmology, although compatible with the standard $\Lambda$CDM predictions, are slightly smaller. Finally, in order to obtain observational constraints on $f(T)$ cosmology, we perform a Monte Carlo fitting using supernova data, involving the most recent union 2.1 data set.

### Binary population synthesis for the core-degenerate scenario of type Ia supernova progenitors [Replacement]

The core-degenerate (CD) scenario has been suggested to be a possible progenitor model of type Ia supernovae (SNe Ia), in which a carbon-oxygen white dwarf (CO WD) merges with the hot CO core of a massive asymptotic giant branch (AGB) star during their common-envelope phase. However, the SN Ia birthrates for this scenario are still uncertain. We conducted a detailed investigation into the CD scenario and then gave the birthrates for this scenario using a detailed Monte Carlo binary population synthesis approach. We found that the delay times of SNe Ia from this scenario are ~70Myrs-1400Myrs, which means that the CD scenario contributes to young SN Ia populations. The Galactic SN Ia birthrates for this scenario are in the range of ~7.4*10^{-5} yr^{-1}-3.7*10^{-4}yr^{-1}, which roughly accounts for ~2-10% of all SNe Ia. This indicates that, under the assumptions made here, the CD scenario only contributes a small portion of all SNe Ia, which is not consistent with the results of Ilkov & Soker (2013).

### Binary population synthesis for the core-degenerate scenario of type Ia supernova progenitors

The core-degenerate (CD) scenario has been suggested to be a possible progenitor model of type Ia supernovae (SNe Ia), in which a carbon-oxygen white dwarf (CO WD) merges with the hot CO core of a massive asymptotic giant branch (AGB) star during their common-envelope phase. However, the SN Ia birthrates for this scenario are still uncertain. We conducted a detailed investigation into the CD scenario and then gave the birthrates for this scenario using a detailed Monte Carlo binary population synthesis approach. We found that the delay times of SNe Ia from this scenario are ~70Myrs-1400Myrs, which means that the CD scenario contributes to young SN Ia populations. The Galactic SN Ia birthrates for this scenario are in the range of ~7.4*10^{-5} yr^{-1}-3.7*10^{-4}yr^{-1}, which roughly accounts for ~2-10% of all SNe Ia. This indicates that, under the assumptions made here, the CD scenario only contributes a small portion of all SNe Ia, which is not consistent with the results of Ilkov & Soker (2013).

### Radiation-hydrodynamical simulations of massive star formation using Monte Carlo radiative transfer: I. Algorithms and numerical methods

We present a set of new numerical methods that are relevant to calculating radiation pressure terms in hydrodynamics calculations, with a particular focus on massive star formation. The radiation force is determined from a Monte Carlo estimator and enables a complete treatment of the detailed microphysics, including polychromatic radiation and anisotropic scattering, in both the free-streaming and optically-thick limits. Since the new method is computationally demanding we have developed two new methods that speed up the algorithm. The first is a photon packet splitting algorithm that enables efficient treatment of the Monte Carlo process in very optically thick regions. The second is a parallelisation method that distributes the Monte Carlo workload over many instances of the hydrodynamic domain, resulting in excellent scaling of the radiation step. We also describe the implementation of a sink particle method that enables us to follow the accretion onto, and the growth of, the protostars. We detail the results of extensive testing and benchmarking of the new algorithms.

### Ab initio lattice results for Fermi polarons in two and three dimensions [Replacement]

We investigate the attractive Fermi polaron problem in two and three dimensions using non-perturbative Monte Carlo simulations. We introduce a new Monte Carlo algorithm called the impurity lattice Monte Carlo method. This algorithm samples the path integral in a computationally efficient manner and has only small sign oscillations for systems with a single impurity. To benchmark the method, we calculate the universal polaron energy in three dimensions in the scale-invariant unitarity limit and find agreement with published results. We then present the first fully non-perturbative calculations of the polaron energy in two dimensions in the limit of zero range interactions. We find evidence for a smooth crossover transition from a fermionic quasiparticle to a molecular state as a function of interaction strength with significant mixing between the two descriptions in the crossover region.

### The tau leptons theory and experimental data: Monte Carlo, fits, software and systematic errors

Status of tau lepton decay Monte Carlo generator TAUOLA is reviewed. Recent efforts on development of new hadronic currents are presented. Multitude new channels for anomalous tau decay modes and parametrization based on defaults used by BaBar collaboration are introduced. Also parametrization based on theoretical considerations are presented as an alternative. Lesson from comparison and fits to the BaBar and Belle data is recalled. It was found that as in the past, in particular at a time of comparisons with CLEO and ALEPH data, proper fitting, to as detailed as possible representation of the experimental data, is essential for appropriate developments of models of tau decays. In the later part of the presentation, use of the TAUOLA program for phenomenology of W,Z,H decays at LHC is adressed. Some new results, relevant for QED bremsstrahlung in such decays are presented as well.

### Monte Carlo error analyses of Spearman's rank test

Spearman’s rank correlation test is commonly used in astronomy to discern whether a set of two variables are correlated or not. Unlike most other quantities quoted in astronomical literature, the Spearman’s rank correlation coefficient is generally quoted with no attempt to estimate the errors on its value. This is a practice that would not be accepted for those other quantities, as it is often regarded that an estimate of a quantity without an estimate of its associated uncertainties is meaningless. This manuscript describes a number of easily implemented, Monte Carlo based methods to estimate the uncertainty on the Spearman’s rank correlation coefficient, or more precisely to estimate its probability distribution.

### Monte Carlo error analyses of Spearman's rank test [Replacement]

Spearman’s rank correlation test is commonly used in astronomy to discern whether a set of two variables are correlated or not. Unlike most other quantities quoted in astronomical literature, the Spearman’s rank correlation coefficient is generally quoted with no attempt to estimate the errors on its value. This is a practice that would not be accepted for those other quantities, as it is often regarded that an estimate of a quantity without an estimate of its associated uncertainties is meaningless. This manuscript describes a number of easily implemented, Monte Carlo based methods to estimate the uncertainty on the Spearman’s rank correlation coefficient, or more precisely to estimate its probability distribution.

### Cosmographic transition redshift in $f(R)$ gravity

We propose a strategy to infer the transition redshift $z_{da}$, which characterizes the passage through the universe decelerated to accelerated phases, in the framework $f(R)$ gravities. To this end, we numerically reconstruct $f(z)$, i.e. the corresponding $f(R)$ function re-expressed in terms of the redshift $z$ and we show how to match $f(z)$ with cosmography. In particular, we relate $f(z)$ and its derivatives to the cosmographic coefficients, i.e. $H_0, q_0$ and $j_0$ and demonstrate that its corresponding evolution may be framed by means of an effective logarithmic dark energy term $\Omega_X$, slightly departing from the case of a pure cosmological constant. Afterwards, we show that our model predicts viable transition redshift constraints, which agree with $\Lambda$CDM. To do so, we compute the corresponding $z_{da}$ in terms of cosmographic outcomes and find that $z_{da}\leq1$. Finally, we reproduce an effective $f(z)$ and show that this class of models is fairly well compatible with present-time data. To do so, we get numerical constraints employing Monte Carlo fits with the Union 2.1 supernova survey and with the Hubble measurement data set.

### Cosmographic transition redshift in $f(R)$ gravity [Cross-Listing]

We propose a strategy to infer the transition redshift $z_{da}$, which characterizes the passage through the universe decelerated to accelerated phases, in the framework $f(R)$ gravities. To this end, we numerically reconstruct $f(z)$, i.e. the corresponding $f(R)$ function re-expressed in terms of the redshift $z$ and we show how to match $f(z)$ with cosmography. In particular, we relate $f(z)$ and its derivatives to the cosmographic coefficients, i.e. $H_0, q_0$ and $j_0$ and demonstrate that its corresponding evolution may be framed by means of an effective logarithmic dark energy term $\Omega_X$, slightly departing from the case of a pure cosmological constant. Afterwards, we show that our model predicts viable transition redshift constraints, which agree with $\Lambda$CDM. To do so, we compute the corresponding $z_{da}$ in terms of cosmographic outcomes and find that $z_{da}\leq1$. Finally, we reproduce an effective $f(z)$ and show that this class of models is fairly well compatible with present-time data. To do so, we get numerical constraints employing Monte Carlo fits with the Union 2.1 supernova survey and with the Hubble measurement data set.

### Cosmographic transition redshift in $f(R)$ gravity [Cross-Listing]

We propose a strategy to infer the transition redshift $z_{da}$, which characterizes the passage through the universe decelerated to accelerated phases, in the framework $f(R)$ gravities. To this end, we numerically reconstruct $f(z)$, i.e. the corresponding $f(R)$ function re-expressed in terms of the redshift $z$ and we show how to match $f(z)$ with cosmography. In particular, we relate $f(z)$ and its derivatives to the cosmographic coefficients, i.e. $H_0, q_0$ and $j_0$ and demonstrate that its corresponding evolution may be framed by means of an effective logarithmic dark energy term $\Omega_X$, slightly departing from the case of a pure cosmological constant. Afterwards, we show that our model predicts viable transition redshift constraints, which agree with $\Lambda$CDM. To do so, we compute the corresponding $z_{da}$ in terms of cosmographic outcomes and find that $z_{da}\leq1$. Finally, we reproduce an effective $f(z)$ and show that this class of models is fairly well compatible with present-time data. To do so, we get numerical constraints employing Monte Carlo fits with the Union 2.1 supernova survey and with the Hubble measurement data set.

### Power Counting to Better Jet Observables

Optimized jet substructure observables for identifying boosted topologies will play an essential role in maximizing the physics reach of the Large Hadron Collider. Ideally, the design of discriminating variables would be informed by analytic calculations in perturbative QCD. Unfortunately, explicit calculations are often not feasible due to the complexity of the observables used for discrimination, and so many validation studies rely heavily, and solely, on Monte Carlo. In this paper we show how methods based on the parametric power counting of the dynamics of QCD, familiar from effective theory analyses, can be used to design, understand, and make robust predictions for the behavior of jet substructure variables. As a concrete example, we apply power counting for discriminating boosted Z bosons from massive QCD jets using observables formed from the n-point energy correlation functions. We show that power counting alone gives a definite prediction for the observable that optimally separates the background-rich from the signal-rich regions of phase space. Power counting can also be used to understand effects of phase space cuts and the effect of contamination from pile-up, which we discuss. As these arguments rely only on the parametric scaling of QCD, the predictions from power counting must be reproduced by any Monte Carlo, which we verify using Pythia8 and Herwig++. We also use the example of quark versus gluon discrimination to demonstrate the limits of the power counting technique.

### Conformal field theories at non-zero temperature: operator product expansions, Monte Carlo, and holography [Replacement]

We compute the non-zero temperature conductivity of conserved flavor currents in conformal field theories (CFTs) in 2+1 spacetime dimensions. At frequencies much greater than the temperature, $\hbar\omega>> k_B T$, the $\omega$ dependence can be computed from the operator product expansion (OPE) between the currents and operators which acquire a non-zero expectation value at T > 0. Such results are found to be in excellent agreement with quantum Monte Carlo studies of the O(2) Wilson-Fisher CFT. Results for the conductivity and other observables are also obtained in vector 1/N expansions. We match these large $\omega$ results to the corresponding correlators of holographic representations of the CFT: the holographic approach then allows us to extrapolate to small $\hbar \omega/(k_B T)$. Other holographic studies implicitly only used the OPE between the currents and the energy-momentum tensor, and this yields the correct leading large $\omega$ behavior for a large class of CFTs. However, for the Wilson-Fisher CFT a relevant "thermal" operator must also be considered, and then consistency with the Monte Carlo results is obtained without a previously needed ad hoc rescaling of the T value. We also establish sum rules obeyed by the conductivity of a wide class of CFTs.

### Conformal field theories at non-zero temperature: operator product expansions, Monte Carlo, and holography [Replacement]

We compute the non-zero temperature conductivity of conserved flavor currents in conformal field theories (CFTs) in 2+1 spacetime dimensions. At frequencies much greater than the temperature, $\hbar\omega>> k_B T$, the $\omega$ dependence can be computed from the operator product expansion (OPE) between the currents and operators which acquire a non-zero expectation value at T > 0. Such results are found to be in excellent agreement with quantum Monte Carlo studies of the O(2) Wilson-Fisher CFT. Results for the conductivity and other observables are also obtained in vector 1/N expansions. We match these large $\omega$ results to the corresponding correlators of holographic representations of the CFT: the holographic approach then allows us to extrapolate to small $\hbar \omega/(k_B T)$. Other holographic studies implicitly only used the OPE between the currents and the energy-momentum tensor, and this yields the correct leading large $\omega$ behavior for a large class of CFTs. However, for the Wilson-Fisher CFT a relevant "thermal" operator must also be considered, and then consistency with the Monte Carlo results is obtained without a previously needed ad hoc rescaling of the T value. We also establish sum rules obeyed by the conductivity of a wide class of CFTs.

### Conformal field theories at non-zero temperature: operator product expansions, Monte Carlo, and holography [Cross-Listing]

We compute the non-zero temperature conductivity of conserved flavor currents in conformal field theories (CFTs) in 2+1 spacetime dimensions. At frequencies much greater than the temperature, $\hbar\omega>> k_B T$, the $\omega$ dependence can be computed from the operator product expansion (OPE) between the currents and operators which acquire a non-zero expectation value at T > 0. Such results are found to be in excellent agreement with quantum Monte Carlo studies of the O(2) Wilson-Fisher CFT. Results for the conductivity and other observables are also obtained in vector 1/N expansions. We match these large $\omega$ results to the corresponding correlators of holographic representations of the CFT: the holographic approach then allows us to extrapolate to small $\hbar \omega/(k_B T)$. Other holographic studies implicitly only used the OPE between the currents and the energy-momentum tensor, and this yields the correct leading large $\omega$ behavior for a large class of CFTs. However, for the Wilson-Fisher CFT a relevant "thermal" operator must also be considered, and then consistency with the Monte Carlo results is obtained without a previously needed ad hoc rescaling of the T value. We also establish sum rules obeyed by the conductivity of a wide class of CFTs.

### Conformal field theories at non-zero temperature: operator product expansions, Monte Carlo, and holography [Cross-Listing]

We compute the non-zero temperature conductivity of conserved flavor currents in conformal field theories (CFTs) in 2+1 spacetime dimensions. At frequencies much greater than the temperature, $\hbar\omega>> k_B T$, the $\omega$ dependence can be computed from the operator product expansion (OPE) between the currents and operators which acquire a non-zero expectation value at T > 0. Such results are found to be in excellent agreement with quantum Monte Carlo studies of the O(2) Wilson-Fisher CFT. Results for the conductivity and other observables are also obtained in vector 1/N expansions. We match these large $\omega$ results to the corresponding correlators of holographic representations of the CFT: the holographic approach then allows us to extrapolate to small $\hbar \omega/(k_B T)$. Other holographic studies implicitly only used the OPE between the currents and the energy-momentum tensor, and this yields the correct leading large $\omega$ behavior for a large class of CFTs. However, for the Wilson-Fisher CFT a relevant "thermal" operator must also be considered, and then consistency with the Monte Carlo results is obtained without a previously needed ad hoc rescaling of the T value. We also establish sum rules obeyed by the conductivity of a wide class of CFTs.

### Signal Attenuation Curve for Different Surface Detector Arrays

Modern cosmic ray experiments consisting of large array of particle detectors measure the signals of electromagnetic or muon components or their combination. The correction for an amount of atmosphere passed is applied to the surface detector signal before its conversion to the shower energy. Either Monte Carlo based approach assuming certain composition of primaries or indirect estimation using real data and assuming isotropy of arrival directions can be used. Toy surface arrays of different sensitivities to electromagnetic and muon components are assumed in MC simulations to study effects imposed on attenuation curves for varying composition or possible high energy anisotropy. The possible sensitivity of the attenuation curve to the mass composition is also tested for different array types focusing on a future apparatus that can separate muon and electromagnetic component signals.

### Heavy Ions Collision evolution modeling with ECHO-QGP

We present a numerical code modeling the evolution of the medium formed in relativistic heavy ion collisions, ECHO-QGP. The code solves relativistic hydrodynamics in $(3+1)-$D, with dissipative terms included within the framework of Israel-Stewart theory; it can work both in Minkowskian and in Bjorken coordinates. Initial conditions are provided through an implementation of the Glauber model (both Optical and Monte Carlo), while freezeout and particle generation are based on the Cooper-Frye prescription. The code is validated against several test problems and shows remarkable stability and accuracy with the combination of a conservative (shock-capturing) approach and the high-order methods employed. In particular it beautifully agrees with the semi-analytic solution known as Gubser flow, both in the ideal and in the viscous Israel-Stewart case, up to very large times and without any ad hoc tuning of the algorithm.

### Properties of Long Gamma Ray Burst Progenitors in Cosmological Simulations

We study the nature of long gamma ray burst (LGRB) progenitors using cosmological simulations of structure formation and galactic evolution. LGRBs are potentially excellent tracers of stellar evolution in the early universe. We developed a Monte Carlo numerical code which generates LGRBs coupled to cosmological simulations. The simulations allows us to follow the ormation of galaxies self-consistently. We model the detectability of LGRBs and their host galaxies in order to compare results with observational data obtained by high-energy satellites. Our code also includes stochastic effects in the observed rate of LGRBs.

### The density of states approach to dense quantum systems [Cross-Listing]

We develop a first-principle generalised density of state method for studying numerically quantum field theories with a complex action. As a proof of concept, we show that with our approach we can solve numerically the strong sign problem of the $Z_3$ spin model at finite density. Our results are confirmed by standard simulations of the theory dual to the considered model, which is free from a sign problem. Our method opens new perspectives on ab initio simulations of cold dense quantum systems, and in particular of Yang-Mills theories with matter at finite densities, for which Monte Carlo based importance sampling are unable to produce sufficiently accurate results.

### The density of states approach to dense quantum systems [Cross-Listing]

We develop a first-principle generalised density of state method for studying numerically quantum field theories with a complex action. As a proof of concept, we show that with our approach we can solve numerically the strong sign problem of the $Z_3$ spin model at finite density. Our results are confirmed by standard simulations of the theory dual to the considered model, which is free from a sign problem. Our method opens new perspectives on ab initio simulations of cold dense quantum systems, and in particular of Yang-Mills theories with matter at finite densities, for which Monte Carlo based importance sampling are unable to produce sufficiently accurate results.

### Small-$x$ dynamics in forward-central dijet decorrelations at the LHC [Replacement]

We provide a description, within the High Energy Factorization formalism, of central-forward dijet decorrelation data measured by the CMS experiment and the predictions for nuclear modification ratio~$R_{pA}$ in p+Pb collisions. In our study, we use the unintegrated gluon density derived from the BFKL and BK equations supplemented with subleading corrections and a hard scale dependence. The latter is introduced at the final step of the calculation by reweighting the Monte Carlo generated events using suitable Sudakov form factors, without changing the total cross section. We achieve a good description of data in the whole region of the azimuthal angle.

### Monte Carlo Radiation Hydrodynamics with Implicit Methods

We explore the application of Monte Carlo transport methods to solving coupled radiation-hydrodynamics problems. We use a time-dependent, frequency-dependent, 3-dimensional radiation transport code, that is special relativistic and includes some detailed microphysical interactions such as resonant line scattering. We couple the transport code to two different 1-dimensional (non-relativistic) hydrodynamics solvers: a spherical Lagrangian scheme and a Eulerian Godunov solver. The gas-radiation energy coupling is treated implicitly, allowing us to take hydrodyanimcal time-steps that are much longer than the radiative cooling time. We validate the code and assess its performance using a suite of radiation hydrodynamical test problems, including ones in the radiation energy dominated regime. We also develop techniques that reduce the noise of the Monte Carlo estimated radiation force by using the spatial divergence of the radiation pressure tensor. The results suggest that Monte Carlo techniques hold promise for simulating the multi-dimensional radiation hydrodynamics of astrophysical systems.

### Monte Carlo Radiation Hydrodynamics with Implicit Methods [Replacement]

We explore the application of Monte Carlo transport methods to solving coupled radiation-hydrodynamics problems. We use a time-dependent, frequency-dependent, 3-dimensional radiation transport code, that is special relativistic and includes some detailed microphysical interactions such as resonant line scattering. We couple the transport code to two different 1-dimensional (non-relativistic) hydrodynamics solvers: a spherical Lagrangian scheme and a Eulerian Godunov solver. The gas-radiation energy coupling is treated implicitly, allowing us to take hydrodyanimcal time-steps that are much longer than the radiative cooling time. We validate the code and assess its performance using a suite of radiation hydrodynamical test problems, including ones in the radiation energy dominated regime. We also develop techniques that reduce the noise of the Monte Carlo estimated radiation force by using the spatial divergence of the radiation pressure tensor. The results suggest that Monte Carlo techniques hold promise for simulating the multi-dimensional radiation hydrodynamics of astrophysical systems.

### (MC)**3 -- a Multi-Channel Markov Chain Monte Carlo algorithm for phase-space sampling

A new Monte Carlo algorithm for phase-space integration, named (MC)**3, is presented. It is based on Markov Chain Monte Carlo techniques but at the same time incorporates prior knowledge about the target distribution in the form of suitable phase-space mappings from a corresponding Multi-Channel importance sampling Monte Carlo. The combined approach inherits the benefits of both techniques while typical drawbacks of either solution get ameliorated.

### (MC)**3 -- a Multi-Channel Markov Chain Monte Carlo algorithm for phase-space sampling [Replacement]

A new Monte Carlo algorithm for phase-space sampling, named (MC)**3, is presented. It is based on Markov Chain Monte Carlo techniques but at the same time incorporates prior knowledge about the target distribution in the form of suitable phase-space mappings from a corresponding Multi-Channel importance sampling Monte Carlo. The combined approach inherits the benefits of both techniques while typical drawbacks of either solution get ameliorated.

### (MC)**3 -- a Multi-Channel Markov Chain Monte Carlo algorithm for phase-space sampling [Replacement]

A new Monte Carlo algorithm for phase-space sampling, named (MC)**3, is presented. It is based on Markov Chain Monte Carlo techniques but at the same time incorporates prior knowledge about the target distribution in the form of suitable phase-space mappings from a corresponding Multi-Channel importance sampling Monte Carlo. The combined approach inherits the benefits of both techniques while typical drawbacks of either solution get ameliorated.

### (MC)**3 -- a Multi-Channel Markov Chain Monte Carlo algorithm for phase-space sampling [Replacement]

A new Monte Carlo algorithm for phase-space sampling, named (MC)**3, is presented. It is based on Markov Chain Monte Carlo techniques but at the same time incorporates prior knowledge about the target distribution in the form of suitable phase-space mappings from a corresponding Multi-Channel importance sampling Monte Carlo. The combined approach inherits the benefits of both techniques while typical drawbacks of either solution get ameliorated.