Context. The magnetic field in spiral galaxies is known to have a large-scale spiral structure along the galactic disk and is observed as X-shaped in the halo of some galaxies. While the disk field can be well explained by dynamo action, the 3-dimensional structure of the halo field and its physical nature is still unclear. Aims. As first steps towards understanding the halo fields, we want to clarify whether the observed X-shaped field is a wide-spread pattern in the halos of spiral galaxies and whether these halo fields are just turbulent fields ordered by compression or shear (anisotropic turbulent fields), or have a large-scale regular structure. Methods. The analysis of the Faraday rotation in the halo is the tool to discern anisotropic turbulent fields from large-scale magnetic fields. This, however, has been challenging until recently because of the faint halo emission in linear polarization. Our sensitive VLA broadband observations C-band and L-band of 35 spiral galaxies seen edge-on (called CHANG-ES) allowed us to perform RM-synthesis in their halos and to analyze the results. We further accomplished a stacking of the observed polarization maps of 28 CHANG-ES galaxies at C-band. Results. Though the stacked edge-on galaxies were of different Hubble types, star formation and interaction activities, the stacked image clearly reveals an X-shaped structure of the apparent magnetic field. We detected a large-scale (coherent) halo field in all 16 galaxies that have extended polarized intensity in their halos. We detected large-scale field reversals in all of their halos. In six galaxies they are along lines about vertical to the galactic midplane (vertical RMTL) with about 2 kpc separation. Only in NGC 3044 and possibly in NGC 3448 we observed vertical giant magnetic ropes (GMRs) similar to those detected recently in NGC 4631.

It is well known that the equation of state (EoS) of compact objects like neutron and quark stars is not determined despite there are several sophisticated models to describe it. From the electromagnetic observations, summarized in \cite{Lattimer01}, and the recent observation of gravitational waves from binary neutron star inspiral GW170817 \cite{Abbott2017_etal} and GW190425 \cite{Abbott2019}, it is possible to make an estimation of the range of masses and so constraint the mass of the neutron and quark stars, determining not only the best approximation for the EoS, but which kind of stars we would be observing. In this paper we explore several configurations of neutron stars assuming a simple polytropic equation of state, using a single layer model without crust. In particular, when the EoS depends on the mass rest density, $p=K \rho_{0}^{\Gamma}$, and when it depends on the energy density $p=K \rho^{\Gamma}$, considerable differences in the mass-radius relationships are found. On the other hand, we also explore quark stars models using the MIT bag EoS for different values of the vacuum energy density $B$.

NGC 4631 is an interacting galaxy which exhibits one of the largest gaseous halos observed among edge-on galaxies. We aim to examine the synchrotron and polarization properties of its disk and halo emission with new radio continuum data. Radio continuum observations of NGC 4631 were performed with the Karl G. Jansky Very Large Array at C-band (5.99 GHz) in the C & D array configurations, and at L-band (1.57 GHz) in the B, C, & D array configurations. The Rotation Measure Synthesis algorithm was utilized to derive the polarization properties. We detected linearly polarized emission at C-band and L-band. The magnetic field in the halo is characterized by strong vertical components above and below the central region of the galaxy. The magnetic field in the disk is only clearly seen in the eastern side of NGC 4631, where it is parallel to the plane of the major axis of the galaxy. We detected for the first time a large-scale, smooth Faraday depth pattern in a halo of an external spiral galaxy, which implies the existence of a regular (coherent) magnetic field. A quasi-periodic pattern in Faraday depth with field reversals was found in the northern halo of the galaxy. The field reversals in the northern halo of NGC 4631, together with the observed polarization angles, indicate giant magnetic ropes (GMRs) with alternating directions. To our knowledge, this is the first time such reversals are observed in an external galaxy.

