Showing 1–43 of 43 results for author: Muñoz, D

HATS-38 b and WASP-139 b join a growing group of eccentric hot Neptunes on polar orbits

Authors: Juan I. Espinoza-Retamal, Guðmundur Stefánsson, Cristobal Petrovich, Rafael Brahm, Andrés Jordán, Elyar Sedaghati, Jennifer P. Lucero, Marcelo Tala Pinto, Diego J. Muñoz, Gavin Boyle, Rodrigo Leiva, Vincent Suc

Abstract: We constrain the sky-projected obliquities of two low-density hot Neptune planets HATS-38 b and WASP-139 b orbiting nearby G and K stars using Rossiter-McLaughlin (RM) observations with VLT/ESPRESSO, yielding $λ= -108_{-16}^{+11}$ deg and $-85.6_{-4.2}^{+7.7}$ deg, respectively. To model the RM effect, we use a new publicly available code, \texttt{ironman}, which is capable of jointly fitting tran… ▽ More We constrain the sky-projected obliquities of two low-density hot Neptune planets HATS-38 b and WASP-139 b orbiting nearby G and K stars using Rossiter-McLaughlin (RM) observations with VLT/ESPRESSO, yielding $λ= -108_{-16}^{+11}$ deg and $-85.6_{-4.2}^{+7.7}$ deg, respectively. To model the RM effect, we use a new publicly available code, \texttt{ironman}, which is capable of jointly fitting transit photometry, Keplerian radial velocities, and RM effects. The two planets have residual eccentricities ($e= 0.112_{-0.070}^{+0.072}$, and $0.103_{-0.041}^{+0.050}$, respectively), and together with the obliquity constraints, we show that they join a growing group of eccentric hot and low-density Neptunes on polar orbits. We use long-term radial velocities to rule out companions with masses $\sim 0.3-50$ $M_J$ within $\sim$10 au. We show that the orbital architectures of the two Neptunes disfavor an origin from primordial disk misalignment and/or in-situ secular interactions with distant companions and instead favor high-eccentricity migration from $\gtrsim 2$ au driven by a distant companion. Finally, we performed a hierarchical Bayesian modeling of the true obliquity distribution of Neptunes and found suggestive evidence for a higher preponderance of polar orbits of hot Neptunes compared to Jupiters. However, we note that the exact distribution is sensitive to the choice of priors, highlighting the need for additional obliquity measurements of Neptunes to robustly compare the hot Neptune obliquity distribution to Jupiters. △ Less

Submitted 26 June, 2024; originally announced June 2024.

Comments: 16 pages, 5 figures. Submitted to AAS Journals

Gaia Focused Product Release: Sources from Service Interface Function image analysis -- Half a million new sources in omega Centauri

Authors: Gaia Collaboration, K. Weingrill, A. Mints, J. Castañeda, Z. Kostrzewa-Rutkowska, M. Davidson, F. De Angeli, J. Hern��ndez, F. Torra, M. Ramos-Lerate, C. Babusiaux, M. Biermann, C. Crowley, D. W. Evans, L. Lindegren, J. M. Martín-Fleitas, L. Palaversa, D. Ruz Mieres, K. Tisanić, A. G. A. Brown, A. Vallenari, T. Prusti, J. H. J. de Bruijne, F. Arenou, A. Barbier , et al. (378 additional authors not shown)

Abstract: Gaia's readout window strategy is challenged by very dense fields in the sky. Therefore, in addition to standard Gaia observations, full Sky Mapper (SM) images were recorded for nine selected regions in the sky. A new software pipeline exploits these Service Interface Function (SIF) images of crowded fields (CFs), making use of the availability of the full two-dimensional (2D) information. This ne… ▽ More Gaia's readout window strategy is challenged by very dense fields in the sky. Therefore, in addition to standard Gaia observations, full Sky Mapper (SM) images were recorded for nine selected regions in the sky. A new software pipeline exploits these Service Interface Function (SIF) images of crowded fields (CFs), making use of the availability of the full two-dimensional (2D) information. This new pipeline produced half a million additional Gaia sources in the region of the omega Centauri ($ω$ Cen) cluster, which are published with this Focused Product Release. We discuss the dedicated SIF CF data reduction pipeline, validate its data products, and introduce their Gaia archive table. Our aim is to improve the completeness of the {\it Gaia} source inventory in a very dense region in the sky, $ω$ Cen. An adapted version of {\it Gaia}'s Source Detection and Image Parameter Determination software located sources in the 2D SIF CF images. We validated the results by comparing them to the public {\it Gaia} DR3 catalogue and external Hubble Space Telescope data. With this Focused Product Release, 526\,587 new sources have been added to the {\it Gaia} catalogue in $ω$ Cen. Apart from positions and brightnesses, the additional catalogue contains parallaxes and proper motions, but no meaningful colour information. While SIF CF source parameters generally have a lower precision than nominal {\it Gaia} sources, in the cluster centre they increase the depth of the combined catalogue by three magnitudes and improve the source density by a factor of ten. This first SIF CF data publication already adds great value to the {\it Gaia} catalogue. It demonstrates what to expect for the fourth {\it Gaia} catalogue, which will contain additional sources for all nine SIF CF regions. △ Less

Submitted 8 November, 2023; v1 submitted 10 October, 2023; originally announced October 2023.

Journal ref: A&A 680, A35 (2023)

Gaia Focused Product Release: A catalogue of sources around quasars to search for strongly lensed quasars

Authors: Gaia Collaboration, A. Krone-Martins, C. Ducourant, L. Galluccio, L. Delchambre, I. Oreshina-Slezak, R. Teixeira, J. Braine, J. -F. Le Campion, F. Mignard, W. Roux, A. Blazere, L. Pegoraro, A. G. A. Brown, A. Vallenari, T. Prusti, J. H. J. de Bruijne, F. Arenou, C. Babusiaux, A. Barbier, M. Biermann, O. L. Creevey, D. W. Evans, L. Eyer, R. Guerra , et al. (376 additional authors not shown)

Abstract: Context. Strongly lensed quasars are fundamental sources for cosmology. The Gaia space mission covers the entire sky with the unprecedented resolution of $0.18$" in the optical, making it an ideal instrument to search for gravitational lenses down to the limiting magnitude of 21. Nevertheless, the previous Gaia Data Releases are known to be incomplete for small angular separations such as those ex… ▽ More Context. Strongly lensed quasars are fundamental sources for cosmology. The Gaia space mission covers the entire sky with the unprecedented resolution of $0.18$" in the optical, making it an ideal instrument to search for gravitational lenses down to the limiting magnitude of 21. Nevertheless, the previous Gaia Data Releases are known to be incomplete for small angular separations such as those expected for most lenses. Aims. We present the Data Processing and Analysis Consortium GravLens pipeline, which was built to analyse all Gaia detections around quasars and to cluster them into sources, thus producing a catalogue of secondary sources around each quasar. We analysed the resulting catalogue to produce scores that indicate source configurations that are compatible with strongly lensed quasars. Methods. GravLens uses the DBSCAN unsupervised clustering algorithm to detect sources around quasars. The resulting catalogue of multiplets is then analysed with several methods to identify potential gravitational lenses. We developed and applied an outlier scoring method, a comparison between the average BP and RP spectra of the components, and we also used an extremely randomised tree algorithm. These methods produce scores to identify the most probable configurations and to establish a list of lens candidates. Results. We analysed the environment of 3 760 032 quasars. A total of 4 760 920 sources, including the quasars, were found within 6" of the quasar positions. This list is given in the Gaia archive. In 87\% of cases, the quasar remains a single source, and in 501 385 cases neighbouring sources were detected. We propose a list of 381 lensed candidates, of which we identified 49 as the most promising. Beyond these candidates, the associate tables in this Focused Product Release allow the entire community to explore the unique Gaia data for strong lensing studies further. △ Less

Submitted 10 October, 2023; originally announced October 2023.

Comments: 35 pages, 60 figures, accepted for publication by Astronomy and Astrophysics

Journal ref: A&A 685, A130 (2024)

Gaia Focused Product Release: Radial velocity time series of long-period variables

Authors: Gaia Collaboration, Gaia Collaboration, M. Trabucchi, N. Mowlavi, T. Lebzelter, I. Lecoeur-Taibi, M. Audard, L. Eyer, P. García-Lario, P. Gavras, B. Holl, G. Jevardat de Fombelle, K. Nienartowicz, L. Rimoldini, P. Sartoretti, R. Blomme, Y. Frémat, O. Marchal, Y. Damerdji, A. G. A. Brown, A. Guerrier, P. Panuzzo, D. Katz, G. M. Seabroke, K. Benson , et al. (382 additional authors not shown)

Abstract: The third Gaia Data Release (DR3) provided photometric time series of more than 2 million long-period variable (LPV) candidates. Anticipating the publication of full radial-velocity (RV) in DR4, this Focused Product Release (FPR) provides RV time series for a selection of LPVs with high-quality observations. We describe the production and content of the Gaia catalog of LPV RV time series, and the… ▽ More The third Gaia Data Release (DR3) provided photometric time series of more than 2 million long-period variable (LPV) candidates. Anticipating the publication of full radial-velocity (RV) in DR4, this Focused Product Release (FPR) provides RV time series for a selection of LPVs with high-quality observations. We describe the production and content of the Gaia catalog of LPV RV time series, and the methods used to compute variability parameters published in the Gaia FPR. Starting from the DR3 LPVs catalog, we applied filters to construct a sample of sources with high-quality RV measurements. We modeled their RV and photometric time series to derive their periods and amplitudes, and further refined the sample by requiring compatibility between the RV period and at least one of the $G$, $G_{\rm BP}$, or $G_{\rm RP}$ photometric periods. The catalog includes RV time series and variability parameters for 9\,614 sources in the magnitude range $6\lesssim G/{\rm mag}\lesssim 14$, including a flagged top-quality subsample of 6\,093 stars whose RV periods are fully compatible with the values derived from the $G$, $G_{\rm BP}$, and $G_{\rm RP}$ photometric time series. The RV time series contain a mean of 24 measurements per source taken unevenly over a duration of about three years. We identify the great most sources (88%) as genuine LPVs, with about half of them showing a pulsation period and the other half displaying a long secondary period. The remaining 12% consists of candidate ellipsoidal binaries. Quality checks against RVs available in the literature show excellent agreement. We provide illustrative examples and cautionary remarks. The publication of RV time series for almost 10\,000 LPVs constitutes, by far, the largest such database available to date in the literature. The availability of simultaneous photometric measurements gives a unique added value to the Gaia catalog (abridged) △ Less

Submitted 9 October, 2023; originally announced October 2023.

Comments: 36 pages, 38 figures

The Aligned Orbit of the Eccentric Proto Hot Jupiter TOI-3362b

Authors: Juan I. Espinoza-Retamal, Rafael Brahm, Cristobal Petrovich, Andrés Jordán, Guðmundur Stefánsson, Elyar Sedaghati, Melissa J. Hobson, Diego J. Muñoz, Gavin Boyle, Rodrigo Leiva, Vincent Suc

