NASA's Interface Region Imaging Spectrograph (IRIS) provides high resolution observations of the solar atmosphere through UV spectroscopy and imaging. Since the launch of IRIS in June 2013, we have conducted systematic observation campaigns in coordination with the Swedish 1-m Solar Telescope (SST) on La Palma. The SST provides complementary high-resolution observations of the photosphere and chromosphere. The SST observations include spectro-polarimetric imaging in photospheric Fe I lines and spectrally-resolved imaging in the chromospheric Ca II 8542 A, H-alpha, and Ca II K lines. We present a database of co-aligned IRIS and SST datasets that is open for analysis to the scientific community. The database covers a variety of targets including active regions, sunspots, plage, quiet Sun, and coronal holes.

We present an efficient numerical code and conduct, for the first time, a null and model-independent CMB test of statistical isotropy using Multipole Vectors (MVs) at all scales. Because MVs are insensitive to the angular power spectrum $C_\ell$, our results are independent from the assumed cosmological model. We avoid a posteriori choices and use pre-defined ranges of scales $\ell\in[2,30]$, $\ell\in[2,600]$ and $\ell\in[2,1500]$ in our analyses. We find that all four masked Planck maps, from both 2015 and 2018 releases, are in agreement with statistical isotropy for $\ell\in[2,30]$, $\ell\in[2,600]$. For $\ell\in[2,1500]$ we detect anisotropies but this is indicative of simply the anisotropy in the noise: there is no anisotropy for $\ell < 1300$ and an increasing level of anisotropy at higher multipoles. Our findings of no large-scale anisotropies seem to be a consequence of avoiding \emph{a posteriori} statistics. We also find that the degree of anisotropy in the full sky (i.e. unmasked) maps vary enormously (between less than 5 and over 1000 standard deviations) among the different mapmaking procedures and data releases.

Within cosmological perturbation theory, the cosmic microwave background anisotropies are usually computed from a Boltzmann hierarchy coupled to the perturbed Einstein equations. In this setup, one set of multipoles describes the temperature anisotropies, while two other sets, of electric and magnetic types, describe the polarization anisotropies. In order to reduce the number of multipoles types needed for polarization, and thus to speed up the numerical resolution, an optimal hierarchy has been proposed in the literature for Einstein-Boltzmann codes. However, it has been recently shown that the separability between directional and orbital eigenfunctions employed in the optimal hierarchy is not correct in the presence of spatial curvature. We investigate how the assumption of separability affects the optimal hierarchy, and show that it introduces relative errors of order $\Omega_K$ with respect to the full hierarchy. Despite of that, we show that the optimal hierarchy still gives extremely good results for temperature and polarization angular spectra, with relative errors that are much smaller than cosmic variance even for curvatures as large as $|\Omega_K|=0.1$. Still, we find that the polarization angular spectra from tensor perturbations are significantly altered when using the optimal hierarchy, leading to errors that are typically of order $50 |\Omega_K| \%$ on that component.

The SunPy Project developed a 13-question survey to understand the software and hardware usage of the solar physics community. 364 members of the solar physics community, across 35 countries, responded to our survey. We found that 99$\pm$0.5% of respondents use software in their research and 66% use the Python scientific software stack. Students are twice as likely as faculty, staff scientists, and researchers to use Python rather than Interactive Data Language (IDL). In this respect, the astrophysics and solar physics communities differ widely: 78% of solar physics faculty, staff scientists, and researchers in our sample uses IDL, compared with 44% of astrophysics faculty and scientists sampled by Momcheva and Tollerud (2015). 63$\pm$4% of respondents have not taken any computer-science courses at an undergraduate or graduate level. We also found that most respondents utilize consumer hardware to run software for solar-physics research. Although 82% of respondents work with data from space-based or ground-based missions, some of which (e.g. the Solar Dynamics Observatory and Daniel K. Inouye Solar Telescope) produce terabytes of data a day, 14% use a regional or national cluster, 5% use a commercial cloud provider, and 29% use exclusively a laptop or desktop. Finally, we found that 73$\pm$4% of respondents cite scientific software in their research, although only 42$\pm$3% do so routinely.

We investigate a corner of the Bianchi models that has not received much attention: "extended FLRW models" (eFLRW) defined as a cosmological model with underlying anisotropic Bianchi geometry that nevertheless expands isotropically and can be mapped onto a reference FLRW model with the same expansion history. In order to investigate the stability and naturalness of such models in a dynamical systems context, we consider spatially homogeneous models that contain a massless scalar field $\varphi$ and a non-tilted perfect fluid obeying an equation of state $p=w\rho$. Remarkably, we find that matter anisotropies and geometrical anisotropies tend to cancel out dynamically. Hence, the expansion is asymptotically isotropic under rather general conditions. Although extended FLRW models require a special matter sector with anisotropies that are 'fine-tuned" relative to geometrical anisotropies, our analysis shows that such solutions are dynamically preferred attractors in general relativity. Specifically, we prove that all locally rotationally symmetric Bianchi type III universes with space-like $\nabla_\mu\varphi$ are asymptotically shear-free, for all $w\in[-1,1]$. Moreover, all shear-free equilibrium sets with anisotropic spatial curvature are proved to be stable with respect to all homogeneous perturbations for $w\geq -1/3$.

A powerful result in theoretical cosmology states that a subset of anisotropic Bianchi models can be seen as the homogeneous limit of (standard) linear cosmological perturbations. Such models are precisely those leading to Friedmann spacetimes in the limit of zero anisotropy. Building on previous works, we give a comprehensive exposition of this result, and perform the detailed identification between anisotropic degrees of freedom and their corresponding scalar, vector, and tensor perturbations of standard perturbation theory. In particular, we find that anisotropic models very close to open (i.e., negatively curved) Friedmann spaces correspond to some type of super-curvature perturbations. As a consequence, provided anisotropy is mild, its effects on all types of cosmological observables can always be computed as simple extensions of the standard techniques used in relativistic perturbation theory around Friedmann models. This fact opens the possibility to consistently constrain, for all cosmological observables, the presence of large scale anisotropies on the top of the stochastic fluctuations.

Context. The solar chromosphere and the lower transition region is believed to play a crucial role in the heating of the solar corona. Models that describe the chromosphere (and the lower transition region), accounting for its highly dynamic and structured character are, so far, found to be lacking. This is partly due to the breakdown of complete frequency redistribution in the chromospheric layers and also because of the difficulty in obtaining complete sets of observations that adequately constrain the solar atmosphere at all relevant heights. Aims. We aim to obtain semi-empirical model atmospheres that reproduce the features of the Mg II h&k line profiles that sample the middle chromosphere with focus on a sunspot. Methods. We use spectropolarimetric observations of the Ca II 8542 A spectra obtained with the Swedish 1-m Solar Telescope (SST) and use NICOLE inversions to obtain semi-empirical model atmospheres for different features in and around a sunspot. These are used to synthesize Mg II h&k spectra using RH1.5D code, which we compare with observations taken with the Interface Region Imaging Spectrograph (IRIS). Results. Comparison of the synthetic profiles with IRIS observations reveals that there are several areas, especially in the penumbra of the sunspot, where most of the observed Mg II h&k profiles are very well reproduced. In addition, we find that supersonic hot downflows, present in our collection of models in the umbra, lead to synthetic profiles that agree well with the IRIS Mg II h&k profiles, with the exception of the line core. Conclusions. We put forward and make available four semi-empirical model atmospheres. Two for the penumbra, reflecting the range of temperatures obtained for the chromosphere, one for umbral flashes, and a model representative of the quiet surroundings of a sunspot. These are available in electronic as well as in table formats.

Ellerman bombs (EBs) and UV bursts are both brightenings related to flux emergence regions and specifically to magnetic flux of opposite polarity that meet in the photosphere. These two reconnection-related phenomena, nominally formed far apart, occasionally occur in the same location and at the same time, thus challenging our understanding of reconnection and heating of the lower solar atmosphere. We consider the formation of an active region, including long fibrils and hot and dense coronal plasma. The emergence of a untwisted magnetic flux sheet, injected $2.5$~Mm below the photosphere, is studied as it pierces the photosphere and interacts with the preexisting ambient field. Specifically, we aim to study whether EBs and UV bursts are generated as a result of such flux emergence and examine their physical relationship. The Bifrost radiative magnetohydrodynamics code was used to model flux emerging into a model atmosphere that contained a fairly strong ambient field, constraining the emerging field to a limited volume wherein multiple reconnection events occur as the field breaks through the photosphere and expands into the outer atmosphere. Synthetic spectra of the different reconnection events were computed using the $1.5$D RH code and the fully 3D MULTI3D code. The formation of UV bursts and EBs at intensities and with line profiles that are highly reminiscent of observed spectra are understood to be a result of the reconnection of emerging flux with itself in a long-lasting current sheet that extends over several scale heights through the chromosphere. Synthetic diagnostics suggest that there are no compelling reasons to assume that UV bursts occur in the photosphere. Instead, EBs and UV bursts are occasionally formed at opposite ends of a long current sheet that resides in an extended bubble of cool gas.

