Showing 1–50 of 112 results for author: D'Souza, R

First resolved stellar halo kinematics of a MW-mass galaxy outside the Local Group: A flat counter-rotating halo in NGC 4945

Authors: Camila Beltrand, Antonela Monachesi, Richard D'Souza, Eric F. Bell, Roelof S. de Jong, Facundo A. Gomez, Jeremy Bailin, In Sung Jang, Adam Smercina

Abstract: Stellar halos of galaxies, primarily formed through the accretion of smaller objects, are important to understand the hierarchical mass assembly of galaxies. However, the inner regions of stellar halos in disk galaxies are predicted to have an in-situ component that is expected to be prominent along the major axis. Kinematic information is crucial to disentangle the contribution of the in-situ com… ▽ More Stellar halos of galaxies, primarily formed through the accretion of smaller objects, are important to understand the hierarchical mass assembly of galaxies. However, the inner regions of stellar halos in disk galaxies are predicted to have an in-situ component that is expected to be prominent along the major axis. Kinematic information is crucial to disentangle the contribution of the in-situ component from the accreted stellar halos. The low surface brightness of stellar halos makes it inaccessible with traditional integrated light spectroscopy. In this work, using a novel technique, we study the kinematics of the stellar halo of the edge-on galaxy NGC 4945. We couple new deep Multi Unit Spectroscopic Explorer spectroscopic observations with existing Hubble Space Telescope imaging data to spectroscopically measure the line-of-sight (LOS) heliocentric velocity and velocity dispersion in two fields at a galactocentric distance of 12.2 kpc (outer disk field) and 34.6 kpc (stellar halo field) along NGC 4945 major axis, by stacking individual spectra of red giant branch and asymptotic giant branch stars. We obtain a LOS velocity and dispersion of 673+/-11 km/s and 73+/-14 km/s, respectively, for the outer disk field. This is consistent with the mean HI velocity of the disk at that distance. For the halo field we obtain a LOS velocity and dispersion of 519+/-12 km/s and 42+/-22 km/s. The halo fields' velocity measurement is within ~40 km/s from the systemic LOS velocity of NGC 4945, which is 563 km/s, suggesting that its stellar halo at 34.6 kpc along the major axis is counter-rotating and is of likely accretion origin. This provides the first ever kinematic measurement of the stellar halo of a Milky Way-mass galaxy outside the Local Group from its resolved stellar population, and establishes a powerful technique for measuring the velocity field of the stellar halos of nearby galaxies. △ Less

Submitted 25 June, 2024; originally announced June 2024.

Comments: Accepted to A&A. 23 pages, 12 figures, 2 appendices

The VISTA Variables in the Vía Láctea eXtended (VVVX) ESO public survey: Completion of the observations and legacy

Authors: R. K. Saito, M. Hempel, J. Alonso-García, P. W. Lucas, D. Minniti, S. Alonso, L. Baravalle, J. Borissova, C. Caceres, A. N. Chené, N. J. G. Cross, F. Duplancic, E. R. Garro, M. Gómez, V. D. Ivanov, R. Kurtev, A. Luna, D. Majaess, M. G. Navarro, J. B. Pullen, M. Rejkuba, J. L. Sanders, L. C. Smith, P. H. C. Albino, M. V. Alonso , et al. (121 additional authors not shown)

Abstract: The ESO public survey VISTA Variables in the Vía Láctea (VVV) surveyed the inner Galactic bulge and the adjacent southern Galactic disk from $2009-2015$. Upon its conclusion, the complementary VVV eXtended (VVVX) survey has expanded both the temporal as well as spatial coverage of the original VVV area, widening it from $562$ to $1700$ sq. deg., as well as providing additional epochs in… ▽ More The ESO public survey VISTA Variables in the Vía Láctea (VVV) surveyed the inner Galactic bulge and the adjacent southern Galactic disk from $2009-2015$. Upon its conclusion, the complementary VVV eXtended (VVVX) survey has expanded both the temporal as well as spatial coverage of the original VVV area, widening it from $562$ to $1700$ sq. deg., as well as providing additional epochs in $JHK_{\rm s}$ filters from $2016-2023$. With the completion of VVVX observations during the first semester of 2023, we present here the observing strategy, a description of data quality and access, and the legacy of VVVX. VVVX took $\sim 2000$ hours, covering about 4% of the sky in the bulge and southern disk. VVVX covered most of the gaps left between the VVV and the VISTA Hemisphere Survey (VHS) areas and extended the VVV time baseline in the obscured regions affected by high extinction and hence hidden from optical observations. VVVX provides a deep $JHK_{\rm s}$ catalogue of $\gtrsim 1.5\times10^9$ point sources, as well as a $K_{\rm s}$ band catalogue of $\sim 10^7$ variable sources. Within the existing VVV area, we produced a $5D$ map of the surveyed region by combining positions, distances, and proper motions of well-understood distance indicators such as red clump stars, RR Lyrae, and Cepheid variables. In March 2023 we successfully finished the VVVX survey observations that started in 2016, an accomplishment for ESO Paranal Observatory upon 4200 hours of observations for VVV+VVVX. The VVV+VVVX catalogues complement those from the Gaia mission at low Galactic latitudes and provide spectroscopic targets for the forthcoming ESO high-multiplex spectrographs MOONS and 4MOST. △ Less

Submitted 24 June, 2024; originally announced June 2024.

Comments: 17 pages, 11 figures (+ appendix). Accepted for publication in Astronomy and Astrophysics in section 14: Catalogs and data

Star formation beyond galaxies: widespread in-situ formation of intra-cluster stars in TNG50

Authors: Niusha Ahvazi, Laura V. Sales, Julio F. Navarro, Andrew Benson, Alessandro Boselli, Richard D'Souza

Abstract: We study the fraction of the intra-cluster light (ICL) formed in-situ in the three most massive clusters of the TNG50 simulation, with virial masses $\sim 10^{14}$ M$_{\odot}$. We find that a significant fraction of ICL stars ($8\%$-$28\%$) are born in-situ. This amounts to a total stellar mass comparable to the central galaxy itself. Contrary to simple expectations, only a sub-dominant fraction o… ▽ More We study the fraction of the intra-cluster light (ICL) formed in-situ in the three most massive clusters of the TNG50 simulation, with virial masses $\sim 10^{14}$ M$_{\odot}$. We find that a significant fraction of ICL stars ($8\%$-$28\%$) are born in-situ. This amounts to a total stellar mass comparable to the central galaxy itself. Contrary to simple expectations, only a sub-dominant fraction of these in-situ ICL stars are born in the central regions and later re-distributed to more energetic orbits during mergers. Instead, many in-situ ICL stars form directly hundreds of kiloparsecs away from the central galaxy, in clouds condensing out of the circum-cluster medium. The simulations predict a present-date diffuse star formation rate of $\sim$1 $\mathrm{M}_{\odot}$/yr, with higher rates at higher redshifts. The diffuse star forming component of the ICL is filamentary in nature, extends for hundreds of kiloparsecs and traces the distribution of neutral gas in the cluster host halo. We discuss briefly how numerical details of the baryonic treatment in the simulation may play a role in this result and conclude that a sensitivity of $1.6 \times 10^{-19} - 2.6 \times 10^{-18}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$ in H$_α$ flux (beyond current observational capabilities) would be necessary to detect this diffuse star-forming component in galaxy clusters. △ Less

Submitted 7 March, 2024; originally announced March 2024.

Comments: 10 pages, 8 figures - comments are welcome

The progenitors of the intra-cluster light and intra-cluster globular clusters in galaxy groups and clusters

Authors: Niusha Ahvazi, Laura V. Sales, Jessica E. Doppel, Andrew Benson, Richard D'Souza, Vicente Rodriguez-Gomez

Abstract: We use the IllustrisTNG50 cosmological hydrodynamical simulation, complemented by a catalog of tagged globular clusters, to investigate the properties and build up of two extended luminous components: the intra-cluster light (ICL) and the intra-cluster globular clusters (ICGC). We select the 39 most massive groups and clusters in the box, spanning the range of virial masses… ▽ More We use the IllustrisTNG50 cosmological hydrodynamical simulation, complemented by a catalog of tagged globular clusters, to investigate the properties and build up of two extended luminous components: the intra-cluster light (ICL) and the intra-cluster globular clusters (ICGC). We select the 39 most massive groups and clusters in the box, spanning the range of virial masses $5 \times 10^{12} < \rm M_{200}/\rm M_{\odot} < 2 \times 10^{14}$. We find good agreement between predictions from the simulations and current observational estimates of the fraction of mass in the ICL and its radial extension. The stellar mass of the ICL is only $\sim10\%-20\%$ of the stellar mass in the central galaxy but encodes useful information on the assembly history of the group or cluster. About half the ICL in all our systems is brought in by galaxies in a narrow stellar mass range, $M_*=10^{10}-10^{11}$ $\rm M_{\odot}$. However, the contribution of low-mass galaxies ($M_*<10^{10}$ $\rm M_{\odot}$) to the build-up of the ICL varies broadly from system to system, $\sim 5\%-45\%$, a feature that might be recovered from the observable properties of the ICL at $z=0$. At fixed virial mass, systems where the accretion of dwarf galaxies plays an important role have shallower metallicity profiles, less metal content and a lower stellar mass in the ICL than systems where the main contributors are more massive galaxies. We show that intra-cluster GCs are also good tracers of this history, representing a valuable alternative when diffuse light is not detectable. △ Less

Submitted 23 August, 2023; originally announced August 2023.

Constraining the assembly time of the stellar haloes of nearby Milky Way-mass galaxies through AGB populations

Authors: Benjamin Harmsen, Eric F. Bell, Richard D'Souza, Antonela Monachesi, Roelof S. de Jong, Adam Smercina, In Sung Jang, Benne W. Holwerda

Abstract: The star formation histories (SFHs) of galactic stellar haloes offer crucial insights into the merger history of the galaxy and the effects of those mergers on their hosts. Such measurements have revealed that while the Milky Way's most important merger was 8-10 Gyr ago, M31's largest merger was more recent, within the last few Gyr. Unfortunately, the required halo SFH measurements are extremely o… ▽ More The star formation histories (SFHs) of galactic stellar haloes offer crucial insights into the merger history of the galaxy and the effects of those mergers on their hosts. Such measurements have revealed that while the Milky Way's most important merger was 8-10 Gyr ago, M31's largest merger was more recent, within the last few Gyr. Unfortunately, the required halo SFH measurements are extremely observationally expensive outside of the Local Group. Here we use asymptotic giant branch (AGB) stars brighter than the tip of the red giant branch (RGB) to constrain stellar halo SFHs. Both stellar population models and archival datasets show that the AGB/RGB ratio constrains the time before which 90% of the stars formed, $t_{90}$. We find AGB stars in the haloes of three highly-inclined roughly Milky Way-mass galaxies with resolved star measurements from the Hubble Space Telescope; this population is most prominent in the stellar haloes of NGC 253 and NGC 891, suggesting that their stellar haloes contain stars born at relatively late times, with inferred $t_{90}\sim 6\pm1.5$Gyr. This ratio also varies from region to region, tending towards higher values along the major axis and in tidal streams or shells. By combining our measurements with previous constraints, we find a tentative anticorrelation between halo age and stellar halo mass, a trend that exists in models of galaxy formation but has never been elucidated before, i.e, the largest stellar haloes of Milky-Way mass galaxies were assembled more recently. △ Less

Submitted 22 August, 2023; originally announced August 2023.

