Showing 1–45 of 45 results for author: Frankel, N

The Age-Dependent Vertical Actions of Young Stars in the Galaxy

Authors: D. N. Garzon, Neige Frankel, Eleonora Zari, Maosheng Xiang, Hans-Walter Rix

Abstract: Stars in the Galactic disk are born on cold, nearly circular orbits with small vertical excursions. After their birth, their orbits evolve, driven by small- or large-scale perturbations in the Galactic disk's gravitational potential. Here, we study the vertical motions of young stars over their first few orbital periods, using a sample of OBA stars from \textit{Gaia} E/DR3, which includes radial v… ▽ More Stars in the Galactic disk are born on cold, nearly circular orbits with small vertical excursions. After their birth, their orbits evolve, driven by small- or large-scale perturbations in the Galactic disk's gravitational potential. Here, we study the vertical motions of young stars over their first few orbital periods, using a sample of OBA stars from \textit{Gaia} E/DR3, which includes radial velocities and ages $τ$ from LAMOST. We constructed a parametric model for the time evolution of the stellar orbits' mean vertical actions $J_z$ as a function of Galactocentric radius, $R_{\mathrm{GC}}$. Accounting for data uncertainties, we use Markov Chain Monte Carlo (MCMC) analysis in annuli of Galactocentric radius to constrain the model parameters. Our best-fit model shows a remarkably linear increase of vertical actions with age across all Galactocentric radii examined. Orbital \textit{heating} by random scattering could offer a straightforward interpretation for this trend. However, various other dynamical aspects of the Galactic disk, such as stars being born in a warped disk, might offer alternative explanations that could be tested in the future. △ Less

Submitted 9 July, 2024; originally announced July 2024.

Comments: 10 pages, 5 figures. Draft version July 9, 2024. Typeset using LATEX default style in AASTeX631

Iron Snails: non-equilibrium dynamics and spiral abundance patterns

Authors: Neige Frankel, David W. Hogg, Scott Tremaine, Adrian Price-Whelan, Jeff Shen

Abstract: Galaxies are not in a dynamical steady state. They continually undergo perturbations, e.g., from infalling dwarf galaxies and dark-matter substructure. After a dynamical perturbation, stars phase mix towards a new steady state; in so doing they generally form spiral structures, such as spiral density waves in galaxy disks and the Gaia Snail observed in the vertical phase-space density in the solar… ▽ More Galaxies are not in a dynamical steady state. They continually undergo perturbations, e.g., from infalling dwarf galaxies and dark-matter substructure. After a dynamical perturbation, stars phase mix towards a new steady state; in so doing they generally form spiral structures, such as spiral density waves in galaxy disks and the Gaia Snail observed in the vertical phase-space density in the solar neighborhood. Structures in phase-space density can be hard to measure accurately, because spatially varying selection effects imprint their own patterns on the density. However, stellar labels such as metallicity, or other element abundances, or stellar masses and ages, can be measured even in the face of complex or unknown spatial selection functions. We show that if the equilibrium galaxy has phase-space gradients in these labels, any perturbation that could raise a spiral wave in the phase-space density will raise a spiral wave in the distribution of labels as well. We work out the relationship between the spiral patterns in the density and in the labels. As an example, we analyze the Gaia Snail and show that its amplitude and dynamical age as derived from elemental abundances (mainly [Mg/Fe]) follow similar patterns to those derived from the phase-space density. Our best model dates the Snail's perturbation to about 400 Myr ago although we find significant variations with angular momentum in the best-fit age. Conceptually, the ideas presented here are related to Orbital Torus Imaging, chemical tagging, and other methods that use stellar labels to trace dynamics. △ Less

Submitted 9 July, 2024; originally announced July 2024.

Comments: Submitted to ApJ, comments are welcome

The Extremely Metal Rich Knot of Stars at the Heart of the Galaxy

Authors: Hans-Walter Rix, Vedant Chandra, Gail Zasowski, Annalisa Pillepich, Sergey Khoperskov, Sofia Feltzing, Rosemary F. Wyse, Neige Frankel, Danny Horta, Juna Kollmeier, Keivan G. Stassun, Melissa Ness, Jonathan C. Bird, David L. Nidever, Jose G. Fernandez, João A. Amarante, Chervin F. Laporte, Jianhui Lian

