It is commonly assumed that the energy density of the Universe was dominated by radiation between reheating after inflation and the onset of matter domination 54,000 years later. While the abundance of light elements indicates that the Universe was radiation dominated during Big Bang Nucleosynthesis (BBN), there is scant evidence that the Universe was radiation dominated prior to BBN. It is therefore possible that the cosmological history was more complicated, with deviations from the standard radiation domination during the earliest epochs. Indeed, several interesting proposals regarding various topics such as the generation of dark matter, matter-antimatter asymmetry, gravitational waves, primordial black holes, or microhalos during a nonstandard expansion phase have been recently made. In this paper, we review various possible causes and consequences of deviations from radiation domination in the early Universe - taking place either before or after BBN - and the constraints on them, as they have been discussed in the literature during the recent years.

We present a complementarity study of gravitational waves and double Higgs production in the $4b$ channel, exploring the gauge singlet scalar extension of the SM. This new physics extension serves as a simplified benchmark model that realizes a strongly first-order electroweak phase transition necessary to generate the observed baryon asymmetry in the universe. In calculating the signal-to-noise ratio of the gravitational waves, we incorporate the effect of the recently discovered significant suppression of the gravitational wave signals from sound waves for strong phase transitions, make sure that supercooled phase transitions do complete and adopt a bubble wall velocity that is consistent with a successful electroweak baryogenesis by solving the velocity profiles of the plasma. The high-luminosity LHC sensitivity to the singlet scalar extension of the SM is estimated using a shape-based analysis of the invariant $4b$ mass distribution. We find that while the region of parameter space giving detectable gravitational waves is shrunk due to the new gravitational wave simulations, the qualitative complementary role of gravitational waves and collider searches remain unchanged.

Supernovae and cooling neutron stars have long been used to constrain the properties of axions, such as their mass and interactions with nucleons and other Standard Model particles. We investigate the prospects of using neutron star mergers as a similar location where axions can be probed in the future. We examine the impact axions would have on mergers, considering both the possibility that they free-stream through the dense nuclear matter and the case where they are trapped. We calculate the mean free path of axions in merger conditions, and find that they would free-stream through the merger in all thermodynamic conditions. In contrast to previous calculations, we integrate over the entire phase space while using a relativistic treatment of the nucleons, assuming the matrix element is momentum-independent. In particular, we use a relativistic mean field theory to describe the nucleons, taking into account the precipitous decrease in the effective mass of the nucleons as density increases above nuclear saturation density. We find that within current constraints on the axion-neutron coupling, axions could cool nuclear matter on timescales relevant to neutron star mergers. Our results may be regarded as first steps aimed at understanding how axions affect merger simulations and potentially interface with observations.

Axion-like-particles (ALPs) emitted from the core of magnetars can convert to photons in the strong magnetic field of the magnetosphere. We study such emissions in the soft gamma-ray range from 300 keV to 1 MeV, where the ALP spectrum peaks and astrophysical backgrounds from resonant Compton upscattering are expected to be suppressed. Using published quiescent soft-gamma flux upper limits in 6 Magnetars obtained with $\textit{CGRO}$ COMPTEL, $\textit{INTEGRAL}$ SPI/IBIS/ISGRI and the $\textit{Fermi}$ Gamma Ray Burst Monitor (GBM), we put limits on the ALP-photon coupling obtained from conversions, assuming that ALP emission from the core is just sub-dominant to bounds from neutrino cooling. For core temperatures of $10^9$ K, the constraints on the ALP-photon coupling coming from magnetars 1E 2259+586 and J170849.0-400910 are better than the current limits obtained from the CAST experiment. We provide a detailed study of the dependence of our results on the magnetar core temperature. Our results motivate a program of studying quiescent soft gamma ray emission from magnetars in the 300 keV - 1 MeV band with the $\textit{Fermi}$-GBM.

