BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Rory Smith

Rory Smith appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.GAastro-ph.COastro-ph.HEgr-qcastro-ph.IMastro-ph.SRhep-phMachine Learningphysics.comp-phphysics.ins-detphysics.optics

Catalog footprint

What is connected

53works
11topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

53 published item(s)

preprint2026arXiv

Analysis of GWTC-3 with fully precessing numerical relativity surrogate models

The third Gravitational-Wave Transient Catalog (GWTC-3) contains 90 binary coalescence candidates detected by the LIGO-Virgo-KAGRA Collaboration (LVK). We provide a re-analysis of binary black hole (BBH) events using a recently developed numerical relativity (NR) waveform surrogate model, NRSur7dq4, that includes all $\ell \leq 4$ spin-weighted spherical harmonic modes as well as the complete physical effects of precession. Properties of the remnant black holes' (BH's) mass, spin vector, and kick vector are found using an associated remnant surrogate model NRSur7dq4Remnant. Both NRSur7dq4 and NRSur7dq4Remnant models have errors comparable to numerical relativity simulations and allow for high-accuracy parameter estimates. We restrict our analysis to 47 BBH events that fall within the regime of validity of NRSur7dq4 (mass ratios greater than 1/6 and total masses greater than $60 M_{\odot}$). While for most of these events our results match the LVK analyses that were obtained using the semi-analytical models such as IMRPhenomXPHM and SEOBNRv4PHM, we find that for more than 20\% of events the NRSur7dq4 model recovers noticeably different measurements of black hole properties like the masses and spins, as well as extrinsic properties like the binary inclination and distance. For instance, GW150914_095045 exhibits noticeable differences in spin precession and spin magnitude measurements. Other notable findings include one event (GW191109_010717) that constrains the effective spin $χ_{eff}$ to be negative at a 99.3\% credible level and two events (GW191109_010717 and GW200129_065458) with well-constrained kick velocities. Furthermore, compared to the models used in the LVK analyses, NRSur7dq4 recovers a larger signal-to-noise ratio and/or Bayes factors for several events.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A New Method to Constrain the Appearance and Disappearance of Observed Jellyfish Galaxy Tails

We present a new approach to observationally constrain where the tails of Jellyfish (JF) galaxies in groups and clusters first appear and how long they remain visible with respect to the moment of their orbital pericenter. This is accomplished by measuring the distribution of their tail directions with respect to their host's center, and their distribution in a projected velocity-radius phase-diagram. We then model these observed distributions using a fast and flexible approach where JF tails are painted onto dark matter halos according to a simple parameterised prescription, and perform a Bayesian analysis to estimate the parameters. We demonstrate the effectiveness of our approach using observational mocks, and then apply it to a known observational sample of 106 JF galaxies with radio continuum tails located inside 68 hosts such as groups and clusters. We find that, typically, the radio continuum tails become visible on first infall when the galaxy reaches roughly three quarters of r$_{200}$, and the tails remain visible for a few hundred Myr after pericenter passage. Lower mass galaxies in more massive hosts tend to form visible tails further out and their tails disappear more quickly after pericenter. We argue that this indicates they are more sensitive to ram pressure stripping. With upcoming large area surveys of JF galaxies in progress, this is a promising new method to constrain the environmental conditions in which visible JF tails exist.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A novel approach for calculating galaxy rotation curves using spaxel cross-correlation and iterative smoothing

Precise measurements of the internal dynamics of galaxies have proven of great importance for understanding the internal dark matter distribution of galaxies. We present a novel method for measuring the line-of-sight (LOS) velocities across the face of galaxies by cross-correlation of spectral pixels (spaxels) and an iterative method of smoothing. On simulated data the method can accurately recover the input LOS velocities for different types of spectra (absorption line dominated, emission line dominated, and differing shapes of the continuum), and can handle stellar population radial gradients. Most important of all, it continues to provide reliable measurements of LOS velocities with reasonable uncertainties even when the spectra are very low signal-to-noise (approaching $\sim 1$), which is a challenge for traditional template-fitting approaches. We apply our method to data from a real MaNGA galaxy as a demonstration and find promising results with good precision. This novel approach can be complementary to existing methods primarily based on template fitting.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Brought to Light III: Colors of Disk and Clump Substructures in Dwarf Early-Type Galaxies of the Fornax Cluster

It has been well established that dwarf early-type galaxies (ETGs) can often exhibit a complex morphology, whereby faint spiral arms, bars, edge-on disks or clumps are embedded in their main, brighter diffuse body. In our first paper (Brought to Light I: Michea et al. 2021), we developed a new method for robustly identifying and extracting substructures in deep imaging data of dwarf ETGs in the Virgo galaxy cluster. Here we apply our method to a sample of 23 dwarf ETGs in the Fornax galaxy cluster, out of which 9 have disk-like and 14 have clump-like substructures. According to Fornax Deep Survey (FDS) data, our sample constitutes $12\%$ of all dwarf ETGs in Fornax brighter than $\text{M}_{r}=-13$ mag, and contains all cases that unequivocally exhibit substructure features. We use $g$ and $r$-band FDS images to measure the relative contribution of the substructures to the total galaxy light and to estimate their $g-r$ colors. We find that substructures typically contribute $8.7\%$ and $5.3\%$ of the total galaxy light in the $g$ and $r$ bands, respectively, within two effective radii. Disk substructures are usually found in dwarf ETGs with redder global colors, and they can be either as red as or bluer than their galaxy's diffuse component. In contrast, clump substructures are found in comparatively bluer dwarf ETGs, and they are always bluer than their galaxy's diffuse component. These results provide further evidence that dwarf ETGs can hide diverse complex substructures, with stellar populations that can greatly differ from those of the dominant diffuse light in which they are embedded.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Comparison of spatial distributions of Intracluster light and Dark Matter

In a galaxy cluster, the relative spatial distributions of dark matter, member galaxies, gas, and intracluster light (ICL) may connote their mutual interactions over the cluster evolution. However, it is a challenging problem to provide a quantitative measure for the shape matching between two multi-dimensional scalar distributions. We present a novel methodology, named the {\em Weighted Overlap Coefficient (WOC)}, to quantify the similarity of 2-dimensional spatial distributions. We compare the WOC with a standard method known as the Modified Hausdorff Distance (MHD). We find that our method is robust, and performs well even with the existence of multiple sub-structures. We apply our methodology to search for a visible component whose spatial distribution resembled with that of dark matter. If such a component could be found to trace the dark matter distribution with high fidelity for more relaxed galaxy clusters, then the similarity of the distributions could also be used as a dynamical stage estimator of the cluster. We apply the method to six galaxy clusters at different dynamical stages simulated within the GRT simulation, which is an N-body simulation using the galaxy replacement technique. Among the various components (stellar particles, galaxies, ICL), the ICL+ brightest cluster galaxy (BCG) component most faithfully trace the dark matter distribution. Among the sample galaxy clusters, the relaxed clusters show stronger similarity in the spatial distribution of the dark matter and ICL+BCG than the dynamically young clusters. While the MHD results show weaker trend with the dynamical stages.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Is Abell 2261 a fossil galaxy cluster in a transitional dynamical state?

Abell 2261 (A2661) is a well-studied fossil cluster, but previous studies give contradictory results on its dynamical states, such as its X-ray central entropy and magnitude gap. To improve our understanding of its dynamical state, we conduct multi-object spectroscopic observations with Hectospec on the MMT, covering an area out to 5 virial radii from the cluster center, and get improved completeness and membership. Using this new data, we calculate multiple dynamical indicators including gaussianity, distance offset, and velocity offset. These indicators suggest that A2261 is moderately relaxed. However, Dressler-Shectman test reveals a group candidate to the south, at a projected distance that is near the virial radius, and that overlaps with an X-ray tail-like feature. One of the galaxies associated with that group would be sufficiently bright to reduce the fossil magnitude gap. This raises the possibility that A2261 could have recently transit in fossil status, if the group had previously crossed the cluster and is only now found outside. In the cluster outskirts, we see an extended feature of galaxies located on the opposite side of the cluster to the group candidate. On even larger scales, we find this feature connects, both on the sky and in velocity space, with a long (~4.4Mpc) filamentary structure in the Sloan Digital Sky Survey data. This could support the idea that a group was fed into the cluster through the filament, temporarily breaking the fossil status, and resulting in a minor merger that weakly disturbed the intracluster medium of the cluster.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