Abstract: High-eccentricity tidal migration predicts the existence of highly eccentric proto-hot Jupiters on the "tidal circularization track," meaning that they might eventually become hot Jupiters, but that their migratory journey remains incomplete. Having experienced moderate amounts of the tidal reprocessing of their orbital elements, proto-hot Jupiters systems can be powerful test beds for the underly… ▽ More High-eccentricity tidal migration predicts the existence of highly eccentric proto-hot Jupiters on the "tidal circularization track," meaning that they might eventually become hot Jupiters, but that their migratory journey remains incomplete. Having experienced moderate amounts of the tidal reprocessing of their orbital elements, proto-hot Jupiters systems can be powerful test beds for the underlying mechanisms of eccentricity growth. Notably, they may be used for discriminating between variants of high-eccentricity migration, each predicting a distinct evolution of misalignment between the star and the planet's orbit. We constrain the spin-orbit misalignment of the proto-hot Jupiter TOI-3362b with high-precision radial velocity observations using ESPRESSO at VLT. The observations reveal a sky-projected obliquity $λ= 1.2_{-2.7}^{+2.8}$ deg and constrain the orbital eccentricity to $e=0.720 \pm 0.016$, making it one of the most eccentric gas giants for which the obliquity has been measured. The large eccentricity and the striking orbit alignment of the planet suggest that ongoing coplanar high-eccentricity migration driven by a distant companion is a possible explanation for the system's architecture. This distant companion would need to reside beyond 5 au at 95% confidence to be compatible with the available radial velocity observations. △ Less

Submitted 29 November, 2023; v1 submitted 6 September, 2023; originally announced September 2023.

Comments: ApJL 958 L20 (9 pages, 6 figures)

Journal ref: ApJL 958 L20 (2023)

Orbital alignment of the eccentric warm Jupiter TOI-677 b

Authors: Elyar Sedaghati, Andrés Jordán, Rafael Brahm, Diego J. Muñoz, Cristobal Petrovich, Melissa J. Hobson

Abstract: Warm Jupiters lay out an excellent laboratory for testing models of planet formation and migration. Their separation from the host star makes tidal reprocessing of their orbits ineffective, which preserves the orbital architectures that result from the planet-forming process. Among the measurable properties, the orbital inclination with respect to the stellar rotational axis, stands out as a cruci… ▽ More Warm Jupiters lay out an excellent laboratory for testing models of planet formation and migration. Their separation from the host star makes tidal reprocessing of their orbits ineffective, which preserves the orbital architectures that result from the planet-forming process. Among the measurable properties, the orbital inclination with respect to the stellar rotational axis, stands out as a crucial diagnostic for understanding the migration mechanisms behind the origin of close-in planets. Observational limitations have made the procurement of spin-orbit measurements heavily biased toward hot Jupiter systems. In recent years, however, high-precision spectroscopy has begun to provide obliquity measurements for planets well into the warm Jupiter regime. In this study, we present Rossiter-McLaughlin (RM) measurements of the projected obliquity angle for the warm Jupiter TOI-677 b using ESPRESSO at the VLT. TOI-677 b exhibits an extreme degree of alignment ($λ= 0.3 \pm 1.3$ deg), which is particularly puzzling given its significant eccentricity ($e \approx 0.45$). TOI-677 b thus joins a growing class of close-in giants that exhibit large eccentricities and low spin-orbit angles, which is a configuration not predicted by existing models. We also present the detection of a candidate outer brown dwarf companion on an eccentric, wide orbit ($e \approx 0.4$ and $P \approx 13$ yr). Using simple estimates, we show that this companion is unlikely to be the cause of the unusual orbit of TOI-677 b. Therefore, it is essential that future efforts prioritize the acquisition of RM measurements for warm Jupiters. △ Less

Submitted 22 August, 2023; v1 submitted 14 July, 2023; originally announced July 2023.

Comments: 15 pages, 6 figures, 5 tables. Accepted for publication in the Astronomical Journal

Journal ref: The Astronomical Journal, 166:130 (12pp), 2023 September

Three long period transiting giant planets from TESS

Authors: Rafael Brahm, Solène Ulmer-Moll, Melissa J. Hobson, Andrés Jordán, Thomas Henning, Trifon Trifonov, Matías I. Jones, Martin Schlecker, Nestor Espinoza, Felipe I. Rojas, Pascal Torres, Paula Sarkis, Marcelo Tala, Jan Eberhardt, Diana Kossakowski, Diego J. Muñoz, Joel D. Hartman, Gavin Boyle, Vincent Suc, François Bouchy, Adrien Deline, Guillaume Chaverot, Nolan Grieves, Monika Lendl, Olga Suarez , et al. (30 additional authors not shown)

Abstract: We report the discovery and orbital characterization of three new transiting warm giant planets. These systems were initially identified as presenting single transit events in the light curves generated from the full frame images of the Transiting Exoplanet Survey Satellite (TESS). Follow-up radial velocity measurements and additional light curves were used to determine the orbital periods and con… ▽ More We report the discovery and orbital characterization of three new transiting warm giant planets. These systems were initially identified as presenting single transit events in the light curves generated from the full frame images of the Transiting Exoplanet Survey Satellite (TESS). Follow-up radial velocity measurements and additional light curves were used to determine the orbital periods and confirm the planetary nature of the candidates. The planets orbit slightly metal-rich late F- and early G-type stars. We find that TOI 4406b has a mass of $M_P$= 0.30 $\pm$ 0.04 $M_J$ , a radius of $R_P$= 1.00 $\pm$ 0.02 $R_J$ , and a low eccentricity orbit (e=0.15 $\pm$ 0.05) with a period of P= 30.08364 $\pm$ 0.00005 d . TOI 2338b has a mass of $M_P$= 5.98 $\pm$ 0.20 $M_J$ , a radius of $R_P$= 1.00 $\pm$ 0.01 $R_J$ , and a highly eccentric orbit (e= 0.676 $\pm$ 0.002 ) with a period of P= 22.65398 $\pm$ 0.00002 d . Finally, TOI 2589b has a mass of $M_P$= 3.50 $\pm$ 0.10 $M_J$ , a radius of $R_P$= 1.08 $\pm$ 0.03 $R_J$ , and an eccentric orbit (e = 0.522 $\pm$ 0.006 ) with a period of P= 61.6277 $\pm$ 0.0002 d . TOI 4406b and TOI 2338b are enriched in metals compared to their host stars, while the structure of TOI 2589b is consistent with having similar metal enrichment to its host star. △ Less

Submitted 4 April, 2023; originally announced April 2023.

Comments: 24 pages, 16 figures, accepted in AJ

Circumbinary Accretion: From Binary Stars to Massive Binary Black Holes

Authors: Dong Lai, Diego J. Muñoz

Abstract: We review recent works on the dynamics of circumbinary accretion, including time variability, angular momentum transfer between the disk and the binary, and the secular evolution of accreting binaries. These dynamics can impact stellar binary formation/evolution, circumbinary planet formation/migration, and the evolution of (super)massive black-hole binaries. We discuss the dynamics and evolution… ▽ More We review recent works on the dynamics of circumbinary accretion, including time variability, angular momentum transfer between the disk and the binary, and the secular evolution of accreting binaries. These dynamics can impact stellar binary formation/evolution, circumbinary planet formation/migration, and the evolution of (super)massive black-hole binaries. We discuss the dynamics and evolution of inclined/warped circumbinary disks and connect with recent observations of protoplanetary disks. A special kind of circumbinary accretion involves binaries embedded in "big" disks, which may contribute to the mergers of stellar-mass black holes in AGN disks. Highlights include: $\bullet$ Circumbinary accretion is highly variable, being modulated at $P_{\rm b}$ (the binary period) or $\sim 5P_{\rm p}$, depending on the binary eccentricity $e_{\rm b}$ and mass ratio $q_{\rm b}$. $\bullet$ The inner region of the circumbinary disk can develop coherent eccentric structure, which may modulate the accretion and affect the physical processes (e.g. planet migration) taking place in the disk. $\bullet$ Over long timescales, circumbinary accretion steers binaries toward equal masses, and it does not always lead to binary orbital decay, as is commonly assumed. The secular orbital evolution depends on the binary parameters ($e_{\rm b}$ and $q_{\rm b}$), and on the thermodynamic properties of the accreting gas. $\bullet$ A misaligned disk around a low-eccentricity binary tends to evolve toward coplanarity due to viscous dissipation. But when $e_{\rm b}$ is significant, the disk can evolve toward "polar alignment", with the disk plane perpendicular to the binary plane. △ Less

Submitted 24 January, 2023; v1 submitted 31 October, 2022; originally announced November 2022.

Comments: Draft version of 2023 Annual Reviews article, comments are welcome. 39 pages including references, 16 figures

Eccentric Mergers of Intermediate-Mass Black Holes from Evection Resonances in AGN Disks

Authors: Diego J. Muñoz, Nicholas C. Stone, Cristobal Petrovich, Frederic A. Rasio

Abstract: We apply the theory of nonlinear resonance capture to the problem of a black hole binary (BHB) orbiting a supermassive black hole (SMBH) while embedded in the accretion disk of an active galactic nucleus (AGN). If successful, resonance capture can trigger dramatic growth in the BHB eccentricity, with important consequences for the BHB merger timescale, as well as for the gravitational wave (GW) si… ▽ More We apply the theory of nonlinear resonance capture to the problem of a black hole binary (BHB) orbiting a supermassive black hole (SMBH) while embedded in the accretion disk of an active galactic nucleus (AGN). If successful, resonance capture can trigger dramatic growth in the BHB eccentricity, with important consequences for the BHB merger timescale, as well as for the gravitational wave (GW) signature of such an eccentric merger. This resonance capture may occur when the orbital period around the SMBH (the "outer binary") and the apsidal precession of the BHB (the "inner binary") are in a 1:1 commensurability. This effect is analogous to the phenomenon of lunar evection resonance in the early Sun-Earth-Moon system, with the distinction that in the present case, the BHB apsidal precession is due to general relativity, rather than rotationally-induced distortion. In contrast to the case of lunar evection, however, the inner binary undergoes orbital decay driven by GW emission, rather than orbital expansion driven by tidal dissipation. This distinction fundamentally alters the three-body dynamics, forbidding resonance capture, and limiting eccentricity growth. However, if the BHB migrates through of a gaseous AGN disk, the change in the outer binary can counterbalance the suppressing effect of BHB decay, permitting evection resonance capture and the production of eccentric BHB mergers. We compute the likelihood of resonance capture assuming an agnostic distribution of parameters for the three bodies involved and for the properties of the AGN disk. We find that intermediate-mass ratio BHBs (involving an intermediate-mass black hole and a stellar-mass black hole) are the most likely to be captured into evection resonance and thus undergo an eccentric merger. We also compute the GW signature of these mergers, showing that they can enter the LISA band while eccentric. △ Less

Submitted 12 April, 2022; originally announced April 2022.

Comments: 15 pages (before references and appendices), 11 figures, comments welcome

Preferential Accretion and Circumbinary Disk Precession in Eccentric Binary Systems

Authors: Magdalena Siwek, Rainer Weinberger, Diego Munoz, Lars Hernquist

