How do you think has Ronnie's character has shaped?
I think as an actor I'm immensely pleased with how the character has evolved. In just 60 episodes I have portrayed such emotions. I have been friends , I have been in love and then been sarcastic and finally coming to terms with the truth. I'm really happy with the character.
Ronnie and Kabir are two most important men in Ananya's life. Will there be a face of?
I wish this was to reach the channel or the production house. I would love for there to be a scene like this. So far , we have not shot any such scenes.
Fans feel that Ronnie and Richa will make a great couple. Can we see romance kindling there?
Would you personally want Ronnie to have a love interest in the show?
Like I said it is not possible since it is a finite show. It will be out of character for Ronnie to move on so quickly. He still loves Ananya but she friend-zoned him. He is a supportive friend and has made peace with Ananya's decision to date Kabir.
Do the low TRP's affect you as an actor?
Initially of course we all felt bad and slightly disheartened. Now that we meet the fans and we get feedback on the Internet we realize it is a well loved show. We have a very dedicated fan base and they are loving the show. As an actor it is such a high to have fans and viewers appreciate your performance.
Ankita Bhargava was expected to join the cast last week but she has not so far. Do TRP's affect the storyline?
No, in this case it is not true. It all depends on the screenplay. At times it stretches for two episodes and alters the joining of a new actor. However, we have the screenplay now with Shreya and she joins soon.
Any specific fan encounter that you can recall?
Oh yea! This one is very interesting. I had gone with friends to Monkey Bar and a girl walked up to me. She told me how she loved me as Ronnie and then added that if I ever got into trouble I should give her a call as she is a practicing lawyer. She shared her number with me !
Who is the prankster on the sets?
We are just a mad group to hang out with it. We are constantly upto something. As such we don't play pranks. However, I pull out my phone whenever someone is upto something and quickly record it and post it on our personal group. Everyone knows if something goes wrong on the sets, Puru is to be blamed!
What have your favorite scenes been so far?
I loved the proposal scene. It was my first rain sequence and was a fantastic experience. I like all my emotional scenes with Kritika. We are both from Delhi and get along really well!
How similar are you to your on screen character ?
I feel we both have a good heart and can be sarcastic when we want to be. Other than that I feel I'm very different from Ronnie. If I loved a girl I wouldn't have given up so easily. I would have fought for her and probably would get her as well.
What has been your most expensive purchase?
I'm a spendthrift. I don't think I have ever stuck to the budget when I go out shopping. I don't know but I always end up liking the most expensive thing in the shop. My mother recently surprised me with a car and she got it to the sets of Reporters.
Do you think finite shows have a future in India?
I'm hoping they work. As an actor I have been associated with mainly finite shows and they are so rewarding. Crispy story lines and well etched out characters! It will take time but I feel people should keep trying. I know so many actors who want to be a part of Reporters and a similar series.
Don't forget to catch him and rest of the 'Reporters' on Sony TV at 9pm Monday to Friday!