LAAT: Locally Aligned Ant Technique for discovering multiple faint low dimensional structures of varying density

Dimensionality reduction and clustering are often used as preliminary steps for many complex machine learning tasks. The presence of noise and outliers can deteriorate the performance of such preprocessing and therefore impair the subsequent analysis tremendously. In manifold learning, several studies indicate solutions for removing background noise or noise close to the structure when the density is substantially higher than that exhibited by the noise. However, in many applications, including astronomical datasets, the density varies alongside manifolds that are buried in a noisy background. We propose a novel method to extract manifolds in the presence of noise based on the idea of Ant colony optimization. In contrast to the existing random walk solutions, our technique captures points that are locally aligned with major directions of the manifold. Moreover, we empirically show that the biologically inspired formulation of ant pheromone reinforces this behavior enabling it to recover multiple manifolds embedded in extremely noisy data clouds. The algorithm performance in comparison to state-of-the-art approaches for noise reduction in manifold detection and clustering is demonstrated, on several synthetic and real datasets, including an N-body simulation of a cosmological volume.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Linking the rates of neutron star binaries and short gamma-ray bursts

Short gamma-ray bursts are believed to be produced by both binary neutron star (BNS) and neutron star-black hole (NSBH) mergers. We use current estimates for the BNS and NSBH merger rates to calculate the fraction of observable short gamma-ray bursts produced through each channel. This allows us to constrain merger rates of BNS to $\mathcal{R}_{\rm{BNS}}=384^{+431}_{-213}{\rm{Gpc}^{-3} \rm{yr}^{-1}}$ ($90\%$ credible interval), a $16\%$ decrease in the rate uncertainties from the second LIGO--Virgo Gravitational-Wave Transient Catalog, GWTC-2. Assuming a top-hat emission profile with a large Lorentz factor, we constrain the average opening angle of gamma-ray burst jets produced in BNS mergers to $\approx 15^\circ$. We also measure the fraction of BNS and NSBH mergers that produce an observable short gamma-ray burst to be $0.02^{+0.02}_{-0.01}$ and $0.01 \pm 0.01$, respectively and find that $\gtrsim 40\%$ of BNS mergers launch jets (90\% confidence). We forecast constraints for future gravitational-wave detections given different modelling assumptions, including the possibility that BNS and NSBH jets are different. With $24$ BNS and $55$ NSBH observations, expected within six months of the LIGO-Virgo-KAGRA network operating at design sensitivity, it will be possible to constrain the fraction of BNS and NSBH mergers that launch jets with $10\%$ precision. Within a year of observations, we can determine whether the jets launched in NSBH mergers have a different structure than those launched in BNS mergers and rule out whether $\gtrsim 80\%$ of binary neutron star mergers launch jets. We discuss the implications of future constraints on understanding the physics of short gamma-ray bursts and binary evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Orbital dynamics and extreme scattering event properties from long-term scintillation observations of PSR J1603-7202

We model long-term variations in the scintillation of binary pulsar PSR J1603$-$7202, observed by the 64 m Parkes radio telescope (Murriyang) between 2004 and 2016. We find that the time variation in the scintillation arc curvature is well-modelled by scattering from an anisotropic thin screen of plasma between the Earth and the pulsar. Using our scintillation model, we measure the inclination angle and longitude of ascending node of the orbit, yielding a significant improvement over the constraints from pulsar timing. From our measurement of the inclination angle, we place a lower bound on the mass of J1603$-$7202's companion of $\gtrsim 0.5\,\text{M}_\odot$ assuming a pulsar mass of $\gtrsim1.2\,\text{M}_\odot$. We find that the scintillation arcs are most pronounced when the electron column density along the line of sight is increased, and that arcs are present during a known extreme scattering event. We measure the distance to the interstellar plasma and its velocity, and we discuss some structures seen in individual scintillation arcs within the context of our model.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Post-starburst galaxies in the centers of intermediate redshift clusters

We present results from MUSE spatially-resolved spectroscopy of 21 post-starburst galaxies in the centers of 8 clusters from $z\sim0.3$ to $z\sim0.4$. We measure spatially resolved star-formation histories (SFHs), the time since quenching ($t_Q$) and the fraction of stellar mass assembled in the past 1.5 Gyr ($μ_{1.5}$). The SFHs display a clear enhancement of star-formation prior to quenching for 16 out of 21 objects, with at least 10% (and up to $>50$%) of the stellar mass being assembled in the past 1.5 Gyr and $t_Q$ ranging from less than 100 Myrs to $\sim800$ Myrs. By mapping $t_Q$ and $μ_{1.5}$, we analyze the quenching patterns of the galaxies. Most galaxies in our sample have quenched their star-formation from the outside-in or show a side-to-side/irregular pattern, both consistent with quenching by ram-pressure stripping. Only three objects show an inside-out quenching pattern, all of which are at the high-mass end of our sample. At least two of them currently host an active galactic nucleus. In two post-starbursts, we identify tails of ionized gas indicating that these objects had their gas stripped by ram pressure very recently. Post-starburst features are also found in the stripped regions of galaxies undergoing ram-pressure stripping in the same clusters, confirming the link between these classes of objects. Our results point to ram-pressure stripping as the main driver of fast quenching in these environments, with active galactic nuclei playing a role at high stellar masses.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Spatial distribution of dark matter in and around galaxy clusters traced by galaxies, gas and intracluster stars in a simulated universe

To understand how well galaxies, gas and intracluster stars trace dark matter in and around galaxy clusters, we use the IllustrisTNG cosmological hydrodynamical simulation and compare the spatial distribution of dark matter with those of baryonic components in clusters. To quantify the global morphology of the density distribution of each component in clusters, we fit an ellipse to the density contour of each component and derive shape parameters at different radii. We find that ellipticity of dark matter is better correlated with that of galaxy mass-weighted number density, rather than with that of galaxy number density or galaxy velocity dispersion. We thus use the galaxy mass-weighted number density map as a representative of the galaxy maps. Among three different density maps from galaxies, gas, and intracluster stars, the ellipticity of dark matter is best reproduced by that of the galaxy map over the entire radii. The 'virialized' galaxy clusters show a better correlation of spatial distribution between dark matter and other components than the 'unvirialized' clusters, suggesting that it requires some time for each component to follow the spatial distribution of dark matter after merging events. Our results demonstrate that galaxies are still good tracers of dark matter distribution even in the non-linear regime corresponding to the scales in and around galaxy clusters, being consistent with the case where galaxies trace well the matter distribution in cosmologically large scales.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The First Fall is the Hardest: The Importance of Peculiar Galaxy Dynamics at Infall Time for Tidal Stripping Acting at the Centers of Groups and Clusters

Using dark matter only N-body cosmological simulations, we measure the pericentre distance of dark matter halos on their first infall into group and cluster halos. We find that the pericentre distance (R$_{\rm{peri}}$) is an important parameter as it significantly affects the strength of tidal mass loss in dense environments, and likely other environmental mechanisms as well. We examine what determines the R$_{\rm{peri}}$ value and find that, for most infallers, the dominant parameter is V$_{\rm{\perp}}$, the tangential component of the orbital velocity as the halo enters the group/cluster halo for the first time. This means that the strength of tidal stripping acting near the cores of groups/clusters are strongly influenced by the external peculiar velocity field of the large scale structure surrounding them, which differs between clusters, and is sensitive to the mass ratio of infaller to host. We find that filament feeding also partially contributes to feeding in low V$_{\rm{\perp}}$ halos. Dynamical friction can also play a role in reducing R$_{\rm{peri}}$ but this is only significant for those few relatively massive infallers ($>$10\% of the mass of their host). These results highlight how the response of galaxies to dense environments will sensitively depend on dynamics inherited from far outside those dense environments.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Unexpected dancing partners: Tracing the coherence between the spin and motion of dark matter halos