Comments: 20 Pages, 10 Figures

Dragon kings in self-organized criticality systems

Authors: Guram Mikaberidze, Arthur Plaud, Raissa M. D'Souza

Abstract: The spontaneous emergence of scale invariance, called self-organized criticality (SOC), is often attributed to a second-order absorbing-state phase transition (ASPT). Many real-world systems display SOC, yet extreme events are often overrepresented, causing significant disruption, and are called dragon kings (DK). We show analytically that the tradeoff between driving impulse and dissipation rate… ▽ More The spontaneous emergence of scale invariance, called self-organized criticality (SOC), is often attributed to a second-order absorbing-state phase transition (ASPT). Many real-world systems display SOC, yet extreme events are often overrepresented, causing significant disruption, and are called dragon kings (DK). We show analytically that the tradeoff between driving impulse and dissipation rate can create DKs in a second-order ASPT. This establishes that DKs exist in SOC systems, reveals a taxonomy of DKs, and shows that larger dissipation and smoother driving lower risk of extreme events. △ Less

Submitted 4 August, 2023; originally announced August 2023.

Comments: 6 pages, 5 figures

Unified framework for hybrid percolation transitions based on microscopic dynamics

Authors: Hoyun Choi, Y. S. Cho, Raissa D'Souza, János Kertész, B. Kahng

Abstract: A hybrid percolation transition (HPT) exhibits both discontinuity of the order parameter and critical behavior at the transition point. Such dynamic transitions can occur in two ways: by cluster pruning with suppression of loop formation of cut links or by cluster merging with suppression of the creation of large clusters. While the microscopic mechanism of the former is understood in detail, a si… ▽ More A hybrid percolation transition (HPT) exhibits both discontinuity of the order parameter and critical behavior at the transition point. Such dynamic transitions can occur in two ways: by cluster pruning with suppression of loop formation of cut links or by cluster merging with suppression of the creation of large clusters. While the microscopic mechanism of the former is understood in detail, a similar framework is missing for the latter. By studying two distinct cluster merging models, we uncover the universal mechanism of the features of HPT-s at a microscopic level. We find that these features occur in three steps: (i) medium-sized clusters accumulate due to the suppression rule hindering the growth of large clusters, (ii) those medium size clusters eventually merge and a giant cluster increases rapidly, and (iii) the suppression effect becomes obsolete and the kinetics is governed by the Erdős-Rényi type of dynamics. We show that during the second and third period, the growth of the largest component must proceed in the form of a Devil's staircase. We characterize the critical behavior by two sets of exponents associated with the order parameter and cluster size distribution, which are related to each other by a scaling relation. Extensive numerical simulations are carried out to support the theory where a specific method is applied for finite-size scaling analysis to enable handling the large fluctuations of the transition point. Our results provide a unified theoretical framework for the HPT. △ Less

Submitted 7 July, 2023; originally announced July 2023.

Comments: 15 pages, 14 figures, 4 tables

Report number: 2307.03584

Origins of the Evil Eye: M64's Stellar Halo Reveals the Recent Accretion of an SMC-mass Satellite

Authors: Adam Smercina, Eric F. Bell, Paul A. Price, Jeremy Bailin, Julianne J. Dalcanton, Roelof S. de Jong, Richard D'Souza, Katya Gozman, In Sung Jang, Antonela Monachesi, David Nidever, Colin T. Slater

Abstract: M64, often called the "Evil Eye" galaxy, is unique among local galaxies. Beyond its dramatic, dusty nucleus, it also hosts an outer gas disk that counter-rotates relative to its stars. The mass of this outer disk is comparable to the gas content of the Small Magellanic Cloud (SMC), prompting the idea that it was likely accreted in a recent minor merger. Yet, detailed follow-up studies of M64's out… ▽ More M64, often called the "Evil Eye" galaxy, is unique among local galaxies. Beyond its dramatic, dusty nucleus, it also hosts an outer gas disk that counter-rotates relative to its stars. The mass of this outer disk is comparable to the gas content of the Small Magellanic Cloud (SMC), prompting the idea that it was likely accreted in a recent minor merger. Yet, detailed follow-up studies of M64's outer disk have shown no evidence of such an event, leading to other interpretations, such as a "flyby" interaction with the distant diffuse satellite Coma P. We present Subaru Hyper Suprime-Cam observations of M64's stellar halo, which resolve its stellar populations and reveal a spectacular radial shell feature, oriented $\sim$30$^{\circ}$ relative to the major axis and along the rotation axis of the outer gas disk. The shell is $\sim$45 kpc southeast of M64, while a similar but more diffuse plume to the northwest extends to $>$100 kpc. We estimate a stellar mass and metallicity for the southern shell of $M_{\star} {=} 1.80~{\pm}~0.54{\times}10^8~M_{\odot}$ and [M/H] $=$ $-$1.0, respectively, and a similar mass of $1.42~{\pm}~0.71{\times}10^8 M_{\odot}$ for the northern plume. Taking into account the accreted material in M64's inner disk, we estimate a total stellar mass for the progenitor satellite of $M_{\rm \star,prog}~{\simeq}~5{\times}10^8~M_{\odot}$. These results suggest that M64 is in the final stages of a minor merger with a gas-rich satellite strikingly similar to the SMC, in which M64's accreted counter-rotating gas originated, and which is responsible for the formation of its dusty inner star-forming disk. △ Less

Submitted 26 May, 2023; originally announced May 2023.

Comments: 12 pages, 5 figures, accepted for publication in ApJ Letters

Saying Hallo to M94's Stellar Halo: Investigating the Accretion History of the Largest Pseudobulge Host in the Local Universe

Authors: Katya Gozman, Eric F. Bell, Adam Smercina, Paul Price, Jeremy Bailin, Roelof S. de Jong, Richard D'Souza, In Sung Jang, Antonela Monachesi, Colin Slater

Abstract: It is not yet settled how the combination of secular processes and merging gives rise to the bulges and pseudobulges of galaxies. The nearby ($D\sim$ 4.2 Mpc) disk galaxy M94 (NGC 4736) has the largest pseudobulge in the local universe, and offers a unique opportunity for investigating the role of merging in the formation of its pseudobulge. We present a first ever look at M94's stellar halo, whic… ▽ More It is not yet settled how the combination of secular processes and merging gives rise to the bulges and pseudobulges of galaxies. The nearby ($D\sim$ 4.2 Mpc) disk galaxy M94 (NGC 4736) has the largest pseudobulge in the local universe, and offers a unique opportunity for investigating the role of merging in the formation of its pseudobulge. We present a first ever look at M94's stellar halo, which we expect to contain a fossil record of M94's past mergers. Using Subaru's Hyper Suprime-Cam, we resolve and identify red giant branch (RGB) stars in M94's halo, finding two distinct populations. After correcting for completeness through artificial star tests, we can measure the radial profile of each RGB population. The metal-rich RGB stars show an unbroken exponential profile to a radius of 30 kpc that is a clear continuation of M94's outer disk. M94's metal poor stellar halo is detectable over a wider area and clearly separates from its metal-rich disk. By integrating the halo density profile, we infer a total accreted stellar mass of $\sim 2.8 \times 10^8 M_\odot$, with a median metallicity of [M/H] $=-$1.4. This indicates that M94's most-massive past merger was with a galaxy similar to, or less massive than, the Small Magellanic Cloud. Few nearby galaxies have had such a low-mass dominant merger; therefore we suggest that M94's pseudobulge was not significantly impacted by merging. △ Less

Submitted 17 April, 2023; originally announced April 2023.

Comments: 19 pages, 12 figures, 1 table; Published in ApJ on April 14, 2023

Ultrafaint Dwarf Galaxy Candidates in the M81 Group: Signatures of Group Accretion

Authors: Eric F. Bell, Adam Smercina, Paul A. Price, Richard D'Souza, Jeremy Bailin, Roelof S. de Jong, Katya Gozman, In Sung Jang, Antonela Monachesi, Oleg Y. Gnedin, Colin T. Slater

Abstract: The faint and ultrafaint dwarf galaxies in the Local Group form the observational bedrock upon which our understanding of small-scale cosmology rests. In order to understand whether this insight generalizes, it is imperative to use resolved-star techniques to discover similarly faint satellites in nearby galaxy groups. We describe our search for ultrafaint galaxies in the M81 group using deep grou… ▽ More The faint and ultrafaint dwarf galaxies in the Local Group form the observational bedrock upon which our understanding of small-scale cosmology rests. In order to understand whether this insight generalizes, it is imperative to use resolved-star techniques to discover similarly faint satellites in nearby galaxy groups. We describe our search for ultrafaint galaxies in the M81 group using deep ground-based resolved-star data sets from Subaru's Hyper Suprime-Cam. We present one new ultrafaint dwarf galaxy in the M81 group and identify five additional extremely low surface brightness candidate ultrafaint dwarfs that reach deep into the ultrafaint regime to $M_V \sim -6$ (similar to current limits for Andromeda satellites). These candidates' luminosities and sizes are similar to known Local Group dwarf galaxies Tucana B, Canes Venatici I, Hercules, and Boötes I. Most of these candidates are likely to be real, based on tests of our techniques on blank fields. Intriguingly, all of these candidates are spatially clustered around NGC 3077, which is itself an M81 group satellite in an advanced state of tidal disruption. This is somewhat surprising, as M81 itself and its largest satellite M82 are both substantially more massive than NGC 3077 and by virtue of their greater masses, would have been expected to host as many or more ultrafaint candidates. These results lend considerable support to the idea that satellites of satellites are an important contribution to the growth of satellite populations around Milky Way-mass galaxies. △ Less

Submitted 13 September, 2022; originally announced September 2022.

Comments: The Astrophysical Journal Letters; in press. 11 pages, 4 figures, 1 table

New Globular Cluster Candidates in the M81 group

Authors: Jiaming Pan, Eric F. Bell, Adam Smercina, Paul Price, Colin T. Slater, Jeremy Bailin, Roelof S. de Jong, Richard D'Souza, In Sung Jang, Antonela Monachesi

Abstract: The study of outer halo globular cluster (GC) populations can give insight into galaxy merging, globular cluster accretion and the origin of GCs. We use archival Subaru Hyper Suprime-Cam (HSC) data in concert with space-based GALEX, IRAC and Gaia EDR3 data to select candidate Globular clusters (GCs) in the outer halo of the M81 group for confirmation and future study. We use a small sample of prev… ▽ More The study of outer halo globular cluster (GC) populations can give insight into galaxy merging, globular cluster accretion and the origin of GCs. We use archival Subaru Hyper Suprime-Cam (HSC) data in concert with space-based GALEX, IRAC and Gaia EDR3 data to select candidate Globular clusters (GCs) in the outer halo of the M81 group for confirmation and future study. We use a small sample of previously-discovered GCs to tune our selection criteria, finding that bright already-known GCs in the M81 group have sizes that are typically slightly larger than the Subaru PSF in our fields. In the optical bands, GCs appear to have colours that are only slightly different from stars. The inclusion of archival IRAC data yields dramatic improvements in colour separation, as the long wavelength baseline aids somewhat in the separation from stars and clearly separates GCs from many compact background galaxies. We show that some previously spectroscopically-identified GCs in the M81 group are instead foreground stars or background galaxies. GCs close to M82 have radial velocities suggesting that they fell into the M81 group along with M82. The overall M81 GC luminosity function is similar to the Milky Way and M31. M81's outer halo GCs are similar to the Milky Way in their metallicities and numbers, and much less numerous than M31's more metal-rich outer halo GC population. These properties reflect differences in the three galaxies' merger histories, highlighting the possibility of using outer halo GCs to trace merger history in larger samples of galaxies. △ Less

