We explore a new kind of field of nonlinear electrodynamics(NLED) which acts as a source of gravity and can accelerate the universe during the inflationary era. We propose a new type of NLED lagrangian which is charecterized by two paremetrs $\alpha$ and $\beta$. We investigate the classical stability and causality aspects of this model by demanding that the speed ($C_s = \frac{dP}{d\rho}$) of the sound wave $C_s^2 > 0$ and $C_s \le 1$ and find that $0 < C_s^2 < 1$ corresponds to $0.25 \le \alpha \le 0.4$ and $0.6 \le \beta B^2 \le 1$. A study of the deceleration parameter($q = \frac{1}{2}(1 + 3 \omega)$, $\omega = P/\rho$ being the equation of state parameter) suggests that the value $q < 0$ (i.e. $\omega < -1/3$ and $\ddot{a}(t) > 0$ ( the accelerating universe)) requires $\beta B^2 \ge 0.13$. During inflation, the energy density $\rho_B$ is found to be maximum and is given by $\rho_B^{max} = 0.65/\beta$ corresponding to $\alpha = 0.3$. The magnetic field necessary to trigger the inflation, is found to be $B (= B_{max}) \simeq \sqrt{\frac{0.4 \rho_B^{max}}{0.65}} = 4 \times 10^{51}~{\rm Gauss}$, where $\rho_{B}^{max}( =10^{64}~{\rm GeV}^4 $) is the energy density of the universe during inflation. The model also predict the e-fold number $N = 71$, which agrees with the experimental result.

We consider the size distribution of superbubbles in a star forming galaxy. Previous studies have tried to explain the distribution by using adiabatic self-similar evolution of wind driven bubbles, assuming that bubbles stall when pressure equilibrium is reached. We show, with the help of hydrodynamical numerical simulations, that this assumption is not valid. We also include radiative cooling of shells. In order to take into account non-thermal pressure in the ambient medium, we assume an equivalent higher temperature than implied by thermal pressure alone. Assuming that bubbles stall when the outer shock speed becomes comparable to the ambient sound speed (which includes non-thermal components), we recover the size distribution with a slope of $\sim -2.7$ for typical values of ISM pressure in Milky Way, which is consistent with observations. Our simulations also allow us to follow the evolution of size distribution in the case of different values of non-thermal pressure, and we show that the size distribution steepens with lower pressure, to slopes intermediate between only-growing and only-stalled cases.

Selection functions are vital for understanding the observational biases of spectroscopic surveys. With the wide variety of multi-object spectrographs currently in operation and becoming available soon, we require easily generalisable methods for determining the selection functions of these surveys. Previous work, however, has largely been focused on generating individual, tailored selection functions for every data release of each survey. Moreover, no methods for combining these selection functions to be used for joint catalogues have been developed. We have developed a Poisson likelihood estimation method for calculating selection functions in a Bayesian framework, which can be generalised to any multi-object spectrograph. We include a robust treatment of overlapping fields within a survey as well as selection functions for combined samples with overlapping footprints. We also provide a method for transforming the selection function that depends on the sky positions, colour, and apparent magnitude of a star to one that depends on the galactic location, metallicity, mass, and age of a star. This `intrinsic' selection function is invaluable for chemodynamical models of the Milky Way. We demonstrate that our method is successful at recreating synthetic spectroscopic samples selected from a mock galaxy catalogue.

We exploit the [Mg/Mn]-[Al/Fe] chemical abundance plane to help identify nearby halo stars in the 14th data release from the APOGEE survey that have been accreted on to the Milky Way. Applying a Gaussian Mixture Model, we find a `blob' of 856 likely accreted stars, with a low disc contamination rate of ~7%. Cross-matching the sample with the second data release from Gaia gives us access to parallaxes and apparent magnitudes, which place constraints on distances and intrinsic luminosities. Using a Bayesian isochrone pipeline, this enables us to estimate new ages for the accreted stars, with typical uncertainties of ~20%. Our new catalogue is further supplemented with estimates of orbital parameters. The blob stars span a metallicities between -0.5 to -2.5, and [Mg/Fe] between -0.1 to 0.5. They constitute ~30% of the metal-poor ([Fe/H] < -0.8) halo at metallicities of ~-1.4. Our new ages are mainly range between 8 to 13 Gyr, with the oldest stars the metal-poorest, and with the highest [Mg/Fe] abundance. If the blob stars are assumed to belong to a single progenitor, the ages imply that the system merged with our Milky Way around 8 Gyr ago and that star formation proceeded for ~5 Gyr. Dynamical arguments suggest that such a single progenitor would have a total mass of ~1011Msun, similar to that found by other authors using chemical evolution models and simulations. Comparing the scatter in the [Mg/Fe]-[Fe/H] plane of the blob stars to that measured for stars belonging to the Large Magellanic Cloud suggests that the blob does indeed contain stars from only one progenitor.

