Researcher profile

Rory Smith

Rory Smith contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
29works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

29 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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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&#39; 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&#39;s stellar and gas kinematics and calls for further studies using a larger sample with higher signal-to-noise.

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.

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 &#34;delayed-then-rapid&#34; 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 &#34;$ex\text{-}situ$ quenching phase&#34; is roughly 3 Gyr with a strong inverse mass dependence. (ii) The pace of quenching is maintained roughly for 2 Gyr (&#34;delay time&#34;) 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.

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.

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.