Submitted 14 June, 2022; originally announced June 2022.

Temperature induced band convergence, intervalley scattering and thermoelectric transport in p-type PbTe

Authors: Ransell D'Souza, José D. Querales-Flores, Jiang Cao, Stephen Fahy, Ivana Savić

Abstract: Achieving high valley degeneracy (i.e. "band convergence") in a material usually results in considerably enhanced thermoelectric properties. However, it is still unclear why this strategy of designing efficient thermoelectric materials is so successful, since the benefit of increased density of states may be severely degraded by intervalley scattering. Using first-principles calculations, we inves… ▽ More Achieving high valley degeneracy (i.e. "band convergence") in a material usually results in considerably enhanced thermoelectric properties. However, it is still unclear why this strategy of designing efficient thermoelectric materials is so successful, since the benefit of increased density of states may be severely degraded by intervalley scattering. Using first-principles calculations, we investigate these effects in $p$-type PbTe, where temperature induces alignment of the $L$ and $Σ$ valleys at $\sim$~620~K. We explicitly show that the thermoelectric power factor and figure of merit peak near the band convergence temperature. The figure of merit maximum is larger than those of the individual $L$ and $Σ$ valleys. Surprisingly, intervalley scattering does not considerably affect the figure of merit near the band convergence temperature and optimal doping conditions, although it reduces the power factor by almost a factor of 2. Our results suggest that band convergence will significantly increase the figure of merit if intervalley scattering is roughly proportional to the density of states and the lattice thermal conductivity is considerably lower than the electronic thermal conductivity, even if intervalley scattering is strong. △ Less

Submitted 27 May, 2022; originally announced May 2022.

Comments: 15 pages, 6 figures, Journal publication

The Global Dynamical Atlas of the Milky Way mergers: Constraints from Gaia EDR3 based orbits of globular clusters, stellar streams and satellite galaxies

Authors: Khyati Malhan, Rodrigo A. Ibata, Sanjib Sharma, Benoit Famaey, Michele Bellazzini, Raymond G. Carlberg, Richard D'Souza, Zhen Yuan, Nicolas F. Martin, Guillaume F. Thomas

Abstract: The Milky Way halo was predominantly formed by the merging of numerous progenitor galaxies. However, our knowledge of this process is still incomplete, especially in regard to the total number of mergers, their global dynamical properties and their contribution to the stellar population of the Galactic halo. Here, we uncover the Milky Way mergers by detecting groupings of globular clusters, stella… ▽ More The Milky Way halo was predominantly formed by the merging of numerous progenitor galaxies. However, our knowledge of this process is still incomplete, especially in regard to the total number of mergers, their global dynamical properties and their contribution to the stellar population of the Galactic halo. Here, we uncover the Milky Way mergers by detecting groupings of globular clusters, stellar streams and satellite galaxies in action ($\mathbf{J}$) space. While actions fully characterize the orbits, we additionally use the redundant information on their energy ($\textit{E}$) to enhance the contrast between groupings. For this endeavour, we use $\textit{Gaia}$ EDR3 based measurements of $170$ globular clusters, $41$ streams and $46$ satellites to derive their $\mathbf{J}$ and $\textit{E}$. To detect groups, we use the $\texttt{ENLINK}$ software, coupled with a statistical procedure that accounts for the observed phase-space uncertainties of these objects. We detect a total of $N=6$ groups, including the previously known mergers $\textit{Sagittarius}$, $\textit{Cetus}$, $\textit{Gaia-Sausage/Enceladus}$, $\textit{LMS-1/Wukong}$, $\textit{Arjuna/Sequoia/I'itoi}$ and one new merger that we call $\textit{Pontus}$. All of these mergers, together, comprise $62$ objects ($\approx 25\%$ of our sample). We discuss their members, orbital properties and metallicity distributions. We find that the three most metal-poor streams of our Galaxy -- "C-19" ([Fe/H]$=-3.4$ dex), "Sylgr" ([Fe/H]$=-2.9$ dex) and "Phoenix" ([Fe/H]$=-2.7$ dex) -- are associated with $\textit{LMS-1/Wukong}$; showing it to be the most metal-poor merger. The global dynamical atlas of Milky Way mergers that we present here provides a present-day reference for galaxy formation models. △ Less

Submitted 15 February, 2022; originally announced February 2022.

Comments: Published in ApJ. Key results in Table 3. The paper contains 38 pages, 17 figures, 4 Tables

Uncertainties associated with the backward integration of dwarf satellites using simple parametric potentials

Authors: Richard D'Souza, Eric F. Bell

Abstract: In order to backward integrate the orbits of Milky Way (MW) dwarf galaxies, much effort has been invested in recent years to constrain their initial phase-space coordinates. Yet equally important are the assumptions on the potential that the dwarf galaxies experience over time, especially given the fact that the MW is currently accreting the Large Magellanic Cloud (LMC). In this work, using a dark… ▽ More In order to backward integrate the orbits of Milky Way (MW) dwarf galaxies, much effort has been invested in recent years to constrain their initial phase-space coordinates. Yet equally important are the assumptions on the potential that the dwarf galaxies experience over time, especially given the fact that the MW is currently accreting the Large Magellanic Cloud (LMC). In this work, using a dark matter-only zoom-in simulation, we test whether the use of common parametric forms of the potential is adequate to successfully backward integrate the orbits of the subhaloes from their present-day positions. We parametrise the recovered orbits and compare them with those from the simulations. We find that simple symmetric parametric forms of the potential fail to capture the complexities and the inhomogeneities of the true potential experienced by the subhaloes. More specifically, modelling a recent massive accretion like that of the LMC as a sum of two spherical parametric potentials leads to substantial errors in the recovered parameters of the orbits. These errors rival those caused due to a) a 30\% uncertainty in the virial mass of the MW and b) not modelling the potential of the recently accreted massive satellite. Our work suggests that i) the uncertainties in the parameters of the recovered orbits of some MW dwarfs may be under-estimated and that ii) researchers should characterise the uncertainties inherent to their choice of integration techniques and assumptions of the potential against cosmological zoom-in simulations of the MW, which include a recently-accreted LMC. △ Less

Submitted 11 February, 2022; originally announced February 2022.

Comments: Accepted for publication in MNRAS

Sandpile cascades on oscillator networks: the BTW model meets Kuramoto

Authors: Guram Mikaberidze, Raissa M. D'Souza

Abstract: Cascading failures abound in complex systems and the BTW sandpile model provides a theoretical underpinning for their analysis. Yet, it does not account for the possibility of nodes having oscillatory dynamics such as in power grids and brain networks. Here we consider a network of Kuramoto oscillators upon which the BTW model is unfolding, enabling us to study how the feedback between the oscilla… ▽ More Cascading failures abound in complex systems and the BTW sandpile model provides a theoretical underpinning for their analysis. Yet, it does not account for the possibility of nodes having oscillatory dynamics such as in power grids and brain networks. Here we consider a network of Kuramoto oscillators upon which the BTW model is unfolding, enabling us to study how the feedback between the oscillatory and cascading dynamics can lead to new emergent behaviors. We assume that the more out-of-sync a node is with its neighbors the more vulnerable it is and lower its load-carrying capacity accordingly. And when a node topples and sheds load, its oscillatory phase is reset at random. This leads to novel cyclic behavior at an emergent, long timescale. The system spends the bulk of its time in a synchronized state where load builds up with minimal cascades. Yet, eventually the system reaches a tipping point where a large cascade triggers a "cascade of larger cascades," which can be classified as a Dragon King event. The system then undergoes a short transient back to the synchronous, build-up phase. The coupling between capacity and synchronization gives rise to endogenous cascade seeds in addition to the standard exogenous ones, and we show their respective roles. We establish the phenomena from numerical studies and develop the accompanying mean-field theory to locate the tipping point, calculate the load in the system, determine the frequency of the long-time oscillations and find the distribution of cascade sizes during the build-up phase. △ Less

Submitted 28 March, 2022; v1 submitted 30 November, 2021; originally announced December 2021.

Comments: 13 pages, 11 figures

Journal ref: Chaos 32, 053121 (2022)

An Algorithm for Local Transverse Feedback Linearization

Authors: Rollen S. D'Souza, Christopher Nielsen

Abstract: Given a multi-input, nonlinear, time-invariant, control-affine system and a controlled invariant, closed, embedded submanifold $\mathsf{N}$, the local transverse feedback linearization (TFL) problem seeks a coordinate and feedback transformation such that, in transformed coordinates, the dynamics governing the system's transverse evolution with respect to $\mathsf{N}$ are linear, time-invariant an… ▽ More Given a multi-input, nonlinear, time-invariant, control-affine system and a controlled invariant, closed, embedded submanifold $\mathsf{N}$, the local transverse feedback linearization (TFL) problem seeks a coordinate and feedback transformation such that, in transformed coordinates, the dynamics governing the system's transverse evolution with respect to $\mathsf{N}$ are linear, time-invariant and controllable. The transformed system is said to be in the TFL normal form. Checkable necessary and sufficient conditions for this problem to be solvable are known, but, unfortunately, the literature does not present a prescription that constructs the required transformation for multi-input systems. In this article we present an algorithm that produces a virtual output of suitable vector relative degree that, using input-output feedback linearization, puts the system into TFL normal form. The procedure is based on dual conditions for TFL and is fundamentally different from existing methods, such as the GS and Blended algorithms, because of the "desired" zero dynamics manifold $\mathsf{N}$. The proposed algorithm is the first to take into consideration the desired zero dynamics manifold. △ Less

Submitted 28 April, 2022; v1 submitted 1 September, 2021; originally announced September 2021.

Comments: 25 pages, 2 figures, 1 algorithm. Resubmitted to SIAM Journal on Control and Optimization

MSC Class: 37N35; 93B10; 93B27

Analyzing states beyond full synchronization on hypergraphs requires methods beyond projected networks

Authors: Anastasiya Salova, Raissa M. D'Souza

Abstract: A common approach for analyzing hypergraphs is to consider the projected adjacency or Laplacian matrices for each order of interactions (e.g., dyadic, triadic, etc.). However, this method can lose information about the hypergraph structure and is not universally applicable for studying dynamical processes on hypergraphs, which we demonstrate through the framework of cluster synchronization. Specif… ▽ More A common approach for analyzing hypergraphs is to consider the projected adjacency or Laplacian matrices for each order of interactions (e.g., dyadic, triadic, etc.). However, this method can lose information about the hypergraph structure and is not universally applicable for studying dynamical processes on hypergraphs, which we demonstrate through the framework of cluster synchronization. Specifically, we show that the projected network does not always correspond to a unique hypergraph structure. This means the projection does not always properly predict the true dynamics unfolding on the hypergraph. Additionally, we show that the symmetry group consisting of permutations that preserve the hypergraph structure can be distinct from the symmetry group of its projected matrix. Thus, considering the full hypergraph is required for analyzing the most general types of dynamics on hypergraphs. We show that a formulation based on node clusters and the corresponding edge clusters induced by the node partitioning, enables the analysis of admissible patterns of cluster synchronization and their effective dynamics. Additionally, we show that the coupling matrix projections corresponding to each edge cluster synchronization pattern, and not just to each order of interactions, are necessary for understanding the structure of the Jacobian matrix and performing the linear stability calculations efficiently. △ Less

Submitted 28 July, 2021; originally announced July 2021.