NGC 4631 is an interacting galaxy that exhibits one of the largest, gaseous halos observed among edge-on galaxies. We aim to examine the synchrotron and cosmic-ray propagation properties of its disk and halo emission with new radio continuum data. Radio continuum observations of NGC 4631 were performed with the Karl G. Jansky Very Large Array at C-band (5.99 GHz) in the C and D array configurations, and at L-band (1.57 GHz) in the B, C, and D array configurations. Complementary observations of NGC 4631 with the Effelsberg telescope were performed at 1.42 and 4.85 GHz. The interferometric total intensity data were combined with the single-dish Effelsberg data in order to recover the missing large-scale total power emission. The thermal and nonthermal components of the total radio emission were separated by estimating the thermal contribution through the extinction-corrected H$\alpha$ emission. The H$\alpha$ radiation was corrected for extinction using a linear combination of the observed H$\alpha$ and 24 $\mu$m data. NGC 4631 has a global thermal fraction at 5.99 (1.57) GHz of 14$\pm$3% (5.4$\pm$1.1%). The mean scale heights of the total emission in the radio halo (thick disk) at 5.99 (1.57) GHz are $1.79\pm0.54$ kpc ($1.75\pm0.27$ kpc) and have about the same values for the synchrotron emission. The total magnetic field of NGC 4631 has a mean strength of $\rm{\langle B_{eq}\rangle} \simeq 9 \rm{\mu G}$ in the disk, and a mean strength of $\rm{\langle B_{eq}\rangle}~\simeq 7~\rm{\mu G}$ in the halo. We also studied a double-lobed background radio galaxy southwest of NGC 4631, which is an FR~II radio galaxy according to the distribution of spectral index across the lobes. From the halo scale heights we estimated that the radio halo is escape-dominated with convective cosmic ray propagation, and conclude that there is a galactic wind in the halo of NGC 4631.

We present a systematic study of the extraplanar gas (EPG) in a sample of 15 nearby late-type galaxies at intermediate inclinations using publicly available, deep interferometric HI data from the HALOGAS survey. For each system we mask the HI emission coming from the regularly rotating disc and use synthetic datacubes to model the leftover "anomalous" HI flux. Our model consists of a smooth, axisymmetric thick component described by 3 structural and 4 kinematical parameters, which are fit to the data via a Bayesian MCMC approach. We find that extraplanar HI is nearly ubiquitous in disc galaxies, as we fail to detect it in only two of the systems with the poorest spatial resolution. The EPG component encloses ~5-25% of the total HI mass, with a mean value of 14%, and has a typical thickness of a few kpc, incompatible with expectations based on hydrostatic equilibrium models. The EPG kinematics is remarkably similar throughout the sample, and consists of a lagging rotation with typical vertical gradients of about -10 km/s/kpc, a velocity dispersion of 15-30 km/s and, for most galaxies, a global inflow in both the vertical and radial directions with speeds of 20-30 km/s. The EPG HI masses are in excellent agreement with predictions from simple models of the galactic fountain powered by stellar feedback. The combined effect of photo-ionisation and interaction of the fountain material with the circumgalactic medium can qualitatively explain the kinematics of the EPG, but dynamical models of the galactic fountain are required to fully test this framework.

We present new narrow-band H-alpha imaging for 24 nearby edge-on galaxies in the CHANG-ES survey. We use the images in conjunction with WISE 22 micron imaging of the sample to estimate improved star formation rates (SFRs) using the updated recipe from Vargas et al. (2018). We explore correlations between the updated star formation properties and radio continuum scale heights, scale lengths, and diameters, measured in Krause et al. (2018). We find a newly discovered correlation between SFR and radio scale height that did not exist using mid-IR only SFR calibrations. This implies that a mid-IR extinction correction should be applied to SFR calibrations when used in edge-on galaxies, due to attenuation by dust. The updated SFR values also show newly discovered correlations with radio scale length and radio diameter, implying that the previously-measured relationship between radio scale height and radio diameter originates from star formation within the disk. We also identify a region of star formation located at extreme distance from the disk of NGC 4157, possibly ionized by a single O5.5 V star. This region is spatially coincident with an XUV disk feature, as traced by GALEX NUV imaging. We theorize that the star formation feature arose due to gravitational instability within gas from an accretion event. New H-alpha images from this work can be found at the CHANG-ES data release web site, https://www.queensu.ca/changes.