We derive model-independent astrophysical constraints on leptophilic dark matter (DM), considering its thermal production in a supernova core and taking into account core temperature fluctuations within the framework of $q$-deformed Tsallis statistics. In an effective field theory approach, where the DM fermions interact with the Standard Model via dimension-six operators of either scalar-pseudoscalar, vector-axial vector, or tensor-axial tensor type, we obtain bounds on the effective cut-off scale $\Lambda$ from supernova cooling and free-streaming of DM from supernova core, and from thermal relic density considerations, depending on the DM mass and the $q$-deformation parameter. Using Raffelt's criterion on the energy loss rate from SN1987A, we obtain a lower bound on $\Lda \gtrsim 3$ (12) TeV corresponding to $q = 1.0~(1.1)$ and an average supernova core temperature of $T_{\rm SN}=30$ MeV. From the optical depth criterion on the free-streaming of DM fermions from the outer 10\% of the SN1987A core, the cooling bound is restricted to $\Lda \gtrsim 1$ TeV. Both cooling and free-streaming bounds are insensitive to the DM mass $m_\chi$ for $m_\chi\lesssim T_{\rm SN}$, whereas for $m_\chi\gg T_{\rm SN}$, the bounds weaken significantly due to the Boltzmann-suppression of the DM number density. We also calculate the thermal relic density of the DM particles in this setup and find that it imposes an upper bound on $\Lambda^4/m_\chi^2$, which together with the cooling/free-streaming bound significantly constrains light leptophilic DM.

We present a new approach (MADE) that generates mass, age, and distance estimates of red giant stars from a combination of astrometric, photometric, and spectroscopic data. The core of the approach is a Bayesian artificial neural network (ANN) that learns from and completely replaces stellar isochrones. The ANN is trained using a sample of red giant stars with mass estimates from asteroseismology. A Bayesian isochrone pipeline uses the astrometric, photometric, spectroscopic, and asteroseismology data to determine posterior distributions for the training outputs: mass, age, and distance. Given new inputs, posterior predictive distributions for the outputs are computed, taking into account both input uncertainties, and uncertainties in the ANN parameters. We apply MADE to $\sim10\,000$ red giants in the overlap between the 14$^{\mathrm{th}}$ data release from the APO Galactic Evolution Experiment (APOGEE, Abolfathi et al. 2018) and the Tycho-Gaia astrometric solution (TGAS, Michalik et al. 2015). The ANN is able to reduce the uncertainty on mass, age, and distance estimates for training-set stars with high output uncertainties allocated through the Bayesian isochrone pipeline. The fractional uncertainties on mass are $<10\%$ and on age are between $10$ to $25\%$. Moreover, the time taken for our ANN to predict masses, ages, and distances for the entire catalogue of APOGEE-TGAS stars is of a similar order of the time taken by the Bayesian isochrone pipeline to run on a handful of stars. Our resulting catalogue clearly demonstrates the expected thick and thin disc components in the [M/H]-[$\alpha$/M] plane, when examined by age.

We present a catalogue of distances, masses and ages for $\sim3$ million stars in the second Gaia data release with spectroscopic parameters available from the large spectroscopic surveys: APOGEE, Gaia-ESO, GALAH, LAMOST, RAVE and SEGUE. We use a Bayesian framework to characterise the probability density functions of distance, mass and age using photometric, spectroscopic and astrometric information, supplemented with spectroscopic masses where available for giant stars. Furthermore, we provide posterior extinction estimates ($A_V$) to every star using published extinction maps as a prior input. We provide an appendix with extinction coefficients for Gaia photometry derived from stellar models, which account for variation with intrinsic colour and total extinction. Our pipeline provides output estimates of the spectroscopic parameters, which can be used to inform improved spectroscopic analysis. We complement our catalogues with Galactocentric coordinates and actions with associated uncertainties. As a demonstration of the power of our catalogue, we produce velocity dispersion profiles of the disc separated by age and Galactocentric radius (between $3$ and $15\,\mathrm{kpc}$ from the Galactic centre). This suggests that the velocity dispersion profiles flatten with radius in the outer Galaxy ($>8\,\mathrm{kpc}$) and that at all radii the velocity dispersion follows the smooth power law with age observed in the solar neighbourhood.