Abstract: We present a suite of high resolution hydrodynamic simulations of binaries immersed in circumbinary accretion disks (CBDs). For the first time, we investigate the preferential accretion rate as a function of both eccentricity $e_{\rm b}$ and mass ratio $q_{\rm b}$ in a densely sampled parameter space, finding that when compared with circular binaries, mass ratios grow more efficiently in binaries… ▽ More We present a suite of high resolution hydrodynamic simulations of binaries immersed in circumbinary accretion disks (CBDs). For the first time, we investigate the preferential accretion rate as a function of both eccentricity $e_{\rm b}$ and mass ratio $q_{\rm b}$ in a densely sampled parameter space, finding that when compared with circular binaries, mass ratios grow more efficiently in binaries on moderately eccentric orbits ($0.0 \lesssim e_{\rm b} \lesssim 0.4$), and high eccentricities ($e_{\rm b} \gtrsim 0.6 $) suppress mass ratio growth. We suggest that this non-monotonic preferential accretion behaviour may produce an observable shift in the mass ratio distributions of stellar binaries and massive black hole binaries. We further find that the response of a CBD can be divided into three regimes, depending on eccentricity and mass ratio: (i) CBDs around circular binaries always precess freely, whereas CBDs around eccentric binaries either (ii) undergo forced precession or (iii) remain locked at an angle with respect to the binary periapsis. Forced precession in eccentric binaries is associated with strong modulation of individual accretion rates on the precession timescale, a potentially observable signature in accreting binaries with short orbital periods. We provide CBD locking angles and precession rates as a function of $e_{\rm b}$ and $q_{\rm b}$ for our simulation suite. △ Less

Submitted 21 November, 2022; v1 submitted 4 March, 2022; originally announced March 2022.

Comments: 15 pages, 14 figures. Published in Monthly Notices of the Royal Astronomical Society

Outward Migration of Super-Jupiters

Authors: Adam M. Dempsey, Diego J. Muñoz, Yoram Lithwick

Abstract: Recent simulations show that giant planets of about one Jupiter mass migrate inward at a rate that differs from the Type II prediction. Here we show that at higher masses, planets migrate outward. Our result differs from previous ones because of our longer simulation times, lower viscosity, and our boundary conditions that allow the disk to reach viscous steady state. We show that, for planets on… ▽ More Recent simulations show that giant planets of about one Jupiter mass migrate inward at a rate that differs from the Type II prediction. Here we show that at higher masses, planets migrate outward. Our result differs from previous ones because of our longer simulation times, lower viscosity, and our boundary conditions that allow the disk to reach viscous steady state. We show that, for planets on circular orbits, the transition from inward to outward migration coincides with the known transition from circular to eccentric disks that occurs for planets more massive than a few Jupiters. In an eccentric disk, the torque on the outer disk weakens due to two effects: the planet launches weaker waves, and those waves travel further before damping. As a result, the torque on the inner disk dominates, and the planet pushes itself outward. Our results suggest that the many super-Jupiters observed by direct-imaging at large distances from the star may have gotten there by outward migration. △ Less

Submitted 3 September, 2021; v1 submitted 11 May, 2021; originally announced May 2021.

Comments: 11 pages, 5 figures. ApJL accepted

Diurnal temperature variation as the source of the preferential direction of fractures on asteroids: theoretical model for the case of Bennu

Authors: D. Uribe, M. Delbo, P. -O. Bouchard, D. Pino Muñoz

Abstract: It has been shown that temperature cycles on airless bodies of our Solar System can cause damaging of surface materials. Nevertheless, propagation mechanisms in the case of space objects are still poorly understood. Present work combines a thermoelasticity model together with linear elastic fracture mechanics theory to predict fracture propagation in the presence of thermal gradients generated by… ▽ More It has been shown that temperature cycles on airless bodies of our Solar System can cause damaging of surface materials. Nevertheless, propagation mechanisms in the case of space objects are still poorly understood. Present work combines a thermoelasticity model together with linear elastic fracture mechanics theory to predict fracture propagation in the presence of thermal gradients generated by diurnal temperature cycling and under conditions similar to those existing on the asteroid Bennu. The crack direction is computed using the maximal strain energy release rate criterion, which is implemented using finite elements and the so-called G$θ$ method (Uribe-Suárez et al. 2020. Eng. Fracture Mech. 227:106918). Using the implemented methodology, crack propagation direction for an initial crack tip in different positions and for different orientations is computed. It is found that cracks preferentially propagate in the North to South (N-S), in the North-East to South-West (NE-SW) and in the North-West to South-East (NW-SE) directions. Finally, thermal fatigue analysis was performed in order to estimate the crack growth rate. Computed value is in good agreement with available experimental evidence. △ Less

Submitted 3 February, 2021; originally announced February 2021.

Comments: To be published in Icarus-Journal

Gaia Early Data Release 3: The Galactic anticentre

Authors: Gaia Collaboration, T. Antoja, P. McMillan, G. Kordopatis, P. Ramos, A. Helmi, E. Balbinot, T. Cantat-Gaudin, L. Chemin, F. Figueras, C. Jordi, S. Khanna, M. Romero-Gomez, G. Seabroke, A. G. A. Brown, A. Vallenari, T. Prusti, J. H. J. de Bruijne, C. Babusiaux, M. Biermann, O. L. Creevey, D. W. Evans, L. Eyer, A. Hutton, F. Jansen , et al. (395 additional authors not shown)

Abstract: We aim to demonstrate the scientific potential of the Gaia Early Data Release 3 (EDR3) for the study of the Milky Way structure and evolution. We used astrometric positions, proper motions, parallaxes, and photometry from EDR3 to select different populations and components and to calculate the distances and velocities in the direction of the anticentre. We explore the disturbances of the current d… ▽ More We aim to demonstrate the scientific potential of the Gaia Early Data Release 3 (EDR3) for the study of the Milky Way structure and evolution. We used astrometric positions, proper motions, parallaxes, and photometry from EDR3 to select different populations and components and to calculate the distances and velocities in the direction of the anticentre. We explore the disturbances of the current disc, the spatial and kinematical distributions of early accreted versus in-situ stars, the structures in the outer parts of the disc, and the orbits of open clusters Berkeley 29 and Saurer 1. We find that: i) the dynamics of the Galactic disc are very complex with vertical asymmetries, and new correlations, including a bimodality with disc stars with large angular momentum moving vertically upwards from below the plane, and disc stars with slightly lower angular momentum moving preferentially downwards; ii) we resolve the kinematic substructure (diagonal ridges) in the outer parts of the disc for the first time; iii) the red sequence that has been associated with the proto-Galactic disc that was present at the time of the merger with Gaia-Enceladus-Sausage is currently radially concentrated up to around 14 kpc, while the blue sequence that has been associated with debris of the satellite extends beyond that; iv) there are density structures in the outer disc, both above and below the plane, most probably related to Monoceros, the Anticentre Stream, and TriAnd, for which the Gaia data allow an exhaustive selection of candidate member stars and dynamical study; and v) the open clusters Berkeley~29 and Saurer~1, despite being located at large distances from the Galactic centre, are on nearly circular disc-like orbits. We demonstrate how, once again, the Gaia are crucial for our understanding of the different pieces of our Galaxy and their connection to its global structure and history. △ Less

Submitted 26 April, 2021; v1 submitted 14 January, 2021; originally announced January 2021.

Comments: Gaia EDR3 performance verification paper, version 2 closer to published version in A&A, complete list of authors

Journal ref: A&A 649, A8 (2021)

Gaia Early Data Release 3: The Gaia Catalogue of Nearby Stars

Authors: Gaia Collaboration, R. L. Smart, L. M. Sarro, J. Rybizki, C. Reylé, A. C. Robin, N. C. Hambly, U. Abbas, M. A. Barstow, J. H. J. de Bruijne, B. Bucciarelli, J. M. Carrasco, W. J. Cooper, S. T. Hodgkin, E. Masana, D. Michalik, J. Sahlmann, A. Sozzetti, A. G. A. Brown, A. Vallenari, T. Prusti, C. Babusiaux, M. Biermann, O. L. Creevey, D. W. Evans , et al. (398 additional authors not shown)

Abstract: We produce a clean and well-characterised catalogue of objects within 100\,pc of the Sun from the \G\ Early Data Release 3. We characterise the catalogue through comparisons to the full data release, external catalogues, and simulations. We carry out a first analysis of the science that is possible with this sample to demonstrate its potential and best practices for its use. The selection of obj… ▽ More We produce a clean and well-characterised catalogue of objects within 100\,pc of the Sun from the \G\ Early Data Release 3. We characterise the catalogue through comparisons to the full data release, external catalogues, and simulations. We carry out a first analysis of the science that is possible with this sample to demonstrate its potential and best practices for its use. The selection of objects within 100\,pc from the full catalogue used selected training sets, machine-learning procedures, astrometric quantities, and solution quality indicators to determine a probability that the astrometric solution is reliable. The training set construction exploited the astrometric data, quality flags, and external photometry. For all candidates we calculated distance posterior probability densities using Bayesian procedures and mock catalogues to define priors. Any object with reliable astrometry and a non-zero probability of being within 100\,pc is included in the catalogue. We have produced a catalogue of \NFINAL\ objects that we estimate contains at least 92\% of stars of stellar type M9 within 100\,pc of the Sun. We estimate that 9\% of the stars in this catalogue probably lie outside 100\,pc, but when the distance probability function is used, a correct treatment of this contamination is possible. We produced luminosity functions with a high signal-to-noise ratio for the main-sequence stars, giants, and white dwarfs. We examined in detail the Hyades cluster, the white dwarf population, and wide-binary systems and produced candidate lists for all three samples. We detected local manifestations of several streams, superclusters, and halo objects, in which we identified 12 members of \G\ Enceladus. We present the first direct parallaxes of five objects in multiple systems within 10\,pc of the Sun. △ Less

Submitted 3 December, 2020; originally announced December 2020.

Comments: 45 Pages, 39 figures in main part and 18 in appendix, tables on CDS

Journal ref: A&A 649, A6 (2021)

Gaia Early Data Release 3: Acceleration of the solar system from Gaia astrometry

Authors: Gaia Collaboration, S. A. Klioner, F. Mignard, L. Lindegren, U. Bastian, P. J. McMillan, J. Hernández, D. Hobbs, M. Ramos-Lerate, M. Biermann, A. Bombrun, A. de Torres, E. Gerlach, R. Geyer, T. Hilger, U. Lammers, H. Steidelmüller, C. A. Stephenson, A. G. A. Brown, A. Vallenari, T. Prusti, J. H. J. de Bruijne, C. Babusiaux, O. L. Creevey, D. W. Evans , et al. (392 additional authors not shown)

Abstract: Context. Gaia Early Data Release 3 (Gaia EDR3) provides accurate astrometry for about 1.6 million compact (QSO-like) extragalactic sources, 1.2 million of which have the best-quality five-parameter astrometric solutions. Aims. The proper motions of QSO-like sources are used to reveal a systematic pattern due to the acceleration of the solar system barycentre with respect to the rest frame of the… ▽ More Context. Gaia Early Data Release 3 (Gaia EDR3) provides accurate astrometry for about 1.6 million compact (QSO-like) extragalactic sources, 1.2 million of which have the best-quality five-parameter astrometric solutions. Aims. The proper motions of QSO-like sources are used to reveal a systematic pattern due to the acceleration of the solar system barycentre with respect to the rest frame of the Universe. Apart from being an important scientific result by itself, the acceleration measured in this way is a good quality indicator of the Gaia astrometric solution. Methods. The effect of the acceleration is obtained as a part of the general expansion of the vector field of proper motions in Vector Spherical Harmonics (VSH). Various versions of the VSH fit and various subsets of the sources are tried and compared to get the most consistent result and a realistic estimate of its uncertainty. Additional tests with the Gaia astrometric solution are used to get a better idea on possible systematic errors in the estimate. Results. Our best estimate of the acceleration based on Gaia EDR3 is $(2.32 \pm 0.16) \times 10^{-10}$ m s${}^{-2}$ (or $7.33 \pm 0.51$ km s$^{-1}$ Myr${}^{-1}$) towards $α= 269.1^\circ \pm 5.4^\circ$, $δ= -31.6^\circ \pm 4.1^\circ$, corresponding to a proper motion amplitude of $5.05 \pm 0.35$ $μ$as yr${}^{-1}$. This is in good agreement with the acceleration expected from current models of the Galactic gravitational potential. We expect that future Gaia data releases will provide estimates of the acceleration with uncertainties substantially below 0.1 $μ$as yr${}^{-1}$. △ Less

