In this work we analyse the structural and photometric properties of 21 barred simulated galaxies from the Auriga Project. These consist of Milky Way-mass magneto-hydrodynamical simulations in a $\Lambda$CDM cosmological context. In order to compare with observations, we generate synthetic SDSS-like broad-band images from the numerical data at z = 0 with different inclinations (from face-on to edge-on). Ellipse fits are used to determine the bar lengths, and 2D bulge/disc/bar decompositions with galfit are also performed, modelling the bar component with the modified Ferrer profile. We find a wide range of bar sizes and luminosities in the sample, and their structural parameters are in good agreement with the observations. All bulges present low S\'ersic indexes, and are classified as pseudobulges. In regard to the discs, the same breaks in the surface brightness profiles observed in real galaxies are found, and the radii at which these take place are in agreement with the observations. Also, from edge-on unsharp-masked images at z = 0, boxy or peanut-shaped (B/P) structures are clearly identified in the inner part of 4 bars, and also 2 more bars are found in buckling phase. The sizes of the B/P match fairly well with those obtained from observations. We thus conclude that the observed photometric and structural properties of galaxies with bars, which are the main drivers of secular evolution, can be developed in present state-of-the-art $\Lambda$CDM cosmological simulations.

Context. Most of our knowledge on the stellar component of galaxies is based on the analysis of distant systems and comes from integrated light data. It is important to test whether the results of the star formation histories (SFH) obtained with standard full-spectrum fitting methods are in agreement with those obtained through colour-magnitude diagram (CMD) fitting (usually considered the most reliable approach). Aims. We compare SFHs recovered from both techniques in Leo~A, a Local Group dwarf galaxy whose majority of stars formed during the last 8 Gyrs. This complements our previous findings in a field in the Large Magellanic Cloud bar, where star formation has been on-going since early epochs though at varying rates. Methods. We have used GTC/OSIRIS in long-slit mode to obtain a high-quality integrated light spectrum by scanning a selected region within Leo~A, for which a CMD reaching the old main mequence turn-off (oMSTO) is available from HST. We compared the SFH obtained from the two datasets, using state-of-art methods of integrated light ({\tt STECKMAP}) and resolved stellar population analysis. In the case of the CMD, we computed the SFH both from a deep CMD (observed with HST/ACS), and from a shallower one (archival data from HST/WFPC2). Results. The agreement between the SFHs recovered from the oMSTO CMD and from full spectrum fitting is remarkable, particularly regarding the time evolution of the star formation rate. The overall extremely low metallicity of Leo~A is recovered up to the last 2 Gyrs, when some discrepancies appear. A relatively high metallicity found for the youngest stars from the integrated data is a recurring feature that might indicate that the current models or synthesis codes should be revised, but that can be significantly mitigated using a more restrictive metallicity range... [Abridged]

We study a sample of 28 S0 galaxies extracted from the integral-field spectroscopic (IFS) survey CALIFA. We combine an accurate two-dimensional (2D) multi-component photometric decomposition with the IFS kinematic properties of their bulges to understand their formation scenario. Our final sample is representative of S0s with high stellar masses ($M_{star}/M_{sun} > 10^{10}$). They lay mainly on the red sequence and live in relatively isolated environments similar to that of the field and loose groups. We use our 2D photometric decomposition to define the size and photometric properties of the bulges, as well as their location within the galaxies. We perform mock spectroscopic simulations mimicking our observed galaxies to quantify the impact of the underlying disc on our bulge kinematic measurements ($\lambda$ and $v/\sigma$). We compare our bulge corrected kinematic measurements with the results from Schwarzschild dynamical modelling. The good agreement confirms the robustness of our results and allows us to use bulge reprojected values of $\lambda$ and $v/\sigma$. We find that the photometric ($n$ and $B/T$) and kinematic ($v/\sigma$ and $\lambda$) properties of our field S0 bulges are not correlated. We demonstrate that this morpho-kinematic decoupling is intrinsic to the bulges and it is not due to projection effects. We conclude that photometric diagnostics to separate different types of bulges (disc-like vs classical) might not be useful for S0 galaxies. The morpho-kinematics properties of S0 bulges derived in this paper suggest that they are mainly formed by dissipation processes happening at high redshift, but dedicated high-resolution simulations are necessary to better identify their origin.

Aims: The physics driving features such as breaks observed in galaxy surface brightness (SB) profiles remains contentious. Here, we assess the importance of stellar radial motions in shaping their characteristics. Methods: We use the simulated Milky Way-mass, cosmological discs, from the Ramses Disc Environment Study (RaDES) to characterise the radial redistribution of stars in galaxies displaying type I (pure exponentials), II (downbending), and III (upbending) SB profiles. We compare radial profiles of the mass fractions and the velocity dispersions of different sub-populations of stars according to their birth and current locations. Results: Radial redistribution of stars is important in all galaxies regardless of their light profiles. Type II breaks seem to be a consequence of the combined effects of outward-moving and accreted stars. The former produces shallower inner profiles (lack of stars in the inner disc) and accumulate material around the break radius and beyond, strengthening the break; the latter can weaken or even convert the break into a pure exponential. Further accretion from satellites can concentrate material in the outermost parts, leading to type III breaks that can coexist with type II breaks, but situated further out. Type III galaxies would be the result of an important radial redistribution of material throughout the entire disc, as well as a concentration of accreted material in the outskirts. In addition, type III galaxies display the most efficient radial redistribution and the largest number of accreted stars, followed by type I and II systems, suggesting that type I galaxies may be an intermediate case between types II and III. In general, the velocity dispersion profiles of all galaxies tend to flatten or even ncrease around the locations where the breaks are found. The age and metallicity profiles are also affected, exhibiting...[abridged]

We estimate the current extinction-corrected H$\alpha$ star formation rate (SFR) of the different morphological components that shape galaxies (bulges, bars, and disks). We use a multi-component photometric decomposition based on SDSS imaging to CALIFA Integral Field Spectroscopy datacubes for a sample of 219 galaxies. This analysis reveals an enhancement of the central SFR and specific SFR (sSFR $=$ SFR/$M_{\star}$) in barred galaxies. Along the Main Sequence, we find more massive galaxies in total have undergone efficient suppression (quenching) of their star formation, in agreement with many studies. We discover that more massive disks have had their star formation quenched as well. We evaluate which mechanisms might be responsible for this quenching process. The presence of type-2 AGNs plays a role at damping the sSFR in bulges and less efficiently in disks. Also, the decrease in the sSFR of the disk component becomes more noticeable for stellar masses around 10$^{10.5}$ M$_{\odot}$; for bulges, it is already present at $\sim$10$^{9.5}$ M$_{\odot}$. The analysis of the line-of-sight stellar velocity dispersions ($\sigma$) for the bulge component and of the corresponding Faber-Jackson relation shows that AGNs tend to have slightly higher $\sigma$ values than star-forming galaxies for the same mass. Finally, the impact of environment is evaluated by means of the projected galaxy density, $\Sigma$$_{5}$. We find that the SFR of both bulges and disks decreases in intermediate-to-high density environments. This work reflects the potential of combining IFS data with 2D multi-component decompositions to shed light on the processes that regulate the SFR.