Comments: 17 pages, 11 figures

Relating the Diverse Merger Histories and Satellite Populations of Nearby Galaxies

Authors: Adam Smercina, Eric F. Bell, Jenna Samuel, Richard D'Souza

Abstract: We investigate whether the considerable diversity in the satellite populations of nearby Milky Way (MW)-mass galaxies is connected with the diversity in their host's merger histories. Analyzing 8 nearby galaxies with extensive observations of their satellite populations and stellar halos, we characterize each galaxy's merger history using the metric of its most dominant merger,… ▽ More We investigate whether the considerable diversity in the satellite populations of nearby Milky Way (MW)-mass galaxies is connected with the diversity in their host's merger histories. Analyzing 8 nearby galaxies with extensive observations of their satellite populations and stellar halos, we characterize each galaxy's merger history using the metric of its most dominant merger, $M_{\rm \star,Dom}$, defined as the greater of either its total accreted stellar mass or most massive current satellite. We find an unexpectedly tight relationship between these galaxies' number of $M_{V}\,{<}\,{-}9$ satellites within 150 kpc ($N_{\rm Sat}$) and $M_{\rm \star,Dom}$. This relationship remains even after accounting for differences in galaxy mass. Using the star formation and orbital histories of satellites around the MW and M81, we demonstrate that both likely evolved along the $M_{\rm\star,Dom}{-}N_{\rm Sat}$ relation during their current dominant mergers with the LMC and M82, respectively. We investigate the presence of this relation in galaxy formation models, including using the FIRE simulations to directly compare to the observations. We find no relation between $M_{\rm\star,Dom}$ and $N_{\rm Sat}$ in FIRE, and a universally large scatter in $N_{\rm Sat}$ with $M_{\rm \star,Dom}$ across models $-$ in direct contrast with the tightness of the empirical relation. This acute difference in the observed and predicted scaling relation between two fundamental galaxy properties signals that current simulations do not sufficiently reproduce diverse merger histories and their effects on satellite populations. Explaining the emergence of this relation is therefore essential for obtaining a complete understanding of galaxy formation. △ Less

Submitted 10 March, 2022; v1 submitted 9 July, 2021; originally announced July 2021.

Comments: 26 pages, 8 figures, 3 tables; Accepted for publication in ApJ

The infall of dwarf satellite galaxies are influenced by their host's massive accretions

Authors: Richard D'Souza, Eric F. Bell

Abstract: Recent progress in constraining the massive accretions (>1:10) experienced by the Milky Way (MW) and the Andromeda galaxy (M31) offers an opportunity to understand the dwarf galaxy population of the Local Group. Using zoom-in dark matter-only simulations of MW-mass haloes and concentrating on subhaloes that are thought to be capable of hosting dwarf galaxies, we demonstrate that the infall of a ma… ▽ More Recent progress in constraining the massive accretions (>1:10) experienced by the Milky Way (MW) and the Andromeda galaxy (M31) offers an opportunity to understand the dwarf galaxy population of the Local Group. Using zoom-in dark matter-only simulations of MW-mass haloes and concentrating on subhaloes that are thought to be capable of hosting dwarf galaxies, we demonstrate that the infall of a massive progenitor is accompanied with the accretion and destruction of a large number of subhaloes. Massive accreted progenitors do not increase the total number of infalling subhaloes onto a MW-mass host, but instead focus surrounding subhaloes onto the host causing a clustering in the infall time of subhaloes. This leads to a temporary elevation in the number of subhaloes as well as changes in their cumulative radial profile within the virial radius of the host. Surviving associated subhaloes with a massive progenitor have a large diversity in their orbits. We find that the star formation quenching times of Local Group dwarf spheroidal galaxies ($10^{5} \mathrm{M_{\odot}} \lesssim \mathrm{M}_{*} \lesssim 10^{7} \mathrm{M_{\odot}}$) are clustered around the times of the most massive accretions suffered by the MW and M31. Our results imply that a) the quenching time of dwarf spheroidals is a good proxy of their infall time and b) the absence of recently quenched satellites around M31 suggests that M33 is not on its first infall and was accreted much earlier. △ Less

Submitted 27 April, 2021; originally announced April 2021.

Comments: Accepted for publication in MNRAS. Accompanying movie available at https://radsouza.github.io/simulations.html

Uncovering near-wall blood flow from sparse data with physics-informed neural networks

Authors: Amirhossein Arzani, Jian-Xun Wang, Roshan M. D'Souza

Abstract: Near-wall blood flow and wall shear stress (WSS) regulate major forms of cardiovascular disease, yet they are challenging to quantify with high fidelity. Patient-specific computational and experimental measurement of WSS suffers from uncertainty, low resolution, and noise issues. Physics-informed neural networks (PINN) provide a flexible deep learning framework to integrate mathematical equations… ▽ More Near-wall blood flow and wall shear stress (WSS) regulate major forms of cardiovascular disease, yet they are challenging to quantify with high fidelity. Patient-specific computational and experimental measurement of WSS suffers from uncertainty, low resolution, and noise issues. Physics-informed neural networks (PINN) provide a flexible deep learning framework to integrate mathematical equations governing blood flow with measurement data. By leveraging knowledge about the governing equations (herein, Navier-Stokes), PINN overcomes the large data requirement in deep learning. In this study, it was shown how PINN could be used to improve WSS quantification in diseased arterial flows. Specifically, blood flow problems where the inlet and outlet boundary conditions were not known were solved by assimilating very few measurement points. Uncertainty in boundary conditions is a common feature in patient-specific computational fluid dynamics models. It was shown that PINN could use sparse velocity measurements away from the wall to quantify WSS with very high accuracy even without full knowledge of the boundary conditions. Examples in idealized stenosis and aneurysm models were considered demonstrating how partial knowledge about the flow physics could be combined with partial measurements to obtain accurate near-wall blood flow data. The proposed hybrid data-driven and physics-based deep learning framework has high potential in transforming high-fidelity near-wall hemodynamics modeling in cardiovascular disease. △ Less

Submitted 7 June, 2021; v1 submitted 16 April, 2021; originally announced April 2021.

Integrating multi-fidelity blood flow data with reduced-order data assimilation

Authors: Milad Habibi, Roshan M. D'Souza, Scott T. M. Dawson, Amirhossein Arzani

Abstract: High-fidelity patient-specific modeling of cardiovascular flows and hemodynamics is challenging. Direct blood flow measurement inside the body with in-vivo measurement modalities such as 4D flow magnetic resonance imaging (4D flow MRI) suffer from low resolution and acquisition noise. In-vitro experimental modeling and patient-specific computational fluid dynamics (CFD) models are subject to uncer… ▽ More High-fidelity patient-specific modeling of cardiovascular flows and hemodynamics is challenging. Direct blood flow measurement inside the body with in-vivo measurement modalities such as 4D flow magnetic resonance imaging (4D flow MRI) suffer from low resolution and acquisition noise. In-vitro experimental modeling and patient-specific computational fluid dynamics (CFD) models are subject to uncertainty in patient-specific boundary conditions and model parameters. Furthermore, collecting blood flow data in the near-wall region (e.g., wall shear stress) with experimental measurement modalities poses additional challenges. In this study, a computationally efficient data assimilation method called reduced-order modeling Kalman filter (ROM-KF) was proposed, which combined a sequential Kalman filter with reduced-order modeling using a linear model provided by dynamic mode decomposition (DMD). The goal of ROM-KF was to overcome low resolution and noise in experimental and uncertainty in CFD modeling of cardiovascular flows. The accuracy of the method was assessed with 1D Womersley flow, 2D idealized aneurysm, and 3D patient-specific cerebral aneurysm models. Synthetic experimental data were used to enable direct quantification of errors using benchmark datasets. The accuracy of ROM-KF in reconstructing near-wall hemodynamics was assessed by applying the method to problems where near-wall blood flow data were missing in the experimental dataset. The ROM-KF method provided blood flow data that were more accurate than the computational and synthetic experimental datasets and improved near-wall hemodynamics quantification. △ Less

Submitted 9 June, 2021; v1 submitted 5 April, 2021; originally announced April 2021.

Quantifying the Global Support Network for Non-State Armed Groups (NAGs)

Authors: Weiran Cai, Belgin San-Akca, Jordan Snyder, Grayson Gordon, Zeev Maoz, Raissa M. D'Souza

Abstract: Human history has been shaped by armed conflicts. Rather than large-scale interstate wars, low-intensity attacks have been more prevalent in the post-World War era. These attacks are often carried out by non-state armed groups (NAGs), which are supported by host states (HSs). We analyze the global bipartite network of NAG-HS support and its evolution over the period of 1945-2010. We find striking… ▽ More Human history has been shaped by armed conflicts. Rather than large-scale interstate wars, low-intensity attacks have been more prevalent in the post-World War era. These attacks are often carried out by non-state armed groups (NAGs), which are supported by host states (HSs). We analyze the global bipartite network of NAG-HS support and its evolution over the period of 1945-2010. We find striking parallels to ecological networks such as mutualistic and parasitic forms of support, and a nested and modular network architecture. The nestedness emerges from preferential behaviors: highly connected players are more likely to both gain and lose connections. Long-persisting major modules are identified, reflecting both regional and trans-regional interests, which show significant turnover in their membership, contrary to the transitory ones. Revealing this architecture further enables the identification of actor's roles and provide insights for effective intervention strategies. △ Less

Submitted 31 January, 2021; originally announced February 2021.

Comments: 13 pages, 6 figures

Cluster synchronization on hypergraphs

Authors: Anastasiya Salova, Raissa M. D'Souza

Abstract: Full synchronization of dynamical elements coupled via hypergraphs can be analyzed with the hypergraph projection onto dyadic matrices, but this is not sufficient for analyzing cluster synchronization. Here we develop the necessary formalism. We introduce the notion of edge clusters and show how node and edge partitions allow us to verify admissible states and simplify their linear stability calcu… ▽ More Full synchronization of dynamical elements coupled via hypergraphs can be analyzed with the hypergraph projection onto dyadic matrices, but this is not sufficient for analyzing cluster synchronization. Here we develop the necessary formalism. We introduce the notion of edge clusters and show how node and edge partitions allow us to verify admissible states and simplify their linear stability calculations. This provides a principled way to track dynamics on hypergraphs, and the projected Laplacian matrices based on each edge cluster are essential to linear stability analysis and its dimensionality reduction. This work goes beyond full synchronization and beyond dyadic interactions. △ Less

Submitted 17 March, 2022; v1 submitted 14 January, 2021; originally announced January 2021.

Comments: 11 pages, 5 figures; updated content, figures, and references

Origins of Algorithmic Instabilities in Crowdsourced Ranking

Authors: Keith Burghardt, Tad Hogg, Raissa M. D'Souza, Kristina Lerman, Marton Posfai

Abstract: Crowdsourcing systems aggregate decisions of many people to help users quickly identify high-quality options, such as the best answers to questions or interesting news stories. A long-standing issue in crowdsourcing is how option quality and human judgement heuristics interact to affect collective outcomes, such as the perceived popularity of options. We address this limitation by conducting a con… ▽ More Crowdsourcing systems aggregate decisions of many people to help users quickly identify high-quality options, such as the best answers to questions or interesting news stories. A long-standing issue in crowdsourcing is how option quality and human judgement heuristics interact to affect collective outcomes, such as the perceived popularity of options. We address this limitation by conducting a controlled experiment where subjects choose between two ranked options whose quality can be independently varied. We use this data to construct a model that quantifies how judgement heuristics and option quality combine when deciding between two options. The model reveals popularity-ranking can be unstable: unless the quality difference between the two options is sufficiently high, the higher quality option is not guaranteed to be eventually ranked on top. To rectify this instability, we create an algorithm that accounts for judgement heuristics to infer the best option and rank it first. This algorithm is guaranteed to be optimal if data matches the model. When the data does not match the model, however, simulations show that in practice this algorithm performs better or at least as well as popularity-based and recency-based ranking for any two-choice question. Our work suggests that algorithms relying on inference of mathematical models of user behavior can substantially improve outcomes in crowdsourcing systems. △ Less

