We build a theoretical picture of how the light from galaxies evolves across cosmic time. In particular, we predict the evolution of the galaxy spectral energy distribution (SED) by carefully integrating the star formation and metal enrichment histories of semi-analytic model (SAM) galaxies and combining these with stellar population synthesis models which we call mentari. Our SAM combines prescriptions to model the interplay between gas accretion, star formation, feedback process, and chemical enrichment in galaxy evolution. From this, the SED of any simulated galaxy at any point in its history can be constructed and compared with telescope data to reverse engineer the various physical processes that may have led to a particular set of observations. The synthetic SEDs of millions of simulated galaxies from mentari can cover wavelengths from the far UV to infrared, and thus can tell a near complete story of the history of galaxy evolution. \keywords{galaxies: evolution - galaxies: stellar content - galaxies.}

We study the dust evolution in galaxies by implementing a detailed dust prescription in the SAGE semi-analytical model for galaxy formation. The new model, called Dusty SAGE, follows the condensation of dust in the ejecta of type II supernovae and asymptotic giant branch (AGB) stars, grain growth in the dense molecular clouds, destruction by supernovae shocks, and the removal of dust from the ISM by star formation, reheating, inflows and outflows. Our model successfully reproduces the observed dust mass function at redshift z = 0 and the observed scaling relations for dust across a wide range of redshifts. We find that the dust mass content in the present Universe is mainly produced via grain growth in the interstellar medium (ISM). By contrast, in the early Universe, the primary production mechanism for dust is the condensation in stellar ejecta. The shift of the significant production channel for dust characterises the scaling relations of dust-to-gas (DTG) and dust-to-metal (DTM) ratios. In galaxies where the grain growth dominates, we find positive correlations for DTG and DTM ratios with both metallicity and stellar mass. On the other hand, in galaxies where dust is produced primarily via condensation, we find negative or no correlation for DTM and DTG ratios with either metallicity or stellar mass. In agreement with observation showing that the circumgalactic medium (CGM) contains more dust than the ISM, our model also shows the same trend for z < 4. Our semi-analytic model is publicly available at https: //github.com/dptriani/dusty-sage.

We compute the bispectra of the 21cm signal during the Epoch of Reionization for three different reionization scenarios that are based on a dark matter N-body simulation combined with a self-consistent, semi-numerical model of galaxy evolution and reionization. Our reionization scenarios differ in their trends of ionizing escape fractions ($f_\mathrm{esc}$) with the underlying galaxy properties and cover the physically plausible range, i.e. $f_\mathrm{esc}$ effectively decreasing, being constant, or increasing with halo mass. We find the 21cm bispectrum to be sensitive to the resulting ionization topologies that significantly differ in their size distribution of ionized and neutral regions throughout reionization. From squeezed to stretched triangles, the 21cm bispectra features a change of sign from negative to positive values, with ionized and neutral regions representing below-average and above-average concentrations contributing negatively and positively, respectively. The position of the change of sign provides a tracer of the size distribution of the ionized and neutral regions, and allows us to identify three major regimes that the 21cm bispectrum undergoes during reionization. In particular the regime during the early stages of reionization, where the 21cm bispectrum tracks the peak of the size distribution of the ionized regions, provides exciting prospects for pinning down reionization with the forthcoming Square Kilometre Array.

emcee is a Python library implementing a class of affine-invariant ensemble samplers for Markov chain Monte Carlo (MCMC). This package has been widely applied to probabilistic modeling problems in astrophysics where it was originally published, with some applications in other fields. When it was first released in 2012, the interface implemented in emcee was fundamentally different from the MCMC libraries that were popular at the time, such as PyMC, because it was specifically designed to work with "black box" models instead of structured graphical models. This has been a popular interface for applications in astrophysics because it is often non-trivial to implement realistic physics within the modeling frameworks required by other libraries. Since emcee's release, other libraries have been developed with similar interfaces, such as dynesty (Speagle 2019). The version 3.0 release of emcee is the first major release of the library in about 6 years and it includes a full re-write of the computational backend, several commonly requested features, and a set of new "move" implementations.

Correlation functions are widely used in extra-galactic astrophysics to extract insights into how galaxies occupy dark matter halos and in cosmology to place stringent constraints on cosmological parameters. A correlation function fundamentally requires computing pair-wise separations between two sets of points and then computing a histogram of the separations. Corrfunc is an existing open-source, high-performance software package for efficiently computing a multitude of correlation functions. In this paper, we will discuss the SIMD AVX512F kernels within Corrfunc, capable of processing 16 floats or 8 doubles at a time. The latest manually implemented Corrfunc AVX512F kernels show a speedup of up to $\sim 4\times$ relative to compiler-generated code for double-precision calculations. The AVX512F kernels show $\sim 1.6\times$ speedup relative to the AVX kernels and compare favorably to a theoretical maximum of $2\times$. In addition, by pruning pairs with too large of a minimum possible separation, we achieve a $\sim 5-10\%$ speedup across all the SIMD kernels. Such speedups highlight the importance of programming explicitly with SIMD vector intrinsics for complex calculations that can not be efficiently vectorized by compilers. Corrfunc is publicly available at https://github.com/manodeep/Corrfunc/.

The two-point correlation function (2PCF) is the most widely used tool for quantifying the spatial distribution of galaxies. Since the distribution of galaxies is determined by galaxy formation physics as well as the underlying cosmology, fitting an observed correlation function yields valuable insights into both. The calculation for a 2PCF involves computing pair-wise separations and consequently, the computing time scales quadratically with the number of galaxies. The next-generation galaxy surveys are slated to observe many millions of galaxies, and computing the 2PCF for such surveys would be prohibitively time-consuming. Additionally, modern modelling techniques require the 2PCF to be calculated thousands of times on simulated galaxy catalogues of {\em at least} equal size to the data and would be completely unfeasible for the next generation surveys. Thus, calculating the 2PCF forms a substantial bottleneck in improving our understanding of the fundamental physics of the universe, and we need high-performance software to compute the correlation function. In this paper, we present Corrfunc --- a suite of highly optimised, OpenMP parallel clustering codes. The improved performance of Corrfunc arises from both efficient algorithms as well as software design that suits the underlying hardware of modern CPUs. Corrfunc can compute a wide range of 2-D and 3-D correlation functions in either simulation (Cartesian) space or on-sky coordinates. Corrfunc runs efficiently in both single- and multi-threaded modes and can compute a typical 2-point projected correlation function ($w_p(r_p)$) for ~1 million galaxies within a few seconds on a single thread. Corrfunc is designed to be both user-friendly and fast and is publicly available at https://github.com/manodeep/Corrfunc.