We present a galaxy classification system for 238 (E1-Sdm) CALIFA (Calar Alto Legacy Integral Field Area) galaxies based on the shapes and amplitudes of their circular velocity curves (CVCs). We infer the CVCs from the de-projected surface brightness of the galaxies, after scaling by a constant mass-to-light ratio based on stellar dynamics - solving axisymmetric Jeans equations via fitting the second velocity moment $V_{\mathrm{rms}}=\sqrt{V^2+\sigma^2}$ of the stellar kinematics. We use principal component analysis (PCA) applied to the CVC shapes to find characteristic features and use a $k$-means classifier to separate circular curves into classes. This objective classification method identifies four different classes, which we name slow-rising (SR), flat (FL), round-peaked (RP) and sharp-peaked (SP) circular curves. SR are typical for low-mass, late-type (Sb-Sdm), young, faint, metal-poor and disc-dominated galaxies. SP are typical for high-mass, early-type (E1-E7), old, bright, metal-rich and bulge-dominated galaxies. FL and RP appear presented by galaxies with intermediate mass, age, luminosity, metallicity, bulge-to-disk ratio and morphologies (E4-S0a, Sa-Sbc). The discrepancy mass factor, $f_d=1-M_{*}/M_{dyn}$, have the largest value for SR and SP classes ($\sim$ 74 per cent and $\sim$ 71 per cent, respectively) in contrast to the FL and RP classes (with $\sim$ 59 per cent and $\sim$ 61 per cent, respectively). Circular curve classification presents an alternative to typical morphological classification and appears more tightly linked to galaxy evolution.

Boxy/peanut bulges are considered to be part of the same stellar structure as bars and both could be linked through the buckling instability. The Milky Way is our closest example. The goal of this letter is determining if the mass assembly of the different components leaves an imprint in their stellar populations allowing to estimate the time of bar formation and its evolution. To this aim we use integral field spectroscopy to derive the stellar age distributions, SADs, along the bar and disc of NGC 6032. The analysis shows clearly different SADs for the different bar areas. There is an underlying old (>=12 Gyr) stellar population for the whole galaxy. The bulge shows star formation happening at all times. The inner bar structure shows stars of ages older than 6 Gyrs with a deficit of younger populations. The outer bar region presents a SAD similar to that of the disc. To interpret our results, we use a generic numerical simulation of a barred galaxy. Thus, we constrain, for the first time, the epoch of bar formation, the buckling instability period and the posterior growth from disc material. We establish that the bar of NGC 6032 is old, formed around 10 Gyr ago while the buckling phase possibly happened around 8 Gyr ago. All these results point towards bars being long-lasting even in the presence of gas.

Spiral arms are the most singular features in disc galaxies. These structures can exhibit different patterns, namely grand design and flocculent arms, with easily distinguishable characteristics. However, their origin and the mechanisms shaping them are unclear. The overall role of spirals in the chemical evolution of disc galaxies is another unsolved question. In particular, it has not been fully explored if the \hii\,regions of spiral arms present different properties from those located in the interarm regions. Here we analyse the radial oxygen abundance gradient of the arm and interarm star forming regions of 63 face-on spiral galaxies using CALIFA Integral Field Spectroscopy data. We focus the analysis on three characteristic parameters of the profile: slope, zero-point, and scatter. The sample is morphologically separated into flocculent versus grand design spirals and barred versus unbarred galaxies. We find subtle but statistically significant differences between the arm and interarm distributions for flocculent galaxies, suggesting that the mechanisms generating the spiral structure in these galaxies may be different to those producing grand design systems, for which no significant differences are found. We also find small differences in barred galaxies, not observed in unbarred systems, hinting that bars may affect the chemical distribution of these galaxies but not strongly enough as to be reflected in the overall abundance distribution. In light of these results, we propose bars and flocculent structure as two distinct mechanisms inducing differences in the abundance distribution between arm and interarm star forming regions.

Context: According to numerical simulations, stars are not always kept at their birth galactocentric distances but migrate. The importance of this radial migration in shaping galactic light distributions is still unclear. However, if it is indeed important, galaxies with different surface brightness (SB) profiles must display differences in their stellar population properties. Aims: We investigate the role of radial migration on the light distribution and the radial stellar content by comparing the inner colour, age and metallicity gradients for galaxies with different SB profiles. We define these inner parts avoiding the bulge and bar regions and up to around three disc scale-lengths (type I, pure exponential) or the break radius (type II, downbending; type III, upbending). Methods: We analyse 214 spiral galaxies from the CALIFA survey covering different SB profiles. We make use of GASP2D and SDSS data to characterise their light distribution and obtain colour profiles. The stellar age and metallicity profiles are computed using a methodology based on full-spectrum fitting techniques (pPXF, GANDALF, and STECKMAP) to the IFS CALIFA data. Results: The distributions of the colour, stellar age and stellar metallicity gradients in the inner parts for galaxies displaying different SB profiles are unalike as suggested by Kolmogorov-Smirnov and Anderson-Darling tests. We find a trend in which type II galaxies show the steepest profiles of all and type III the shallowest, with type I galaxies displaying an intermediate behaviour. Conclusions: These results are consistent with a scenario in which radial migration is more efficient for type III galaxies than for type I systems with type II galaxies presenting the lowest radial migration efficiency. In such scenario, radial migration mixes the stellar content flattening the radial stellar properties and... [abriged]

