-
One Stream or Two -- Exploring Andromeda's North West Stream
Authors:
Janet Preston,
Denis Erkal,
Michelle L. M. Collins,
R. Michael Rich,
Rodrigo Ibata,
Maxime Delorme
Abstract:
We present results of our dynamical stream modelling for the North West Stream in the outer halo of the Andromeda galaxy (M31). Comprising two main segments, the North West Stream was thought to be a single structured arching around M31. However, recent evidence suggests that it is two separate, unrelated, streams. To test this hypothesis we use observational data from 6 fields associated with the…
▽ More
We present results of our dynamical stream modelling for the North West Stream in the outer halo of the Andromeda galaxy (M31). Comprising two main segments, the North West Stream was thought to be a single structured arching around M31. However, recent evidence suggests that it is two separate, unrelated, streams. To test this hypothesis we use observational data from 6 fields associated with the upper segment of the North West Stream together with 8 fields and 5 globular clusters associated with the lower segment to constrain model orbits. We fit both segments of the stream using a fixed potential model for M31 and an orbit integrator to compare orbits with the observed streams. We measure the central tracks and predict proper motions for for the upper segment (lower segment) finding ${μ^*_α}$ = 0.078$^{+0.015}_{-0.012}$ (0.085$^{+0.001}_{-0.002}$) mas/yr and ${μ_δ}$ = $-$0.05$^{+0.008}_{-0.009}$ ($-$0.095$^{+0.003}_{-0.005}$) mas/yr. Our results support the hypothesis that the dwarf spheroidal galaxy Andromeda XXVII is the progenitor of the upper segment of the North West Stream and that the upper and lower segments do not comprise a single structure. We propose that the upper segment, which appears to be on an infall trajectory with M31, be renamed the "Andromeda XXVII Stream" and the lower segment, also apparently infalling towards M31, retain the name "North West Stream".
△ Less
Submitted 11 September, 2024;
originally announced September 2024.
-
EDGE: The shape of dark matter haloes in the faintest galaxies
Authors:
Matthew D. A. Orkney,
Ethan Taylor,
Justin I. Read,
Martin P. Rey,
Andrew Pontzen,
Oscar Agertz,
Stacy Y. Kim,
Maxime Delorme
Abstract:
Collisionless Dark Matter Only (DMO) structure formation simulations predict that Dark Matter (DM) haloes are prolate in their centres and triaxial towards their outskirts. The addition of gas condensation transforms the central DM shape to be rounder and more oblate. It is not clear, however, whether such shape transformations occur in `ultra-faint' dwarfs, which have extremely low baryon fractio…
▽ More
Collisionless Dark Matter Only (DMO) structure formation simulations predict that Dark Matter (DM) haloes are prolate in their centres and triaxial towards their outskirts. The addition of gas condensation transforms the central DM shape to be rounder and more oblate. It is not clear, however, whether such shape transformations occur in `ultra-faint' dwarfs, which have extremely low baryon fractions. We present the first study of the shape and velocity anisotropy of ultra-faint dwarf galaxies that have gas mass fractions of $f_{\rm gas}(r<R_{\rm half}) < 0.06$. These dwarfs are drawn from the Engineering Dwarfs at Galaxy formation's Edge (EDGE) project, using high resolution simulations that allow us to resolve DM halo shapes within the half light radius ($\sim 100\,$pc). We show that gas-poor ultra-faints ($M_{\rm 200c} \leqslant 1.5\times10^9\,$M$_\odot$; $f_{\rm gas} < 10^{-5}$) retain their pristine prolate DM halo shape even when gas, star formation and feedback are included. This could provide a new and robust test of DM models. By contrast, gas-rich ultra-faints ($M_{\rm 200c} > 3\times10^9\,$M$_\odot$; $f_{\rm gas} > 10^{-4}$) become rounder and more oblate within $\sim 10$ half light radii. Finally, we find that most of our simulated dwarfs have significant radial velocity anisotropy that rises to $\tildeβ > 0.5$ at $R \gtrsim 3 R_{\rm half}$. The one exception is a dwarf that forms a rotating gas/stellar disc because of a planar, major merger. Such strong anisotropy should be taken into account when building mass models of gas-poor ultra-faints.
△ Less
Submitted 5 September, 2023; v1 submitted 24 February, 2023;
originally announced February 2023.
-
EDGE: the puzzling ellipticity of Eridanus II's star cluster and its implications for dark matter at the heart of an ultra-faint dwarf
Authors:
Matthew D. A. Orkney,
Justin I. Read,
Oscar Agertz,
Andrew Pontzen,
Martin P. Rey,
Alex Goater,
Ethan Taylor,
Stacy Y. Kim,
Maxime Delorme
Abstract:
The Eridanus II (EriII) 'ultra-faint' dwarf has a large ($15\,\text{pc}$) and low mass ($4.3\times10^3\,\text{M}_\odot$) star cluster (SC) offset from its centre by $23\pm3\,\text{pc}$ in projection. Its size and offset are naturally explained if EriII has a central dark matter core, but such a core may be challenging to explain in a $Λ$CDM cosmology. In this paper, we revisit the survival and evo…
▽ More
The Eridanus II (EriII) 'ultra-faint' dwarf has a large ($15\,\text{pc}$) and low mass ($4.3\times10^3\,\text{M}_\odot$) star cluster (SC) offset from its centre by $23\pm3\,\text{pc}$ in projection. Its size and offset are naturally explained if EriII has a central dark matter core, but such a core may be challenging to explain in a $Λ$CDM cosmology. In this paper, we revisit the survival and evolution of EriII's SC, focussing for the first time on its puzzlingly large ellipticity ($0.31^{+0.05}_{-0.06}$). We perform a suite of 960 direct $N$-body simulations of SCs, orbiting within a range of spherical background potentials fit to ultra-faint dwarf (UFD) galaxy simulations. We find only two scenarios that come close to explaining EriII's SC. In the first, EriII has a low density dark matter core (of size $\sim70\,\text{pc}$ and density $\lesssim2\times10^8\,\text{M}_{\odot}\,\text{kpc}^{-3}$). In this model, the high ellipticity of EriII's SC is set at birth, with the lack of tidal forces in the core allowing its ellipticity to remain frozen in for long times. In the second, EriII's SC orbits in a partial core, with its high ellipticity owing to its imminent tidal destruction. However, this latter model struggles to reproduce the large size of EriII's SC, and it predicts substantial tidal tails around EriII's SC that should have already been seen in the data. This leads us to favour the cored model. We discuss potential caveats to these findings, and the implications of the cored model for galaxy formation and the nature of dark matter.
△ Less
Submitted 1 August, 2022; v1 submitted 31 January, 2022;
originally announced January 2022.
-
EDGE: Two routes to dark matter core formation in ultra-faint dwarfs
Authors:
Matthew D. A. Orkney,
Justin I. Read,
Martin P. Rey,
Imran Nasim,
Andrew Pontzen,
Oscar Agertz,
Stacy Y. Kim,
Maxime Delorme,
Walter Dehnen
Abstract:
In the standard Lambda cold dark matter paradigm, pure dark matter simulations predict dwarf galaxies should inhabit dark matter haloes with a centrally diverging density `cusp'. This is in conflict with observations that typically favour a constant density `core'. We investigate this `cusp-core problem' in `ultra-faint' dwarf galaxies simulated as part of the `Engineering Dwarfs at Galaxy formati…
▽ More
In the standard Lambda cold dark matter paradigm, pure dark matter simulations predict dwarf galaxies should inhabit dark matter haloes with a centrally diverging density `cusp'. This is in conflict with observations that typically favour a constant density `core'. We investigate this `cusp-core problem' in `ultra-faint' dwarf galaxies simulated as part of the `Engineering Dwarfs at Galaxy formation's Edge' (EDGE) project. We find, similarly to previous work, that gravitational potential fluctuations within the central region of the simulated dwarfs kinematically heat the dark matter particles, lowering the dwarfs' central dark matter density. However, these fluctuations are not exclusively caused by gas inflow/outflow, but also by impulsive heating from minor mergers. We use the genetic modification approach on one of our dwarf's initial conditions to show how a delayed assembly history leads to more late minor mergers and, correspondingly, more dark matter heating. This provides a mechanism by which even ultra-faint dwarfs ($M_* < 10^5\,\text{M}_{\odot}$), in which star formation was fully quenched at high redshift, can have their central dark matter density lowered over time. In contrast, we find that late major mergers can regenerate a central dark matter cusp, if the merging galaxy had sufficiently little star formation. The combination of these effects leads us to predict significant stochasticity in the central dark matter density slopes of the smallest dwarfs, driven by their unique star formation and mass assembly histories.
△ Less
Submitted 7 January, 2021;
originally announced January 2021.
-
Formation of the largest galactic cores through binary scouring and gravitational wave recoil
Authors:
Imran Nasim,
Alessia Gualandris,
Justin I. Read,
Fabio Antonini,
Walter Dehnen,
Maxime Delorme
Abstract:
Massive elliptical galaxies are typically observed to have central cores in their projected radial light profiles. Such cores have long been thought to form through `binary scouring' as supermassive black holes (SMBHs), brought in through mergers, form a hard binary and eject stars from the galactic centre. However, the most massive cores, like the ~3kpc core in A2261-BCG, remain challenging to ex…
▽ More
Massive elliptical galaxies are typically observed to have central cores in their projected radial light profiles. Such cores have long been thought to form through `binary scouring' as supermassive black holes (SMBHs), brought in through mergers, form a hard binary and eject stars from the galactic centre. However, the most massive cores, like the ~3kpc core in A2261-BCG, remain challenging to explain in this way. In this paper, we run a suite of dry galaxy merger simulations to explore three different scenarios for central core formation in massive elliptical galaxies: `binary scouring', `tidal deposition' and `gravitational wave (GW) induced recoil'. Using the Griffin code, we self-consistently model the stars, dark matter and SMBHs in our merging galaxies, following the SMBH dynamics through to the formation of a hard binary. We find that we can only explain the large surface brightness core of A2261-BCG with a combination of a major merger that produces a small ~1kpc core through binary scouring, followed by the subsequent GW recoil of its SMBH that acts to grow the core size. We show that this same model can also explain the bright `knots' observed in the core region of A2261-BCG. Key predictions of this scenario are an offset SMBH surrounded by a compact cluster of bound stars and a non-divergent central density profile. We show that the bright `knots' observed in the core region of A2261-BCG are best explained as stalled perturbers resulting from minor mergers, though the brightest may also represent ejected SMBHs surrounded by a stellar cloak of bound stars.
△ Less
Submitted 11 February, 2021; v1 submitted 9 November, 2020;
originally announced November 2020.
-
Integer Linear Programming for the Tutor Allocation Problem: A Practical Case in a British University
Authors:
Giulia Caselli,
Maxence Delorme,
Manuel Iori
Abstract:
In the Tutor Allocation Problem, the objective is to assign a set of tutors to a set of workshops in order to maximize tutors' preferences. The problem is solved every year by many universities, each having its own specific set of constraints. In this work, we study the tutor allocation in the School of Mathematics at the University of Edinburgh, and solve it with an integer linear programming mod…
▽ More
In the Tutor Allocation Problem, the objective is to assign a set of tutors to a set of workshops in order to maximize tutors' preferences. The problem is solved every year by many universities, each having its own specific set of constraints. In this work, we study the tutor allocation in the School of Mathematics at the University of Edinburgh, and solve it with an integer linear programming model. We tested the model on the 2019/2020 case, obtaining a significant improvement with respect to the manual assignment in use. Further tests on randomly created instances show that the model can be used to address cases of broad interest. We also provide meaningful insights on how input parameters, such as the number of workshop locations and the length of the tutors' preference list, might affect the performance of the model and the average number of preferences satisfied.
△ Less
Submitted 19 May, 2020;
originally announced May 2020.
-
Defeating stochasticity: coalescence timescales of massive black holes in galaxy mergers
Authors:
Imran Nasim,
Alessia Gualandris,
Justin Read,
Walter Dehnen,
Maxime Delorme,
Fabio Antonini
Abstract:
The coalescence of massive black hole binaries (BHBs) in galactic mergers is the primary source of gravitational waves (GWs) at low frequencies. Current estimates of GW detection rates for the Laser Interferometer Space Antenna and the Pulsar Timing Array vary by three orders of magnitude. To understand this variation, we simulate the merger of equal-mass, eccentric, galaxy pairs with central mass…
▽ More
The coalescence of massive black hole binaries (BHBs) in galactic mergers is the primary source of gravitational waves (GWs) at low frequencies. Current estimates of GW detection rates for the Laser Interferometer Space Antenna and the Pulsar Timing Array vary by three orders of magnitude. To understand this variation, we simulate the merger of equal-mass, eccentric, galaxy pairs with central massive black holes and shallow inner density cusps. We model the formation and hardening of a central BHB using the Fast Multiple Method as a force solver, which features a $O(N)$ scaling with the number $N$ of particles and obtains results equivalent to direct-summation simulations. At $N \sim 5\times 10^5$, typical for contemporary studies, the eccentricity of the BHBs can vary significantly for different random realisations of the same initial condition, resulting in a substantial variation of the merger timescale. This scatter owes to the stochasticity of stellar encounters with the BHB and decreases with increasing $N$. We estimate that $N \sim 10^7$ within the stellar half-light radius suffices to reduce the scatter in the merger timescale to $\sim 10$\%. Our results suggest that at least some of the uncertainty in low-frequency GW rates owes to insufficient numerical resolution.
△ Less
Submitted 29 April, 2020;
originally announced April 2020.
-
Uncovering the Orbit of the Hercules Dwarf Galaxy
Authors:
Alexandra L. Gregory,
Michelle L. M. Collins,
Denis Erkal,
Erik Tollerud,
Maxime Delorme,
Lewis Hill,
David J. Sand,
Jay Strader,
Beth Willman
Abstract:
We present new chemo--kinematics of the Hercules dwarf galaxy based on Keck II-- DEIMOS spectroscopy. Our 21 confirmed members have a systemic velocity of $v_{\mathrm{Herc}}=46.4\pm1.3$ kms$^{-1}$ and a velocity dispersion $σ_{v,\mathrm{Herc}}=4.4^{+1.4}_{-1.2}$ kms$^{-1}$. From the strength of the Ca II triplet, we obtain a metallicity of [Fe/H]= $-2.48\pm0.19$ dex and dispersion of…
▽ More
We present new chemo--kinematics of the Hercules dwarf galaxy based on Keck II-- DEIMOS spectroscopy. Our 21 confirmed members have a systemic velocity of $v_{\mathrm{Herc}}=46.4\pm1.3$ kms$^{-1}$ and a velocity dispersion $σ_{v,\mathrm{Herc}}=4.4^{+1.4}_{-1.2}$ kms$^{-1}$. From the strength of the Ca II triplet, we obtain a metallicity of [Fe/H]= $-2.48\pm0.19$ dex and dispersion of $σ_{\rm{[Fe/H]}}= 0.63^{+0.18}_{-0.13}$ dex. This makes Hercules a particularly metal--poor galaxy, placing it slightly below the standard mass--metallicity relation. Previous photometric and spectroscopic evidence suggests that Hercules is tidally disrupting and may be on a highly radial orbit. From our identified members, we measure no significant velocity gradient. By cross--matching with the second \textit{Gaia} data release, we determine an uncertainty--weighted mean proper motion of $μ_α^*=μ_α\cos(δ)=-0.153\pm{0.074}$ mas yr$^{-1}$, $μ_δ=-0.397\pm0.063$ mas yr$^{-1}$. This proper motion is slightly misaligned with the elongation of Hercules, in contrast to models which suggest that any tidal debris should be well aligned with the orbital path. Future observations may resolve this tension.
△ Less
Submitted 30 November, 2019;
originally announced December 2019.
-
A Dwarf Disrupting -- Andromeda XXVII and the North West Stream
Authors:
Janet Preston,
Michelle L. M. Collins,
Rodrigo A. Ibata,
Erik J. Tollerud,
R. Michael Rich,
Ana Bonaca,
Alan W. McConnachie,
Dougal Mackey,
Geraint F. Lewis,
Nicolas F. Martin,
Jorge Peñarrubia,
Scott C. Chapman,
Maxime Delorme
Abstract:
We present a kinematic and spectroscopic analysis of 38 red giant branch stars, in 7 fields, spanning the dwarf spheroidal galaxy Andromeda XXVII and the upper segment of the North West Stream. Both features are located in the outer halo of the Andromeda galaxy at a projected radius of 50-80 kpc, with the stream extending for $\sim$3$^{\circ}$ on the sky. Our data is obtained as part of the PAndAS…
▽ More
We present a kinematic and spectroscopic analysis of 38 red giant branch stars, in 7 fields, spanning the dwarf spheroidal galaxy Andromeda XXVII and the upper segment of the North West Stream. Both features are located in the outer halo of the Andromeda galaxy at a projected radius of 50-80 kpc, with the stream extending for $\sim$3$^{\circ}$ on the sky. Our data is obtained as part of the PAndAS survey and enables us to confirm that Andromeda XXVII's heliocentric distance is 827 $\pm$ 47 kpc and spectroscopic metallicity is -2.1$^{+0.4}_{-0.5}$. We also re-derive Andromeda XXVII's kinematic properties, measuring a systemic velocity = -526.1$^{+10.0}_{-11.0}${\kms} and a velocity dispersion that we find to be non-Gaussian but for which we derive a formal value of 27.0$^{+2.2}_{-3.9}${\kms}. In the upper segment of the North West Stream we measure mean values for the metallicity = -1.8$\pm$0.4, systemic velocity = -519.4 $\pm$4.0{\kms} and velocity dispersion = 10.0$\pm$4.0{\kms}. We also detect a velocity gradient of 1.7$\pm$0.3 {\kms} kpc$^{-1}$ on an infall trajectory towards M31. With a similar gradient, acting in the same direction, in the lower segment we suggest that the North West Stream is not a single structure. As the properties of the upper segment of the North West Stream and Andromeda XXVII are consistent within ~90\% confidence limits, it is likely that the two are related and plausible that Andromeda XXVII is the progenitor of this stream.
△ Less
Submitted 20 September, 2019;
originally announced September 2019.
-
Mathematical models for stable matching problems with ties and incomplete lists
Authors:
Maxence Delorme,
Sergio García,
Jacek Gondzio,
Joerg Kalcsics,
David Manlove,
William Pettersson
Abstract:
We present new integer linear programming (ILP) models for NP-hard optimisation problems in instances of the Stable Marriage problem with Ties and Incomplete lists (SMTI) and its many-to-one generalisation, the Hospitals / Residents problem with Ties (HRT). These models can be used to efficiently solve these optimisation problems when applied to (i) instances derived from real-world applications,…
▽ More
We present new integer linear programming (ILP) models for NP-hard optimisation problems in instances of the Stable Marriage problem with Ties and Incomplete lists (SMTI) and its many-to-one generalisation, the Hospitals / Residents problem with Ties (HRT). These models can be used to efficiently solve these optimisation problems when applied to (i) instances derived from real-world applications, and (ii) larger instances that are randomly-generated. In the case of SMTI, we consider instances arising from the pairing of children with adoptive families, where preferences are obtained from a quality measure of each possible pairing of child to family. In this case we seek a maximum weight stable matching. We present new algorithms for preprocessing instances of SMTI with ties on both sides, as well as new ILP models. Algorithms based on existing state-of-the-art models only solve 6 of our 22 real-world instances within an hour per instance, and our new models solve all 22 instances within a mean runtime of 60 seconds. For HRT, we consider instances derived from the problem of assigning junior doctors to foundation posts in Scottish hospitals. Here we seek a maximum size stable matching. We show how to extend our models for SMTI to the HRT case. For the real instances, we reduce the mean runtime from an average of 144 seconds when using state-of-the-art methods, to 3 seconds when using our new ILP-based algorithms. We also show that our models outperform considerably state-of-the-art models on larger randomly-generated instances of SMTI and HRT.
△ Less
Submitted 17 March, 2019; v1 submitted 5 October, 2018;
originally announced October 2018.
-
Generalized arcsine laws for fractional Brownian motion
Authors:
Tridib Sadhu,
Mathieu Delorme,
Kay Jörg Wiese
Abstract:
The three arcsine laws for Brownian motion are a cornerstone of extreme-value statistics. For a Brownian $B_t$ starting from the origin, and evolving during time $T$, one considers the following three observables: (i) the duration $t_+$ the process is positive, (ii) the time $t_{\rm last}$ the process last visits the origin, and (iii) the time $t_{\rm max}$ when it achieves its maximum (or minimum…
▽ More
The three arcsine laws for Brownian motion are a cornerstone of extreme-value statistics. For a Brownian $B_t$ starting from the origin, and evolving during time $T$, one considers the following three observables: (i) the duration $t_+$ the process is positive, (ii) the time $t_{\rm last}$ the process last visits the origin, and (iii) the time $t_{\rm max}$ when it achieves its maximum (or minimum). All three observables have the same cumulative probability distribution expressed as an arcsine function, thus the name of arcsine laws. We show how these laws change for fractional Brownian motion $X_t$, a non-Markovian Gaussian process indexed by the Hurst exponent $H$. It generalizes standard Brownian motion (i.e. $H=\tfrac{1}{2}$). We obtain the three probabilities using a perturbative expansion in $ε= H-\tfrac{1}{2}$. While all three probabilities are different, this distinction can only be made at second order in $ε$. Our results are confirmed to high precision by extensive numerical simulations.
△ Less
Submitted 6 June, 2017;
originally announced June 2017.
-
Probing dark matter with star clusters: a dark matter core in the ultra-faint dwarf Eridanus II
Authors:
Filippo Contenta,
Eduardo Balbinot,
James A. Petts,
Justin I. Read,
Mark Gieles,
Michelle L. M. Collins,
Jorge Peñarrubia,
Maxime Delorme,
Alessia Gualandris
Abstract:
We present a new technique to probe the central dark matter (DM) density profile of galaxies that harnesses both the survival and observed properties of star clusters. As a first application, we apply our method to the `ultra-faint' dwarf Eridanus II (Eri II) that has a lone star cluster ~45 pc from its centre. Using a grid of collisional $N$-body simulations, incorporating the effects of stellar…
▽ More
We present a new technique to probe the central dark matter (DM) density profile of galaxies that harnesses both the survival and observed properties of star clusters. As a first application, we apply our method to the `ultra-faint' dwarf Eridanus II (Eri II) that has a lone star cluster ~45 pc from its centre. Using a grid of collisional $N$-body simulations, incorporating the effects of stellar evolution, external tides and dynamical friction, we show that a DM core for Eri II naturally reproduces the size and the projected position of its star cluster. By contrast, a dense cusped galaxy requires the cluster to lie implausibly far from the centre of Eri II (>1 kpc), with a high inclination orbit that must be observed at a particular orbital phase. Our results, therefore, favour a dark matter core. This implies that either a cold DM cusp was `heated up' at the centre of Eri II by bursty star formation, or we are seeing an evidence for physics beyond cold DM.
△ Less
Submitted 4 April, 2018; v1 submitted 4 May, 2017;
originally announced May 2017.
-
Pickands' constant at first order in an expansion around Brownian motion
Authors:
Mathieu Delorme,
Alberto Rosso,
Kay Jörg Wiese
Abstract:
In the theory of extreme values of Gaussian processes, many results are expressed in terms of the Pickands constant $\mathcal{H}_α$. This constant depends on the local self-similarity exponent $α$ of the process, i.e. locally it is a fractional Brownian motion (fBm) of Hurst index $H=α/2$. Despite its importance, only two values of the Pickands constant are known: ${\cal H}_1 =1$ and…
▽ More
In the theory of extreme values of Gaussian processes, many results are expressed in terms of the Pickands constant $\mathcal{H}_α$. This constant depends on the local self-similarity exponent $α$ of the process, i.e. locally it is a fractional Brownian motion (fBm) of Hurst index $H=α/2$. Despite its importance, only two values of the Pickands constant are known: ${\cal H}_1 =1$ and ${\cal H}_2=1/\sqrtπ$. Here, we extend the recent perturbative approach to fBm to include drift terms. This allows us to investigate the Pickands constant $\mathcal{H}_α$ around standard Brownian motion ($α=1$) and to derive the new exact result $\mathcal{H}_α=1 - (α-1) γ_{\rm E} + \mathcal{O}\!\left( α-1\right)^{2}$.
△ Less
Submitted 26 September, 2016;
originally announced September 2016.
-
Extreme-Value Statistics of Fractional Brownian Motion Bridges
Authors:
Mathieu Delorme,
Kay Jörg Wiese
Abstract:
Fractional Brownian motion is a self-affine, non-Markovian and translationally invariant generalization of Brownian motion, depending on the Hurst exponent $H$. Here we investigate fractional Brownian motion where both the starting and the end point are zero, commonly referred to as bridge processes. Observables are the time $t_+$ the process is positive, the maximum $m$ it achieves, and the time…
▽ More
Fractional Brownian motion is a self-affine, non-Markovian and translationally invariant generalization of Brownian motion, depending on the Hurst exponent $H$. Here we investigate fractional Brownian motion where both the starting and the end point are zero, commonly referred to as bridge processes. Observables are the time $t_+$ the process is positive, the maximum $m$ it achieves, and the time $t_{\rm max}$ when this maximum is taken. Using a perturbative expansion around Brownian motion ($H=\frac12$), we give the first-order result for the probability distribution of these three variables, and the joint distribution of $m$ and $t_{\rm max}$. Our analytical results are tested, and found in excellent agreement, with extensive numerical simulations, both for $H>\frac12$ and $H<\frac12$. This precision is achieved by sampling processes with a free endpoint, and then converting each realization to a bridge process, in generalization to what is usually done for Brownian motion.
△ Less
Submitted 13 May, 2016;
originally announced May 2016.
-
An efficient positive potential-density pair expansion for modelling galaxies
Authors:
Armando Rojas-Niño,
Justin I. Read,
Luis Aguilar,
Maxime Delorme
Abstract:
We present a novel positive potential-density pair expansion for modelling galaxies, based on the Miyamoto-Nagai (MN) disc. By using three sets of such discs, each one of them aligned along each symmetry axis, we are able to reconstruct a broad range of potentials that correspond to density profiles from exponential discs to 3D power law models with varying triaxiality (henceforth simply "twisted"…
▽ More
We present a novel positive potential-density pair expansion for modelling galaxies, based on the Miyamoto-Nagai (MN) disc. By using three sets of such discs, each one of them aligned along each symmetry axis, we are able to reconstruct a broad range of potentials that correspond to density profiles from exponential discs to 3D power law models with varying triaxiality (henceforth simply "twisted" models). We increase the efficiency of our expansion by allowing the scale length parameter of each disc to be negative. We show that, for suitable priors on the scale length and height parameters, these "MNn discs" have just one negative density minimum. This allows us to ensure global positivity by demanding that the total density at the global minimum is positive. We find that at better than 10\% accuracy in our density reconstruction, we can represent a radial and vertical exponential disc over $0.1-10$ scale lengths/heights with 4 MNn discs, an NFW profile over $0.1-10$ scale lengths with 4 MNn discs, and a twisted triaxial NFW profile with 3 MNn discs per symmetry axis. Our expansion is efficient, fully analytic, and well-suited to reproducing the density distribution and gravitational potential of galaxies from discs to ellipsoids.
△ Less
Submitted 13 April, 2016;
originally announced April 2016.
-
Perturbative Expansion for the Maximum of Fractional Brownian Motion
Authors:
Mathieu Delorme,
Kay Jörg Wiese
Abstract:
Brownian motion is the only random process which is Gaussian, stationary and Markovian. Dropping the Markovian property, i.e. allowing for memory, one obtains a class of processes called fractional Brownian motion, indexed by the Hurst exponent $H$. For $H=1/2$, Brownian motion is recovered. We develop a perturbative approach to treat the non-locality in time in an expansion in…
▽ More
Brownian motion is the only random process which is Gaussian, stationary and Markovian. Dropping the Markovian property, i.e. allowing for memory, one obtains a class of processes called fractional Brownian motion, indexed by the Hurst exponent $H$. For $H=1/2$, Brownian motion is recovered. We develop a perturbative approach to treat the non-locality in time in an expansion in $\varepsilon = H-1/2$. This allows us to derive analytic results beyond scaling exponents for various observables related to extreme value statistics: The maximum $m$ of the process and the time $t_{\text{max}}$ at which this maximum is reached, as well as their joint distribution. We test our analytical predictions with extensive numerical simulations for different values of $H$. They show excellent agreement, even for $H$ far from $1/2$.
△ Less
Submitted 2 March, 2016;
originally announced March 2016.
-
Distribution of joint local and total size and of extension for avalanches in the Brownian force model
Authors:
Mathieu Delorme,
Pierre Le Doussal,
Kay Jörg Wiese
Abstract:
The Brownian force model (BFM) is a mean-field model for the local velocities during avalanches in elastic interfaces of internal space dimension $d$, driven in a random medium. It is exactly solvable via a non-linear differential equation. We study avalanches following a kick, i.e. a step in the driving force. We first recall the calculation of the distributions of the global size (total swept ar…
▽ More
The Brownian force model (BFM) is a mean-field model for the local velocities during avalanches in elastic interfaces of internal space dimension $d$, driven in a random medium. It is exactly solvable via a non-linear differential equation. We study avalanches following a kick, i.e. a step in the driving force. We first recall the calculation of the distributions of the global size (total swept area) and of the local jump size for an arbitrary kick amplitude. We extend this calculation to the joint density of local and global sizes within a single avalanche, in the limit of an infinitesimal kick. When the interface is driven by a single point we find new exponents $τ_0=5/3$ and $τ=7/4$, depending on whether the force or the displacement is imposed. We show that the extension of a single avalanche along one internal direction (i.e. the total length in $d=1$) is finite and we calculate its distribution, following either a local or a global kick. In all cases it exhibits a divergence $P(\ell) \sim \ell^{-3}$ at small $\ell$. Most of our results are tested in a numerical simulation in dimension $d=1$.
△ Less
Submitted 19 January, 2016;
originally announced January 2016.
-
The Maximum of a Fractional Brownian Motion: Analytic Results from Perturbation Theory
Authors:
Mathieu Delorme,
Kay Joerg Wiese
Abstract:
Fractional Brownian motion is a non-Markovian Gaussian process $X_t$, indexed by the Hurst exponent $H$. It generalises standard Brownian motion (corresponding to $H=1/2$). We study the probability distribution of the maximum $m$ of the process and the time $t_{\rm max}$ at which the maximum is reached. They are encoded in a path integral, which we evaluate perturbatively around a Brownian, settin…
▽ More
Fractional Brownian motion is a non-Markovian Gaussian process $X_t$, indexed by the Hurst exponent $H$. It generalises standard Brownian motion (corresponding to $H=1/2$). We study the probability distribution of the maximum $m$ of the process and the time $t_{\rm max}$ at which the maximum is reached. They are encoded in a path integral, which we evaluate perturbatively around a Brownian, setting $H=1/2 + \varepsilon$. This allows us to derive analytic results beyond the scaling exponents. Extensive numerical simulations for different values of $H$ test these analytical predictions and show excellent agreement, even for large $\varepsilon$.
△ Less
Submitted 22 July, 2015;
originally announced July 2015.
-
Bulking II: Classifications of Cellular Automata
Authors:
Marianne Delorme,
Jacques Mazoyer,
Nicolas Ollinger,
Guillaume Theyssier
Abstract:
This paper is the second part of a series of two papers dealing with bulking: a way to define quasi-order on cellular automata by comparing space-time diagrams up to rescaling. In the present paper, we introduce three notions of simulation between cellular automata and study the quasi-order structures induced by these simulation relations on the whole set of cellular automata. Various aspects of t…
▽ More
This paper is the second part of a series of two papers dealing with bulking: a way to define quasi-order on cellular automata by comparing space-time diagrams up to rescaling. In the present paper, we introduce three notions of simulation between cellular automata and study the quasi-order structures induced by these simulation relations on the whole set of cellular automata. Various aspects of these quasi-orders are considered (induced equivalence relations, maximum elements, induced orders, etc) providing several formal tools allowing to classify cellular automata.
△ Less
Submitted 24 September, 2010; v1 submitted 29 January, 2010;
originally announced January 2010.