Submitted 3 December, 2020; originally announced December 2020.

Comments: A&A, accepted

Journal ref: A&A 649, A9 (2021)

Gaia Early Data Release 3: Structure and properties of the Magellanic Clouds

Authors: Gaia Collaboration, X. Luri, L. Chemin, G. Clementini, H. E. Delgado, P. J. McMillan, M. Romero-Gómez, E. Balbinot, A. Castro-Ginard, R. Mor, V. Ripepi, L. M. Sarro, M. -R. L. Cioni, C. Fabricius, A. Garofalo, A. Helmi, T. Muraveva, A. G. A. Brown, A. Vallenari, T. Prusti, J. H. J. de, C. Babusiaux, M. Biermann, O. L. Creevey, D. W. Evans , et al. (395 additional authors not shown)

Abstract: We compare the Gaia DR2 and Gaia EDR3 performances in the study of the Magellanic Clouds and show the clear improvements in precision and accuracy in the new release. We also show that the systematics still present in the data make the determination of the 3D geometry of the LMC a difficult endeavour; this is at the very limit of the usefulness of the Gaia EDR3 astrometry, but it may become feasib… ▽ More We compare the Gaia DR2 and Gaia EDR3 performances in the study of the Magellanic Clouds and show the clear improvements in precision and accuracy in the new release. We also show that the systematics still present in the data make the determination of the 3D geometry of the LMC a difficult endeavour; this is at the very limit of the usefulness of the Gaia EDR3 astrometry, but it may become feasible with the use of additional external data. We derive radial and tangential velocity maps and global profiles for the LMC for the several subsamples we defined. To our knowledge, this is the first time that the two planar components of the ordered and random motions are derived for multiple stellar evolutionary phases in a galactic disc outside the Milky Way, showing the differences between younger and older phases. We also analyse the spatial structure and motions in the central region, the bar, and the disc, providing new insights into features and kinematics. Finally, we show that the Gaia EDR3 data allows clearly resolving the Magellanic Bridge, and we trace the density and velocity flow of the stars from the SMC towards the LMC not only globally, but also separately for young and evolved populations. This allows us to confirm an evolved population in the Bridge that is slightly shift from the younger population. Additionally, we were able to study the outskirts of both Magellanic Clouds, in which we detected some well-known features and indications of new ones. △ Less

Submitted 4 January, 2021; v1 submitted 3 December, 2020; originally announced December 2020.

Comments: This paper is part of the "demonstration papers" released with Gaia EDR3: https://www.cosmos.esa.int/web/gaia/earlydr3

Journal ref: A&A 649, A7 (2021)

Gaia Early Data Release 3: Summary of the contents and survey properties

Authors: Gaia Collaboration, A. G. A Brown, A. Vallenari, T. Prusti, J. H. J. de Bruijne, C. Babusiaux, M. Biermann, O. L. Creevey, D. W. Evans, L. Eyer, A. Hutton, F. Jansen, C. Jordi, S. A. Klioner, U. Lammers, L. Lindegren, X. Luri, F. Mignard, C. Panem, D. Pourbaix, S. Randich, P. Sartoretti, C. Soubiran, N. A. Walton, F. Arenou , et al. (401 additional authors not shown)

Abstract: We present the early installment of the third Gaia data release, Gaia EDR3, consisting of astrometry and photometry for 1.8 billion sources brighter than magnitude 21, complemented with the list of radial velocities from Gaia DR2. Gaia EDR3 contains celestial positions and the apparent brightness in G for approximately 1.8 billion sources. For 1.5 billion of those sources, parallaxes, proper motio… ▽ More We present the early installment of the third Gaia data release, Gaia EDR3, consisting of astrometry and photometry for 1.8 billion sources brighter than magnitude 21, complemented with the list of radial velocities from Gaia DR2. Gaia EDR3 contains celestial positions and the apparent brightness in G for approximately 1.8 billion sources. For 1.5 billion of those sources, parallaxes, proper motions, and the (G_BP-G_RP) colour are also available. The passbands for G, G_BP, and G_RP are provided as part of the release. For ease of use, the 7 million radial velocities from Gaia DR2 are included in this release, after the removal of a small number of spurious values. New radial velocities will appear as part of Gaia DR3. Finally, Gaia EDR3 represents an updated materialisation of the celestial reference frame (CRF) in the optical, the Gaia-CRF3, which is based solely on extragalactic sources. The creation of the source list for Gaia EDR3 includes enhancements that make it more robust with respect to high proper motion stars, and the disturbing effects of spurious and partially resolved sources. The source list is largely the same as that for Gaia DR2, but it does feature new sources and there are some notable changes. The source list will not change for Gaia DR3. Gaia EDR3 represents a significant advance over Gaia DR2, with parallax precisions increased by 30 percent, proper motion precisions increased by a factor of 2, and the systematic errors in the astrometry suppressed by 30--40 percent for the parallaxes and by a factor ~2.5 for the proper motions. The photometry also features increased precision, but above all much better homogeneity across colour, magnitude, and celestial position. A single passband for G, G_BP, and G_RP is valid over the entire magnitude and colour range, with no systematics above the 1 percent level. △ Less

Submitted 9 June, 2021; v1 submitted 2 December, 2020; originally announced December 2020.