The cold dark matter fermions with mass MeV scale, pair produced inside the supernova SN1987A core, can freely stream away from the supernovae and hence contributes to its energy loss rate. Similar type of DM fermions(having similar kind of coupling to the standard model photon), produced from some other sources earlier, could have contributed to the relic density of the Universe. Working in a theory with an effective dark matter-photon coupling (inversely proportional to the scale $\Lambda$) in the formalism of Tsallis statistics, we find the dark matter contribution to the relic density and obtain a upper bound on $\Lambda$ using the experimental bound on the relic density for cold non-baryonic matter i.e. $\Omega h^2 = 0.1186 \pm 0.0020 $. The upper bound obtained from the relic density is shown with the lower bound obtained from the Raffelt's criterion on the emissibity rate of the supernovae SN1987A energy loss $\dot{\varepsilon}(e^+ e^- \to \chi \overline{\chi}) \le 10^{19}~\rm{erg~g^{-1}s^{-1}}$ and the optical depth criteria on the free streaming of the dark matter fermion (produced inside the supernovae core). As the deformation parameter $q$ changes from $1.0$ (undeformed scenario) to $1.1$(deformed scenario), the relic density bound on $\Lambda$ is found to vary from $ \sim 4.9 \times 10^7 $ TeV to $1.6 \times 10^8$ TeV for a fermion dark matter($\chi$) of mass $m_\chi = 30~\rm{MeV}$, which is almost $10$ times more than the lower bound obtained from the SN1987A energy loss rate and the optical depth criteria. \noindent {{\bf Keywords}: Dark matter, Relic density, Supernova cooling, Tsallis statistics, free-streaming, } }

Phylogenetics is a widely used concept in evolutionary biology. It is the reconstruction of evolutionary history by building trees that represent branching patterns and sequences. These trees represent shared history, and it is our intention for this approach to be employed in the analysis of Galactic history. In Galactic archaeology the shared environment is the interstellar medium in which stars form and provides the basis for tree-building as a methodological tool. Using elemental abundances of solar-type stars as a proxy for DNA, we built in Jofre et al 2017 such an evolutionary tree to study the chemical evolution of the solar neighbourhood. In this proceeding we summarise these results and discuss future prospects.

In the context of the ESA M5 (medium mission) call we proposed a new satellite mission, Theia, based on relative astrometry and extreme precision to study the motion of very faint objects in the Universe. Theia is primarily designed to study the local dark matter properties, the existence of Earth-like exoplanets in our nearest star systems and the physics of compact objects. Furthermore, about 15 $\%$ of the mission time was dedicated to an open observatory for the wider community to propose complementary science cases. With its unique metrology system and "point and stare" strategy, Theia's precision would have reached the sub micro-arcsecond level. This is about 1000 times better than ESA/Gaia's accuracy for the brightest objects and represents a factor 10-30 improvement for the faintest stars (depending on the exact observational program). In the version submitted to ESA, we proposed an optical (350-1000nm) on-axis TMA telescope. Due to ESA Technology readiness level, the camera's focal plane would have been made of CCD detectors but we anticipated an upgrade with CMOS detectors. Photometric measurements would have been performed during slew time and stabilisation phases needed for reaching the required astrometric precision.

