Search

The Large Magellanic Cloud (LMC) is a Milky Way (MW) satellite that is massive enough to gravitationally attract the MW disc and inner halo, causing significant motion of the inner MW with respect to the outer halo. In this work, we probe this interaction by constructing a sample of 9,866 blue horizontal branch (BHB) stars with radial velocities from the DESI spectroscopic survey out to 120 kpc from the Galactic centre. This is the largest spectroscopic set of BHB stars in the literature to date, and it contains four times more stars with Galactocentric distances beyond 50 kpc than previous BHB catalogues. Using the DESI BHB sample combined with SDSS BHBs, we measure the bulk radial velocity of stars in the outer halo and observe that the velocity in the Southern Galactic hemisphere is different by 3.7$\sigma$ from the North. Modelling the projected velocity field shows that its dipole component is directed at a point 22 degrees away from the LMC along its orbit, which we interpret as the travel direction of the inner MW. The velocity field includes a monopole term that is -24 km/s, which we refer to as compression velocity. This velocity is significantly larger than predicted by the current models of the MW and LMC interaction. This work uses DESI data from its first two years of observations, but we expect that with upcoming DESI data releases, the sample of BHB stars will increase and our ability to measure the MW-LMC interaction will improve significantly., Comment: 22 pages, 19 figures. Submitted to MNRAS

Barred galaxies exhibit boxy/peanut or X-shapes (BP/X) protruding from their disks in edge-on views. Two types of BP/X morphologies exist depending on whether the X-wings meet at the center (CX) or are off-centered (OX). Orbital studies indicate that various orbital types can generate X-shaped structures. Here, we provide a classification approach that identifies the specific orbit families responsible for generating OX and CX-shaped structures. Applying this approach to three different N-body bar models, we show that both OX and CX structures are associated with the x1 orbit family, but OX-supporting orbits possess higher angular momentum (closer to x1 orbits) than orbits in CX structures. Consequently, as the bar slows down the contribution of higher angular momentum OX-supporting orbits decreases and that of lower angular momentum orbits increases resulting in an evolution of the morphology from OX to CX. If the bar does not slow down, the shape of the BP/X structure and the fractions of OX/CX supporting orbits remain substantially unchanged. Bars that do not undergo buckling but that do slow down initially show the OX structure and are dominated by high angular momentum orbits, transitioning to a CX morphology. Bars that buckle exhibit a combination of both OX and CX supporting orbits immediately after the buckling, but become more CX dominated as their pattern speed decreases. This study demonstrates that the evolution of BP/X morphology and orbit populations strongly depends on the evolution of the bar angular momentum., Comment: 16 pages, 12 figures, Accepted for publication in ApJ

Modern artificial intelligence (AI) systems are powered by foundation models. This paper presents a new set of foundation models, called Llama 3. It is a herd of language models that natively support multilinguality, coding, reasoning, and tool usage. Our largest model is a dense Transformer with 405B parameters and a context window of up to 128K tokens. This paper presents an extensive empirical evaluation of Llama 3. We find that Llama 3 delivers comparable quality to leading language models such as GPT-4 on a plethora of tasks. We publicly release Llama 3, including pre-trained and post-trained versions of the 405B parameter language model and our Llama Guard 3 model for input and output safety. The paper also presents the results of experiments in which we integrate image, video, and speech capabilities into Llama 3 via a compositional approach. We observe this approach performs competitively with the state-of-the-art on image, video, and speech recognition tasks. The resulting models are not yet being broadly released as they are still under development.