In this paper the problem of consistency of smoothed particle hydrodynamics (SPH) is solved. A novel error analysis is developed in $n$-dimensional space using the Poisson summation formula, which enables the treatment of the kernel and particle approximation errors in combined fashion. New consistency integral relations are derived for the particle approximation which correspond to the cosine Fourier transform of the classically known consistency conditions for the kernel approximation. The functional dependence of the error bounds on the SPH interpolation parameters, namely the smoothing length $h$ and the number of particles within the kernel support ${\cal{N}}$ is demonstrated explicitly from which consistency conditions are seen to follow naturally. As ${\cal{N}}\to\infty$, the particle approximation converges to the kernel approximation independently of $h$ provided that the particle mass scales with $h$ as $m\propto h^{\beta}$, with $\beta >n$. This implies that as $h\to 0$, the joint limit $m\to 0$, ${\cal{N}}\to\infty$, and $N\to\infty$ is necessary for complete convergence to the continuum, where $N$ is the total number of particles. The analysis also reveals the presence of a dominant error term of the form $(\ln {\cal{N}})^{n}/{\cal{N}}$, which tends asymptotically to $1/{\cal{N}}$ when ${\cal{N}}\gg 1$, as it has long been conjectured based on the similarity between the SPH and the quasi-Monte Carlo estimates.

We analyze the application of star formation rate (SFR) calibrations using H$\alpha$ and 22 micron infrared imaging data in predicting the thermal radio component for a test sample of 3 edge-on galaxies (NGC 891, NGC 3044, and NGC 4631) in the Continuum Halos in Nearby Galaxies -- an EVLA Survey (CHANG-ES). We use a mixture of H$\alpha$ and 24 micron calibration from Calzetti et al. (2007), and a linear 22 micron only calibration from Jarrett et al. (2013) on the test sample. We apply these relations on a pixel-to-pixel basis to create thermal prediction maps in the two CHANG-ES bands: L- and C-band (1.5 GHz and 6.0 GHz, respectively). We analyze the resulting non-thermal spectral index maps, and find a characteristic steepening of the non-thermal spectral index with vertical distance from the disk after application of all methods. We find possible evidence of extinction in the 22 micron data as compared to 70 micron Spitzer Multband Imaging Photometer (MIPS) imaging in NGC 891. We analyze a larger sample of edge-on and face-on galaxy 25 micron to 100 micron flux ratios, and find that the ratios for edge-ons are systematically lower by a factor of 1.36, a result we attribute to excess extinction in the mid-IR in edge-ons. We introduce a new calibration for correcting the H$\alpha$ luminosity for dust when galaxies are edge-on or very dusty.

The vertical halo scale height is a crucial parameter to understand the transport of cosmic-ray electrons (CRE) and their energy loss mechanisms in spiral galaxies. Until now, the radio scale height could only be determined for a few edge-on galaxies because of missing sensitivity at high resolution. We developed a sophisticated method for the scale height determination of edge-on galaxies. With this we determined the scale heights and radial scale lengths for a sample of 13 galaxies from the CHANG-ES radio continuum survey in two frequency bands. The sample average value for the radio scale heights of the halo are 1.1 +/- 0.3kpc in C-band and 1.4 +/- 0.7kpc in L-band. From the frequency dependence analysis of the halo scale heights we found that the wind velocities (estimated using the adiabatic loss time) are above the escape velocity. We found that the halo scale heights increase linearly with the radio diameters. In order to exclude the diameter dependence, we defined a normalized scale height which is quite similar for all sample galaxies at both frequency bands and does not depend on the star formation rate or the magnetic field strength. However, the normalized scale height shows a tight anticorrelation with the mass surface density. The sample galaxies with smaller scale lengths are more spherical in the radio emission, while those with larger scale lengths are flatter. The radio scale height depends mainly on the radio diameter of the galaxy. The sample galaxies are consistent with an escape-dominated radio halo with convective cosmic ray propagation, indicating that galactic winds are a widespread phenomenon in spiral galaxies. While a higher star formation rate or star formation surface density does not lead to a higher wind velocity, we deceleration of CRE outflow, e.g. a lowering of the wind velocity from the galactic disk.