text:
Hortonville man charged, accused of soliciting nude images of girls on Facebook
"We have allegations that, quite frankly, span the globe," Outagamie County Court Commissioner said Friday.
Hortonville man charged, accused of soliciting nude images of girls on Facebook "We have allegations that, quite frankly, span the globe," Outagamie County Court Commissioner said Friday. Check out this story on postcrescent.com: https://post.cr/2GvF0dC
Alison Dirr, USA TODAY NETWORK-Wisconsin Published 5:04 p.m. CT May 18, 2018
If you can help solve a crime, please contact the appropriate police agency. Duke Behnke/USA TODAY NETWORK-Wisconsin
Jason Junek(Photo: Photo courtesy of Outagamie County Sheriff's Department)
A Hortonville man was charged Friday in a case in which investigators say they received "several" cyber tips about a man soliciting and receiving nude images of girls through Facebook, the Outagamie County Sheriff's Department said in a statement Friday.
The tips were received early this year.
The investigation revealed the 41-year-old man used a fake Facebook profile to communicate with girls in Texas, Illinois, Colorado, Nebraska and the United Kingdom, the department said.
"We have allegations that, quite frankly, span the globe," Outagamie County Court Commissioner Brian Figy said Friday in setting a $25,000 cash bond for Jason W. Junek. A preliminary hearing is set for May 25.
He was charged with four counts of sexual exploitation of a child, two counts of attempted sexual exploitation of a child, three counts of exposing a child to harmful material and a single count of soliciting an intimate representation from a minor, according to court records.
According to the criminal complaint:
Tips from the National Center for Missing and Exploited Children led an Outagamie County Sheriff's sergeant to Junek's IP address and ultimately his home address. The investigator determined that Junek was "communicating with and soliciting nude photographs from several juvenile females who lived in various locations around the world," the complaint states.
A January tip reported that a user in the United Kingdom had uploaded an image of suspected child pornography to a user in Hortonville. There were 71 images and 20 videos attached to the cyber tip.
The 13-year-old girl in the images and videos said that the other Facebook user had asked her to send nude images. The girl and the other user exchanged sexually explicit photos and videos in September and October 2017.
The second victim was a 16-year-old Illinois girl who had sent "several" images and videos to an individual in Hortonville with Junek's IP address. She said the user she sent the explicit content that said he was 17 years old. She said he had asked repeatedly for the nude photos and videos and had once sent her an explicit video.
In the Facebook account, the investigator also found a conversation in which a Nebraska girl said she was 14 years old. The user asked for pictures and videos and the girl sent four nude images.
The fourth victim was a 16-year-old girl from Texas, who sent an explicit photo.
Junek told an investigator that he created a fake Facebook profile under the name that had communicated with the girls.
He said he realized he had a problem and had deactivated his Facebook page five or six times.
He did not admit to soliciting nude photographs from the girls.
Read or Share this story: https://post.cr/2GvF0dC
How is the legacy of war shaped?
Investigators: Restaurant fires caused by sushi ingredient
Snodgrass announces state Assembly campaign
Neenah Joint School District considers buying 225 acres, selling other property

redpajama_set_name: RedPajamaCommonCrawl


text:
\section{Introduction}

In the current cosmological paradigm, the Universe is composed of a filamentary network of structures shaped by gravity. In this framework, dark matter haloes correspond to overdense regions that evolve by gravitational instability due to mergers and interactions with other haloes. Galaxy formation occurs inside haloes where baryons collapse in the gravitational potentials and the condensation of cold gas allows the formation of stars and the evolution of galaxies \citep{White78}. A detailed description of the halo-galaxy connection enables using the galaxies to constrain the cosmological model.

The evolution of dark matter haloes can be followed, to high accuracy, using N-body simulations which use a set of cosmological parameters as inputs. In contrast, the evolution of galaxies in  haloes involves many physical processes that are still poorly understood. The fate of baryons within dark matter haloes has been modelled using different approaches, such as hydrodynamical simulations that provide an insight into the formation and evolution of galaxies \citep[e.g.][]{Vogelsberger14,Schaye15}. However, these models are computationally expensive and cannot be run over the large volumes needed for cosmological studies. Alternatively, the effect of baryons can be probed in such large volumes using semi-analytical models (SAMs) of galaxy formation. These start from haloes extracted from a large volume dark matter only simulation and use simplified physical models of the processes that shape the evolution of baryons \citep{cole2000,Baugh06-rv,Benson10,Sommerville15}. Hence, SAMs make predictions for the abundance and clustering of galaxies that can be compared and tested with large surveys.

Another way to describe the galaxy population is with the halo occupation distribution (HOD) framework \citep{Benson00, Peacokc00,Scoccimarro01,Yang03}. This is an empirical approach that provides a relation between the mass of haloes and the number of galaxies hosted by them. This is expressed as the probability distribution $P(N|M_h)$ that a halo of virial mass $M_h$ hosts $N$ galaxies which satisfy some selection criteria. This approach provides insight into the halo-galaxy connection and can be used to study galaxy clustering \citep{Berlind02,Zheng05,Conroy06,Zehavi11,Wechsler18}. Furthermore, the HOD parameters can be tuned in detail because they only aim to reproduce a limited set of observables such as the galaxy number density and clustering. Thus, HOD modelling is one of the most efficient ways to populate very large volumes or to produce many realizations required for, e.g., estimating covariance matrices using mock galaxy catalogues \citep[e.g.][]{Norberg:2009, Manera13}. These mock catalogues can then be used to test and develop new algorithms that will be used for the next generation of surveys.

The study of star forming emission line galaxies (ELGs) has gained interest over the last decade as they will be targeted by surveys such as Euclid and the Dark Energy Spectroscopic Instrument (DESI) surveys \citep{Laureijs:2011, DESI16}. The luminosity of an emission line depends on a number of factors, including the star formation rate (SFR), gas metallicity and the conditions in the HII regions \citep[e.g.][]{Orsi14}. Even though ELGs samples are related to star formation, they are not the same as SFR-selected samples. Still, a similar HOD approach can be used to study both galaxy populations \citep{Geach12, Cochrane17, Cochrane18}. In particular, the shape of the HOD in SFR selected samples is more complex than the case of the more widely studied stellar mass selected samples \citep[e.g.][]{contreras13,gp18}. For example, the occupation function of central galaxies in ELG samples does not follow the canonical step-like form. Accurate modelling of the HOD will provide the more realistic mock catalogues needed for the analysis of future observational samples.

Here, we use the HOD formalism to test three different ways to populate dark matter haloes with galaxies. The prescriptions of these models aim to replicate as accurately as possible the target galaxy populations of a SAM sample. The comparison between the galaxy population in the mock catalogues and SAM samples is done via the analysis of their two-point correlation function (2PCF), which is related to the power spectra of density fluctuations and is sensitive to cosmology \citep[e.g.][]{DeRose19}. We also include the scatter of the HOD of satellites in our modelling, and quantify the impact of using this additional parameter on the clustering.

The outline of this paper is as follows. The definition of galaxy samples used and the basic properties of the N-body simulation and the SAM are given in Section~\ref{sec: simulation data}. The correlation functions and the HODs of the samples are presented in Section~\ref{sec: Characterization of samples}. In Section~\ref{sec: building a mock} we introduce the HOD models used to build the mock catalogues and the recipes employed to perform this procedure. The main results and analysis are discussed in Section~\ref{sec: results} while in Section~\ref{sec: conclusions} we present our conclusions. Appendix~\ref{section:1-HODoccupationfunctions} shows the predicted occupation functions for a particular HOD model.
\section{Simulation data}
\label{sec: simulation data}

In this section we give a brief overview of the galaxy formation model used (\S~2.1) and the N-body simulation in which it is implemented (\S~2.2).

\subsection{Galaxy formation model}

A galaxy formation model needs to take into account a variety of physical processes such as radiative cooling of gas; AGN, supernovae and photoionisation feedback; chemical evolution; star formation; disc instabilities; collapse and merging of dark matter haloes; and galaxy mergers. These affect the fate of baryons in haloes which lead to the formation and evolution of galaxies. Several physical processes such as star formation and gas cooling are not fully understood due to their complexity. As a consequence, a set of free parameters are used in the equations that model these processes. These free parameters are tuned in order to reproduce observations such as the luminosity functions, colours and the distribution of morphological types. In this context, different SAMs usually have their own implementations to model these physical processes, predicting different galaxy populations. Here we use the outputs at $z=0$ from the SAM of \citet{Guo13} (hereafter G13) which is a version of the {\tt  L-GALAXIES} code from the Munich group \citep{DeLucia04, Croton06, DeLucia07, guo11, Henriques13}. The outputs are publicly available from the Millennium Archive\footnote{\url{http://gavo.mpa-garching.mpg.de/Millennium/}}.

The samples used here are defined according to three different number densities where we rank the galaxies in the SAM by their stellar mass or SFR in a decreasing way (hereafter the SAM samples). These samples are useful in order to compare with observational catalogues with similar space densities. Table~\ref{tab: density samples} shows the three number densities and the cuts in stellar mass and SFR used in each case.

\subsection{The Millennium simulation}
\label{sec:simulation data}

The distribution of dark matter haloes used in this work is drawn from the Millennium-WMAP7 simulation \citet{Guo13} which is identical with the Millennium Simulation \citet{Springel05}, but with updated cosmological parameters that match the results from the WMAP7 observations. This version assumes a flat $\Lambda$CDM universe considering $\Omega_m = 0.27$, $\Omega_{\Lambda} = 0.73$, $h = H_0/(\rm 100\ km\ s^{-1}\ Mpc^{-1}) = 0.704$ and $\sigma_{8} = 0.81$. The simulation was carried out in a box-size of $500\ h^{-1} \rm Mpc$ following $2160^3$ particles of mass $9.31\times 10^{8}h^{-1}\rm M_{\odot}$.  The run produced 61 simulation snapshots from $z=50$ up to $z=0$.  G13 use a friends-of-friends group finding algorithm ({\tt FOF}) to identify dark matter haloes in each snapshot \citep{Davis85} and then run {\tt SUBFIND} to identify the subhaloes \citep{Springel01}. Halo merger trees are constructed for each output and track the evolution of haloes through cosmic time. These trees are the starting points for the SAM.

\begin{table}
    \centering
    \begin{tabular}{ccc}
        \hline
         $n/h^{3}\rm Mpc^{-3}$ & $ {\rm M_{min}^*}/ h^{-1} {\rm M_{\odot}}$ & $\rm SFR_{min}/yr^{-1}M_{\odot}$ \\
        \hline
        $\rm 10^{-3.0}$ & $5.95\times 10^{10}$ & $5.25$ \\
        $\rm 10^{-2.5}$ & $3.38\times 10^{10}$ & $2.53$ \\
        $\rm 10^{-2.0}$ & $1.25\times 10^{10}$ & $0.70$\\
        \hline
    \end{tabular}
    \caption{The first column shows the abundance of galaxies in the three density samples used here. The second and third columns show the cuts applied to G13 galaxies in stellar mass and star formation rate, respectively, to achieve these abundances.}
    \label{tab: density samples}
\end{table}

\section{Characterization of the SAM galaxy samples}
\label{sec: Characterization of samples}

This section introduces the statistics used to characterize the distribution of galaxies, starting with the measurement of the correlation function (\S~3.1), the form of the HOD predicted by the SAM (\S~3.2) and the scatter in the HOD (\S~3.3).

\subsection{Clustering measurement: two-point galaxy correlation function}
\label{subsec: xi}

The spatial two-point correlation function, $\xi(r)$,  measures the excess probability of finding a pair of galaxies at a given separation with respect to a random distribution. We compute the 2PCF of the galaxy samples with the {\tt Corrfunc} code \citep{Sinha17}.

Fig.~\ref{fig: xi from SAMs} shows the 2PCF of the stellar mass (top) and SFR (bottom) selected samples for the different space densities. For the former, the amplitude of the clustering increases with decreasing number density, as we consider more massive galaxies. The impact of the inclusion of these massive galaxies is stronger at small scales and is weaker at large scales. In contrast, for the SFR selected galaxies the amplitude of the 2PCF for the different samples remains largely unchanged except for small scales where the satellite-satellite pairs make an important contribution to the clustering amplitude. For both selections, the satellite fraction increases with increasing number density.

In the 2PCF, we can distinguish between the contribution from galaxy pairs in the same halo and from different haloes. The former are the main contributors to the amplitude of the 2PCF on small scales, namely the one-halo term which dominates up to $\sim 1\ h^{-1} {\rm Mpc}$, while galaxy pairs between different haloes contribute mostly to the two-halo term which determines the clustering on large scales.
In this regime the total number of galaxies in the halo, regardless of whether they are satellites or the central,  drives the amplitude of the clustering, acting as a weighting for the bias of each halo in computing an overall ``effective'' bias for the sample (see e.g. \citealt{Baugh:1999}).
The one-halo term is sensitive not only to the number of satellites, but also depends on their spatial distribution.

\begin{figure}
    \centering
    \includegraphics[width=\columnwidth]{xi_SAM.pdf}
    \caption{Two-point correlation functions ($\xi(r)$) of different galaxy samples from G13 and defined in Table~\ref{tab: density samples}. ({\it Top}) stellar mass and SFR selected samples ({\it bottom}). Colors indicate each sample as labeled in the bottom panel. The fraction of satellites in each sample is shown in both panels, with the color indicating the  sample number density.}
    \label{fig: xi from SAMs}
\end{figure}

\subsection{The halo occupation function predicted by the SAMs}
\label{subsec: HODs from SAMs}

The galaxy populations in the SAMs depend on the choices adopted for the modelling of the baryonic processes. Hence, depending on the SAM employed, different galaxy catalogues with different luminosity functions, stellar mass functions or correlation functions can be obtained for the same dark matter simulation. For example, \citet{contreras13} studied the effects on the clustering predicted from different SAMs and found some differences particularly in galaxy samples selected by SFR and cold gas mass. Moreover, they show that the shapes of the HODs are model-dependent which reflects  the differences in the implementation of physical processes in each SAM. For example, the specific modelling of dynamical friction affects the satellite population in SAMs. Here we are not interested in the detailed shape of the HOD predicted by a particular SAM, but on how best to use the occupation functions to populate dark matter haloes with galaxies to produce a similar spatial distribution to that resulting from a SAM.

The HOD is usually broken down into the contribution from central and satellite galaxies. Fig.~\ref{fig: HODs from SAMs} shows these two components for stellar mass and SFR selected samples with the same number density for the G13 SAM. Here, each HOD is computed in bins of width 0.08 dex in the logarithm of the halo mass where the position of each $\langle N \rangle$ value is plotted at the median value within each bin. The striking difference in the shape of the HOD of centrals between the two selections is due to the different galaxies that are included. Massive centrals tend to be red galaxies hosted by massive dark haloes. Such centrals are included in stellar mass selected samples but not when selecting by SFR. The galaxies in the SFR samples correspond mainly to blue star-forming galaxies excluding luminous red galaxies with high stellar mass but low SFR. It is noteworthy that the fraction of halos that contain a central passing the SFR selection never reaches unity for the sample plotted in Fig.~\ref{fig: HODs from SAMs}. These features of  the HOD of SFR galaxies have been noted in SAMs before  \citep[e.g.][]{contreras13,contreras19,gp18} and inferred for blue galaxies in the SDSS  \citep[e.g.][]{Zehavi11}. This shows that a significant number of haloes in the SFR selected samples do not host a central as their SFR is below the threshold. The same situation is found in the other number density samples. Note that, in observational samples, the ranking of galaxies in order of their emission line luminosity may not correspond to the ranking in SFR due to dust attenuation, which means that the highest SFR galaxies may not necessarily have the brightest emission lines.

We estimate the uncertainties of the HOD values using jackknife resampling \citep{Norberg:2009}, dividing the simulation volume into 10 slices. We use the position of the centre of the potential of haloes to classify the galaxies within each halo. The resulting errors are shown as the shaded regions in Fig.~\ref{fig: HODs from SAMs}, and they are negligible for all halo masses except at the high mass end and for the HOD of centrals selected by SFR.

Because of the simple relation between halo mass and occupation number, the HOD represents a useful approach for the construction of mock galaxy catalogues. Here we have described the first moment of the HOD, the main ingredient to building-mocks recipes. Nevertheless, it is important to also consider the second moment, i.e the dispersion in the HOD of satellites.

\begin{figure}
    \centering
    \includegraphics[width=\columnwidth]{HODs_SAM.pdf}
    \caption{The HOD predicted by G13 for stellar mass ({\it top}) and SFR selected samples ({\it bottom}), for a number density of $10^{-2.5} h^3 \rm Mpc^{-3}$. Black lines show the HOD for the full sample and red and blue indicate the HOD for central and satellite galaxies, respectively.  The red and blue shaded regions represent jackknife errors calculated using 10 subsamples.
    The horizontal black dotted line shows an average occupation value of unity.}
    \label{fig: HODs from SAMs}
\end{figure}

\subsection{The predicted dispersion in the halo occupation number}
\label{subsec: HOD dispersion}

When the simplest HOD approach is used to build mock catalogues, the mean of the distribution is the main parameter. Central galaxies are assumed to follow a nearest-integer distribution where the mean $\rm \langle N_{cen} \rangle$ is between zero and one. For satellites, a Poisson distribution with mean $\rm \langle N_{sat} \rangle$ is the most widely assumed distribution \citep[e.g][]{Kravtsov04,Zheng05}.

In G13, satellites are classified as type-1 if they are hosted by a resolved subhalo, and type-2 or orphans if the subhalo has been destroyed by tidal effects and is no longer identified. \citet{BK10} found that the number of low mass subhaloes in main haloes in the Millennium-II Simulation \citep{BK09} is well described by a negative binomial distribution which corresponds to a super-Poissonian statistic as its scatter is larger than a Poisson distribution. This suggests that the type-1 satellite population can also be described by this distribution. Based on the outputs from the SAM presented in \citet{Jian16-I}, and using the Bolshoi \citep{Klypin11} and MultiDark \citep{Prada12} simulations, \citet{Jiang16} showed that ignoring this non-Poissonity in the HOD of subhaloes results in systematic errors in the predicted clustering of galaxies. Here we extend the application of the negative binomial distribution by checking whether the HOD of G13 satellites, that is including type-1 and type-2, is well described by this statistic. We expect that the HOD scatter is model-dependent because of the different treatments of dynamical friction. Moreover, as some galaxy properties, such as SFR and stellar mass, have a model-dependent scatter it is reasonable to assume the same for HODs. For example, \citet{Guo16} showed that different galaxy formation models do not have  the same dispersion in the stellar mass-halo mass relation. Therefore our results are specific to the G13 model. It is likely that a different SAM would require an adjustment to to the value of $\beta$ to describe the scatter of the satellite HOD. Nevertheless we expect our general results to hold for any SAM, and that the satellite distribution displays more scatter than Poisson.

The Poisson and negative binomial distributions differ in their shapes so it is useful to parametrise the departure from the Poisson scatter. We use the parameter $\beta$ (defined below) to denote this departure. For a Poisson distribution the variance is given by the mean value of the random variable, namely $\rm \langle N_{sat}\rangle$, with the standard deviation given by $\sigma = \sqrt{\rm \langle N_{sat}\rangle}$. The negative binomial distribution has the same mean as the Poisson distribution, but a larger scatter which can be expressed as

\begin{equation}
     \sigma_{\rm NB} = \sigma + \beta \sigma.
\end{equation}

\noindent where $0<\beta<1$. Then, $\beta$ indicates the fractional change in the variance with respect to the Poisson standard deviation $\sigma$. Under this definition, when $\beta=0$ the distribution is Poissonian and if $\beta =1$ the standard deviation is twice that from a Poisson distribution.

The probability function of the negative binomial distribution is given by

\begin{equation}
 P(N\left| r,p \right.) = \frac{\Gamma(N+r)}{\Gamma(r)\Gamma(N+1)}p^r (1-p)^N.
\end{equation}

\noindent Here $\Gamma(x) = (x+1)!$ is the gamma function. The parameters $r$ and $p$ are determined by the first moment $\left< N\right>$ and second moment $\sigma^2$ of the distribution,

\begin{equation} \label{eq. parameters of NB}
 p = \frac{\left< N\right>}{\sigma_{\rm NB}^2}, \hspace{2mm} r = \frac{\left< N\right>^2}{\sigma_{\rm NB}^2 - \left< N\right>}.
\end{equation}

Thus, we can control the width of the negative binomial distribution through the parameter $\beta$ and compute the value of $\sigma_{\rm NB}^2$.

\section{Generating HOD mock catalogues}
\label{sec: building a mock}

We now describe the procedure followed to build HOD mock galaxy catalogues using the HODs of the SAM samples.  Section~\ref{subsec: HOD models} presents the three methods we use to populate haloes with galaxies. In Section~{\ref{subsec: scatter in the HOD}} we specify the treatment of the scatter in the HOD of satellites. Section~{\ref{subsec: distribution of satellites}} explains how we impose a standard Navarro-Frenk-White (NFW) density profile for satellites. Section~\ref{subsec: AB and shuffle} presents the impact of assembly bias in the SAM samples and explains why it must be removed from the SAM in order to compare with the HOD mock catalogues. Finally in Section~4.5 we discuss the treatment of the radial distribution of satellite galaxies within haloes.

\subsection{The HOD models used to build mocks}
\label{subsec: HOD models}

We test three different HOD schemes of increasing complexity. This helps us to understand the level of complexity needed to obtain accurate clustering predictions. Each model uses occupation functions obtained from linear interpolations of the HOD values in each bin, rather than fitting a parametric form to the SAM HOD. The distribution of galaxies can be nearest-integer (centrals only) and Poisson or negative binomial (satellites).

\subsubsection{1-HOD}

The 1-HOD model builds mock catalogues using the HOD of all galaxies from the SAM sample (black solid lines in Fig.~\ref{fig: HODs from SAMs}) including both centrals and satellites. The model assumes either a Poisson or negative binomial distribution for the occupation number. We adopt a Monte Carlo approach to obtain the final number of galaxies.

This approach does not distinguish between centrals and satellites. If the model predicts that $\rm N \geq 1$ we assume that this halo hosts a central and $\rm N_{sat}=N-1$. Because of this, the number of centrals and satellites in the 1-HOD mock catalogues can be notably different with respect to the SAM samples where there are haloes with satellites but no central. Moreover, the HODs of these two separate components in the mock catalogues are completely different with respect to the HODs of the SAM samples (see Appendix~\ref{section:1-HODoccupationfunctions}). However, the total number of galaxies in these mock catalogues is essentially the same as in the SAM samples.

\subsubsection{2-HOD}

The 2-HOD model uses the HOD of centrals and satellites separately, i.e, the red and blue solid lines in Fig.~\ref{fig: HODs from SAMs}, respectively. Thus, a particular distribution can be assumed for each component and the modelling is done independently for each one. For centrals, we use the nearest-integer distribution and for satellites the Poisson or negative binomial distribution.

This scheme predicts practically the same number of central and satellites as the SAM samples. Note that in a non-negligible number of realizations it is possible to get haloes without a central. This is more likely for haloes with masses for which $\rm \left\langle N_{cen}\right\rangle <  1$ which is more frequently found in SFR-selected samples.

\subsubsection{4-HOD}

This model contains more information about the galaxy population of the SAM sample than the 2-HOD model. The 4-HOD requires us to store the number of haloes that host a central ($\rm N_{cen}$) and the number of haloes that do not host a central ($\rm N_{nocen}$) as a function of halo mass. Under this definition the total number of haloes in the volume is the sum of both quantities. Furthermore, the 4-HOD also needs knowledge of the number of satellites in haloes with a central ($\rm N_{sat\_cen}$) and without a central ($\rm N_{sat\_nocen}$). Thus, the total number of satellites is the sum of these two quantities. With these definitions, we build new HODs for satellites that take into account the population of centrals in the SAM samples. The SAM samples contain haloes with satellites but no centrals. This is more common in SFR selected samples. Indeed the HOD of centrals in these samples indicates that a large number of haloes do not host a central (see Fig~\ref{fig: HODs from SAMs}), and the 4-HOD takes this feature into account.

We then define the satellite occupation functions conditioned on whether or not haloes host a central. With the four quantities explained above, we can define the conditional HODs,

\begin{align}
\rm \langle N_{\rm sat\_cen}(M_h) \rangle  &= \rm \frac{N_{sat\_cen}}{N_{cen}}(M_h) \label{eq. sat_cen} \\
\rm \langle N_{\rm sat\_nocen}(M_h) \rangle  &= \rm \frac{N_{sat\_nocen}}{N_{nocen}}(M_h) \label{eq. sat_nocen}
\end{align}

Fig.~\ref{fig: 4-HOD} shows the conditional HODs where the main differences are observed at low halo masses. Even though the ratio between these two HODs is close to unity, it is the galaxies hosted by these haloes ($\approx 10^{12} h^{-1} {\rm M_{\odot}}$) that dominate the amplitude of clustering. The conditional HODs are well fitted by a negative binomial distribution, including the HODs of the other number density samples. The 4-HOD method uses a Monte Carlo approach to decide if a halo hosts a central galaxy. Depending on this outcome, one of the two conditional HODs is then chosen to obtain the number of satellites.

\begin{figure}
    \centering
    \includegraphics[width=\columnwidth]{4HODplot.pdf}
    \caption{Conditional HODs from the 4-HOD method for an SFR selected sample with number density ${\rm 10^{-2.5}} h^3 \,{\rm Mpc^{-3}}$. {\it Top:} Average number of satellites in haloes with a central (cyan) and without a central (magenta). {\it Bottom:} Ratio of the two HODs shown in the upper panel. Shaded regions represent jackknife errors calculated using 10 subsamples.}
    \label{fig: 4-HOD}
\end{figure}

\subsection{Treatment of scatter in the HOD of satellites}
\label{subsec: scatter in the HOD}

A Poisson distribution is fully described by its first moment. In the case of satellites this is $\rm \langle N_{sat} \rangle$.  If the distribution of the number of satellites follows instead a negative binomial distribution, an additional parameter $\beta$ is needed which specifies the increase in the scatter with respect to a Poisson distribution (see Eq.~1). We fix the $\beta$ value so that we reproduce as closely as possible the scatter of the HOD of satellites in a given SAM sample. Fig.~\ref{fig: beta definition} shows the scatter of the HOD of satellites in SAMs and 2-HOD mock catalogues for two illustrative $\beta$ values in a SFR selected sample. This shows that a small but non-zero $\beta$ is required to reproduce the HOD scatter of the SAM sample. The same is found for the other number density samples, and for the conditional HODs. We do not perform this analysis for the 1-HOD model as satellites are not treated independently in this case.

It is not possible to replicate the HOD scatter in SAMs more closely as this would require $\beta$ to be a function of mass. Instead, we assume a constant scatter for the HOD of satellites by using the same $\beta$ for all halo mass bins. The accuracy of the $\beta$ values used are judged by checking the quality of the resulting mocks via comparison of their 2PCFs with the clustering of the shuffled-NFW samples (see \S~4.5 below for the definition of this catalogue). We show in Section~\ref{sec: results}, that when the scatter of the SAM and HOD mocks are matched up to  $M_h \lesssim 10^{13.5} M_{\odot}h^{-1}$ ($\beta=0.05$ for SFR-selected samples), we obtain the most accurate clustering predictions. In contrast, using larger values for $\beta$ worsens the predictions (as does using $\beta=0$, which corresponds to Poisson scatter).

\begin{figure}
	\centering
	\includegraphics[width=\columnwidth]{HODs_scatter.pdf}
    \caption{The HOD of satellites in a SAM sample (dashed blue) and a 2-HOD mock catalogue (solid red) contrasting two values of the parameter $\beta$ that controls the scatter (see Eq.~1): $\beta=0.05$ ({\it top}) and $\beta=0.18$ ({\it bottom}). The shaded regions show the HOD scatter  and the red dotted lines correspond to the  scatter in the HOD mocks. The subpanels show the ratios between the HOD scatter of the mocks and SAM sample. Note that it is not possible to visually distinguish a Poisson scatter from the $\beta$ scaled versions plotted in the main panels, but this choice would lead to a larger ratio of variances than the range plotted in the lower subpanels.}
    \label{fig: beta definition}
\end{figure}

The satellite HOD is well described by the negative binomial distribution for a wide range of halo masses. Fig.~\ref{fig: histograms} shows the satellite PDF in a particular mass bin for a stellar mass and a SFR selected sample. We show negative binomial distributions defined by $\beta=0.08$ and $\beta=0.05$. In order to compute the satellite distributions, we split satellites according to whether or not their haloes host a central galaxy, which is relevant for the 4-HOD model. The satellite distribution matches with the negative binomial when most of haloes in the bin are included. A similar close match is found when comparing with Poisson distributions ($\beta=0$). Note that in the SFR selection case most of the haloes do not host a central galaxy, as the HOD of centrals in that bin suggests. The opposite behaviour is observed when selecting by stellar mass.

\begin{figure*}
  \begin{minipage}{\textwidth}
   \centering
      \includegraphics[width=0.8\textwidth]{Histograms.pdf}
      \caption{Probability distributions of satellites (cyan histograms) that are hosted by haloes with masses within the blue shaded (vertical) mass range  of the HODs shown in the insets. The galaxies are selected by stellar mass ({\it top}) and SFR ({\it bottom}), with a number density of ${\rm 10^{-2.5}} h^3 {\rm Mpc^{-3}}$. {\it Left:} Haloes without centrals (i.e its central did not enter to the cut). {\it Right} Same as left panel but considering haloes with centrals. Note the high probability to find haloes that do not host centrals in the SFR sample, which is expected by the low value of $\langle \rm  N_{cen} \rangle $ in the mass range analyzed and as is shown in the inset. The distributions are well described by negative binomial distributions (magenta) in the cases when most of haloes in the bin are included (i.e top right panel and bottom left panel). The negative binomial distributions shown here are obtained using a scatter that is $5\%$ larger than that from a Poisson distribution in the SFR selected sample, and $8\%$ larger in the stellar mass selection case}.
     \label{fig: histograms}
  \end{minipage}
\end{figure*}

\subsection{The radial distribution of satellite galaxies in halos}
\label{subsec: distribution of satellites}

The number of satellites in the mock catalogues is obtained from the adopted HOD model (see Sec~\ref{subsec: HOD models}). Their positions in haloes are set according the standard NFW density profile \citep{Navarro96} which requires two parameters, the concentration and scale radius. The former depends on halo mass and the latter is a function of the virial radius. For simplicity, we assume that all haloes in the simulation volume have the same concentration parameter $c=13.98$ which corresponds to the concentration of a halo at redshift $\rm z=0$ and mass $M_h = 10^{12.5}{\rm M_{\odot}}h^{-1}$. We do not use a more realistic model for concentration as we are interested in comparing the HOD models rather than obtain a realistic redistribution of satellites. We impose that the maximum distance from a satellite to the halo center is two virial radii which depends on the halo mass. This defines the NFW mass profile used to obtain the satellite distances by a Monte Carlo approach. We modify the SAM output to impose a similar satellite distribution as described below (\S~4.5).

\subsection{Removing assembly bias from the SAM output}
\label{subsec: AB and shuffle}

In order to determine the best methodology to produce HOD mock catalogues, we aim to compare them with the clustering of the original SAM samples via their 2PCF. Before making this comparison it is necessary to remove assembly bias from the SAM samples.

The clustering of dark matter haloes depends on additional properties besides mass. For example, \citet{Gao05} showed that the clustering of low mass haloes depends on their formation redshift and other works have found dependencies on concentration and subhalo occupation number \citep[e.g.][]{Wechsler06} among other secondary properties. This additional contribution to the clustering is commonly known as assembly bias and potentially changes the galaxy clustering amplitude on large scales.

The standard HOD approach considers only halo mass as the variable regulating the galaxy population. SAMs include assembly bias because they follow the evolution of baryons in halo merger histories that are shaped by the large-scale environment in the N-body simulation. Namely,  SAMs include a dependence on secondary halo properties as these affect the halo merger history and the evolution of galaxies that live within them. Thus, in order to compare the clustering between SAM samples and HOD mocks which use only halo mass as input, it is necessary to remove the assembly bias signal in the former samples.

Assembly bias can be eliminated from SAM samples through the shuffling technique introduced by \citet{Croton07}. This consists of randomly exchanging the galaxy populations between haloes of the same mass, thus removing any connection to the assembly history of the haloes. This procedure does not change the distances from satellites to their central galaxy in each halo. In clustering terms, the one-halo term of this ``shuffled'' catalogue is the same as the original SAM sample but its two-halo term is different because assembly bias is not present in the shuffled sample. If the SAM samples did not have assembly bias, we would measure the same 2PCFs for their shuffled samples as measured for the original output.

Fig.~\ref{fig: xi_SAM_Shuffle} shows the correlation functions of a SAM sample and its shuffled version, for both the stellar mass and the SFR selected samples. The assembly bias signature, shown in the middle panels, is evident in the clustering differences between these two catalogues at large separations. We also show the 2PCF of a modified shuffled sample that will be introduced below. The assembly bias signatures remain unchanged for the other samples, but they are noisier for the lowest number density samples as they contain fewer galaxies.

\begin{figure*}
 \begin{minipage}{\textwidth}
     \centering
     \includegraphics[width=0.95\textwidth]{xi_SAM_Shuffle.pdf}
    \caption{ {\it Top:} Correlation functions of SAM samples (dotted black in the main panel, solid black in the subpanel), and their modifications: the shuffled (dashed magenta) and the shuffled-NFW samples (dashed cyan). The galaxy samples are selected by stellar mass ({\it left}) and SFR ({\it right}), with a number density of ${\rm 10^{-2.5}} h^3{\rm Mpc^{-3}}$. {\it Middle:} Ratios between the 2PCF of the SAM and shuffled samples. Differences at large scales are signatures of assembly bias. {\it Bottom:} Ratios of the 2PCF with respect to the 2PCF from the shuffled-NFW samples. The differences in the one-halo term  below ${\rm 1\ Mpc/}h$ indicates the departure of the satellite profiles from a NFW.}
    \label{fig: xi_SAM_Shuffle}
\end{minipage}
\end{figure*}

It can be seen that assembly bias increases the clustering for stellar mass selected samples, as was shown by \citet{Zehavi18}. SFR selected samples, on the other hand, show a decreased clustering amplitude. For the intermediate galaxy density sample, the assembly bias enhances the two-halo term by $\sim 12\%$ for the stellar mass selected sample and suppresses the amplitude in the SFR selection case by $\sim 4\%$. The enhance of clustering amplitude for the other stellar mass selections remains similar. For the SFR selections, we observe that suppression on clustering becomes weaker for higher density samples. Indeed, assembly bias can enhance the amplitude if the density of the sample if very high,  as shown in \citep{contreras19}.

\subsection{The shuffled-NFW target catalogue: changing the satellite distribution in the SAM}

The radial profile of satellites in the G13 SAM deviates from the standard NFW profile of dark matter within halos because the SAM associates galaxies with subhalos (or a proxy, such as the most bound particle, in the case of subhalos which are no longer resolved). The radial profile of subhalos is different from that of the dark matter (see e.g. \citealt{Angulo:2009}). The choice of which subhalos (and former subhalos) are associated with galaxies is driven by the galaxy formation model, which determines the luminosity of any galaxy associated with a subhalo and whether or not it has merged due to dynamical friction (only type 2 satellites, those that no longer have a resolved subhalo associated with them, are considered as candidates for galaxy mergers).

The final step before testing the accuracy of the HOD models is to modify the shuffled SAM catalogue to force the satellites in each halo to follow an NFW profile.  We call the result the shuffled-NFW catalogue. Because satellite galaxies in the SAMs and shuffled samples do not follow an NFW profile, the one-halo term of their 2PCFs are different to the one-halo term of the shuffled-NFW sample, as shown in the bottom panel of Fig.~\ref{fig: xi_SAM_Shuffle}. The shuffled-NFW catalogue does not contain assembly bias, and the satellites follow the {\it same} NFW profile as adopted in the HOD mocks. We note that the shuffled-NFW is not intended to be the ``best'' prediction of galaxy clustering but rather is the target sample for the reconstructions using the HOD mocks which has a controlled 1-halo clustering pattern to facilitate testing.

\section{{\bf Testing the accuracy of the HOD models}}
\label{sec: results}

\subsection{Satellite radial distributions and clustering of HOD mocks}

In Fig.~\ref{fig: rad_dist} we show the satellite profiles in the SAM samples, for stellar mass and SFR selected samples separated into the contributions from type-1 satellites and from orphan galaxies. To examine the departure from NFW, we produce a SAM-NFW catalogue where satellites in SAM are forced to follow the same NFW profile as used in the HOD mocks. For this catalogue, we update the satellite positions in the SAM samples according the same NFW density profile used to produce the HOD mocks. Note that this SAM-NFW is different from the shuffled and shuffled-NFW catalogues mentioned above.

It can be seen that the NFW profile is different from  the profile of type-1 satellites and orphans, particularly for the SFR selected sample. The profiles in the 2-HOD and 4-HOD mock catalogues are also shown, and they match with the NFW as expected from the construction of the HOD mock. Both models also reproduce the NFW profile in the other number density samples.

The masses of host haloes of 1-HOD satellites do not correspond with the masses in the original SFR selected samples (see Fig.~\ref{fig: 1HOD results}). Thus, the virial radii of these haloes define NFW density profiles that are different from the profiles in the other models. This has an impact on the positions of satellites generating the striking difference with NFW in Fig.~\ref{fig: rad_dist} for the SFR selection. The same occurs for the stellar mass selections but it is less extreme than the SFR case.

\begin{figure}
    \centering
    \includegraphics[width=\columnwidth]{sat_profiles.pdf}
    \caption{Profile of satellites hosted by subhaloes (dashed magenta) and orphans (dashed cyan) in a stellar mass ({\it top}) and SFR selected sample ({\it bottom}). The lines show the SAM with an NFW imposed for all satellites (solid black), and HOD mocks built by the 1-HOD (solid red), 2-HOD (solid blue) and 4-HOD models (solid green) where the HODs of satellites are described by the Poisson distribution.}
    \label{fig: rad_dist}
\end{figure}

The HOD models predict different galaxy populations for the G13 SAM samples. Table~\ref{table: sat_fractions} shows the satellite fraction of the SAM samples and HOD mocks built by the three different models, assuming a Poisson distribution for the HOD of satellites. Note that the 2-HOD and 4-HOD predict almost the same satellite fraction as in the SAM samples because of the separate HOD modelling of centrals and satellites.

\begin{table}
\caption{Satellite fractions of the galaxy samples used. The first column indicates their number densities. Column 2,3,4 and 5 show the satellite fraction in the SAM samples and in the HOD mock built using the 1-HOD, 2-HOD and 4-HOD models, respectively.}
 \begin{tabular}{c c c c c }
  \hline
 & & Stellar Mass & &  \\
 \hline
 $n/h^3 {\rm Mpc^{-3}}$  & SAM & 1-HOD & 2-HOD & 4-HOD \\
  $10^{-3.0}$ & $0.280$ & $0.317$ & $0.279$ & $0.280$ \\
  $10^{-2.5}$ & $0.322$ & $0.405$  & $0.324$ & $0.324$ \\
  $10^{-2.0}$ & $0.381$ & $0.517$ & $0.383$ & $0.383$\\
 \hline
  & & SFR & &  \\
 \hline
 $n/h^3 {\rm Mpc^{-3}}$ & SAM & 1-HOD & 2-HOD & 4-HOD \\
  $10^{-3.0}$ & $0.171$ & $0.084$ & $0.175$ & $0.172$ \\
  $10^{-2.5}$ & $0.195$ & $0.192$ & $0.197$ & $0.197$ \\
  $10^{-2.0}$ & $0.244$ & $0.334$ & $0.246$ & $0.246$ \\
  \hline
  \end{tabular}

\label{table: sat_fractions}
\end{table}

We check the accuracy of each HOD model by comparing the HOD mocks with the shuffled-NFW sample via their 2PCFs. Fig.~\ref{fig: xi_d3} shows the clustering of the shuffled-NFW and the HOD mocks built using the HOD models described in Sec.~\ref{subsec: HOD models}. These particular models assume a Poisson distribution for the HOD of satellites. It can be seen that the three schemes produce accurate clustering predictions on large scales. On small scales, the 2-HOD and 4-HOD models produce similar accurate results while the 1-HOD shows striking differences. These deviations come from the overprediction of the number of satellites in the stellar mass selected samples. For the SFR selection cases, the difference is due to the notably different occupation function of central and satellites in the 1-HOD mock (see Appendix~\ref{section:1-HODoccupationfunctions}). The inaccuracy of the 1-HOD modelling is also present for the other number density samples too, whereas the 2-HOD and 4-HOD models produce similar quality results to the one shown here. As the 2-HOD and 4-HOD models are clearly the best, we drop the 1-HOD model henceforth.

\begin{figure}
    \centering
    \includegraphics[width=\columnwidth]{xi_d3.pdf}
    \caption{Clustering of HOD mock catalogues of stellar mass ({\it top}) and SFR selected samples ({\it bottom}) for a number density of ${\rm n=10^{-2.5}}/h^{-3}\rm Mpc^3$. The mocks are built using the 1-HOD (red), 2-HOD (blue) and 4-HOD (green) models, assuming a Poisson distribution for the HOD of satellites. The clustering of the shuffled catalogue with an NFW profile is shown as the dotted line. Subpanels show the ratios of the 2PCF of the mocks with respect to the 2PCF of the shuffle-NFW catalogue}
    \label{fig: xi_d3}
\end{figure}

\subsection{Impact of the assumed HOD scatter}

To study the impact of the scatter of the HOD on the clustering, we consider different dispersions for the negative binomial distribution in the construction of HOD mocks. Fig.~\ref{fig: xi_bump} shows the 2PCF of HOD mocks using different $\beta$ values. For the stellar mass selected samples, the scatter of the HOD does not have a significant impact.

For the SFR selection, the amplitude of the clustering on small scales is very sensitive to the scatter in the number of satellites. We find that increasing $\beta$ changes the amplitude of the one-halo term. When we split the contribution to the clustering from low and high mass haloes, we observe that the scatter mainly impacts the one-halo term of low mass haloes. This feature is present in HOD mocks built using the 2-HOD and 4-HOD methods. This is reproduced by both HOD models, indicating that this is a feature particular to the SFR selected samples.

\begin{figure*}
  \begin{minipage}{\textwidth}
      \centering
      \includegraphics[width=0.9\textwidth]{xi_bump.pdf}
      \caption{2PCF of HOD mocks built using the 2-HOD ({\it right}) and 4-HOD ({\it left}) models, for stellar mass ({\it top}) and SFR selected samples ({\it bottom}). The HOD models are used to build mock catalogues assuming a Poisson distribution (solid red) and negative binomial distributions of $\beta=0.05$ (solid blue) and $\beta=0.18$ (solid green) for the HOD of satellites. The clustering of the shuffled-NFW samples (dotted black) is shown in each panel. The ratios between the clustering of HOF mocks and shuffled-NFW are shown in the subpanels.}
      \label{fig: xi_bump}
  \end{minipage}
\end{figure*}

The most accurate clustering reconstructions for the G13 samples are obtained when we use the 2-HOD or 4-HOD to build mock catalogues assuming a negative binomial distribution for the HOD of satellites. Note that clustering predictions from both models do not show significant differences.

Fig~\ref{fig: 4HOD_results} shows the particular results from the 4-HOD modelling for all the space density samples. For the G13 SFR selected samples, the 4-HOD (and the 2-HOD) modelling produces the best results when $\beta=0.05$, which corresponds to a distribution slightly wider than a Poisson distribution. For the case of stellar mass selected samples, the best reproduction is obtained with $\beta=0.08$. Using instead the Poisson distribution (i.e, $\beta=0$) produces worse results for both selections particularly in the one halo regime.

For SFR selections, when using $\beta=0.05$ and $\beta=0$, the departures from the shuffled-NFW catalogues are below $\sim 8\%$ and $\sim 15\%$, respectively. It can be seen that the dispersion of the 2PCFs becomes important in the lowest number density sample. However, the assumption of the negative binomial distribution still produces better results, especially in the transition from the one- to the two-halo term.

For the stellar mass selection cases, the impact on clustering when using different $\beta$ values is much less significant. Indeed, the weak relation between clustering and HOD scatter, shown in Fig.~\ref{fig: xi_bump}, suggests that it is not necessary to include additional scatter in the construction of HOD mock catalogues for stellar mass selections. To compare with the Poisson distribution, we show also the clustering prediction for the stellar mass samples using $\beta=0.08$, for which we obtain the best result.

\begin{figure*}
  \begin{minipage}{\textwidth}
      \centering
      \includegraphics[width=0.9\textwidth]{4HOD_results.pdf}
      \caption{Ratios between the 2PCFs of mock catalogues, constructed with the 4-HOD method, and the shuffled-NFW catalogue. The HOD of satellites in the HOD mocks follows either the negative binomial (green) or Poisson distributions (cyan), with the color indicating the value of $\beta$. We show results for stellar mass ({\it top}) and SFR selected samples ({\it bottom}). Number densities increase from left to the right as labelled at the top of each column. The shaded regions represent jackknife errors calculated  using  10 subsamples.}
     \label{fig: 4HOD_results}
  \end{minipage}
\end{figure*}

Satellites in the G13 SAM are well described by a non-Poisson distribution. This is consistent with the HOD of subhaloes found in \citet{BK10}. The recipes that build mock catalogues of SFR selected samples using the HOD approach must undertake an analysis of the scatter of the HOD of satellites as it impacts the clustering. This analysis will provide the best $\beta$ to construct a HOD mock of a particular sample. For stellar mass samples, the HOD scatter has a weak impact on clustering, so the same analysis is not necessary in the context of HOD mock catalogues.

\section{Conclusion}
\label{sec: conclusions}

The next generation of surveys will measure the clustering of the galaxy distribution over a wide range of redshifts. Mock catalogues have proven to be important tools in preparation for this because of their multiple applications including error analysis, data interpretation and survey planning. SAMs are a physical approach to obtain such mocks, but
sometimes their direct application to a simulation is not possible due to the limited resolution of the halo merger trees \citep{Angulo:2014} or the trees may not be available, as in the case of the Euclid flagship simulation \citep{Potter:2017}. Even if the trees were available at the required resolution, the sheer number of halos in a giga-parsec side N-body simulation may preclude a direct calculation with a SAM.

The HOD model provides a simple yet efficient way to construct mock catalogues. The modelling consists of using a probability distribution to obtain the number of galaxies hosted in a halo of a particular mass. This simple method allows us to create large sets of mock catalogues for huge cosmological simulations. This is useful in the context of ELGs, as they are targets in current and coming surveys.

To determine the level of complexity needed to produce accurate mock catalogues, we test different HOD models. The 1-HOD uses the HOD of all galaxies making no distinction between centrals and satellites while the 2-HOD uses the HOD of these two components separately. The 4-HOD stores additional information about whether or not haloes host a central, and it constructs conditional HODs for satellites taking into account this information.

Because SAMs include assembly bias by construction, and in their simplest form HOD mocks do not, we remove the assembly bias from the G13 SAM samples by shuffling the galaxy populations among haloes of the same mass creating the shuffled catalogue. This allows us to make a direct comparison between the clustering of our mocks and the SAMs from which we extract the HOD measurements. For example, We find that, for the intermediate galaxy density sample in the G13 SAM, the assembly bias affects the 2-halo term of the 2PCF of stellar mass selected galaxies increasing the amplitude by $\sim 12\%$. For the SFR selected galaxies, in contrast, the assembly bias suppresses the clustering by  $\sim 4\%$.  We also impose the standard NFW profile for satellites in the shuffled catalogue as is done for satellites in the HOD mock catalogues. We can then check the accuracy of the HOD models through a comparison between the 2PCFs of the HOD mocks and the shuffled-NFW catalogues.

The 2-HOD and 4-HOD produce the best mock catalogues as their 2PCFs are in close agreement with the clustering of the shuffled-NFW sample. We obtain the best results using a negative binomial distribution for the (conditional) HOD (see Eq.~1); in previous works this was commonly considered to be a Poisson distribution. This is consistent with the subhalo HOD found in \citet{BK10} using the Millennium-II simulation. Furthermore, we found that the assumption of this non-Poisonnian HOD changes the galaxy clustering. Previously, \citet{Jiang16} found a similar result using subhaloes from the Bolshoi, the MultiDark N-body simulations, and the SAM presented in \citet{Jian16-I}.

The scatter of the HOD of satellites in G13 is reproduced by  a negative binomial distribution up to halo masses of ${\rm M_h \lesssim 10^{13.5}M_{\odot}}h^{-1}$. The galaxies in this halo mass range dominate the amplitude of the 2PCF.
We quantify the departure from the Poisson distribution with the parameter $\beta$ (see Eq. \ref{eq. parameters of NB}). We obtain the best clustering predictions for SFR selected samples using $\beta=0.05$ and $\beta = 0.08$ for stellar mass selected samples. These correspond to negative binomial distributions slightly wider than Poisson. Because of the specific modelling of different SAMs, we expect that the best $\beta$ values for each sample are model-dependent. For stellar mass selected samples, we find that the HOD scatter has a weak impact on clustering, making the addition of this additional parameter unnecessary in the context of mock catalogues.

The analysis of the HOD of satellites is important because the width of the distribution (determined by the $\beta$ parameter) has a large impact on the one-halo term of the 2PCF of mock catalogues that emulate SFR-selected sample and ELG samples. If we consider the Poisson distribution for the HOD of satellites ($\beta=0$) the 2PCF of the mock catalogues is underestimated with respect to the clustering of the shuffled-NFW. In contrast, using the negative binomial distribution increases the amplitude of clustering in the one-halo regime. If we assume a value of $\beta$ larger than the one present in the distribution of number of satellites, the clustering on small scales is further overestimated. We highlight the importance to perform a careful analysis of the satellite HOD if the HOD framework is used to produce mock catalogues, for ELGs or star forming galaxies, following a particular model or observation.

\section*{Acknowledgements}

This work was made possible by the efforts of Gerard Lemson and colleagues at the German Astronomical Virtual Observatory in setting up the Millennium Simulation database in Garching.
EJ, SC, NP \& IZ acknowledge the hospitality of the ICC at Durham University.
EJ acknowledges support from ``Centro de Astronom\'{i}a y Tecnolog\'{i}as Afines'' BASAL 170002.
NP acknowledges support from Fondecyt Regular 1191813.
IZ acknowledges support by NSF grant AST-1612085.
This project has received funding from the European Union's Horizon 2020 Research and Innovation Programme
under the Marie Sk\l{}odowska-Curie grant agreement No 734374. The calculations for this paper were performed
on the Geryon computer at the Center for Astro-Engineering UC, part of the BASAL PFB-06, which received additional
funding from QUIMAL 130008 and Fondequip AIC-57 for upgrades.

redpajama_set_name: RedPajamaArXiv


text:
Supermarket in Town Centre, Stockton-On-Tees, County Durham
Displaying 1 - 30 of 30 Map: Show Map
Darlington Lane, Town Centre, Stockton-On-Tees
24 Wellington Square, Town Centre, Stockton-On-Tees
John Walker Square, Town Centre, Stockton-On-Tees
Heron Frozen Foods
Castle Way, Town Centre, Stockton-On-Tees
5 Bath Lane, Town Centre, Stockton-On-Tees
Portrack Lane, Town Centre, Stockton-On-Tees
Urlay Nook Road, Town Centre, Stockton-On-Tees
Durham Road, Town Centre, Stockton-On-Tees
The Glebe, Town Centre, Stockton-On-Tees
Sainsbury's Local
1-3 Marske Parade, Town Centre, Stockton-On-Tees
Norton Road, Town Centre, Stockton-On-Tees
The Co-operative Food
24-28 High Street, Town Centre, Stockton-On-Tees
11 Birkdale Road, Town Centre, Stockton-On-Tees
Tesco Express
6 Premier Parade, Town Centre, Stockton-On-Tees
1 Sunningdale Drive, Town Centre, Stockton-On-Tees
2-6 Durham Lane, Town Centre, Stockton-On-Tees
63 Redhill Road, Town Centre, Stockton-On-Tees
Bishopton Road, Town Centre, Stockton-On-Tees
10 Butterwick Road, Town Centre, Stockton-On-Tees
10 Harper Parade, Town Centre, Stockton-On-Tees
10 Harper Terrace, Town Centre, Stockton-On-Tees
1 White House Drive, Town Centre, Stockton-On-Tees
Front Street, Town Centre, Stockton-On-Tees
318 Bishopton Road, Town Centre, Stockton-On-Tees
Pearson Way, Teesside Park, Stockton-On-Tees
59 Westbury Street, Teesside Park, Stockton-On-Tees
29-31 Havelock Street, Teesside Park, Stockton-On-Tees
Newmarket Avenue, Teesside Park, Stockton-On-Tees
Teesside Park
Coulby Newham


Narrative documents are typically a narrative type. Narrative essay matters aren't easy and it's a really depth focused and unbelievably particular style of creating. Following are a few really great bigotry article subjects that will definitely allow you to come up along with a very good written piece. This really is because a literacy story isn't just about authoring your own favourite literature, it's more private than that. This is understood as free wheeling composition creating. Narrative composition, in precisely the same manner as some other class article, ought to be restricted in range. The straightforward belief here will be to write an essay that's important. It essay writing service certainly is great to select the perfect narrative essay topic for your own function. Narrative creating is really a kind of creating that doesn't require the writer to run any type of background study.
With no community present, jack unleashes his inner animal and becomes dangerously violent.


text:
A Case of Anti-Caspr2 Autoimmune Encephalitis Associated with Adenocarcinoma of the Lung
Fábio Pé D'Arca Barbosa1, Francisco Vara-Luiz2, Ana Albuquerque1, Ana Valada Marques1, Vanda Spencer1
1 Internal Medicine Department, Hospital Garcia de Orta, Almada, Portugal
2 Gastroenterology Department, Hospital Garcia de Orta, Almada, Portugal
How to cite this article: Pé D'Arca Barbosa F, Vara-Luiz F, Albuquerque A, Valada Marques A, Spencer V. A case of anti-Caspr2 autoimmune encephalitis associated with adenocarcinoma of the lung. EJCRIM 2022;9 doi:10.12890/2022_003544.
Patient consent: Informed consent was obtained from the patient for publication of this case report
Autoimmune encephalitis (AIE) is an inflammatory brain condition with multiple aetiologies but is mainly associated with paraneoplastic syndromes. Several antibodies described in AIE are being investigated in relation to different cancers, including antibodies against Contactin-associated protein-like 2 (Caspr2), which have been associated with thymoma but very rarely with lung cancer. The authors present the case of a 72-year-old man with cardiovascular risk factors, who presented with a 3-week history of left hemichorea following a first unprovoked seizure the week before, with no other signs or symptoms. The patient was submitted to extensive aetiological investigation, testing positive for anti-Caspr2 antibodies in the cerebrospinal fluid. AIE was diagnosed after other possible causes had been excluded. A PET scan showed signs of abnormal metabolism, with a lung biopsy confirming lung adenocarcinoma. This case highlights a very rare association and the importance of a thorough aetiological investigation for neurological complaints.
Neurological complaints require a thorough aetiological investigation, and positive findings on imaging studies must be carefully examined.
Anti-Caspr2 antibodies are classically associated with Morvan syndrome, but other neurological presentations must also be considered.
Anti-Caspr2 autoimmune encephalitis may be the primary presentation of lung adenocarcinoma.
Autoimmune encephalitis, paraneoplastic syndrome, anti-Caspr2, lung adenocarcinoma
The presence of autoantibodies to Contactin-associated protein-like 2 (Caspr2)—a cell adhesion molecule—can manifest with central and peripheral nervous system symptoms [1]. Depending on the region involved, there can be deficits in motor control, sensation, behaviour, memory, awareness, and cognitive and autonomic function.
The characteristic constellation of symptoms consists of neuromyotonia, encephalitis, dysautonomia and insomnia (Morvan syndrome), but there may be other manifestations [1]. Throughout the course of the disease, cognitive disturbance and epilepsy are present in up to 80% and 50% of patients, respectively [2]. Anti-Caspr2 antibodies are mostly associated with thymoma, but are very rarely associated with other solid tumours, namely lung adenocarcinoma [2–4].
The diagnosis of autoimmune encephalitis (AIE) is challenging and requires the exclusion of other causes. The presence of anti-neuronal antibodies aids diagnosis and is crucial for determining management and guiding aetiological investigation [5]. Treatment consists of immunosuppression and, when indicated, tumour resection. Without treatment, neurological deterioration and death can occur [6].
A 72-year-old man with a medical history of hypertension, dyslipidaemia, acute myocardial infarction and transitory ischaemic stroke, was admitted to the emergency department with a 3-week history of left hemichorea, preceded a week earlier by a first episode of nocturnal generalized seizure, with no other signs or symptoms.
Computed tomography and magnetic resonance imaging revealed the presence of chronic vascular lesions in the basal ganglia bilaterally, as well as signs of chronic microangiopathic leukoencephalopathy with findings of T2 hyperintensities around the ventricles, in subcortical areas, basal ganglia and brainstem, compatible with enlarged peri-vascular spaces considering attenuation in T2-FLAIR. The electroencephalogram showed slow and abrupt waves in the right temporal lobe. In light of these findings, the diagnosis of epilepsy was assumed and the patient was started on levetiracetam. He was then submitted to a vascular study of the neck, intracranial and supra-aortic vessels and cardiac imaging, with no relevant abnormalities being detected. Blood tests revealed secondary hyperparathyroidism due to vitamin D deficiency and positivity for antiphospholipid antibodies, without other serum autoimmune markers, nutritional imbalances, organ dysfunction or electrolyte disturbances detected. Treatment for vitamin D deficiency was started promptly and the patient was put on anticoagulation with warfarin.
A lumbar puncture was performed and anti-Caspr2 antibodies were seen in both cerebrospinal fluid and serum, confirming the diagnosis of AIE. An investigation for occult primary neoplasm was carried out with a whole-body CT scan and positron emission tomography scan, which showed an area of hypermetabolism within the lung. The patient was then submitted to lung biopsy which confirmed the diagnosis of lung adenocarcinoma. He was treated with chemotherapy and immunotherapy.
AIE is an inflammatory brain disease characterized by compromised brain function due to antibodies directed against neuronal cell-surface proteins, ion channels or receptors [7]. Depending on the brain area affected, patients can present with a wide variety of symptoms. Anti-Caspr2 antibody-mediated AIE usually evolves over a period of 4 months, with some cases developing over a year [2]. Our patient's symptoms, which are present in up to a third of cases at first evaluation, prompted several follow-up investigations early in the course of the disease, allowing for a diagnosis approximately 1 month after symptom onset. During the investigation, a diagnosis of idiopathic antiphospholipid syndrome (aPL) was made according to the Sydney Classification Criteria, in light of the presence of brain vascular lesions associated with positive antiphospholipid antibodies. Several neurological disorders have been reported in aPL-positive patients, including epilepsy, psychosis, chorea, hemiballismus, transverse myelopathy, sensorineural hearing loss and migraine. However, a strong association has not been established between these manifestations and aPL [8]. The neuroimaging abnormalities in our patient were not typical for AIE and because the presence of anti-neuronal antibodies alone is not diagnostic of AIE, the positive finding of anti-Caspr2 antibodies was only taken into account because the distribution of the vascular lesions did not fully explain the clinical presentation and because other possible causes had been ruled out, namely neurologically relevant vitamin deficiencies, toxic/metabolic disturbances, infections and other autoimmune and neurodegenerative disorders.
A wide panel of anti-neuronal antibodies should be used in patients with symptoms consistent with AIE which has not yet been definitively diagnosed. Correlation between AIE and tumours is variable [9]. In patients with lung cancer, most neurological paraneoplastic manifestations occur with small-cell carcinoma, with very few cases reporting an association with lung adenocarcinoma [4, 9]. The diagnosis of lung adenocarcinoma in this patient allowed the underlying condition to be treated, thus improving prognosis. This clinical report highlights the association between non-small-cancer carcinomas and lesser-known antibodies, namely anti-Caspr2, as well as the importance of an extensive search for tumours in patients with antibody-mediated AIE. In view of the clinical variables of this case and the few reports associating anti-Caspr2 with lung adenocarcinoma, further studies are warranted to better clarify these findings.
Patterson K, Dalmau J, Lancaster E. Mechanisms of Caspr2 antibodies in autoimmune encephalitis and neuromyotonia. Ann Neurol 2018;83(1):40–51.
van Sonderen A, Ariño H, Petit-Pedrol M, Leypoldt F, Körtvélyessy P, Wandinger K, et al. The clinical spectrum of Caspr2 antibody-associated disease. Neurology 2016;87(5):521–528.
Qin X, Yang H, Zhu F, Wang Q, Shan W. Clinical character of CASPR2 autoimmune encephalitis: a multiple center retrospective study. Front Immunol 2021;12:652864.
Boyko M, Au K, Casault C, de Robles P, Pfeffer G. Systematic review of the clinical spectrum of CASPR2 antibody syndrome. J Neurol 2020;267(4):1137–1146.
Graus F, Vogrig A, Muñiz-Castrillo S, Antoine J, Desestret V, Dubey D, et al. Updated diagnostic criteria for paraneoplastic neurologic syndromes. Neurol Neuroimmunol Neuroinflamm 2021;8(4):e1014.
Abboud H, Probasco J, Irani S, Ances B, Benavides D, Bradshaw M, et al. Autoimmune encephalitis: proposed best practice recommendations for diagnosis and acute management. J Neurol Neurosurg Psychiatry 2021;92(7):757–768.
Dalmau J, Geis C, Graus F. Autoantibodies to synaptic receptors and neuronal cell surface proteins in autoimmune diseases of the central nervous system. Physiol Rev 2017;97(2):839–887.
Cervera R, Piette J, Font J, Khamashta M, Shoenfeld Y, Camps M, et al. Antiphospholipid syndrome: clinical and immunologic manifestations and patterns of disease expression in a cohort of 1,000 patients. Arthritis Rheum 2002;46(4):1019–1027.
Seluk L, Taliansky A, Yonath H, Gilburd B, Amital H, Shoenfeld Y, et al. A large screen for paraneoplastic neurological autoantibodies; diagnosis and predictive values. Clin Immunol 2019;199:29–36.