The Southern Photometric Local Universe Survey (S-PLUS) is a project to map $\sim9300$ sq deg of the sky using twelve bands (seven narrow and five broadbands). Observations are performed with the T80-South telescope, a robotic telescope located at the Cerro Tololo Observatory in Chile. The survey footprint consists of several large contiguous areas, including fields at high and low galactic latitudes, and towards the Magellanic Clouds. S-PLUS uses fixed exposure times to reach point source depths of about $21$ mag in the $griz$ and $20$ mag in the $u$ and the narrow filters. This paper describes the S-PLUS Data Release 4 (DR4), which includes calibrated images and derived catalogues for over 3000 sq deg, covering the aforementioned area. The catalogues provide multi-band photometry performed with the tools \texttt{DoPHOT} and \texttt{SExtractor} -- point spread function (\PSF) and aperture photometry, respectively. In addition to the characterization, we also present the scientific potential of the data. We use statistical tools to present and compare the photometry obtained through different methods. Overall we find good agreement between the different methods, with a slight systematic offset of 0.05\,mag between our \PSF and aperture photometry. We show that the astrometry accuracy is equivalent to that obtained in previous S-PLUS data releases, even in very crowded fields where photometric extraction is challenging. The depths of main survey (MS) photometry for a minimum signal-to-noise ratio $S/N = 3$ reach from $\sim19.5$ for the bluer bands to $\sim21.5$ mag on the red. The range of magnitudes over which accurate \PSF photometry is obtained is shallower, reaching $\sim19$ to $\sim20.5$ mag depending on the filter. Based on these photometric data, we provide star-galaxy-quasar classification and photometric redshift for millions of objects., Comment: 26 pages, 17 figures, 14 tables, accepted for A&A

A tracer sample in a gravitational potential, starting from a generic initial condition, phase-mixes towards a stationary state. This evolution is accompanied by an entropy increase, and the final state is characterized by a distribution function (DF) that depends only on integrals of motion (Jeans theorem). We present a method to constrain a gravitational potential where a sample is stationary by minimizing the entropy the sample would have if it were allowed to phase-mix in trial potentials. This method avoids assuming a known DF, and is applicable to any sets of integrals. We provide expressions for the entropy of DFs depending on energy, $f(E)$, energy and angular momentum, $f(E,L)$, or three actions, $f(\vec{J})$, and investigate the bias and fluctuations in their estimates. We show that the method correctly recovers the potential parameters for spherical and axisymmetric models. We also present a methodology to characterize the posterior probability distribution of the parameters with an Approximate Bayesian Computation, indicating a pathway for application to observational data. Using $N=10^4$ tracers with $20\%$-uncertainties in the 6D coordinates, we recover the flattening parameter $q$ of an axisymmetric potential with $\sigma_q/q\sim 10\%$., Comment: Significant changes in appendix with mathematical results. Final expressions and numerical results are the same. A couple of discussions subsections added, and a new appendix on the possibility of maximizing an entropy in angle-space. Submitted to ApJ

We present ~ 126 new spectroscopically identified members of the GD-1 tidal stream obtained with the 5000-fiber Dark Energy Spectroscopic Instrument (DESI). We confirm the existence of a ``cocoon'' which is broad (FWHM~2.932deg~460pc) and kinematically hot (velocity dispersion, sigma~5-8km/s) component that surrounds a narrower (FWHM~0.353deg~55pc) and colder (sigma~ 2.2-2.6km/s) thin stream component (based on a median per star velocity precision of 2.7km/s). The cocoon extends over at least a ~ 20deg segment of the stream observed by DESI. The thin and cocoon components have similar mean values of [Fe/H]: -2.54+/- 0.04dex and -2.45+/-0.06dex suggestive of a common origin. The data are consistent with the following scenarios for the origin of the cocoon. The progenitor of the GD-1 stream was an accreted globular cluster (GC) and: (a) the cocoon was produced by pre-accretion tidal stripping of the GC while it was still inside its parent dwarf galaxy; (b) the cocoon is debris from the parent dwarf galaxy; (c) an initially thin GC tidal stream was heated by impacts from dark subhalos in the Milky Way; (d) an initially thin GC stream was heated by a massive Sagittarius dwarf galaxy; or a combination of some these. In the first two cases the velocity dispersion and mean metallicity are consistent with the parent dwarf galaxy having a halo mass of ~0^9\msun. Future DESI spectroscopy and detailed modeling may enable us to distinguish between these possible origins., Comment: Submitted to ApJ, 23 pages, 13 figures 4 tables