We use new deep 21 cm HI observations of the moderately inclined galaxy NGC 4559 in the HALOGAS survey to investigate the properties of extra-planar gas. We use TiRiFiC to construct simulated data cubes to match the HI observations. We find that a thick disk component of scale height $\sim\,2\,\mathrm{kpc}$, characterized by a negative vertical gradient in its rotation velocity (lag) of $\sim13 \pm 5$ km s$^{-1}$ kpc$^{-1}$ is an adequate fit to extra-planar gas features. The tilted ring models also present evidence for a decrease in the magnitude of the lag outside of $R_{25}$, and a radial inflow of $\sim 10$ km s$^{-1}$. We extracted lagging extra-planar gas through Gaussian velocity profile fitting. From both the 3D models and and extraction analyses we conclude that $\sim10-20\%$ of the total {\HI} mass is extra-planar. Most of the extra-planar gas is spatially coincident with regions of star formation in spiral arms, as traced by H$\alpha$ and GALEX FUV images, so it is likely due to star formation processes driving a galactic fountain. We also find the signature of a filament of a kinematically "forbidden" HI, containing $\sim 1.4\times 10^{6}$ M$_{\odot}$ of HI, and discuss its potential relationship to a nearby HI hole. We discover a previously undetected dwarf galaxy in HI located $\sim 0.4^{\circ}$ ($\sim 58$ kpc) from the center of NGC 4559, containing $\sim 4\times10^{5}$ M$_{\odot}$. This dwarf has counterpart sources in SDSS with spectra typical of HII regions, and we conclude it is two merging blue compact dwarf galaxies.

The problem of consistency of smoothed particle hydrodynamics (SPH) has demanded considerable attention in the past few years due to the ever increasing number of applications of the method in many areas of science and engineering. A loss of consistency leads to an inevitable loss of approximation accuracy. In this paper, we revisit the issue of SPH kernel and particle consistency and demonstrate that SPH has a limiting second-order convergence rate. Numerical experiments with suitably chosen test functions validate this conclusion. In particular, we find that when using the root mean square error as a model evaluation statistics, well-known corrective SPH schemes, which were thought to converge to second, or even higher order, are actually first-order accurate, or at best close to second order. We also find that observing the joint limit when $N\to\infty$, $h\to 0$, and $n\to\infty$, as was recently proposed by Zhu et al., where $N$ is the total number of particles, $h$ is the smoothing length, and $n$ is the number of neighbor particles, standard SPH restores full $C^{0}$ particle consistency for both the estimates of the function and its derivatives and becomes insensitive to particle disorder.

Based on the Lema\^itre-Tolman-Bondi (LTB) metric we consider two flat inhomogeneous big-bang models. We aim at clarifying, as far as possible analytically, basic features of the dynamics of the simplest inhomogeneous models and to point out the potential usefulness of exact inhomogeneous solutions as generalizations of the homogeneous configurations of the cosmological standard model. We discuss explicitly partial successes but also potential pitfalls of these simplest models. Although primarily seen as toy models, the relevant free parameters are fixed by best-fit values using the Joint Light-curve Analysis (JLA)-sample data. On the basis of a likelihood analysis we find that a local hump with an extension of almost 2 Gpc provides a better description of the observations than a local void for which we obtain a best-fit scale of about 30 Mpc. Future redshift-drift measurements are discussed as a promising tool to discriminate between inhomogeneous configurations and the $\Lambda$CDM model.