We present a two-dimensional multi-component photometric decomposition of 404 galaxies from the CALIFA Data Release 3. They represent all possible galaxies with no clear signs of interaction and not strongly inclined in the final CALIFA data release. Galaxies are modelled in the g, r, and i SDSS images including, when appropriate, a nuclear point source, bulge, bar, and an exponential or broken disc component. We use a human-supervised approach to determine the optimal number of structures to be included in the fit. The dataset, including the photometric parameters of the CALIFA sample, is released together with statistical errors and a visual analysis of the quality of each fit. The analysis of the photometric components reveals a clear segregation of the structural composition of galaxies with stellar mass. At high masses (log(Mstar/Msun)>11), the galaxy population is dominated by galaxies modelled with a single Sersic or a bulge+disc with a bulge-to-total (B/T) luminosity ratio B/T>0.2. At intermediate masses (9.5<log(Mstar/Msun)<11), galaxies described with bulge+disc but B/T < 0.2 are preponderant, whereas, at the low mass end (log(Mstar/Msun)<9.5), the prevailing population is constituted by galaxies modelled with either pure discs or nuclear point sources+discs (i.e., no discernible bulge). We obtain that 57% of the volume corrected sample of disc galaxies in the CALIFA sample host a bar. This bar fraction shows a significant drop with increasing galaxy mass in the range 9.5<log(Mstar/Msun)<11.5. The analyses of the extended multi-component radial profile result in a volume-corrected distribution of 62%, 28%, and 10% for the so-called Type I, Type II, and Type III disc profiles, respectively. These fractions are in discordance with previous findings. We argue that the different methodologies used to detect the breaks are the main cause for these differences.

We compute and investigate four types of imprint of a stochastic background of primordial magnetic fields (PMFs) on the cosmic microwave background (CMB) anisotropies: the impact of PMFs on the CMB spectra; the effect on CMB polarization induced by Faraday rotation; the impact of PMFs on the ionization history; magnetically-induced non-Gaussianities; and the magnetically-induced breaking of statistical isotropy. Overall, Planck data constrain the amplitude of PMFs to less than a few nanogauss. In particular, individual limits coming from the analysis of the CMB angular power spectra, using the Planck likelihood, are $B_{1\,\mathrm{Mpc}}< 4.4$ nG (where $B_{1\,\mathrm{Mpc}}$ is the comoving field amplitude at a scale of 1 Mpc) at 95% confidence level, assuming zero helicity, and $B_{1\,\mathrm{Mpc}}< 5.6$ nG for a maximally helical field.For nearly scale-invariant PMFs we obtain $B_{1\,\mathrm{Mpc}}<2.0$ nG and $B_{1\,\mathrm{Mpc}}<0.9$ nG if the impact of PMFs on the ionization history of the Universe is included. From the analysis of magnetically-induced non-Gaussianity we obtain three different values, corresponding to three applied methods, all below 5 nG. The constraint from the magnetically-induced passive-tensor bispectrum is $B_{1\,\mathrm{Mpc}}< 2.8$ nG. A search for preferred directions in the magnetically-induced passive bispectrum yields $B_{1\,\mathrm{Mpc}}< 4.5$ nG, whereas the the compensated-scalar bispectrum gives $B_{1\,\mathrm{Mpc}}< 3$ nG. The analysis of the Faraday rotation of CMB polarization by PMFs uses the Planck power spectra in $EE$ and $BB$ at 70 GHz and gives $B_{1\,\mathrm{Mpc}}< 1380$ nG. In our final analysis, we consider the harmonic-space correlations produced by Alfv\'en waves, finding no significant evidence for the presence of these waves. Together, these results comprise a comprehensive set of constraints on possible PMFs with Planck data.

We measured the gas abundance profiles in a sample of 122 face-on spiral galaxies observed by the CALIFA survey and included all spaxels whose line emission was consistent with star formation. This type of analysis allowed us to improve the statistics with respect to previous studies, and to properly estimate the oxygen distribution across the entire disc to a distance of up to 3-4 disc effective radii (r$_e$). We confirm the results obtained from classical HII region analysis. In addition to the general negative gradient, an outer flattening can be observed in the oxygen abundance radial profile. An inner drop is also found in some cases. There is a common abundance gradient between 0.5 and 2.0 r$_e$ of $\alpha_{O/H} = -\,0.075\,\rm{dex}/r_e$ with a scatter of $\sigma = 0.016\,\rm{dex}/r_e$ when normalising the distances to the disc effective radius. By performing a set of Kolmogorov-Smirnov tests, we determined that this slope is independent of other galaxy properties, such as morphology, absolute magnitude, and the presence or absence of bars. In particular, barred galaxies do not seem to display shallower gradients, as predicted by numerical simulations. Interestingly, we find that most of the galaxies in the sample with reliable oxygen abundance values beyond $\sim 2$ effective radii (57 galaxies) present a flattening of the abundance gradient in these outer regions. This flattening is not associated with any morphological feature, which suggests that it is a common property of disc galaxies. Finally, we detect a drop or truncation of the abundance in the inner regions of 27 galaxies in the sample; this is only visible for the most massive galaxies.

Aims: We study the effects of the cosmological assembly history on the chemical and dynamical properties of the discs of spiral galaxies as a function of radius. Methods: We make use of the simulated Milky-Way mass, fully-cosmological discs, from {\tt RaDES} (Ramses Disc Environment Study). We analyse their assembly history by examining the proximity of satellites to the galactic disc, instead of their merger trees, to better gauge which satellites impact the disc. We present stellar age and metallicity profiles, Age-Metallicity Relation (AMR), Age-Velocity dispersion Relation (AVR), and Stellar Age Distribution (SAD) in several radial bins for the simulated galaxies. Results: Assembly histories can be divided into three different stages: i) a merger dominated phase, when a large number of mergers with mass ratios of $\sim$1:1 take place (lasting $\sim$3.2$\pm$0.4 Gyr on average); ii) a quieter phase, when $\sim$1:10 mergers take place (lasting $\sim$4.4$\pm$2.0 Gyr) - these two phases are able to kinematically heat the disc and produce a disc that is chemically mixed over its entire radial extension; iii) a "secular" phase where the few mergers that take place have mass ratios below 1:100, and which do not affect the disc properties (lasting $\sim$5.5$\pm$2.0 Gyr). Phase ii ends with a final merger event (at time $t_\mathrm{jump}$) marking the onset of important radial differences in the AMR, AVR, and SAD. Conclusions: Inverted AMR trends in the outer parts of discs, for stars younger than $t_\mathrm{jump}$, are found as the combined effect of radial motions and star formation in satellites temporarily located in these outer parts. "U-shaped" stellar age profiles change to an old plateau ($\sim$10 Gyr) in the outer discs for the entire {\tt RaDES} sample. This shape is a consequence of inside-out growth of the disc, radial motions of disc stars ... [abridged]