Submitted 23 October, 2020; originally announced October 2020.

Comments: 12 pages, 20 figures

Journal ref: Proc. ACM Hum.-Comput. Interact., Vol. 4, No. CSCW2, Article 166. 2020

Network and Phase Symmetries Reveal That Amplitude Dynamics Stabilize Decoupled Oscillator Clusters

Authors: J. Emenheiser, A. Salova, J. Snyder, J. P. Crutchfield, R. M. D'Souza

Abstract: Oscillator networks display intricate synchronization patterns. Determining their stability typically requires incorporating the symmetries of the network coupling. Going beyond analyses that appeal only to a network's automorphism group, we explore synchronization patterns that emerge from the phase-shift invariance of the dynamical equations and symmetries in the nodes. We show that these nonstr… ▽ More Oscillator networks display intricate synchronization patterns. Determining their stability typically requires incorporating the symmetries of the network coupling. Going beyond analyses that appeal only to a network's automorphism group, we explore synchronization patterns that emerge from the phase-shift invariance of the dynamical equations and symmetries in the nodes. We show that these nonstructural symmetries simplify stability calculations. We analyze a ring-network of phase-amplitude oscillators that exhibits a "decoupled" state in which physically-coupled nodes appear to act independently due to emergent cancellations in the equations of dynamical evolution. We establish that this state can be linearly stable for a ring of phase-amplitude oscillators, but not for a ring of phase-only oscillators that otherwise require explicit long-range, nonpairwise, or nonphase coupling. In short, amplitude-phase interactions are key to stable synchronization at a distance. △ Less

Submitted 10 December, 2020; v1 submitted 18 October, 2020; originally announced October 2020.

Comments: 8 pages, 4 figures; http://csc.ucdavis.edu/~cmg/compmech/pubs/gamma_sync.htm

Electron-phonon scattering and thermoelectric transport in $p$-type PbTe from first principles

Authors: Ransell D'Souza, Jiang Cao, José D. Querales-Flores, Stephen Fahy, Ivana Savić

Abstract: We present a first principles based model of electron-phonon scattering mechanisms and thermoelectric transport at the L and $Σ$ valleys in <200b>$p$-type PbTe, accounting for their thermally induced shifts. Our calculated values of all thermoelectric transport parameters at room temperature are in very good agreement with experiments for a wide range of doping concentrations. Scattering due to lo… ▽ More We present a first principles based model of electron-phonon scattering mechanisms and thermoelectric transport at the L and $Σ$ valleys in <200b>$p$-type PbTe, accounting for their thermally induced shifts. Our calculated values of all thermoelectric transport parameters at room temperature are in very good agreement with experiments for a wide range of doping concentrations. Scattering due to longitudinal optical phonons is the main scattering mechanism in $p$-type PbTe, while scattering due to transverse optical modes is the weakest. The L valleys contribute most to thermoelectric transport at 300 K due to the sizeable energy difference between the L and $Σ$ valleys. We show that both scattering between the L and $Σ$ valleys and additional transport channels of the $Σ$ valleys are beneficial for the overall thermoelectric performance of $p$-type PbTe at 300 K. Our findings thus support the idea that materials with high valley degeneracy may be good thermoelectrics. △ Less

Submitted 24 August, 2020; originally announced August 2020.

Comments: 16 pages, 6 figures

Journal ref: Phys. Rev. B 102, 115204 (2020)

Is NGC 300 a pure exponential disk galaxy?

Authors: In Sung Jang, Roelof S. de Jong, I. Minchev, Eric F. Bell, Antonela Monachesi, Benne W. Holwerda, Jeremy Bailin, Adam Smercina, Richard D'Souza