To date, no compelling evidence has been found that the universe has non-trivial spatial topology. Meanwhile, anomalies in the observed CMB temperature map, such as the lack of correlations at large angular separations, remain observationally robust. We show that if our universe is flat and has one compact dimension of appropriate size (slab topology), this would suppress large-angle temperature correlations while maintaining a low-$\ell$ angular power spectrum consistent with observations. The optimal length appears to be $1.4$ times the conformal radius of the CMB's last scattering surface ($\chi_{rec}$). We construct the probability distribution function of the statistic $S_{1/2}$ using simulated Sachs-Wolf-only skies for each of several values of $L_z/\chi_{rec}$. For $L_z\simeq1.4\chi_{rec}$ the $p$-value of four standard masked Planck maps is $p\simeq0.15$, compared to $p\lesssim0.003$ for the conventional topologically trivial space. The mean angular power spectrum $\langle C_{\ell} \rangle$ of the $L_z=1.4\chi_{rec}$ slab space matches the observed power spectrum at $2\leq\ell\lesssim6$ -- including a substantially suppressed quadrupole $C_2$, a slightly suppressed octopole $C_3$, and unsuppressed higher multipoles. It does not predict other low-$\ell$ CMB anomalies, and does not take account of normally sub-dominant Integrated Sachs Wolfe contributions. An $L_z=1.4\chi_{rec}$ slab topology is consistent with published limits from the Planck maps ($L_z\gtrsim1.12\chi_{rec}$). It is within the 95% confidence range $1.2\leq L_z/\chi_{rec}\leq2.1$ inferred using the covariance-matrix of temperature fluctuations. However, it violates published circles-in-the-sky limits from WMAP and related unpublished limits from Planck ($L_z/\chi_{rec}\gtrsim1.9$). We remark on the possibility to satisfy these limits, and "postdict" other large-angle anomalies, with closely related topologies.

The interface between the bright solar surface and the million-degree corona continues to hold the key to many unsolved problems in solar physics. Advances in instrumentation now allow us to observe the dynamic structures of the solar chromosphere down to less than 0.1" with cadences of just a few seconds and in multiple polarisation states. Such observational progress has been matched by the ever-increasing sophistication of numerical models, which have become necessary to interpret the complex observations. With an emphasis on the quiet Sun, I will review recent progress in the observation and modelling of the chromosphere. Models have come a long way from 1D static atmospheres, but their predictions still fail to reproduce several key observed features. Nevertheless, they have given us invaluable insight into the physical processes that energise the atmosphere. With more physics being added to models, the gap between predictions and observations is narrowing. With the next generation of solar observatories just around the corner, the big question is: will they close the gap?

High precision astrometry now enables to measure the time drift of astrophysical observables in real time, hence providing new ways to probe different cosmological models. This article presents a general derivation of the redshift and direction drifts for general observers. It is then applied to the standard cosmological framework of a Friedmann-Lemaitre space- time including all effects at first order in the cosmological perturbations, as well as in the class of spatially anisotropic universe models of the Bianchi I family. It shows that for a general observer, the direction drift splits into a parallax and an aberration drifts and order of magnitude estimates of these two components are provided. The multipolar decomposition of the redshift and aberration drifts is also derived and shows that the observer's peculiar velocity contributes only as a dipole whereas the anisotropic shear contributes as a quadrupole.

We present the results of 1D hydrodynamic simulations of coronal loops which are subject to nanoflares, caused by either in-situ thermal heating, or non-thermal electrons (NTE) beams. The synthesized intensity and Doppler shifts can be directly compared with IRIS and AIA observations of rapid variability in the transition region (TR) of coronal loops, associated with transient coronal heating. We find that NTE with high enough low-energy cutoff (E$_\textrm{C}$) deposit energy in the lower TR and chromosphere causing blueshifts (up to~$\sim$~20 km/s) in the \emph{IRIS} \siiv~lines, which thermal conduction cannot reproduce. The E$_\textrm{C}$ threshold value for the blueshifts depends on the total energy of the events ($\approx$~5 keV for 10$^{24}$ ergs, up to 15 keV for 10$^{25}$ ergs). The observed footpoint emission intensity and flows, combined with the simulations, can provide constraints on both the energy of the heating event and E$_\textrm{C}$. The response of the loop plasma to nanoflares depends crucially on the electron density: significant \siiv~intensity enhancements and flows are observed only for initially low-density loops ($<$~10$^{9}$~cm$^{-3}$). This provides a possible explanation of the relative scarcity of observations of significant moss variability. While the TR response to single heating episodes can be clearly observed, the predicted coronal emission (AIA 94\AA) for single strands is below current detectability, and can only be observed when several strands are heated closely in time. Finally, we show that the analysis of the IRIS \mgii~chromospheric lines can help further constrain the properties of the heating mechanisms.

Low-lying loops have been discovered at the solar limb in transition region temperatures by the Interface Region Imaging Spectrograph (IRIS). They do not appear to reach coronal temperatures, and it has been suggested that they are the long-predicted unresolved fine structures (UFS). These loops are dynamic and believed to be visible during both heating and cooling phases. Making use of coordinated observations between IRIS and the Swedish 1-m Solar Telescope, we study how these loops impact the solar chromosphere. We show for the first time that there is indeed a chromospheric signal of these loops, seen mostly in the form of strong Doppler shifts and a conspicuous lack of chromospheric heating. In addition, we find that several instances have a inverse Y-shaped jet just above the loop, suggesting that magnetic reconnection is driving these events. Our observations add several puzzling details to the current knowledge of these newly discovered structures; this new information must be considered in theoretical models.

We present results from an observing campaign in support of the VAULT2.0 sounding rocket launch on September 30, 2014. VAULT2.0 is a Ly$\alpha$ (1216 \AA) spectroheliograph capable of providing spectroheliograms at high cadence. Ly$\alpha$ observations are highly complementary to the IRIS observations of the upper chromosphere and the low transition region (TR) but have previously been unavailable. The VAULT2.0 data provide new constraints on upper-chromospheric conditions for numerical models. The observing campaign was closely coordinated with the IRIS mission. Taking advantage of this simultaneous multi-wavelength coverage of target AR 12172 and by using state-of-the-art radiative-MHD simulations of spicules, we investigate in detail a type-II spicule associated with a fast (300 km s$^{-1}$) network jet recorded in the campaign observations. Our analysis suggests that spicular material exists suspended high in the atmosphere but in lower temperatures (seen in Ly$\alpha$) until it is heated and becomes visible in TR temperatures as a network jet. The heating begins lower in the spicule and propagates upwards as a rapidly propagating thermal front. The front is then observed as fast, plane-of-the-sky motion typical of a network jet, but contained inside the pre-existing spicule. This work supports that the high speeds reported in network jets should not be taken as real mass upflows but only as apparent speeds of a rapidly propagating heating front along the pre-existing spicule.

In addition to their jet-like dynamic behaviour, spicules usually exhibit strong transverse speeds, multi-stranded structure and heating from chromospheric to transition region temperatures. In this work we first analyse \textit{Hinode} \& \textit{IRIS} observations of spicules and find different behaviours in terms of their Doppler velocity evolution and collective motion of their sub-structure. Some have a Doppler shift sign change that is rather fixed along the spicule axis, and lack coherence in the oscillatory motion of strand-like structure, matching rotation models or long wavelength torsional Alfv\'en waves. Others exhibit a Doppler shift sign change at maximum displacement and coherent motion of their strands, suggesting a collective MHD wave. By comparing with an idealised 3-D MHD simulation combined with radiative transfer modelling, we analyse the role of transverse MHD waves and associated instabilities in spicule-like features. We find that Transverse Wave Induced Kelvin-Helmholtz (TWIKH) rolls lead to coherence of strand-like structure in imaging and spectral maps, as seen in some observations. The rapid transverse dynamics and the density and temperature gradients at the spicule boundary lead to ring-shaped \ion{Mg}{2} k and \ion{Ca}{2} H source functions in the transverse cross-section, potentially allowing IRIS to capture the KHI dynamics. Twists and currents propagate along the spicule at Alfv\'enic speeds, and the temperature variations within TWIKH rolls produce sudden appearance / disappearance of strands seen in Doppler velocity and in \ion{Ca}{2} H intensity. However, only a mild intensity increase in higher temperature lines is obtained, suggesting there is an additional heating mechanism at work in spicules.