We study the stellar properties of 44 face-on spiral galaxies from the Calar Alto Legacy Integral Field Area survey via full spectrum fitting techniques. We compare the age profiles with the surface brightness distribution in order to highlight differences between profile types (type I, exponential profile; and II, down-bending profile). We observe an upturn ("U-shape") in the age profiles for 17 out of these 44 galaxies with reliable stellar information up to their outer parts. This "U-shape" is not a unique feature for type II galaxies but can be observed in type I as well. These findings suggest that the mechanisms shaping the surface brightness and stellar population distributions are not directly coupled. This upturn in age is only observable in the light-weighted profiles while it flattens out in the mass-weighted profiles. Given recent results on the outer parts of nearby systems and the results presented in this Letter, one of the most plausible explanations for the age upturn is an early formation of the entire disc ($\sim$~10~Gyr ago) followed by an inside-out quenching of the star formation.

Bar-induced gas inflows towards the galaxy centres are recognized as a key agent for the secular evolution of galaxies. One immediate consequence is the accumulation of gas in the centre of galaxies where it can form stars and alter the chemical and physical properties. We use a sample of nearby face--on disc galaxies with available SDSS spectra to study whether the properties of the ionised gas in the central parts (radii <~0.6-2.1 kpc) of barred galaxies are altered by the presence of a bar, and whether the bar effect is related to bar and/or parent galaxy properties. The distributions of all parameters analysed are different for barred and unbarred galaxies, except for the R23 metallicity tracer and the oxygen abundance (from photoionisation models). The median values point towards (marginally) larger dust content, star formation rate per unit area, electron density and ionisation parameter in the centres of barred galaxies than in the unbarred counterpart. The most remarkable barred/unbarred difference appears in the [NII]6583/Ha line ratio, which is on average ~25% larger in barred galaxies, due to a larger N/O in the centres of these galaxies. We observe an enhancement of the central gas differences in later-type galaxies or galaxies with less massive bulges. However the bar seems to have a lower impact on the central gas properties for galaxies with more massive bulges (M_bulge > 10^10 M_sun) or galaxies with total stellar mass above ~ 10^10.8 M_sun. In conclusion, we find observational evidence that the presence of a galactic bar affects the central ionised gas properties of disc galaxies, where the most striking effect is an enhancement in the N/O abundance ratio, which can be qualitatively interpreted as due to a different origin or evolutionary processes for less and more massive bulges, with the gaseous phase of the former having currently a closer relation with bars.

We study, for the first time in a statistically significant and well-defined sample, the relation between the outer-disk ionized-gas metallicity gradients and the presence of breaks in the surface brightness profiles of disk galaxies. SDSS g'- and r'-band surface brightness, (g'- r') color, and ionized-gas oxygen abundance profiles for 324 galaxies within the CALIFA survey are used for this purpose. We perform a detailed light-profile classification finding that 84% of our disks show down- or up-bending profiles (Type II and Type III, respectively) while the remaining 16% are well fitted by one single exponential (Type I). The analysis of the color gradients at both sides of this break shows a U-shaped profile for most Type II galaxies with an average minimum (g'- r') color of ~0.5 mag and a ionized-gas metallicity flattening associated to it only in the case of low-mass galaxies. More massive systems show a rather uniform negative metallicity gradient. The correlation between metallicity flattening and stellar mass results in p-values as low as 0.01. Independently of the mechanism having shaped the outer light profiles of these galaxies, stellar migration or a previous episode of star formation in a shrinking star-forming disk, it is clear that the imprint in their ionized-gas metallicity was different for low- and high-mass Type II galaxies. In the case of Type III disks, a positive correlation between the change in color and abundance gradient is found (the null hypothesis is ruled out with a p-value of 0.02), with the outer disks of Type III galaxies with masses $\leq$10$^{10}$ M$_{\odot}$ showing a weak color reddening or even a bluing. This is interpreted as primarily due to a mass down-sizing effect on the population of Type III galaxies having recently experienced an enhanced inside-out growth.

Accurate star formation histories (SFHs) of galaxies are fundamental for understanding the build-up of their stellar content. However, the most accurate SFHs - those obtained from colour-magnitude diagrams (CMDs) of resolved stars reaching the oldest main sequence turnoffs (oMSTO) - are presently limited to a few systems in the Local Group. It is therefore crucial to determine the reliability and range of applicability of SFHs derived from integrated light spectroscopy, as this affects our understanding of unresolved galaxies from low to high redshift. To evaluate the reliability of current full spectral fitting techniques in deriving SFHs from integrated light spectroscopy by comparing SFHs from integrated spectra to those obtained from deep CMDs of resolved stars. We have obtained a high signal--to--noise (S/N $\sim$ 36.3 per \AA) integrated spectrum of a field in the bar of the Large Magellanic Cloud (LMC) using EFOSC2 at the 3.6 meter telescope at La Silla Observatory. For this same field, resolved stellar data reaching the oMSTO are available. We have compared the star formation rate (SFR) as a function of time and the age-metallicity relation (AMR) obtained from the integrated spectrum using {\tt STECKMAP}, and the CMD using the IAC-star/MinnIAC/IAC-pop set of routines. For the sake of completeness we also use and discuss other synthesis codes ({\tt STARLIGHT} and {\tt ULySS}) to derive the SFR and AMR from the integrated LMC spectrum. We find very good agreement (average differences $\sim$ 4.1 $\%$) between the SFR(t) and the AMR obtained using {\tt STECKMAP} on the integrated light spectrum, and the CMD analysis. {\tt STECKMAP} minimizes the impact of the age-metallicity degeneracy and has the advantage of preferring smooth solutions to recover complex SFHs by means of a penalized $\chi^2$. [abridged]

The European Space Agency's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14~May 2009 and scanned the microwave and submillimetre sky continuously between 12~August 2009 and 23~October 2013. In February~2015, ESA and the Planck Collaboration released the second set of cosmology products based on data from the entire Planck mission, including both temperature and polarization, along with a set of scientific and technical papers and a web-based explanatory supplement. This paper gives an overview of the main characteristics of the data and the data products in the release, as well as the associated cosmological and astrophysical science results and papers. The science products include maps of the cosmic microwave background (CMB), the thermal Sunyaev-Zeldovich effect, and diffuse foregrounds in temperature and polarization, catalogues of compact Galactic and extragalactic sources (including separate catalogues of Sunyaev-Zeldovich clusters and Galactic cold clumps), and extensive simulations of signals and noise used in assessing the performance of the analysis methods and assessment of uncertainties. The likelihood code used to assess cosmological models against the Planck data are described, as well as a CMB lensing likelihood. Scientific results include cosmological parameters deriving from CMB power spectra, gravitational lensing, and cluster counts, as well as constraints on inflation, non-Gaussianity, primordial magnetic fields, dark energy, and modified gravity.