Abstract: NGC 300 is a low-mass disk galaxy in the Sculptor group. In the literature, it has been identified as a pure exponential disk galaxy, as its luminosity profile could be well fitted with a single exponential law over many disk scale lengths (Type I). We investigate the stellar luminosity distribution of NGC 300 using $Hubble$ $Space$ $Telescope$ (HST) archive data, reaching farther and deeper than… ▽ More NGC 300 is a low-mass disk galaxy in the Sculptor group. In the literature, it has been identified as a pure exponential disk galaxy, as its luminosity profile could be well fitted with a single exponential law over many disk scale lengths (Type I). We investigate the stellar luminosity distribution of NGC 300 using $Hubble$ $Space$ $Telescope$ (HST) archive data, reaching farther and deeper than any other previous studies. Color magnitude diagrams show a significant population of old red giant branch (RGB) stars in all fields out to $R\sim19$ kpc ($32'$), as well as younger populations in the inner regions. We construct the density profiles of the young, intermediate-aged, and old stellar populations. We find two clear breaks in the density profiles of the old RGB and intermediate-aged stars: one down-bending (Type II) at $R\sim5.9$ kpc, and another up-bending (Type III) at $R\sim8.3$ kpc. Moreover, the old RGB stars exhibit a negative radial color gradient with an up-bending at $R\sim8$~kpc, beyond which the stellar populations are uniformly old ($>$7~Gyr) and metal-poor ($\rm[Fe/H] = -1.6^{+0.2}_{-0.4}$ dex). The outer stellar component at $R\gtrapprox8$ kpc is, therefore, well separated from the inner disk in terms of the stellar density and stellar populations. While our results cast doubt on the currently established wisdom that NGC\,300 is a pure exponential disk galaxy, a more detailed survey should be carried out to identify the outskirts as either a disk or a stellar halo. △ Less

Submitted 27 July, 2020; originally announced July 2020.

Comments: 6 pages, 5 figures, 1 appendix, accepted for publication in A&A Letters

Decoupled synchronized states in networks of linearly coupled limit cycle oscillators

Authors: Anastasiya Salova, Raissa M. D'Souza

Abstract: Networks of limit cycle oscillators can show intricate patterns of synchronization such as splay states and cluster synchronization. Here we analyze dynamical states that display a continuum of seemingly independent splay clusters. Each splay cluster is a block splay state consisting of sub-clusters of fully synchronized nodes with uniform amplitudes. Phases of nodes within a splay cluster are equ… ▽ More Networks of limit cycle oscillators can show intricate patterns of synchronization such as splay states and cluster synchronization. Here we analyze dynamical states that display a continuum of seemingly independent splay clusters. Each splay cluster is a block splay state consisting of sub-clusters of fully synchronized nodes with uniform amplitudes. Phases of nodes within a splay cluster are equally spaced, but nodes in different splay clusters have an arbitrary phase difference that can be fixed or evolve linearly in time. Such coexisting splay clusters form a decoupled state in that the dynamical equations become effectively decoupled between oscillators that can be physically coupled. We provide the conditions that allow the existence of particular decoupled states by using the eigendecomposition of the coupling matrix. Additionally, we provide an algorithm to search for admissible decoupled states using the external equitable partition and orbital partition considerations combined with symmetry groupoid formalism. Unlike previous studies, our approach is applicable when existence does not follow from symmetries alone and also illustrates the differences between adjacency and Laplacian coupling. We show that the decoupled state can be linearly stable for a substantial range of parameters using a simple eight-node cube network and its modifications as an example. We also demonstrate how the linear stability analysis of decoupled states can be simplified by taking into account the symmetries of the Jacobian matrix. Some network structures can support multiple decoupled patterns. To illustrate that, we show the variety of qualitatively different decoupled states that can arise on two-dimensional square and hexagonal lattices. △ Less

Submitted 30 June, 2020; v1 submitted 10 June, 2020; originally announced June 2020.

Comments: 19 pages, 11 figures

Journal ref: Phys. Rev. Research 2, 043261 (2020)

Correlated structural evolution within multiplex networks

Authors: Haochen Wu, Ryan G. James, James P. Crutchfield, Raissa M. D'Souza

Abstract: Many natural, engineered, and social systems can be represented using the framework of a layered network, where each layer captures a different type of interaction between the same set of nodes. The study of such multiplex networks is a vibrant area of research. Yet, understanding how to quantify the correlations present between pairs of layers, and more so present in their co-evolution, is lackin… ▽ More Many natural, engineered, and social systems can be represented using the framework of a layered network, where each layer captures a different type of interaction between the same set of nodes. The study of such multiplex networks is a vibrant area of research. Yet, understanding how to quantify the correlations present between pairs of layers, and more so present in their co-evolution, is lacking. Such methods would enable us to address fundamental questions involving issues such as function, redundancy and potential disruptions. Here we show first how the edge-set of a multiplex network can be used to construct an estimator of a joint probability distribution describing edge existence over all layers. We then adapt an information-theoretic measure of general correlation called the conditional mutual information, which uses the estimated joint probability distribution, to quantify the pairwise correlations present between layers. The pairwise comparisons can also be temporal, allowing us to identify if knowledge of a certain layer can provide additional information about the evolution of another layer. We analyze datasets from three distinct domains---economic, political, and airline networks---to demonstrate how pairwise correlation in structure and dynamical evolution between layers can be identified and show that anomalies can serve as potential indicators of major events such as shocks. △ Less

Submitted 9 May, 2020; originally announced May 2020.

Comments: Accepted to Journal of Complex Networks

Degree-targeted cascades in modular, degree-heterogeneous networks

Authors: Jordan Snyder, Weiran Cai, Raissa M. D'Souza

Abstract: The dynamics of cascading activation, such as rapid changes in public opinion and the outbreak of disease epidemics, have a crucial dependence on the connectivity patterns among the agents. We study cascading dynamics in modular, degree-heterogeneous networks, and consider the impact of intra-module seeding strategy on inter-module spread. Specifically, we establish that although activating the hi… ▽ More The dynamics of cascading activation, such as rapid changes in public opinion and the outbreak of disease epidemics, have a crucial dependence on the connectivity patterns among the agents. We study cascading dynamics in modular, degree-heterogeneous networks, and consider the impact of intra-module seeding strategy on inter-module spread. Specifically, we establish that although activating the highest-degree nodes is more effective than random selection at growing a cascade locally, there is a critical level of inter-module connectivity required for a cascade to cross from one module to another, irrespective of the seeding strategy. We present an analytical proof of this statement for the case that each module has the same degree distribution and all module pairs have the same inter-module connectivity, while our simulation results suggest its validity for more general situations, including a ring of modules. Interestingly, we find that on a network comprised of two modules, this critical level is primarily determined by the degree distribution of the \emph{alter} module, as opposed to the seed module. Our analytical approach extends a method developed by Gleeson, but is able to capture different seeding strategies using only one dynamical variable per module, namely the conditional exposure probability. Our work shows that the possibility of a global cascade depends sensitively on inter-module connectivity, and less on the intra-module seeding strategy. This suggests, for example, that slight changes to inter-module connectivity can be a feasible intervention strategy to promote or inhibit global cascades. △ Less

Submitted 11 January, 2022; v1 submitted 16 April, 2020; originally announced April 2020.

Strain induced large enhancement of thermoelectric figure-of-merit ($ZT\sim 2$) in transition metal dichalcogenide monolayers ZrX$_2$ (X= S, Se, Te)

Authors: Ransell D'Souza, Sugata Mukherjee, Sohail Ahmad

Abstract: Two-dimensional group IV transition-metal dichalcogenides have encouraging thermoelectric applications since their electronic and lattice properties can be manipulated with strain. In this paper, we report the thermoelectric parameters such as electrical conductivity, Seebeck coefficients, electrical relaxation times, and the mode dependent contributions to the lattice thermal conductivity of ZrX… ▽ More Two-dimensional group IV transition-metal dichalcogenides have encouraging thermoelectric applications since their electronic and lattice properties can be manipulated with strain. In this paper, we report the thermoelectric parameters such as electrical conductivity, Seebeck coefficients, electrical relaxation times, and the mode dependent contributions to the lattice thermal conductivity of ZrX$_2$ (X = S, Se, Te) from first principles methods. Our calculations indicate that due to tensile strain, the powerfactor increases while simultaneously decreasing the lattice thermal conductivity thus enhancing the thermoelectric figure of merit. Tensile strain widens the bandgap which corresponds to higher powerfactor. The lattice thermal conductivity decreases due to the stiffening of the out-of-plane phonon modes thus reducing the anharmonic scattering lifetimes and increasing the thermoelectric figure-of-merit. △ Less

Submitted 28 January, 2020; originally announced January 2020.

Comments: 10 pages, 8 figures

Journal ref: Journal of Applied Physics 126, 214302 (2019)

Homophily based on few attributes can impede structural balance

Authors: Piotr J. Górski, Klavdiya Bochenina, Janusz A. Hołyst, Raissa M. D'Souza

Abstract: Two complementary mechanisms are thought to shape social groups: homophily between agents and structural balance in connected triads. Here we consider $N$ fully connected agents, where each agent has $G$ underlying attributes, and the similarity between agents in attribute space (i.e., homophily) is used to determine the link weight between them. To incorporate structural balance we use a triad-up… ▽ More Two complementary mechanisms are thought to shape social groups: homophily between agents and structural balance in connected triads. Here we consider $N$ fully connected agents, where each agent has $G$ underlying attributes, and the similarity between agents in attribute space (i.e., homophily) is used to determine the link weight between them. To incorporate structural balance we use a triad-updating rule where only one attribute of one agent is changed intentionally in each update, but this also leads to accidental changes in link weights and even link polarities. The link weight dynamics in the limit of large $G$ is described by a Fokker-Planck equation from which the conditions for a phase transition to a fully balanced state with all links positive can be obtained. This "paradise state" of global cooperation is, however, difficult to achieve requiring $G > O(N^2)$ and $p>0.5$, where the parameter $p$ captures a willingness to consensus. Allowing edge weights to be a consequence of attributes naturally captures homophily and reveals that many real-world social systems would have a subcritical number of attributes necessary to achieve structural balance. △ Less

Submitted 17 January, 2020; originally announced January 2020.

Comments: 6 pages, 5 figures

Journal ref: Phys. Rev. Lett. 125, 078302 (2020)

The Saga of M81: Global View of a Massive Stellar Halo in Formation

Authors: Adam Smercina, Eric F. Bell, Paul A. Price, Colin T. Slater, Richard D'Souza, Jeremy Bailin, Roelof S. de Jong, In Sung Jang, Antonela Monachesi, David Nidever

Abstract: Recent work has shown that Milky Way-mass galaxies display an incredible range of stellar halo properties, yet the origin of this diversity is unclear. The nearby galaxy M81 $-$ currently interacting with M82 and NGC 3077 $-$ sheds unique light on this problem. We present a Subaru Hyper Suprime-Cam survey of the resolved stellar populations around M81, revealing M81's stellar halo in never-before-… ▽ More Recent work has shown that Milky Way-mass galaxies display an incredible range of stellar halo properties, yet the origin of this diversity is unclear. The nearby galaxy M81 $-$ currently interacting with M82 and NGC 3077 $-$ sheds unique light on this problem. We present a Subaru Hyper Suprime-Cam survey of the resolved stellar populations around M81, revealing M81's stellar halo in never-before-seen detail. We resolve the halo to unprecedented $V$-band equivalent surface brightnesses of 33 mag arcsec$^{-2}$, and produce the first-ever global stellar mass density map for a Milky Way-mass stellar halo outside of the Local Group. Using the minor axis, we confirm M81's halo as one of the lowest mass and metal-poorest known ($M_{\star} \simeq 1.16{\times}10^9 M_{\odot}$, [Fe/H] $\simeq {-}1.2$) $-$ indicating a relatively quiet prior accretion history. Yet, our global halo census finds that tidally unbound material from M82 and NGC 3077 provides a substantial infusion of metal-rich material ($M_{\star} \simeq 5.4{\times}10^8$ $M_{\odot}$, [Fe/H] $\simeq {-}$0.9). We further show that, following the accretion of its massive satellite M82 (and the LMC-like NGC 3077), M81 will host one of the most massive and metal-rich stellar halos in the nearby universe. Thus, the saga of M81: following a passive history, M81's merger with M82 will completely transform its halo from a low-mass, anemic halo rivaling the MW, to a metal-rich behemoth rivaled only by systems such as M31. This dramatic transformation indicates that the observed diversity in stellar halo properties is primarily driven by diversity in the largest mergers these galaxies have experienced. △ Less

Submitted 26 October, 2020; v1 submitted 31 October, 2019; originally announced October 2019.

Comments: 25 pages, 15 figures, 3 tables, accepted to ApJ

Explosive Phenomena in Complex Networks

Authors: Raissa M. D'Souza, Jesus Gómez-Gardeñes, Jan Nagler, Alex Arenas

Abstract: The emergence of large-scale connectivity and synchronization are crucial to the structure, function and failure of many complex socio-technical networks. Thus, there is great interest in analyzing phase transitions to large-scale connectivity and to global synchronization, including how to enhance or delay the onset. These phenomena are traditionally studied as second-order phase transitions wher… ▽ More The emergence of large-scale connectivity and synchronization are crucial to the structure, function and failure of many complex socio-technical networks. Thus, there is great interest in analyzing phase transitions to large-scale connectivity and to global synchronization, including how to enhance or delay the onset. These phenomena are traditionally studied as second-order phase transitions where, at the critical threshold, the order parameter increases rapidly but continuously. In 2009, an extremely abrupt transition was found for a network growth process where links compete for addition in attempt to delay percolation. This observation of "explosive percolation" was ultimately revealed to be a continuous transition in the thermodynamic limit, yet with very atypical finite-size scaling, and it started a surge of work on explosive phenomena and their consequences. Many related models are now shown to yield discontinuous percolation transitions and even hybrid transitions. Explosive percolation enables many other features such as multiple giant components, modular structures, discrete scale invariance and non-self-averaging, relating to properties found in many real phenomena such as explosive epidemics, electric breakdowns and the emergence of molecular life. Models of explosive synchronization provide an analytic framework for the dynamics of abrupt transitions and reveal the interplay between the distribution in natural frequencies and the network structure, with applications ranging from epileptic seizures to waking from anesthesia. Here we review the vast literature on explosive phenomena and synthesize the fundamental connections between models and survey the application areas. We attempt to classify explosive phenomena based on underlying mechanisms and to provide a coherent overview and perspective for future research to address the many vital questions that remained unanswered. △ Less

Submitted 23 July, 2019; originally announced July 2019.

Comments: Accepted to Advances in Physics, Review article, 82 pages, 44 figures

Diversity of structural controllability of complex networks with given degree sequence

Authors: Abdorasoul Ghasemi, Márton Pósfai, Raissa M. D'Souza

Abstract: We investigate to what extent the degree sequence of a directed network constrains the number of driver nodes. We develop a pair of algorithms that take a directed degree sequence as input and aim to output a network with the maximum or minimum number of driver nodes. We find an upper bound for the maximum and a lower bound for the minimum, and show that the algorithms achieve these bounds for all… ▽ More We investigate to what extent the degree sequence of a directed network constrains the number of driver nodes. We develop a pair of algorithms that take a directed degree sequence as input and aim to output a network with the maximum or minimum number of driver nodes. We find an upper bound for the maximum and a lower bound for the minimum, and show that the algorithms achieve these bounds for all real and model networks, with few exceptions characterized by tiny system size and heterogeneous degree distributions. Applying these algorithms to a broad range of real networks, we show the gap between the upper and lower bounds can vary dramatically across different degree sequences. Thus, we introduce the notion of \bl{structural control} complexity to capture how much more difficult it is to control a specific network beyond what is required given its degree sequence, suggesting additional structure is present. Using model networks, we numerically and analytically investigate how typical features of the degree distribution affect the range of required driver nodes. We find that the minimum is determined by the number of sources or sinks, while the maximum is strongly affected by the presence of hubs. △ Less

Submitted 23 March, 2020; v1 submitted 16 May, 2019; originally announced May 2019.

Comments: 12 pages, 8 figures

Koopman Operator and its Approximations for Systems with Symmetries

Authors: Anastasiya Salova, Jeffrey Emenheiser, Adam Rupe, James P. Crutchfield, Raissa M. D'Souza

Abstract: Nonlinear dynamical systems with symmetries exhibit a rich variety of behaviors, including complex attractor-basin portraits and enhanced and suppressed bifurcations. Symmetry arguments provide a way to study these collective behaviors and to simplify their analysis. The Koopman operator is an infinite dimensional linear operator that fully captures a system's nonlinear dynamics through the linear… ▽ More Nonlinear dynamical systems with symmetries exhibit a rich variety of behaviors, including complex attractor-basin portraits and enhanced and suppressed bifurcations. Symmetry arguments provide a way to study these collective behaviors and to simplify their analysis. The Koopman operator is an infinite dimensional linear operator that fully captures a system's nonlinear dynamics through the linear evolution of functions of the state space. Importantly, in contrast with local linearization, it preserves a system's global nonlinear features. We demonstrate how the presence of symmetries affects the Koopman operator structure and its spectral properties. In fact, we show that symmetry considerations can also simplify finding the Koopman operator approximations using the extended and kernel dynamic mode decomposition methods (EDMD and kernel DMD). Specifically, representation theory allows us to demonstrate that an isotypic component basis induces block diagonal structure in operator approximations, revealing hidden organization. Practically, if the data is symmetric, the EDMD and kernel DMD methods can be modified to give more efficient computation of the Koopman operator approximation and its eigenvalues, eigenfunctions, and eigenmodes. Rounding out the development, we discuss the effect of measurement noise. △ Less

Submitted 25 April, 2019; originally announced April 2019.

Astro2020 Science White Paper: Construction of an L* Galaxy: the Transformative Power of Wide Fields for Revealing the Past, Present and Future of the Great Andromeda System

Authors: Karoline M. Gilbert, Erik J. Tollerud, Jay Anderson, Rachael L. Beaton, Eric F. Bell, Alyson Brooks, Thomas M. Brown, James Bullock, Jeffrey L. Carlin, Michelle Collins, Andrew Cooper, Denija Crnojevic, Julianne Dalcanton, Andres del Pino, Richard D'Souza, Ivanna Escala, Mark Fardal, Andreea Font, Marla Geha, Puragra Guhathakurta, Evan Kirby, Geraint F. Lewis, Jennifer L. Marshall, Nicolas F. Martin, Kristen McQuinn , et al. (12 additional authors not shown)

Abstract: The Great Andromeda Galaxy (M31) is the nexus of the near-far galaxy evolution connection and a principal data point for near-field cosmology. Due to its proximity (780 kpc), M31 can be resolved into individual stars like the Milky Way (MW). Unlike the MW, we have the advantage of a global view of M31, enabling M31 to be observed with techniques that also apply to more distant galaxies. Moreover,… ▽ More The Great Andromeda Galaxy (M31) is the nexus of the near-far galaxy evolution connection and a principal data point for near-field cosmology. Due to its proximity (780 kpc), M31 can be resolved into individual stars like the Milky Way (MW). Unlike the MW, we have the advantage of a global view of M31, enabling M31 to be observed with techniques that also apply to more distant galaxies. Moreover, recent evidence suggests that M31 may have survived a major merger within the last several Gyr, shaping the morphology of its stellar halo and triggering a starburst, while leaving the stellar disk largely intact. The MW and M31 thus provide complementary opportunities for in-depth studies of the disks, halos, and satellites of L* galaxies. Our understanding of the M31 system will be transformed in the 2020s if they include wide field facilities for both photometry (HST-like sensitivity and resolution) and spectroscopy (10-m class telescope, >1 sq. deg. field, highly multiplexed, R~ 3000 to 6000). We focus here on the power of these facilities to constrain the past, present, and future merger history of M31, via chemo-dynamical analyses and star formation histories of phase-mixed stars accreted at early times, as well as stars in surviving tidal debris features, M31's extended disk, and intact satellite galaxies that will eventually be tidally incorporated into the halo. This will yield an unprecedented view of the hierarchical formation of the M31 system and the subhalos that built it into the L* galaxy we observe today. △ Less

Submitted 1 April, 2019; originally announced April 2019.

Comments: Submitted as a science white paper to the Astro2020 Decadal Survey

Why understanding multiplex social network structuring processes will help us better understand the evolution of human behavior

Authors: Curtis Atkisson, Piotr J. Górski, Matthew O. Jackson, Janusz A. Hołyst, Raissa M. D'Souza

Abstract: Social scientists have long appreciated that relationships between individuals cannot be described from observing a single domain, and that the structure across domains of interaction can have important effects on outcomes of interest (e.g., cooperation).1 One debate explicitly about this surrounds food sharing. Some argue that failing to find reciprocal food sharing means that some process other… ▽ More Social scientists have long appreciated that relationships between individuals cannot be described from observing a single domain, and that the structure across domains of interaction can have important effects on outcomes of interest (e.g., cooperation).1 One debate explicitly about this surrounds food sharing. Some argue that failing to find reciprocal food sharing means that some process other than reciprocity must be occurring, whereas others argue for models that allow reciprocity to span domains in the form of trade.2 Multilayer networks, high-dimensional networks that allow us to consider multiple sets of relationships at the same time, are ubiquitous and have consequences, so processes giving rise to them are important social phenomena. The analysis of multi-dimensional social networks has recently garnered the attention of the network science community.3 Recent models of these processes show how ignoring layer interdependencies can lead one to miss why a layer formed the way it did, and/or draw erroneous conclusions.6 Understanding the structuring processes that underlie multiplex networks will help understand increasingly rich datasets, giving more accurate and complete pictures of social interactions. △ Less

Submitted 27 May, 2020; v1 submitted 26 March, 2019; originally announced March 2019.

Consensus ranking for multi-objective interventions in multiplex networks

Authors: Márton Pósfai, Niklas Braun, Brianne A. Beisner, Brenda McCowan, Raissa M. D'Souza

Abstract: High-centrality nodes have disproportionate influence on the behavior of a network; therefore controlling such nodes can efficiently steer the system to a desired state. Existing multiplex centrality measures typically rank nodes assuming the layers are qualitatively similar. Many real systems, however, are comprised of networks heterogeneous in nature, for example, social networks may have both a… ▽ More High-centrality nodes have disproportionate influence on the behavior of a network; therefore controlling such nodes can efficiently steer the system to a desired state. Existing multiplex centrality measures typically rank nodes assuming the layers are qualitatively similar. Many real systems, however, are comprised of networks heterogeneous in nature, for example, social networks may have both agnostic and affiliative layers. Here, we use rank aggregation methods to identify intervention targets in multiplex networks when the structure, the dynamics, and our intervention goals are qualitatively different for each layer. Our approach is to rank the nodes separately in each layer considering their different function and desired outcome, and then we use Borda count or Kemeny aggregation to identify a consensus ranking - top nodes in the consensus ranking are expected to effectively balance the competing goals simultaneously among all layers. To demonstrate the effectiveness of consensus ranking, we apply our method to a degree-based node removal procedure such that we aim to destroy the largest component in some layers, while maintaining large-scale connectivity in others. For any multi-objective intervention, optimal targets only exist in the Pareto-sense; we, therefore, use a weighted generalization of consensus ranking to investigate the trade-off between the competing objectives. We use a collection of model and real networks to systematically investigate how this trade-off is affected by multiplex network structure. We use the copula representation of the multiplex centrality distributions to generate model multiplex networks with given rank correlations. This allows use to separately manipulate the marginal centrality distribution of each layer and the interdependence between the layers and to independently investigate the role of the two using both analytical and numerical methods. △ Less

Submitted 7 June, 2019; v1 submitted 5 March, 2019; originally announced March 2019.

Journal ref: New J. Phys., 21:055001 (2019)

Competitive percolation strategies for network recovery

Authors: Andrew M. Smith, Márton Pósfai, Martin Rohden, Andrés D. Gonzáles, Leonardo Dueńas-Osorio, Raissa M. D'Souza

Abstract: Restoring operation of critical infrastructure systems after catastrophic events is an important issue, inspiring work in multiple fields, including network science, civil engineering, and operations research. We consider the problem of finding the optimal order of repairing elements in power grids and similar infrastructure. Most existing methods either only consider system network structure, pot… ▽ More Restoring operation of critical infrastructure systems after catastrophic events is an important issue, inspiring work in multiple fields, including network science, civil engineering, and operations research. We consider the problem of finding the optimal order of repairing elements in power grids and similar infrastructure. Most existing methods either only consider system network structure, potentially ignoring important features, or incorporate component level details leading to complex optimization problems with limited scalability. We aim to narrow the gap between the two approaches. Analyzing realistic recovery strategies, we identify over- and undersupply penalties of commodities as primary contributions to reconstruction cost, and we demonstrate traditional network science methods, which maximize the largest connected component, are cost inefficient. We propose a novel competitive percolation recovery model accounting for node demand and supply, and network structure. Our model well approximates realistic recovery strategies, suppressing growth of the largest connected component through a process analogous to explosive percolation. Using synthetic power grids, we investigate the effect of network characteristics on recovery process efficiency. We learn that high structural redundancy enables reduced total cost and faster recovery, however, requires more information at each recovery step. We also confirm that decentralized supply in networks generally benefits recovery efforts. △ Less

Submitted 2 March, 2019; originally announced March 2019.

Comments: 14 pages, 6 figures

Journal ref: Scientific Reports 9, 11843, 2019

Cascading failures in scale-free interdependent networks

Authors: Malgorzata Turalska, Keith Burghardt, Martin Rohden, Ananthram Swami, Raissa M. D'Souza

Abstract: Large cascades are a common occurrence in many natural and engineered complex systems. In this paper we explore the propagation of cascades across networks using realistic network topologies, such as heterogeneous degree distributions, as well as intra- and interlayer degree correlations. We find that three properties, scale-free degree distribution, internal network assortativity, and cross-netwo… ▽ More Large cascades are a common occurrence in many natural and engineered complex systems. In this paper we explore the propagation of cascades across networks using realistic network topologies, such as heterogeneous degree distributions, as well as intra- and interlayer degree correlations. We find that three properties, scale-free degree distribution, internal network assortativity, and cross-network hub-to-hub connections, are all necessary components to significantly reduce the size of large cascades in the Bak-Tang-Wiesenfeld sandpile model. We demonstrate that correlations present in the structure of the multilayer network influence the dynamical cascading process and can prevent failures from spreading across connected layers. These findings highlight the importance of internal and cross-network topology in optimizing stability and robustness of interconnected systems. △ Less

Submitted 19 February, 2019; originally announced February 2019.

Journal ref: Phys. Rev. E 99, 032308 (2019)

An Application of Manifold Learning in Global Shape Descriptors

Authors: Fereshteh S. Bashiri, Reihaneh Rostami, Peggy Peissig, Roshan M. D'Souza, Zeyun Yu

Abstract: With the rapid expansion of applied 3D computational vision, shape descriptors have become increasingly important for a wide variety of applications and objects from molecules to planets. Appropriate shape descriptors are critical for accurate (and efficient) shape retrieval and 3D model classification. Several spectral-based shape descriptors have been introduced by solving various physical equat… ▽ More With the rapid expansion of applied 3D computational vision, shape descriptors have become increasingly important for a wide variety of applications and objects from molecules to planets. Appropriate shape descriptors are critical for accurate (and efficient) shape retrieval and 3D model classification. Several spectral-based shape descriptors have been introduced by solving various physical equations over a 3D surface model. In this paper, for the first time, we incorporate a specific group of techniques in statistics and machine learning, known as manifold learning, to develop a global shape descriptor in the computer graphics domain. The proposed descriptor utilizes the Laplacian Eigenmap technique in which the Laplacian eigenvalue problem is discretized using an exponential weighting scheme. As a result, our descriptor eliminates the limitations tied to the existing spectral descriptors, namely dependency on triangular mesh representation and high intra-class quality of 3D models. We also present a straightforward normalization method to obtain a scale-invariant descriptor. The extensive experiments performed in this study show that the present contribution provides a highly discriminative and robust shape descriptor under the presence of a high level of noise, random scale variations, and low sampling rate, in addition to the known isometric-invariance property of the Laplace-Beltrami operator. The proposed method significantly outperforms state-of-the-art algorithms on several non-rigid shape retrieval benchmarks. △ Less

Submitted 8 January, 2019; originally announced January 2019.

A Dynamic Niche Model for the Emergence and Evolution of Mutualistic Network Structures

Authors: Weiran Cai, Jordan Snyder, Alan Hastings, Raissa M. D'Souza

Abstract: Mutualistic interactions are vital constituents of ecological and socio-economic systems. Empirical studies have found that the patterns of reciprocal relations among the participants often shows the salient features of being simultaneously nested and modular. Whether and how these two structural properties of mutualistic networks can emerge out of a common mechanism however remains unclear. We pr… ▽ More Mutualistic interactions are vital constituents of ecological and socio-economic systems. Empirical studies have found that the patterns of reciprocal relations among the participants often shows the salient features of being simultaneously nested and modular. Whether and how these two structural properties of mutualistic networks can emerge out of a common mechanism however remains unclear. We propose a unified dynamic model based on the adaptation of niche relations that gives rise to both structural features. We apply Hutchinson's concept of niche interaction to networked cooperative species. Their niche relation evolves under the assumption of fitness maximization. Modularity and nestedness emerge concurrently through the accumulated local advantages in the structural and demographic distribution. A rich ensemble of key dynamical behaviors are unveiled in the dynamical framework. We demonstrate that mutualism can exhibit either a stabilizing or destabilizing effect on the evolved network, which undergoes a drastic transition with the overall competition level. Most strikingly, the adaptive network may exhibit a profound nature of history-dependency in response to environmental changes, allowing it to be found in alternative stable structures. The adaptive nature of niche interactions, as captured in our framework, can underlie a broad class of ecological relations and also socio-economic networks that engage in bipartite cooperation. △ Less

Submitted 9 December, 2018; originally announced December 2018.

Comments: 18 pages, 4 figures

Enhancement of thermoelectric figure-of-merit of Graphene upon BN-doping and sample length reduction

Authors: Ransell D'Souza, Sugata Mukherjee

Abstract: Using first-principles density functional perturbation theory based calculations of length-dependent lattice thermal conductivity (\k{appa} L ) and using our previously calculated results (Phys Rev B 95 085435 (2017)) of electrical transport, we report results of thermoelectric figure-of-merit (ZT ) of monolayer and bilayer Graphene. We find nearly ten-fold increase in ZT for the graphene sample d… ▽ More Using first-principles density functional perturbation theory based calculations of length-dependent lattice thermal conductivity (\k{appa} L ) and using our previously calculated results (Phys Rev B 95 085435 (2017)) of electrical transport, we report results of thermoelectric figure-of-merit (ZT ) of monolayer and bilayer Graphene. We find nearly ten-fold increase in ZT for the graphene sample doped with boron nitride and reduced sample length. We also compare \k{appa} L calculated using the iterative real space method with conventional analytical Callaway-Klemens method and obtain the flexural (ZA) phonon modes to be dominant in thermal transport unlike in the latter method. Our calculations are in good agreement with available experimental data. △ Less

Submitted 4 September, 2018; originally announced September 2018.

Comments: 11 pages, 9 figures. Recently accepted in Journal of Applied Physics

Talent and experience shape competitive social hierarchies

Authors: Márton Pósfai, Raissa M. D'Souza

Abstract: Hierarchy of social organization is a ubiquitous property of animal and human groups, linked to resource allocation, collective decisions, individual health, and even to social instability. Experimental evidence shows that both intrinsic abilities of individuals and social reinforcement processes impact hierarchies; existing mathematical models, however, focus on the latter. Here, we develop a rig… ▽ More Hierarchy of social organization is a ubiquitous property of animal and human groups, linked to resource allocation, collective decisions, individual health, and even to social instability. Experimental evidence shows that both intrinsic abilities of individuals and social reinforcement processes impact hierarchies; existing mathematical models, however, focus on the latter. Here, we develop a rigorous model that incorporates both features and explore their synergistic effect on stability and the structure of hierarchy. For pairwise interactions, we show that there is a trade-off between relationship stability and having the most talented individuals in the highest ranks. Extending this to open societies, where individuals enter and leave the population, we show that important societal effects arise from the interaction between talent and social processes: (i) despite positive global correlation between talent and rank, paradoxically, local correlation is negative, and (ii) the removal of an individual can induce a series of rank reversals. We show that the mechanism underlying the latter is the removal of an older individual of limited talent, who nonetheless was able to suppress the rise of younger, more talented individuals. △ Less

Submitted 12 August, 2018; originally announced August 2018.

Comments: Accepted to Physical Review E Rapid Communications

Journal ref: Márton Pósfai and Raissa M. D'Souza. Phys. Rev. E 98, 020302(R) (2018)

The Andromeda galaxy's most important merger about 2 billion years ago as M32's likely progenitor

Authors: Richard D'Souza, Eric F. Bell

Abstract: Although the proximity of the Andromeda galaxy (M31) offers an opportunity to understand how mergers affect galaxies, uncertainty remains about M31's most important mergers. Previous studies focused individually on the giant stellar stream or the impact of M32 on M31's disk, thereby suggesting many substantial satellite interactions. Yet models of M31's disk heating and the similarity between the… ▽ More Although the proximity of the Andromeda galaxy (M31) offers an opportunity to understand how mergers affect galaxies, uncertainty remains about M31's most important mergers. Previous studies focused individually on the giant stellar stream or the impact of M32 on M31's disk, thereby suggesting many substantial satellite interactions. Yet models of M31's disk heating and the similarity between the stellar populations of different tidal substructures in M31's outskirts both suggested a single large merger. M31's stellar halo (its outer low-surface-brightness regions) is built up from the tidal debris of satellites and provides information about its important mergers. Here we use cosmological models of galaxy formation to show that M31's massive and metal-rich stellar halo, containing intermediate-age stars, dramatically narrows the range of allowed interactions, requiring a single dominant merger with a large galaxy (with stellar mass about 2.5 x 10^10 solar masses, the third largest member of the Local Group) about 2 Gyr ago. This single event explains many observations that were previously considered separately: M31's compact and metal-rich satellite M32 is likely to be the stripped core of the disrupted galaxy, its rotating inner stellar halo contains most of the merger debris, and the giant stellar stream is likely to have been thrown out during the merger. This interaction may explain M31's global burst of star formation about 2 Gyr ago in which approximately a fifth of its stars were formed. Moreover, M31's disk and bulge were already in place, suggesting that mergers of this magnitude need not dramatically affect galaxy structure. △ Less

Submitted 23 July, 2018; originally announced July 2018.

Comments: Published in Nature Astronomy. Available here https://www.nature.com/articles/s41550-018-0533-x

A Lonely Giant: The Sparse Satellite Population of M94 Challenges Galaxy Formation

Authors: Adam Smercina, Eric F. Bell, Paul A. Price, Richard D'Souza, Colin T. Slater, Jeremy Ballin, Antonela Monachesi, David Nidever

Abstract: The dwarf satellites of `giant' Milky Way (MW)-mass galaxies are our primary probes of low-mass dark matter halos. The number and velocities of the satellite galaxies of the MW and M31 initially puzzled galaxy formation theorists, but are now reproduced well by many models. Yet, are the MW's and M31's satellites representative? Were galaxy formation models `overfit'? These questions motivate deep… ▽ More The dwarf satellites of `giant' Milky Way (MW)-mass galaxies are our primary probes of low-mass dark matter halos. The number and velocities of the satellite galaxies of the MW and M31 initially puzzled galaxy formation theorists, but are now reproduced well by many models. Yet, are the MW's and M31's satellites representative? Were galaxy formation models `overfit'? These questions motivate deep searches for satellite galaxies outside the Local Group. We present a deep survey of the `classical' satellites ($M_{\star}$$\geqslant$4$\times$10$^5 M_{\odot}$) of the MW-mass galaxy M94 out to 150 kpc projected distance. We find $only\ two$ satellites, each with $M_{\star}{\sim}10^6 M_{\odot}$, compared with 6-12 such satellites in the four other MW-mass systems with comparable data (MW, M31, M81, M101). Using a `standard' prescription for occupying dark matter halos (taken from the fully hydrodynamical EAGLE simulation) with galaxies, we find that such a sparse satellite population occurs in $<$0.2% of MW-mass systems $-$ a $<$1% probability among a sample of five (known systems + M94). In order to produce an M94-like system more frequently we make satellite galaxy formation much more stochastic than is currently predicted by dramatically increasing the slope and scatter of the stellar mass-halo mass (SMHM) relation. Surprisingly, the SMHM relation must be altered even for halos masses up to 10$^{11}M_{\odot} -$ significantly above the mass scales predicted to have increased scatter from current hydrodynamical models. The sparse satellite population of this `lonely giant' thus advocates for an important modification to ideas of how the satellites around MW-mass galaxies form. △ Less

Submitted 20 August, 2018; v1 submitted 10 July, 2018; originally announced July 2018.

Comments: 9 pages, 4 figures, 1 table, published + additional reference

Journal ref: The Astrophysical Journal, 2018, 863, 152

The Prevalence of Type III Disc Breaks in HI-rich and Low-spin Galaxies

Authors: Jing Wang, Zheng Zheng, Richard D'Souza, Houjun Mo, Gyula Józsa, Cheng Li, Peter Kamphuis, Barbara Catinella, Li Shao, Claudia del P. Lagos, Min Du, Zhizheng Pan

Abstract: We investigate the origin of the type III (up-bending) discs based on a sample of 1808 galaxies from SDSS and a sub-sample of 286 galaxies with HI data from ALFALFA. We examine how the type III fraction f3, the fraction of disc galaxies which host up-bending disc breaks, depends on other galactic properties. We confirm that f3 strongly depends on the stellar concentration of galaxies. We find that… ▽ More We investigate the origin of the type III (up-bending) discs based on a sample of 1808 galaxies from SDSS and a sub-sample of 286 galaxies with HI data from ALFALFA. We examine how the type III fraction f3, the fraction of disc galaxies which host up-bending disc breaks, depends on other galactic properties. We confirm that f3 strongly depends on the stellar concentration of galaxies. We find that HI-rich galaxies with low spins tend to have significantly more type III disc breaks than control galaxies, which are matched in concentration and stellar mass. This trend is independent of the existence of strong bars or environment of the galaxies. This result is broadly consistent with predictions from theoretical simulations, and indicates in-situ star formation fueled by gas accretion to be an important mechanism that builds the outer discs of type III galaxies. △ Less

Submitted 25 June, 2018; originally announced June 2018.

Comments: 17 pages, key results in Figure 5, accepted for publication at MNRAS

Length dependent lattice thermal conductivity of single \& multi layered hexagonal boron nitride: A first-principles study using the Callaway-Klemens \& real space super cell methods

Authors: Ransell D'Souza, Sugata Mukherjee

Abstract: The phonon dispersion, density of states, Grüneisen parameters, and the lattice thermal conductivity of single- and multi-layered boron nitride were calculated using first-principles methods. For the bulk {\it h}-BN we also report the two-phonon density of states. We also present simple analytical solutions to the acoustic vibrational mode-dependent lattice thermal conductivity. Moreover, computat… ▽ More The phonon dispersion, density of states, Grüneisen parameters, and the lattice thermal conductivity of single- and multi-layered boron nitride were calculated using first-principles methods. For the bulk {\it h}-BN we also report the two-phonon density of states. We also present simple analytical solutions to the acoustic vibrational mode-dependent lattice thermal conductivity. Moreover, computations based on the elaborate Callaway-Klemens and the real space super cell methods are presented to calculate the sample length and temperature dependent lattice thermal conductivity of single- and multi-layered hexagonal boron nitride which shows good agreement with experimental data. △ Less

Submitted 7 November, 2017; originally announced November 2017.

Comments: 13 pages, 7 figures

Journal ref: Phys. Rev. B 96, 205422 (2017)

Influence of Interface Geometry on Phase Stability and Bandgap Engineering in Boron Nitride substituted Graphene: A Combined First-principles and Monte Carlo Study

Authors: Ransell D'Souza, Sugata Mukherjee, Tanusri Saha-Dasgupta

Abstract: Using combination of Density Functional Theory and Monte Carlo simulation, we study the phase stability and electronic properties of two dimensional hexagonal composites of boron nitride and graphene, with a goal to uncover the role of the interface geometry formed between the two. Our study highlights that preferential creation of extended armchair interfaces may facilitate formation of solid sol… ▽ More Using combination of Density Functional Theory and Monte Carlo simulation, we study the phase stability and electronic properties of two dimensional hexagonal composites of boron nitride and graphene, with a goal to uncover the role of the interface geometry formed between the two. Our study highlights that preferential creation of extended armchair interfaces may facilitate formation of solid solution of boron nitride and graphene within a certain temperature range. We further find that for band-gap engineering, armchair interfaces or patchy interfaces with mixed geometry are most suitable. Extending the study to nanoribbon geometry shows that reduction of dimensionality makes the tendency to phase segregation of the two phases even stronger. Our thorough study should form an useful database in designing boron nitride-graphene composites with desired properties. △ Less

Submitted 10 July, 2017; originally announced July 2017.

Comments: 8 pages, 6 figures