Magnetic reconnection is thought to drive a wide variety of dynamic phenomena in the solar atmosphere. Yet the detailed physical mechanisms driving reconnection are difficult to discern in the remote sensing observations that are used to study the solar atmosphere. In this paper we exploit the high-resolution instruments Interface Region Imaging Spectrograph (IRIS) and the new CHROMIS Fabry-Perot instrument at the Swedish 1-m Solar Telescope (SST) to identify the intermittency of magnetic reconnection and its association with the formation of plasmoids in so-called UV bursts in the low solar atmosphere. The Si IV 1403A UV burst spectra from the transition region show evidence of highly broadened line profiles with often non-Gaussian and triangular shapes, in addition to signatures of bidirectional flows. Such profiles had previously been linked, in idealized numerical simulations, to magnetic reconnection driven by the plasmoid instability. Simultaneous CHROMIS images in the chromospheric Ca II K 3934A line now provide compelling evidence for the presence of plasmoids, by revealing highly dynamic and rapidly moving brightenings that are smaller than 0.2 arcsec and that evolve on timescales of order seconds. Our interpretation of the observations is supported by detailed comparisons with synthetic observables from advanced numerical simulations of magnetic reconnection and associated plasmoids in the chromosphere. Our results highlight how subarcsecond imaging spectroscopy sensitive to a wide range of temperatures combined with advanced numerical simulations that are realistic enough to compare with observations can directly reveal the small-scale physical processes that drive the wide range of phenomena in the solar atmosphere.

In the lower solar atmosphere, the chromosphere is permeated by jets known as spicules, in which plasma is propelled at speeds of 50 to 150 kilometers per second into the corona. The origin of the spicules is poorly understood, although they are expected to play a role in heating the million-degree corona and are associated with Alfv\'enic waves that help drive the solar wind. We compare magnetohydrodynamic simulations of spicules with observations from the Interface Region Imaging Spectrograph and the Swedish 1-m Solar Telescope. Spicules are shown to occur when magnetic tension is amplified and transported upward through interactions between ions and neutrals or ambipolar diffusion. The tension is impulsively released to drive flows, heat plasma (through ambipolar diffusion), and generate Alfv\'enic waves.

Spectral observations of the solar transition region (TR) and corona show broadening of spectral lines beyond what is expected from thermal and instrumental broadening. The remaining non-thermal broadening is significant (5-30 km/s) and correlated with the intensity. Here we study spectra of the TR Si IV 1403 Angstrom line obtained at high resolution with the Interface Region Imaging Spectrograph (IRIS). We find that the large improvement in spatial resolution (0.33 arcsec) of IRIS compared to previous spectrographs (2 arcsec) does not resolve the non-thermal line broadening which in most regions remains at pre-IRIS levels of about 20 km/s. This invariance to spatial resolution indicates that the processes behind the broadening occur along the line-of-sight (LOS) and/or on spatial scales (perpendicular to the LOS) smaller than 250 km. Both effects appear to play a role. Comparison with IRIS chromospheric observations shows that, in regions where the LOS is more parallel to the field, magneto-acoustic shocks driven from below impact the TR and can lead to significant non-thermal line broadening. This scenario is supported by MHD simulations. While these do not show enough non-thermal line broadening, they do reproduce the long-known puzzling correlation between non-thermal line broadening and intensity. This correlation is caused by the shocks, but only if non-equilibrium ionization is taken into account. In regions where the LOS is more perpendicular to the field, the prevalence of small-scale twist is likely to play a significant role in explaining the invariance and correlation with intensity.

A spectacular manifestation of solar activity is the appearance of transient brightenings in the far wings of the H$\alpha$ line, known as Ellerman bombs (EBs). Recent observations obtained by the Interface Region Imaging Spectrograph (IRIS) have revealed another type of plasma "bombs" (UV bursts) with high temperatures of perhaps up to $8\times 10^4$ K within the cooler lower solar atmosphere. Realistic numerical modeling showing such events is needed to explain their nature. Here, we report on 3D radiative magneto-hydrodynamic simulations of magnetic flux emergence in the solar atmosphere. We find that ubiquitous reconnection between emerging bipolar magnetic fields can trigger EBs in the photosphere, UV bursts in the mid/low chromosphere and small (nano-/micro-) flares ($10^6$ K) in the upper chromosphere. These results provide new insights on the emergence and build up of the coronal magnetic field and the dynamics and heating of the solar surface and lower atmosphere.

Our ability to model the shapes and strengths of iron lines in the solar spectrum is a critical test of the accuracy of the solar iron abundance, which sets the absolute zero-point of all stellar metallicities. We use an extensive 463-level Fe atom with new photoionisation cross-sections for FeI as well as quantum mechanical calculations of collisional excitation and charge transfer with neutral hydrogen; the latter effectively remove a free parameter that has hampered all previous line formation studies of Fe in non-local thermodynamic equilibrium (NLTE). For the first time, we use realistic 3D NLTE calculations of Fe for a quantitative comparison to solar observations. We confront our theoretical line profiles with observations taken at different viewing angles across the solar disk with the Swedish 1-m Solar Telescope. We find that 3D modelling well reproduces the observed centre-to-limb behaviour of spectral lines overall, but highlight aspects that may require further work, especially cross-sections for inelastic collisions with electrons. Our inferred solar iron abundance is log(eps(Fe))=7.48+-0.04.

Two-point correlation functions are ubiquitous tools of modern cosmology, appearing in disparate topics ranging from cosmological inflation to late-time astrophysics. When the background spacetime is maximally symmetric, invariance arguments can be used to fix the functional dependence of this function as the invariant distance between any two points. In this paper we introduce a novel formalism which fixes this functional dependence directly from the isometries of the background metric, thus allowing one to quickly assess the overall features of Gaussian correlators without resorting to the full machinery of perturbation theory. As an application we construct the CMB temperature correlation function in one inhomogeneous (namely, an off-center LTB model) and two spatially flat and anisotropic (Bianchi) universes, and derive their covariance matrices in the limit of almost Friedmannian symmetry. We show how the method can be extended to arbitrary N -point correlation functions and illustrate its use by constructing three-point correlation functions in some simple geometries.

This article details the computation of the two-point correlators of the convergence, $E$- and $B$-modes of the cosmic shear induced by the weak-lensing by large scale structure assuming that the background spacetime is spatially homogeneous and anisotropic. After detailing the perturbation equations and the general theory of weak-lensing in an anisotropic universe, it develops a weak shear approximation scheme in which one can compute analytically the evolution of the Jacobi matrix. It allows one to compute the angular power spectrum of the $E$- and $B$-modes. In the linear regime, the existence of $B$-modes is a direct tracer of a late time anisotropy and their angular power spectrum scales as the square of the shear. It is then demonstrated that there must also exist off-diagonal correlations between the $E$-modes, $B$-modes and convergence that are linear in the geometrical shear and allow one to reconstruct the eigendirections of expansion. These spectra can be measured in future large scale surveys, such as Euclid and SKA, and offer a new tool to test the isotropy of the expansion of the universe at low redshift.

We compute the angular power spectrum of the $B$-modes of the weak-lensing shear in a spatially anisotropic spacetime. We find that there must also exist off-diagonal correlations between the $E$-modes, $B$-modes, and convergence that allow one to reconstruct the eigendirections of expansion. Focusing on future surveys such as Euclid and SKA, we show that observations can constrain the geometrical shear in units of the Hubble rate at the percent level, or even better, offering a new and powerful method to probe our cosmological model.

We present the scientific case for a future space-based telescope aimed at very high spatial and temporal resolution imaging of the solar photosphere and chromosphere. Previous missions (e.g., HINODE, SUNRISE) have demonstrated the power of observing the solar photosphere and chromosphere at high spatial resolution without contamination from Earth's atmosphere. We argue here that increased spatial resolution (from currently 70 km to 25 km in the future) and high temporal cadence of the observations will vastly improve our understanding of the physical processes controlling solar magnetism and its characteristic scales. This is particularly important as the Sun's magnetic field drives solar activity and can significantly influence the Sun-Earth system. At the same time a better knowledge of solar magnetism can greatly improve our understanding of other astrophysical objects.

Solar spicules are chromospheric fibrils that appear everywhere on the Sun, yet their origin is not understood. Using high resolution observations of spicules obtained with the Swedish 1-m Solar Telescope, we aim to understand how spicules appear in filtergrams and Dopplergrams, how they compare in Ca II H and H-alpha, and what can make them appear and disappear. We find thatspicules display a rich and detailed spatial structure, and show a distribution of transverse velocities that when aligned with the line of sight can make them appear at different H-alpha wing positions. They become more abundant at positions closer to the line core, reflecting a distribution of Doppler shifts and widths. In H-alpha width maps they stand out as bright features both on disk and off limb, reflecting their large Doppler motions and possibly higher temperatures than in the typical H-alpha formation region. Spicule lifetimes measured from narrowband images at only a few positions will be an underestimate because Doppler shifts can make them disappear prematurely from such images; for such cases, width maps are a more robust tool. In H-alpha and Ca II H filtergrams, off limb spicules essentially have the same properties, appearance, and evolution. We find that the sudden appearance of spicules can be explained by Doppler shifts from their transverse motions, and does not require other convoluted explanations.