The Square Kilometre Array (SKA) is called to revolutionise essentially all areas of Astrophysics. With a collecting area of about a square kilometre, the SKA will be a transformational instrument, and its scientific potential will go beyond the interests of astronomers. Its technological challenges and huge cost requires a multinational effort, and Europe has recognised this by putting the SKA on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI). The Spanish SKA White Book is the result of the coordinated effort of 120 astronomers from 40 different research centers. The book shows the enormous scientific interest of the Spanish astronomical community in the SKA and warrants an optimum scientific exploitation of the SKA by Spanish researchers, if Spain enters the SKA project.

This paper describes the Second Public Data Release (DR2) of the Calar Alto Legacy Integral Field Area (CALIFA) survey. The data for 200 objects are made public, including the 100 galaxies of the First Public Data Release (DR1). Data were obtained with the integral-field spectrograph PMAS/PPak mounted on the 3.5 m telescope at the Calar Alto observatory. Two different spectral setups are available for each galaxy, (i) a low-resolution V500 setup covering the wavelength range 3745-7500 \AA\ with a spectral resolution of 6.0 \AA\ (FWHM), and (ii) a medium-resolution V1200 setup covering the wavelength range 3650-4840 \AA\ with a spectral resolution of 2.3 \AA\ (FWHM). The sample covers a redshift range between 0.005 and 0.03, with a wide range of properties in the Color-Magnitude diagram, stellar mass, ionization conditions, and morphological types. All released cubes were reduced with the latest pipeline, including improved spectrophotometric calibration, spatial registration and spatial resolution. The spectrophotometric calibration is better than 6% and the median spatial resolution is 2.5". Altogether the second data release contains over 1.5 million spectra. It is available at http://califa.caha.es/DR2.

We present optical integral field spectroscopy (IFS) observations of the Mice, a major merger between two massive (>10^11Msol) gas-rich spirals NGC4676A and B, observed between first passage and final coalescence. The spectra provide stellar and gas kinematics, ionised gas properties and stellar population diagnostics, over the full optical extent of both galaxies. The Mice provide a perfect case study highlighting the importance of IFS data for improving our understanding of local galaxies. The impact of first passage on the kinematics of the stars and gas has been significant, with strong bars likely induced in both galaxies. The barred spiral NGC4676B exhibits a strong twist in both its stellar and ionised gas disk. On the other hand, the impact of the merger on the stellar populations has been minimal thus far: star formation induced by the recent close passage has not contributed significantly to the global star formation rate or stellar mass of the galaxies. Both galaxies show bicones of high ionisation gas extending along their minor axes. In NGC4676A the high gas velocity dispersion and Seyfert-like line ratios at large scaleheight indicate a powerful outflow. Fast shocks extend to ~6.6kpc above the disk plane. The measured ram pressure and mass outflow rate (~8-20Msol/yr) are similar to superwinds from local ULIRGs, although NGC4676A has only a moderate infrared luminosity of 3x10^10Lsol. Energy beyond that provided by the mechanical energy of the starburst appears to be required to drive the outflow. We compare the observations to mock kinematic and stellar population maps from a merger simulation. The models show little enhancement in star formation during and following first passage, in agreement with the observations. We highlight areas where IFS data could help further constrain the models.

We present the results on the study of the global and local M-Z relation based on the first data available from the CALIFA survey (150 galaxies). This survey provides integral field spectroscopy of the complete optical extent of each galaxy (up to 2-3 effective radii), with enough resolution to separate individual HII regions and/or aggregations. Nearly $\sim$3000 individual HII regions have been detected. The spectra cover the wavelength range between [OII]3727 and [SII]6731, with a sufficient signal-to-noise to derive the oxygen abundance and star-formation rate associated with each region. In addition, we have computed the integrated and spatially resolved stellar masses (and surface densities), based on SDSS photometric data. We explore the relations between the stellar mass, oxygen abundance and star-formation rate using this dataset. We derive a tight relation between the integrated stellar mass and the gas-phase abundance, with a dispersion smaller than the one already reported in the literature ($\sigma_{\Delta{\rm log(O/H)}}=$0.07 dex). Indeed, this dispersion is only slightly larger than the typical error derived for our oxygen abundances. However, we do not find any secondary relation with the star-formation rate, other than the one induced due to the primary relation of this quantity with the stellar mass. We confirm the result using the $\sim$3000 individual HII regions, for the corresponding local relations. Our results agree with the scenario in which gas recycling in galaxies, both locally and globally, is much faster than other typical timescales, like that of gas accretion by inflow and/or metal loss due to outflows. In essence, late-type/disk dominated galaxies seem to be in a quasi-steady situation, with a behavior similar to the one expected from an instantaneous recycling/closed-box model.

This work investigates the effect of the aperture size on derived galaxy properties for which we have spatially-resolved optical spectra. We focus on some indicators of star formation activity and dust attenuation for spiral galaxies that have been widely used in previous work on galaxy evolution. We have used 104 spiral galaxies from the CALIFA survey for which 2D spectroscopy with complete spatial coverage is available. From the 3D cubes we have derived growth curves of the most conspicuous Balmer emission lines (Halpha, Hbeta) for circular apertures of different radii centered at the galaxy's nucleus after removing the underlying stellar continuum. We find that the Halpha flux (f(Halpha)) growth curve follows a well defined sequence with aperture radius showing low dispersion around the median value. From this analysis, we derive aperture corrections for galaxies in different magnitude and redshift intervals. Once stellar absorption is properly accounted for, the f(Halpha)/f(Hbeta) ratio growth curve shows a smooth decline, pointing towards the absence of differential dust attenuation as a function of radius. Aperture corrections as a function of the radius are provided in the interval [0.3,2.5]R_50. Finally, the Halpha equivalent width (EW(Halpha)) growth curve increases with the size of the aperture and shows a very large dispersion for small apertures. This large dispersion prevents the use of reliable aperture corrections for this quantity. In addition, this result suggests that separating star-forming and quiescent galaxies based on observed EW(Halpha) through small apertures is likely to result in low EW(Halpha) star-forming galaxies begin classified as quiescent.