This paper presents the first public data release of the S-PLUS Ultra-Short Survey (USS), a photometric survey with short exposure times, covering approximately 9300 deg$^{2}$ of the Southern sky. The USS utilizes the Javalambre 12-band magnitude system, including narrow and medium-band and broad-band filters targeting prominent stellar spectral features. The primary objective of the USS is to identify bright, extremely metal-poor (EMP; [Fe/H] $\leq -3$) and ultra metal-poor (UMP; [Fe/H] $\leq -4$) stars for further analysis using medium- and high-resolution spectroscopy.}{This paper provides an overview of the survey observations, calibration method, data quality, and data products. Additionally, it presents the selection of EMP and UMP candidates.}{The data from the USS were reduced and calibrated using the same methods as presented in the S-PLUS DR2. An additional step was introduced, accounting for the offset between the observed magnitudes off the USS and the predicted magnitudes from the very low-resolution Gaia XP spectra.}{This first release contains data for 163 observed fields totaling $\sim$324 deg$^{2}$ along the Celestial Equator. The magnitudes obtained from the USS are well-calibrated, showing a difference of $\sim 15$ mmag compared to the predicted magnitudes by the GaiaXPy toolkit. By combining colors and magnitudes, 140 candidates for EMP or UMP have been identified for follow-up studies.}{The S-PLUS USS DR1 is an important milestone in the search for bright metal-poor stars, with magnitudes in the range 10 $ < r \leq 14$. The USS is an ongoing survey; in the near future, it will provide many more bright metal-poor candidate stars for spectroscopic follow-up., Comment: Accepted for publication in A&A. 17 pages 6 figures. Long table at the end of the paper

We study azimuthal variations in the mean stellar metallicity, <[Fe/H]>, in a self-consistent, isolated simulation in which all stars form out of gas. We find <[Fe/H]> variations comparable to those observed in the Milky Way and which are coincident with the spiral density waves. The azimuthal variations are present in young and old stars and therefore are not a result of recently formed stars. Similar variations are present in the mean age and alpha-abundance. We measure the pattern speeds of the <[Fe/H]>-variations and find that they match those of the spirals, indicating that they are at the origin of the metallicity patterns. Because younger stellar populations are not just more [Fe/H]-rich and alpha-poor but also dynamically cooler, we expect them to more strongly support spirals, which is indeed the case in the simulation. However, if we measure the radial action, J_R, using the Stackel axisymmetric approximation, we find that the spiral ridges are traced by regions of high J_R, contrary to expectations. Assuming that the passage of stars through the spirals leads to unphysical variations in the measured J_R, we obtain an improved estimate of J_R by averaging over a 1 Gyr time interval. This time-averaged J_R is a much better tracer of the spiral structure, with minima at the spiral ridges. We conclude that the errors incurred by the axisymmetric approximation introduce correlated deviations large enough to render the instantaneous radial action inadequate for tracing spirals., Comment: 21 pages. Submitted to MNRAS. Comments welcomed

The classification of phases and the detection of phase transitions are central and challenging tasks in diverse fields. Within physics, it relies on the identification of order parameters and the analysis of singularities in the free energy and its derivatives. Here, we propose an alternative framework to identify quantum phase transitions. Using the axial next-nearest neighbor Ising (ANNNI) model as a benchmark, we show how machine learning can detect three phases (ferromagnetic, paramagnetic, and a cluster of the antiphase with the floating phase). Employing supervised learning, we demonstrate the feasibility of transfer learning. Specifically, a machine trained only with nearest-neighbor interactions can learn to identify a new type of phase occurring when next-nearest-neighbor interactions are introduced. We also compare the performance of common classical machine learning methods with a version of the quantum nearest neighbors (QNN) algorithm.

The Milky Way stellar disk has both a thin and a thick component. The thin disk is composed mostly of younger stars ($\lesssim$8 Gyr) with a lower abundance of $\alpha$ elements, while the thick disk contains predominantly older stars ($\gtrsim$8--12 Gyr) with a higher $\alpha$ abundance, giving rise to an $\alpha$-bimodality most prominent at intermediate metallicities. A proposed explanation for the bimodality is an episode of clumpy star formation, where high-$\alpha$ stars form in massive clumps that appear in the first few Gyrs of the Milky Way's evolution, while low-$\alpha$ stars form throughout the disk and over a longer time span. To better understand the evolution of clumps, we track them and their constituent stars in two clumpy Milky Way simulations that reproduce the $\alpha$-abundance bimodality, one with 10% and the other with 20% supernova feedback efficiency. We investigate the paths that these clumps take in the chemical space ([O/Fe]--[Fe/H]) as well as their mass, star formation rate (SFR), formation location, lifetime, and merger history. The clumps in the simulation with lower feedback last longer on average, with several lasting hundreds of Myr. Some of the clumps do not reach high-$\alpha$, but the ones that do on average had a higher SFR, longer lifetime, greater mass, and form closer to the galactic center than the ones that do not. Most clumps that reach high-$\alpha$ merge with others and eventually spiral into the galactic center, but shed stars along the way to form most of the thick disk component., Comment: 14 pages, 14 figures, accepted for publication in ApJ

As our nearest large neighbor, the Andromeda Galaxy provides a unique laboratory for investigating galaxy formation and the distribution and substructure properties of dark matter in a Milky Way-like galaxy. Here, we propose an initial 2-epoch ($\Delta t\approx 5$yr), 2-band Roman survey of the entire halo of Andromeda, covering 500 square degrees, which will detect nearly every red giant star in the halo (10$\sigma$ detection in F146, F062 of 26.5, 26.1AB mag respectively) and yield proper motions to $\sim$25 microarcsec/year (i.e., $\sim$90 km/s) for all stars brighter than F146 $\approx 23.6$ AB mag (i.e., reaching the red clump stars in the Andromeda halo). This survey will yield (through averaging) high-fidelity proper motions for all satellites and compact substructures in the Andromeda halo and will enable statistical searches for clusters in chemo-dynamical space. Adding a third epoch during the extended mission will improve these proper motions by $\sim t^{-1.5}$, to $\approx 11$ km/s, but this requires obtaining the first epoch in Year 1 of Roman operations. In combination with ongoing and imminent spectroscopic campaigns with ground-based telescopes, this Roman survey has the potential to yield full 3-d space motions of $>$100,000 stars in the Andromeda halo, including (by combining individual measurements) robust space motions of its entire globular cluster and most of its dwarf galaxy satellite populations. It will also identify high-velocity stars in Andromeda, providing unique information on the processes that create this population. These data offer a unique opportunity to study the immigration history, halo formation, and underlying dark matter scaffolding of a galaxy other than our own., Comment: Submitted in response to the call for Roman Space Telescope Core Community Survey white papers

The Nancy Grace Roman Space Telescope is capable of delivering an unprecedented all-sky, high-spatial resolution, multi-epoch infrared map to the astronomical community. This opportunity arises in the midst of numerous ground- and space-based surveys that will provide extensive spectroscopy and imaging together covering the entire sky (such as Rubin/LSST, Euclid, UNIONS, SPHEREx, DESI, SDSS-V, GALAH, 4MOST, WEAVE, MOONS, PFS, UVEX, NEO Surveyor, etc.). Roman can uniquely provide uniform high-spatial-resolution (~0.1 arcsec) imaging over the entire sky, vastly expanding the science reach and precision of all of these near-term and future surveys. This imaging will not only enhance other surveys, but also facilitate completely new science. By imaging the full sky over two epochs, Roman can measure the proper motions for stars across the entire Milky Way, probing 100 times fainter than Gaia out to the very edge of the Galaxy. Here, we propose NANCY: a completely public, all-sky survey that will create a high-value legacy dataset benefiting innumerable ongoing and forthcoming studies of the universe. NANCY is a pure expression of Roman's potential: it images the entire sky, at high spatial resolution, in a broad infrared bandpass that collects as many photons as possible. The majority of all ongoing astronomical surveys would benefit from incorporating observations of NANCY into their analyses, whether these surveys focus on nearby stars, the Milky Way, near-field cosmology, or the broader universe., Comment: Submitted to the call for white papers for the Roman Core Community Survey (June 16th, 2023), and to the Bulletin of the AAS

Using N-body simulations we explore the effects of growing a supermassive black hole (SMBH) prior to or during the formation of a stellar bar. Keeping the final mass and growth rate of the SMBH fixed, we show that if it is introduced before or while the bar is still growing, the SMBH does not cause a decrease in bar amplitude. Rather, in most cases, it is strengthened. In addition early growing SMBHs always either decreases the buckling amplitude, delay buckling, or both. This weakening of buckling is caused by an increase in the disk vertical velocity dispersion at radii well beyond the nominal black hole sphere-of-influence. While we find considerable stochasticity and sensitivity to initial conditions, the only case where the SMBH causes a decrease in bar amplitude is when it is introduced after the bar has attained a steady state. In this case we confirm previous findings that the decrease in bar strength is a result of scattering of bar-supporting orbits with small pericenter radii. By heating the inner disk both radially and vertically, an early growing SMBH increases the fraction of stars that can be captured by the Inner Lindblad Resonance (ILR) and the vertical ILR, thereby strengthening both the bar and the boxy peanut shaped bulge. Using orbital frequency analysis of star particles, we show that when an SMBH is introduced early and the bar forms around it, the bar is populated by different families of regular bar-supporting orbits than when the bar forms without an SMBH., Comment: 26 pages, 15 figures. Updated figures and some arguments in Sec. 5.3,5.4, Results Unchanged. Updates reflect format and style of final published version

Metal-poor stars are key to our understanding of the early stages of chemical evolution in the Universe. New multi-filter surveys, such as the Southern Photometric Local Universe Survey (S-PLUS), are greatly advancing our ability to select low-metallicity stars. In this work, we analyse the chemodynamical properties and ages of 522 metal-poor candidates selected from the S-PLUS data release 3. About 92% of these stars were confirmed to be metal-poor ([Fe/H] $\leq -1$) based on previous medium-resolution spectroscopy. We calculated the dynamical properties of a subsample containing 241 stars, using the astrometry from Gaia Data Release 3. Stellar ages are estimated by a Bayesian isochronal method formalized in this work. We analyse the metallicity distribution of these metal-poor candidates separated into different subgroups of total velocity, dynamical properties, and ages. Our results are used to propose further restrictions to optimize the selection of metal-poor candidates in S-PLUS. The proposed astrometric selection ($\mathrm{parallax}>0.85$ mas) is the one that returns the highest fraction of extremely metal-poor stars (16.3% have [Fe/H] $\leq -3$); the combined selection provides the highest fraction of very metal-poor stars (91.0% have [Fe/H] $\leq -2$), whereas the dynamical selection (eccentricity > 0.35 and diskness < 0.75) is better for targetting metal-poor (99.5% have [Fe/H] $\leq -1$). Using only S-PLUS photometric selections, it is possible to achieve selection fractions of 15.6%, 88.5% and 98.3% for metallicities below $-$3, $-$2 and $-$1, respectively. We also show that it is possible to use S-PLUS to target metal-poor stars in halo substructures such as Gaia-Sausage/Enceladus, Sequoia, Thamnos and the Helmi stream., Comment: 18 pages, 13 figures. To be published in MNRAS main journal (accepted 15-may-2023)

In Paper I we showed that clumps in high-redshift galaxies, having a high star formation rate density (\Sigma_SFR), produce disks with two tracks in the [Fe/H]-[\alpha/Fe] chemical space, similar to that of the Milky Way's (MW's) thin + thick disks. Here we investigate the effect of clumps on the bulge's chemistry. The chemistry of the MW's bulge is comprised of a single track with two density peaks separated by a trough. We show that the bulge chemistry of an N-body + smoothed particle hydrodynamics clumpy simulation also has a single track. Star formation within the bulge is itself in the high-\Sigma_SFR clumpy mode, which ensures that the bulge's chemical track follows that of the thick disk at low [Fe/H] and then extends to high [Fe/H], where it peaks. The peak at low metallicity instead is comprised of a mixture of in-situ stars and stars accreted via clumps. As a result, the trough between the peaks occurs at the end of the thick disk track. We find that the high-metallicity peak dominates near the mid-plane and declines in relative importance with height, as in the MW. The bulge is already rapidly rotating by the end of the clump epoch, with higher rotation at low [\alpha/Fe]. Thus clumpy star formation is able to simultaneously explain the chemodynamic trends of the MW's bulge, thin + thick disks and the Splash., Comment: 25 pages, in press at ApJ