If the spacetime metric has anisotropic spatial curvature, one can afford to expand the universe isotropically, provided that the energy-momentum tensor satisfy a certain con- straint. This leads to the so-called shear-free metrics, which have the interesting property of violating the cosmological principle while still preserving the isotropy of the cosmic mi- crowave background (CMB) radiation. In this work we show that shear-free cosmologies correspond to an attractor solution in the space of models with anisotropic spatial curva- ture. Through a rigorous definition of linear perturbation theory in these spacetimes, we show that shear-free models represent a viable alternative to describe the large-scale evo- lution of the universe, leading, in particular, to a kinematically equivalent Sachs-Wolfe effect. Alternatively, we discuss some specific signatures that shear-free models would imprint on the temperature spectrum of CMB.

Despite its continued observational successes, there is a persistent (and growing) interest in extending cosmology beyond the standard model, $\Lambda$CDM. This is motivated by a range of apparently serious theoretical issues, involving such questions as the cosmological constant problem, the particle nature of dark matter, the validity of general relativity on large scales, the existence of anomalies in the CMB and on small scales, and the predictivity and testability of the inflationary paradigm. In this paper, we summarize the current status of $\Lambda$CDM as a physical theory, and review investigations into possible alternatives along a number of different lines, with a particular focus on highlighting the most promising directions. While the fundamental problems are proving reluctant to yield, the study of alternative cosmologies has led to considerable progress, with much more to come if hopes about forthcoming high-precision observations and new theoretical ideas are fulfilled.

A primordial inflationary phase allows one to erase any possible anisotropic expansion thanks to the cosmic no-hair theorem. If there is no global anisotropic stress, then the anisotropic expansion rate tends to decrease. What are the observational consequences of a possible early anisotropic phase? We first review the dynamics of anisotropic universes and report analytic approximations. We then discuss the structure of dynamical equations for perturbations and the statistical properties of observables, as well as the implication of a primordial anisotropy on the quantization of these perturbations during inflation. Finally we briefly review models based on primordial vector field which evade the cosmic no-hair theorem.

The C II 133.5 nm multiplet has been observed by NASA's Interface Region Imaging Spectrograph (IRIS) in unprecedented spatial resolution. The aims of this work are to characterize these new observations of the C II lines, place them in context with previous work, and to identify any additional value the C II lines bring when compared with other spectral lines. We make use of wide, long exposure IRIS rasters covering the quiet Sun and an active region. Line properties such as velocity shift and width are extracted from individual spectra and analyzed. The lines have a variety of shapes (mostly single-peak or double-peak), are strongest in active regions and weaker in the quiet Sun. The ratio between the 133.4 nm and 133.5 nm components is always less than 1.8, indicating that their radiation is optically thick in all locations. Maps of the C II line widths are a powerful new diagnostic of chromospheric structures, and their line shifts are a robust velocity diagnostic. Compared with earlier quiet Sun observations, we find similar absolute intensities and mean line widths, but smaller red shifts; this difference can perhaps be attributed to differences in spectral resolution and spatial coverage. The C II intensity maps are somewhat similar to those of transition region lines, but also share some features with chromospheric maps such as those from the Mg II k line, indicating that they are formed between the upper chromosphere and transition region. C II intensity, width, and velocity maps can therefore be used to gather additional information about the upper chromosphere.

When searching for deviations of statistical isotropy in CMB, a popular strategy is to write the two-point correlation function (2pcf) as the most general function of four spherical angles (i.e., two unit vectors) in the celestial sphere. Then, using a basis of bipolar spherical harmonics, statistical anisotropy will show up if and only if any coefficient of the expansion with non-trivial bipolar momentum is detected -- although this detection will not in general elucidate the origin of the anisotropy. In this work we show that two new sets of four angles and basis functions exist which completely specifies the 2pcf, while, at the same time, offering a clearer geometrical interpretation of the mechanisms generating the signal. Since the coefficients of these expansions are zero if and only if isotropy holds, they act as a simple and geometrically motivated null test of statistical isotropy, with the advantage of allowing cosmic variance to be controlled in a systematic way. We report the results of the application of these null tests to the latest temperature data released by the Planck collaboration.

If one is willing to give up the cherished hypothesis of spatial isotropy, many interesting cosmological models can be developed beyond the simple anisotropically expanding scenarios. One interesting possibility is presented by shear-free models in which the anisotropy emerges at the level of the curvature of the homogeneous spatial sections, whereas the expansion is dictated by a single scale factor. We show that such models represent viable alternatives to describe the large-scale structure of the inflationary universe, leading to a kinematically equivalent Sachs-Wolfe effect. Through the definition of a complete set of spatial eigenfunctions we compute the two-point correlation function of scalar perturbations in these models. In addition, we show how such scenarios would modify the spectrum of the CMB assuming that the observations take place in a small patch of a universe with anisotropic curvature.

Spicules are ubiquitous, fast moving jets observed off-limb in chromospheric spectral lines. Combining the recently-launched Interface Region Imaging Spectrograph with the Solar Dynamics Observatory and Hinode, we have a unique opportunity to study spicules simultaneously in multiple passbands and from a seeing free environment. This makes it possible to study their thermal evolution over a large range of temperatures. A recent study showed that spicules appear in several chromospheric and transition region spectral lines, suggesting that spicules continue their evolution in hotter passbands after they fade from Ca II H. In this follow-up paper we answer some of the questions that were raised in the introductory study. In addition, we study spicules off-limb in C II 1330 {\AA} for the first time. We find that Ca II H spicules are more similar to Mg II 2796 {\AA} spicules than initially reported. For a sample of 54 spicules, we find that 44% of Si IV 1400 {\AA} spicules are brighter toward the top; 56% of the spicules show an increase in Si IV emission when the Ca II H component fades. We find several examples of spicules that fade from passbands other than Ca II H, and we observe that if a spicule fades from a passband, it also generally fades from the passbands with lower formation temperatures. We discuss what these new, multi-spectral results mean for the classification of type I and II spicules.

In this work, we analyse coordinated observations spanning chromospheric, TR and coronal temperatures at very high resolution which reveal essential characteristics of thermally unstable plasmas. Coronal rain is found to be a highly multi-thermal phenomenon with a high degree of co-spatiality in the multi-wavelength emission. EUV darkening and quasi-periodic intensity variations are found to be strongly correlated to coronal rain showers. Progressive cooling of coronal rain is observed, leading to a height dependence of the emission. A fast-slow two-step catastrophic cooling progression is found, which may reflect the transition to optically thick plasma states. The intermittent and clumpy appearance of coronal rain at coronal heights becomes more continuous and persistent at chromospheric heights just before impact, mainly due to a funnel effect from the observed expansion of the magnetic field. Strong density inhomogeneities on spatial scales of 0.2"-0.5" are found, in which TR to chromospheric temperature transition occurs at the lowest detectable scales. The shape of the distribution of coronal rain widths is found to be independent of temperature with peaks close to the resolution limit of each telescope, ranging from 0.2" to 0.8". However we find a sharp increase of clump numbers at the coolest wavelengths and especially at higher resolution, suggesting that the bulk of the rain distribution remains undetected. Rain clumps appear organised in strands in both chromospheric and TR temperatures, suggesting an important role of thermal instability in the shaping of fundamental loop substructure. We further find structure reminiscent of the MHD thermal mode. Rain core densities are estimated to vary between 2x10^{10} cm^{-3} and 2.5x10^{11} cm^{-3} leading to significant downward mass fluxes per loop of 1-5x10^{9} g s^{-1}, suggesting a major role in the chromosphere-corona mass cycle.

A triplet of subordinate lines of Mg II exists in the region around the h&k lines. In solar spectra these lines are seen mostly in absorption, but in some cases can become emission lines. The aim of this work is to study the formation of this triplet, and investigate any diagnostic value they can bring. Using 3D radiative magnetohydrodynamic simulations of quiet Sun and flaring flux emergence, we synthesize spectra and investigate how spectral features respond to the underlying atmosphere. We find that emission in the lines is rare and is typically caused by a steep temperature increase in the lower chromosphere (above 1500 K, with electron densities above 10$^{18}$ m$^{-3}$). In both simulations the lines are sensitive to temperature increases taking place at column masses >= 5e-4 g cm$^{-2}$. Additional information can also be inferred from the peak-to-wing ratio and shape of the line profiles. Using observations from NASA's Interface Region Imaging Spectrograph we find both absorption and emission line profiles with similar shapes to the synthetic spectra, which suggests that these lines represent a useful diagnostic that complements the MgII h&k lines.