We present the first public data release of the CALIFA survey. It consists of science-grade optical datacubes for the first 100 of eventually 600 nearby (0.005<z<0.03) galaxies, obtained with the integral-field spectrograph PMAS/PPak mounted on the 3.5m telescope at the Calar Alto observatory. The galaxies in DR1 already cover a wide range of properties in color-magnitude space, morphological type, stellar mass, and gas ionization conditions. This offers the potential to tackle a variety of open questions in galaxy evolution using spatially resolved spectroscopy. Two different spectral setups are available for each galaxy, (i) a low-resolution V500 setup covering the nominal wavelength range 3745-7500A with a spectral resolution of 6.0A (FWHM), and (ii) a medium-resolution V1200 setup covering the nominal wavelength range 3650-4840A with a spectral resolution of 2.3A (FWHM). We present the characteristics and data structure of the CALIFA datasets that should be taken into account for scientific exploitation of the data, in particular the effects of vignetting, bad pixels and spatially correlated noise. The data quality test for all 100 galaxies showed that we reach a median limiting continuum sensitivity of 1.0x10^-18erg/s/cm^2/A/arcsec^2 at 5635A and 2.2x10^-18erg/s/cm^2/A/arcsec^2 at 4500A for the V500 and V1200 setup respectively, which corresponds to limiting r and g band surface brightnesses of 23.6mag/arcsec^2 and 23.4mag/arcsec^2, or an unresolved emission-line flux detection limit of roughly 1x10^-17erg/s/cm^2/arcsec^2 and 0.6x10^-17erg/s/cm^2/arcsec^2, respectively. The median spatial resolution is 3.7", and the absolute spectrophotometric calibration is better than 15% (1sigma). We also describe the available interfaces and tools that allow easy access to this first public CALIFA data.

This is the third paper of a series devoted to study the properties of bars from long slit spectroscopy to understand their formation, evolution and their influence on the evolution of disk galaxies. In this work we aim to determine the gas metallicity distribution of a sample of 20 barred early-type galaxies. We compare the nebular and stellar metallicity distributions to conclude about the origin of the warm gas. We compare the results of nebular emission metallicities using different semi-empirical methods. We carry out AGN diagnostic diagrams along the radius to determine the radius of influence of the AGN and the nuclei nature of the studied galaxies. We then derive the gas metallicities along the bars and compare the results to the distribution of stellar metallicities in the same regions. Most of the gas emission is centrally concentrated, although 15 galaxies also show emission along the bar. In the central regions, gas oxygen abundances are in the range 12+$\log$(O/H)= 8.4-9.1. The nebular metallicity gradients are very shallow in the bulge and bar regions. For three galaxies (one of them a LINER), the gas metallicities lie well below the stellar ones in the bulge region. These results do not depend on the choice of the semi-empirical calibration used to calculate the abundances. We see that the galaxies with the lowest abundances are those with the largest rotational velocities. The presence of gas of significantly lower metallicity than the stellar abundances in three of our galaxies, points to an external origin as the source of the gas that fuels the present star formation in the centre of some early-type barred galaxies. The fact that the bar/disk nebular metallicities are higher than the central ones might be indicating that the gas could be accreted via cooling flows instead of radial accretion from gas sitting in the outer parts of the disk.

Recent observations of the rotation curve of M31 show a rise of the outer part that can not be understood in terms of standard dark matter models or perturbations of the galactic disc by M31's satellites. Here, we propose an explanation of this dynamical feature based on the influence of the magnetic field within the thin disc. We have considered standard mass models for the luminous mass distribution, a NFW model to describe the dark halo, and we have added up the contribution to the rotation curve of a magnetic field in the disc, which is described by an axisymmetric pattern. Our conclusion is that a significant improvement of the fit in the outer part is obtained when magnetic effects are considered. The best-fit solution requires an amplitude of ~4 microG with a weak radial dependence between 10 and 38 kpc.

Context: Warps occurring in galactic discs have been studied extensively in HI and in the optical, but rarely in the near-infrared (NIR) bands that trace the older stellar populations. Aims: We provide NIR data of nearby edge-on galaxies, combined with optical observations, for direct comparison of the properties of galactic warps as a function of wavelength, and calculate warp curves for each galaxy and obtain the characteristic warp parameters. We discuss these properties as possible constraints to the different mechanisms that have been proposed for the development and persistence of galactic warps. Methods: We observed 20 galaxies that were selected from a statistically complete diameter-limited subsample of edge-on disc galaxies. We used the Cerro Tololo Infrared Imager (CIRIM) at the CTIO 1.5m Ritchey-Chretien telescope to acquire the NIR data. We used the 1.54m Danish and 0.92m Dutch telescopes at the European Southern Observatory's La Silla site for our optical observations. Results: Our results show that 13 of our 20 sample galaxies are warped, with the warp more pronounced in the optical than at NIR wavelengths. In the remaining seven galaxies, no warp is apparent within the limitations of our automated detection method. The transition between the unperturbed inner disc and the outer, warped region is rather abrupt. S0 galaxies exhibit very small or no warps. The magnetic model remains one of a number of interesting formation scenarios.

CONTEXT. Hydrodynamical cosmological simulations predict flows of the intergalactic medium along the radial vector of the voids, approximately in the direction of the infall of matter at the early stages of the galaxy formation. AIMS. These flows might be detected by analysing the dependence of the warp amplitude on the inclination of the galaxies at the shells of the voids with respect to the radial vector of the voids. This analysis will be the topic of this paper. METHODS. We develop a statistical method of analysing the correlation of the amplitude of the warp and the inclination of the galaxy at the void surface. This is applied to a sample of 97 edge-on galaxies from the Sloan Digital Sky Survey. Our results are compared with the theoretical expectations, which are also derived in this paper. RESULTS. Our results allow us to reject the null hypothesis (i.e., the non-correlation of the warp amplitude and the inclination of the galaxy with respect to the void surface) at 94.4% C. L., which is not conclusive. The absence of the radial flows cannot be excluded at present, although we can put a constraint on the maximum average density of baryonic matter of the radial flows of <rho_b> <~ 4 Omega_b rho_crit.