Fundamental coding and software development skills are increasingly necessary for success in nearly every aspect of astronomical and astrophysical research as large surveys and high resolution simulations become the norm. However, professional training in these skills is inaccessible or impractical for many members of our community. Students and professionals alike have been expected to acquire these skills on their own, apart from formal classroom curriculum or on-the-job training. Despite the recognized importance of these skills, there is little opportunity to develop them - even for interested researchers. To ensure a workforce capable of taking advantage of the computational resources and the large volumes of data coming in the next decade, we must identify and support ways to make software development training widely accessible to community members, regardless of affiliation or career level. To develop and sustain a technology capable astronomical and astrophysical workforce, we recommend that agencies make funding and other resources available in order to encourage, support and, in some cases, require progress on necessary training, infrastructure and policies. In this white paper, we focus on recommendations for how funding agencies can lead in the promotion of activities to support the astronomy and astrophysical workforce in the 2020s.

One of the most important parameters in characterizing the Epoch of Reionization, the escape fraction of ionizing photons, $f_\mathrm{esc}$, remains unconstrained both observationally and theoretically. With recent work highlighting the impact of galaxy-scale feedback on the instantaneous value of $f_\mathrm{esc}$, it is important to develop a model in which reionization is self-consistently coupled to galaxy evolution. In this work, we present such a model and explore how physically motivated functional forms of $f_\mathrm{esc}$ affect the evolution of ionized hydrogen within the intergalactic medium. Using the $21$cm power spectrum evolution, we investigate the likelihood of observationally distinguishing between a constant $f_\mathrm{esc}$ and other models that depend upon different forms of galaxy feedback. We find that changing the underlying connection between $f_\mathrm{esc}$ and galaxy feedback drastically alters the large-scale $21$cm power. The upcoming Square Kilometre Array Low Frequency instrument possesses the sensitivity to differentiate between our models at a fixed optical depth, requiring only $200$ hours of integration time focused on redshifts $z = 7.5-8.5$. Generalizing these results to account for a varying optical depth will require multiple $800$ hour observations spanning redshifts $z = 7-10$. This presents an exciting opportunity to observationally constrain one of the most elusive parameters during the Epoch of Reionization.

Software is a critical part of modern research, and yet there are insufficient mechanisms in the scholarly ecosystem to acknowledge, cite, and measure the impact of research software. The majority of academic fields rely on a one-dimensional credit model whereby academic articles (and their associated citations) are the dominant factor in the success of a researcher's career. In the petabyte era of astronomical science, citing software and measuring its impact enables academia to retain and reward researchers that make significant software contributions. These highly skilled researchers must be retained to maximize the scientific return from petabyte-scale datasets. Evolving beyond the one-dimensional credit model requires overcoming several key challenges, including the current scholarly ecosystem and scientific culture issues. This white paper will present these challenges and suggest practical solutions for elevating the role of software as a product of the research enterprise.

We have built PRISM, a "Probabilistic Regression Instrument for Simulating Models". PRISM uses the Bayes linear approach and history matching to construct an approximation ('emulator') of any given model, by combining limited model evaluations with advanced regression techniques, covariances and probability calculations. It is designed to easily facilitate and enhance existing Markov chain Monte Carlo (MCMC) methods by restricting plausible regions and exploring parameter space efficiently. However, PRISM can additionally be used as a standalone alternative to MCMC for model analysis, providing insight into the behavior of complex scientific models. With PRISM, the time spent on evaluating a model is minimized, providing developers with an advanced model analysis for a fraction of the time required by more traditional methods. This paper provides an overview of the different techniques and algorithms that are used within PRISM. We demonstrate the advantage of using the Bayes linear approach over a full Bayesian analysis when analyzing complex models. Our results show how much information can be captured by PRISM and how one can combine it with MCMC methods to significantly speed up calibration processes (>15 times faster). PRISM is an open-source Python package that is available under the BSD 3-Clause License (BSD-3) at https://github.com/1313e/PRISM and hosted at https://prism-tool.readthedocs.io. PRISM has also been reviewed by "The Journal of Open Source Software" (https://doi.org/10.21105/joss.01229).

We investigate the role of angular momentum in the clustering of dark matter haloes. We make use of data from two high-resolution N-body simulations spanning over four orders of magnitude in halo mass, from $10^{9.8}$ to $10^{14}\ h^{-1}\ \text{M}_\odot$. We explore the hypothesis that mass accretion in filamentary environments alters the angular momentum of a halo, thereby driving a correlation between the spin parameter $\lambda$ and the strength of clustering. However, we do not find evidence that the distribution of matter on large scales is related to the spin of haloes. We find that a halo's spin is correlated with its age, concentration, sphericity, and mass accretion rate. Removing these correlations strongly affects the strength of secondary spin bias at low halo masses. We also find that high spin haloes are slightly more likely to be found near another halo of comparable mass. These haloes that are found near a comparable mass neighbour - a \textit{twin} - are strongly spatially biased. We demonstrate that this \textit{twin bias}, along with the relationship between spin and mass accretion rates, statistically accounts for halo spin secondary bias.

We introduce the THE THREE HUNDRED project, an endeavour to model 324 large galaxy clusters with full-physics hydrodynamical re-simulations. Here we present the data set and study the differences to observations for fundamental galaxy cluster properties and scaling relations. We find that the modelled galaxy clusters are generally in reasonable agreement with observations with respect to baryonic fractions and gas scaling relations at redshift z = 0. However, there are still some (model-dependent) differences, such as central galaxies being too massive, and galaxy colours (g - r) being bluer (about 0.2 dex lower at the peak position) than in observations. The agreement in gas scaling relations down to 10^{13} h^{-1} M_{\odot} between the simulations indicates that particulars of the sub-grid modelling of the baryonic physics only has a weak influence on these relations. We also include - where appropriate - a comparison to three semi-analytical galaxy formation models as applied to the same underlying dark-matter-only simulation. All simulations and derived data products are publicly available.