Rather general considerations of the string theory landscape imply a mild statistical draw towards large soft SUSY breaking terms tempered by the requirement of proper electroweak symmetry breaking where SUSY contributions to the weak scale are not too far from m(weak)~ 100 GeV. Such a picture leads to the prediction that m_h~ 125 GeV while most sparticles are beyond current LHC reach. Here we explore the possibility that the magnitude of the Peccei-Quinn (PQ) scale f_a is also set by string landscape considerations within the framework of a compelling SUSY axion model. First, we examine the case where the PQ symmetry arises as an accidental approximate global symmetry from a more fundamental gravity-safe Z(24)^R symmetry and where the SUSY mu parameter arises from a Kim-Nilles operator. The pull towards large soft terms then also pulls the PQ scale as large as possible. Unless this is tempered by rather severe (unknown) cosmological or anthropic bounds on the density of dark matter, then we would expect a far greater abundance of dark matter than is observed. This conclusion cannot be negated by adopting a tiny axion misalignment angle theta_i because WIMPs are also overproduced at large f_a. Hence, we conclude that setting the PQ scale via anthropics is highly unlikely. Instead, requiring soft SUSY breaking terms of order the gravity-mediation scale m_{3/2}~ 10-100 TeV places the mixed axion-neutralino dark matter abundance into the intermediate scale sweet zone where f_a~ 10^{11}-10^{12} GeV. We compare our analysis to the more general case of a generic SUSY DFSZ axion model with uniform selection on theta_i but leading to the measured dark matter abundance: this approach leads to a preference for f_a~ 10^{12} GeV.

This white paper highlights compact object and fundamental physics science opportunities afforded by high-throughput broadband (0.1-60 keV) X-ray polarization observations. X-ray polarimetry gives new observables with geometric information about stellar remnants which are many orders of magnitude too small for direct imaging. The X-ray polarimetric data also reveal details about the emission mechanisms and the structure of the magnetic fields in and around the most extreme objects in the Universe. Whereas the Imaging X-ray Polarimetry Explorer (IXPE) to be launched in 2021 will obtain first results for bright objects, a follow-up mission could be one order of magnitude more sensitive and would be able to use a broader bandpass to perform physics type experiments for representative samples of sources.

We propose to use gravitational waves from extreme mass ratio inspirals (EMRI), composed of a boson star and a supermassive black hole in the center of galaxies, as a new method to search for boson stars. Gravitational waves from EMRI have the advantage of being long-lasting within the frequency band of future space-based interferometer gravitational wave detectors and can accumulate large signal-to-noise ratio (SNR) for very sub-solar mass boson stars. Compared to gravitational waves from boson star binaries, which fall within the LIGO band, we find that much larger ranges of the mass and compactness of boson stars, as well as the underlying particle physics parameter space, can be probed by EMRI. We take tidal disruption of the boson stars into account and distinguish those which dissolve before the inner-most-stable-circular-orbit (ISCO) and those which dissolve after it. Due to the large number of cycles recorded, EMRIs can lead to a very precise mass determination of the boson star and distinguish it from standard astrophysical compact objects in the event of a signal. Tidal effects in inspiralling binary systems, as well as possible correlated electromagnetic signals, can also serve as potential discriminants.

We present a dedicated complementarity study of gravitational wave and collider measurements of the simplest extension of the Higgs sector: the singlet scalar augmented Standard Model. We study the following issues: (i) the electroweak phase transition patterns admitted by the model, and the proportion of parameter space for each pattern; (ii) the regions of parameter space that give detectable gravitational waves at future space-based detectors; and (iii) the current and future collider measurements of di-Higgs production, as well as searches for a heavy weak diboson resonance, and how these searches interplay with regions of parameter space that exhibit strong gravitational wave signals. We carefully investigate the behavior of the normalized energy released during the phase transition as a function of the model parameters, address subtle issues pertaining to the bubble wall velocity, and provide a description of different fluid velocity profiles. On the collider side, we identify the subset of points that are most promising in terms of di-Higgs and weak diboson production studies while also giving detectable signals at LISA, setting the stage for future benchmark points that can be used by both communities.

Axion-like-particles (ALPs) produced in the core of a magnetar can convert to photons in the magnetosphere, giving rise to novel features in the X-ray spectrum. Since ALPs only mix with the parallel mode of the photon, the polarization of the soft and hard X-ray spectra is predicted to have an O-mode component, in addition to the mainly X-mode component given by most astrophysical models. The relative strength of the O-mode component depends on the intensity of ALPs produced in the core and the probability of conversion. We quantify our results by considering X-ray emission produced both by astrophysical processes and by ALP-photon conversion, in an uncorrelated fashion, and in different relative proportions, which we parametrize by the angle $\chi_0$. We then define a normalized astrophysics-subtracted Stokes parameter $R$ which only acquires non-zero values in the presence of ALP-photon conversion. We find, remarkably, that the parameter $R$ factorizes into a product of the ALP-to-photon conversion probability and $\cos(2\chi_0)$ and display $R$, as well as the usual Stokes parameter $Q$, as a function of the photon energy and relative fractions of ALP and photon intensities. For benchmark points currently allowed by the CAST experiment, the O-mode prediction can be tested in future X-ray polarimeters and used either to constrain ALPs or find evidence for them.