Abstract: We show with Gaia XP spectroscopy that extremely metal-rich stars in the Milky Way (EMR; $[M/H]_{XP} > 0.5$) - but only those - are largely confined to a tight "knot" at the center of the Galaxy. This EMR knot is round in projection, has a fairly abrupt edge near $\sim 1.5$kpc, and is a dynamically hot system. This central knot also contains very metal-rich (VMR; $+0.2\le [M/H]_{XP} \le +0.4$) sta… ▽ More We show with Gaia XP spectroscopy that extremely metal-rich stars in the Milky Way (EMR; $[M/H]_{XP} > 0.5$) - but only those - are largely confined to a tight "knot" at the center of the Galaxy. This EMR knot is round in projection, has a fairly abrupt edge near $\sim 1.5$kpc, and is a dynamically hot system. This central knot also contains very metal-rich (VMR; $+0.2\le [M/H]_{XP} \le +0.4$) stars. However, in contrast to EMR stars, the bulk of VMR stars form an extended, highly flattened distribution in the inner Galaxy ($R_{\mathrm{GC}}\lesssim 5$ kpc). We draw on TNG50 simulations of Milky Way analogs for context and find that compact, metal-rich knots confined to $<1.5$kpc are a universal feature. In typical simulated analogs, the top 5-10% most metal-rich stars are confined to a central knot; however, in our Milky Way data this fraction is only 0.1%. Dust-penetrating wide-area near-infrared spectroscopy, such as SDSS-V, will be needed for a rigorous estimate of the fraction of stars in the Galactic EMR knot. Why in our Milky Way only EMR giants are confined to such a central knot remains to be explained. Remarkably, the central few kiloparsecs of the Milky Way harbor both the highest concentration of metal-poor stars (the `poor old heart') and almost all EMR stars. This highlights the stellar population diversity at the bottom of galactic potential wells. △ Less

Submitted 3 June, 2024; originally announced June 2024.

Comments: 11 pages, 7 figures, submitted to ApJ

A path towards constraining the evolution of the interstellar medium and outflows in the Milky Way using APOGEE

Authors: Piyush Sharda, Yuan-Sen Ting, Neige Frankel

Abstract: In recent years, the study of the Milky Way has significantly advanced due to extensive spectroscopic surveys of its stars, complemented by astroseismic and astrometric data. However, it remains disjoint from recent advancements in understanding the physics of the Galactic interstellar medium (ISM). This paper introduces a new model for the chemical evolution of the Milky Way that can be constrain… ▽ More In recent years, the study of the Milky Way has significantly advanced due to extensive spectroscopic surveys of its stars, complemented by astroseismic and astrometric data. However, it remains disjoint from recent advancements in understanding the physics of the Galactic interstellar medium (ISM). This paper introduces a new model for the chemical evolution of the Milky Way that can be constrained on stellar data, because it combines a state-of-the-art ISM model with a Milky Way stellar disc model. Utilizing a dataset of red clump stars from APOGEE, known for their precise ages and metallicities, we concentrate on the last 6 billion years -- a period marked by Milky Way's secular evolution. We examine the oxygen abundance in the low-$α$ disc stars relative to their ages and birth radii, validating or constraining critical ISM parameters that remain largely unexplored in extragalactic observations. The models that successfully reproduce the radius -- metallicity distribution and the age -- metallicity distribution of stars without violating existing ISM observations indicate a need for modest differential oxygen enrichment in Galactic outflows, meaning that the oxygen abundance of outflows is higher than the local ISM abundance, irrespective of outflow mass loading. The models also suggest somewhat elevated ISM gas velocity dispersion levels over the past 6 billion years compared to galaxies of similar mass. The extra turbulence necessary could result from energy from gas accretion onto the Galaxy, supernovae clustering in the ISM, or increased star formation efficiency per freefall time. This work provides a novel approach to constraining the Galactic ISM and outflows, leveraging the detailed insights available from contemporary Milky Way surveys. △ Less

Submitted 29 May, 2024; v1 submitted 28 May, 2024; originally announced May 2024.

Comments: 20 pages, 9 figures. Accepted by MNRAS

What Does the Large Magellanic Cloud Look Like? It Depends on [M/H] and Age

Authors: Neige Frankel, Rene Andrae, Hans-Walter Rix, Joshua Povick, Vedant Chandra

Abstract: We offer a new way to look at the Large Magellanic Cloud through stellar mono-abundance and mono-age-mono-abundance maps. These maps are based on $\gtrsim 500\,000$ member stars with photo-spectroscopic [M/H] and age estimates from Gaia DR3 data, and they are the first area-complete, metallicity- and age-differentiated stellar maps of any disk galaxy. Azimuthally averaged, these maps reveal a surp… ▽ More We offer a new way to look at the Large Magellanic Cloud through stellar mono-abundance and mono-age-mono-abundance maps. These maps are based on $\gtrsim 500\,000$ member stars with photo-spectroscopic [M/H] and age estimates from Gaia DR3 data, and they are the first area-complete, metallicity- and age-differentiated stellar maps of any disk galaxy. Azimuthally averaged, these maps reveal a surprisingly simple picture of the Milky Way's largest satellite galaxy. For any [M/H] below -0.1 dex, the LMC's radial profile is well described by a simple exponential, but with a scale length that steadily shrinks towards higher metallicities, from nearly 2.3~kpc at [M/H]$=-1.8$ to only 0.75~kpc at [M/H]$=-0.25$. The prominence of the bar decreases dramatically with [M/H], making it barely discernible at [M/H]$\lesssim -1.5$. Yet, even for metal-rich populations, the bar has little impact on the azimuthally averaged profile of the mono-abundance components. Including ages, we find that the scale length is a greater function of age than of metallicity, with younger populations far more centrally concentrated. At old ages, the scale length decreases with increasing metallicity; at young ages, the scale-length is independent of metallicity. These findings provide quantitative support for a scenario where the LMC built its stellar structure effectively outside in. △ Less

Submitted 13 March, 2024; originally announced March 2024.

Comments: Submitted to ApJ; constructive comments would be appreciated

Coherent evolution of superexchange interaction in seconds long optical clock spectroscopy

Authors: William R. Milner, Stefan Lannig, Mikhail Mamaev, Lingfeng Yan, Anjun Chu, Ben Lewis, Max N. Frankel, Ross B. Hutson, Ana Maria Rey, Jun Ye

Abstract: Measurement science now connects strongly with engineering of quantum coherence, many-body states, and entanglement. To scale up the performance of an atomic clock using a degenerate Fermi gas loaded in a three-dimensional optical lattice, we must understand complex many-body Hamiltonians to ensure meaningful gains for metrological applications. In this work, we use a near unity filled Sr 3D latti… ▽ More Measurement science now connects strongly with engineering of quantum coherence, many-body states, and entanglement. To scale up the performance of an atomic clock using a degenerate Fermi gas loaded in a three-dimensional optical lattice, we must understand complex many-body Hamiltonians to ensure meaningful gains for metrological applications. In this work, we use a near unity filled Sr 3D lattice to study the effect of a tunable Fermi-Hubbard Hamiltonian. The clock laser introduces a spin-orbit coupling spiral phase and breaks the isotropy of superexchange interactions, changing the Heisenberg spin model into one exhibiting XXZ-type spin anisotropy. By tuning the lattice confinement and applying imaging spectroscopy we map out favorable atomic coherence regimes. With weak transverse confinement, both s- and p-wave interactions contribute to decoherence and atom loss, and their contributions can be balanced. At deep transverse confinement, we directly observe coherent superexchange interactions, tunable via on-site interaction and site-to-site energy shift, on the clock Ramsey fringe contrast over timescales of multiple seconds. This study provides a groundwork for using a 3D optical lattice clock to probe quantum magnetism and spin entanglement △ Less

Submitted 20 February, 2024; originally announced February 2024.

A comparison of the Milky Way's recent star formation revealed by dust thermal emission and high-mass stars

Authors: J. D. Soler, E. Zari, D. Elia, S. Molinari, C. Mininni, E. Schisano, A. Traficante, R. S. Klessen, S. C. O. Glover, P. Hennebelle, T. Colman, N. Frankel, T. Wenger

Abstract: We present a comparison of the Milky Way's star formation rate (SFR) surface density ($Σ_{\rm SFR}$) obtained with two independent state-of-the-art observational methods. The first method infers $Σ_{\rm SFR}$ from observations of the dust thermal emission from interstellar dust grains in far-infrared wavelengths registered in the Herschel infrared Galactic Plane Survey (Hi-GAL). The second method… ▽ More We present a comparison of the Milky Way's star formation rate (SFR) surface density ($Σ_{\rm SFR}$) obtained with two independent state-of-the-art observational methods. The first method infers $Σ_{\rm SFR}$ from observations of the dust thermal emission from interstellar dust grains in far-infrared wavelengths registered in the Herschel infrared Galactic Plane Survey (Hi-GAL). The second method determines $Σ_{\rm SFR}$ by modeling the current population of O-, B-, and A-type stars in a 6 kpc $\times$ 6 kpc area around the Sun. We find an agreement between the two methods within a factor of two for the mean SFRs and the SFR surface density profiles. Given the broad differences between the observational techniques and the independent assumptions in the methods for computing the SFRs, this agreement constitutes a significant advance in our understanding of the star formation of our Galaxy and implies that the local SFR has been roughly constant over the past 10\,Myr. △ Less

Submitted 18 September, 2023; v1 submitted 2 August, 2023; originally announced August 2023.

Comments: Accepted for publication at Astronomy & Astrophysics (11SEP2023). 8 pages, 7 figures

Journal ref: A&A 678, A95 (2023)

Decoding the age-chemical structure of the Milky Way disk: An application of Copulas and Elicitable Maps

Authors: Aarya A. Patil, Jo Bovy, Sebastian Jaimungal, Neige Frankel, Henry W. Leung

Abstract: In the Milky Way, the distribution of stars in the $[α/\mathrm{Fe}]$ vs. $[\mathrm{Fe/H}]$ and $[\mathrm{Fe/H}]$ vs. age planes holds essential information about the history of star formation, accretion, and dynamical evolution of the Galactic disk. We investigate these planes by applying novel statistical methods called copulas and elicitable maps to the ages and abundances of red giants in the A… ▽ More In the Milky Way, the distribution of stars in the $[α/\mathrm{Fe}]$ vs. $[\mathrm{Fe/H}]$ and $[\mathrm{Fe/H}]$ vs. age planes holds essential information about the history of star formation, accretion, and dynamical evolution of the Galactic disk. We investigate these planes by applying novel statistical methods called copulas and elicitable maps to the ages and abundances of red giants in the APOGEE survey. We find that the low- and high-$α$ disk stars have a clean separation in copula space and use this to provide an automated separation of the $α$ sequences using a purely statistical approach. This separation reveals that the high-$α$ disk ends at the same [$α$/Fe] and age at high $[\mathrm{Fe/H}]$ as the low-$[\mathrm{Fe/H}]$ start of the low-$α$ disk, thus supporting a sequential formation scenario for the high- and low-$α$ disks. We then combine copulas with elicitable maps to precisely obtain the correlation between stellar age $τ$ and metallicity $[\mathrm{Fe/H}]$ conditional on Galactocentric radius $R$ and height $z$ in the range $0 < R < 20$ kpc and $|z| < 2$ kpc. The resulting trends in the age-metallicity correlation with radius, height, and [$α$/Fe] demonstrate a $\approx 0$ correlation wherever kinematically-cold orbits dominate, while the naively-expected negative correlation is present where kinematically-hot orbits dominate. This is consistent with the effects of spiral-driven radial migration, which must be strong enough to completely flatten the age-metallicity structure of the low-$α$ disk. △ Less

Submitted 15 June, 2023; originally announced June 2023.

Comments: 18 pages, 16 figures, Submitted to MNRAS

Disk flaring with TNG50: diversity across Milky Way and M31 analogs

Authors: Diego Sotillo-Ramos, Martina Donnari, Annalisa Pillepich, Neige Frankel, Dylan Nelson, Volker Springel, Lars Hernquist

Abstract: We use the sample of 198 Milky Way (MW) and Andromeda (M31) analogs from TNG50 to quantify the level of disk flaring predicted by a modern, high-resolution cosmological hydrodynamical simulation. Disk flaring refers to the increase of vertical stellar disk height with galactocentric distance. The TNG50 galaxies are selected to have stellar disky morphology, a stellar mass in the range of… ▽ More We use the sample of 198 Milky Way (MW) and Andromeda (M31) analogs from TNG50 to quantify the level of disk flaring predicted by a modern, high-resolution cosmological hydrodynamical simulation. Disk flaring refers to the increase of vertical stellar disk height with galactocentric distance. The TNG50 galaxies are selected to have stellar disky morphology, a stellar mass in the range of $M_* = 10^{10.5 - 11.2}~\rm{M_{\odot}}$, and a MW-like Mpc-scale environment at $z=0$. The stellar disks of such TNG50 MW/M31 analogs exhibit a wide diversity of structural properties, including a number of galaxies with disk scalelength and thin and thick disk scaleheights that are comparable to those measured or inferred for the Galaxy and Andromeda. With one set of physical ingredients, TNG50 returns a large variety of flaring flavours and amounts, also for mono-age stellar populations. With this paper, we hence propose a non-parametric characterization of flaring. The typical MW/M31 analogs exhibit disk scaleheights that are $1.5-2$ times larger in the outer than in the inner regions of the disk for both old and young stellar populations, but with a large galaxy-to-galaxy variation. Which stellar population flares more, and by how much, also varies from galaxy to galaxy. TNG50 de facto brackets existing observational constraints for the Galaxy and all previous numerical findings. A link between the amount of flaring and the $z=0$ global galaxy structural properties or merger history is complex. However, a connection between the scaleheights and the local stellar vertical kinematics and gravitational potential is clearly in place. △ Less

Submitted 28 March, 2023; originally announced March 2023.

Comments: Submitted to MNRAS. Main figures: 5, 13. See presentation and data release of TNG50 MW/M31 analogs by Pillepich et al. and see also Ramesh et al. on astro-ph today

The Origin of Stars in the Inner 500 Parsecs in TNG50 Galaxies

Authors: Alina Boecker, Nadine Neumayer, Annalisa Pillepich, Neige Frankel, Rahul Ramesh, Ryan Leaman, Lars Hernquist

Abstract: We investigate the origin of stars in the innermost $500\,\mathrm{pc}$ of galaxies spanning stellar masses of $5\times10^{8-12}\,\mathrm{M}_{\odot}$ at $\mathrm{z=0}$ using the cosmological magnetohydrodynamical TNG50 simulation. Three different origins of stars comprise galactic centers: 1) in-situ (born in the center), 2) migrated (born elsewhere in the galaxy and ultimately moved to the center)… ▽ More We investigate the origin of stars in the innermost $500\,\mathrm{pc}$ of galaxies spanning stellar masses of $5\times10^{8-12}\,\mathrm{M}_{\odot}$ at $\mathrm{z=0}$ using the cosmological magnetohydrodynamical TNG50 simulation. Three different origins of stars comprise galactic centers: 1) in-situ (born in the center), 2) migrated (born elsewhere in the galaxy and ultimately moved to the center), 3) ex-situ (accreted from other galaxies). In-situ and migrated stars dominate the central stellar mass budget on average with 73% and 23% respectively. The ex-situ fraction rises above 1% for galaxies $\gtrsim10^{11}\,\mathrm{M}_{\odot}$. Yet, only 9% of all galaxies exhibit no ex-situ stars in their centers and the scatter of ex-situ mass is significant ($4-6\,\mathrm{dex}$). Migrated stars predominantly originate closely from the center ($1-2\,\mathrm{kpc}$), but if they travelled together in clumps distances reach $\sim10\,\mathrm{kpc}$. Central and satellite galaxies possess similar amounts and origins of central stars. Star forming galaxies ($\gtrsim10^{10}\,\mathrm{M}_{\odot}$) have on average more ex-situ mass in their centers than quenched ones. We predict readily observable stellar population and dynamical properties: 1) migrated stars are distinctly young ($\sim2\,\mathrm{Gyr}$) and rotationally supported, especially for Milky Way mass galaxies, 2) in-situ stars are most metal-rich and older than migrated stars, 3) ex-situ stars are on random motion dominated orbits and typically the oldest, most metal-poor and $α$-enhanced population. We demonstrate that the interaction history with other galaxies leads to diverse pathways of building up galaxy centers in a $Λ$CDM universe. Our work highlights the necessity for cosmological context in formation scenarios of central galactic components and the potential to use galaxy centers as tracers of overall galaxy assembly. △ Less

Submitted 27 January, 2023; originally announced January 2023.

Comments: 24 pages, 13 Figures, published in MNRAS

Vertical motion in the Galactic disc: unwinding the Snail

Authors: Neige Frankel, Jo Bovy, Scott Tremaine, David W. Hogg

Abstract: The distribution of stars in the Milky Way disc shows a spiral structure--the Snail--in the space of velocity and position normal to the Galactic mid-plane. The Snail appears as straight lines in the vertical frequency--vertical phase plane when effects from sample selection are removed. Their slope has the dimension of inverse time, with the simplest interpretation being the inverse age of the Sn… ▽ More The distribution of stars in the Milky Way disc shows a spiral structure--the Snail--in the space of velocity and position normal to the Galactic mid-plane. The Snail appears as straight lines in the vertical frequency--vertical phase plane when effects from sample selection are removed. Their slope has the dimension of inverse time, with the simplest interpretation being the inverse age of the Snail. Here, we devise and fit a simple model in which the spiral starts as a lopsided perturbation from steady state, that winds up into the present-day morphology. The winding occurs because the vertical frequency decreases with vertical action. We use data from stars in Gaia EDR3 that have measured radial velocities, pruned by simple distance and photometric selection functions. We divide the data into boxels of dynamical invariants (radial action, angular momentum); our model fits the data well in many of the boxels. The model parameters have physical interpretations: one, $A$, is a perturbation amplitude, and one, $t$, is interpretable in the simplest models as the time since the event that caused the Snail. We find trends relating the strength and age to angular momentum: (i) the amplitude $A$ is small at low angular momentum ($<1\,600\mathrm{\,kpc\,km\,s}^{-1}$ or guiding-centre radius $< 7.3\,$kpc), and over a factor of three larger, with strong variations, in the outer disc; (ii) there is no single well-defined perturbation time, with $t$ varying between 0.2 and 0.6 Gyr. Residuals between the data and the model display systematic trends, implying that the data call for more complex models. △ Less

Submitted 22 December, 2022; originally announced December 2022.

Comments: 9 pages, 8 figures. Submitted to MNRAS (constructive comments are welcome)

The origin and fate of the Gaia phase-space snail

Authors: Scott Tremaine, Neige Frankel, Jo Bovy

Abstract: The Gaia snail is a spiral feature in the distribution of solar-neighbourhood stars in position and velocity normal to the Galactic midplane. The snail probably arises from phase mixing of gravitational disturbances that perturbed the disc in the distant past. The most common hypothesis is that the strongest disturbance resulted from a passage of the Sagittarius dwarf galaxy close to the solar nei… ▽ More The Gaia snail is a spiral feature in the distribution of solar-neighbourhood stars in position and velocity normal to the Galactic midplane. The snail probably arises from phase mixing of gravitational disturbances that perturbed the disc in the distant past. The most common hypothesis is that the strongest disturbance resulted from a passage of the Sagittarius dwarf galaxy close to the solar neighbourhood. In this paper we investigate the alternative hypothesis that the snail is created by many small disturbances rather than one large one, that is, by Gaussian noise in the gravitational potential. Probably most of this noise is due to substructures in the dark-matter halo. We show that this hypothesis naturally reproduces most of the properties of the snail. In particular it predicts correctly, with no free parameters, that the apparent age of the snail will be $ \sim 0.5$ Gyr. An important ingredient of this model is that any snail-like feature in the solar neighbourhood, whatever its cause, is erased by scattering from giant molecular clouds or other small-scale structure on a time-scale $\lesssim 1$ Gyr. △ Less

Submitted 22 December, 2022; originally announced December 2022.

Comments: 9 pages, 6 figures. Submitted to MNRAS

Did the Milky Way just light up? The recent star formation history of the Galactic disc

Authors: E. Zari, N. Frankel, H. -W. Rix

Abstract: [Abridged] We map the stellar age distribution ($\lesssim 1$ Gyr) across a 6kpc$\,\times\,$6kpc area of the Galactic disc to constrain our Galaxy's recent star-formation history. Our modelling draws on the sample of Zari et al. (2021) that encompasses all presumed disc OBA stars ($\sim 500,000$ sources) with $G<16$. To be less sensitive to reddening, we do not forward model the detailed CMD distri… ▽ More [Abridged] We map the stellar age distribution ($\lesssim 1$ Gyr) across a 6kpc$\,\times\,$6kpc area of the Galactic disc to constrain our Galaxy's recent star-formation history. Our modelling draws on the sample of Zari et al. (2021) that encompasses all presumed disc OBA stars ($\sim 500,000$ sources) with $G<16$. To be less sensitive to reddening, we do not forward model the detailed CMD distribution of these stars, but instead the K-band absolute magnitude distribution, $n(M_K)$, among stars with $M_K<0$ and $T_{\mathrm{eff}} > 7000$K at a certain positions $\vec{x}$ in the disc as a step function with five age bins, $b(τ~|~\vec{x}, \vecα)$, logarithmically-spaced in age from $τ= 5$ Myr to $τ= 1$ Gyr. Given a set of isochrones and a Kroupa (2001) initial mass function, we sample $b(τ\,|\,\vec{x}, \vecα)$ to maximize the likelihood of the data, accounting for the selection function. This results in a set of mono-age stellar density maps across a sizeable portion of the Galactic disc. These maps show that some, but not all, spiral arms are reflected in overdensities of stars younger than 500 Myr. The maps of the youngest stars ($<10$ Myr) trace major star forming regions. The maps at all ages exhibit an outward density gradient and distinct spiral-like spatial structure, which is qualitatively similar on large scales among all age intervals. When summing over the maps' area and extrapolating to the whole disc, we find an effective star-formation rate over the last 10 Myr of $\approx 3.3 \mathrm{M_{\odot}/yr}$, higher than previously published estimates. Remarkably, our stellar age distribution implies that the star-formation rate has been three times lower throughout most of the last Gyr, having risen distinctly only in the very recent past. △ Less

Submitted 6 June, 2022; originally announced June 2022.

Comments: 11 pages, 9 figures. Submitted to A&A: comments welcome

Journal ref: A&A 669, A10 (2023)

High-contrast, high-angular resolution view of the GJ 367 exoplanet system

Authors: Wolfgang Brandner, Per Calissendorff, Neige Frankel, Faustine Cantalloube

Abstract: We search for additional companions in the GJ 367 exoplanet system, and aim at better constraining its age and evolutionary status. We analyse high contrast direct imaging observations obtained with HST/NICMOS, VLT/NACO, and VLT/SPHERE. We investigate and critically discuss conflicting age indicators based on theoretical isochrones and models for Galactic dynamics. A comparison of GAIA EDR3 parall… ▽ More We search for additional companions in the GJ 367 exoplanet system, and aim at better constraining its age and evolutionary status. We analyse high contrast direct imaging observations obtained with HST/NICMOS, VLT/NACO, and VLT/SPHERE. We investigate and critically discuss conflicting age indicators based on theoretical isochrones and models for Galactic dynamics. A comparison of GAIA EDR3 parallax and photometric measurements with theoretical isochrones suggest a young age $\le$60 Myr for GJ 367. The star's Galactic kinematics exclude membership to any nearby young moving group or stellar stream. Its highly eccentric Galactic orbit, however, is atypical for a young star. Age estimates considering Galactic dynamical evolution are most consistent with an age of 1 to 8 Gyr. We find no evidence for a significant mid-infrared excess in the WISE bands, suggesting the absence of warm dust in the GJ 367 system. The direct imaging data provide significantly improved detection limits compared to previous studies. At 530 mas (5 au) separation, the SPHERE data achieve a 5 sigma contrast of $2.6 \times 10^{-6}$. The data exclude the presence of a stellar companion at projected separations $\ge$0.4 au. At projected separations $\ge$5 au we can exclude substellar companions with a mass $\ge$ 1.5 M$_{\rm Jup}$ for an age of 50 Myr, and $\ge$ 20 M$_{\rm Jup}$ for an age of 5 Gyr. By applying the stellar parameters corresponding to the 50 Myr isochrone, we derive a bulk density of $ρ_{\rm planet} = 6.2$ g/cm$^3$ for GJ 367b, which is 25% smaller than a previous estimate. △ Less

Submitted 6 April, 2022; originally announced April 2022.

Comments: 9 pages, 8 figures, accepted by MNRAS

Journal ref: 2022, MNRAS 513, 661

Simulated Bars May Be Shorter But Are Not Slower Than Observed: TNG50 vs. MaNGA

Authors: Neige Frankel, Annalisa Pillepich, Hans-Walter Rix, Vicente Rodriguez-Gomez, Jason Sanders, Jo Bovy, Juna Kollmeier, Norm Murray, Ted Mackereth

Abstract: Galactic bars are prominent dynamical structures within disk galaxies whose size, formation time, strength, and pattern speed influence the dynamical evolution of their hosts galaxies. Yet, their formation and evolution in a cosmological context is not well understood, as cosmological simulation studies have been limited by the classic trade off between simulation volume and resolution. Here we an… ▽ More Galactic bars are prominent dynamical structures within disk galaxies whose size, formation time, strength, and pattern speed influence the dynamical evolution of their hosts galaxies. Yet, their formation and evolution in a cosmological context is not well understood, as cosmological simulation studies have been limited by the classic trade off between simulation volume and resolution. Here we analyze barred disk galaxies in the cosmological magneto-hydrodynamical simulation TNG50 and quantitatively compare the distributions of bar size and pattern speed to those from MaNGA observations at $z=0$. TNG50 galaxies are selected to match the stellar mass and size distributions of observed galaxies, to account for observational selection effects. We find that the high-resolution of TNG50 yields bars with a wide range of pattern speeds (including those with $\geq 40~\mathrm{km}\,\mathrm{s}^{-1}$\,$\mathrm{kpc}^{-1}$) and a mean value of $\sim36~\mathrm{km}\,\mathrm{s}^{-1}\,\mathrm{kpc}$ consistent with observations within $6\,\mathrm{km}\,\mathrm{s}^{-1}$\,$\mathrm{kpc}^{-1}$, in contrast with previous lower-resolution cosmological simulations that produced bars that were too slow. We find, however, that bars in TNG50 are on average $\sim 35\%$ shorter than observed, although this discrepancy may partly reflect remaining inconsistencies in the simulation-data comparison. This leads to higher values of $\mathcal{R} = R_\mathrm{corot}/R_\mathrm{bar}$ in TNG50, but points to simulated bars being `too short' rather than `too slow'. After repeating the analysis on the lower-resolution run of the same simulation (with the same physical model), we qualitatively reproduce the results obtained in previous studies: this implies that, along with physical model variations, numerical resolution effects may explain the previously found `slowness' of simulated bars. △ Less

Submitted 20 January, 2022; originally announced January 2022.

Comments: 14 pages, key figures are Fig. 5 and 6. Submitted to ApJ, comments are welcome

Gravitational wave sources in our Galactic backyard: Predictions for BHBH, BHNS and NSNS binaries detectable with LISA

Authors: Tom Wagg, Floor S. Broekgaarden, Selma E. de Mink, Lieke A. C. van Son, Neige Frankel, Stephen Justham

Abstract: Future searches for gravitational waves from space will be sensitive to double compact objects (DCOs) in our Milky Way. We present new simulations of the populations of double black holes (BHBHs), black hole neutron stars (BHNSs) and double neutron stars (NSNSs) that will be detectable by the planned space-based gravitational wave detector LISA. For our estimates, we use an empirically-informed mo… ▽ More Future searches for gravitational waves from space will be sensitive to double compact objects (DCOs) in our Milky Way. We present new simulations of the populations of double black holes (BHBHs), black hole neutron stars (BHNSs) and double neutron stars (NSNSs) that will be detectable by the planned space-based gravitational wave detector LISA. For our estimates, we use an empirically-informed model of the metallicity dependent star formation history of the Milky Way. We populate it using an extensive suite of binary population-synthesis predictions for varying assumptions relating to mass transfer, common-envelope, supernova kicks, remnant masses and wind mass loss physics. For a 4(10)-year LISA mission, we predict between 30-370(50-550) detections over these variations, out of which 6-154(9-238) are BHBHs, 2-198(3-289) are BHNSs and 3-35(4-57) are NSNSs. We discuss how the variations in the physics assumptions alter the distribution of properties of the detectable systems, even when the detection rates are unchanged. In particular we discuss the observable characteristics such as the chirp mass, eccentricity and sky localisation and how the BHBH, BHNS and NSNS populations can be distinguished, both from each other and from the more numerous double white dwarf population. We further discuss the possibility of multi-messenger observations of pulsar populations with the Square Kilometre Array (SKA) and assess the benefits of extending the LISA mission. △ Less

Submitted 26 November, 2021; originally announced November 2021.

Comments: Submitted. All code for reproducing results and figures is available https://github.com/TomWagg/detecting-DCOs-in-LISA. Data is available https://doi.org/10.5281/zenodo.4699712

Mapping Luminous Hot Stars in the Galaxy

Authors: E. Zari, H. -W. Rix, N. Frankel, M. Xiang, E. Poggio, R. Drimmel, A. Tkachenko

Abstract: [Abridged] Luminous hot stars dominate the stellar energy input to the interstellar medium throughout cosmological time, they are laboratories to test theories of stellar evolution and multiplicity, and they serve as luminous tracers of star formation in the Milky Way and other galaxies. Massive stars occupy well-defined loci in colour-colour and colour-magnitude spaces, enabling selection based o… ▽ More [Abridged] Luminous hot stars dominate the stellar energy input to the interstellar medium throughout cosmological time, they are laboratories to test theories of stellar evolution and multiplicity, and they serve as luminous tracers of star formation in the Milky Way and other galaxies. Massive stars occupy well-defined loci in colour-colour and colour-magnitude spaces, enabling selection based on the combination of Gaia EDR3 astrometry and photometry and 2MASS photometry, even in the presence of substantive dust extinction. In this paper we devise an all-sky sample of such luminous OBA-type stars, designed to be quite complete rather than very pure, to serve as targets for spectroscopic follow-up with the SDSS-V survey. We estimate "astro-kinematic" distances by combining parallaxes and proper motions with a model for the expected velocity and density distribution of young stars; we show that this adds useful constraints on the stars' distances, and hence luminosities. With these distances we map the spatial distribution of a more stringently selected sub-sample across the Galactic disc, and find it to be highly structured, with distinct over- and under-densities. The most evident over-densities can be associated with the presumed spiral arms of the Milky Way, in particular the Sagittarius-Carina and Scutum-Centaurus arms. Yet, the spatial picture of the Milky Way's young disc structure emerging in this study is complex, and suggests that most young stars in our Galaxy ($t_{age}<t_{dyn}$) are not neatly organised into distinct spiral arms. The combination of the comprehensive spectroscopy to come from SDSS-V (yielding velocities, ages, etc..) with future Gaia data releases will be crucial to reveal the dynamical nature of the spiral arm themselves. △ Less

Submitted 17 February, 2021; originally announced February 2021.

Comments: 19 pages, 21 figures; accepted for publication in A&A

Journal ref: A&A 650, A112 (2021)

The GALAH survey: tracing the Galactic disk with Open Clusters

Authors: Lorenzo Spina, Yuan-Sen Ting, Gayandhi M. De Silva, Neige Frankel, Sanjib Sharma, Tristan Cantat-Gaudin, Meridith Joyce, Dennis Stello, Amanda I. Karakas, Martin B. Asplund, Thomas Nordlander, Luca Casagrande, Valentina D'Orazi, Andrew R. Casey, Peter Cottrell, Thor Tepper-García, Martina Baratella, Janez Kos, Klemen Čotar, Joss Bland-Hawthorn, Sven Buder, Ken C. Freeman, Michael R. Hayden, Geraint F. Lewis, Jane Lin , et al. (6 additional authors not shown)

Abstract: Open clusters are unique tracers of the history of our own Galaxy's disk. According to our membership analysis based on \textit{Gaia} astrometry, out of the 226 potential clusters falling in the footprint of GALAH or APOGEE, we find that 205 have secure members that were observed by at least one of the survey. Furthermore, members of 134 clusters have high-quality spectroscopic data that we use to… ▽ More Open clusters are unique tracers of the history of our own Galaxy's disk. According to our membership analysis based on \textit{Gaia} astrometry, out of the 226 potential clusters falling in the footprint of GALAH or APOGEE, we find that 205 have secure members that were observed by at least one of the survey. Furthermore, members of 134 clusters have high-quality spectroscopic data that we use to determine their chemical composition. We leverage this information to study the chemical distribution throughout the Galactic disk of 21 elements, from C to Eu. The radial metallicity gradient obtained from our analysis is $-$0.076$\pm$0.009 dex kpc$^{-1}$, which is in agreement with previous works based on smaller samples. Furthermore, the gradient in the [Fe/H] - guiding radius (r$_{\rm guid}$) plane is $-$0.073$\pm$0.008 dex kpc$^{-1}$. We show consistently that open clusters trace the distribution of chemical elements throughout the Galactic disk differently than field stars. In particular, at given radius, open clusters show an age-metallicity relation that has less scatter than field stars. As such scatter is often interpreted as an effect of radial migration, we suggest that these differences are due to the physical selection effect imposed by our Galaxy: clusters that would have migrated significantly also had higher chances to get destroyed. Finally, our results reveal trends in the [X/Fe]$-$r$_{\rm guid}$$-$age space, which are important to understand production rates of different elements as a function of space and time. △ Less

Submitted 16 February, 2021; v1 submitted 4 November, 2020; originally announced November 2020.

Comments: 18 pages, 10 figures, accepted for publication MNRAS

MSC Class: 85-02

Orbital and spectral characterization of the benchmark T-type brown dwarf HD 19467B

Authors: A. -L. Maire, K. Molaverdikhani, S. Desidera, T. Trifonov, P. Mollière, V. D'Orazi, N. Frankel, J. -L. Baudino, S. Messina, A. Müller, B. Charnay, A. Cheetham, P. Delorme, R. Ligi, M. Bonnefoy, W. Brandner, D. Mesa, F. Cantalloube, R. Galicher, T. Henning, B. A. Biller, J. Hagelberg, A. -M. Lagrange, B. Lavie, E. Rickman , et al. (20 additional authors not shown)

Abstract: Context. Detecting and characterizing substellar companions for which the luminosity, mass, and age can be determined independently is of utter importance to test and calibrate the evolutionary models due to uncertainties in their formation mechanisms. HD 19467 is a bright and nearby star hosting a cool brown dwarf companion detected with RV and imaging, making it a valuable object for such studie… ▽ More Context. Detecting and characterizing substellar companions for which the luminosity, mass, and age can be determined independently is of utter importance to test and calibrate the evolutionary models due to uncertainties in their formation mechanisms. HD 19467 is a bright and nearby star hosting a cool brown dwarf companion detected with RV and imaging, making it a valuable object for such studies. Aims. We aim to further characterize the orbital, spectral, and physical properties of the HD 19467 system. Methods. We present new high-contrast imaging data with the SPHERE and NaCo instruments. We also analyze archival data from HARPS, NaCo, HIRES, UVES, and ASAS. We also use proper motion data of the star from Hipparcos and Gaia. Results. We refine the properties of the host star and derive an age of 8.0$^{+2.0}_{-1.0}$ Gyr based on isochrones, gyrochronology, and chemical and kinematic arguments. This estimate is slightly younger than previous estimates of ~9-11 Gyr. No orbital curvature is seen in the current imaging, RV, and astrometric data. From a joint fit of the data, we refine the orbital parameters for HD 19467B: period 398$^{+95}_{-93}$ yr, inclination 129.8$^{+8.1}_{-5.1}$ deg, eccentricity 0.56$\pm$0.09, longitude of the ascending node 134.8$\pm$4.5 deg, and argument of the periastron 64.2$^{+5.5}_{-6.3}$ deg. We assess a dynamical mass of 74$^{+12}_{-9}$ MJ. The fit with atmospheric models of the spectrophotometric data of HD 19467B indicates an atmosphere without clouds or with very thin clouds, an effective temperature of 1042$^{+77}_{-71}$ K, and a large surface gravity of 5.34$^{+0.08}_{-0.09}$ dex. The comparison to model predictions of the bolometric luminosity and dynamical mass of HD 19467B, assuming our system age estimate, indicates a better agreement with the Burrows et al. models; whereas the other evolutionary models used tend to underestimate its cooling rate. △ Less

Submitted 4 June, 2020; v1 submitted 20 May, 2020; originally announced May 2020.

Comments: Accepted for publication in A&A, 25 pages, 23 figures, 9 tables. Abstract slightly abridged to match arXiv requirements. Updated after language editing

Journal ref: A&A 639, A47 (2020)

The Strength of the Dynamical Spiral Perturbation in the Galactic Disk

Authors: Anna-Christina Eilers, David W. Hogg, Hans-Walter Rix, Neige Frankel, Jason A. S. Hunt, Jean-Baptiste Fouvry, Tobias Buck

Abstract: The mean Galactocentric radial velocities $\langle v_{R}\rangle(R,\varphi)$ of luminous red giant stars within the mid-plane of the Milky Way reveal a spiral signature, which could plausibly reflect the response to a non-axisymmetric perturbation of the gravitational potential in the Galactic disk. We apply a simple steady-state toy model of a logarithmic spiral to interpret these observations, an… ▽ More The mean Galactocentric radial velocities $\langle v_{R}\rangle(R,\varphi)$ of luminous red giant stars within the mid-plane of the Milky Way reveal a spiral signature, which could plausibly reflect the response to a non-axisymmetric perturbation of the gravitational potential in the Galactic disk. We apply a simple steady-state toy model of a logarithmic spiral to interpret these observations, and find a good qualitative and quantitative match. Presuming that the amplitude of the gravitational potential perturbation is proportionate to that in the disk's surface mass density, we estimate the surface mass density amplitude to be $Σ_{\rm max} (R_{\odot})\approx 5.5\,\rm M_{\odot}\,pc^{-2}$ at the solar radius when choosing a fixed pattern speed of $Ω_{\mathrm p}=12\,\rm km\,s^{-1}\,kpc^{-1}$. Combined with the local disk density, this implies a surface mass density contrast between the arm and inter-arm regions of approximately $\pm 10\%$ at the solar radius, with an increases towards larger radii. Our model constrains the pitch angle of the dynamical spiral arms to be approximately $12^{\circ}$. △ Less

Submitted 3 August, 2020; v1 submitted 2 March, 2020; originally announced March 2020.

Comments: Accepted for publication in ApJ

Keeping it Cool: Much Orbit Migration, yet Little Heating, in the Galactic Disk

Authors: Neige Frankel, Jason Sanders, Yuan-Sen Ting, Hans-Walter Rix

Abstract: A star in the Milky Way's disk can now be at a Galactocentric radius quite distant from its birth radius for two reasons: either its orbit has become eccentric through radial heating, which increases its radial action $J_R$ (`blurring'); or merely its angular momentum $L_z$ has changed and thereby its guiding radius (`churning'). We know that radial orbit migration is strong in the Galactic low-… ▽ More A star in the Milky Way's disk can now be at a Galactocentric radius quite distant from its birth radius for two reasons: either its orbit has become eccentric through radial heating, which increases its radial action $J_R$ (`blurring'); or merely its angular momentum $L_z$ has changed and thereby its guiding radius (`churning'). We know that radial orbit migration is strong in the Galactic low-$α$ disk and set out to quantify the relative importance of these two effects, by devising and applying a parameterized model for the distribution $p(L_z, J_R, τ, \mathrm[Fe/H])$ in the stellar disk. This model describes the orbit evolution for stars of age $τ$ and metallicity [Fe/H], presuming coeval stars were initially born on (near-)circular orbits, and with a unique [Fe/H] at a given birth angular momentum and age. We fit this model to APOGEE red clump stars, accounting for the complex selection function of the survey. The best fit model implies changes of angular momentum of $\sqrt{\langle ΔL_z \rangle^2} \approx 619\, \mathrm{kpc~km/s~}(τ/\mathrm{6~Gyr})^{0.5}$, and changes of radial action as $\sqrt{\langle ΔJ_R \rangle^2} \approx 63\, \mathrm{kpc~km/s~} (τ/\mathrm{6~Gyr})^{0.6}$ at 8 kpc. This suggests that the secular orbit evolution of the disk is dominated by diffusion in angular momentum, with radial heating being an order of magnitude lower. △ Less

Submitted 11 February, 2020; originally announced February 2020.

Comments: Submitted to ApJ, comments are welcome

The Inside-out Growth of the Galactic Disk

Authors: Neige Frankel, Jason Sanders, Hans-Walter Rix, Yuan-Sen Ting, Melissa Ness

Abstract: We quantify the inside-out growth of the Milky Way's low-alpha stellar disk, modelling the ages, metallicities and Galactocentric radii of APOGEE red clump stars with 6 < R < 13 kpc. The current stellar distribution differs significantly from that expected from the star formation history due to the redistribution of stars through radial orbit mixing. We propose and fit a global model for the Milky… ▽ More We quantify the inside-out growth of the Milky Way's low-alpha stellar disk, modelling the ages, metallicities and Galactocentric radii of APOGEE red clump stars with 6 < R < 13 kpc. The current stellar distribution differs significantly from that expected from the star formation history due to the redistribution of stars through radial orbit mixing. We propose and fit a global model for the Milky Way disk, specified by an inside-out star formation history, radial orbit mixing, and an empirical, parametric model for its chemical evolution. We account for the spatially complex survey selection function, and find that the model fits all data well. We find distinct inside-out growth of the Milky Way disk; the best fit model implies that the half-mass radius of the Milky Way disk has grown by 43\% over the last 7 Gyr. Yet, such inside-out growth still results in present-day age gradient weaker than 0.1 Gyr/kpc. Our model predicts the half-mass and half-light sizes of the Galactic disk at earlier epochs, which can be compared to the observed redshift -size relations of disk galaxies. We show that radial orbit migration can reconcile the distinct disk-size evolution with redshift, also expected from cosmological simulations, with the modest present-day age gradients seen in the Milky Way and other galaxies. △ Less

Submitted 16 September, 2019; originally announced September 2019.

Comments: Accepted in ApJ

Spatial variations in the Milky Way disc metallicity-age relation

Authors: Diane K. Feuillet, Neige Frankel, Karin Lind, Peter M. Frinchaboy, D. A. Garcia-Hernandez, Richard R. Lane, Christian Nitschelm, Alexandre Roman-Lopes

Abstract: Stellar ages are a crucial component to studying the evolution of the Milky Way. Using Gaia DR2 distance estimates, it is now possible to estimate stellar ages for a larger volume of evolved stars through isochrone matching. This work presents [M/H]-age and [$α$/M]-age relations derived for different spatial locations in the Milky Way disc. These relations are derived by hierarchically modelling t… ▽ More Stellar ages are a crucial component to studying the evolution of the Milky Way. Using Gaia DR2 distance estimates, it is now possible to estimate stellar ages for a larger volume of evolved stars through isochrone matching. This work presents [M/H]-age and [$α$/M]-age relations derived for different spatial locations in the Milky Way disc. These relations are derived by hierarchically modelling the star formation history of stars within a given chemical abundance bin. For the first time, we directly observe that significant variation is apparent in the [M/H]-age relation as a function of both Galactocentric radius and distance from the disc mid-plane. The [M/H]-age relations support claims that radial migration has a significant effect in the plane of the disc. Using the [M/H] bin with the youngest mean age at each radial zone in the plane of the disc, the present-day metallicity gradient is measured to be $-0.059 \pm 0.010$ dex kpc$^{-1}$, in agreement with Cepheids and young field stars. We find a vertically flared distribution of young stars in the outer disc, confirming predictions of models and previous observations. The mean age of the [M/H]-[$α$/M] distribution of the solar neighborhood suggests that the high-[M/H] stars are not an evolutionary extension of the low-$α$ sequence. Our observational results are important constraints to Galactic simulations and models of chemical evolution. △ Less

Submitted 7 August, 2019; originally announced August 2019.

Comments: 12 pages, 7 figures, accepted to MNRAS

SanD primes and numbers

Authors: Freeman J. Dyson, Norman E. Frankel, Anthony J. Guttmann

Abstract: We define S(um)anD(ifference) numbers as ordered pairs $(m,\, m+Δ)$ such that the digital-sum $DS(m(m+Δ))=Δ.$ We consider both the decimal and the binary case. If both $m$ and $m+Δ$ are prime numbers, we refer to SanD {\em primes}. We show that the number of (decimal-based) SanD numbers less than $x$ grows as $c1\cdot x,$ where $c1 = 2/3,$ while the number of SanD primes less than $x$ grows as… ▽ More We define S(um)anD(ifference) numbers as ordered pairs $(m,\, m+Δ)$ such that the digital-sum $DS(m(m+Δ))=Δ.$ We consider both the decimal and the binary case. If both $m$ and $m+Δ$ are prime numbers, we refer to SanD {\em primes}. We show that the number of (decimal-based) SanD numbers less than $x$ grows as $c1\cdot x,$ where $c1 = 2/3,$ while the number of SanD primes less than $x$ grows as $c2\cdot x/\log^2{x},$ where $c2 = 3/4.$ Due to the quasi-fractal nature of the digital-sum function, convergence is both slow and erratic compared to twin primes, which, apart from the constant, have the same leading asymptotics. △ Less

Submitted 2 March, 2020; v1 submitted 6 April, 2019; originally announced April 2019.

Comments: 22 pages, 2 figures. Revised version corrects one proof and fixes some typos

MSC Class: 11A41; 11A63; 11Y55; 11Y60

Measuring Radial Orbit Migration in the Milky Way Disk

Authors: Neige Frankel, Hans-Walter Rix, Yuan-Sen Ting, Melissa K. Ness, David W. Hogg

Abstract: We develop and apply a model to quantify the global efficiency of radial orbit migration among stars in the Milky Way disk. This model parameterizes the possible star formation and enrichment histories, radial birth profiles, and combines them with a migration model that relates present-day orbital radii to birth radii through a Gaussian probability, broadening with age $τ$ as… ▽ More We develop and apply a model to quantify the global efficiency of radial orbit migration among stars in the Milky Way disk. This model parameterizes the possible star formation and enrichment histories, radial birth profiles, and combines them with a migration model that relates present-day orbital radii to birth radii through a Gaussian probability, broadening with age $τ$ as $σ_\mathrm{RM8}~\sqrt{τ/8{~\mathrm{Gyr}}}$. Guided by observations, we assume that stars are born with an initially tight age--metallicity relation at given radius, which becomes subsequently scrambled by radial orbit migration, thereby providing a direct observational constraint on radial orbit migration strength $σ_\mathrm{RM8}$. We fit this model with MCMC to the observed age--metallicity distribution of low-$α$ red clump stars with Galactocentric radii between 5 and 14 kpc from APOGEE DR12, sidestepping the complex spatial selection function and accounting for the considerable age uncertainties. This simple model reproduces well the observed data, and we find a global (in radius and time) radial orbit migration efficiency in the Milky Way of $σ_\mathrm{RM8}=3.6\pm 0.1$ kpc when marginalizing over all other model aspects. This shows that radial orbit migration in the Milky Way's main disk is indeed rather strong, in line with theoretical expectations: stars migrate by about a half-mass radius over the age of the disk. The model finds the Sun's birth radius at $\sim 5.2$ kpc. If such strong radial orbit migration is typical, this mechanism plays indeed an important role in setting the structural regularity of disk galaxies. △ Less

Submitted 23 May, 2018; originally announced May 2018.

Comments: Submitted to ApJ, comments are welcome

Non-genetic diversity modulates population performance

Authors: Adam James Waite, Nicholas W. Frankel, Yann S. Dufour, Jessica F. Johnston, Junjiajia Long, Thierry Emonet

Abstract: Biological functions are typically performed by groups of cells that express predominantly the same genes, yet display a continuum of phenotypes. While it is known how one genotype can generate such non-genetic diversity, it remains unclear how different phenotypes contribute to the performance of biological function at the population level. We developed a microfluidic device to simultaneously mea… ▽ More Biological functions are typically performed by groups of cells that express predominantly the same genes, yet display a continuum of phenotypes. While it is known how one genotype can generate such non-genetic diversity, it remains unclear how different phenotypes contribute to the performance of biological function at the population level. We developed a microfluidic device to simultaneously measure the phenotype and chemotactic performance of tens of thousands of individual, freely-swimming Escherichia coli as they climbed a gradient of attractant. We discovered that spatial structure spontaneously emerged from initially well-mixed wild type populations due to non-genetic diversity. By manipulating the expression of key chemotaxis proteins, we established a causal relationship between protein expression, non- genetic diversity, and performance that was theoretically predicted. This approach generated a complete phenotype-to-performance map, in which we found a nonlinear regime. We used this map to demonstrate how changing the shape of a phenotypic distribution can have as large of an effect on collective performance as changing the mean phenotype, suggesting that selection could act on both during the process of adaptation. △ Less

Submitted 3 January, 2017; originally announced January 2017.

Journal ref: Molecular Systems Biology, 12(12), 895 (2016)

Adaptability of non-genetic diversity in bacterial chemotaxis

Authors: Nicholas W Frankel, William Pontius, Yann S Dufour, Junjiajia Long, Luis Hernandez- Nunez, Thierry Emonet

Abstract: Bacterial chemotaxis systems are as diverse as the environments that bacteria inhabit, but how much environmental variation can cells tolerate with a single system? Diversification of a single chemotaxis system could serve as an alternative, or even evolutionary stepping-stone, to switching between multiple systems. We hypothesized that mutations in gene regulation could lead to heritable control… ▽ More Bacterial chemotaxis systems are as diverse as the environments that bacteria inhabit, but how much environmental variation can cells tolerate with a single system? Diversification of a single chemotaxis system could serve as an alternative, or even evolutionary stepping-stone, to switching between multiple systems. We hypothesized that mutations in gene regulation could lead to heritable control of chemotactic diversity. By simulating foraging and colonization of E. coli using a single-cell chemotaxis model, we found that different environments selected for different behaviors. The resulting trade-offs show that populations facing diverse environments would ideally diversify behaviors when time for navigation is limited. We show that advantageous diversity can arise from changes in the distribution of protein levels among individuals, which could occur through mutations in gene regulation. We propose experiments to test our prediction that chemotactic diversity in a clonal population could be a selectable trait that enables adaptation to environmental variability. △ Less

Submitted 25 November, 2015; originally announced November 2015.

Comments: Journal link: http://elifesciences.org/content/3/e03526

Journal ref: eLife 3, e03526 (2014)

Scaling the universe

Authors: Norman E. Frankel

Abstract: A model is presented for the origin of the large scale structure of the universe and their Mass-Radius scaling law; a fractal power law, $M \propto R^D$, with dimension $D=2$, most significantly. The physics is conventional, orthodox, but it is used to fashion a highly unorthodox model of the origin of the galaxies, their groups, clusters, super-clusters, and great walls. The scaling law fits the… ▽ More A model is presented for the origin of the large scale structure of the universe and their Mass-Radius scaling law; a fractal power law, $M \propto R^D$, with dimension $D=2$, most significantly. The physics is conventional, orthodox, but it is used to fashion a highly unorthodox model of the origin of the galaxies, their groups, clusters, super-clusters, and great walls. The scaling law fits the observational results and the model offers new suggestions and predictions. These include a largest, a supreme, cosmic structure, and possible implications for the recently observed pressing cosmological anomalies. △ Less

Submitted 1 November, 2014; v1 submitted 26 January, 2014; originally announced January 2014.

Comments: In v2 and v3 a paragraph has been added at the end of the paper wrt the rippling of the space time fabric; v4/v5 has new parts in Abstract, Introduction and close of Conclusion

Lehmer's Interesting Series

Authors: F. J. Dyson, N. E. Frankel, M. L. Glasser

Abstract: The series $$S_k(z)=\sum_{m=1}^{\infty}\frac{m^kz^m}{(\{array}{c} 2m m \{array})}$$ is evaluated in non-recursive closed and analytically continued beyond its domain of convergence $0\le |z|<4$ for $k=0,1,2,\...$. From this we provide a firm basis for Lehmer's observation that $π$ emerges from the limiting behavior of $S_k(2)$ as $k\rightarrow\infty$. The series $$S_k(z)=\sum_{m=1}^{\infty}\frac{m^kz^m}{(\{array}{c} 2m m \{array})}$$ is evaluated in non-recursive closed and analytically continued beyond its domain of convergence $0\le |z|<4$ for $k=0,1,2,\...$. From this we provide a firm basis for Lehmer's observation that $π$ emerges from the limiting behavior of $S_k(2)$ as $k\rightarrow\infty$. △ Less

Submitted 28 March, 2011; v1 submitted 22 September, 2010; originally announced September 2010.

Comments: 14 pages

Edge effects in some perturbations of the GUE

Authors: K. E. Bassler, P. J. Forrester, N. E. Frankel

Abstract: A bordering of GUE matrices is considered, in which the bordered row consists of zero mean complex Gaussians N$[0,σ/2] + i {\rm N}[0,σ/2]$ off the diagonal, and the real Gaussian N$[μ,σ/\sqrt{2}]$ on the diagonal. We compute the explicit form of the eigenvalue probability function for such matrices, as well as that for matrices obtained by repeating the bordering. The correlations are in general d… ▽ More A bordering of GUE matrices is considered, in which the bordered row consists of zero mean complex Gaussians N$[0,σ/2] + i {\rm N}[0,σ/2]$ off the diagonal, and the real Gaussian N$[μ,σ/\sqrt{2}]$ on the diagonal. We compute the explicit form of the eigenvalue probability function for such matrices, as well as that for matrices obtained by repeating the bordering. The correlations are in general determinantal, and in the single bordering case the explicit form of the correlation kernel is computed. In the large $N$ limit it is shown that $μ$ and/or $σ$ can be tuned to induce a separation of the largest eigenvalue. This effect is shown to be controlled by a single parameter, universal correlation kernel. △ Less

Submitted 18 May, 2010; originally announced May 2010.

Comments: 16 pages, no figures

Eigenvalue Separation in Some Random Matrix Models

Authors: Kevin E. Bassler, Peter J. Forrester, Norman E. Frankel

Abstract: The eigenvalue density for members of the Gaussian orthogonal and unitary ensembles follows the Wigner semi-circle law. If the Gaussian entries are all shifted by a constant amount c/Sqrt(2N), where N is the size of the matrix, in the large N limit a single eigenvalue will separate from the support of the Wigner semi-circle provided c > 1. In this study, using an asymptotic analysis of the secul… ▽ More The eigenvalue density for members of the Gaussian orthogonal and unitary ensembles follows the Wigner semi-circle law. If the Gaussian entries are all shifted by a constant amount c/Sqrt(2N), where N is the size of the matrix, in the large N limit a single eigenvalue will separate from the support of the Wigner semi-circle provided c > 1. In this study, using an asymptotic analysis of the secular equation for the eigenvalue condition, we compare this effect to analogous effects occurring in general variance Wishart matrices and matrices from the shifted mean chiral ensemble. We undertake an analogous comparative study of eigenvalue separation properties when the size of the matrices are fixed and c goes to infinity, and higher rank analogues of this setting. This is done using exact expressions for eigenvalue probability densities in terms of generalized hypergeometric functions, and using the interpretation of the latter as a Green function in the Dyson Brownian motion model. For the shifted mean Gaussian unitary ensemble and its analogues an alternative approach is to use exact expressions for the correlation functions in terms of classical orthogonal polynomials and associated multiple generalizations. By using these exact expressions to compute and plot the eigenvalue density, illustrations of the various eigenvalue separation effects are obtained. △ Less

Submitted 8 October, 2008; originally announced October 2008.

Comments: 25 pages, 9 figures included

Journal ref: J.Math.Phys.50:033302,2009

Analytic solutions of the 1D finite coupling delta function Bose gas

Authors: P. J. Forrester, N. E. Frankel, M. I. Makin

Abstract: An intensive study for both the weak coupling and strong coupling limits of the ground state properties of this classic system is presented. Detailed results for specific values of finite $N$ are given and from them results for general $N$ are determined. We focus on the density matrix and concomitantly its Fourier transform, the occupation numbers, along with the pair correlation function and c… ▽ More An intensive study for both the weak coupling and strong coupling limits of the ground state properties of this classic system is presented. Detailed results for specific values of finite $N$ are given and from them results for general $N$ are determined. We focus on the density matrix and concomitantly its Fourier transform, the occupation numbers, along with the pair correlation function and concomitantly its Fourier transform, the structure factor. These are the signature quantities of the Bose gas. One specific result is that for weak coupling a rational polynomial structure holds despite the transcendental nature of the Bethe equations. All these new results are predicated on the Bethe ansatz and are built upon the seminal works of the past. △ Less

Submitted 11 September, 2006; v1 submitted 21 July, 2006; originally announced July 2006.

Comments: 23 pages, 0 figures, uses rotate.sty. A few lines added. Accepted by Phys. Rev. A

Short Wavelength Analysis of the Evolution of Perturbations in a Two-component Cosmological Fluid

Authors: R. M. Gailis, N. E. Frankel

Abstract: The equations describing a two-component cosmological fluid with linearized density perturbations are investigated in the small wavelength or large $k$ limit. The equations are formulated to include a baryonic component, as well as either a hot dark matter (HDM) or cold dark matter (CDM) component. Previous work done on such a system in static spacetime is extended to reveal some interesting phy… ▽ More The equations describing a two-component cosmological fluid with linearized density perturbations are investigated in the small wavelength or large $k$ limit. The equations are formulated to include a baryonic component, as well as either a hot dark matter (HDM) or cold dark matter (CDM) component. Previous work done on such a system in static spacetime is extended to reveal some interesting physical properties, such as the Jeans wavenumber of the mixture, and resonant mode amplitudes. A WKB technique is then developed to study the expanding universe equations in detail, and to see whether such physical properties are also of relevance in this more realistic scenario. The Jeans wavenumber of the mixture is re-interpreted for the case of an expanding background spacetime. The various modes are obtained to leading order, and the amplitudes of the modes are examined in detail to compare to the resonances observed in the static spacetime results. It is found that some conclusions made in the literature about static spacetime results cannot be carried over to an expanding cosmology. △ Less

Submitted 23 July, 2006; originally announced July 2006.

Comments: 42 pages, 12 figures

Journal ref: J.Math.Phys. 47 (2006) 062506

Two-Component Cosmological Fluids with Gravitational Instabilities

Authors: R. M. Gailis, N. E. Frankel

Abstract: A survey of linearized cosmological fluid equations with a number of different matter components is made. To begin with, the one-component case is reconsidered to illustrate some important mathematical and physical points rarely discussed in the literature. The work of some previous studies of two-component systems are examined and re-analyzed to point out some deficiencies of solutions, and fur… ▽ More A survey of linearized cosmological fluid equations with a number of different matter components is made. To begin with, the one-component case is reconsidered to illustrate some important mathematical and physical points rarely discussed in the literature. The work of some previous studies of two-component systems are examined and re-analyzed to point out some deficiencies of solutions, and further solutions and physical interpretation are then presented. This leads into a general two-component model with variable velocity dispersion parameters and mass density fractions of each component. The equations, applicable to both hot dark matter (HDM) and cold dark matter (CDM) universes are solved in the long wavelength limit. This region is of interest, because some modes in this range of wavenumbers are Jeans unstable. The mixture Jeans wavenumber of the two-component system is introduced and interpreted, and the solutions are discussed, particularly in comparison to analogous solutions previously derived for plasma modes. This work is applicable to that region in the early Universe ($20 < z < 140$), where large scale structure formation is thought to have occurred. △ Less

Submitted 23 July, 2006; originally announced July 2006.

Comments: 36 pages, 2 figures

Journal ref: J.Math.Phys. 47 (2006) 062505

Asymptotic form of the density profile for Gaussian and Laguerre random matrix ensembles with orthogonal and symplectic symmetry

Authors: P. J. Forrester, N. E. Frankel, T. M. Garoni

Abstract: In a recent study we have obtained correction terms to the large N asymptotic expansions of the eigenvalue density for the Gaussian unitary and Laguerre unitary ensembles of random N by N matrices, both in the bulk and at the soft edge of the spectrum. In the present study these results are used to similarly analyze the eigenvalue density for Gaussian and Laguerre random matrix ensembles with or… ▽ More In a recent study we have obtained correction terms to the large N asymptotic expansions of the eigenvalue density for the Gaussian unitary and Laguerre unitary ensembles of random N by N matrices, both in the bulk and at the soft edge of the spectrum. In the present study these results are used to similarly analyze the eigenvalue density for Gaussian and Laguerre random matrix ensembles with orthogonal and symplectic symmetry. As in the case of unitary symmetry, a matching is exhibited between the asymptotic expansion of the bulk density, expanded about the edge, and the asymptotic expansion of the edge density, expanded into the bulk. In addition, aspects of the asymptotic expansion of the smoothed density, which involves delta functions at the endpoints of the support, are interpreted microscopically. △ Less

Submitted 15 August, 2005; originally announced August 2005.

Comments: 28 pages, 1 figure

Asymptotic corrections to the eigenvalue density of the GUE and LUE

Authors: T. M. Garoni, P. J. Forrester, N. E. Frankel

Abstract: We obtain correction terms to the large N asymptotic expansions of the eigenvalue density for the Gaussian unitary and Laguerre unitary ensembles of random N by N matrices, both in the bulk of the spectrum and near the spectral edge. This is achieved by using the well known orthogonal polynomial expression for the kernel to construct a double contour integral representation for the density, to w… ▽ More We obtain correction terms to the large N asymptotic expansions of the eigenvalue density for the Gaussian unitary and Laguerre unitary ensembles of random N by N matrices, both in the bulk of the spectrum and near the spectral edge. This is achieved by using the well known orthogonal polynomial expression for the kernel to construct a double contour integral representation for the density, to which we apply the saddle point method. The main correction to the bulk density is oscillatory in N and depends on the distribution function of the limiting density, while the corrections to the Airy kernel at the soft edge are again expressed in terms of the Airy function and its first derivative. We demonstrate numerically that these expansions are very accurate. A matching is exhibited between the asymptotic expansion of the bulk density, expanded about the edge, and the asymptotic expansion of the edge density, expanded into the bulk. △ Less

Submitted 16 April, 2005; originally announced April 2005.

Comments: 14 pages, 4 figures

Applications and generalizations of Fisher-Hartwig asymptotics

Authors: P. J. Forrester, N. E. Frankel

Abstract: Fisher-Hartwig asymptotics refers to the large $n$ form of a class of Toeplitz determinants with singular generating functions. This class of Toeplitz determinants occurs in the study of the spin-spin correlations for the two-dimensional Ising model, and the ground state density matrix of the impenetrable Bose gas, amongst other problems in mathematical physics. We give a new application of the… ▽ More Fisher-Hartwig asymptotics refers to the large $n$ form of a class of Toeplitz determinants with singular generating functions. This class of Toeplitz determinants occurs in the study of the spin-spin correlations for the two-dimensional Ising model, and the ground state density matrix of the impenetrable Bose gas, amongst other problems in mathematical physics. We give a new application of the original Fisher-Hartwig formula to the asymptotic decay of the Ising correlations above $T_c$, while the study of the Bose gas density matrix leads us to generalize the Fisher-Hartwig formula to the asymptotic form of random matrix averages over the classical groups and the Gaussian and Laguerre unitary matrix ensembles. Another viewpoint of our generalizations is that they extend to Hankel determinants the Fisher-Hartwig asymptotic form known for Toeplitz determinants. △ Less

Submitted 6 January, 2004; originally announced January 2004.

Comments: 25 pages

Random matrix averages and the impenetrable Bose gas in Dirichlet and Neumann boundary conditions

Authors: P. J. Forrester, N. E. Frankel, T. M. Garoni

Abstract: The density matrix for the impenetrable Bose gas in Dirichlet and Neumann boundary conditions can be written in terms of $<\prod_{l=1}^n| \cosφ_1-\cosθ_l| |\cosφ_2-\cosθ_l|>$, where the average is with respect to the eigenvalue probability density function for random unitary matrices from the classical groups $Sp(n)$ and $O^+(2n)$ respectively. In the large $n$ limit log-gas considerations imp… ▽ More The density matrix for the impenetrable Bose gas in Dirichlet and Neumann boundary conditions can be written in terms of $<\prod_{l=1}^n| \cosφ_1-\cosθ_l| |\cosφ_2-\cosθ_l|>$, where the average is with respect to the eigenvalue probability density function for random unitary matrices from the classical groups $Sp(n)$ and $O^+(2n)$ respectively. In the large $n$ limit log-gas considerations imply that the average factorizes into the product of averages of the form $<\prod_{l=1}^n|\cosφ-\cosθ_l>$. By changing variables this average in turn is a special case of the function of $t$ obtained by averaging $\prod_{l=1}^n| t-x_l|^{2q}$ over the Jacobi unitary ensemble from random matrix theory. The latter task is accomplished by a duality formula from the theory of Selberg correlation integrals, and the large $n$ asymptotic form is obtained. The corresponding large $n$ asymptotic form of the density matrix is used, via the exact solution of a particular integral equation, to compute the asymptotic form of the low lying effective single particle states and their occupations, which are proportional to $\sqrt{N}$. △ Less

Submitted 21 May, 2003; v1 submitted 30 January, 2003; originally announced January 2003.

Comments: typos corrected; ref. 7, Phys Rev A67,043607 (2003); accepted for publication in J Math Phys

On the ubiquity of the Levy integral, its relationship with the generalised Euler-Jacobi series, and their asymptotics beyond all orders

Authors: T. M. Garoni, N. E. Frankel

Abstract: We present here an overview of the history, applications and important properties of a function which we refer to as the Levy integral. For certain values of its characteristic parameter the Levy integral defines the symmetric Levy stable probability density function. As we discuss however the Levy integral has applications to a number of other fields besides probability, including random matrix… ▽ More We present here an overview of the history, applications and important properties of a function which we refer to as the Levy integral. For certain values of its characteristic parameter the Levy integral defines the symmetric Levy stable probability density function. As we discuss however the Levy integral has applications to a number of other fields besides probability, including random matrix theory, number theory and asymptotics beyond all orders. We exhibit a direct relationship between the Levy integral and a number theoretic series which we refer to as the generalised Euler-Jacobi series. The complete asymptotic expansions for all natural values of its parameter are presented, and in particular it is pointed out that the intricate exponentially small series become dominant for certain parameter values. △ Less

Submitted 7 November, 2002; originally announced November 2002.

Comments: 7 pages, no figures

Journal ref: Fractals 11, 93-100 (2003)

Finite one dimensional impenetrable Bose systems: Occupation numbers

Authors: P. J. Forrester, N. E. Frankel, T. M. Garoni, N. S. Witte

Abstract: Bosons in the form of ultra cold alkali atoms can be confined to a one dimensional (1d) domain by the use of harmonic traps. This motivates the study of the ground state occupations $λ_i$ of effective single particle states $φ_i$, in the theoretical 1d impenetrable Bose gas. Both the system on a circle and the harmonically trapped system are considered. The $λ_i$ and $φ_i$ are the eigenvalues an… ▽ More Bosons in the form of ultra cold alkali atoms can be confined to a one dimensional (1d) domain by the use of harmonic traps. This motivates the study of the ground state occupations $λ_i$ of effective single particle states $φ_i$, in the theoretical 1d impenetrable Bose gas. Both the system on a circle and the harmonically trapped system are considered. The $λ_i$ and $φ_i$ are the eigenvalues and eigenfunctions respectively of the one body density matrix. We present a detailed numerical and analytic study of this problem. Our main results are the explicit scaled forms of the density matrices, from which it is deduced that for fixed $i$ the occupations $λ_i$ are asymptotically proportional to $\sqrt{N}$ in both the circular and harmonically trapped cases. △ Less

Submitted 6 November, 2002; originally announced November 2002.

Comments: 22 pages, 8 figures (.eps), uses REVTeX4

Painlevé transcendent evaluations of finite system density matrices for 1d impenetrable Bosons

Authors: P. J. Forrester, N. E. Frankel, T. M. Garoni, N. S. Witte

Abstract: The recent experimental realisation of a one-dimensional Bose gas of ultra cold alkali atoms has renewed attention on the theoretical properties of the impenetrable Bose gas. Of primary concern is the ground state occupation of effective single particle states in the finite system, and thus the tendency for Bose-Einstein condensation. This requires the computation of the density matrix. For the… ▽ More The recent experimental realisation of a one-dimensional Bose gas of ultra cold alkali atoms has renewed attention on the theoretical properties of the impenetrable Bose gas. Of primary concern is the ground state occupation of effective single particle states in the finite system, and thus the tendency for Bose-Einstein condensation. This requires the computation of the density matrix. For the impenetrable Bose gas on a circle we evaluate the density matrix in terms of a particular Painlevé VI transcendent in $σ$-form, and furthermore show that the density matrix satisfies a recurrence relation in the number of particles. For the impenetrable Bose gas in a harmonic trap, and with Dirichlet or Neumann boundary conditions, we give a determinant form for the density matrix, a form as an average over the eigenvalues of an ensemble of random matrices, and in special cases an evaluation in terms of a transcendent related to Painlevé V and VI. We discuss how our results can be used to compute the ground state occupations. △ Less

Submitted 26 January, 2003; v1 submitted 1 July, 2002; originally announced July 2002.

Comments: AMSLatex, requires Commun. Math. Phys. style macros

MSC Class: 81Hxx;15A52;33E17

The Princess and the Pea

Authors: Neil Cornish, Norm Frankel

Abstract: Like a fairy-tale princess, trajectories around black holes can be sensitive to small disturbances. We describe how a small disturbance can lead to erratic orbits and an increased production of gravitational waves. Like a fairy-tale princess, trajectories around black holes can be sensitive to small disturbances. We describe how a small disturbance can lead to erratic orbits and an increased production of gravitational waves. △ Less

Submitted 28 January, 1997; v1 submitted 11 June, 1996; originally announced June 1996.

Comments: 5 pages, RevTeX, 6 figures, revised and expanded version

Report number: UM-P-96/45

Journal ref: Phys.Rev. D56 (1997) 1903-1907

Chaos and Fractals around Black Holes

Authors: C. Dettmann, N. Frankel, N. Cornish

Abstract: Fractal basin boundaries provide an important means of characterizing chaotic systems. We apply these ideas to general relativity, where other properties such as Lyapunov exponents are difficult to define in an observer independent manner. Here we discuss the difficulties in describing chaotic systems in general relativity and investigate the motion of particles in two- and three-black-hole spacet… ▽ More Fractal basin boundaries provide an important means of characterizing chaotic systems. We apply these ideas to general relativity, where other properties such as Lyapunov exponents are difficult to define in an observer independent manner. Here we discuss the difficulties in describing chaotic systems in general relativity and investigate the motion of particles in two- and three-black-hole spacetimes. We show that the dynamics is chaotic by exhibiting the basins of attraction of the black holes which have fractal boundaries. Overcoming problems of principle as well as numerical difficulties, we evaluate Lyapunov exponents numerically and find that some trajectories have a positive exponent. △ Less

Submitted 6 February, 1995; originally announced February 1995.

Comments: To appear in "Fractals" March issue (World Scientific), 20 figures available by request, also available from SLAC's gr-qc postscript archive

Report number: UTPT-94-22, UM-P-94/87

Journal ref: Fractals 3 (1995) 161

Fractal basins and chaotic trajectories in multi-black hole space-times

Authors: N. J. Cornish, C. P. Dettmann, N. E. Frankel

Abstract: We investigate the phase-space for trajectories in multi-black hole spacetimes. We find that complete, chaotic geodesics are well described by Lyapunov exponents, and that the attractor basin boundary scales as a fractal in a diffeomorphism invariant manner. We investigate the phase-space for trajectories in multi-black hole spacetimes. We find that complete, chaotic geodesics are well described by Lyapunov exponents, and that the attractor basin boundary scales as a fractal in a diffeomorphism invariant manner. △ Less

Submitted 14 February, 1994; originally announced February 1994.

Comments: 9 pages, RevTeX, 4 figures (appended as a compressed tar file), UTPT-93-32, UM-P-94/18

Journal ref: Phys.Rev.D50:618-621,1994

Gravitation vs. Rotation in 2+1 Dimensions

Authors: N. J. Cornish, N. E. Frankel

Abstract: We investigate rotation and rotating structures in (2+1)-dimensional Einstein gravity. We show that rotation generally leads to pathological physical situations. We investigate rotation and rotating structures in (2+1)-dimensional Einstein gravity. We show that rotation generally leads to pathological physical situations. △ Less

Submitted 10 September, 1993; v1 submitted 8 June, 1993; originally announced June 1993.

Comments: 11 pages, UM-P-93/43 & UTPT-93-10

Journal ref: Class.Quant.Grav.11:723,1994