In the present work we constrain three different profiles of a Lema\^itre-Tolman-Bondi model using supernovae type Ia and baryon acoustic oscillation data. We improve common practice in the literature by carefully calibrating the supernovae in the appropriate inhomogeneous background dynamics. In addition, we address subtle issues in order to propagate the primordial BAO scale to present epoch. The combined analysis of BAO+SNIa offers a stringent test for these models. We use two distinct parameter estimation approaches, namely, the $\chi^2$ and the complete likelihood functional. It has been argued that these two approaches are not equivalent and indeed our analysis shows a specific example of their departure.

How a galaxy regulates its SNe energy into different interstellar/circumgalactic medium components strongly affects galaxy evolution. Based on the JVLA D-configuration C- (6 GHz) and L-band (1.6 GHz) continuum observations, we perform statistical analysis comparing multi-wavelength properties of the CHANG-ES galaxies. The high-quality JVLA data and edge-on orientation enable us for the first time to include the halo into the energy budget for a complete radio-flux-limited sample. We find tight correlations of $L_{\rm radio}$ with the mid-IR-based SFR. The normalization of our $I_{\rm 1.6GHz}/{\rm W~Hz^{-1}}-{\rm SFR}$ relation is $\sim$2-3 times of those obtained for face-on galaxies, probably a result of enhanced IR extinction at high inclination. We also find tight correlations between $L_{\rm radio}$ and the SNe energy injection rate $\dot{E}_{\rm SN(Ia+CC)}$, indicating the energy loss via synchrotron radio continuum accounts for $\sim0.1\%$ of $\dot{E}_{\rm SN}$, comparable to the energy contained in CR electrons. The integrated C-to-L-band spectral index is $\alpha\sim0.5-1.1$ for non-AGN galaxies, indicating a dominance by the diffuse synchrotron component. The low-scatter $L_{\rm radio}-{\rm SFR}$/$L_{\rm radio}-\dot{E}_{\rm SN (Ia+CC)}$ relationships have super-linear logarithmic slopes at $\sim2~\sigma$ in L-band ($1.132\pm0.067$/$1.175\pm0.102$) while consistent with linear in C-band ($1.057\pm0.075$/$1.100\pm0.123$). The super-linearity could be naturally reproduced with non-calorimeter models for galaxy disks. Using Chandra halo X-ray measurements, we find sub-linear $L_{\rm X}-L_{\rm radio}$ relations. These results indicate that the observed radio halo of a starburst galaxy is close to electron calorimeter, and a galaxy with higher SFR tends to distribute an increased fraction of SNe energy into radio emission (than X-ray).

We use the Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey (CANDELS) GOODS-S multi-wavelength catalog to identify counterparts for 20 Lya Emitting (LAE) galaxies at z=2.1. We build several types of stacked Spectral Energy Distributions (SEDs) of these objects. We combine photometry to form average and median flux-stacked SEDs, and postage stamp images to form average and median image-stacked SEDs. We also introduce scaled flux stacks that eliminate the influence of variation in overall brightness. We use the SED fitting code SpeedyMC to constrain the physical properties of individual objects and stacks. Our LAEs at z = 2.1 have stellar masses ranging from 2x10^7 Msun - 8x10^9 Msun (median = 3x10^8 Msun), ages ranging from 4 Myr to 500 Myr (median =100 Myr), and E(B-V) between 0.02 and 0.24 (median = 0.12). We do not observe strong correlations between Lya equivalent width (EW) and stellar mass, age, or E(B-V). The Lya radiative transfer (q) factors of our sample are predominantly close to one and do not correlate strongly with EW or E(B-V), implying that Lya radiative transfer prevents Lya photons from resonantly scattering in dusty regions. The SED parameters of the flux stacks match the average and median values of the individual objects, with the flux-scaled median SED performing best with reduced uncertainties. Median image-stacked SEDs provide a poor representation of the median individual object, and none of the stacking methods captures the large dispersion of LAE properties.