Axion-like particles (ALPs) produced in the core of a magnetar will convert to photons in the magnetosphere, leading to possible signatures in the hard X-ray band. We perform a detailed calculation of the ALP-to-photon conversion probability in the magnetosphere, recasting the coupled differential equations that describe ALP-photon propagation into a form that is efficient for large scale numerical scans. We show the dependence of the conversion probability on the ALP energy, mass, ALP-photon coupling, magnetar radius, surface magnetic field, and the angle between the magnetic field and direction of propagation. Along the way, we develop an analytic formalism to perform similar calculations in more general $n$-state oscillation systems. Assuming ALP emission rates from the core that are just subdominant to neutrino emission, we calculate the resulting constraints on the ALP mass versus ALP-photon coupling space, taking SGR 1806-20 as an example. In particular, we take benchmark values for the magnetar radius and core temperature, and constrain the ALP parameter space by the requirement that the luminosity from ALP-to-photon conversion should not exceed the total observed luminosity from the magnetar. The resulting constraints are competitive with constraints from helioscope experiments in the relevant part of ALP parameter space.

The string theory landscape of vacua solutions provides physicists with some understanding as to the magnitude of the cosmological constant. Similar reasoning can be applied to the magnitude of the soft SUSY breaking terms in supersymmetric models of particle physics: there appears to be a statistical draw towards large soft terms which is tempered by the anthropic requirement of the weak scale lying not too far from ~100 GeV. For a mild statistical draw of m_{soft}^n with n=1 (as expected from SUSY breaking due to a single F term) then the light Higgs mass is preferred at ~125 GeV while sparticles are all pulled beyond LHC bounds. We confront a variety of LHC and WIMP dark matter search limits with the statistical expectations from a fertile patch of string theory landscape. The end result is that LHC and WIMP dark matter detectors see exactly that which is expected from the string theory landscape: a Standard Model-like Higgs boson of mass 125 GeV but as yet no sign of sparticles or WIMP dark matter. SUSY from the n=1 landscape is most likely to emerge at LHC in the soft opposite-sign dilepton plus jet plus MET channel. Multi-ton noble liquid WIMP detectors should be able to completely explore the n=1 landscape parameter space.

The density spike of dark matter (DM) in the subparsec region near the supermassive black hole at the Galactic Center can provide potentially observable gamma-ray signals coming from DM annihilations. Taking Fermi-$LAT$ data for the gamma-ray flux from the point source 3FGL J1745.6-2859c (Sgr A*), we calculate the resulting constraints on generic models of DM, allowing for the possibility of a non-negligible velocity-dependent component of the annihilation cross section. We consider a variety of selections for the astrophysical parameters that describe the spike profile and find that the gamma-ray flux is strongly dependent on these selections, particularly the modelling of spike depletion effects due to gravitational interactions with baryons, which affect the spike radius and steepness profile. We calculate constraints on the DM parameter space for both attenuated and idealized spikes, considering different choices for the steepness profiles in each case. We find that for the most conservative selection of parameters, corresponding to a depleted spike with an NFW cusp profile, the gamma-ray flux for a 100 GeV thermal relic is lower than current observational constraints by several orders of magnitude. For more optimistic choices of parameters corresponding to spikes that have not been attenuated, bounds on DM masses can be obtained for a variety of choices of steepness profiles. We then specialize to a class of simplified models of fermionic DM that annihilate dominantly through the $t-$channel exchange of two scalar mediators with arbitrary mixing angle $\alpha$, and calculate the indirect detection constraints coming from the DM spike. Along the way, we discuss constraints on the astrophysical parameters describing the DM spike, taking the Galactic Center excess as a signal of DM annihilation.