A recent study conducted using CALIFA survey data (Lee et al, 2019b) has found that the orbital motions of neighbor galaxies are coherent with the spin direction of a target galaxy on scales of many Megaparsecs. We study this so called `large-scale coherence' phenomena using N-body cosmological simulations. We confirm a strong coherence signal within 1 Mpc/h of a target galaxy, reaching out to 6 Mpc/h. We divide the simulation halos into subsamples based on mass, spin, merger history and local halo number density for both target and neighbor halos. We find a clear dependency on the mass of the target halo only. Another key parameter was the local number density of both target and neighbor halos, with high density regions such as clusters and groups providing the strongest coherence signals, rather than filaments or lower density regions. However we do not find a clear dependency on halo spin or time since last major merger. The most striking result we find is that the signal can be detected up to~15Mpc/h from massive halos. These result provide valuable lessons for how observational studies could more clearly detect coherence, and we discuss the implications of our results for the origins of large-scale coherence.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Exploring the AGN-ram pressure stripping connection in local clusters

Ram-pressure stripping by the intracluster medium (ICM) is one of the most advocated mechanisms that affect the properties of cluster galaxies. A recent study based on a small sample has found that many galaxies showing strong signatures of ram-pressure stripping also possess an active galactic nucleus (AGN), suggesting a possible correlation between the two phenomena. This result has not been confirmed by a subsequent study. Building upon previous findings, here we combine MUSE observations conducted within the GASP program and a general survey of the literature to robustly measure the AGN fraction in cluster's ram pressure stripped galaxies using BPT emission line diagrams. Considering a sample of 115 ram pressure stripped galaxies with stellar masses $\geq 10^9 \, M_{\odot}$, we find an AGN fraction of $\sim27\%$. This fraction strongly depends on stellar mass: it raises to 51% when only ram-pressure stripped galaxies of masses $M_* \geq 10^{10} \, M_{\odot}$ are considered. We then investigate whether the AGN incidence is in excess in ram pressure stripped galaxies compared to non-stripped galaxies, using as comparison a sample of non-cluster galaxies observed by the survey MaNGA. Considering mass-matched samples, we find that the incidence of AGN activity is significantly higher (at a confidence level $>99.95\%$) when ram-pressure stripping is on act, supporting the hypothesis of an AGN-ram pressure connection.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

GW190521 as a merger of Proca stars: a potential new vector boson of $8.7 \times 10^{-13}$ eV

Advanced LIGO-Virgo reported a short gravitational-wave signal (GW190521) interpreted as a quasi-circular merger of black holes, one populating the pair-instability supernova gap, forming a remnant black hole of $M_f\sim 142 M_\odot$ at a luminosity distance of $d_L \sim 5.3$ Gpc. With barely visible pre-merger emission, however, GW190521 merits further investigation of the pre-merger dynamics and even of the very nature of the colliding objects. We show that GW190521 is consistent with numerically simulated signals from head-on collisions of two (equal mass and spin) horizonless vector boson stars (aka Proca stars), forming a final black hole with $M_f = 231^{+13}_{-17}\,M_\odot$, located at a distance of $d_L = 571^{+348}_{-181}$ Mpc. The favoured mass for the ultra-light vector boson constituent of the Proca stars is $μ_{\rm V}= 8.72^{+0.73}_{-0.82}\times10^{-13}$ eV. This provides the first demonstration of close degeneracy between these two theoretical models, for a real gravitational-wave event. Confirmation of the Proca star interpretation, which we find statistically slightly preferred, would provide the first evidence for a long sought dark matter particle.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Multiwavelength Analysis of A1240, the Double Radio Relic Merging Galaxy Cluster Embedded in a ~80 Mpc-long Cosmic Filament

We present a multiwavelength study of the double radio relic cluster A1240 at z=0.195. Our Subaru-based weak lensing analysis detects three mass clumps forming a ~4 Mpc filamentary structure elongated in the north-south orientation. The northern ($M_{200}=2.61_{-0.60}^{+0.51}\times10^{14} M_{\odot}$) and middle ($M_{200}=1.09_{-0.43}^{+0.34}\times10^{14} M_{\odot}$) mass clumps separated by ~1.3 Mpc are associated with A1240 and co-located with the X-ray peaks and cluster galaxy overdensities revealed by Chandra and MMT/Hectospec observations, respectively. The southern mass clump ($M_{200}=1.78_{-0.55}^{+0.44}\times10^{14} M_{\odot}$), ~1.5 Mpc to the south of the middle clump, coincides with the galaxy overdensity in A1237, the A1240 companion cluster at z=0.194. Considering the positions, orientations, and polarization fractions of the double radio relics measured by the LOFAR study, we suggest that A1240 is a post-merger binary system in the returning phase with the time-since-collision ~1.7 Gyr. With the SDSS DR16 data analysis, we also find that A1240 is embedded in the much larger-scale (~80 Mpc) filamentary structure whose orientation is in remarkable agreement with the hypothesized merger axis of A1240.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Temperature dependent appearance of exotic matter makes nascent neutron stars spin faster

Neutron stars offer the opportunity to study the behaviour of matter at densities and temperatures inaccessible to terrestrial experiments. Gravitational-wave observations of binary neutron star coalescences can constrain the neutron-star equation of state before and after merger. After the neutron star binary merges, hyperons can form in the remnant, changing the behaviour of the neutron-star equation of state. In this study, we use finite-entropy equations of state to show that a post-merger remnant can spin up due to cooling. The magnitude of the spin-up depends on the neutron-star equation of state. If hyperons are present, the post-merger spin-up changes the peak gravitational-wave frequency by $\sim 540$ Hz, when the entropy per baryon drops from $s=2$ $k_B$ to $s=0$ $k_B$. If hyperons are not present, the post-merger spin-up changes by $\sim 360$ Hz, providing a gravitational-wave signature for exotic matter. We expect the same qualitative behaviour whenever temperature dependent phase transitions are triggered.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The Galaxy Replacement Technique (GRT): a New Approach to Study Tidal Stripping and Formation of Intracluster Light in a Cosmological Context

We introduce the Galaxy Replacement Technique (GRT) that allows us to model tidal stripping of galaxies with very high-mass (m$_{\rm{star}}=5.4\times10^4$~M$_\odot$/h) and high-spatial resolution (10 pc/h), in a fully cosmological context, using an efficient and fast technique. The technique works by replacing multiple low-resolution DM halos in the base cosmological simulation with high-resolution models, including a DM halo and stellar disk. We apply the method to follow the hierarchical build-up of a cluster since redshift $\sim8$ to now, through the hierarchical accretion of galaxies, individually or in substructures such as galaxy groups. We find we can successfully reproduce the observed total stellar masses of observed clusters since redshift $\sim$1. The high resolution allows us to accurately resolve the tidal stripping process and well describe the formation of ultra-low surface brightness features in the cluster ($μ_{V}<32$ mag arcsec$^{-2}$) such as the intracluster light (ICL), shells and tidal streams. We measure the evolution of the fraction of light in the ICL and brightest cluster galaxy (BCG) using several different methods. While their broad response to the cluster mass growth history is similar, the methods show systematic differences, meaning we must be careful when comparing studies that use distinct methods. The GRT represents a powerful new tool for studying tidal effects on galaxies and exploring the formation channels of the ICL in a fully cosmological context and with large samples of simulated groups and clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

An analysis and visualization of the output mode-matching requirements for squeezing in Advanced LIGO and future gravitational wave detectors

The sensitivity of ground-based gravitational wave (GW) detectors will be improved in the future via the injection of frequency-dependent squeezed vacuum. The achievable improvement is ultimately limited by losses of the interferometer electromagnetic field that carries the GW signal. The analysis and reduction of optical loss in the GW signal chain will be critical for optimal squeezed light-enhanced interferometry. In this work we analyze a strategy for reducing output-side losses due to spatial mode mismatch between optical cavities with the use of adaptive optics. Our goal is not to design a detector from the top down, but rather to minimize losses within the current design. Accordingly, we consider actuation on optics already present and one transmissive optic to be added between the signal recycling mirror and the output mode cleaner. The results of our calculation show that adaptive mode-matching with the current Advanced LIGO design is a suitable strategy for loss reduction that provides less than 2% mean output mode-matching loss. The range of actuation required is +47 uD on SR3, +140 mD on OM1 and OM2, +50 mD on the SRM substrate, and -50 mD on the added new transmissive optic. These requirements are within the demonstrated ranges of real actuators in similar or identical configurations to the proposed implementation. We also present a novel technique that graphically illustrates the matching of interferometer modes and allows for a quantitative comparison of different combinations of actuators.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Inferring the population properties of binary black holes from unresolved gravitational waves

The vast majority of compact binary mergers in the Universe produce gravitational waves that are too weak to yield unambiguous detections; they are unresolved. We present a method to infer the population properties of compact binaries -- such as their merger rates, mass spectrum, and spin distribution -- using both resolved and unresolved gravitational waves. By eliminating entirely the distinction between resolved and unresolved signals, we eliminate bias from selection effects. To demonstrate this method, we carry out a Monte Carlo study using an astrophysically motivated population of binary black holes. We show that some population properties of compact binaries are well constrained by unresolved signals after about one week of observation with Advanced LIGO at design sensitivity.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Measuring the neutron star equation of state with gravitational waves: the first forty binary neutron star mergers

Gravitational waves from binary neutron star coalescences contain rich information about matter at supranuclear densities encoded by the neutron star equation of state. We can measure the equation of state by analyzing the tidal interactions between neutron stars, which is quantified by the tidal deformability. Multiple merger events are required to probe the equation of state over a range of neutron star masses. The more events included in the analysis, the stronger the constraints on the equation of state. In this paper, we build on previous work to explore the constraints that LIGO and Virgo are likely to place on the neutron star equation of state by combining the first forty binary neutron star detections, a milestone we project to be reached during the first year of accumulated design-sensitivity data. We carry out Bayesian inference on a realistic mock dataset of binaries to obtain posterior distributions for neutron star tidal parameters. In order to combine posterior samples from multiple observations, we employ a random forest regressor, which allows us to efficiently interpolate the likelihood distribution. Assuming a merger rate of 1540 Gpc$^{-3}$ yr$^{-1}$ and a LIGO-Virgo detector network operating for one year at the sensitivity of the third-observation run, plus an additional eight months of design sensitivity, we find that the radius of a 1.4 $M_\odot$ neutron star can be constrained to $\sim 10$% at 90% confidence. At the same time, the pressure at twice the nuclear saturation density can be constrained to $\sim 45$ % at 90% confidence. Finally, we add an appendix following publication of the paper in the journal, showing the posterior distribution of the maximum neutron star mass allowed by the equation of state. We find that the maximum mass can be constrained to $\sim 0.3$ $M_\odot$ at 90% confidence.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Revealing the formation mechanism of the shell galaxy NGC474 with MUSE

Stellar shells around galaxies could provide precious insights into their assembly history. However, their formation mechanism remains poorly empirically constrained, in particular the type of galaxy collisions at their origin. We present MUSE@VLT data of the most prominent outer shell of NGC 474, to constrain its formation history. The stellar shell spectrum is clearly detected, with a signal-to-noise ratio of around 65 pix$^{-1}$. We use a full spectral fitting method to determine the line-of-sight velocity and the age and metallicity of the shell and associated point-like sources within the MUSE field of view. We detect six GC candidates and eight PN candidates which are all kinematically associated to the stellar shell. We show that the shell has an intermediate metallicity, [M/H] = $-0.83^{+0.12}_{-0.12}$ and a possible $α$-enrichment, [$α$/Fe] ~ 0.3. Assuming the material of the shell comes from a lower mass companion, and that the latter had no initial metallicity gradient, such a stellar metallicity would constrain the mass of the progenitor to be around 7.4 x 10^8 M$_\odot$, implying a merger mass ratio of about 1:100. However our census of PNs and earlier photometry of the shell would suggest a much higher ratio, around 1:20. Given the uncertainties, this difference is significant only at the ~1sigma level. We discuss the characteristics of the progenitor, in particular whether the progenitor could also be composed of stars from the low metallicity outskirts from a more massive galaxy. Ultimately, the presented data does not allow us to put a firm constraint on the progenitor mass. We show that at least two globular cluster candidates possibly associated with the shell are quite young, with ages below 1.5~Gyr. We also note the presence of a young (~1Gyr) stellar population in the center of NGC 474. The two may have resulted from the same event.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Blue Compact Dwarf Galaxy VCC 848 Formed by Dwarf-Dwarf Merging: HI Gas, Star Formation and Numerical Simulations

A clear link between a dwarf-dwarf merger event and enhanced star formation (SF) in the recent past was recently identified in the gas-dominated merger remnant VCC 848, offering by far the clearest view of a gas-rich late-stage dwarf-dwarf merger. We present a joint analysis of JVLA HI emission-line mapping, optical imaging and numerical simulations of VCC 848, in order to examine the impact of the merger on the stellar and gaseous distributions. VCC 848 has less than 30% of its HI gas concentrated within the central high-surface-brightness star-forming region, while the remaining HI is entrained in outlying tidal features. Particularly, a well-defined tidal arm reaches N(HI) comparable to the galaxy center but lacks SF. The molecular gas mass inferred from the current SF rate (SFR) dominates over the atomic gas mass in the central ~ 1.5 kpc. VCC 848 is consistent with being a main-sequence star-forming galaxy for its current stellar mass and SFR. The HII region luminosity distribution largely agrees with that of normal dwarf irregulars with similar luminosities, except that the brightest HII region is extraordinarily luminous. Our N-body/hydrodynamical simulations imply that VCC 848 is a merger between a gas-dominated primary progenitor and a gas-bearing star-dominated secondary. The progenitors had their first passage on a near-radial non-coplanar orbit more than 1 Gyr ago. The merger did not build up a core as compact as typical compact dwarfs with centralized starburst, which may be partly ascribed to the star-dominated nature of the secondary, and in a general sense, a negative stellar feedback following intense starbursts triggered at early stages of the merger.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Blue Compact Dwarf Galaxy VCC848 Formed by Dwarf-Dwarf Merging

It has long been speculated that many starburst or compact dwarf galaxies are resulted from dwarf-dwarf galaxy merging, but unequivocal evidence for this possibility has rarely been reported in the literature. We present the first study of deep optical broadband images of a gas-dominated blue compact dwarf galaxy (BCD) VCC848 (Mstar=2e8Msun) which hosts extended stellar shells and thus is confirmed to be a dwarf-dwarf merger. VCC848 is located in the outskirts of the Virgo Cluster. By analyzing the stellar light distribution, we found that VCC848 is the result of a merging between two dwarf galaxies with a primary-to-secondary mass ratio < ~ 5 for the stellar components and < ~ 2 for the presumed dark matter halos. The secondary progenitor galaxy has been almost entirely disrupted. The age-mass distribution of photometrically selected star cluster candidates in VCC848 implies that the cluster formation rate (CFR, proportional to star formation rate) was enhanced by a factor of ~ 7 - 10 during the past 1 Gyr. The merging-induced enhancement of CFR peaked near the galactic center a few hundred Myr ago and has started declining in the last few tens of Myr. The current star formation activities, as traced by the youngest clusters, mainly occur at large galactocentric distances (> ~ 1 kpc). The fact that VCC848 is still (atomic) gas-dominated after the period of most violent collision suggests that gas-rich dwarf galaxy merging can result in BCD-like remnants with extended atomic gas distribution surrounding a blue compact center, in general agreement with previous numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The SAMI Galaxy Survey: Kinematics of stars and gas in brightest group galaxies; the role of group dynamics

We study the stellar and gas kinematics of the brightest group galaxies (BGGs) in dynamically relaxed and unrelaxed galaxy groups for a sample of 154 galaxies in the SAMI galaxy survey. We characterize the dynamical state of the groups using the luminosity gap between the two most luminous galaxies and the BGG offset from the luminosity centroid of the group. We find that the misalignment between the rotation axis of gas and stellar components is more frequent in the BGGs in unrelaxed groups, although with quite low statistical significance. Meanwhile galaxies whose stellar dynamics would be classified as `regular rotators' based on their kinemetry are more common in relaxed groups. We confirm that this dependency on group dynamical state remains valid at fixed stellar mass and Sersic index. The observed trend could potentially originate from a differing BGG accretion history in virialised and evolving groups. Amongst the halo relaxation probes, the group BGG offset appears to play a stronger role than the luminosity gap on the stellar kinematic differences of the BGGs. However, both the group BGG offset and luminosity gap appear to roughly equally drive the misalignment between the gas and stellar component of the BGGs in one direction. This study offers the first evidence that the dynamical state of galaxy groups may influence the BGG's stellar and gas kinematics and calls for further studies using a larger sample with higher signal-to-noise.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Sami Galaxy Survey: stellar populations of passive spiral galaxies in different environment