We explore the connection between the atomic gas fraction, f_atm, and `global disc stability' parameter, q, of galaxies within a fully cosmological context by examining galaxies in the Dark Sage semi-analytic model. The q parameter is determined by the ratio of disc specific angular momentum to mass. Dark Sage is well suited to our study, as it includes the numerical evolution of one-dimensional disc structure, making both j_disc and q predicted quantities. We show that Dark Sage produces a clear correlation between gas fraction and j_disc at fixed disc mass, in line with recent results from observations and hydrodynamic simulations. This translates to a tight q--f_atm sequence for star-forming central galaxies, which closely tracks the analytic prediction of Obreschkow et al. The scatter in this sequence is driven by the probability distribution function of mass as a function of j (PDF of j) within discs, specifically where it peaks. We find that halo mass is primarily responsible for the peak location of the PDF of j, at least for low values of q. Two main mechanisms of equal significance are then identified for disconnecting f_atm from q. Mergers in the model can trigger quasar winds, with the potential to blow out most of the gas disc, while leaving the stellar disc relatively unharmed. Ram-pressure stripping of satellite galaxies has a similar effect, where f_atm can drop drastically with only a minimal effect to q. We highlight challenges associated with following these predictions up with observations.

Galactic conformity is the phenomenon whereby galaxy properties exhibit excess correlations across distance than that expected if these properties only depended on halo mass. We perform a comprehensive study of conformity at low redshift using a galaxy group catalogue from the SDSS DR7 spectroscopic sample. We study correlations both between central galaxies and their satellites (1-halo), and between central galaxies in separate haloes (2-halo). We use the quenched fractions and the marked correlation function (MCF), to probe for conformity in three galaxy properties, $(g-r)$ colour, specific star formation rate (sSFR), and morphology. We assess the statistical significance of conformity signals with a suite of mock galaxy catalogues that have no built-in conformity, but contain the same group-finding and mass assignment errors as the real data. In the case of 1-halo conformity, quenched fractions show strong signals at all group masses. However, these signals are equally strong in mock catalogues, indicating that the conformity signal is spurious and likely entirely caused by group-finding systematics, calling into question previous claims of 1-halo conformity detection. The MCF reveals a significant detection of radial segregation within massive groups, but no evidence of conformity. In the case of 2-halo conformity, quenched fractions show no significant evidence of conformity in colour or sSFR once compared with mock catalogues, but a clear signal using morphology. In contrast, the MCF reveals a small, yet highly significant signal for all three properties in low mass groups and scales of $0.8-4\ h^{-1}\textrm{Mpc}$, possibly representing the first robust detection of 2-halo conformity.

Interpreting the small-scale clustering of galaxies with halo models can elucidate the connection between galaxies and dark matter halos. Unfortunately, the modelling is typically not sufficiently accurate for ruling out models statistically. It is thus difficult to use the information encoded in small scales to test cosmological models or probe subtle features of the galaxy-halo connection. In this paper, we attempt to push halo modelling into the "accurate" regime with a fully numerical mock-based methodology and careful treatment of statistical and systematic errors. With our forward-modelling approach, we can incorporate clustering statistics beyond the traditional two-point statistics. We use this modelling methodology to test the standard $\Lambda\mathrm{CDM}$ + halo model against the clustering of SDSS DR7 galaxies. Specifically, we use the projected correlation function, group multiplicity function and galaxy number density as constraints. We find that while the model fits each statistic separately, it struggles to fit them simultaneously. Adding group statistics leads to a more stringent test of the model and significantly tighter constraints on model parameters. We explore the impact of varying the adopted halo definition and cosmological model and find that changing the cosmology makes a significant difference. The most successful model we tried (Planck cosmology with Mvir halos) matches the clustering of low luminosity galaxies, but exhibits a 2.3$\sigma$ tension with the clustering of luminous galaxies, thus providing evidence that the "standard" halo model needs to be extended. This work opens the door to adding interesting freedom to the halo model and including additional clustering statistics as constraints.

The exact role of quasars during the Epoch of Reionization remains uncertain. With consensus leaning towards quasars producing a negligible amount of ionizing photons, we pose an alternate question: Can quasars indirectly contribute to reionization by allowing ionizing photons from stars to escape more easily? Using the Semi-Analytic Galaxy Evolution model to evolve a galaxy population through cosmic time, we construct an idealized scenario in which the escape fraction of stellar ionizing photons ($f_\mathrm{esc}$) is boosted following quasar wind events, potentially for several dynamical times. We find that under this scenario, the mean value of $f_\mathrm{esc}$ as a function of galaxy stellar mass peaks for intermediate mass galaxies. This mass dependence will have consequences for the 21cm power spectrum, enhancing power at small scales and suppressing it at large scales. This hints that whilst quasars may not directly contribute to the ionizing photon budget, they could influence reionization indirectly by altering the topology of ionized regions.

Standard present day large-scale structure (LSS) analyses make a major assumption in their Bayesian parameter inference --- that the likelihood has a Gaussian form. For summary statistics currently used in LSS, this assumption, even if the underlying density field is Gaussian, cannot be correct in detail. We investigate the impact of this assumption on two recent LSS analyses: the Beutler et al. (2017) power spectrum multipole ($P_\ell$) analysis and the Sinha et al. (2017) group multiplicity function ($\zeta$) analysis. Using non-parametric divergence estimators on mock catalogs originally constructed for covariance matrix estimation, we identify significant non-Gaussianity in both the $P_\ell$ and $\zeta$ likelihoods. We then use Gaussian mixture density estimation and Independent Component Analysis on the same mocks to construct likelihood estimates that approximate the true likelihood better than the Gaussian $pseudo$-likelihood. Using these likelihood estimates, we accurately estimate the true posterior probability distribution of the Beutler et al. (2017) and Sinha et al. (2017) parameters. Likelihood non-Gaussianity shifts the $f\sigma_8$ constraint by $-0.44\sigma$, but otherwise, does not significantly impact the overall parameter constraints of Beutler et al. (2017). For the $\zeta$ analysis, using the pseudo-likelihood significantly underestimates the uncertainties and biases the constraints of Sinha et al. (2017) halo occupation parameters. For $\log M_1$ and $\alpha$, the posteriors are shifted by $+0.43\sigma$ and $-0.51\sigma$ and broadened by $42\%$ and $66\%$, respectively. The divergence and likelihood estimation methods we present provide a straightforward framework for quantifying the impact of likelihood non-Gaussianity and deriving more accurate parameter constraints.

We explore the phenomenon commonly known as halo assembly bias, whereby dark matter halos of the same mass are found to be more or less clustered when a second halo property is considered, for halos in the mass range $3.7 \times 10^{11} \; h^{-1} \mathrm{M_{\odot}} - 5.0 \times 10^{13} \; h^{-1} \mathrm{M_{\odot}}$. Using the Large Suite of Dark Matter Simulations (LasDamas) we consider nine commonly used halo properties and find that a clustering bias exists if halos are binned by mass or by any other halo property. This secondary bias implies that no single halo property encompasses all the spatial clustering information of the halo population. The mean values of some halo properties depend on their halo's distance to a more massive neighbor. Halo samples selected by having high values of one of these properties therefore inherit a neighbor bias such that they are much more likely to be close to a much more massive neighbor. This neighbor bias largely accounts for the secondary bias seen in halos binned by mass and split by concentration or age. However, halos binned by other mass-like properties still show a secondary bias even when the neighbor bias is removed. The secondary bias of halos selected by their spin behaves differently than that for other halo properties, suggesting that the origin of the spin bias is different than of other secondary biases.

We present the public release of the MultiDark-Galaxies: three distinct galaxy catalogues derived from one of the Planck cosmology MultiDark simulations (i.e. MDPL2, with a volume of (1 Gpc/$h$)$^{3}$ and mass resolution of $1.5 \times 10^{9} M_{\odot}/h$) by applying the semi-analytic models GALACTICUS, SAG, and SAGE to it. We compare the three models and their conformity with observational data for a selection of fundamental properties of galaxies like stellar mass function, star formation rate, cold gas fractions, and metallicities - noting that they sometimes perform differently reflecting model designs and calibrations. We have further selected galaxy subsamples of the catalogues by number densities in stellar mass, cold gas mass, and star formation rate in order to study the clustering statistics of galaxies. We show that despite different treatment of orphan galaxies, i.e. galaxies that lost their dark-matter host halo due to the finite mass resolution of the N-body simulation or tidal stripping, the clustering signal is comparable, and reproduces the observations in all three models - in particular when selecting samples based upon stellar mass. Our catalogues provide a powerful tool to study galaxy formation within a volume comparable to those probed by on-going and future photometric and redshift surveys. All model data consisting of a range of galaxy properties - including broad-band SDSS magnitudes - are publicly available.

We present the first stable release of Halotools (v0.2), a community-driven Python package designed to build and test models of the galaxy-halo connection. Halotools provides a modular platform for creating mock universes of galaxies starting from a catalog of dark matter halos obtained from a cosmological simulation. The package supports many of the common forms used to describe galaxy-halo models: the halo occupation distribution (HOD), the conditional luminosity function (CLF), abundance matching, and alternatives to these models that include effects such as environmental quenching or variable galaxy assembly bias. Satellite galaxies can be modeled to live in subhalos, or to follow custom number density profiles within their halos, including spatial and/or velocity bias with respect to the dark matter profile. The package has an optimized toolkit to make mock observations on a synthetic galaxy population, including galaxy clustering, galaxy-galaxy lensing, galaxy group identification, RSD multipoles, void statistics, pairwise velocities and others, allowing direct comparison to observations. Halotools is object-oriented, enabling complex models to be built from a set of simple, interchangeable components, including those of your own creation. Halotools has an automated testing suite and is exhaustively documented on http://halotools.readthedocs.io, which includes quickstart guides, source code notes and a large collection of tutorials. The documentation is effectively an online textbook on how to build and study empirical models of galaxy formation with Python.

MeV-scale magnetic fields in the interiors of magnetars suppress the pairing of neutrons and protons in the $S$-wave state. In the case of a neutron condensate the suppression is the consequence of the Pauli-paramagnetism of the neutron gas, i.e., the alignment of the neutron spins along the magnetic field. The proton $S$-wave pairing is suppressed because of the Landau diamagnetic currents of protons induced by the field. The Ginzburg-Landau and BCS theories of the critical magnetic fields for unpairing are reviewed. The macrophysical implications of the suppression (unpairing) of the condensates are discussed for the rotational crust-core coupling in magnetars and the neutrino-dominated cooling era of their thermal evolution.

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.

Recent observations of AGN activity in massive galaxies (log Mstar / Msun > 10.4) show that: 1) at z < 1, AGN-hosting galaxies do not show enhanced merger signatures compared to normal galaxies, 2) also at z < 1, most AGNs are hosted by quiescent galaxies; and 3) at z > 1, percentage of AGNs in star forming galaxies increases and becomes comparable to AGN percentage in quiescent galaxies at z ~ 2. How can major mergers explain AGN activity in massive quiescent galaxies which have no merger features and no star formation to indicate recent galaxy merger? By matching merger events in a cosmological N-body simulation to the observed AGN incidence probability in the COSMOS survey, we show that major merger triggered AGN activity is consistent with the observations. By distinguishing between "peak" AGNs (recently merger triggered and hosted by star forming galaxies) and "faded" AGNs (merger triggered a long time ago and now residing in quiescent galaxies), we show that the AGN occupation fraction in star forming and quiescent galaxies simply follows the evolution of the galaxy merger rate. Since the galaxy merger rate drops dramatically at z < 1, the only AGNs left to be observed are the ones triggered by old mergers and are now in the declining phase of their nuclear activity, hosted by quiescent galaxies. As we go toward higher redshifts the galaxy merger rate increases and the percentages of "peak" AGNs and "faded" AGNs become comparable.

We use subhalo abundance matching (SHAM) to model the stellar mass function (SMF) and clustering of the Baryon Oscillation Spectroscopic Survey (BOSS) "CMASS" sample at $z\sim0.5$. We introduce a novel method which accounts for the stellar mass incompleteness of CMASS as a function of redshift, and produce CMASS mock catalogs which include selection effects, reproduce the overall SMF, the projected two-point correlation function $w_{\rm p}$, the CMASS $dn/dz$, and are made publicly available. We study the effects of assembly bias above collapse mass in the context of "age matching" and show that these effects are markedly different compared to the ones explored by Hearin et al. (2013) at lower stellar masses. We construct two models, one in which galaxy color is stochastic ("AbM" model) as well as a model which contains assembly bias effects ("AgM" model). By confronting the redshift dependent clustering of CMASS with the predictions from our model, we argue that that galaxy colors are not a stochastic process in high-mass halos. Our results suggest that the colors of galaxies in high-mass halos are determined by other halo properties besides halo peak velocity and that assembly bias effects play an important role in determining the clustering properties of this sample.

In the second paper of this series, we present results from cosmological simulations on the demographics of flyby interactions to gauge their potential impact on galaxy evolution. In a previous paper, we demonstrated that flybys -- an interaction where two independent halos inter-penetrate but detach at a later time and do not merge -- occur much more frequently than previously believed. In particular, we found that the frequency of flybys increases at low redshift and is comparable to or even greater than the frequency of mergers for halos $\gtrsim 10^{11} M_\odot/h$. In this paper, we classify flybys according to their orbits and the level of perturbation exacted on both the halos involved. We find that the majority of flybys penetrate deeper than $\sim R_{half}$ of the primary and have an initial relative speed $\sim 1.6\times V_{vir}$ of the primary. The typical flyby mass-ratio is $\sim 0.1$ at high $z$ for all halos, while at low $z$, massive primary halos undergo flybys with small secondary halos. We estimate the perturbation from the flyby on both the primary and the secondary and find that a typical flyby is mostly non-perturbative for the primary halo. However, since a massive primary experiences so many flybys at any given time, they are nearly continually a victim of a perturbative event. In particular, we find flybys that cause $\sim 1\%$ change in the binding energy of a primary halo occurs $\gtrsim 1 $ Gyr$^{-1}$ for halos $> 10^{10} M_\odot/h$ for $z \lesssim 4$. Secondary halos, on the other hand, are highly perturbed by the typical encounter, experiencing a change in binding energy of nearly order unity. Our results imply that flybys can drive a significant part of galaxy transformation at moderate to lower redshifts ($z \lesssim 4$). We touch on implications for observational surveys, mass-to-light ratios, and galaxy assembly bias.

We present and test TesseRACt, a non-parametric technique for recovering the concentration of simulated dark matter halos using Voronoi tessellation. TesseRACt is tested on idealized N-body halos that are axisymmetric, triaxial, and contain substructure and compared to traditional least-squares fitting as well as two non-parametric techniques that assume spherical symmetry. TesseRACt recovers halo concentrations within 0.3% of the true value regardless of whether the halo is spherical, axisymmetric, or triaxial. Traditional fitting and non-parametric techniques that assume spherical symmetry can return concentrations that are systematically off by as much as 10% from the true value for non-spherical halos. TesseRACt also performs significantly better when there is substructure present outside $0.5R_{200}$. Given that cosmological halos are rarely spherical and often contain substructure, we discuss implications for studies of halo concentration in cosmological N-body simulations including how choice of technique for measuring concentration might bias scaling relations.

We describe the microphysics, phenomenology, and astrophysical implication of a $B$-field induced unpairing effect that may occur in magnetars, if the local $B$-field in the core of a magnetar exceeds a critical value $H_{c2}$. Using the Ginzburg-Landau theory of superconductivity, we derive the $H_{c2}$ field for proton condensate taking into the correction ($\le 30\%$) which arises from its coupling to the background neutron condensate. The density dependence of pairing of proton condensate implies that $H_{c2}$ is maximal at the crust-core interface and decreases towards the center of the star. As a consequence, magnetar cores with homogenous constant fields will be partially superconducting for "medium-field" magnetars ($10^{15}\le B\le 5 \times 10^{16}$ G) whereas "strong-field" magnetars ($B>5\times 10^{16}$ G) will be void of superconductivity. The neutrino emissivity of a magnetar's core changes in a twofold manner: (i)~the $B$-field assisted direct Urca process is enhanced by orders of magnitude, because of the unpairing effect in regions where $B\ge H_{c2}$; (ii)~the Cooper-pair breaking processes on protons vanish in these regions and the overall emissivity by the pair-breaking processes is reduced by a factor of only a few.

Compact stars having strong magnetic fields (magnetars) have been observationally determined to have surface magnetic fields of order of 10^14-10^15 G, the implied internal field strength being several orders larger. We study the equation of state and composition of hypernuclear matter and quark matter - two forms of dense matter in strong magnetic fields. We find that the magnetic field has substantial influence on the properties of hypernuclear matter and quark matter for magnetic field B \ge 10^17 G and B \ge 10^18 G respectively. In particular the matter properties become anisotropic. Moreover, above a critical field B_cr, both hypernuclear and quark matter show instability, although the values of B_cr are different for two kinds of matter.

We study the structure and evolution of dark matter halos from z = 300 to z = 6 for two cosmological N-body simulation initialization techniques. While the second order Lagrangian perturbation theory (2LPT) and the Zel'dovich approximation (ZA) both produce accurate present day halo mass functions, earlier collapse of dense regions in 2LPT can result in larger mass halos at high redshift. We explore the differences in dark matter halo mass and concentration due to initialization method through three 2LPT and three ZA initialized cosmological simulations. We find that 2LPT induces more rapid halo growth, resulting in more massive halos compared to ZA. This effect is most pronounced for high mass halos and at high redshift. Halo concentration is, on average, largely similar between 2LPT and ZA, but retains differences when viewed as a function of halo mass. For both mass and concentration, the difference between typical individual halos can be very large, highlighting the shortcomings of ZA-initialized simulations for high-z halo population studies.

We construct equilibrium models of compact stars using a realistic equation of state and obtain the density range occupied by the proton superconductor in strong $B$-fields. We do so by combining the density profiles of our models with microscopic calculations of proton pairing gaps and the critical unpairing field $H_{c2}$ above which the proton type-II superconductivity is destroyed. We find that magnetars with interior homogeneous field within the range $0.1 \le B_{16}\le 2$, where $B_{16} = B/10^{16}$ G, are partially superconducting, whereas those with $B_{16} > 2$ are void of superconductivity. We briefly discuss the neutrino emissivity and superfluid dynamics of magnetars in the light of their (non)-superconductivity.

Recently, both simulations and observations have revealed that flybys - fast, one-time interactions between two galaxy halos - are surprisingly common, nearing/comparable to galaxy mergers. Since these are rapid, transient events with the closest approach well outside the galaxy disk, it is unclear if flybys can transform the galaxy in a lasting way. We conduct collisionless N-body simulations of three co-planer flyby interactions between pure-disk galaxies to take a first look at the effects flybys have on disk structure, with particular focus on stellar bar formation. We find that some flybys are capable of inciting a bar with bars forming in both galaxies during our 1:1 interaction and in the secondary during our 10:1 interaction. The bars formed have ellipticities >0.5, sizes on the order of the host disk's scale length, and persist to the end of our simulations, ~5 Gyr after pericenter. The ability of flybys to incite bar formation implies that many processes associated with secular bar evolution may be more closely tied with interactions than previously though.

We construct an equation of state of strange quark matter in strong magnetic field within a confining model. The confinement is modeled by means of the Richardson potential for quark-quark interaction modified suitably to account for strong magnetic field. We compare our results for the equation of state and magnetization of matter to those derived within the MIT bag model. The differences between these models arise mainly due to the momentum dependence of the strong interaction between quarks in the Richardson model. Specifically, we find that the magnetization of strange quark matter in this model has much more pronounced de Haas-van Alf\'{v}en oscillations than in the MIT bag model, which is the consequence of the (static) gluon-exchange structure of the confining potential.

Hierarchical structure formation theory is based on the notion that mergers drive galaxy evolution, so a considerable framework of semi-analytic models and N-body simulations has been constructed to calculate how mergers transform a growing galaxy. However, galaxy mergers are only one type of major dynamical interaction between halos -- another class of encounter, a close flyby, has been largely ignored. We analyze a 50 Mpc/h, $1024^3$ collisionless cosmological simulation and find that the number of close flyby interactions is comparable to, or even surpasses, the number of mergers for halo masses $\ga 10^{11}\,{h^{-1} M_\odot}$ at $z \la 2$. Halo flybys occur so frequently to high mass halos that they are continually perturbed, unable to reach a dynamical equilibrium. We also find tentative evidence that at high redshift, $z \ga 14$, flybys are as frequent as mergers. Our results suggest that close halo flybys can play an important role in the evolution of the earliest dark matter halos and their galaxies, and can still influence galaxy evolution at the present epoch.

Compact stars with strong magnetic fields (magnetars) have been observationally determined to have surface magnetic fields of order of 10^{14}-10^{15} G, the implied internal field strength being several orders larger. We study the equation of state and composition of dense hypernuclear matter in strong magnetic fields in a range expected in the interiors of magnetars. Within the non-linear Boguta-Bodmer-Walecka model we find that the magnetic field has sizable influence on the properties of matter for central magnetic field B \ge 10^{17} G, in particular the matter properties become anisotropic. Moreover, for the central fields B \ge 10^{18} G, the magnetized hypernuclear matter shows instability, which is signaled by the negative sign of the derivative of the pressure parallel to the field with respect to the density, and leads to vanishing parallel pressure at the critical value B_{\rm cr} \simeq 10^{19} G. This limits the range of admissible homogeneously distributed fields in magnetars to fields below the critical value B_{\rm cr}.

The magnetars are believed to be highly magnetized neutron stars having surface magnetic field 10^{14} - 10^{15} G. It is believed that at the center, the magnetic field may be higher than that at the surface. We study the effect of the magnetic field on the neutron star matter. We model the nuclear matter with the relativistic mean field approach considering the possibility of appearance of hyperons at higher density. We find that the effect of magnetic field on the matter of neutron stars and hence on the mass-radius relation is important, when the central magnetic field is atleast of the order of 10^{17} G. Very importantly, the effect of strong magnetic field reveals anisotropy to the system. Moreover, if the central field approaches 10^{19} G, then the matter becomes unstable which limits the maximum magnetic field at the center of magnetars.

Observations of the Galactic Center (GC) have accumulated a multitude of "forensic" evidence indicating that several million years ago the center of the Milky Way galaxy was teaming with starforming and accretion-powered activity -- this paints a rather different picture from the GC as we understand it today. We examine a possibility that this epoch of activity could have been triggered by the infall of a satellite galaxy into the Milky Way which began at the redshift of 10 and ended few million years ago with a merger of the Galactic supermassive black hole with an intermediate mass black hole brought in by the inspiralling satellite.

We explore structure formation in the dark ages ($z\sim 30-6$) using two well-known methods for initializing cosmological $N$-body simulations. Overall, both the Zel'dovich approximation (\za) and second order Lagrangian perturbation theory (\lpt) are known to produce accurate present-day dark matter halo mass functions. However, since the \lpt method drives more rapid evolution of dense regions, it increases the occurrence of rare massive objects -- an effect that is most pronounced at high redshift. We find that \lpt produces more halos that could harbor Population III stars and their black hole remnants, and they produce them earlier. Although the differences between the \lpt and \za mass functions are nearly erased by $z=6$, this small boost to the number and mass of black holes more than doubles the reionized volume of the early Universe. We discuss the implications for reionization and massive black hole growth.

We present a comparative study of recent works on merger-timescales with dynamical friction and find a strong contrast between idealized/isolated mergers (Boylan-Kolchin et al. 2008) and mergers from a cosmological volume (Jiang et al. 2008). Our study measures the duration of mergers in a cosmological N-body simulation of dark matter, with emphasis on higher redshifts (z < 10) and a lower mass range. In our analysis we consider and compare two merger definitions; tidal disruption and coalescence. We find that the merger-time formula proposed by Jiang et al. (2008) describes our results well and conclude that cosmologically motivated merger-time formulae provide a more versatile and statistically robust approximation for practical applications such as semi-analytic/hybrid models.

Hierarchical structure formation theory is based on the notion that mergers drive galaxy evolution, so a considerable framework of semi-analytic models and N-body simulations has been constructed to calculate how mergers transform a growing galaxy. However, galaxy mergers are only one type of major dynamical interaction between halos -- another class of encounter, a close flyby, has been largely ignored. We use cosmological N-body simulations to reconstruct the entire dynamical interaction history of dark matter halos. We present a careful method of identifying and tracking a dark matter halo which resolves the typical classes of anomalies that occur in N-body data. This technique allows us to robustly follow halos and several hierarchical levels of subhalos as they grow, dissolve, merge, and flyby one another -- thereby constructing both a census of the dynamical interactions in a volume and an archive of the dynamical evolution of an individual halo. In addition to a census of mergers, our tool characterizes the frequency of close flyby interactions in the Universe. We find that the number of close flyby interactions is comparable to, or even surpasses, the number of mergers for halo masses $\gtrsim 10^{11}\,\Msun$ at $z \lesssim 2$. Halo flybys occur so frequently to high mass halos that they are continually perturbed, unable to reach a dynamical equilibrium. In particular, we find that Milky Way type halos undergo a similar number of flybys as mergers irrespective of mass-ratio for $z\lesssim 2$. We also find tentative evidence that at high redshift, $z \gtrsim 14$, flybys are as frequent as mergers. Our results suggest that close halo flybys can play an important role in the evolution of the earliest dark matter halos and their galaxies, and can still influence galaxy evolution at the present epoch. [Abridged]

Recent results and data suggests that high magnetic field in neutron stars (NS) strongly affects the characteristic (radius, mass) of the star. They are even separated as a class known as magnetars, for whom the surface magnetic field are greater than $10^{14}$ G. In this work we discuss the effect of such high magnetic field on the phase transition of NS to quark star (QS). We study the effect of magnetic field on the transition from NS to QS including the magnetic field effect in equation of state (EoS). The inclusion of the magnetic field increases the range of baryon number density, for which the flow velocities of the matter in the respective phase are finite. The magnetic field helps in initiation of the conversion process. The velocity of the conversion front however decreases due to the presence of magnetic field, as the presence of magnetic field reduces the effective pressure (P). The magnetic field of the star gets decreased by the conversion process, and the resultant QS has lower magnetic field than that of the initial NS.

Observations indicate that roughly 60% of the baryons may exist in a Warm-Hot Intergalactic Medium (WHIM) at low redshifts. Following up on previous results showing that gas is released through galaxy mergers, we use a semi-analytic technique to estimate the fraction of gas mass lost from haloes solely due to mergers. We find that up to ~25% of the gas in a halo can unbind over the course of galaxy assembly. This process does not act preferentially on smaller mass haloes; bigger haloes \emph{always} release larger amounts of gas in a given volume of the Universe. However, if we include multi-phase gas accretion onto haloes, we find that only a few percent is unbound. We conclude that either non-gravitational processes may be in play to heat up the gas in the galaxies prior to unbinding by mergers or most of the baryons in the WHIM have never fallen into virialised dark matter haloes. We present a budget for stocking the WHIM compiled from recent work.

We investigate the transition of a radiatively inefficient phase of a viscous two temperature accreting flow to a cooling dominated phase and vice versa around black holes. Based on a global sub-Keplerian accretion disc model in steady state, including explicit cooling processes self-consistently, we show that general advective accretion flow passes through various phases during its infall towards a black hole. Bremsstrahlung, synchrotron and inverse Comptonization of soft photons are considered as possible cooling mechanisms. Hence the flow governs a much lower electron temperature ~10^8 - 10^{9.5}K compared to the hot protons of temperature ~10^{10.2} - 10^{11.8}K in the range of the accretion rate in Eddington units 0.01 - 100. Therefore, the solutions may potentially explain the hard X-rays and the gamma-rays emitted from AGNs and X-ray binaries. We finally compare the solutions for two different regimes of viscosity and conclude that a weakly viscous flow is expected to be cooling dominated compared to its highly viscous counterpart which is radiatively inefficient. The flow is successfully able to reproduce the observed luminosities of the under-fed AGNs and quasars (e.g. Sgr A*), ultra-luminous X-ray sources (e.g. SS433), as well as the highly luminous AGNs and ultra-luminous quasars (e.g. PKS 0743-67) at different combinations of the mass accretion rate and ratio of specific heats.

We study the bulk viscosity of neutron star matter including $\Lambda$ hyperons in the presence of quantizing magnetic fields. Relaxation time and bulk viscosity due to both the non-leptonic weak process involving $\Lambda$ hyperons and direct Urca processes are calculated here. In the presence of a strong magnetic field of $10^{17}$ G, the hyperon bulk viscosity coefficient is reduced whereas bulk viscosity coefficients due to direct Urca processes are enhanced compared with their field free cases when many Landau levels are populated by protons, electrons and muons.

Galaxy merger simulations have explored the behaviour of gas within the galactic disk, yet the dynamics of hot gas within the galaxy halo has been neglected. We report on the results of high-resolution hydrodynamic simulations of colliding galaxies with hot halo gas. We explore a range of mass ratios, gas fractions and orbital configurations to constrain the shocks and gas dynamics within the progenitor haloes. We find that : (i) A strong shock is produced in the galaxy haloes before the first passage, increasing the temperature of the gas by almost an order of magnitude to $T\sim 10^{6.3}$ K. (ii) The X-ray luminosity of the shock is strongly dependent on the gas fraction; it is $\gtrsim 10^{39}$ erg/s for halo gas fractions larger than 10%. (iii) The hot diffuse gas in the simulation produces X-ray luminosities as large as $10^{42}$ erg/s. This contributes to the total X-ray background in the Universe. (iv) We find an analytic fit to the maximum X-ray luminosity of the shock as a function of merger parameters. This fit can be used in semi-analytic recipes of galaxy formation to estimate the total X-ray emission from shocks in merging galaxies. (v) $\sim$ 10-20% of the initial gas mass is unbound from the galaxies for equal-mass mergers, while $3-5%$ of the gas mass is released for the 3:1 and 10:1 mergers. This unbound gas ends up far from the galaxy and can be a feasible mechanism to enrich the IGM with metals.

High energy density ($\eps$) and temperature (T) links general relativity and hydrodynamics leading to a lower bound for the ratio of shear viscosity ($\eta$) and entropy density ($s$). We get the interesting result that the bound is saturated in the simple model for quark matter that we use for strange stars at the surface for $T \sim 80 MeV$. At this $T$ we have the possibility of cosmic separation of phases. At the surface of the star where the pressure is zero - the density $\eps$ has a fixed value for all stars of various masses with correspondingly varying central energy density $\eps_c$. Inside the star where this density is higher, the ratio of $\eta/s$ is larger and are like the known results found for perturbative QCD. This serves as a check of our calculation. The deconfined quarks at the surface of the strange star at $T = 80 MeV$ seem to constitute the most perfect interacting fluid permitted by nature.

We study the consequences of CPT and lepton number violation in neutrino sector. For CPT violation we take gravity with which neutrino and antineutrino couple differently. Gravity mixes neutrino and antineutrino in an unequal ratio to give two mass eigenstates. Lepton number violation interaction together with CPT violation gives rise to neutrino-antineutrino oscillation. Subsequently, we study the neutrino flavor mixing and oscillation under the influence of gravity. It is found that gravity changes flavor oscillation significantly which influences the relative abundance of different flavors in present universe. We show that the neutrinoless double beta decay rate is modified due to presence of gravity- the origin of CPT violation, as the mass of the flavor state is modified.

At finite temperature (T) there is a link with general relativity and hydrodynamics that leads to a lower bound for the ratio of shear viscosity and entropy density (\eta/s). We find that the bound is saturated in the simple model for quark matter that we use for strange stars at T = 80 MeV, at the surface of a strange star. At this T we have the possibility of cosmic separation of phases. We find that, although strongly correlated, the quark matter at the surface of strange stars constitute the most perfect interacting fluid permitted by nature. At the centre of the star, however, the density is higher and conditions are more like the results found for perturbative QCD.

The pions decouple from quarks while their coupling to nucleons $g_{\pi NN}$ remains unchanged at high density. This will have consequences for some recent papers.

It appears that there is a genuine shortage of radio pulsars with surface magnetic fields significantly smaller than $\sim 10^8$ Gauss. We propose that the pulsars with very low magnetic fields are actually strange stars locked in a state of minimum free energy and therefore at a limiting value of the magnetic field which can not be lowered by the system spontaneously.

We suggest that superbursts from some low mass X-ray binaries may be due to breaking and re-formation of diquark pairs, on the surface of realistic strange stars. Diquarks are expected to break up due to the explosion and shock of the thermonuclear process. The subsequent production of copious diquark pairing may produce sufficient energy to produce the superbursts.

Careful assessment of four good superburst candidates for GX 17+2 reveals that superburst is possible at near Eddington mass accretion rates. For the other seven stars, where superburst is found, there is the standard model of burning accumulated carbon from normal type I bursts of the accreting stars. However, there is the need for carbon, nitrogen and oxygen mass fraction (Z_{CNO}) which must be larger than Z_{CNO,\odot}, where the latter refers to the standard value found in the sun. Also it is very difficult to incorporate GX 17+2 into the standard picture of superbursts. In case of superbursts from strange stars, arising from broken quark pairs going over to diquarks at the surface of the star, these problems do not arise. Furthermore there is a natural explanation for the large value of ~ 1000 for \alpha which is defined in the literature as the ratio of energy released between normal bursts to the energy released during the normal burst. In the scenario for superbursts in strange stars it may be argued that the relatively smaller value of \alpha of ~ 440 indicates frequent recurrence of superbursts which is reflected in 4U 1636-53.

Surface tension (S) is due to the inward force experienced by particles at the surface and usually gravitation does not play an important role in this force. But in compact stars the gravitational force on the particles is very large and S is found to depend not only on the interactions in the strange quark matter, but also on the structure of the star, i.e. on its mass and radius. Indeed, it has been claimed recently that 511 keV photons observed by the space probe INTEGRAL from the galactic bulge may be due to electron-positron annihilation, and their source may be the positron cloud outside of an antiquark star. Such stars, if they exist, may also go a long way towards explaining away the antibaryon deficit of the universe. For that to happen S must be high enough to allow for survival of quark/antiquark stars born in early stages of the formation of the universe. High value of S may also assist explanation of delayed gamma-ray burst after a supernova explosion, as conversion from normal matter to strange matter takes place. The possibility of some implications from formation of surface waves are also discussed.

Chiral symmetry is restored at high density, quarks become nearly massless and pion, the Goldstone of the symmetry breaking decouples from the quarks. What happens at high density is important for finding the density dependence of Strange Quark Matter (SQM), - which in turn is relevant for understanding the structure of compact stars.

A realistic Equation of State (EOS) leads to strange stars (ReSS) which are compact in the mass radius plot, close to the Schwarzchild limiting line (Dey et al. 1998). We carry out a stability analysis under radial oscillations and compare with the EOS of other SS models. We find that the ReSS is stable and an M-R region can be identified to that effect.

Strange stars calculated from a realistic equation of state (EOS), that incorporate chiral symmetry restoration as well as deconfinement at high density show compact objects in the mass radius curve. We compare our calculations of incompressibility for this EOS with that of nuclear matter. One of the nuclear matter EOS has a continuous transition to ud-matter at about five times normal density. Another nuclear matter EOS incorporates density dependent coupling constants. From a look at the consequent velocity of sound, it is found that the transition to ud-matter seems necessary.

Skin vibration of ReSS and consequent resonance absorption can account for the absorption lines in the spectrum of X-ray emission from many compact stellar objects and in particular, the stars J1210$-$5226 and RXJ1856$-$3754. Observations of the X-ray spectrum of these stars is difficult to explain, if they are neutron stars.

X-ray bursts from compact stars is believed to be due to type I thermonuclear processes which are short lived, typically ~ 10 to 100 s. There are some low mass X-ray binaries (LMXB) like 4U 1820-30, 4U 1636-53, KS 1731-260 and Serpens X-1, known as Super Bursters (SB) which emit X-rays close to the Eddington luminosity limit for long periods of several hours. Recently there are reports of some long bursters (LB), which have bursts lasting 6-25 minutes, whereas the 4U 1735-44 has a burst period of 86 minutes. We suggest that these bursts from SB and LB may be due to breaking and re-formation of diquark pairs, on the surface of realistic strange quark stars. We use the beta equilibrated u, d and s quark model of Dey et al. (1998) and Li et al. (1999a and 1999b) and allow for spin dependent hyperfine interaction between quarks. The interaction produces pairing of specific colour-spin diquarks, leading to further lowering of energy by several MeV-s for each pair, on the average. Diquarks are expected to break up due to the explosion and shock of the TN process. The subsequent production of copious diquark pairing may produce sufficient energy to produce the very long bursts seen in SB or LB. The estimated total energy liberated, 10^{42} ergs, can be explained in our model with the calculated pair density ~ 0.275/fm^3 and a surface thickness of only half a micron, if the entire surface is involved. The depth of the surface involved in the process may be only few microns if the process is restricted to small part of the surface near the equator as suggested by Bildsten.

A realistic equation of state (EOS) leads to realistic strange stars (ReSS) which are compact in the mass radius plot, close to the Schwarzchild limiting line (Dey et al 1998). Many of the observed stars fit in with this kind of compactness, irrespective of whether they are X-ray pulsars, bursters or soft $\gamma$ repeaters or even radio pulsars. We point out that a change in the radius of a star can be small or large, when its mass is increasing and this depends on the position of a particular star on the mass radius curve. We carry out a stability analysis against radial oscillations and compare with the EOS of other strange star (SS) models. We find that the ReSS is stable and an M-R region can be identified to that effect.