The results obtained in paper I are used to study possible relationships between the truncation radius of stellar discs in the NIR and structural parameters of the galaxies. The NIR truncation radius is larger for brighter galaxies, being proportional to $V_m^c$ with $c \approx 3/2$, and with $V_m$ being the asymptotic rotation velocity at large radii (when the rotation curve becomes flat), and is lower for higher wavelengths. When it is normalized to the scalelength, the truncation is an increasing function of the central surface brightness and is lower for late type galaxies, although these correlations are weaker. These relations are in agreement with the scenario of magnetically driven truncations.

We present NIR surface photometry of 11 edge-on galaxies obtained in the course of a long term project aimed at analysing the occurrence and type of the truncation of the outer disks. Observations were carried out at the 1.5 m CST (Carlos S\'anchez Telescope) in Tenerife (Spain) using the CAIN infrared camera. 7 galaxies exhibit clear truncation on their disk profiles and 4 galaxies were observed to be clearly untruncated within observational limits. We describe the truncations as real, smooth and complete (as suggested by extrapolation and in the sense that the measured truncation curve goes into the noise at a truncation radius $R_{tr}$), following a decline proportional to $(R_{tr} -R)^{-n}$ (where $R$ is the radius). Despite its deep photometric reach, the data presented do not permit a detailed exploration of the region where optical data show a second slope. Special care was taken concerning the surface brightness deprojection of edge-on galaxies, which was carried out by two methods, one comprising the inversion of Abel's integral equation and the other following a numerical method. These methods gave nearly identical results. NIR observations of truncations could differ from observations in the optical, since the two domains trace different stellar populations.

Following the model of magnetically supported rotation of spiral galaxies, the inner disk rotation is dominated by gravity but magnetism is not negligible at radii where the rotation curve becomes flat, and indeed becomes dominant at very large radii. Values of the order of 1 $\mu$G, or even less, produce a centripetal force when the absolute value of the slope of the curve [$B_\phi$, R] (azimuthal field strength versus radius) is less than $R^{-1}$. The $R^{-1}$-profile is called the critical profile. From this hypothesis, the following is to be expected: at large radii, a ``subcritical'' profile (slope flatter than $R^{-1}$); at still larger radii a $B_\phi$-profile becoming asymptotically critical as the density becomes asymptotically vanishing. Recent observations of magnetic fields in NGC 6946 and the Milky Way are in very good agreement with these predictions. This magnetic alternative requires neither galactic dark matter (DM) nor modification of fundamental laws of physics, but it is not in conflict with these hypotheses, especially with the existence of cosmological cold dark matter (CDM).

We present integral field spectroscopy of NGC 5033, a low luminosity Seyfert galaxy. The observations were made with INTEGRAL, a fiber based system operating at the WHT. The intensity map of the H$\beta$ emission line represents a spiral or ring-like pattern of HII regions. On the contrary, the [OIII] intensity map morphology is markedly anisotropic. The strong morphological differences imply that the [OIII] emitters represent highly ionized gas illuminated by the central source. The [OIII] map morphology is compatible with a biconical structure of ionization induced by strong extinction in the galaxy disc that also obscures half of the spheroidal stellar bulge. We identify the spectrum corresponding to the Seyfert 1 nucleus from the presence of H$\beta$ broad emission lines. This spectrum is located in a region where strong extinction is expected but exhibits the bluest spectral energy distribution. The Seyfert 1 nucleus seems to be offcenter with respect to the stellar rotation center. This result has been also found in other Seyfert galaxies and interpreted in terms of a past merger. The offcentering could indicate the presence of nonsymmetric departures in the gravitational potential which could be fueling the active nucleus. The kinematics of the [OIII] emitters show important deviations at a kpc scale with respect to the stellar velocity field and show features related to the asymmetrical morphology of the high ionization region.

New Integral Field Spectroscopy of the central region of NGC 7331 reveals strong H$\alpha$ emission in the well-known CO and HI ring of NGC 7331. The [NII]/H$\alpha$ ratio indicates that a large scale stellar formation process is taken place at the ring in agreement with previous hypothesis about the exhaustion of gas in the inner to the ring region. The dynamics of stars and gas are not coupled. There is a ring of peculiar velocities in the ionized gas velocity map. These peculiar velocities can be well interpreted by the presence of an axisymmetric inflow of 40 km/s at the inner boundary of the large-scale gaseous ring. We infer an inwards total flux of 1.6 M$_{\odot}$yr$^{-1}$. This value is typical of the accretion rates in hypothetical {\bf large} nuclear black holes. Despite the large differences in the scales of the nucleus and the gas ring of NGC 7331, we suggest that this inwards flux is feeding the nucleus.

We propose a hypothesis of the truncation of stellar discs based on the magnetic model of the rotation curve of spiral galaxies. Once the disc had formed and acquired its present structure, approximately, three balanced forces were acting on the initial gas: gravity and magnetic forces, inwards, and centrifugal force. When stars are formed from this gas, the magnetic force is suddenly suppressed. Gravitation alone cannot retain the newly-formed stars and at birth places beyond a certain galactocentric radius they escape to intergalactic space. This radius is the so-called ``truncation radius'', which is predicted to be at about 4-5 disc radial scale lengths, in promising agreement with observations.

We present a first study of truncation of the stellar disks of spiral galaxies in the near infrared. Observations of NGC4013, NGC4217, NGC6504 and NGC5981 were made with the CAIN NIR camera on the CST in Tenerife. This wavelength range provides the best description of the phenomenon, not only because extinction effects are minimized, but also because the distribution of the old stellar population is directly obtained. The four galaxies are edge-on and an inversion method was developed to obtain the deprojected profiles. We did not assume any model of the different galactic components. The ``truncation curve'', i.e. T(R)= mu(R)- mu_{D}(R), where mu is the actual surface brightness in mag/arcsec^2 and mu_{D} the exponential disk surface brightness, has been obtained with unprecedented precision. It is suggested that T(R) is proportional to (R_{t}-R)^{-1}, where R_{t} is the truncation radius, i.e. the radius beyond which no star is observed.

We review the topic of rotation curves of spiral galaxies emphasizing the standard interpretation as evidence for the existence of dark matter halos. Galaxies other than spirals and late-type dwarfs may also possess great amounts of dark matter, and therefore ellipticals, dwarf spirals, lenticulars and polar ring galaxies are also considered. Furthermore, other methods for determining galactic dark matter, such as those provided by binaries, satellites or globular clusters, have to be included. Cold dark matter hierarchical models constitute the standard way to explain rotation curves, and thus the problem becomes just one aspect of a more general theory explaining structure and galaxy formation. Alternative theories also are included. In the magnetic model, rotation curves could also be a particular aspect of the whole history of cosmic magnetism during different epochs of the Universe. Modifications of Newtonian Dynamics provide another interesting possibility which is discussed here.

Magnetic fields are considered to be dominant when $\epsilon_{B}\geq\epsilon_{K}$, being $\epsilon_{B}=B^{2}/8\pi$ the magnetic energy density and $\epsilon_{K}=1/2 \rho\theta^{2}$ the rotation energy density, for a conventional moderate B= 1 $\mu$G. They are considered to be negligible when $\epsilon_{B}<\epsilon_{K}$ for $B\sim 10 \mu$G. With no assumption and no theoretical calculation, we show that magnetic fields cannot be ignored in the outer parts of a galaxy like the Milky Way and in the whole disk of a dwarf galaxy.

The large-scale structure of the Universe has been observed to be characterized by long filaments, forming polyhedra, with a remarkable 100-200 Mpc periodicity, suggesting a regular network. The introduction of magnetic fields into the physics of the evolution of structure formation provides some clues to understanding this unexpected lattice structure. A relativistic treatment of the evolution of pre-recombination inhomogeneities, including magnetic fields, is presented to show that equivalent-to-present field strengths of the order of $10^{-8}$ G could have played an important role. Primordial magnetic tubes generated at inflation, at scales larger than the horizon before recombination, could have produced filamentary density structures, with comoving lengths larger than about 10 Mpc. Structures shorter than this would have been destroyed by diffusion due to the small pre-recombination conductivity. If filaments constitute a lattice, the primordial magnetic field structures that produced the post-recombination structures of matter, impose several restrictions on the lattice. The simplest lattice compatible with these restrictions is a network of octahedra contacting at their vertexes, which is indeed identifiable in the observed distribution of superclusters.

The distribution of superclusters in the Local Supercluster neighbourhood presents such a remarkable periodicity that some kind of network must fit the observed large scale structure. A three dimension chessboard has been suggested$^{1}$. The existence of this network is really a challenge for currently-suggested theoretical models. For instance, CDM models of the formation of the large scale structure predict a random distribution of superclusters$^{2}$. If the filaments of matter that are now observed building up the network are fossil relics of over-dense regions of magnetic field energy before Recombination, then it has been shown$^{3}$ that the simplest network compatible with magnetic field constraints is made up of octahedra contacting at their vertexes. This suggests a set of superimposed egg-carton structures. Our aim in this paper is to show that the real large-scale structure is actually fitted by the theoretical octahedron structure.

We present a simplified model in which we suggest that two important galactic problems -the magnetic field configuration at large scales and the flat rotation curve- may be simultaneously explained. A highly convective disc produces a high turbulent magnetic diffusion in the vertical direction, stablishing a merging of extragalactic and galactic magnetic fields. The outer disc may then adquire a magnetic energy gradient very close to the gradient required to explain the rotation curve, without the hypothesis of galactic dark matter. Our model predicts symmetries of the galactic field in noticeable agreement with the large scale structure of our galaxy.

The truncation of stellar discs is not abrupt but characterized by a continuous distancing from the exponential profile. There exists a truncation curve, $t(r)$, ending at a truncation radius, $r_t$. We present here a theoretical model in which it is assumed that the magnetic hypothesis explaining the flat rotation curve also explains the truncation. Once stars are born, the centripetal magnetic force previously acting on the progenitor gas cloud is suddenly interrupted, and stars must move to larger orbits or escape. The agreement between theoretical and observed truncation curves is very satisfactory. Parameters defining the disc gas rotation curve should therefore be related to those defining the truncation. It is predicted that rotation curves that quickly reach the asymptotic value $\theta_0 = \theta (r=\infty)$ would have small truncation radii. On the contrary, $r_t$ and $\theta_0$ itself, would be uncorrelated quantities.

In paper I, we obtained an equation for the evolution of density inhomogeneities in a radiation dominated universe when they are affected by magnetic fields. In this second paper we apply this equation to the case in which the subjacent magnetic configuration is a flux tube. For scales of the order of 1 Mpc or less the differential equation is elliptical. To solve it, we have used the numerical method based on "Simultaneous Over Relaxation", SOR, with Chebyshev acceleration and we have treated the problem as a boundary value problem, which restricts the prediction ability of the integration. For large-scale flux tubes, much larger than 1 Mpc, the equation can be analytically integrated and no assumption about the final shape or magnitude of the inhomogeneity is required. In both cases we obtain an evolution which does not differ very much from linear in time. The inhomogeneity in the density becomes filamentary. Large scale structures ($\ge$ 10 Mpc) are probably unaffected by damping, non-linear and amplification mechanisms after Equality, so that this model provides a tool to interpret the present observed large scale structure. Filaments are very frequently found in the large-scale structure in the Universe. It is suggested here that they could arise from primordial magnetic flux tubes, thus providing an alternative hypothesis for its interpretation; in particular we consider the case of the Coma-A1367 supercluster, where the magnetic field is known to be high.

We provide a new tool to interpret the large scale structure of the Universe. As suggested in Paper II, energy density filaments could have been produced by subjacent magnetic flux tubes when the Universe was dominated by radiation. In more recent time epochs, small scale filaments have evolved in a complicated way, but large scale filaments have probably survived and should be identified with present observed matter filaments of clusters. A primordial magnetic field lattice made up of magnetic flux tubes would then have originated a matter lattice made up of cluster filaments. Taking into account some restrictions in the magnetic field configuration we speculate what type of matter network this would have generated. The simplest lattice is made up of octohedra connected at their vertexes, reminding us of a structure of superimposed egg cartons. The vertexes would correspond to large superclusters, which would therefore have eight filaments emerging from them. Some observational evidence is found in the Local Supercluster, which is spider-like, with eight legs.

We consider that no mean magnetic field exists during this epoch, but that there is a mean magnetic energy associated with large-scale magnetic inhomogeneities. We study the evolution of these inhomogeneities and their influence on the large scale density structure, by introducing linear perturbations in Maxwell equations, the conservation of momentum-energy equation, and in Einstein field equations. The primordial magnetic field structure is time independent in the linear approximation, only being diluted by the general expansion, so that $\vec{B}R^2$ is conserved in comoving coordinates. Magnetic fields have a strong influence on the formation of large-scale structure. Firstly, relatively low fields are able to generate density structures even if they were inexistent at earlier times. Second, magnetic fields act anisotropically more recently, modifying the evolution of individual density clouds. Magnetic flux tubes have a tendency to concentrate photons in filamentary patterns.