We investigate the stellar populations of passive spiral galaxies as a function of mass and environment, using integral field spectroscopy data from the Sydney-AAO Multi-object Integral field spectrograph Galaxy Survey. Our sample consists of $52$ cluster passive spirals and $18$ group/field passive spirals, as well as a set of S0s used as a control sample. The age and [Z/H] estimated by measuring Lick absorption line strength indices both at the center and within $1R_{\rm e}$ do not show a significant difference between the cluster and the field/group passive spirals. However, the field/group passive spirals with log(M$_\star$/M$_\odot)\gtrsim10.5$ show decreasing [$α$/Fe] along with stellar mass, which is $\sim0.1$ dex smaller than that of the cluster passive spirals. We also compare the stellar populations of passive spirals with S0s. In the clusters, we find that passive spirals show slightly younger age and lower [$α$/Fe] than the S0s over the whole mass range. In the field/group, stellar populations show a similar trend between passive spirals and S0s. In particular, [$α$/Fe] of the field/group S0s tend to be flattening with increasing mass above log(M$_\star$/M$_\odot)\gtrsim10.5$, similar to the field/group passive spirals. We relate the age and [$α$/Fe] of passive spirals to their mean infall time in phase-space; we find a positive correlation, in agreement with the prediction of numerical simulations. We discuss the environmental processes that can explain the observed trends. The results lead us to conclude that the formation of the passive spirals and their transformation into S0s may significantly depend on their environments.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