We argue that interactions in the dark sector may have a crucial impact on the cosmological dynamics. In particular, the future cosmic evolution may be very different from that predicted by the LCDM model. An example is a scenario in which the currently observed accelerated expansion is an interaction-induced transient phenomenon. We discuss such type of behavior on the basis of a two-fluid toy model.

Lyman Alpha Emitting (LAE) galaxies are thought to be progenitors of present-day L* galaxies. Clustering analyses have suggested that LAEs at z ~ 3 might evolve into LAEs at z ~ 2, but it is unclear whether the physical nature of these galaxies is compatible with this hypothesis. Several groups have investigated the properties of LAEs using spectral energy distribution (SED) fitting, but direct comparison of their results is complicated by inconsistencies in the treatment of the data and in the assumptions made in modeling the stellar populations, which are degenerate with the effects of galaxy evolution. By using the same data analysis pipeline and SED fitting software on two stacked samples of LAEs at z = 3.1 and z = 2.1, and by eliminating several systematic uncertainties that might cause a discrepancy, we determine that the physical properties of these two samples of galaxies are dramatically different. LAEs at z = 3.1 are found to be old (age ~ 1 Gyr) and metal-poor (Z < 0.2 Z_Sun), while LAEs at z = 2.1 appear to be young (age ~ 50 Myr) and metal-rich (Z > Z_Sun). The difference in the observed stellar ages makes it very unlikely that z = 3.1 LAEs evolve directly into z = 2.1 LAEs. Larger samples of galaxies, studies of individual objects and spectroscopic measurements of metallicity at these redshifts are needed to confirm this picture, which is difficult to reconcile with the effects of 1 Gyr of cosmological evolution.

According to the standard $\Lambda$CDM model, the accelerated expansion of the Universe will go on forever. Motivated by recent observational results, we explore the possibility of a finite phase of acceleration which asymptotically approaches another period of decelerated expansion. Extending an earlier study on a corresponding homogeneous and isotropic dynamics, in which interactions between dark matter and dark energy are crucial, the present paper also investigates the dynamics of the matter perturbations both on the Newtonian and GR levels and quantifies the potential relevance of perturbations of the dark-energy component. In the background, the model is tested against the Supernova type Ia (SNIa) data of the Constitution set and on the perturbative level against growth rate data and the data of the 2dFGRS project. Our results indicate, that a transient phase of accelerated expansion is not excluded by current observations.

This paper describes the Hubble Space Telescope imaging data products and data reduction procedures for the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). This survey is designed to document the evolution of galaxies and black holes at $z\sim1.5-8$, and to study Type Ia SNe beyond $z>1.5$. Five premier multi-wavelength sky regions are selected, each with extensive multiwavelength observations. The primary CANDELS data consist of imaging obtained in the Wide Field Camera 3 / infrared channel (WFC3/IR) and UVIS channel, along with the Advanced Camera for Surveys (ACS). The CANDELS/Deep survey covers \sim125 square arcminutes within GOODS-N and GOODS-S, while the remainder consists of the CANDELS/Wide survey, achieving a total of \sim800 square arcminutes across GOODS and three additional fields (EGS, COSMOS, and UDS). We summarize the observational aspects of the survey as motivated by the scientific goals and present a detailed description of the data reduction procedures and products from the survey. Our data reduction methods utilize the most up to date calibration files and image combination procedures. We have paid special attention to correcting a range of instrumental effects, including CTE degradation for ACS, removal of electronic bias-striping present in ACS data after SM4, and persistence effects and other artifacts in WFC3/IR. For each field, we release mosaics for individual epochs and eventual mosaics containing data from all epochs combined, to facilitate photometric variability studies and the deepest possible photometry. A more detailed overview of the science goals and observational design of the survey are presented in a companion paper.