Dwarf spheroidal galaxies (dSphs) are promising targets for the indirect detection of dark matter through gamma-ray emission due to their proximity, lack of astrophysical backgrounds and high dark matter density. They are often used to place restrictive bounds on the dark matter annihilation cross section. In this paper, we analyze six years of {\it Fermi}-LAT gamma-ray data from 19 dSphs that are satellites of the Milky Way, and derive from a stacked analysis of 15 dSphs, robust 95\% confidence level lower limits on the dark matter lifetime for several decay channels and dark matter masses between $\sim 1$GeV and $10$TeV. Our findings are based on a bin-by-bin maximum likelihood analysis treating the J-factor as a nuisance parameter using PASS 8 event-class. Our constraints from this ensemble are among the most stringent and solid in the literature, and competitive with existing ones coming from the extragalactic gamma-ray background, galaxy clusters, AMS-02 cosmic ray data, Super-K and ICECUBE neutrino data, while rather insensitive to systematic uncertainties. In particular, among gamma-ray searches, we improve existing limits for dark matter decaying into $\bar{b}b$, ($\mu^+\mu^-$) for DM masses below $\sim 30\, (200)$~GeV, demonstrating that dSphs are compelling targets for constraining dark matter decay lifetimes.

In standard cosmology, the growth of structure becomes significant following matter-radiation equality. In non-thermal histories, where an effectively matter-dominated phase occurs due to scalar oscillations prior to Big Bang Nucleosynthesis, a new scale at smaller wavelengths appears in the matter power spectrum. Density perturbations that enter the horizon during the matter-dominated phase grow linearly with the scale factor prior to the onset of radiation domination, which leads to enhanced inhomogeneity on small scales if dark matter thermally and kinetically decouples during the matter-dominated phase. The microhalos that form from these enhanced perturbations significantly boost the self-annihilation rate for dark matter. This has important implications for indirect detection experiments: the larger annihilation rate will result in observable signals from dark matter candidates that are usually deemed untestable. As a proof of principle, we consider Binos in heavy supersymmetry with an intermediate extended Higgs sector and all other superpartners decoupled. We find that these isolated Binos, which lie under the neutrino floor, can account for the dark matter relic density while also leading to observable predictions for Fermi-LAT. Current limits on the annihilation cross section from Fermi-LAT's observations of dwarf spheroidal galaxies may already constrain Bino dark matter up to masses $\mathcal{O}(300)$ GeV, depending on the internal structure of the microhalos. More extensive constraints are possible with improved gamma-ray bounds and boost calculations from $N$-body simulations.

In this note we propose a model independent framework for inflationary (p)reheating. Our approach is analogous to the Effective Field Theory of Inflation, however here the inflaton oscillations provide an additional source of (discrete) symmetry breaking. Using the Goldstone field that non-linearly realizes time diffeormorphism invariance we construct a model independent action for both the inflaton and reheating sectors. Utilizing the hierarchy of scales present during the reheating process we are able to recover known results in the literature in a simpler fashion, including the presence of oscillations in the primordial power spectrum. We also construct a class of models where the shift symmetry of the inflaton is preserved during reheating, which helps alleviate past criticisms of (p)reheating in models of Natural Inflation. Extensions of our framework suggest the possibility of analytically investigating non-linear effects (such as rescattering and back-reaction) during thermalization without resorting to lattice methods. By construction, the EFT relates the strength of many of these interactions to other operators in the theory, including those responsible for the efficiency of (p)reheating. We conclude with a discussion of the limitations and challenges for our approach.

We critically review the role of cosmological moduli in determining the post-inflationary history of the universe. Moduli are ubiquitous in string and M-theory constructions of beyond the Standard Model physics, where they parametrize the geometry of the compactification manifold. For those with masses determined by supersymmetry breaking this leads to their eventual decay slightly before Big Bang Nucleosynthesis (without spoiling its predictions). This results in a matter dominated phase shortly after inflation ends, which can influence baryon and dark matter genesis, as well as observations of the Cosmic Microwave Background and the growth of large-scale structure. Given progress within fundamental theory, and guidance from dark matter and collider experiments, non-thermal histories have emerged as a robust and theoretically well-motivated alternative to a strictly thermal one. We review this approach to the early universe and discuss both the theoretical challenges and the observational implications.

We exploit the complementarity among supersymmetry, inflation, axions, Big Bang Nucleosynthesis (BBN) and Cosmic Microwave Background Radiation (CMB) to constrain supersymmetric axion models in the light of the recent Planck and BICEP results. In particular, we derive BBN bounds coming from altering the light element abundances by taking into account hadronic and electromagnetic energy injection, and CMB constraints from black-body spectrum distortion. Lastly, we outline the viable versus excluded region of these supersymetric models that might account for the mild dark radiation observed.