Using 17 chemical elements as a proxy for stellar DNA, we present a full phylogenetic study of stars in the solar neighbourhood. This entails applying a clustering technique that is widely used in molecular biology to construct an evolutionary tree from which three branches emerge. These are interpreted as stellar populations which separate in age and kinematics and can be thus attributed to the thin disk, the thick disk, and an intermediate population of probable distinct origin. We further find six lone stars of intermediate age that could not be assigned to any population with enough statistical significance. Combining the ages of the stars with their position on the tree, we are able to quantify the mean rate of chemical enrichment of each of the populations, and thus show in a purely empirical way that the star formation rate in the thick disk is much higher than in the thin disk. We are also able to estimate the relative contribution of dynamical processes such as radial migration and disk heating to the distribution of chemical elements in the solar neighbourhood. Our method offers an alternative approach to chemical tagging methods with the advantage of visualising the behaviour of chemical elements in evolutionary trees. This offers a new way to search for `common ancestors' that can reveal the origin of solar neighbourhood stars.

Light dark matter($\simeq 1-30~\rm{MeV}$) particles pair produced in electron-positron annihilation $ e^-e^+ \stackrel{\gamma}{\longrightarrow} \chi \bar{\chi}$ inside the supernova core can take away the energy released in the supernova SN1987A explosion. Working within the formalism of $q$-deformed statistics (with the average value of the supernovae core temperature(fluctuating) being $T_{SN} = 30~\rm{MeV}$) and using the Raffelt's criterion on the emissivity for any new channel $\dot{\varepsilon}(e^+ e^- \to \chi \overline{\chi}) \le 10^{19}~{erg~g^{-1}s^{-1}}$, we find that as the deformation parameter $q$ changes from $1.0$ (undeformed scenario) to $1.1$(deformed scenario), the lower bound on the scale $\Lambda$ of the dark matter effective theory varies from $3.3\times 10^6$ TeV to $3.2 \times 10^7$ TeV for a dark matter fermion of mass $m_\chi = 30~\rm{MeV}$. Using the optical depth criteria on the free streaming of the dark matter fermion, we find the lower bound on $\Lambda \sim 10^{8}~\rm{TeV}$ for $m_\chi = 30~\rm{MeV}$. In a scenerio,where the dark matter fermions are pair produced in the outermost sector of the supernova core (with radius $0.9 R_c \le r \le R_c$, $R_c (=10~\rm{km})$ being the supernova core radius or the radius of proto-neutron star), we find that the bound on $\Lambda$ ($\sim 3 \times 10^7$ TeV)obtained from SN cooling criteria (Raffelt's criteria) is comparable with the bound obtained from free streaming (optical depth criterion) for light fermion dark matter of mass $m_{\chi}=10 - 30$ MeV.

The distribution of Milky Way halo blue horizontal-branch (BHB) stars is examined using action-based extended distribution functions (EDFs) that describe the locations of stars in phase space, metallicity, and age. The parameters of the EDFs are fitted using stars observed in the Sloan Extension for Galactic Understanding and Exploration-II (SEGUE-II) survey that trace the phase-space kinematics and chemistry out to ~70 kpc. A maximum a posteriori probability (MAP) estimate method and a Markov Chain Monte Carlo method are applied, taking into account the selection function in positions, distance, and metallicity for the survey. The best-fit EDF declines with actions less steeply at actions characteristic of the inner halo than at the larger actions characteristic of the outer halo, and older ages are found at smaller actions than at larger actions. In real space, the radial density profile steepens smoothly from -2 at ~2 kpc to -4 in the outer halo, with an axis ratio ~0.7 throughout. There is no indication for rotation in the BHBs, although this is highly uncertain. A moderate level of radial anisotropy is detected, with $\beta_s$ varying from isotropic to between ~0.1 and ~0.3 in the outer halo depending on latitude. The BHB data are consistent with an age gradient of -0.03 Gyr kpc$^{-1}$, with some uncertainty in the distribution of the larger ages. These results are consistent with a scenario in which older, larger systems contribute to the inner halo, whilst the outer halo is primarily comprised of younger, smaller systems.

We fit an Extended Distribution Function (EDF) to K giants in the Sloan Extension for Galactic Understanding and Exploration (SEGUE) survey. These stars are detected to radii ~80 kpc and span a wide range in [Fe/H]. Our EDF, which depends on [Fe/H] in addition to actions, encodes the entanglement of metallicity with dynamics within the Galaxy's stellar halo. Our maximum-likelihood fit of the EDF to the data allows us to model the survey's selection function. The density profile of the K giants steepens with radius from a slope ~-2 to ~-4 at large radii. The halo's axis ratio increases with radius from 0.7 to almost unity. The metal-rich stars are more tightly confined in action space than the metal-poor stars and form a more flattened structure. A weak metallicity gradient ~-0.001 dex/kpc, a small gradient in the dispersion in [Fe/H] of ~0.001 dex/kpc, and a higher degree of radial anistropy in metal-richer stars result. Lognormal components with peaks at ~-1.5 and ~-2.3 are required to capture the overall metallicity distribution, suggestive of the existence of two populations of K giants. The spherical anisotropy parameter varies between 0.3 in the inner halo to isotropic in the outer halo. If the Sagittarius stream is included, a very similar model is found but with a stronger degree of radial anisotropy throughout.

In this paper, we extend our earlier study of spontaneous CP violating scattering of quarks and anti-quarks from QCD $Z\left(3\right)$ domain walls for the situation when these walls have asymmetric profiles of the Polyakov loop order parameter $l(x)$. Dynamical quarks lead to explicit breaking of $Z(3)$ symmetry, which lifts the degeneracy of the $Z(3)$ vacua arising from spontaneous breaking of the $Z(3)$ symmetry in the quark-gluon plasma (QGP) phase. Resulting domain walls have asymmetric profile of $l(x)$ (under reflection $x \rightarrow -x$ for a domain wall centered at the origin). We calculate the background gauge field profile $A_0$ associated with this domain wall profile. Interestingly, even with the asymmetric $l(x)$ profile, quark-antiquark scattering from the corresponding gauge field configuration does not reflect this asymmetry. We show that the expected asymmetry in scattering arises when we include the effect of asymmetric profile of $l(x)$ on the effective mass of quarks and antiquarks and calculate resultant scattering. We discuss the effects of such asymmetric Z(3) walls in generating quark and antiquark density fluctuations in cosmology, and in relativistic heavy-ion collisions e.g. event-by-event baryon fluctuations.

In the Randall-Sundrum model where the Standard Model fields are confined to the TeV brane located at the orbifold point $\theta = \pi$ and the gravity peaks at the Planck brane located at $\theta = 0$, the stabilized modulus (radion) field is required to stabilize the size of the fifth spatial dimension. It can be produced copiously inside the supernova core due to nucleon-nucleon bremstrahlung, electron-positron and plasmon-plasmon annihilations, which then subsequently decays to neutrino-antineutrino pair and take away the energy released in SN1987A explosion. Assuming that the supernovae cooling rate $\dot{\varepsilon} \le 7.288\times 10^{-27} \rm{GeV}$, we find the lower bound on the radion vev $\vphi \sim 9.0$ TeV, 2.2 TeV and 0.9 TeV corresponding to the radion mass $m_\phi = 5$ GeV, 20 GeV and 50 GeV, respectively.

We create dynamical models of the massive elliptical galaxy, NGC 4649, using the N-body made-to-measure code, NMAGIC, and kinematic constraints from long-slit and planetary nebula (PN) data. We explore a range of potentials based on previous determinations from X-ray observations and a dynamical model fitting globular cluster (GC) velocities and a stellar density profile. The X-ray mass distributions are similar in the central region but have varying outer slopes, while the GC mass profile is higher in the central region and on the upper end of the range further out. Our models cannot differentiate between the potentials in the central region, and therefore if non-thermal pressures or multi-phase components are present in the hot gas, they must be smaller than previously inferred. In the halo, we find that the PN velocities are sensitive tracers of the mass, preferring a less massive halo than that derived from the GC mass profile, but similar to one of the mass distributions derived from X-rays. Our results show that the GCs may form a dynamically distinct system, and that the properties of the hot gas derived from X-rays in the outer halo have considerable uncertainties that need to be better understood. Estimating the mass in stars using photometric information and a stellar population mass-to-light ratio, we infer a dark matter mass fraction in NGC 4649 of ~0.39 at 1Re (10.5 kpc) and ~0.78 at 4Re. We find that the stellar orbits are isotropic to mildly radial in the central ~6 kpc depending on the potential assumed. Further out, the orbital structure becomes slightly more radial along R and more isotropic along z, regardless of the potential assumed. In the equatorial plane, azimuthal velocity dispersions dominate over meridional velocity dispersions, implying that meridional velocity anisotropy is the mechanism for flattening the stellar system.

Using a slitless spectroscopy method with (a) the 8.2 m Subaru telescope and its FOCAS Cassegrain spectrograph, and (b) the ESO Very Large Telescope (VLT) unit 1 (Antu) and its FORS2 Cassegrain spectrograph, we have detected 326 planetary nebulae (PNs) in the giant Virgo elliptical galaxy NGC 4649 (M 60), and we have measured their radial velocities. After rejecting some PNs more likely to belong to the companion galaxy NGC 4647, we have built a catalog with kinematic information for 298 PNs in M 60. Using these radial velocities we have concluded that they support the presence of a dark matter halo around M 60. Based on an isotropic, two-component Hernquist model, we estimate the dark matter halo mass within 3$R_{\rm e}$ to be 4$\times10^{11} M_{\odot}$, which is almost one half of the total mass of about $10^{12} M_{\odot}$ within 3$R_{\rm e}$. This total mass is similar to that estimated from globular cluster, XMM-Newton and Chandra observations. The dark matter becomes dominant outside. More detailed dynamical modeling of the PN data is being published in a companion paper. We have also measured the $m$(5007) magnitudes of many of these PNs, and built a statistically complete sample of 218 PNs. The resulting PN luminosity function (PNLF) was used to estimate a distance modulus of 30.7$\pm$0.2 mag, equivalent to 14$\pm$1 Mpc. This confirms an earlier PNLF distance measurement, based on a much smaller sample. The PNLF distance modulus remains smaller than the surface brightness fluctuation (SBF) distance modulus by 0.4 mag.

To investigate the origins of S0 galaxies, we present a new method of analyzing their stellar kinematics from discrete tracers such as planetary nebulae. This method involves binning the data in the radial direction so as to extract the most general possible non-parametric kinematic profiles, and using a maximum likelihood fit within each bin in order to make full use of the information in the discrete kinematic tracers. Both disk and spheroid kinematic components are fitted, with a two-dimensional decomposition of imaging data used to attribute to each tracer a probability of membership in the separate components. Likelihood clipping also allows us to identify objects whose properties are not consistent with the adopted model, rendering the technique robust against contaminants and able to identify additional kinematic features. The method is first tested on an N-body simulated galaxy to assess possible sources of systematic error associated with the structural and kinematic decomposition, which are found to be small. It is then applied to the S0 system NGC~1023, for which a planetary nebula catalogue has already been released and analyzed by (Noordermeer et al., 2008). The correct inclusion of the spheroidal component allows us to show that, contrary to previous claims, the stellar kinematics of this galaxy are indistinguishable from those of a normal spiral galaxy, indicating that it may have evolved directly from such a system via gas stripping or secular evolution. The method also successfully identifies a population of outliers whose kinematics are different from those of the main galaxy; these objects can be identified with a stellar stream associated with the companion galaxy NGC~1023A.

As part of our current programme to test LCDM predictions for dark matter (DM) haloes using extended kinematical observations of early-type galaxies, we present a dynamical analysis of the bright elliptical galaxy NGC 4374 (M84) based on ~450 Planetary Nebulae (PNe) velocities from the PN.Spectrograph, along with extended long-slit stellar kinematics. This is the first such analysis of a galaxy from our survey with a radially constant velocity dispersion profile. We find that the spatial and kinematical distributions of the PNe agree with the field stars in the region of overlap. The velocity kurtosis is consistent with zero at almost all radii. We construct a series of Jeans models, fitting both velocity dispersion and kurtosis to help break the mass-anisotropy degeneracy. Our mass models include DM halos either with shallow cores or with central cusps as predicted by cosmological simulations - along with the novel introduction in this context of adiabatic halo contraction from baryon infall. Both classes of models confirm a very massive dark halo around NGC 4374, demonstrating that PN kinematics data are well able to detect such haloes when present. Considering the default cosmological mass model, we confirm earlier suggestions that bright galaxies tend to have halo concentrations higher than LCDM predictions, but this is found to be solved if either a Salpeter IMF or adiabatic contraction with a Kroupa IMF is assumed. Thus for the first time a case is found where the PN dynamics may well be consistent with a standard dark matter halo. A cored halo can also fit the data, and prefers a stellar mass consistent with a Salpeter IMF. The less dramatic dark matter content found in lower-luminosity "ordinary" ellipticals suggests a bimodality in the halo properties which may be produced by divergent baryonic effects during their assembly histories.

Using a counter-dispersed slitless spectroscopy technique, we detect and measure the line-of-sight velocities of 187 planetary nebulae (PNe) around one of the nearest cD galaxies, NGC 1399, with FORS1 on the VLT. We describe the method for identifying and classifying the emission-line sources and the procedure for computing their J2000 coordinates and velocities. The number of PN detections and the errors in the velocity measurements (37 km/s indicate that this technique is comparable to other methods, such as that described by Teodorescu et al. (2005). We present the spatial distribution of the PNe and a basic analysis of their velocities. The PN two-dimensional velocity field shows marginal rotation consistent with other studies. We also find a low-velocity substructure in the halo and a flatter velocity-dispersion profile compared to previous observations that extends to ~400 arcsec. The detection of a low-velocity subcomponent underscores the importance of discrete velocity tracers for the detection of un-mixed components. The new velocity-dispersion profile is in good agreement with revised velocity dispersions for the red globular clusters in NGC 1399, using the data of Schuberth et al. (2009). The outer parts of this profile are consistent with one of the dynamical models of Kronawitter et al. (2000), which corresponds to a circular velocity of ~340 km/s and a rescaled B-band mass-to-light ratio of ~20 at 7' radius. These measurements trace the kinematics of the outer halo and disentangle the heterogenous populations in the Fornax Cluster core. The new data set the stage for a revised dynamical model of the outer halo of NGC 1399.

We use a new non-parametric Bayesian approach to obtain the most probable mass distributions and circular velocity curves along with their confidence ranges, given deprojected density and temperature profiles of the hot gas surrounding X-ray bright elliptical galaxies. For a sample of six X-ray bright ellipticals, we find that all circular velocity curves are rising in the outer parts due to a combination of a rising temperature profile and a logarithmic pressure gradient that increases in magnitude. Comparing the circular velocity curves we obtain from X-rays to those obtained from dynamical models, we find that the former are often lower in the central ~10 kpc. This is probably due to a combination of: i) Non-thermal contributions of up to ~35% in the pressure (with stronger effects in NGC 4486), ii) multiple-temperature components in the hot gas, iii) incomplete kinematic spatial coverage in the dynamical models, and iv) mass profiles that are insufficiently general in the dynamical modelling. Complementing the total mass information from the X-rays with photometry and stellar population models to infer the dark matter content, we find evidence for massive dark matter haloes with dark matter mass fractions of ~35-80% at 2Re, rising to a maximum of 80-90% at the outermost radii. We also find that the six galaxies follow a Tully-Fisher relation with slope ~4 and that their circular velocities at 1Re correlate strongly with the velocity dispersion of the local environment. As a result, the galaxy luminosity at 1Re also correlates with the velocity dispersion of the environment. These relations suggest a close link between the properties of central X-ray bright elliptical galaxies and their environments (abridged).

We analyze six X-ray bright elliptical galaxies, observed with Chandra and XMM-Newton, and approximate their gravitational potentials by isothermal spheres phi(r)=v_c^2 ln(r) over a range of radii from ~0.5 to ~25 kpc. We then compare the circular speed v_c derived from X-ray data with the estimators available from optical data. In particular we discuss two simple and robust procedures for evaluating the circular speed of the galaxy using the observed optical surface brightness and the line-of-sight velocity dispersion profiles. The best fitting relation between the circular speeds derived from optical observations of stars and X-ray observations of hot gas is v_{c,opt}~ \eta * v_{c,X}, where \eta=1.10-1.15 (depending on the method), suggesting, albeit with large statistical and systematic uncertainties, that non-thermal pressure on average contributes ~20-30% of the gas thermal pressure.

We present high resolution FLAMES/VLT spectroscopy of intracluster planetary nebula (ICPN) candidates, targeting three new fields in the Virgo cluster core with surface brightness down to mu_B = 28.5. Based on the projected phase space information we separate the old and 12 newly-confirmed PNs into galaxy and intracluster components. The M87 PNs are confined to the extended stellar envelope of M87, within a projected radius of ~ 160 kpc, while the ICPNs are scattered across the whole surveyed region between M87 and M86. The velocity dispersions determined from the M87 PNs at projected radii of 60 kpc and 144 kpc show that the galaxy's velocity dispersion profile decreases in the outer halo, down to 78 +/- 25 km/s. A Jeans model for the M87 halo stars in the gravitational potential traced by the X-ray emission fits the observed velocity dispersion profile only if the stellar orbits are strongly radially anisotropic (beta ~= 0.4 at r ~= 10 kpc increasing to 0.8 at the outer edge), and if additionally the stellar halo is truncated at ~= 150 kpc average elliptical radius. From the spatial and velocity distribution of the ICPNs we infer that M87 and M86 are falling towards each other and that we may be observing them just before the first close pass. The inferred luminosity-specific PN numbers for the M87 halo and the ICL are in the range of values observed for old (> 10 Gyr) stellar populations (abridged).

We present new planetary nebulae (PNe) positions, radial velocities, and magnitudes for 6 early-type galaxies obtained with the Planetary Nebulae Spectrograph, their two-dimensional velocity and velocity dispersion fields. We extend this study to include an additional 10 early-type galaxies with PNe radial velocity measurements available from the literature, to obtain a broader description of the outer-halo kinematics in early-type galaxies. These data extend the information derived from stellar kinematics to typically up to ~8 Re. The combination of photometry, stellar and PNe kinematics shows: i) good agreement between the PNe number density and the stellar surface brightness in the region where the two data sets overlap; ii) good agreement between PNe and stellar kinematics; iii) that the mean rms velocity profiles fall into two groups: with of the galaxies characterized by slowly decreasing profiles and the remainder having steeply falling profiles; iv) a larger variety of velocity dispersion profiles; v) that twists and misalignments in the velocity fields are more frequent at large radii, including some fast rotators; vi) that outer haloes are characterised by more complex radial profiles of the specific angular momentum-related lambda_R parameter than observed within 1Re; vii) that many objects are more rotationally dominated at large radii than in their central parts; and viii) that the halo kinematics are correlated with other galaxy properties, such as total luminosity, isophotal shape, total stellar mass, V/sigma, and alpha parameter, with a clear separation between fast and slow rotators.

We present first results of a study of the halo kinematics for a sample of early type galaxies using planetary nebulae (PNe) as kinematical tracers. PNe allow to extend up to several effective radii (Re) the information from absorption line kinematics (confined to within 1 or 2 Re), providing valuable information and constraints for merger simulations and galaxy formation models. We find that the specific angular momentum per unit mass has a more complex radial dependence when the halo region is taken into account and that the halo velocity dispersion is related to the total galaxy luminosity, isophotal shape, and number of PNe per unit of luminosity

We present new Planetary Nebula Spectrograph observations of the ordinary elliptical galaxy NGC 4494, resulting in positions and velocities of 255 PNe out to 7 effective radii (25 kpc). We also present new wide-field surface photometry from MMT/Megacam, and long-slit stellar kinematics from VLT/FORS2. The spatial and kinematical distributions of the PNe agree with the field stars in the region of overlap. The mean rotation is relatively low, with a possible kinematic axis twist outside 1 Re. The velocity dispersion profile declines with radius, though not very steeply, down to ~70 km/s at the last data point. We have constructed spherical dynamical models of the system, including Jeans analyses with multi-component LCDM-motivated galaxies as well as logarithmic potentials. These models include special attention to orbital anisotropy, which we constrain using fourth-order velocity moments. Given several different sets of modelling methods and assumptions, we find consistent results for the mass profile within the radial range constrained by the data. Some dark matter (DM) is required by the data; our best-fit solution has a radially anisotropic stellar halo, a plausible stellar mass-to-light ratio, and a DM halo with an unexpectedly low central density. We find that this result does not substantially change with a flattened axisymmetric model. Taken together with other results for galaxy halo masses, we find suggestions for a puzzling pattern wherein most intermediate-luminosity galaxies have very low concentration halos, while some high-mass ellipticals have very high concentrations. We discuss some possible implications of these results for DM and galaxy formation.

We examine the dark matter properties of nearby early-type galaxies using planetary nebulae (PNe) as mass probes. We have designed a specialised instrument, the Planetary Nebula Spectrograph (PN.S) operating at the William Herschel telescope, with the purpose of measuring PN velocities with best efficiency. The primary scientific objective of this custom-built instrument is the study of the PN kinematics in 12 ordinary round galaxies. Preliminary results showing a dearth of dark matter in ordinary galaxies (Romanowsky et al. 2003) are now confirmed by the first complete PN.S datasets. On the other hand early-type galaxies with a "regular" dark matter content are starting to be observed among the brighter PN.S target sample, thus confirming a correlation between the global dark-to-luminous mass virial ratio (f_DM=M_DM/M_star) and the galaxy luminosity and mass.

Individual pulses from pulsars have intensity-phase profiles that differ widely from pulse to pulse, from the average profile, and from phase to phase within a pulse. Widely accepted explanations do not exist for this variability or for the mechanism producing the radiation. The variability corresponds to the field statistics, particularly the distribution of wave field amplitudes, which are predicted by theories for wave growth in inhomogeneous media. This paper shows that the field statistics of the Vela pulsar (PSR B0833-45) are well-defined and vary as a function of pulse phase, evolving from Gaussian intensity statistics off-pulse to approximately power-law and then lognormal distributions near the pulse peak to approximately power-law and eventually Gaussian statistics off-pulse again. Detailed single-component fits confirm that the variability corresponds to lognormal statistics near the peak of the pulse profile and Gaussian intensity statistics off-pulse. The lognormal field statistics observed are consistent with the prediction of stochastic growth theory (SGT) for a purely linear system close to marginal stability. The simplest interpretations are that the pulsar's variability is a direct manifestation of an SGT state and the emission mechanism is linear (either direct or indirect), with no evidence for nonlinear mechanisms like modulational instability and wave collapse which produce power-law field statistics. Stringent constraints are placed on nonlinear mechanisms: they must produce lognormal statistics when suitably ensemble-averaged. Field statistics are thus a powerful, potentially widely applicable tool for understanding variability and constraining mechanisms and source characteristics of coherent astrophysical and space emissions.

The variability of the Vela pulsar (PSR B0833-45) corresponds to well-defined field statistics that vary with pulsar phase, ranging from Gaussian intensity statistics off-pulse to approximately power-law statistics in a transition region and then lognormal statistics on-pulse, excluding giant micropulses. These data are analyzed here in terms of two superposed wave populations, using a new calculation for the amplitude statistics of two vectorially-combined transverse fields. Detailed analyses show that the approximately power-law and lognormal distributions observed are fitted well at essentially all on-pulse phases by Gaussian-lognormal and double-lognormal combinations, respectively. These good fits, plus the smooth but significant variations in fit parameters across the source, provide strong evidence that the approximately power-law statistics observed in the transition region are not intrinsic. Instead, the data are consistent with normal pulsar emission having lognormal statistics at all phases. This is consistent with generation in an inhomogeneous source obeying stochastic growth theory (SGT) and with the emission mechanism being purely linear (either direct or indirect). A nonlinear mechanism is viable only if it produces lognormal statistics when suitably ensemble-averaged. Variations in the SGT fit parameters with phase imply that the radiation is relatively more variable near the pulse edges than near the center, as found in earlier work. In contrast, Vela's giant micropulses come from a very restricted phase range and have power-law statistics with indices ($6.7 \pm 0.6$) not inconsistent with nonlinear wave collapse. These results imply that normal pulses have a different source and generation mechanism than giant micropulses, as suggested previously on other grounds.

Individual pulses from pulsars have intensity-phase profiles that differ widely from pulse to pulse, from the average profile, and from phase to phase within a pulse. Widely accepted explanations for pulsar radio emission and its time variability do not exist. Here, by analysing data near the peak of the Vela pulsar's average profile, we show that Vela's variability corresponds to lognormal field statistics, consistent with the prediction of stochastic growth theory (SGT) for a purely linear system close to marginal stability. Vela's variability is therefore a direct manifestation of an SGT state and the field statistics constrain the emission mechanism to be linear (either direct or indirect), ruling out nonlinear mechanisms like wave collapse. Field statistics are thus a powerful, potentially widely applicable tool for understanding variability and constraining mechanisms and source characteristics of coherent astrophysical and space emissions.

Low energy string theory predicts the existence of an axion field which can lead to cosmic birefringence. We solve the electromagnetic wave equations in the presence of such an axion and a dilaton field in order to determine their effect on the polarization of light. We find that the presence of dilaton field leads to a nontrivial modification of the final result. We comment on the possibility of discovering such an effect by observations of radio wave polarizations from distant radio galaxies and quasars. We have also determined the limits on the string theory parameters that are imposed by the current polarization data from distant quasars and radio galaxies.