The cyclicity problem, crucial in analyzing planar vector fields, consists in estimating the number of limit cycles emanating from monodromic singularities. Traditionally, this estimation relies on Lyapunov coefficients. However, in nonsmooth systems, besides the limit cycles bifurcating by varying the Lyapunov coefficients, monodromic singularities on the switching curve can always be split apart yielding, under suitable conditions, a sliding region and an additional limit cycle surrounding it. This bifurcation phenomenon, known as pseudo-Hopf bifurcation, has enhanced lower-bound cyclicity estimations for monodromic singularities in Filippov systems. In this study, we push beyond the pseudo-Hopf bifurcation, demonstrating that the destruction of $(2k,2k)$-monodromic tangential singularities yields at least $k$ limit cycles surrounding sliding segments. This new bifurcation phenomenon expands our understanding of limit cycle bifurcations in nonsmooth systems and, in addition to the theoretical significance, has practical relevance in various applied models involving switches and abrupt processes.

Many barred galaxies exhibit upturns (shoulders) in their bar major-axis density profile. Simulation studies have suggested that shoulders are supported by looped $x_1$ orbits, occur in growing bars, and can appear after bar-buckling. We investigate the orbital support and evolution of shoulders via frequency analyses of orbits in simulations. We confirm that looped orbits are shoulder-supporting, and can remain so, to a lesser extent, after being vertically thickened. We show that looped orbits appear at the resonance $(\Omega_\varphi - \Omega_\mathrm{P})/\Omega_R=1/2$ (analogous to the classical Inner Lindblad Resonance, and here called ILR) with vertical-to-radial frequency ratios $1 \lesssim\Omega_z/\Omega_R \lesssim 3/2$ (vertically warm orbits). Cool orbits at the ILR (those with $\Omega_z/\Omega_R > 3/2$) are vertically thin and have no loops, contributing negligibly to shoulders. As bars slow and thicken, either secularly or by buckling, they populate warm orbits at the ILR. Further thickening carries these orbits towards crossing the vertical ILR [vILR, $(\Omega_\varphi - \Omega_\mathrm{P})/\Omega_z=1/2$], where they convert in-plane to vertical motion, become chaotic, kinematically hotter and less shoulder-supporting. Hence, persistent shoulders require bars to trap new stars, consistent with the need for a growing bar. Since buckling speeds up trapping on warm orbits at the ILR, it can be followed by shoulder formation, as seen in simulations. This sequence supports the recent observational finding that shoulders likely precede the emergence of BP-bulges. The python module for the frequency analysis, naif, is made available., Comment: Published in ApJ; A few additional figures in response to referee's request. Conclusions unchanged