A detection of primordial B-modes has been heralded not only as a smoking gun for the existence of inflation, but also as a way to establish the scale at which inflation took place. In this paper we critically reinvestigate the connection between a detection of primordial gravity waves and the scale of inflation. We consider whether the presence of additional fields and non-adiabaticity during inflation may have provided an additional source of primordial B-modes competitive with those of the quasi-de Sitter vacuum. In particular, we examine whether the additional sources could provide the dominant signal, which could lead to a misinterpretation of the scale of inflation. In light of constraints on the level of non-Gaussianity coming from Planck we find that only hidden sectors with strictly gravitationally strength couplings provide a feasible mechanism. The required model building is somewhat elaborate, and so we discuss possible UV completions in the context of Type IIB orientifold compactifications with RR axions. We find that an embedding is possible and that dangerous sinusoidal corrections can be suppressed through the compactification geometry. Our main result is that even when additional sources of primordial gravity waves are competitive with the inflaton, a positive B-mode detection would still be a relatively good indicator of the scale of inflation. This conclusion will be strengthened by future constraints on both non-Gaussianity and CMB polarization.

We study the viability of pseudo Nambu-Goldstone bosons (Majorons) arising in see-saw models as dark matter candidates. Interestingly the stability of the Majoron as dark matter is related to the scale that sets the see-saw and leptogenesis mechanisms, while its annihilation and scattering cross section off nuclei can be set through the Higgs portal. For $\mathcal{O}(GeV) - \mathcal{O}(TeV)$ Majorons, we compute observables such as the abundance, scattering cross section, Higgs invisible decay width, and emission lines and compare with current data in order to outline the excluded versus still viable parameter space regions. We conclude that the simplest Majoron dark matter models coupling through the Higgs portal, except at the Higgs resonance, are excluded by current direct detection data for Majorons lighter than $225$~GeV and future runnings are expected to rule out decisively the 1GeV-1TeV window. Lastly, we point out that light keV-scale Majorons whose relic density is set by thermal freeze-in from sterile neutrinos can account for the keV line observed by XMM-Newton observatory in the spectrum of 73 galaxy clusters, within a see-saw model with a triplet Higgs.

Following our predictions on detectability, very weak lines of the SH molecules have been reported as identified in the photospheric spectrum of the Sun. This could be the first detection of a mercapto radical in the solar spectrum, thus placing confidence in the solar abundance of Sulfur. The observations additionally tested the utilized theoretical band oscillator strength. Sunspots being cooler than the photosphere are the hosts to a large number of molecular species. However, owing to photospheric radiation scattered into the observed umbra, the discoverers could not detect the lines of SH in the sunspot spectrum where they are expected to show larger than the photospheric equivalent widths (EWs in short). Detection of weak lines in photospheric spectrum coupled with a no-detection of the relatively strong same lines in the sunspot spectrum might cause doubts on the detection itself. The above problem is investigated here in detail with a choice on photospheric and sunspot models and micro-turbulence values. The new results indicate that the predicted sunspot lines are about half intense than reported before and the lines remain of detectable strengths. Further, a need for a laboratory determination of oscillator strengths is felt.

Reheating in string compactifications is generically driven by the decay of the lightest modulus which produces Standard Model particles, dark matter and light hidden sector degrees of freedom that behave as dark radiation. This common origin allows us to find an interesting correlation between dark matter and dark radiation. By combining present upper bounds on the effective number of neutrino species N_eff with lower bounds on the reheating temperature as a function of the dark matter mass m_DM from Fermi data, we obtain strong constraints on the (N_eff,m_DM)-plane. Most of the allowed region in this plane corresponds to non-thermal scenarios with Higgsino-like dark matter. Thermal dark matter can be allowed only if N_eff tends to its Standard Model value. We show that the above situation is realised in models with perturbative moduli stabilisation where the production of dark radiation is unavoidable since bulk closed string axions remain light and do not get eaten up by anomalous U(1)s.