YZiCS: Unveiling Quenching History of Cluster Galaxies Using Phase-space Analysis

We used the time since infall (TSI) of galaxies, obtained from the Yonsei Zoom-in Cluster Simulation, and the star formation rate (SFR) from the Sloan Digital Sky Survey (SDSS) Data Release 10 to study how quickly star formation of disk galaxies is quenched in cluster environments. We first confirm that both simulated and observed galaxies are consistently distributed in phase space. We then hypothesize that the TSI and SFR are causally connected; thus, both the TSI and SFR of galaxies at each position of phase space can be associated through abundance matching. Using a flexible model, we derive the star formation history (SFH) of cluster galaxies that best reproduces the relationship between the TSI and SFR at $z\sim 0.08$. According to this SFH, we find that the galaxies with $M_{*} > 10^{9.5} M_{\odot}$ generally follow the so-called "delayed-then-rapid" quenching pattern. Our main results are as following: (i) Part of the quenching takes place outside clusters through mass quenching and pre-processing. The e-folding timescale of this "$ex\text{-}situ$ quenching phase" is roughly 3 Gyr with a strong inverse mass dependence. (ii) The pace of quenching is maintained roughly for 2 Gyr ("delay time") during the first crossing time into the cluster. During the delay time, quenching remains gentle probably because gas loss happens primarily on hot and neutral gases. (iii) Quenching becomes more dramatic (e-folding timescale of roughly 1 Gyr) after delay time, probably because ram pressure stripping is strongest near the cluster center. Counter-intuitively, more massive galaxies show shorter quenching timescales mainly because they enter their clusters with lower gas fractions due to $ex\text{-}situ$ quenching.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The Gas Accretion History of Low Mass Halos within the Cosmic Web from Cosmological Simulations

Using high resolution hydrodynamical cosmological simulations, we study the gas accretion history of low mass halos located in a field-like, low density environment. We track their evolution individually from the early, pre-reionization era, through reionization, and beyond until $z=0$. Before reionization, low mass halos accrete cool cosmic web gas at a very rapid rate, often reaching the highest gas mass they will ever have. But when reionization occurs, we see that almost all halos lose significant quantities of their gas content, although some respond less quickly than others. We find that the response rate is influenced by halo mass first, and secondarily by their internal gas density at the epoch of reionization. Reionization also fully ionises the cosmic web gas by z$\sim$6. As a result, the lowest mass halos (M$\sim$10$^6~h^{-1}$M$_\odot$ at $z=6$) can never again re-accrete gas from the cosmic web, and by $z\sim5$ have lost all their internal gas to ionisation, resulting in a halt in star formation at this epoch. However, more massive halos can recover from their gas mass loss, and re-accrete ionised cosmic web gas. We find the efficiency of this re-accretion is a function of halo mass first, followed by local surrounding gas density. Halos that are closer to the cosmic web structure can accrete denser gas more rapidly. We find that our lower mass halos have a sweet spot for rapid, dense gas accretion at distances of roughly 1-5 virial radii from the most massive halos in our sample ($>$10$^8~h^{-1}$M$_\odot$), as these tend to be embedded deeply within the cosmic web.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The impact of the dynamical state of galaxy groups on the stellar populations of central galaxies

We study the stellar populations of the brightest group galaxies (BGGs) in groups with different dynamical states, using GAMA survey data. We use two independent, luminosity dependent indicators to probe the relaxedness of their groups; the magnitude gap between the two most luminous galaxies ($ΔM_{12}$), and offset between BGG and the luminosity center ($D_{offset}$) of the group. Combined, these two indicators were previously found useful for identifying relaxed and unrelaxed groups. We find that the BGGs of unrelaxed groups have significantly bluer NUV-r colours than in relaxed groups. This is also true at the fixed sersic index. We find the bluer colours cannot be explained away by differing dust fraction, suggesting there are real differences in their stellar populations. SFRs derived from SED-fitting tend to be higher in unrelaxed systems. This is in part because of a greater fraction of BGGs with non-elliptical morphology, but also because unrelaxed systems have larger numbers of mergers, some of which may bring fuel for star formation. The SED-fitted stellar metallicities of BGGs in unrelaxed systems also tend to be higher by around 0.05 dex, perhaps because their building blocks were more massive. We find that the $ΔM_{12}$ parameter is the most important parameter behind the observed differences in the relaxed/unrelaxed groups, in contrast with the previous study of Trevisan et al. (2017). We also find that groups selected to be unrelaxed using our criteria tend to have higher velocity offsets between the BGG and their group.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

BUDHIES III: The fate of HI and the quenching of galaxies in evolving environments

In a hierarchical Universe clusters grow via the accretion of galaxies from the field, groups and even other clusters. As this happens, galaxies can lose their gas reservoirs via different mechanisms, eventually quenching their star-formation. We explore the diverse environmental histories of galaxies through a multi-wavelength study of the combined effect of ram-pressure stripping and group "processing" in Abell 963, a massive growing cluster at $z=0.2$ from the Blind Ultra Deep HI Environmental Survey (BUDHIES). We incorporate hundreds of new optical redshifts (giving a total of 566 cluster members), as well as Subaru and XMM-Newton data from LoCuSS, to identify substructures and evaluate galaxy morphology, star-formation activity, and HI content (via HI deficiencies and stacking) out to $3\times R_{200}$. We find that Abell 963 is being fed by at least 7 groups, that contribute to the large number of passive galaxies outside the cluster core. More massive groups have a higher fraction of passive and HI-poor galaxies, while low-mass groups host younger (often interacting) galaxies. For cluster galaxies not associated with groups we corroborate our previous finding that HI gas (if any) is significantly stripped via ram-pressure during their first passage through the intra-cluster medium, and find mild evidence for a starburst associated with this event. In addition, we find an overabundance of morphologically peculiar and/or star-forming galaxies near the cluster core. We speculate that these arise as groups pass through the cluster (post-processing). Our study highlights the importance of environmental quenching and the complexity added by evolving environments.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Could Segue 1 be a destroyed star cluster? - a dynamical perspective

We attempt to find a progenitor for the ultra-faint object Segue 1 under the assumption that it formed as a dark matter free star cluster in the past. We look for orbits, using the elongation of Segue 1 on the sky as a tracer of its path. Those orbits are followed backwards in time to find the starting points of our N-body simulations. The successful orbit, with which we can reproduce Segue 1 has a proper motion of mu_alpha = -0.19 mas/yr and mu_delta = -1.9 mas/yr, placing Segue 1 near its apo-galacticon today. Our best fitting model has an initial mass of 6224 Msun and an initial scale-length of 5.75 pc.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Fast and Accurate Inference on Gravitational Waves from Precessing Compact Binaries

