A conceptual model of the Plio-Pleistocene glacial cycles is developed based on the Budyko-Sellers type energy balance model. The model is shown to admit a phenomenon like the Mid-Pleistocene transition, capturing the essence of the albedo and the temperature precipitation feedbacks.

In the standard model (SM), lepton flavor violating (LFV) Higgs decay is absent at renormalizable level and thus it is a good probe to new physics. In this article we study a type of new physics that could lead to large LFV Higgs decay, i.e., a lepton-flavored dark matter (DM) model which is specified by a Majorana DM and scalar lepton mediators. Different from other similar models with similar setup, we introduce both left-handed and right-handed scalar leptons. They allow large LFV Higgs decay and thus may explain the tentative Br$(h\ra\tau\mu)\sim1\%$ experimental results from the LHC. In particular, we find that the stringent bound from $\tau\ra\mu\gamma$ can be naturally evaded. One reason, among others, is a large chirality violation in the mediator sector. Aspects of relic density and especially radiative direct detection of the leptonic DM are also investigated, stressing the difference from previous lepton-flavored DM models.

We present a scalar dark matter (DM) model where DM ($X_I$) is stabilized by a local $Z_2$ symmetry originating from a spontaneously broken local dark $U(1)_X$. Compared with the usual scalar DM with a global $Z_2$ symmetry, the local $Z_2$ model possesses three new extra fields, dark photon $Z^{'}$, dark Higgs $\phi$ and the excited partner of scalar DM ($X_R$), with the kinetic mixing and Higgs portal interactions dictated by local dark gauge invariance. The resulting model can accommodate thermal relic density of scalar DM without conflict with the invisible Higgs branching ratio and the bounds from DM direct detections, thanks to the newly opened channels, $X_I X_I \rightarrow Z^{'} Z^{'}, \phi\phi$. In particular, due to the new particles, the ${\rm GeV}$ scale $\gamma$-ray excess from the Galactic Center (GC) can be originated from the decay of dark Higgs boson which is produced in DM annihilations.

We show that the 't Hooft-Polyakov monopole model in the hidden sector with Higgs portal interaction makes a viable dark matter model, where monopole and massive vector dark matter (VDM) are stable due to topological conservation and the unbroken subgroup $U(1)_X$. We show that, even though observed CMB data requires the dark gauge coupling to be quite small, a right amount of VDM thermal relic can be obtained via $s$-channel resonant annihilation for the mass of VDM close to or smaller than the half of SM higgs mass, thanks to Higgs portal interaction. Monopole relic density turns out to be several orders of magnitude smaller than observed dark matter relic density. Direct detection experiments, particularly, the projected XENON1T experiment, may probe the parameter space where the dark Higgs is lighter than $\lesssim 60 {\rm GeV}$. In addition, the dark photon associated with unbroken $U(1)_X$ contributes to the radiation energy density at present, giving $\Delta N_{\rm eff}^\nu \sim 0.1$ as the extra relativistic neutrino species.

The correlation between the invisible Higgs branching ratio ($B_h^{\rm inv} $) vs. dark matter (DM) direct detection ($\sigma_p^{\rm SI}$) in Higgs portal DM models is usually presented in the effective field theory (EFT) framework. This is fine for singlet scalar DM, but not in the singlet fermion DM (SFDM) or vector DM (VDM) models. In this paper, we derive the explicit expressions for this correlation within UV completions of SFDM and VDM models with Higgs portals, and discuss the limitation of the EFT approach. We show that there are at least two additional hidden parameter in $\sigma_p^{\rm SI}$ in the UV completions: the singlet-like scalar mass $m_2$ and its mixing angle $\alpha$ with the SM Higgs boson ($h$). In particular, if the singlet-like scalar is lighter than the SM Higgs boson ($m_2 < m_h \cos \alpha / \sqrt{1 + \cos^2 \alpha}$), the collider bound becomes weaker than the one based on EFT.

Using a suite of detailed numerical simulations we estimate the level of anisotropy generated by the time evolution along the light cone of the 21cm signal from the epoch of reionization. Our simulations include the physics necessary to model the signal during both the late emission regime and the early absorption regime, namely X-ray and Lyman-band 3D radiative transfer in addition to the usual dynamics and ionizing UV transfer. The signal is analysed using correlation functions perpendicular and parallel to the line of sight (LOS). We reproduce general findings from previous theoretical studies: the overall amplitude of the correlations and the fact that the light cone anisotropy is visible only on large scales (100 cMpc). However, the detailed behaviour is different. At 3 different epochs, the amplitude of the correlations along and perpendicular to the LOS differ from each other, indicating anisotropy. These 3 epochs are associated with 3 events of the global reionization history: the overlap of ionized bubbles, the onset of mild heating by X-rays in regions around the sources, and the onset of efficient Lyman-alpha coupling in regions around the sources. A 20x20 deg^2 survey area may be necessary to mitigate sample variance when we use the directional correlation functions. On a 100 cMpc scale the light cone anisotropy dominates over the anisotropy generated by peculiar velocity gradients computed in the linear regime. By modelling instrumental noise and limited resolution, we find that the anisotropy should be easily detectable by the SKA, assuming perfect foreground removal, the limiting factor being a large enough survey size. In the case of the LOFAR, it is likely that only first anisotropy episode will fall in the observing frequency range and will be detectable only if sample variance is much reduced (i.e. a larger than 20x20 deg^2 survey, which is not presently planned).

Assuming dark matter is absolutely stable due to unbroken dark gauge symmetry and singlet operators are portals to the dark sector, we present a simple extension of the standard seesaw model that can accommodate all the cosmological observations as well as terrestrial experiments available as of now, including leptogenesis, extra dark radiation of $\sim 0.08$ (resulting in $N_{\rm eff} = 3.130$ the effective number of neutrino species), Higgs inflation, small and large scale structure formation, and current relic density of scalar DM ($X$). The Higgs signal strength is equal to one as in the SM for unbroken $U(1)_X$ case with a scalar dark matter, but it could be less than one independent of decay channels if the dark matter is a dark sector fermion or if $U(1)_X$ is spontaneously broken, because of a mixing with a new neutral scalar boson in the models.

Assuming ({\it i}) the seesaw mechanism works for neutrino masses and mixings, ({\it ii}) dark matter is absolutely stable due to unbroken $U(1)_X$ dark gauge symmetry, and ({\it iii}) the singlet fields ($H^\dagger H$ and $N_R$) are portal to the dark sector, we construct a simple model which is consistent with all the cosmological observations as well as terrestrial experiments available as of now, including leptogenesis, extra dark radiation of $\sim 8 %$ (resulting in $N_{\rm eff} = 3.130$ the effective number of neutrino species), Higgs inflation, small and large scale structure formation, and current relic density of scalar dark matter ($X$). The electroweak vacuum of this model is stable up to Planck scale for $m_H = 125$ GeV without any other new physics. The Higgs signal strength is equal to one as in the standard model for unbroken $U(1)_X$ case with a scalar dark matter, but it could be less than one independent of decay channels if the dark matter is a dark sector fermion or if $U(1)_X$ is spontaneously broken, because of a mixing with a new neutral scalar boson in the models. Detailed study of Higgs properties at the LHC would shed light on the models described in this work.

By combining high resolution, radiative transfer cosmological simulations of z~6 galaxies with a sub-grid multi-phase model of their interstellar medium we derive the expected intensity of several far infrared (FIR) emission lines ([C II] 158 micron, [O I] 63 micron, and [N II] 122 micron) for different values of the gas metallicity, Z. For Z = Z_sun the [C II] spectrum is very complex due to the presence of several emitting clumps of individual size < 3 kpc; the peak is displaced from the galaxy center by ~100 km/s. While the [O I] spectrum is also similarly displaced, the [N II] line comes predominantly from the central ionized regions of the galaxy. When integrated over ~500 km/s, the [C II] line flux is 185 mJy km/s; 95% of such flux originates from the cold (T~250 K) H I phase, and only 5% from the warm (T~5000 K) neutral medium. The [O I] and [N II] fluxes are ~6 and ~90 times lower than the [C II] one, respectively. By comparing our results with observations of Himiko, the most extended and luminous Lyman Alpha Emitter (LAE) at z = 6.6, we find that the gas metallicity in this source must be sub-solar. We conclude that the [C II] line from z~6 galaxies is detectable by the ALMA full array in 1.9 < t_ON < 7.7 hr observing time, depending on Z.

Lyman Alpha Emitters (LAEs) are usually thought to be powered by star formation. It has been recently reported that a fraction of LAEs at redshift z~3-4 hosts an Active Galactic Nuclei (AGN). If an AGN is present it could be obscured and undetectable in X-rays, but yet dominate the Lya luminosity. We examine the properties of these AGN-powered LAEs at high redshift (z>=6) using radiative transfer cosmological simulations and obtain a reliable criterion to identify them from their observed \Lya line and surface brightness. We find that these sources should have: (a) negative line weighted skewness, S_w <0, and (b) surface brightness profiles FWHM >= 1.5". This parameter space cannot be populated by starburst LAEs. Thus, LAEs satisfying this criterion would be strong candidates for the presence of a hidden AGN powering their luminosity.

We present a novel method to investigate cosmic reionization, using joint spectral information on high redshift Lyman Alpha Emitters (LAE) and quasars (QSOs). Although LAEs have been proposed as reionization probes, their use is hampered by the fact their Ly{\alpha} line is damped not only by intergalactic HI but also internally by dust. Our method allows to overcome such degeneracy. First, we carefully calibrate a reionization simulation with QSO absorption line experiments. Then we identify LAEs in two simulation boxes at z=5.7 and z=6.6 and we build synthetic images/spectra of a prototypical LAE. At redshift 5.7, we find that the Ly{\alpha} transmissivity (T_LAE) ~ 0.25, almost independent of the halo mass. This constancy arises from the conspiracy of two effects: (i) the intrinsic Ly{\alpha} line width and (ii) the infall peculiar velocity. At higher redshift, z=6.6, where the transmissivity is instead largely set by the local HI abundance and LAE transmissivity consequently increases with halo mass from 0.15 to 0.3. Although outflows are present, they are efficiently pressure-confined by infall in a small region around the LAE; hence they only marginally affect transmissivity. Finally, we cast LOS originating from background QSOs passing through foreground LAEs at different impact parameters, and compute the quasar transmissivity (T_QSO). At smaller impact parameters, d < 1 cMpc, a positive correlation between T_QSO and halo mass is found at z = 5.7, which tends to become less pronounced (i.e. flatter) at larger distances. By cross-correlating T_LAE and T_QSO, we can obtain a HI density estimate unaffected by dust. At z= 5.7, the cross-correlation is relatively weak,whereas at z = 6.6 we find a clear positive correlation. We conclude by briefly discussing the perspectives for the application of the method to existing and forthcoming data.

It is predicted that sources emitting UV radiation in the Lyman band during the epoch of reionization (EoR) showed a series of discontinuities in their Ly-alpha flux radial profile as a consequence of the thickness of the Lyman line series in the primeval intergalactic medium. Through unsaturated Wouthuysen-Field coupling, these spherical discontinuities are also present in the 21 cm emission of the neutral IGM. In this article, we study the effects these discontinuities have on the differential brightness temperature of the 21 cm signal of neutral hydrogen in a realistic setting including all other sources of fluctuations. We focus on the early phases of the EoR, and we address the question of the detectability by the planned Square Kilometre Array. Such a detection would be of great interest, because these structures could provide an unambiguous diagnostic for the cosmological origin of the signal remaining after the foreground cleaning procedure. Also, they could be used as a new type of standard rulers. We determine the differential brightness temperature of the 21 cm signal in the presence of inhomogeneous Wouthuysen-Field effect using simulations which include (hydro)dynamics and both ionizing and Lyman lines 3D radiative transfer with the code LICORICE. We find that the Lyman horizons are clearly visible on the maps and radial profiles around the first sources of our simulations, but for a limited time interval, typically \Delta z \approx 2 at z \sim 13. Stacking the profiles of the different sources of the simulation at a given redshift results in extending this interval to \Delta z \approx 4. When we take into account the implementation and design planned for the SKA (collecting area, sensitivity, resolution), we find that detection will be challenging. It may be possible with a 10 km diameter for the core, but will be difficult with the currently favored design of a 5 km core.

We present simulations of the 21-cm signal during the epoch of reionization. We focus on properly modeling the absorption regime in the presence of inhomogeneous Wouthuysen-Field effect and X-ray heating. We ran radiative transfer simulations for three bands in the source spectrum (Lyman, UV, and X-ray) to fully account for these processes. We find that the brightness temperature fluctuation of the 21 cm signal has an amplitude greater than 100 mK during the early reionization, up to 10 times greater than the typical amplitude of a few 10 mK obtained during the later emission phase. More importantly, we find that even a rather high contribution from QSO-like sources only damps the absorption regime without erasing it. Heating the IGM with X-ray takes time. Our results show that observations of the early reionization will probably benefit from a higher signal-to-noise value than during later stages. After analyzing the statistical properties of the signal (power spectrum and PDF) we find three diagnostics to constrain the level of X-ray, hence the nature of the first sources.

The development of radiation hydrodynamical methods that are able to follow gas dynamics and radiative transfer self-consistently is key to the solution of many problems in numerical astrophysics. Such fluid flows are highly complex, rarely allowing even for approximate analytical solutions against which numerical codes can be tested. An alternative validation procedure is to compare different methods against each other on common problems, in order to assess the robustness of the results and establish a range of validity for the methods. Previously, we presented such a comparison for a set of pure radiative transfer tests (i.e. for fixed, non-evolving density fields). This is the second paper of the Cosmological Radiative Transfer (RT) Comparison Project, in which we compare 9 independent RT codes directly coupled to gasdynamics on 3 relatively simple astrophysical hydrodynamics problems: (5) the expansion of an H II region in a uniform medium; (6) an ionization front (I-front) in a 1/r^2 density profile with a flat core, and (7), the photoevaporation of a uniform dense clump. Results show a broad agreement between the different methods and no big failures, indicating that the participating codes have reached a certain level of maturity and reliability. However, many details still do differ, and virtually every code has showed some shortcomings and has disagreed, in one respect or another, with the majority of the results. This underscores the fact that no method is universal and all require careful testing of the particular features which are most relevant to the specific problem at hand.

The 21 cm emission of neutral hydrogen is the most promising probe of the epoch of reionization(EoR). In the next few years, the SKA pathfinders will provide statistical measurements of this signal. Numerical simulations predicting these observations are necessary to optimize the design of the instruments. The main difficulty is the computation of the spin temperature of neutral hydrogen which depends on the gas kinetic temperature and on the level of the local Lyman-alpha flux. A T_s >> T_cmb assumption is usual. However, this assumption does not apply early in the reionization history, or even later in the history as long as the sources of X-rays are too weak to heat the intergalactic medium significantly. This work presents the first EoR numerical simulations including, beside dynamics and ionizing continuum radiative transfer, a self-consistent treatment of the Ly-alpha radiative transfer. This allows us to compute the spin temperature more accurately. We use two different box sizes, 20 Mpc/h and 100 Mpc/h, and a star source model. Using the redshift dependence of average quantities, maps, and power spectra, we quantify the effect of using different assumptions to compute the spin temperature and the influence of the box size. The first effect comes from allowing for a signal in absorption. The magnitude of this effect depends on the amount of heating by hydrodynamic shocks and X-rays in the intergalactic medium(IGM). The second effects comes from using the real, local, Lyman-alpha flux. This effect is important for an average ionization fraction of less than 10%: it changes the overall amplitude of the 21 cm signal, and adds its own fluctuations to the power spectrum.

During the epoch of reionization, Ly-alpha photons emitted by the first stars can couple the neutral hydrogen spin temperature to the kinetic gas temperature, providing the opportunity to observe the gas in emission or absorption in the 21-cm line. Given the bright foregrounds, it is of prime importance to determine precisely the fluctuations signature of the signal, to be able to extract it by its correlation power. LICORICE is a Monte-Carlo radiative transfer code, coupled to the dynamics via an adaptative Tree-SPH code. We present here the Ly-alpha part of the implementation, and validate it through three classical tests. Contrary to previous works, we do not assume that P_alpha, the number of scatterings of Ly-alpha photons per atom per second, is proportional to the Ly-alpha background flux, but take into account the scatterings in the Ly-alpha line wings. The latter have the effect to steepen the radial profile of P_alpha around each source, and re-inforce the contrast of the fluctuations. In the particular geometry of cosmic filaments of baryonic matter, Ly-alpha photons are scattered out of the filament, and the large scale structure of P_alpha is significantly anisotropic. This could have strong implications for the possible detection of the 21-cm signal.

We present 3.65''x3.34'' angular-resolution IRAM PdBI observations of the CS J=2-1 line toward the Horsehead Photodissociation Region (PDR), complemented with IRAM-30m single-dish observations of several rotational lines of CS, C34S and HCS+. We analyse the CS and HCS+ photochemistry, excitation and radiative transfer to obtain their abundances and the physical conditions prevailing in the cloud edge. Since the CS abundance scales to that of sulfur, we determine the gas phase sulfur abundance in the PDR, an interesting intermediate medium between translucent clouds (where sulfur remains in the gas phase) and dark clouds (where large depletions have been invoked). A nonlocal non-LTE radiative transfer code including dust and cosmic background illumination adapted to the Horsehead geometry has been developed to carefuly analyse the CS, C34S, HCS+ and C18O rotational line emission. We use this model to consistently link the line observations with photochemical models to determine the CS/HCS+/S/S+ structure of the PDR. The combination of CS, C34S and HCS+ observations together with the inclusion of the most recent CS collisional and chemical rates in our models implies that sulfur depletion invoked to account for CS and HCS+ abundances is much smaller than in previous studies.

The direct detection of neutralino dark matter is analysed in general supergravity scenarios, where non-universal soft scalar and gaugino masses can be present. In particular, the theoretical predictions for the neutralino-nucleon cross section are studied and compared with the sensitivity of dark matter detectors. We take into account the most recent astrophysical and experimental constraints on the parameter space, including the current limit on B(Bs-> mu+ mu-). The latter puts severe limitations on the dark matter scattering cross section, ruling out most of the regions that would be within the reach of present experiments. We show how this constraint can be softened with the help of appropriate choices of non-universal parameters which increase the Higgsino composition of the lightest neutralino and minimise the chargino contribution to the b->s transition.