Comments: Accepted for A&A Special Issue on Gaia EDR3, 21 pages, 2 figures. This version includes the updates in the erratum (https://doi.org/10.1051/0004-6361/202039657e)

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

Kozai Migration Naturally Explains the White Dwarf Planet WD1856b

Authors: Diego J. Muñoz, Cristobal Petrovich

Abstract: The Jovian-sized object WD~1856~b transits a white dwarf (WD) in a compact $1.4$-day orbit. Unlikely to have endured stellar evolution in its current orbit, WD~1856~b is thought to have migrated from much wider separations. Because the WD is old, and a member of a well-characterized hierarchical multiple, the well-known Kozai mechanism provides an effective migration channel for WD~1856~b. Moreove… ▽ More The Jovian-sized object WD~1856~b transits a white dwarf (WD) in a compact $1.4$-day orbit. Unlikely to have endured stellar evolution in its current orbit, WD~1856~b is thought to have migrated from much wider separations. Because the WD is old, and a member of a well-characterized hierarchical multiple, the well-known Kozai mechanism provides an effective migration channel for WD~1856~b. Moreover, the lack of tides in the star allows us to directly connect the current semi-major axis to the pre-migration one, from which we can infer the initial conditions of the system. By further demanding that successful migrators survive all previous phases of stellar evolution, we are able to constrain the mass of WD~1856~b to be $\simeq0.7-3M_{\rm J}$ and its main sequence semi-major axis to be $\simeq 2-2.5$ au. These properties imply that WD~1856~b was born a typical gas giant. We further estimate the occurrence rate of Kozai-migrated planets around WDs to be ${\cal O}(10^{-3}{-}10^{-4})$, suggesting that WD~1856~b is the only one in the {\it TESS} sample, but implying ${\cal O}(10^2)$ future detections by LSST. In a sense, WD~1856~b was an ordinary Jovian planet that underwent an extraordinary dynamical history. △ Less

Submitted 9 October, 2020; originally announced October 2020.

Comments: 9 pages, 5 figuress (submitted)

A disk-driven resonance as the origin of high inclinations of close-in planets

Authors: Cristobal Petrovich, Diego J. Muñoz, Kaitlin M. Kratter, Renu Malhotra

Abstract: The recent characterization of transiting close-in planets has revealed an intriguing population of sub-Neptunes with highly tilted and even polar orbits relative to their host star's equator. Any viable theory for the origin of these close-in, polar planets must explain (1) the observed stellar obliquities, (2) the substantial eccentricities, and (3) the existence of Jovian companions with large… ▽ More The recent characterization of transiting close-in planets has revealed an intriguing population of sub-Neptunes with highly tilted and even polar orbits relative to their host star's equator. Any viable theory for the origin of these close-in, polar planets must explain (1) the observed stellar obliquities, (2) the substantial eccentricities, and (3) the existence of Jovian companions with large mutual inclinations. In this work, we propose a theoretical model that satisfies these requirements without invoking tidal dissipation or large primordial inclinations. Instead, tilting is facilitated by the protoplanetary disk dispersal during the late stage of planet formation, initiating a process of resonance sweeping and parametric instability. This mechanism consists of two steps. First, a nodal secular resonance excites the inclination to large values; then, once the inclination reaches a critical value, a linear eccentric instability is triggered, which detunes the resonance and ends inclination growth. The critical inclination is pushed to high values by general relativistic precession, making polar orbits an inherently post-Newtonian outcome. Our model predicts that polar, close-in sub-Neptunes coexist with cold Jupiters in low stellar obliquity orbits. △ Less

Submitted 19 August, 2020; originally announced August 2020.

Comments: submitted to AAS Journals

Long-Lived Eccentric Modes in Circumbinary Disks

Authors: Diego J. Muñoz, Yoram Lithwick

Abstract: Hydrodynamical simulations show that circumbinary disks become eccentric, even when the binary is circular. Here we demonstrate that, in steady state, the disk's eccentricity behaves as a long-lived free mode trapped by turning points that naturally arise from a continuously truncated density profile. Consequently, both the disk's precession rate and eccentricity profile may be calculated via the… ▽ More Hydrodynamical simulations show that circumbinary disks become eccentric, even when the binary is circular. Here we demonstrate that, in steady state, the disk's eccentricity behaves as a long-lived free mode trapped by turning points that naturally arise from a continuously truncated density profile. Consequently, both the disk's precession rate and eccentricity profile may be calculated via the simple linear theory for perturbed pressure-supported disks. By formulating and solving the linear theory we find that (i) surprisingly, the precession rate is roughly determined by the binary's quadrupole, even when the quadrupole is very weak relative to pressure; (ii) the eccentricity profile is largest near the inner edge of the disk, and falls exponentially outwards; and (iii) the results from linear theory indeed agree with what is found in simulations. Understanding the development of eccentric modes in circumbinary disks is a crucial first step for understanding the long term (secular) exchange of eccentricity, angular momentum and mass between the binary and the gas. Potential applications include the search for a characteristic kinematic signature in disks around candidate binaries and precession-induced modulation of accretion over long timescales. △ Less

Submitted 8 November, 2020; v1 submitted 18 August, 2020; originally announced August 2020.

Comments: 11 pages, 11 figures (accepted to ApJ)

Inner Boundary Condition in Quasi-Lagrangian Simulations of Accretion Disks

Authors: Adam M. Dempsey, Diego Muñoz, Yoram Lithwick

Abstract: In simulations of viscously evolving accretion disks, the inner boundary condition is particularly important. If treated incorrectly, it induces incorrect behavior very quickly, because the viscous time is shortest near the inner boundary. Recent work has determined the correct inner boundary in Eulerian simulations. But in quasi-Lagrangian simulations (e.g., SPH, moving mesh, and mesh-less), wher… ▽ More In simulations of viscously evolving accretion disks, the inner boundary condition is particularly important. If treated incorrectly, it induces incorrect behavior very quickly, because the viscous time is shortest near the inner boundary. Recent work has determined the correct inner boundary in Eulerian simulations. But in quasi-Lagrangian simulations (e.g., SPH, moving mesh, and mesh-less), where the inner boundary is modeled by removing mass within a finite zone, the inner density profile typically becomes anomalously depleted. Here we show how the boundary condition should be applied in such codes, via a simple modification of the usual approach: when one removes mass, one must speed up the remaining material so that the disk's angular momentum is unchanged. We show with both 1D and 2D moving-mesh (AREPO) simulations that this scheme works as desired in viscously evolving disks. It produces no spurious density depletions and is independent of the mass removal rate, provided that the disk is adequately resolved and that the mass removal rate is not so extreme as to trigger instabilities. This "torque-free" mass removal technique permits the use of quasi-Lagrangian codes to simulate viscously evolving disks, while including a variety of additional effects. As an example, we apply our scheme to a 2D simulation of an accretion disk perturbed by a very massive planet, in which the disk is evolved to viscous steady state. △ Less

Submitted 17 March, 2020; v1 submitted 12 February, 2020; originally announced February 2020.

Comments: 10 pages, 5 figures, Accepted to ApJL

Full-field modeling of heat transfer in asteroid regolith: Radiative thermal conductivity of polydisperse particulates

Authors: Andrew J. Ryan, Daniel Pino Muñoz, Marc Bernacki, Marco Delbo

Abstract: Characterizing the surface material of an asteroid is important for understanding its geology and for informing mission decisions, such as the selection of a sample site. Diurnal surface temperature amplitudes are directly related to the thermal properties of the materials on the surface. We describe a numerical model for studying the thermal conductivity of particulate regolith in vacuum. Heat di… ▽ More Characterizing the surface material of an asteroid is important for understanding its geology and for informing mission decisions, such as the selection of a sample site. Diurnal surface temperature amplitudes are directly related to the thermal properties of the materials on the surface. We describe a numerical model for studying the thermal conductivity of particulate regolith in vacuum. Heat diffusion and surface-to-surface radiation calculations are performed using the finite element (FE) method in three-dimensional meshed geometries of randomly packed spherical particles. We validate the model for test cases where the total solid and radiative conductivity values of particulates with monodisperse particle size frequency distributions (SFDs) are determined at steady-state thermal conditions. Then, we use the model to study the bulk radiative thermal conductivity of particulates with polydisperse, cumulative power-law particle SFDs. We show that for each polydisperse particulate geometry tested, there is a corresponding monodisperse geometry with some effective particle diameter that has an identical radiative thermal conductivity. These effective diameters are found to correspond very well to the Sauter mean particle diameter, which is essentially the surface-area-weighted mean. Next, we show that the thermal conductivity of the particle material can have an important effect on the radiative component of the thermal conductivity of particulates, especially if the particle material conductivity is very low or the spheres are relatively large, owing to non-isothermality in each particle. We provide an empirical correlation to predict the effects of non-isothermality on radiative thermal conductivity in both monodisperse and polydisperse particulates. △ Less

Submitted 23 February, 2020; v1 submitted 1 February, 2020; originally announced February 2020.

Comments: Accepted by Journal of Geophysical Research: Planets

Circumbinary Accretion from Finite and Infinite Disks

Authors: Diego Muñoz, Dong Lai, Kaitlin Kratter, Ryan Miranda

Abstract: We carry out 2D viscous hydrodynamics simulations of circumbinary disk (CBD) accretion using {\footnotesize AREPO}. We resolve the accretion flow from a large-scale CBD down to the streamers and disks around individual binary components. Extending our recent studies \citep{mun19}, we consider circular binaries with various mass ratios ($0.1\leq q_{\rm{b}}\leq1$) and study accretion from ``infinite… ▽ More We carry out 2D viscous hydrodynamics simulations of circumbinary disk (CBD) accretion using {\footnotesize AREPO}. We resolve the accretion flow from a large-scale CBD down to the streamers and disks around individual binary components. Extending our recent studies \citep{mun19}, we consider circular binaries with various mass ratios ($0.1\leq q_{\rm{b}}\leq1$) and study accretion from ``infinite'', steady-supply disks and from finite-sized, viscously spreading tori. For ``infinite'' disks, a global steady state can be reached, and the accretion variability has a dominant frequency ${\sim}0.2Ω_{\rm{b}}$ for $q_{\rm{b}}>0.5$ and $Ω_{\rm{b}}$ for $q_{\rm{b}}<0.5$, ($Ω_{\rm{b}}$ is the binary angular frequency). We find that the accretion ``eigenvalue'' $l_0$ -- the net angular momentum transfer from the disk to the binary per unit accreted mass -- is always positive and falls in the range ($0.65$-$0.85)a_{\rm b}^2Ω_{\rm{b}}$ (with $a_{\rm{b}}$ the binary separation), depending weakly on the mass ratio and viscosity. This leads to binary expansion when $q_{\rm{b}}\gtrsim0.3$. Accretion from a finite torus can be separated into two phases: an initial transient phase, corresponding to the filling of the binary cavity, followed by a viscous pseudo-stationary phase, during which the torus viscously spreads and accretes onto the binary. In the viscous phase, the net torque on the binary per unit accreted mass is close to $l_0$, the value derived for ``infinite'' disks. We conclude that similar-mass binaries accreting from CBDs gain angular momentum and expand over long time scales. This result significantly impacts the coalescence of supermassive binary black holes and newly formed binary stars. We offer a word of caution against conclusions drawn from simulations of transient accretion onto empty circumbinary cavities. △ Less

Submitted 16 December, 2019; v1 submitted 10 October, 2019; originally announced October 2019.

Comments: 10 pages, 13 figures, accepted to ApJ

Hydrodynamics of circumbinary accretion: Angular momentum transfer and binary orbital evolution

Authors: Diego J. Muñoz, Ryan Miranda, Dong Lai

Abstract: We carry out 2D viscous hydrodynamical simulations of circumbinary accretion using the AREPO code. We self-consistently compute the accretion flow over a wide range of spatial scales, from the circumbinary disk (CBD) far from the central binary, through accretion streamers, to the disks around individual binary components, resolving the flow down to 2% of the binary separation. We focus on equal-m… ▽ More We carry out 2D viscous hydrodynamical simulations of circumbinary accretion using the AREPO code. We self-consistently compute the accretion flow over a wide range of spatial scales, from the circumbinary disk (CBD) far from the central binary, through accretion streamers, to the disks around individual binary components, resolving the flow down to 2% of the binary separation. We focus on equal-mass binaries with arbitrary eccentricities. We evolve the flow over long (viscous) timescales until a quasi-steady is reached, in which the mass supply rate at large distances $\dot{M}_0$ (assumed constant) equals the time-averaged mass transfer rate across the disk and the total mass accretion rate onto the binary components. This quasi-steady state allows us to compute the secular angular momentum transfer rate onto the binary, $\langle\dot{J}_b\rangle$, and the resulting orbital evolution. Through direct computation of the gravitational and accretion torques on the binary, we find that $\langle\dot{J}_b\rangle$ is consistently positive (i.e., the binary gains angular momentum), with $l_0\equiv\langle\dot{J}_b\rangle/\dot M_0$ in the range of $(0.4-0.8)a_b^2Ω_b$, depending on the binary eccentricity (where $a_b,~Ω_b$ are the binary semi-major axis and angular frequency); we also find that this $\langle\dot{J}_b\rangle$ is equal to the net angular momentum current across the CBD, indicating that global angular momentum balance is achieved in our simulations. We compute the time-averaged rate of change of the binary orbital energy for eccentric binaries, and thus obtain the secular rates $\langle\dot a_b\rangle$ and $\langle \dot{e}_b\rangle$. In all cases, $\langle\dot{a}_b\rangle$ is positive, i.e., the binary expands while accreting. We discuss the implications of our results for the merger of supermassive binary black holes and for the formation of close stellar binaries. △ Less

Submitted 9 January, 2019; v1 submitted 10 October, 2018; originally announced October 2018.

Comments: 18 pages, 13 figures, 1 table. Accepted for publication in ApJ

Small-scale magnetic flux emergence in the quiet Sun

Authors: F. Moreno-Insertis, J. Martinez-Sykora, V. H. Hansteen, D. Muñoz

Abstract: Small bipolar magnetic features are observed to appear in the interior of individual granules in the quiet Sun, signaling the emergence of tiny magnetic loops from the solar interior. We study the origin of those features as part of the magnetoconvection process in the top layers of the convection zone. Two quiet-Sun magnetoconvection models, calculated with the radiation-magnetohydrodynamic (MHD)… ▽ More Small bipolar magnetic features are observed to appear in the interior of individual granules in the quiet Sun, signaling the emergence of tiny magnetic loops from the solar interior. We study the origin of those features as part of the magnetoconvection process in the top layers of the convection zone. Two quiet-Sun magnetoconvection models, calculated with the radiation-magnetohydrodynamic (MHD) Bifrost code and with domain stretching from the top layers of the convection zone to the corona, are analyzed. Using 3D visualization as well as a posteriori spectral synthesis of Stokes parameters, we detect the repeated emergence of small magnetic elements in the interior of granules, as in the observations. Additionally, we identify the formation of organized horizontal magnetic sheets covering whole granules. Our approach is twofold, calculating statistical properties of the system, like joint probability density functions (JPDFs), and pursuing individual events via visualization tools. We conclude that the small magnetic loops surfacing within individual granules in the observations may originate from sites at or near the downflows in the granular and mesogranular levels, probably in the first 1 or 1.5 Mm below the surface. We also document the creation of granule-covering magnetic sheet-like structures through the sideways expansion of a small subphotospheric magnetic concentration picked up, and pulled out of the interior, by a nascent granule. The sheet-like structures we found in the models may match the recent observations of Centeno et al. (2017). △ Less

Submitted 1 June, 2018; originally announced June 2018.

Comments: 9 pages, 5 figures, Published in The Astrophysical Journal Letters

Journal ref: The Astrophysical Journal Letters, 859:L26 (7pp), 2018 June 1

Statistical Trends in the Obliquity Distribution of Exoplanet Systems

Authors: Diego J. Muñoz, Hagai B. Perets

Abstract: Important clues on the formation and evolution of planetary systems can be inferred from the stellar obliquity $ψ$. We study the distribution of obliquities using the California-Kepler Survey and the TEPCat Catalog of Rossiter-MacLaughlin (RM) measurements, from which we extract, respectively, 275 and 118 targets. We infer a "best fit" obliquity distribution in $ψ$ with a single parameter $κ$. Lar… ▽ More Important clues on the formation and evolution of planetary systems can be inferred from the stellar obliquity $ψ$. We study the distribution of obliquities using the California-Kepler Survey and the TEPCat Catalog of Rossiter-MacLaughlin (RM) measurements, from which we extract, respectively, 275 and 118 targets. We infer a "best fit" obliquity distribution in $ψ$ with a single parameter $κ$. Large values of $κ$ imply that $ψ$ is distributed narrowly around zero, while small values imply approximate isotropy. Our findings are: (1) The distribution of $ψ$ in Kepler systems is narrower than found by previous studies and consistent with $κ\sim15$ (mean $\langleψ\rangle\sim19^\circ$ and spread $ψ\sim10^\circ$). (2) The value of $κ$ in Kepler systems does not depend, at a statistically significant level, on planet multiplicity, stellar multiplicity or stellar age; on the other hand, metal-rich hosts, small planet hosts and long-period planet hosts tend to be more oblique than the general sample (at a $\sim$2.5-$σ$ significance level). (3) Hot Jupiter (HJ) systems with RM measurements are consistent with $κ\sim2$, more broadly distributed than the general Kepler population. (4) A separation of the RM sample into cooler ($T_{\rm eff}\lesssim$6250 K) and hotter ($T_{\rm eff}\gtrsim$6250 K) HJ hosts results in two distinct distributions, $κ_{\rm cooler}\sim4$ and $κ_{\rm hotter}\sim1$ (4-$σ$ significance), both more oblique than the Kepler sample. We hypothesize that the total mass in planets may be behind the increasing obliquity with metallicity and planet radius, and that the period dependence could be due to primordial disk alignment rather than realignment of stellar spin. △ Less

Submitted 19 October, 2018; v1 submitted 9 May, 2018; originally announced May 2018.

Comments: 16 pages, 11 figures, accepted to AJ

An outflow in the Seyfert ESO 362-G18 revealed by Gemini-GMOS/IFU Observations

Authors: Pedro K. Humire, Neil M. Nagar, Carolina Finlez, Verónica Firpo, Roy Slater, Davide Lena, Pamela R. Soto, Dania Muñoz, Rogemar A. Riffel, H. R. Schmitt, S. B. Kraemer, Allan Schnorr-Müller, T. C. Fischer, Andrew Robinson, Thaisa Storchi-Bergmann, Mike Crenshaw, Martin S Elvis

Abstract: We present two-dimensional stellar and gaseous kinematics of the inner 0.7 $\times$ 1.2 kpc$^{2}$ of the Seyfert galaxy ESO 362-G18, derived from optical spectra obtained with the GMOS/IFU on the Gemini South telescope at a spatial resolution of $\approx$170 pc and spectral resolution of 36 km s$^{-1}$. ESO 362-G18 is a strongly perturbed galaxy of morphological type Sa or S0/a, with a minor merge… ▽ More We present two-dimensional stellar and gaseous kinematics of the inner 0.7 $\times$ 1.2 kpc$^{2}$ of the Seyfert galaxy ESO 362-G18, derived from optical spectra obtained with the GMOS/IFU on the Gemini South telescope at a spatial resolution of $\approx$170 pc and spectral resolution of 36 km s$^{-1}$. ESO 362-G18 is a strongly perturbed galaxy of morphological type Sa or S0/a, with a minor merger approaching along the NE direction. Previous studies have shown that the [OIII] emission shows a fan-shaped extension of $\approx$ 10\arcsec\ to the SE. We detect the [OIII] doublet, [NII] and H$α$ emission lines throughout our field of view. The stellar kinematics is dominated by circular motions in the galaxy plane, with a kinematic position angle of $\approx$137$^{\circ}$. The gas kinematics is also dominated by rotation, with kinematic position angles ranging from 122$^{\circ}$ to 139$^{\circ}$. A double-Gaussian fit to the [OIII]$λ$5007 and H$α$ lines, which have the highest signal to noise ratios of the emission lines, reveal two kinematic components: (1) a component at lower radial velocities which we interpret as gas rotating in the galactic disk; and (2) a component with line of sight velocities 100-250 km s$^{-1}$ higher than the systemic velocity, interpreted as originating in the outflowing gas within the AGN ionization cone. We estimate a mass outflow rate of 7.4 $\times$ 10$^{-2}$ M$_{\odot}$ yr$^{-1}$ in the SE ionization cone (this rate doubles if we assume a biconical configuration), and a mass accretion rate on the supermassive black hole (SMBH) of 2.2 $\times$ 10$^{-2}$ M$_{\odot}$ yr$^{-1}$. The total ionized gas mass within $\sim$84 pc of the nucleus is 3.3 $\times$ 10$^{5}$ M$_{\odot}$; infall velocities of $\sim$34 km s$^{-1}$ in this gas would be required to feed both the outflow and SMBH accretion. △ Less

Submitted 9 March, 2018; v1 submitted 21 February, 2018; originally announced February 2018.

Comments: 18 pages, 14 figures

Journal ref: A&A 614, A94 (2018)

Viscous Hydrodynamics Simulations of Circumbinary Accretion Discs: Variability, Quasi-Steady State, and Angular Momentum Transfer

Authors: Ryan Miranda, Diego Muñoz, Dong Lai

Abstract: We carry out numerical simulations of circumbinary discs, solving the viscous hydrodynamics equations on a polar grid covering an extended disc outside the binary co-orbital region. We use carefully controlled outer boundary conditions and long-term integrations to ensure that the disc reaches a quasi-steady state, in which the time-averaged mass accretion rate onto the binary,… ▽ More We carry out numerical simulations of circumbinary discs, solving the viscous hydrodynamics equations on a polar grid covering an extended disc outside the binary co-orbital region. We use carefully controlled outer boundary conditions and long-term integrations to ensure that the disc reaches a quasi-steady state, in which the time-averaged mass accretion rate onto the binary, $\langle\dot{M}\rangle$, matches the mass supply rate at the outer disc. We focus on binaries with comparable masses and a wide range of eccentricities ($e_\mathrm{B}$). For $e_\mathrm{B} \lesssim 0.05$, the mass accretion rate of the binary is modulated at about $5$ times the binary period; otherwise it is modulated at the binary period. The inner part of the circumbinary disc ($r \lesssim 6 a_\mathrm{B}$) generally becomes coherently eccentric. For low and high $e_\mathrm{B}$, the disc line of apsides precesses around the binary, but for intermediate $e_\mathrm{B}$ ($0.2 - 0.4$), it instead becomes locked with that of the binary. By considering the balance of angular momentum transport through the disc by advection, viscous stress, and gravitational torque, we determine the time-averaged net angular momentum transfer rate to the binary, $\langle\dot{J}\rangle$. The specific angular momentum, $l_0 = \langle\dot{J}\rangle/\langle\dot{M}\rangle$, depends non-monotonically on $e_\mathrm{B}$. Contrary to previous claims, we find that $l_0$ is positive for most $e_\mathrm{B}$, implying that the binary receives net angular momentum, which may cause its separation to grow with time. The minimum $l_0$ occurs at intermediate $e_\mathrm{B}$ ($0.2 - 0.4$), corresponding to the regime where the inner eccentric disc is apsidally aligned with the binary. △ Less

Submitted 6 December, 2016; v1 submitted 23 October, 2016; originally announced October 2016.

Comments: 24 pages, 20 figures, accepted for publication in MNRAS

Planetary engulfment as a trigger for white dwarf pollution

Authors: Cristobal Petrovich, Diego J. Muñoz

Abstract: The presence of a planetary system can shield a planetesimal disk from the secular gravitational perturbations due to distant outer massive objects (planets or stellar companions). As the host star evolves off the main sequence to become a white dwarf, these planets can be engulfed, triggering secular instabilities and leading to the tidal disruptions of small rocky bodies. These disrupted bodies… ▽ More The presence of a planetary system can shield a planetesimal disk from the secular gravitational perturbations due to distant outer massive objects (planets or stellar companions). As the host star evolves off the main sequence to become a white dwarf, these planets can be engulfed, triggering secular instabilities and leading to the tidal disruptions of small rocky bodies. These disrupted bodies can feed the white dwarfs with rocky material and possibly explain the high-metallicity material in their atmospheres. We illustrate how this mechanism can operate when the gravitational perturbations are due to the Kozai-Lidov mechanism from a stellar binary companion. We show that this mechanism can explain the observed levels of accretion if: (1) the planetary engulfment happens fast compared to the secular timescale, which is generally the case for wide binaries ($>100$ AU) and planetary engulfment during the Asymptotic Giant Branch; (2) the planetesimal disk has a total mass of $\sim10^{-4}-10^{-2}M_\oplus$. We show that this new mechanism can provide a steady supply of material throughout the entire life of the white dwarfs for all cooling ages and can account for a large fraction (up to nearly half) of the observed polluted WDs. △ Less

Submitted 21 October, 2016; v1 submitted 17 July, 2016; originally announced July 2016.

Comments: 14 pages, 4 figures. Revised version after referee's comments. Major changes include a new Figure 1 and section 4.4

Pulsed Accretion onto Eccentric and Circular Binaries

Authors: Diego J. Muñoz, Dong Lai

Abstract: We present numerical simulations of circumbinary accretion onto eccentric and circular binaries using the moving-mesh code AREPO. This is the first set of simulations to tackle the problem of binary accretion using a finite-volume scheme on a freely moving mesh, which allows for accurate measurements of accretion onto individual stars for arbitrary binary eccentricity. While accretion onto a circu… ▽ More We present numerical simulations of circumbinary accretion onto eccentric and circular binaries using the moving-mesh code AREPO. This is the first set of simulations to tackle the problem of binary accretion using a finite-volume scheme on a freely moving mesh, which allows for accurate measurements of accretion onto individual stars for arbitrary binary eccentricity. While accretion onto a circular binary shows bursts with period of ~5 times the binary period P_b,accretion onto an eccentric binary is predominantly modulated at the period ~1P_b. For an equal-mass circular binary, the accretion rates onto individual stars are quite similar to each other, following the same variable pattern in time. By contrast, for eccentric binaries, one of the stars can accrete at a rate 10-20 times larger than its companion. This "symmetry breaking" between the stars, however, alternates over timescales of order 200 P_b, and can be attributed to a slowly precessing, eccentric circumbinary disk. Over longer timescales, the net accretion rates onto individual stars are the same, reaching a quasi-steady state with the circumbinary disk. These results have important implications for the accretion behavior of binary T-Tauri stars and supermassive binary black holes. △ Less

Submitted 31 March, 2016; originally announced April 2016.

Comments: 7 pages, 5 figures, submitted

The formation efficiency of close-in planets via Lidov-Kozai migration: analytic calculations

Authors: Diego J. Muñoz, Dong Lai, Bin Liu

Abstract: Lidov-Kozai oscillations of planets in stellar binaries, combined with tidal dissipation, can lead to the formation of hot Jupiters (HJs) or tidal disruption of planets. Recent population synthesis studies have found that the fraction of systems resulting in HJs (F_HJ) depends strongly on the planet mass, host stellar type and tidal dissipation strength, while the total migration fraction F_mig =… ▽ More Lidov-Kozai oscillations of planets in stellar binaries, combined with tidal dissipation, can lead to the formation of hot Jupiters (HJs) or tidal disruption of planets. Recent population synthesis studies have found that the fraction of systems resulting in HJs (F_HJ) depends strongly on the planet mass, host stellar type and tidal dissipation strength, while the total migration fraction F_mig = F_ HJ + F_dis (including both HJ formation and tidal disruption) exhibits much weaker dependence. We present an analytical method for calculating F_HJ and F_mig in the Lidov-Kozai migration scenario. The key ingredient of our method is to determine the critical initial planet-binary inclination angle that drives the planet to reach sufficiently large eccentricity for efficient tidal dissipation or disruption. This calculation includes the effects of the octupole potential and short-range forces on the planet. Our analytical method reproduces the resulting planet migration/disruption fractions from population synthesis, and can be easily implemented for various planet, stellar/companion types, and for different distributions of initial planetary semi-major axes, binary separations and eccentricities. We extend our calculations to planets in the super-Earth mass range and discuss the conditions for such planets to survive Lidov-Kozai migration and form close-in rocky planets. △ Less

Submitted 21 January, 2016; originally announced January 2016.

Comments: 8 pages, 5 figures. Submitted

Survival of Planets Around Shrinking Stellar Binaries

Authors: Diego J. Muñoz, Dong Lai

Abstract: The discovery of transiting circumbinary planets by the Kepler mission suggests that planets can form efficiently around binary stars. None of the stellar binaries currently known to host planets has a period shorter than 7 days, despite the large number of eclipsing binaries found in the Kepler target list with periods shorter than a few days. These compact binaries are believed to have evolved f… ▽ More The discovery of transiting circumbinary planets by the Kepler mission suggests that planets can form efficiently around binary stars. None of the stellar binaries currently known to host planets has a period shorter than 7 days, despite the large number of eclipsing binaries found in the Kepler target list with periods shorter than a few days. These compact binaries are believed to have evolved from wider orbits into their current configurations via the so-called Lidov-Kozai migration mechanism, in which gravitational perturbations from a distant tertiary companion induce large-amplitude eccentricity oscillations in the binary, followed by orbital decay and circularization due to tidal dissipation in the stars. Here we explore the orbital evolution of planets around binaries undergoing orbital decay by this mechanism. We show that planets may survive and become misaligned from their host binary, or may develop erratic behavior in eccentricity, resulting in their consumption by the stars or ejection from the system as the binary decays. Our results suggest that circumbinary planets around compact binaries could still exist, and we offer predictions as to what their orbital configurations should be like. △ Less

Submitted 2 July, 2015; v1 submitted 20 May, 2015; originally announced May 2015.

Comments: Accepted for publication. 10 pages, 7 figures in main text; 4 pages, 4 figures in supplementary material

Improving the convergence properties of the moving-mesh code AREPO

Authors: Ruediger Pakmor, Volker Springel, Andreas Bauer, Philip Mocz, Diego J. Munoz, Sebastian T. Ohlmann, Kevin Schaal, Chenchong Zhu

Abstract: Accurate numerical solutions of the equations of hydrodynamics play an ever more important role in many fields of astrophysics. In this work, we reinvestigate the accuracy of the moving-mesh code \textsc{Arepo} and show how its convergence order can be improved for general problems. In particular, we clarify that for certain problems \textsc{Arepo} only reaches first-order convergence for its orig… ▽ More Accurate numerical solutions of the equations of hydrodynamics play an ever more important role in many fields of astrophysics. In this work, we reinvestigate the accuracy of the moving-mesh code \textsc{Arepo} and show how its convergence order can be improved for general problems. In particular, we clarify that for certain problems \textsc{Arepo} only reaches first-order convergence for its original formulation. This can be rectified by simple modifications we propose to the time integration scheme and the spatial gradient estimates of the code, both improving the accuracy of the code. We demonstrate that the new implementation is indeed second-order accurate under the $L^1$ norm, and in particular substantially improves conservation of angular momentum. Interestingly, whereas these improvements can significantly change the results of smooth test problems, we also find that cosmological simulations of galaxy formation are unaffected, demonstrating that the numerical errors eliminated by the new formulation do not impact these simulations. In contrast, simulations of binary stars followed over a large number of orbital times are strongly affected, as here it is particularly crucial to avoid a long-term build up of errors in angular momentum conservation. △ Less

Submitted 21 December, 2015; v1 submitted 2 March, 2015; originally announced March 2015.

Comments: 11 pages, 9 figures, accepted for publication

Stellar orbit evolution in close circumstellar disc encounters

Authors: Diego J. Muñoz, Kaitlin M. Kratter, Mark Vogelsberger, Lars Hernquist, Volker Springel

Abstract: The formation and early evolution of circumstellar discs often occurs within dense, newborn stellar clusters. For the first time, we apply the moving-mesh code AREPO, to circumstellar discs in 3-D, focusing on disc-disc interactions that result from stellar fly-bys. Although a small fraction of stars are expected to undergo close approaches, the outcomes of the most violent encounters might leave… ▽ More The formation and early evolution of circumstellar discs often occurs within dense, newborn stellar clusters. For the first time, we apply the moving-mesh code AREPO, to circumstellar discs in 3-D, focusing on disc-disc interactions that result from stellar fly-bys. Although a small fraction of stars are expected to undergo close approaches, the outcomes of the most violent encounters might leave an imprint on the discs and host stars that will influence both their orbits and their ability to form planets. We first construct well-behaved 3-D models of self-gravitating discs, and then create a suite of numerical experiments of parabolic encounters, exploring the effects of pericenter separation r_p, disc orientation and disc-star mass ratio (M_d/M_*) on the orbital evolution of the host stars. Close encounters (2r_p<~ disc radius) can truncate discs on very short time scales. If discs are massive, close encounters facilitate enough orbital angular momentum extraction to induce stellar capture. We find that for realistic primordial disc masses M_d<~0.1M_*, non-colliding encounters induce minor orbital changes, which is consistent with analytic calculations of encounters in the linear regime. The same disc masses produce entirely different results for grazing/colliding encounters. In the latter case, rapidly cooling discs lose orbital energy by radiating away the energy excess of the shock-heated gas, thus causing capture of the host stars into a bound orbit. In rare cases, a tight binary with a circumbinary disc forms as a result of this encounter. △ Less

Submitted 7 January, 2015; v1 submitted 16 October, 2014; originally announced October 2014.

Comments: 20 pages, 14 figures, 1 table. Published by MNRAS

Suppression of extreme orbital evolution in triple systems with short range forces

Authors: Bin Liu, Diego J. Muñoz, Dong Lai

Abstract: The Lidov-Kozai (LK) mechanism plays an important role in the secular evolution of many hierarchical triple systems. The standard LK mechanism consists of large-amplitude oscillations in eccentricity and inclination of a binary subject to the quadrupole potential from an outer perturber. Recent work has shown that when the octupole terms are included in the potential, the inner binary can reach mo… ▽ More The Lidov-Kozai (LK) mechanism plays an important role in the secular evolution of many hierarchical triple systems. The standard LK mechanism consists of large-amplitude oscillations in eccentricity and inclination of a binary subject to the quadrupole potential from an outer perturber. Recent work has shown that when the octupole terms are included in the potential, the inner binary can reach more extreme eccentricities as well as undergo orientation flips. It is known that pericenter precessions due to short-range effects, such as General Relativity and tidal and rotational distortions, can limit the growth of eccentricity and even suppress standard (quadrupolar) LK oscillations, but their effect on the octupole-level LK mechanism has not been fully explored. In this paper, we systematically study how these short-range forces affect the extreme orbital behaviour found in octupole LK cycles. In general, the influence of the octupole potential is confined to a range of initial mutual inclinations Itot centered around 90deg (when the inner binary mass ratio is <<1), with this range expanding with increasing octupole strength. We find that, while the short-range forces do not change the width and location of this "window of influence", they impose a strict upper limit on the maximum achievable eccentricity. This limiting eccentricity can be calculated analytically, and its value holds even for strong octupole potential and for the general case of three comparable masses. Short-range forces also affect orbital flips, progressively reducing the range of Itot within which flips are possible as the intensity of these forces increases. △ Less

Submitted 7 January, 2015; v1 submitted 23 September, 2014; originally announced September 2014.

Comments: 20 pages,15 figures, updated to match version published on MNRAS

Planet-disc interaction on a freely moving mesh

Authors: Diego J. Munoz, Kaitlin M. Kratter, Volker Springel, Lars Hernquist

Abstract: General-purpose, moving-mesh schemes for hydrodynamics have opened the possibility of combining the accuracy of grid-based numerical methods with the flexibility and automatic resolution adaptivity of particle-based methods. Due to their supersonic nature, Keplerian accretion discs are in principle a very attractive system for applying such freely moving mesh techniques. However, the high degree o… ▽ More General-purpose, moving-mesh schemes for hydrodynamics have opened the possibility of combining the accuracy of grid-based numerical methods with the flexibility and automatic resolution adaptivity of particle-based methods. Due to their supersonic nature, Keplerian accretion discs are in principle a very attractive system for applying such freely moving mesh techniques. However, the high degree of symmetry of simple accretion disc models can be difficult to capture accurately by these methods, due to the generation of geometric grid noise and associated numerical diffusion, which is absent in polar grids. To explore these and other issues, in this work we study the idealized problem of two-dimensional planet-disc interaction with the moving-mesh code AREPO. We explore the hydrodynamic evolution of discs with planets through a series of numerical experiments that vary the planet mass, the disc viscosity and the mesh resolution, and compare the resulting surface density, vortensity field and tidal torque with results from the literature. We find that the performance of the moving-mesh code in this problem is in accordance with published results, showing good consistency with grid codes written in polar coordinates. We also conclude that grid noise and mesh distortions do not introduce excessive numerical diffusion. Finally, we show how the moving-mesh approach can naturally increase resolution in regions of high densityaround planets and planetary wakes, while retaining the background flow at low resolution. This provides an alternative to the difficult task of implementing adaptive mesh refinement in conventional polar-coordinate codes. △ Less

Submitted 7 January, 2015; v1 submitted 27 August, 2014; originally announced August 2014.

Comments: 21 pages, 15 figures, 2 tables. Updated to match version published by MNRAS

ALMA Observations of the T Tauri Binary System AS 205: Evidence for Molecular Winds and/or Binary Interactions

Authors: Colette Salyk, Klaus Pontoppidan, Stuartt Corder, Diego Muñoz, Ke Zhang, Geoffrey A. Blake

Abstract: In this study, we present high-resolution millimeter observations of the dust and gas disk of the T Tauri star AS 205 N and its companion, AS 205 S, obtained with the Atacama Large Millimeter Array. The gas disk around AS 205 N, for which infrared emission spectroscopy demonstrates significant deviations from Keplerian motion that has been interpreted as evidence for a disk wind (Pontoppidan et al… ▽ More In this study, we present high-resolution millimeter observations of the dust and gas disk of the T Tauri star AS 205 N and its companion, AS 205 S, obtained with the Atacama Large Millimeter Array. The gas disk around AS 205 N, for which infrared emission spectroscopy demonstrates significant deviations from Keplerian motion that has been interpreted as evidence for a disk wind (Pontoppidan et al. 2011; Bast et al. 2011), also displays significant deviations from Keplerian disk emission in the observations presented here. Detections near both AS 205 N and S are obtained in 1.3 mm continuum, 12CO 2-1, 13CO 2-1 and C18O 2-1. The 12CO emission is extended up to 2 arcsec from AS 205N, and both 12CO and 13CO display deviations from Keplerian rotation at all angular scales. Two possible explanations for these observations hold up best to close scrutiny - tidal interaction with AS 205 S or disk winds (or a combination of the two), and we discuss these possibilities in some detail. △ Less

Submitted 9 July, 2014; originally announced July 2014.

Comments: accepted by The Astrophysical Journal

Multi-Dimensional, Compressible Viscous Flow on a Moving Voronoi Mesh

Authors: Diego Muñoz, Volker Springel, Robert Marcus, Mark Vogelsberger, Lars Hernquist

Abstract: Numerous formulations of finite volume schemes for the Euler and Navier-Stokes equations exist, but in the majority of cases they have been developed for structured and stationary meshes. In many applications, more flexible mesh geometries that can dynamically adjust to the problem at hand and move with the flow in a (quasi) Lagrangian fashion would, however, be highly desirable, as this can allow… ▽ More Numerous formulations of finite volume schemes for the Euler and Navier-Stokes equations exist, but in the majority of cases they have been developed for structured and stationary meshes. In many applications, more flexible mesh geometries that can dynamically adjust to the problem at hand and move with the flow in a (quasi) Lagrangian fashion would, however, be highly desirable, as this can allow a significant reduction of advection errors and an accurate realization of curved and moving boundary conditions. Here we describe a novel formulation of viscous continuum hydrodynamics that solves the equations of motion on a Voronoi mesh created by a set of mesh-generating points. The points can move in an arbitrary manner, but the most natural motion is that given by the fluid velocity itself, such that the mesh dynamically adjusts to the flow. Owing to the mathematical properties of the Voronoi tessellation, pathological mesh-twisting effects are avoided. Our implementation considers the full Navier-Stokes equations and has been realized in the AREPO code both in 2D and 3D. We propose a new approach to compute accurate viscous fluxes for a dynamic Voronoi mesh, and use this to formulate a finite volume solver of the Navier-Stokes equations. Through a number of test problems, including circular Couette flow and flow past a cylindrical obstacle, we show that our new scheme combines good accuracy with geometric flexibility, and hence promises to be competitive with other highly refined Eulerian methods. This will in particular allow astrophysical applications of the AREPO code where physical viscosity is important, such as in the hot plasma in galaxy clusters, or for viscous accretion disk models. △ Less

Submitted 5 March, 2012; originally announced March 2012.

Comments: 26 pages, 21 figures. Submitted to MNRAS

The Circular Polarization of Sagittarius A* at Submillimeter Wavelengths

Authors: Diego J. Muñoz, Daniel P. Marrone, James M. Moran, Ramprasad Rao

Abstract: We report the first detections of circularly polarized emission at submillimeter wavelengths from the compact radio source and supermassive black hole candidate Sgr A* at a level of 1.2\pm0.3% at 1.3 mm wavelength (230 GHz) and 1.6\pm0.3% at 860 microns (345 GHz) with the same handedness as observed at lower frequencies (1.4-15 GHz). The observations, taken with the Submillimeter Array in multiple… ▽ More We report the first detections of circularly polarized emission at submillimeter wavelengths from the compact radio source and supermassive black hole candidate Sgr A* at a level of 1.2\pm0.3% at 1.3 mm wavelength (230 GHz) and 1.6\pm0.3% at 860 microns (345 GHz) with the same handedness as observed at lower frequencies (1.4-15 GHz). The observations, taken with the Submillimeter Array in multiple epochs, also show simultaneous linear polarization (LP) at both wavelengths of about 6%. These properties differ sharply from those at wavelengths longer than 1 cm (frequencies below 30 GHz), where weak circular polarization (CP) (~ 0.5%) dominates over LP, which is not detected at similar fractional limits. We describe an extensive set of tests to ensure the accuracy of our measurements. We find no CP in any other source, including the bright quasar 1924-292, which traces the same path on the sky as Sgr A* and therefore should be subject to identical systematic errors originating in the instrument frame. Since a relativistic synchrotron plasma is expected to produce little CP, the observed CP is probably generated close to the event horizon by the Faraday conversion process. We use a simple model to show that the phase shift associated with Faraday conversion can be nearly independent of frequency, a sufficient condition to make the handedness of CP independent of frequency. Because the size of the tau=1-surface changes by more than an order of magnitude between 1.4 and 345 GHz, the magnetic field must be coherent over such scales to consistently produce left CP. △ Less

Submitted 5 January, 2012; v1 submitted 2 May, 2011; originally announced May 2011.

Comments: 13, pages, 6 figures, 5 tables. Accepted for publication in ApJ

An X-ray, IR, and Submillimeter Flare of Sagittarius A*

Authors: D. P. Marrone, F. K. Baganoff, M. R. Morris, J. M. Moran, A. M. Ghez, S. D. Hornstein, C. D. Dowell, D. J. Munoz, M. W. Bautz, G. R. Ricker, W. N. Brandt, G. P. Garmire, J. R. Lu, K. Matthews, J. -H. Zhao, R. Rao, G. C. Bower

Abstract: Energetic flares are observed in the Galactic supermassive black hole Sagittarius A* from radio to X-ray wavelengths. On a few occasions, simultaneous flares have been detected in IR and X-ray observations, but clear counterparts at longer wavelengths have not been seen. We present a flare observed over several hours on 2006 July 17 with the Chandra X-Ray Observatory, the Keck II telescope, the… ▽ More Energetic flares are observed in the Galactic supermassive black hole Sagittarius A* from radio to X-ray wavelengths. On a few occasions, simultaneous flares have been detected in IR and X-ray observations, but clear counterparts at longer wavelengths have not been seen. We present a flare observed over several hours on 2006 July 17 with the Chandra X-Ray Observatory, the Keck II telescope, the Caltech Submillimeter Observatory, and the Submillimeter Array. All telescopes observed strong flare events, but the submillimeter peak is found to occur nearly 100 minutes after the X-ray peak. Submillimeter polarization data show linear polarization in the excess flare emission, increasing from 9% to 17% as the flare passes through its peak, consistent with a transition from optically thick to thin synchrotron emission. The temporal and spectral behavior of the flare require that the energetic electrons responsible for the emission cool faster than expected from their radiative output. This is consistent with adiabatic cooling in an expanding emission region, with X-rays produced through self-Compton scattering, although not consistent with the simplest model of such expansion. We also present a submillimeter flare that followed a bright IR flare on 2005 July 31. Compared to 2006, this event had a larger peak IR flux and similar submillimeter flux, but it lacked measurable X-ray emission. It also showed a shorter delay between the IR and submillimeter peaks. Based on these events we propose a synchrotron and self-Compton model to relate the submillimeter lag and the variable IR/X-ray luminosity ratio. △ Less

Submitted 14 July, 2008; v1 submitted 18 December, 2007; originally announced December 2007.

Comments: To appear in The Astrophysical Journal 682: 373, 2008 July 20. Corrected in response to referee comments, matches published version

Journal ref: Astrophys. J.682:373-383, 2008

Massive Clumps in the NGC 6334 Star Forming Region

Authors: Diego J. Munoz, Diego Mardones, Guido Garay, David Rebolledo, Kate Brooks, Sylvain Bontemps

Abstract: We report observations of dust continuum emission at 1.2 mm toward the star forming region NGC 6334 made with the SEST SIMBA bolometer array. The observations cover an area of $\sim 2$ square degrees with approximately uniform noise. We detected 181 clumps spanning almost three orders of magnitude in mass (3\Msun$-6\times10^3$ \Msun) and with sizes in the range 0.1--1.0 pc. We find that the clum… ▽ More We report observations of dust continuum emission at 1.2 mm toward the star forming region NGC 6334 made with the SEST SIMBA bolometer array. The observations cover an area of $\sim 2$ square degrees with approximately uniform noise. We detected 181 clumps spanning almost three orders of magnitude in mass (3\Msun$-6\times10^3$ \Msun) and with sizes in the range 0.1--1.0 pc. We find that the clump mass function $dN/d\log M$ is well fit with a power law of the mass with exponent -0.6 (or equivalently $dN/dM \propto M^{-1.6}$). The derived exponent is similar to those obtained from molecular line emission surveys and is significantly different from that of the stellar initial mass function. We investigated changes in the mass spectrum by changing the assumptions on the temperature distribution of the clumps and on the contribution of free-free emission to the 1.2 mm emission, and found little changes on the exponent. The Cumulative Mass Distribution Function is also analyzed giving consistent results in a mass range excluding the high-mass end where a power-law fit is no longer valid. The masses and sizes of the clumps observed in NGC 6334 indicate that they are not direct progenitors of stars and that the process of fragmentation determines the distribution of masses later on or occurs at smaller spatial scales. The spatial distribution of the clumps in NGC 6334 reveals clustering which is strikingly similar to that exhibited by young stars in other star forming regions. A power law fit to the surface density of companions gives $Σ\propto θ^{-0.62}$. △ Less

Submitted 20 June, 2007; originally announced June 2007.

Comments: 16 pages, 11 figures, 4 tables. To appear in the Astrophysical Journal

Exploring a deep meridional flow hypothesis for a circulation dominated solar dynamo model

Authors: G. A. Guerrero, J. D. Muñoz, E M. de Gouveia Dal Pino

Abstract: Circulation-dominated solar dynamo models, which employ a helioseismic rotation profile and a fixed meridional flow, give a good approximation to the large scale solar magnetic phenomena, such as the 11-year cycle or the so called Hale's law of polarities. Nevertheless, the larger amplitude of the radial shear at the high latitudes makes the dynamo to produce a strong toroidal magnetic field at… ▽ More Circulation-dominated solar dynamo models, which employ a helioseismic rotation profile and a fixed meridional flow, give a good approximation to the large scale solar magnetic phenomena, such as the 11-year cycle or the so called Hale's law of polarities. Nevertheless, the larger amplitude of the radial shear at the high latitudes makes the dynamo to produce a strong toroidal magnetic field at high latitudes, in contradiction with the observations of the sunspots (Sporer`s Law). A possible solution was proposed by Nandy and Choudhur (2002) in which a deep meridional flow can conduct the magnetic field inside of a stable layer (the radiative core) and then allow that it erupts just at lower latitudes. We have recently explored this hypothesis in a different approximation, using the magnetic buoyancy mechanism proposed by Dikpaty and Charbonneau (1999) and found that a deep meridional flow pushes the maximum of the toroidal magnetic field towards the solar equator, but, in contrast to Nandy and Choudhuri (2002), a second zone of maximum fields remains at the poles (Guerrero and Munoz (2004)). We here review these results and also discuss a new possible scenario where the tachocline has an ellipsoidal shape, following early helioseismologic observations Charbonneau et. al., (1999), and find that the modification of the geometry of the tachocline can lead to results which are in good agreement with observations and opens the possibility to explore in more detail, through the dynamo model, the place where the magnetic field could be really stored. △ Less

Submitted 4 April, 2005; v1 submitted 4 April, 2005; originally announced April 2005.

Comments: 8 pages, 6 (low resolution) figures, Procs. of Conference Magnetic Fields in the Universe: from laboratory and stars to primordial structures, AIP (NY), eds. E. M. de Gouveia Dal Pino, G. Lugones and A. Lazarian

Kinematic solar dynamo models with a deep meridional flow

Authors: G. A. Guerrero, J. D. Muñoz

Abstract: We develop two different solar dynamo models to verify the hypothesis that a deep meridional flow can restrict the apperance of sunspots below 45 degrees, proposed by Nandy & Choudhuri (2002). In the first one, a single polytropic approximation for the density profile was taken, for both radiative and convective zones. In the second one, two polytropes were used to distinguish between both zones… ▽ More We develop two different solar dynamo models to verify the hypothesis that a deep meridional flow can restrict the apperance of sunspots below 45 degrees, proposed by Nandy & Choudhuri (2002). In the first one, a single polytropic approximation for the density profile was taken, for both radiative and convective zones. In the second one, two polytropes were used to distinguish between both zones Pinzon & Calvo-Mozo (2001). The magnetic buoyancy mechanism proposed by Dikpati & Charbonneau (1999) was chosen in both models. We, actually, have obtained that a deep meridional flow pushes the maxima of toroidal magnetic field toward the solar equator, but in contrast to Nandy & Choudhuri (2002) a second zone of maximal fields remains at the poles. The second model, although closely resembling the solar standard model of Bahcall, Pinsonneault & Wasserbug (1995); Bahcall, Pinsonneault & Basu (2001), gives solar cyles three times longer than observed. △ Less

Submitted 4 February, 2004; originally announced February 2004.

Comments: 7 pages, 6 figures, to be published in MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc. 350 (2004) 317