Inferring astrophysical information from gravitational waves emitted by compact binaries is one of the key science goals of gravitational-wave astronomy. In order to reach the full scientific potential of gravitational-wave experiments we require techniques to mitigate the cost of Bayesian inference, especially as gravitational-wave signal models and analyses become increasingly sophisticated and detailed. Reduced order models (ROMs) of gravitational waveforms can significantly reduce the computational cost of inference by removing redundant computations. In this paper we construct the first reduced order models of gravitational-wave signals that include the effects of spin-precession, inspiral, merger, and ringdown in compact object binaries, and which are valid for component masses describing binary neutron star, binary black hole and mixed binary systems. This work utilizes the waveform model known as "IMRPhenomPv2". Our ROM enables the use of a fast reduced order quadrature (ROQ) integration rule which allows us to approximate Bayesian probability density functions at a greatly reduced computational cost. We find that the ROQ rule can be used to speed up inference by factors as high as 300 without introducing systematic bias. This corresponds to a reduction in computational time from around half a year to a half a day, for the longest duration/lowest mass signals. The ROM and ROQ rule are available with the main inference library of the LIGO Scientific Collaboration, LALInference.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Demographics of galactic bulges in the SDSS database

We present a new database of our two-dimensional bulge-disk decompositions for 14,233 galaxies drawn from SDSS DR12 in order to examine the properties of bulges residing in the local universe ($0.005 < z < 0.05$). We performed decompositions in $g$ and $r$ bands by utilizing the {\sc{galfit}} software. The bulge colors and bulge-to-total ratios are found to be sensitive to the details in the decomposition technique, and hence we hereby provide full details of our method. The $g-r$ colors of bulges derived are almost constantly red regardless of bulge size except for the bulges in the low bulge-to-total ratio galaxies ($B/T_{\rm r} \lesssim 0.3$). Bulges exhibit similar scaling relations to those followed by elliptical galaxies, but the bulges in galaxies with lower bulge-to-total ratios clearly show a gradually larger departure in slope from the elliptical galaxy sequence. The scatters around the scaling relations are also larger for the bulges in galaxies with lower bulge-to-total ratios. Both the departure in slopes and larger scatters are likely originated from the presence of young stars. The bulges in galaxies with low bulge-to-total ratios show signs of a frosting of young stars so substantial that their luminosity-weighted Balmer-line ages are as small as 1\,Gyr in some cases. While bulges seem largely similar in optical properties to elliptical galaxies, they do show clear and systematic departures as a function of bulge-to-total ratio. The stellar properties and perhaps associated formation processes of bulges seem much more diverse than those of elliptical galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Preferential Tidal Stripping of Dark Matter versus Stars in Galaxies

Using high resolution hydrodynamical cosmological simulations, we conduct a comprehensive study of how tidal stripping removes dark matter and stars from galaxies. We find that dark matter is always stripped far more significantly than the stars -- galaxies that lose $\sim$80$\%$ of their dark matter, typically lose only 10$\%$ of their stars. This is because the dark matter halo is initially much more extended than the stars. As such, we find the stellar-to-halo size-ratio (measured using r$_{\rm{eff}}$/r$_{\rm{vir}}$) is a key parameter controlling the relative amounts of dark matter and stellar stripping. We use simple fitting formulae to measure the relation between the fraction of bound dark matter and fraction of bound stars. We measure a negligible dependence on cluster mass or galaxy mass. Therefore these formulae have general applicability in cosmological simulations, and are ideal to improve stellar stripping recipes in semi-analytical models, and/or to estimate the impact that tidal stripping would have on galaxies when only their halo mass evolution is known.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Formation Scenario for the Disk of Satellites: Accretion of Satellites during Mergers

The Disk of Satellites (DoS) observed in the Andromeda galaxy is a thin and extended group of satellites, nearly perpendicular to the disk plane, that share a common direction of rotation about the centre of Andromeda. Although a DoS is also observed in the Milky Way galaxy, the prevalance of such structures in more distant galaxies remains controversial. Explanations for the formation of such DoSs vary widely from filamentary infall, or flattening due to the potential field from large scale structure, to galaxy interactions in a Mondian paradigm. Here we present an alternative scenario -- during a merger, a galaxy may bring its own satellite population when merging with another galaxy. We demonstrate how, under the correct circumstances, during the coalescence of the two galaxies, the satellite population can be spread into an extended, flattened structure, with a common direction of rotation about the merger remnant. We investigate the key parameters of the interaction, and the satellite population, that are required to form a DoS in this scenario.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Accelerated gravitational-wave parameter estimation with reduced order modeling

Inferring the astrophysical parameters of coalescing compact binaries is a key science goal of the upcoming advanced LIGO-Virgo gravitational-wave detector network and, more generally, gravitational-wave astronomy. However, current parameter estimation approaches for such scenarios can lead to computationally intractable problems in practice. Therefore there is a pressing need for new, fast and accurate Bayesian inference techniques. In this letter we demonstrate that a reduced order modeling approach enables rapid parameter estimation studies. By implementing a reduced order quadrature scheme within the LIGO Algorithm Library, we show that Bayesian inference on the 9-dimensional parameter space of non-spinning binary neutron star inspirals can be sped up by a factor of 30 for the early advanced detectors' configurations. This speed-up will increase to about $150$ as the detectors improve their low-frequency limit to 10Hz, reducing to hours analyses which would otherwise take months to complete. Although these results focus on gravitational detectors, the techniques are broadly applicable to any experiment where fast Bayesian analysis is desirable.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

BUDHIES II: A phase-space view of HI gas stripping and star-formation quenching in cluster galaxies

We investigate the effect of ram-pressure from the intracluster medium on the stripping of HI gas in galaxies in a massive, relaxed, X-ray bright, galaxy cluster at z=0.2 from the Blind Ultra Deep HI Environmental Survey (BUDHIES). We use cosmological simulations, and velocity vs. position phase-space diagrams to infer the orbital histories of the cluster galaxies. In particular, we embed a simple analytical description of ram-pressure stripping in the simulations to identify the regions in phase-space where galaxies are more likely to have been sufficiently stripped of their HI gas to fall below the detection limit of our survey. We find a striking agreement between the model predictions and the observed location of HI-detected and non-detected blue (late-type) galaxies in phase-space, strongly implying that ram-pressure plays a key role in the gas removal from galaxies, and that this can happen during their first infall into the cluster. However, we also find a significant number of gas-poor, red (early-type) galaxies in the infall region of the cluster that cannot easily be explained with our model of ram-pressure stripping alone. We discuss different possible additional mechanisms that could be at play, including the pre-processing of galaxies in their previous environment. Our results are strengthened by the distribution of galaxy colours (optical and UV) in phase-space, that suggests that after a (gas-rich) field galaxy falls into the cluster, it will lose its gas via ram-pressure stripping, and as it settles into the cluster, its star formation will decay until it is completely quenched. Finally, this work demonstrates the utility of phase-space diagrams to analyze the physical processes driving the evolution of cluster galaxies, in particular HI gas stripping.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Demographics of Isolated Galaxies Along the Hubble Sequence

Isolated galaxies in low-density regions are significant in the sense that they are least affected by the hierarchical pattern of galaxy growth, and interactions with perturbers, at least for the last few Gyr. To form a comprehensive picture of the star formation history of isolated galaxies, we constructed a catalog of isolated galaxies and their comparison sample in relatively denser environments. The galaxies are drawn from the SDSS DR7 in the redshift range of $0.025<z<0.044$. We performed a visual inspection and classified their morphology following the Hubble classification scheme. For the spectroscopic study, we make use of the OSSY catalog. We confirm most of the earlier understanding on isolated galaxies. The most remarkable additional results are as follows. Isolated galaxies are dominantly late type with the morphology distribution (E: S0: S: Irr) = (9.9: 11.3: 77.6: 1.2)\%. The frequency of elliptical galaxies among isolated galaxies is only a third of that of the comparison sample. Most of the photometric and spectroscopic properties are surprisingly similar between isolated and comparison samples. However, early-type isolated galaxies are less massive by 50\% and younger (by H$β$) by 20\% than their counterparts in the comparison sample. This can be explained as a result of different merger and star formation histories for differing environments in the hierarchical merger paradigm. We provide an on-line catalog for the list and properties of our sample galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Fast Simulation of Gaussian-Mode Scattering for Precision Interferometry