The Interface Region Imaging Spectrograph (IRIS) reveals small-scale rapid brightenings in the form of bright grains all over coronal holes and the quiet sun. These bright grains are seen with the IRIS 1330 \AA, 1400 \AA\ and 2796 \AA\ slit-jaw filters. We combine coordinated observations with IRIS and from the ground with the Swedish 1-m Solar Telescope (SST) which allows us to have chromospheric (Ca II 8542 \AA, Ca II H 3968 \AA, H\alpha, and Mg II k 2796 \AA), and transition region (C II 1334 \AA, Si IV 1402) spectral imaging, and single-wavelength Stokes maps in Fe I 6302 \AA at high spatial (0.33"), temporal and spectral resolution. We conclude that the IRIS slit-jaw grains are the counterpart of so-called acoustic grains, i.e., resulting from chromospheric acoustic waves in a non-magnetic environment. We compare slit-jaw images with spectra from the IRIS spectrograph. We conclude that the grain intensity in the 2796 \AA\ slit-jaw filter comes from both the Mg II k core and wings. The signal in the C II and Si IV lines is too weak to explain the presence of grains in the 1300 and 1400 \AA\ slit-jaw images and we conclude that the grain signal in these passbands comes mostly from the continuum. Even though weak, the characteristic shock signatures of acoustic grains can often be detected in IRIS C II spectra. For some grains, spectral signature can be found in IRIS Si IV. This suggests that upward propagating acoustic waves sometimes reach all the way up to the transition region.

The emergence of three-dimensional magneto-hydrodynamic (MHD) simulations of stellar atmospheres has sparked a need for efficient radiative transfer codes to calculate detailed synthetic spectra. We present RH 1.5D, a massively parallel code based on the RH code and capable of performing Zeeman polarised multi-level non-local thermodynamical equilibrium (NLTE) calculations with partial frequency redistribution for an arbitrary amount of chemical species. The code calculates spectra from 3D, 2D or 1D atmospheric models on a column-by-column basis (or 1.5D). While the 1.5D approximation breaks down in the cores of very strong lines in an inhomogeneous environment, it is nevertheless suitable for a large range of scenarios and allows for faster convergence with finer control over the iteration of each simulation column. The code scales well to at least tens of thousands of CPU cores, and is publicly available. In the present work we briefly describe its inner workings, strategies for convergence optimisation, its parallelism, and some possible applications.

We use coordinated observations with the Interface Region Imaging Spectrograph (IRIS) and the Swedish 1-m Solar Telescope (SST) to identify the disk counterpart of type II spicules in upper-chromospheric and transition region (TR) diagnostics. These disk counterparts were earlier identified through short-lived asymmetries in chromospheric spectral lines: rapid blue- or red-shifted excursions (RBEs or RREs). We find clear signatures of RBEs and RREs in Mg II h & k, often with excursions of the central h3 and k3 absorption features in concert with asymmetries in co-temporal and co-spatial H-alpha spectral profiles. We find spectral signatures for RBEs and RREs in C II 1335 and 1336 A and Si IV 1394 and 1403 A spectral lines and interpret this as a sign that type II spicules are heated to at least TR temperatures, supporting other recent work. These C II and Si IV spectral signals are weaker for a smaller network region than for more extended network regions in our data. A number of bright features around extended network regions observed in IRIS slit-jaw imagery SJI 1330 and 1400, recently identified as network jets, can be clearly connected to H-alpha RBEs and/or RREs in our coordinated data. We speculate that at least part of the diffuse halo around network regions in the IRIS SJI 1330 and 1400 images can be attributed to type II spicules with insufficient opacity in the C II and Si IV lines to stand out as single features in these passbands.

The heating of the outer solar atmospheric layers, i.e., the transition region and corona, to high temperatures is a long standing problem in solar (and stellar) physics. Solutions have been hampered by an incomplete understanding of the magnetically controlled structure of these regions. The high spatial and temporal resolution observations with the Interface Region Imaging Spectrograph (IRIS) at the solar limb reveal a plethora of short, low lying loops or loop segments at transition-region temperatures that vary rapidly, on the timescales of minutes. We argue that the existence of these loops solves a long standing observational mystery. At the same time, based on comparison with numerical models, this detection sheds light on a critical piece of the coronal heating puzzle.

Photon spheres, surfaces where massless particles are confined in closed orbits, are expected to be common astrophysical structures surrounding ultracompact objects. In this paper a semiclassical treatment of a photon sphere is proposed. We consider the quantum Maxwell field and derive its energy spectra. A thermodynamic approach for the quantum photon sphere is developed and explored. Within this treatment, an expression for the spectral energy density of the emitted radiation is presented. Our results suggest that photon spheres, when thermalized with their environment, have nonusual thermodynamic properties, which could lead to distinct observational signatures.

The solar chromosphere and transition region (TR) form an interface between the Sun's surface and its hot outer atmosphere. Here most of the non-thermal energy that powers the solar atmosphere is transformed into heat, although the detailed mechanism remains elusive. High-resolution (0.33-arcsec) observations with NASA's Interface Region Imaging Spectrograph (IRIS) reveal a chromosphere and TR that are replete with twist or torsional motions on sub-arcsecond scales, occurring in active regions, quiet Sun regions, and coronal holes alike. We coordinated observations with the Swedish 1-m Solar Telescope (SST) to quantify these twisting motions and their association with rapid heating to at least TR temperatures. This view of the interface region provides insight into what heats the low solar atmosphere.

We report the results of a statistical analysis performed with the four foreground-cleaned Planck maps by means of a suitably defined {\it local-variance} estimator. Our analysis shows a clear dipolar structure in Planck's {\it variance map} pointing in the direction $(l,b) \,\simeq\, (220^{\circ},-32^{\circ})$, thus consistent with the North-South asymmetry phenomenon. Surprisingly, and contrary to previous findings, removing the CMB quadrupole and octopole makes the asymmetry stronger. Our results show a maximal statistical significance, of $98.1\%$ CL, in the scales ranging from $\ell=4$ to $\ell=500$. Additionally, through exhaustive analyses of the four foreground-cleaned and individual frequency Planck maps, we find unlikely that residual foregrounds could be causing this dipole variance asymmetry. Moreover, we find that the dipole gets lower amplitudes for larger masks, evidencing that most of the contribution to the variance dipole comes from a region near the galactic plane. Finally, our results are robust against different foreground cleaning procedures, different Planck masks, pixelization parameters, and the addition of inhomogeneous real noise.

Solar spicules have eluded modelers and observers for decades. Since the discovery of the more energetic type II, spicules have become a heated topic but their contribution to the energy balance of the low solar atmosphere remains unknown. Here we give a first glimpse of what quiet Sun spicules look like when observed with NASA's recently launched Interface Region Imaging Spectrograph (IRIS). Using IRIS spectra and filtergrams that sample the chromosphere and transition region we compare the properties and evolution of spicules as observed in a coordinated campaign with Hinode and the Atmospheric Imaging Assembly. Our IRIS observations allow us to follow the thermal evolution of type II spicules and finally confirm that the fading of Ca II H spicules appears to be caused by rapid heating to higher temperatures. The IRIS spicules do not fade but continue evolving, reaching higher and falling back down after 500-800 s. Ca II H type II spicules are thus the initial stages of violent and hotter events that mostly remain invisible in Ca II H filtergrams. These events have very different properties from type I spicules, which show lower velocities and no fading from chromospheric passbands. The IRIS spectra of spicules show the same signature as their proposed disk counterparts, reinforcing earlier work. Spectroheliograms from spectral rasters also confirm that quiet Sun spicules originate in bushes from the magnetic network. Our results suggest that type II spicules are indeed the site of vigorous heating (to at least transition region temperatures) along extensive parts of the upward moving spicular plasma.

The Interface Region Imaging Spectrograph (IRIS) small explorer spacecraft provides simultaneous spectra and images of the photosphere, chromosphere, transition region, and corona with 0.33-0.4 arcsec spatial resolution, 2 s temporal resolution and 1 km/s velocity resolution over a field-of-view of up to 175 arcsec x 175 arcsec. IRIS was launched into a Sun-synchronous orbit on 27 June 2013 using a Pegasus-XL rocket and consists of a 19-cm UV telescope that feeds a slit-based dual-bandpass imaging spectrograph. IRIS obtains spectra in passbands from 1332-1358, 1389-1407 and 2783-2834 Angstrom including bright spectral lines formed in the chromosphere (Mg II h 2803 Angstrom and Mg II k 2796 Angstrom) and transition region (C II 1334/1335 Angstrom and Si IV 1394/1403 Angstrom). Slit-jaw images in four different passbands (C II 1330, Si IV 1400, Mg II k 2796 and Mg II wing 2830 Angstrom) can be taken simultaneously with spectral rasters that sample regions up to 130 arcsec x 175 arcsec at a variety of spatial samplings (from 0.33 arcsec and up). IRIS is sensitive to emission from plasma at temperatures between 5000 K and 10 MK and will advance our understanding of the flow of mass and energy through an interface region, formed by the chromosphere and transition region, between the photosphere and corona. This highly structured and dynamic region not only acts as the conduit of all mass and energy feeding into the corona and solar wind, it also requires an order of magnitude more energy to heat than the corona and solar wind combined. The IRIS investigation includes a strong numerical modeling component based on advanced radiative-MHD codes to facilitate interpretation of observations of this complex region. Approximately eight Gbytes of data (after compression) are acquired by IRIS each day and made available for unrestricted use within a few days of the observation.

The Mg II h&k lines are the prime chromospheric diagnostics of NASA's Interface Region Imaging Spectrograph (IRIS). In the previous papers of this series we used a realistic three-dimensional radiative magnetohydrodynamics model to calculate the h&k lines in detail and investigated how their spectral features relate to the underlying atmosphere. In this work, we employ the same approach to investigate how the h&k diagnostics fare when taking into account the finite resolution of IRIS and different noise levels. In addition, we investigate the diagnostic potential of several other photospheric lines and near-continuum regions present in the near-ultraviolet (NUV) window of IRIS and study the formation of the NUV slit-jaw images. We find that the instrumental resolution of IRIS has a small effect on the quality of the h&k diagnostics; the relations between the spectral features and atmospheric properties are mostly unchanged. The peak separation is the most affected diagnostic, but mainly due to limitations of the simulation. The effects of noise start to be noticeable at a signal-to-noise ratio (S/N) of 20, but we show that with noise filtering one can obtain reliable diagnostics at least down to a S/N of 5. The many photospheric lines present in the NUV window provide velocity information for at least eight distinct photospheric heights. Using line-free regions in the h&k far wings we derive good estimates of photospheric temperature for at least three heights. Both of these diagnostics, in particular the latter, can be obtained even at S/Ns as low as 5.

The effect of the scalar spectral index on inflationary super-Hubble waves is to amplify/damp large wavelengths according to whether the spectrum is red ($n_{s}<1$) or blue ($n_{s}>1$). As a consequence, the large-scale temperature correlation function will unavoidably change sign at some angle if our spectrum is red, while it will always be positive if it is blue. We show that this inflationary filtering property also affects our estimates of the size of the homogeneous patch of the universe through the Grishchuk-Zel'dovich effect. Using the recent quadrupole measurement of ESA's Planck mission, we find that the homogeneous patch of universe is at least 87 times bigger than our visible universe if we accept Planck's best fit value $n_{s}=0.9624$. An independent estimation of the size of the universe could be used to independently constrain $n_{s}$, thus narrowing the space of inflationary models.

NASA's Interface Region Imaging Spectrograph (IRIS) small explorer mission will study how the solar atmosphere is energized. IRIS contains an imaging spectrograph that covers the Mg II h&k lines as well as a slit-jaw imager centered at Mg II k. Understanding the observations requires forward modeling of Mg II h&k line formation from 3D radiation-MHD models. We compute the vertically emergent h&k intensity from a snapshot of a dynamic 3D radiation-MHD model of the solar atmosphere, and investigate which diagnostic information about the atmosphere is contained in the synthetic line profiles. We find that the Doppler shift of the central line depression correlates strongly with the vertical velocity at optical depth unity, which is typically located less than 200 km below the transition region (TR). By combining the Doppler shifts of the h and the k line we can retrieve the sign of the velocity gradient just below the TR. The intensity in the central line depression is anticorrelated with the formation height, especially in subfields of a few square Mm. This intensity could thus be used to measure the spatial variation of the height of the transition region. The intensity in the line-core emission peaks correlates with the temperature at its formation height, especially for strong emission peaks. The peaks can thus be exploited as a temperature diagnostic. The wavelength difference between the blue and red peaks provides a diagnostic of the velocity gradients in the upper chromosphere. The intensity ratio of the blue and red peaks correlates strongly with the average velocity in the upper chromosphere. We conclude that the Mg II h&k lines are excellent probes of the very upper chromosphere just below the transition region, a height regime that is impossible to probe with other spectral lines.

NASA's Interface Region Imaging Spectrograph (IRIS) space mission will study how the solar atmosphere is energized. IRIS contains an imaging spectrograph that covers the Mg II h&k lines as well as a slit-jaw imager centered at Mg II k. Understanding the observations will require forward modeling of Mg II h&k line formation from 3D radiation-MHD models. This paper is the first in a series where we undertake this forward modeling. We discuss the atomic physics pertinent to h&k line formation, present a quintessential model atom that can be used in radiative transfer computations and discuss the effect of partial redistribution (PRD) and 3D radiative transfer on the emergent line profiles. We conclude that Mg II h&k can be modeled accurately with a 4-level plus continuum Mg II model atom. Ideally radiative transfer computations should be done in 3D including PRD effects. In practice this is currently not possible. A reasonable compromise is to use 1D PRD computations to model the line profile up to and including the central emission peaks, and use 3D transfer assuming complete redistribution to model the central depression.

We performed a 3D radiative MHD simulation of the solar atmosphere. This simulation shows a jet-like feature that shows similarities to the type II spicules observed for the first time with Hinode. Rapid Blueshifted Events (RBEs) on the solar disk are associated with these spicules. Observational results suggest they may contribute significantly in supplying the corona with hot plasma. We perform a detailed comparison of the properties of the simulated jet with those of type II spicules (observed with Hinode) and RBEs (with ground-based instruments). We analyze variety of synthetic emission and absorption lines from the simulations including chromospheric Ca II and Ha to TR and coronal temperatures (10E4 to several 10E6K). We compare their synthetic intensities, line profiles, Doppler shifts, line widths and asymmetries with observations from Hinode/SOT and EIS, SOHO/SUMER, SST and SDO/AIA. Many properties of the synthetic observables resemble the observations, and we describe in detail the physical processes that lead to these observables. Detailed analysis of the synthetic observables provides insight into how observations should be analyzed to derive information about physical variables in such a dynamic event. For example, we find that LOS superposition in the optically thin atmosphere requires the combination of Doppler shifts and spectral line asymmetry to determine the velocity in the jet. Other properties differ from the observations, especially in the chromospheric lines. The mass density of the part of the spicule with a chromospheric temperature is too low to produce significant opacity in chromospheric lines. These and other discrepancies are described in detail, and we discuss which mechanisms and physical processes may need to be included in the MHD simulations to mimic the thermodynamic processes of the chromosphere and corona, in particular to reproduce type II spicules.

Recently, new solar model atmospheres have been developed to replace classical 1D LTE hydrostatic models and used to for example derive the solar chemical composition. We aim to test various models against key observational constraints. In particular, a 3D model used to derive the solar abundances, a 3D MHD model (with an imposed 10 mT vertical magnetic field), 1D models from the PHOENIX project, the 1D MARCS model, and the 1D semi-empirical model of Holweger & M\"uller. We confront the models with observational diagnostics of the temperature profile: continuum centre-to-limb variations (CLV), absolute continuum fluxes, and the wings of hydrogen lines. We also test the 3D models for the intensity distribution of the granulation and spectral line shapes. The predictions from the 3D model are in excellent agreement with the continuum CLV observations, performing even better than the Holweger & M\"uller model (constructed largely to fulfil such observations). The predictions of the 1D theoretical models are worse, given their steeper temperature gradients. For the continuum fluxes, predictions for most models agree well with the observations. No model fits all hydrogen lines perfectly, but again the 3D model comes ahead. The 3D model also reproduces the observed continuum intensity fluctuations and spectral line shapes very well. The excellent agreement of the 3D model with the observables reinforces the view that its temperature structure is realistic. It outperforms the MHD simulation in all diagnostics, implying that recent claims for revised abundances based on MHD modelling are premature. Several weaknesses in the 1D models are exposed. The differences between the PHOENIX LTE and NLTE models are small. We conclude that the 3D hydrodynamical model is superior to any of the tested 1D models, which gives further confidence in the solar abundance analyses based on it.

This article derives a multipolar hierarchy for the propagation of the weak-lensing shear and convergence in a general spacetime. The origin of B-modes, in particular on large angular scales, is related to the local isotropy of space. Known results assuming a Friedmann-Lema\^itre background are naturally recovered. The example of a Bianchi I spacetime illustrates our formalism and its implications for future observations are stressed.

Spicules have been observed on the sun for more than a century, typically in chromospheric lines such as H-alpha and Ca II H. Recent work has shown that so-called 'type II' spicules may have a role in providing mass to the corona and the solar wind. In chromospheric filtergrams these spicules are not seen to fall back down, and they are shorter-lived and more dynamic than the spicules that have been classically reported in ground-based observations. Observations of type II spicules with Hinode show fundamentally different properties from what was previously measured. In earlier work we showed that these dynamic type II spicules are the most common type, a view that was not properly identified by early observations.The aim of this work is to investigate the effects of spatio-temporal resolution in the classical spicule measurements. Making use of Hinode data degraded to match the observing conditions of older ground-based studies, we measure the properties of spicules with a semi-automated algorithm. These results are then compared to measurements using the original Hinode data. We find that degrading the data has a significant effect on the measured properties of spicules. Most importantly, the results from the degraded data agree well with older studies (e.g. mean spicule duration more than 5 minutes, and upward apparent velocities of about 25 km/s). These results illustrate how the combination of spicule superposition, low spatial resolution and cadence affect the measured properties of spicules, and that previous measurements can be misleading.

The effect of our Galaxy's motion through the Cosmic Microwave Background rest frame, which aberrates and Doppler shifts incoming photons measured by current CMB experiments, has been shown to produce mode-mixing in the multipole space temperature coefficients. However, multipole space determinations are subject to many difficulties, and a real-space analysis can provide a straightforward alternative. In this work we describe a numerical method for removing Lorentz- boost effects from real-space temperature maps. We show that to deboost a map so that one can accurately extract the temperature power spectrum requires calculating the boost kernel at a finer pixelization than one might naively expect. In idealized cases that allow for easy comparison to analytic results, we have confirmed that there is indeed mode mixing among the spherical harmonic coefficients of the temperature. We find that using a boost kernel calculated at Nside=8192 leads to a 1% bias in the binned boosted power spectrum at l~2000, while individual Cls exhibit ~5% fluctuations around the binned average. However, this bias is dominated by pixelization effects and not the aberration and Doppler shift of CMB photons that causes the fluctuations. Performing analysis on maps with galactic cuts does not induce any additional error in the boosted, binned power spectra over the full sky analysis. For multipoles that are free of resolution effects, there is no detectable deviation between the binned boosted and unboosted spectra. This result arises because the power spectrum is a slowly varying function of and does not show that, in general, Lorentz boosts can be neglected for other cosmological quantities such as polarization maps or higher-point functions.

Understanding the dynamic solar chromosphere is fundamental in solar physics. Spicules are an important feature of the chromosphere, connecting the photosphere to the corona, potentially mediating the transfer of energy and mass. The aim of this work is to study the properties of spicules over different regions of the sun. Our goal is to investigate if there is more than one type of spicules, and how spicules behave in the quiet sun, coronal holes, and active regions. We make use of high-cadence and high-spatial resolution Ca II H observations taken by Hinode/SOT. Making use of a semi-automated detection algorithm, we self-consistently track and measure the properties of 519 spicules over different regions. We find clear evidence of two types of spicules. Type I spicules show a rise and fall and have typical lifetimes of 150-400 s and maximum ascending velocities of 15-40 km/s, while type II spicules have shorter lifetimes of 50-150 s, faster velocities of 30-110 km/s, and are not seen to fall down, but rather fade at around their maximum length. Type II spicules are the most common, seen in quiet sun and coronal holes. Type I spicules are seen mostly in active regions. There are regional differences between quiet sun and coronal hole spicules, likely attributable to the different field configurations. The properties of type II spicules are consistent with published results of Rapid Blueshifted Events (RBEs), supporting the hypothesis that RBEs are their disk counterparts. For type I spicules we find the relations between their properties to be consistent with a magnetoacoustic shock wave driver, and with dynamic fibrils as their disk counterpart. The driver of type II spicules remains unclear from limb observations.

Radiative transfer modeling of spectral lines including partial redistribution (PRD) effects requires the evaluation of the ratio of the emission to the absorption profile. This quantity requires a large amount of computational work if one employs the angle-dependent redistribution function, which prohibits its use in 3D radiative transfer computations with model atmospheres containing velocity fields. We aim to provide a method to compute the emission to absorption profile ratio that requires less computational work but retains the effect of angle-dependent scattering in the resulting line profiles. We present a method to compute the profile ratio that employs the angle-averaged redistribution function and wavelength transforms to and from the rest frame of the scattering particles. We compare the emergent line profiles of the \MgII\,k and \Lyalpha\ lines computed with angle-dependent PRD, angle-averaged PRD and our new method in two representative test atmospheres. The new method yields a good approximation of true angle-dependent profile ratio and the resulting emergent line profiles while keeping the computational speed and simplicity of angle-averaged PRD theory.

In this work, the linear and gauge-invariant theory of cosmological perturbations in a class of anisotropic and shear-free spacetimes is developed. After constructing an explicit set of complete eigenfunctions in terms of which perturbations can be expanded, we identify the effective degrees of freedom during a generic slow-roll inflationary phase. These correspond to the anisotropic equivalent of the standard Mukhanov-Sasaki variables. The associated equations of motion present a remarkable resemblance to those found in perturbed Friedmann-Robertson-Walker spacetimes with curvature, apart from the spectrum of the Laplacian, which exhibits the characteristic frequencies of the underlying geometry. In particular, it is found that the perturbations cannot develop arbitrarily large super-Hubble modes.

Context: In studies of the solar spectrum relative to spectra of solar twin stars, it has been found that the chemical composition of the Sun seems to depart systematically from those of the twins. One possible explanation is that the effect is due to the special aspect angle of the Sun when observed from Earth, as compared with the aspect angles of the twins. Thus, a latitude dependence of the solar spectrum, even with the heliocentric angle constant, could lead to effects of the type observed. Aim: We explore a possible variation in the strength of certain spectral lines, used in the comparisons between the composition of the Sun and the twins, at loci on the solar disk with different latitudes but at constant heliocentric angle. Methods: We use the TRIPPEL spectrograph at the Swedish 1-m Solar Telescope on La Palma to record spectra in five spectral regions in order to compare different locations on the solar disk at a heliocentric angle of 45 deg. Equivalent widths and other parameters are measured for fifteen different lines representing nine atomic species. Results: The relative variations in equivalent widths at the equator and at solar latitude 45 deg are found to be less than 1.5 % for all spectral lines studied. Translated to elemental abundances as they would be measured from a terrestrial and a hypothetical pole-on observer, the difference is estimated to be within 0.005 dex in all cases. Conclusion: It is very unlikely that latitude effects could cause the reported abundance difference between the Sun and the solar twins. The accuracy obtainable in measurements of small differences in spectral line strengths between different solar disk positions is very high.

The largest fluctuation in the observed CMB temperature field is the dipole, its origin being usually attributed to the Doppler Effect - the Earth's velocity with respect to the CMB rest frame. The lowest order boost correction to temperature multipolar coefficients appears only as a second order correction in the temperature power spectrum, $C_{\ell}$. Since v/c - 10-3, this effect can be safely ignored when estimating cosmological parameters [4-7]. However, by cutting our galaxy from the CMB sky we induce large-angle anisotropies in the data. In this case, the corrections to the cut-sky $C_{\ell}$s show up already at first order in the boost parameter. In this paper we investigate this issue and argue that this effect might turn out to be important when reconstructing the power spectrum from the cut-sky data.

We study the evolution of inhomogeneous spherical perturbations in the universe in a way that generalizes the spherical top hat collapse in a straightforward manner. For that purpose we derive a dynamical equation for the evolution of the density contrast in the context of a Lemaitre-Tolman-Bondi metric and construct solutions with and without a cosmological constant for the evolution of a spherical perturbation with a given initial radial profile.

We review the basic hypotheses which motivate the statistical framework used to analyze the cosmic microwave background, and how that framework can be enlarged as we relax those hypotheses. In particular, we try to separate as much as possible the questions of gaussianity, homogeneity and isotropy from each other. We focus both on isotropic estimators of non-gaussianity as well as statistically anisotropic estimators of gaussianity, giving particular emphasis on their signatures and the enhanced "cosmic variances" that become increasingly important as our putative Universe becomes less symmetric. After reviewing the formalism behind some simple model-independent tests, we discuss how these tests can be applied to CMB data when searching for large scale "anomalies"

Gaussianity and statistical isotropy of the Universe are modern cosmology's minimal set of hypotheses. In this work we introduce a new statistical test to detect observational deviations from this minimal set. By defining the temperature correlation function over the whole celestial sphere, we are able to independently quantify both angular and planar dependence (modulations) of the CMB temperature power spectrum over different slices of this sphere. Given that planar dependence leads to further modulations of the usual angular power spectrum $C_l$, this test can potentially reveal richer structures in the morphology of the primordial temperature field. We have also constructed an unbiased estimator for this angular-planar power spectrum which naturally generalizes the estimator for the usual $C_l$'s. With the help of a chi-square analysis, we have used this estimator to search for observational deviations of statistical isotropy in WMAP's 5 year release data set (ILC5), where we found only slight anomalies on the angular scales $l=7$ and $l=8$. Since this angular-planar statistic is model-independent, it is ideal to employ in searches of statistical anisotropy (e.g., contaminations from the galactic plane) and to characterize non-Gaussianities.

We search for planar deviations of statistical isotropy in the Wilkinson Microwave Anisotropy Probe (WMAP) data by applying a recently introduced angular-planar statistics both to full-sky and to masked temperature maps, including in our analysis the effect of the residual foreground contamination and systematics in the foreground removing process as sources of error. We confirm earlier findings that full-sky maps exhibit anomalies at the planar ($l$) and angular ($\ell$) scales $(l,\ell)=(2,5),(4,7),$ and $(6,8)$, which seem to be due to unremoved foregrounds since this features are present in the full-sky map but not in the masked maps. On the other hand, our test detects slightly anomalous results at the scales $(l,\ell)=(10,8)$ and $(2,9)$ in the masked maps but not in the full-sky one, indicating that the foreground cleaning procedure (used to generate the full-sky map) could not only be creating false anomalies but also hiding existing ones. We also find a significant trace of an anomaly in the full-sky map at the scale $(l,\ell)=(10,5)$, which is still present when we consider galactic cuts of 18.3% and 28.4%. As regards the quadrupole ($\ell=2$), we find a coherent over-modulation over the whole celestial sphere, for all full-sky and cut-sky maps. Overall, our results seem to indicate that current CMB maps derived from WMAP data do not show significant signs of anisotropies, as measured by our angular-planar estimator. However, we have detected a curious coherence of planar modulations at angular scales of the order of the galaxy's plane, which may be an indication of residual contaminations in the full- and cut-sky maps.

We present results from a series of observational tests to 3D and 1D solar models. In particular, emphasis is given to the line formation of atomic oxygen lines, used to derive the much debated solar oxygen photospheric abundance. Using high-quality observations obtained with the Swedish Solar Telescope (SST) we study the centre-to-limb variation of the O I lines, testing the models and line formation (LTE and non-LTE). For the O I 777 nm triplet, the centre-to-limb variation sets strong constraints in the non-LTE line formation, and is used to derive an empirical correction factor (S_H) to the classical Drawin recipe for neutral hydrogen collisions. Taking advantage of the spatially-resolved character of the SST data, an additional framework for testing the 3D model and line formation is also studied. From the tests we confirm that the employed 3D model is realistic and its predictions agree very well with the observations.

Context: There is a lively debate about the solar oxygen abundance and the role of 3D models in its recent downward revision. The models have been tested using high resolution solar atlases. Further testing can be done using centre-to-limb variations. Aims: Using high quality observations of oxygen lines across the solar surface we seek to test if the 3D and 1D models reproduce their observed centre-to-limb variation (CLV). In particular we seek to assess whether the 3D model is appropriate to derive the solar oxygen abundance. Methods: We use our recent observations of OI 777 nm, OI 615.81 nm, [OI] 630.03 nm and nine lines of other elements for five viewing angles 0.2<mu<1 of the quiet solar disk. We compare them with the predicted line profiles from the 3D and 1D models computed using the most up-to-date line formation codes, line data and allowing for departures of LTE. The CLV of the OI 777 nm lines is also used to obtain an empirical correction for the poorly known efficiency of the inelastic collisions with H. Results: The 3D model generally reproduces the CLV observations of the lines very well, particularly the oxygen lines. From the OI 777 nm lines we find that the classical Drawin recipe slightly overestimates H collisions. The limb observations of the OI 615.82 nm line allow us to identify a previously unknown contribution of molecules for this line, prevalent at the solar limb. A detailed treatment of the [OI] 630.03 nm line shows that the 3D modeling provides an excellent agreement with the observations. The derived oxygen abundances with the 3D model are 8.68 (OI 777 nm), 8.66 ([OI] 630.03 nm) and 8.62 (OI 615.82 nm). Conclusions: These additional tests have reinforced the trustworthiness of the 3D model and line formation for abundance analyses.

Aims: we seek to provide additional tests of the line formation of theoretical 3D solar photosphere models. In particular, we set out to test the spatially-resolved line formation at several viewing angles, from the solar disk-centre to the limb and focusing on atomic oxygen lines. The purpose of these tests is to provide additional information on whether the 3D model is suitable to derive the solar oxygen abundance. We also aim to empirically constrain the NLTE recipes for neutral hydrogen collisions, using the spatially-resolved observations of the OI 777 nm lines. Methods: using the Swedish 1-m Solar Telescope we obtained high-spatial-resolution observations of five atomic oxygen lines (along with lines for other species) for five positions on the solar disk. These observations have a high spatial and spectral resolution, and a continuum intensity contrast up to 9% at 615 nm. The theoretical line profiles were computed using the 3D model, with a full 3D NLTE treatment for oxygen and LTE for the other lines. Results: at disk-centre we find an excellent agreement between predicted and observed line shifts, strengths, FWHM and asymmetries. At other viewing angles the agreement is also good, but the smaller continuum intensity contrast makes a quantitative comparison harder. We use the disk-centre observations we constrain S_H, the scaling factor for the efficiency of collisions with neutral hydrogen. We find that S_H=1 provides the best match to the observations. Conclusions: overall there is a very good agreement between predicted and observed line properties over the solar granulation. This further reinforces the view that the 3D model is realistic and a reliable tool to derive the solar oxygen abundance.

This article investigates the predictions of an inflationary phase starting from a homogeneous and anisotropic universe of the Bianchi I type. After discussing the evolution of the background spacetime, focusing on the number of e-folds and the isotropization, we solve the perturbation equations and predict the power spectra of the curvature perturbations and gravity waves at the end of inflation. The main features of the early anisotropic phase is (1) a dependence of the spectra on the direction of the modes, (2) a coupling between curvature perturbations and gravity waves, and (3) the fact that the two gravity waves polarisations do not share the same spectrum on large scales. All these effects are significant only on large scales and die out on small scales where isotropy is recovered. They depend on a characteristic scale that can, but a priori must not, be tuned to some observable scale. To fix the initial conditions, we propose a procedure that generalises the one standardly used in inflation but that takes into account the fact that the WKB regime is violated at early times when the shear dominates. We stress that there exist modes that do not satisfy the WKB condition during the shear-dominated regime and for which the amplitude at the end of inflation depends on unknown initial conditions. On such scales, inflation loses its predictability. This study paves the way to the determination of the cosmological signature of a primordial shear, whatever the Bianchi I spacetime. It thus stresses the importance of the WKB regime to draw inflationary predictions and demonstrates that when the number of e-folds is large enough, the predictions converge toward those of inflation in a Friedmann-Lemaitre spacetime but that they are less robust in the case of an inflationary era with a small number of e-folds.

This article describes the theory of cosmological perturbations around a homogeneous and anisotropic universe of the Bianchi I type. Starting from a general parameterisation of the perturbed spacetime a la Bardeen, a complete set of gauge invariant variables is constructed. Three physical degrees of freedom are identified and it is shown that, in the case where matter is described by a scalar field, they generalize the Mukhanov-Sasaki variables. In order to show that they are canonical variables, the action for the cosmological perturbations at second order is derived. Two major physical imprints of the primordial anisotropy are identified: (1) a scalar-tensor ``see-saw'' mechanism arising from the fact that scalar, vector and tensor modes do not decouple and (2) an explicit dependence of the statistical properties of the density perturbations and gravity waves on the wave-vector instead of its norm. This analysis extends, but also sheds some light on, the quantization procedure that was developed under the assumption of a Friedmann-Lemaitre background spacetime, and allows to investigate the robustness of the predictions of the standard inflationary scenario with respect to the hypothesis on the symmetries of the background spacetime. These effects of a primordial anisotropy may be related to some anomalies of the cosmic microwave background anisotropies on large angular scales.

We study the CMB constraints on two Dark Energy models described by scalar fields with different Lagrangians, namely a Klein-Gordon and a Born-Infeld field. The speed of sound of field fluctuations are different in these two theories, and therefore the predictions for CMB and structure formation are different. Employing the WMAP data on CMB, we make a likelihood analysis on a grid of theoretical models. We constrain the parameters of the models and compute the probability distribution functions for the equation of state. We show that the effect of the different sound speeds affects the low multipoles of CMB anisotropies, but is at most marginal for the class of models studied here.