Non-thermal cosmological histories are capable of greatly increasing the available parameter space of different particle physics dark matter (DM) models and are well-motivated by the ubiquity of late-decaying gravitationally coupled scalars in UV theories like string theory. A non-thermal DM model is presented in the context of LARGE Volume Scenarios in type IIB string theory. The model is capable of addressing both the moduli-induced gravitino problem as well as the problem of overproduction of axionic dark radiation and/or DM. We show that the right abundance of neutralino DM can be obtained in both thermal under and overproduction cases for DM masses between O(GeV) to O(TeV). In the latter case the contribution of the QCD axion to the relic density is totally negligible, while in the former case it can be comparable to that of the neutralino thus resulting in a multi-component DM scenario.

If the lightest supersymmetric particle (LSP) is Higgsino-like, the thermal relic density is lower than the observed dark matter content for a LSP mass in the sub-TeV region. We outline constraints arising from the Fermi Gamma-Ray Telescope data and LSP production from gravitino decay that must be satisfied by a successful non-thermal Higgsino scenario. We show that in a generic class of models where anomaly and modulus mediated contributions to supersymmetry breaking are of comparable size, Higgsino arises as the only viable sub-TeV dark matter candidate if gravitinos are heavy enough to decay before the onset of big bang nucleosynthesis (BBN). The correct relic density can be obtained via modulus decay in these models. As an explicit example, we consider a modulus sector in effective field theory ($D=4, N=1$ supergravitiy arising from type IIB KKLT compactification). Within this class of mirage mediation models, heaviness of the gravitino forces a sub-TeV Higgsino LSP and gives a Higgs mass around 125 GeV. In this example, the constraints from direct detection experiments are also satisfied.

Studies of the initial conditions for inflation have conflicting predictions from exponential suppression to inevitability. At the level of phase space, this conflict arises from the competing intuitions of CPT invariance and thermodynamics. After reviewing this conflict, we enlarge the ensemble beyond phase space to include scalar potential data. We show how this leads to an important contribution from inflection point inflation, enhancing the likelihood of inflation to an inverse cubic power law. In the process, we emphasize the attractor dynamics of the gravity-scalar system and the existence of universality classes from inflection point inflation. Finally, we comment on the predictivity of inflation in light of these results.

We study inflection point inflation using Singularity Theory, which relates degenerate critical points of functions to their local behavior. This approach illuminates universal features of small-field models and gives analytic control over parametrized families of scalar potentials near inflationary solutions. The behavior of the scalar potential is tied to the number of physical input parameters, which determines a set of universality classes. Within these classes, we obtain universal scaling relations for density perturbations and the scale of inflation. In specific models, we show that the scale of supersymmetry breaking also possesses scaling behavior. We illustrate this general structure with a specific example: the Racetrack Inflation model in type IIB string theory, with the inflaton being the real part of the Kahler modulus, and the input parameters being flux dependent quantities that appear in the 4D, N=1 superpotential.

Recent anomalies in cosmic rays could be due to dark matter annihilation in our galaxy. In order to get the required large cross-section to explain the data while still obtaining the right relic density, we rely on a non standard thermal history between dark matter freeze-out and Big-Bang Nucleosynthesis (BBN). We show that through a reheating phase from the decay of a heavy moduli or even the gravitino, we can produce the right relic density of dark matter if its self-annihilation cross-section is large enough. In addition to fitting the recent data, this scenario solves the cosmological moduli and gravitino problems. We illustrate this mechanism with a specific example in the context of U(1)_{B-L} extended MSSM where supersymmetry is broken via mirage mediation. These string motivated models naturally contain heavy moduli decaying to the gravitino, whose subsequent decay to the LSP can reheat the universe at a low temperature. The right-handed sneutrino and the B-L gaugino can both be viable dark matter candidates with large cross-section. They are leptophilic because of B-L charges. We also show that it is possible to distinguish the non-thermal from the thermal scenario (using Sommerfeld enhancement) in direct detection experiments for certain regions of parameter space.

Identification of Fraunhofer lines with the known atomic and molecular absorbers and predictions leading to such an effort has been a challenging area of study crowned with occasional success. Such studies have also lead, amongst other things to (i) a determination of abundances of elements and that of their isotopes (ii) valuable information on model atmospheres and (iii) use of Sun as a laboratory source. We summarize and review here the work done in the last four decades in the area of solar spectroscopy at Aryabhatta Research Institute of observational sciencES (ARIES in short) with a view to pick up new and interesting areas for future investigations in the light of the tremendous progress made elsewhere in observations of the sun and in the laboratory studies.