Understanding how laser light scatters from realistic mirror surfaces is crucial for the design, com- missioning and operation of precision interferometers, such as the current and next generation of gravitational-wave detectors. Numerical simulations are indispensable tools for this task but their utility can in practice be limited by the computational cost of describing the scattering process. In this paper we present an efficient method to significantly reduce the computational cost of optical simulations that incorporate scattering. This is accomplished by constructing a near optimal representation of the complex, multi-parameter 2D overlap integrals that describe the scattering process (referred to as a reduced order quadrature). We demonstrate our technique by simulating a near-unstable Fabry-Perot cavity and its control signals using similar optics to those installed in one of the LIGO gravitational-wave detectors. We show that using reduced order quadrature reduces the computational time of the numerical simulation from days to minutes (a speed-up of $\approx 2750 \times$) whilst incurring negligible errors. This significantly increases the feasibility of modelling interferometers with realistic imperfections to overcome current limits in state-of-the-art optical systems. Whilst we focus on the Hermite-Gaussian basis for describing the scattering of the optical fields, our method is generic and could be applied with any suitable basis. An implementation of this reduced order quadrature method is provided in the open source interferometer simulation software Finesse.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

MUSE-ings on AM1354-250: Collisions, Shocks and Rings

We present MUSE observations of AM1354-250, confirming its status as a collisional ring galaxy which has recently undergone an interaction, creating its distinctive shape. We analyse the stellar and gaseous emission throughout the galaxy finding direct evidence that the gaseous ring is expanding with a velocity of $\sim$70km.s$^{-1}$ and that star formation is occurring primarily in HII regions associated with the ring. This star formation activity is likely triggered by the interaction. We find evidence for several excitation mechanisms in the gas, including emission consistent with shocked gas in the expanding ring and a region of LINER-like emission in the central core of the galaxy. Evidence of kinematic disturbance in both the stars and gas, possibly also triggered by the interaction, can be seen in all of the velocity maps. The ring galaxy retains weak spiral structure, strongly suggesting the progenitor galaxy was a massive spiral prior to the collision with its companion an estimated $140 \pm 12$ Myr ago.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Origin of the Metallicity Distribution in the Thick Disc

Aims. Using a suite of cosmological chemodynamical disc galaxy simulations, we assess how (a) radial metallicity gradients evolve with scaleheight; (b) the vertical metallicity gradients change through the thick disc; and (c) the vertical gradient of the stellar rotation velocity varies through the disc. We compare with the Milky Way to search for analogous trends. Methods. We analyse five simulated spiral galaxies with masses comparable to the Milky Way. The simulations span a range of star formation and energy feedback strengths and prescriptions, particle- and grid-based hydrodynamical implementations, as well as initial conditions/assembly history. Results. Consistently, we find that the steeper, negative, radial metallicity gradients seen in the mid-plane flatten with increasing height away from the plane. In simulations with stronger (and/or more spatially-extended) feedback, the negative radial gradients invert, becoming positive for heights in excess of 1 kpc. Such behaviour is consistent with that inferred from recent observations. Our measurements of the vertical metallicity gradients show no clear correlation with galactocentric radius, and are in good agreement with those observed in the Milky Way's thick disc (locally). Conclusions. Simulations employing stronger/more extended feedback prescriptions possess radial and vertical metallicity and kinematic gradients more in line with recent observations. The inverted, positive, radial metallicity gradients seen in the simulated thick stellar discs originate from a population of younger, more metal-rich, stars formed in-situ, superimposed upon a background population of older migrators from the inner disc; the contrast provided by the former increases radially, due to the inside-out growth of the disc. A similar behaviour may be responsible for the same flattening seen in the radial gradients with scaleheight in the Milky Way.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Robust parameter estimation for compact binaries with ground-based gravitational-wave observations using the LALInference software library

The Advanced LIGO and Advanced Virgo gravitational wave (GW) detectors will begin operation in the coming years, with compact binary coalescence events a likely source for the first detections. The gravitational waveforms emitted directly encode information about the sources, including the masses and spins of the compact objects. Recovering the physical parameters of the sources from the GW observations is a key analysis task. This work describes the LALInference software library for Bayesian parameter estimation of compact binary signals, which builds on several previous methods to provide a well-tested toolkit which has already been used for several studies. We show that our implementation is able to correctly recover the parameters of compact binary signals from simulated data from the advanced GW detectors. We demonstrate this with a detailed comparison on three compact binary systems: a binary neutron star, a neutron star black hole binary and a binary black hole, where we show a cross-comparison of results obtained using three independent sampling algorithms. These systems were analysed with non-spinning, aligned spin and generic spin configurations respectively, showing that consistent results can be obtained even with the full 15-dimensional parameter space of the generic spin configurations. We also demonstrate statistically that the Bayesian credible intervals we recover correspond to frequentist confidence intervals under correct prior assumptions by analysing a set of 100 signals drawn from the prior. We discuss the computational cost of these algorithms, and describe the general and problem-specific sampling techniques we have used to improve the efficiency of sampling the compact binary coalescence parameter space.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Simple and accurate modelling of the gravitational potential produced by thick and thin exponential disks

We present accurate models of the gravitational potential produced by a radially exponential disk mass distribution. The models are produced by combining three separate Miyamoto-Nagai disks. Such models have been used previously to model the disk of the Milky Way, but here we extend this framework to allow its application to disks of any mass, scalelength, and a wide range of thickness from infinitely thin to near spherical (ellipticities from 0 to 0.9). The models have the advantage of simplicity of implementation, and we expect faster run speeds over a double exponential disk treatment. The potentials are fully analytical, and differentiable at all points. The mass distribution of our models deviates from the radial mass distribution of a pure exponential disk by <0.4% out to 4 disk scalelengths, and <1.9% out to 10 disk scalelengths. We tabulate fitting parameters which facilitate construction of exponential disks for any scalelength, and a wide range of disk thickness (a user-friendly, web-based interface is also available). Our recipe is well suited for numerical modelling of the tidal effects of a giant disk galaxy on star clusters or dwarf galaxies. We consider three worked examples; the Milky Way thin and thick disk, and a disky dwarf galaxy.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Sensitivity of Harassment to Orbit: Mass Loss from Early-Type Dwarfs in Galaxy Clusters

We conduct a comprehensive numerical study of the orbital dependence of harassment on early-type dwarfs consisting of 168 different orbits within a realistic, Virgo-like cluster, varying in eccentricity and pericentre distance. We find harassment is only effective at stripping stars or truncating their stellar disks for orbits that enter deep into the cluster core. Comparing to the orbital distribution in cosmological simulations, we find that the majority of the orbits (more than three quarters) result in no stellar mass loss. We also study the effects on the radial profiles of the globular cluster systems of early-type dwarfs. We find these are significantly altered only if harassment is very strong. This suggests that perhaps most early-type dwarfs in clusters such as Virgo have not suffered any tidal stripping of stars or globular clusters due to harassment, as these components are safely embedded deep within their dark matter halo. We demonstrate that this result is actually consistent with an earlier study of harassment of dwarf galaxies, despite the apparent contradiction. Those few dwarf models that do suffer stellar stripping are found out to the virial radius of the cluster at redshift=0, which mixes them in with less strongly harassed galaxies. However when placed on phase-space diagrams, strongly harassed galaxies are found offset to lower velocities compared to weakly harassed galaxies. This remains true in a cosmological simulation, even when halos have a wide range of masses and concentrations. Thus phase-space diagrams may be a useful tool for determining the relative likelihood that galaxies have been strongly or weakly harassed.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Possible Formation Scenario for Dwarf Spheroidal Galaxies I: Fiducial Model

We use numerical simulations to study a formation scenario for dwarf spheroidal galaxies in which their stellar populations are the products of the dissolution of open star clusters and stellar associations within cosmological dark matter haloes. This paper shows that this process gives rise to objects which resemble the observed dwarf spheroidal satellites of the Milky Way without invoking external influences. The presence of long-lived kinematic substructures within the stellar components of these objects affects their projected velocity dispersions. We find that this in turn affects mass estimates based on the projected velocity dispersion profiles which may over-estimate the actual dark matter halo mass depending on the amount of substructure which is present. Our models make predictions about the detailed kinematic and photometric properties of the dSphs which can be tested using future observations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Rapidly evaluating the compact binary likelihood function via interpolation