We measure and compare the pattern speeds of vertical breathing, vertical bending, and spiral density waves in two isolated N-body+SPH simulations, using windowed Fourier transforms over 1 Gyr time intervals. We show that the pattern speeds of the breathing waves match those of the spirals but are different from those of the bending waves. We also observe matching pattern speeds between the bar and breathing waves. Our results not only strengthen the case that, throughout the disc, breathing motions are driven by spirals but indeed that the breathing motions are part and parcel of the spirals., Comment: 5 pages, 4 figures, published in MNRAS Letters

We describe the Milky Way Survey (MWS) that will be undertaken with the Dark Energy Spectroscopic Instrument (DESI) on the Mayall 4m telescope at the Kitt Peak National Observatory. Over the next 5 yr DESI MWS will observe approximately seven million stars at Galactic latitudes |b|>20 degrees, with an inclusive target selection scheme focused on the thick disk and stellar halo. MWS will also include several high-completeness samples of rare stellar types, including white dwarfs, low-mass stars within 100pc of the Sun, and horizontal branch stars. We summarize the potential of DESI to advance understanding of Galactic structure and stellar evolution. We introduce the final definitions of the main MWS target classes and estimate the number of stars in each class that will be observed. We describe our pipelines for deriving radial velocities, atmospheric parameters, and chemical abundances. We use ~500,000 spectra of unique stellar targets from the DESI Survey Validation program (SV) to demonstrate that our pipelines can measure radial velocities to ~1 km/s and [Fe/H] accurate to ~0.2 dex for typical stars in our main sample. We find the stellar parameter distributions from ~100 sq. deg of SV observations with >90% completeness on our main sample are in good agreement with expectations from mock catalogs and previous surveys., Comment: Accepted for publication in ApJ, 44 pages, 25 figures, 4 tables, one of a suite of 8 papers detailing targeting for DESI; v2 added links to data shown in figures, added citations to other DESI papers, corrected author list and minor typos; v3 fixed minor errors in Fig. 6 and clarified associated text; v4 updated to include minor changes in response to review