Bayesian parameter estimation on gravitational waves from compact binary coalescences (CBCs) typically requires millions of template waveform computations at different values of the parameters describing the binary. Sampling techniques such as Markov chain Monte Carlo and nested sampling evaluate likelihoods, and hence compute template waveforms, serially; thus, the total computational time of the analysis scales linearly with that of template generation. Here we address the issue of rapidly computing the likelihood function of CBC sources with non-spinning components. We show how to efficiently compute the continuous likelihood function on the three-dimensional subspace of parameters on which it has a non-trivial dependence (the chirp mass, symmetric mass ratio and coalescence time) via interpolation. Subsequently, sampling this interpolated likelihood function is a significantly cheaper computational process than directly evaluating the likelihood; we report improvements in computational time of two to three orders of magnitude while keeping likelihoods accurate to < 0.025 %. Generating the interpolant of the likelihood function over a significant portion of the CBC mass space is computationally expensive but highly parallelizable, so the wall time can be very small relative to the time of a full parameter-estimation analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Infant mortality in the hierarchical merging scenario: Dependency on gas expulsion timescales

We examine the effects of gas expulsion on initially sub-structured and out-of-equilibrium star clusters. We perform N-body simulations of the evolution of star clusters in a static background potential before adjusting that potential to model gas expulsion. We investigate the impact of varying the rate at which the gas is removed, and the instant at which gas removal begins. Reducing the rate at which the gas is expelled results in an increase in cluster survival. Quantitatively, this dependency is approximately in agreement with previous studies, despite their use of smooth, and virialised initial stellar distributions. However, the instant at which gas expulsion occurs is found to have a strong effect on cluster response to gas removal. We find if gas expulsion occurs prior to one crossing time, cluster response is poorly described by any global parameters. Furthermore in real clusters the instant of gas expulsion is poorly constrained. Therefore our results emphasis the highly stochastic and variable response of star clusters to gas expulsion.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Leo IV & V - A possible dwarf galaxy pair ?

The last few years have seen the discovery of many faint and ultra-faint dwarf spheroidal galaxies around the Milky Way. Among these is a pair of satellites called Leo IV and Leo V. This pair is found at large distances from the Milky Way (154 and 175 kpc respectively). The rather small difference in radial distance, and the fact that they also show a close projected distance on the sky, has led to the idea that we might be seeing a new pair of bound galaxies - like the Magellanic Clouds. In this paper we investigate this speculation by means of a simple integration code (confirming the results with full N-body simulations). As the luminous mass of both faint dwarfs is far too low to allow them to be bound, we simulate the pair assuming extended dark matter haloes. Our results show that the minimum dark matter mass required for the pair to be bound is rather high - ranging from 1.6 x 10^10 Msun to 5.4 x 10^10 Msun (within the virial radii). Computing the mass of dark matter within a commonly adopted radius of 300 pc shows that our models are well within the predicted range of dark matter content for satellites so faint. We therefore conclude that it could be possible that the two galaxies constitute a bound pair.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Numerical modelling of Auriga's Wheel - a new ring galaxy

We model the formation of Auriga's Wheel - a recently discovered collisional ring galaxy. Auriga's Wheel has a number of interesting features including a bridge of stars linking the neighbouring elliptical to the ring galaxy, and evidence for components of expansion and rotation within the ring. Using N-body/SPH modelling, we study collisions between an elliptical galaxy and a late-type disk galaxy. A near direct collision, with a mildy inclined disk, is found to reasonably reproduce the general system morphology ~50 Myr following the collision. The collision must have a relatively low velocity (initially ~150 km s^{-1}) in order to form the observed bridge, and simultaneously match the galaxies separation. Our best-match model suggests the total disk galaxy is ~5 times more massive than the elliptical. We find that the velocity of expansion of the ring is sensitive to the mass of the elliptical, while insensitive to the encounter velocity. We evolve our simulation beyond the current epoch to study the future destiny of the galaxy pair. In our model, the nucleus moves further away from the plane of the ring in the direction of the stellar bridge. The nucleus eventually merges with the elliptical galaxy ~100 Myr after the present time. The ring continues to expand for ~200 Myr before collapsing back. The low initial relative velocity of the two galaxies will eventually result in a complete merger.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The impact of galaxy harassment on the globular cluster systems of early-type cluster dwarf galaxies

The dynamics of globular cluster systems (GCSs) around galaxies are often used to assess the total enclosed mass, and even to constrain the dark matter distribution. The globular cluster system of a galaxy is typically assumed to be in dynamical equilibrium within the potential of the host galaxy. However cluster galaxies are subjected to a rapidly evolving and, at times, violently destructive tidal field. We investigate the impact of the harassment on the dynamics of GCs surrounding early type cluster dwarfs, using numerical simulations. We find that the dynamical behaviour of the GCS is strongly influenced by the fraction of bound dark matter f_{DM} remaining in the galaxy. Only when f_{DM} falls to ~15%, do stars and GCs begin to be stripped. Still the observed GC velocity dispersion can be used to measure the true enclosed mass to within a factor of 2, even when f_{DM} falls as low as ~3%. This is possible partly because unbound GCs quickly separate from the galaxy body. However even the distribution of {\it{bound}} GCs may spatially expand by a factor of 2-3. Once f_{DM} falls into the <3% regime, the galaxy is close to complete disruption, and GCS dynamics can no longer be used to reliably estimate the enclosed mass. In this regime, the remaining bound GCS may spatially expand by a factor of 4 to 8. It may be possible to test if a galaxy is in this regime by measuring the dynamics of the stellar disk. We demonstrate that if a stellar disk is rotationally supported, it is likely that a galaxy has sufficient dark matter, that the dynamics of the GCS can be used to reliably estimate the enclosed mass.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A New Collisional Ring Galaxy at z = 0.111: Auriga's Wheel

We report the serendipitous discovery of a collision ring galaxy, identified as 2MASX J06470249+4554022, which we have dubbed 'Auriga's Wheel', found in a SUPRIME-CAM frame as part of a larger Milky Way survey. This peculiar class of galaxies is the result of a near head-on collision between typically, a late type and an early type galaxy. Subsequent GMOS-N long-slit spectroscopy has confirmed both the relative proximity of the components of this interacting pair and shown it to have a redshift of 0.111. Analysis of the spectroscopy reveals that the late type galaxy is a LINER class Active Galactic Nuclei while the early type galaxy is also potentially an AGN candidate, this is very uncommon amongst known collision ring galaxies. Preliminary modeling of the ring finds an expansion velocity of ~200 kms^-1 consistent with our observations, making the collision about 50 Myr old. The ring currently has a radius of about 10 kpc and a bridge of stars and gas is also visible connecting the two galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Star Cluster collisions - a formation scenario for the Extended Globular Cluster Scl-dE1 GC1

Recent observations of the dwarf elliptical galaxy Scl-dE1 (Sc22) in the Sculptor group of galaxies revealed an extended globular cluster (Scl-dE1 GC1), which exhibits an extremely large core radius of about 21.2 pc. The authors of the discovery paper speculated on whether this object could reside in its own dark matter halo and/or if it might have formed through the merging of two or more star clusters. In this paper, we present N-body simulations to explore thoroughly this particular formation scenario. We follow the merger of two star clusters within dark matter haloes of a range of masses (as well as in the absence of a dark matter halo). In order to obtain a remnant which resembles the observed extended star cluster, we find that the star formation efficiency has to be quite high (around 33 per cent) and the dark matter halo, if present at all, has to be of very low mass, i.e. raising the mass to light ratio of the object within the body of the stellar distribution by at most a factor of a few. We also find that expansion of a single star cluster following mass loss provides another viable formation path. Finally, we show that future measurements of the velocity dispersion of this system may be able to distinguish between the various scenarios we have explored.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Surviving infant mortality in the hierarchical merging scenario

We examine the effects of gas expulsion on initially sub-structured and out-of-equilibrium star clusters. We perform $N$-body simulations of the evolution of star clusters in a static background potential before removing that potential to model gas expulsion. We find that the initial star formation efficiency is not a good measure of the survivability of star clusters. This is because the stellar distribution can change significantly, causing a large change in the relative importance of the stellar and gas potentials. We find that the initial stellar distribution and velocity dispersion are far more important parameters than the initial star formation efficiency, and that clusters with very low star formation efficiencies can survive gas expulsion. We suggest that it is variations in cluster initial conditions rather than in their star formation efficiencies that cause some clusters to be destroyed while a few survive.

0 citations0 reviewsNo code linkedOpen access