Showing 1–50 of 137 results for author: Do, T

Using Galaxy Evolution as Source of Physics-Based Ground Truth for Generative Models

Authors: Yun Qi Li, Tuan Do, Evan Jones, Bernie Boscoe, Kevin Alfaro, Zooey Nguyen

Abstract: Generative models producing images have enormous potential to advance discoveries across scientific fields and require metrics capable of quantifying the high dimensional output. We propose that astrophysics data, such as galaxy images, can test generative models with additional physics-motivated ground truths in addition to human judgment. For example, galaxies in the Universe form and change ove… ▽ More Generative models producing images have enormous potential to advance discoveries across scientific fields and require metrics capable of quantifying the high dimensional output. We propose that astrophysics data, such as galaxy images, can test generative models with additional physics-motivated ground truths in addition to human judgment. For example, galaxies in the Universe form and change over billions of years, following physical laws and relationships that are both easy to characterize and difficult to encode in generative models. We build a conditional denoising diffusion probabilistic model (DDPM) and a conditional variational autoencoder (CVAE) and test their ability to generate realistic galaxies conditioned on their redshifts (galaxy ages). This is one of the first studies to probe these generative models using physically motivated metrics. We find that both models produce comparable realistic galaxies based on human evaluation, but our physics-based metrics are better able to discern the strengths and weaknesses of the generative models. Overall, the DDPM model performs better than the CVAE on the majority of the physics-based metrics. Ultimately, if we can show that generative models can learn the physics of galaxy evolution, they have the potential to unlock new astrophysical discoveries. △ Less

Submitted 9 July, 2024; originally announced July 2024.

Comments: 20 pages, 14 figures, 1 Table, code: https://github.com/astrodatalab/li2024_public, training data: https://zenodo.org/records/11117528

General Relativistic effects and the NIR variability of Sgr A* II: A systematic approach to temporal asymmetry

Authors: Sebastiano D. von Fellenberg, Gunther Witzel, Michi Bauboeck, Hui-Hsuan Chung, Nicola Marchili, Greg Martinez, Matteo Sadun-Bordoni, Guillaume Bourdarot, Tuan Do, Antonia Drescher, Giovanni Fazio, Frank Eisenhauer, Reinhard Genzel, Stefan Gillessen, Joseph L. Hora, Felix Mang, Thomas Ott, Howard A. Smith, Eduardo Ros, Diogo C. Ribeiro, Felix Widmann, S. P. Willner, J. Anton Zensus

Abstract: A systematic study, based on the third-moment structure function, of Sgr A*'s variability finds an exponential rise time $τ_{1,\rm{obs}}=14.8^{+0.4}_{-1.5}~\mathrm{minutes}$ and decay time $τ_{2,\rm{obs}}=13.1^{+1.3}_{-1.4}~\mathrm{minutes}$. This symmetry of the flux-density variability is consistent with earlier work, and we interpret it as caused by the dominance of Doppler boosting, as opposed… ▽ More A systematic study, based on the third-moment structure function, of Sgr A*'s variability finds an exponential rise time $τ_{1,\rm{obs}}=14.8^{+0.4}_{-1.5}~\mathrm{minutes}$ and decay time $τ_{2,\rm{obs}}=13.1^{+1.3}_{-1.4}~\mathrm{minutes}$. This symmetry of the flux-density variability is consistent with earlier work, and we interpret it as caused by the dominance of Doppler boosting, as opposed to gravitational lensing, in Sgr~A*'s light curve. A relativistic, semi-physical model of Sgr~A* confirms an inclination angle $i<45$ degrees. The model also shows that the emission of the intrinsic radiative process can have some asymmetry even though the observed emission does not. The third-moment structure function, which is a measure of the skewness of the light-curve increments, may be a useful summary statistic in other contexts of astronomy because it senses only temporal asymmetry, i.e., it averages to zero for any temporally symmetric signal. △ Less

Submitted 9 July, 2024; originally announced July 2024.

Comments: Accepted for publication in A&A letters

Photometric Classification of Stars Around the Milky Way's Central Black Hole: I. Central Parsec

Authors: E. Gallego-Cano, T. Fritz, R. Schödel, A. Feldmeier-Krause, T. Do, S. Nishiyama

Abstract: The presence of young massive stars in the Galactic Centre (GC) raises questions about star formation near the black hole Sagittarius A* (Sgr A*). Additionally, the initial mass function (IMF) in this region appears different from the standard Salpeter/Kroupa law. Extreme extinction and crowding limit our understanding of the stellar population, with spectroscopic data available only for selected… ▽ More The presence of young massive stars in the Galactic Centre (GC) raises questions about star formation near the black hole Sagittarius A* (Sgr A*). Additionally, the initial mass function (IMF) in this region appears different from the standard Salpeter/Kroupa law. Extreme extinction and crowding limit our understanding of the stellar population, with spectroscopic data available only for selected bright sources. We aim to improve knowledge about the distribution and IMF of young, massive stars near Sgr A*. Using intermediate band (IB) photometry, we identify candidates for massive young stars through Bayesian inference, a neural network, and a gradient-boosted trees algorithm. We obtained spectral energy distributions for 6590 stars, 1181 of which have been previously classified spectroscopically. We identify 351 stars classified as early types by all three methods, including 155 newly identified candidates. The radial density profiles for late and early-type stars fit broken power laws, with a break radius of 9.2 +- 0.6'' for early-type stars. Late-type stars show a core-like distribution around Sgr A*, while early-type stars' density increases steeply towards the black hole. We infer a top-heavy IMF of young stars near Sgr A* (R < 9''), with a power-law of 1.6 +- 0.1. At greater distances, a standard Salpeter/Kroupa IMF fits the data. IB photometry also constrains the metallicities of late-type stars, estimating metallicities for over 600 stars. The IMF variation with radial distance suggests different star formation mechanisms, with a top-heavy IMF near Sgr A* consistent with disc star formation. △ Less

Submitted 13 June, 2024; v1 submitted 6 June, 2024; originally announced June 2024.

Orbital and Atmospheric Characterization of the 1RXS J034231.8+121622 System Using High-Resolution Spectroscopy Confirms That The Companion is a Low-Mass Star

Authors: Clarissa R. Do Ó, Ben Sappey, Quinn M. Konopacky, Jean-Baptiste Ruffio, Kelly K. O'Neil, Tuan Do, Gregory Martinez, Travis S. Barman, Jayke S. Nguyen, Jerry W. Xuan, Christopher A. Theissen, Sarah Blunt, William Thompson, Chih-Chun Hsu, Ashley Baker, Randall Bartos, Geoffrey A. Blake, Benjamin Calvin, Sylvain Cetre, Jacques-Robert Delorme, Greg Doppmann, Daniel Echeverri, Luke Finnerty, Michael P. Fitzgerald, Julie Inglis , et al. (11 additional authors not shown)

Abstract: The 1RXS J034231.8+121622 system consists of an M dwarf primary and a directly imaged low-mass stellar companion. We use high resolution spectroscopic data from Keck/KPIC to estimate the objects' atmospheric parameters and radial velocities (RVs). Using PHOENIX stellar models, we find that the primary has a temperature of 3460 $\pm$ 50 K a metallicity of 0.16 $\pm$ 0.04, while the secondary has a… ▽ More The 1RXS J034231.8+121622 system consists of an M dwarf primary and a directly imaged low-mass stellar companion. We use high resolution spectroscopic data from Keck/KPIC to estimate the objects' atmospheric parameters and radial velocities (RVs). Using PHOENIX stellar models, we find that the primary has a temperature of 3460 $\pm$ 50 K a metallicity of 0.16 $\pm$ 0.04, while the secondary has a temperature of 2510 $\pm$ 50 K and a metallicity of $0.13\substack{+0.12 \\ -0.11}$. Recent work suggests this system is associated with the Hyades, placing it an older age than previous estimates. Both metallicities agree with current $[Fe/H]$ Hyades measurements (0.11 -- 0.21). Using stellar evolutionary models, we obtain significantly higher masses for the objects, of 0.30 $\pm$ 0.15 $M_\odot$ and 0.08 $\pm$ 0.01 $M_\odot$ (84 $\pm$ 11 $M_{Jup}$) respectively. Using the RVs and a new astrometry point from Keck/NIRC2, we find that the system is likely an edge-on, moderately eccentric ($0.41\substack{+0.27 \\ -0.08}$) configuration. We also estimate the C/O ratio of both objects using custom grid models, obtaining 0.42 $\pm$ 0.10 (primary) and 0.55 $\pm$ 0.10 (companion). From these results, we confirm that this system most likely went through a binary star formation process in the Hyades. The significant changes in this system's parameters since its discovery highlight the importance of high resolution spectroscopy for both orbital and atmospheric characterization of directly imaged companions. △ Less

Submitted 11 April, 2024; originally announced April 2024.

Comments: 30 pages, 18 figures, accepted for publication in AJ

Stability investigations of de Sitter inflationary solutions in power-law extensions of the Starobinsky model

Authors: Tuyen M. Pham, Duy H. Nguyen, Tuan Q. Do, W. F. Kao

Abstract: In this paper, we would like to examine whether stable de Sitter inflationary solutions appear within power-law extensions of the Starobinsky model. In particular, we will address general constraints for the existence along with the stability of de Sitter inflationary solutions in a general case involving not only the Starobinsky $R^2$ term but also an additional power-law $R^n$ one. According to… ▽ More In this paper, we would like to examine whether stable de Sitter inflationary solutions appear within power-law extensions of the Starobinsky model. In particular, we will address general constraints for the existence along with the stability of de Sitter inflationary solutions in a general case involving not only the Starobinsky $R^2$ term but also an additional power-law $R^n$ one. According to the obtained results, we will be able to identify which extension is more suitable for an early inflationary phase rather than a late-time cosmic acceleration phase. To be more specific, we will consider several values of $n$ to see whether the corresponding de Sitter inflationary solutions are stable or not. △ Less

Submitted 4 March, 2024; originally announced March 2024.

Comments: 31 pages, 16 figures. Comments are welcome

An Estimate of the Binary Star Fraction Among Young Stars at the Galactic Center: Possible Evidence of a Radial Dependence

Authors: Abhimat K. Gautam, Tuan Do, Andrea M. Ghez, Devin S. Chu, Matthew W. Hosek Jr., Shoko Sakai, Smadar Naoz, Mark R. Morris, Anna Ciurlo, Zoe Haggard, Jessica R. Lu

Abstract: We present the first estimate of the intrinsic binary fraction of young stars across the central $\approx$ 0.4 pc surrounding the supermassive black hole (SMBH) at the Milky Way Galactic center (GC). This experiment searched for photometric variability in 102 young stars, using 119 nights of 10"-wide adaptive optics imaging observations taken at Keck Observatory over 16 years in the K'- and H-band… ▽ More We present the first estimate of the intrinsic binary fraction of young stars across the central $\approx$ 0.4 pc surrounding the supermassive black hole (SMBH) at the Milky Way Galactic center (GC). This experiment searched for photometric variability in 102 young stars, using 119 nights of 10"-wide adaptive optics imaging observations taken at Keck Observatory over 16 years in the K'- and H-bands. We photometrically detected three binary stars, all of which are situated more than 1" (0.04 pc) from the SMBH and one of which, S2-36, is newly reported here with spectroscopic confirmation. To convert the observed binary fraction into an estimate of the underlying binary fraction, we determined experiment sensitivity through detailed light curve simulations, incorporating photometric effects of eclipses, irradiation, and tidal distortion in binaries. The simulations assumed a population of young binaries, with stellar ages (4 Myr) and masses matched to the most probable values measured for the GC young star population and underlying binary system parameters similar to those of local massive stars. The detections and simulations imply young, massive stars in the GC have a stellar binary fraction $\geq$ 71% (68% confidence), or $\geq$ 42% (95% confidence). This inferred GC young star binary fraction is consistent with that typically seen in young stellar populations in the solar neighborhood. Furthermore, our measured binary fraction is significantly higher than that recently reported by Chu et al. (2023) based on RV measurements of young stars <~1" of the SMBH. Constrained with these two studies, the probability that the same underlying young binary fraction extends across the entire region is <1.4%. This tension provides support for a radial dependence of the binary star fraction and, therefore, for the dynamical predictions of binary merger and evaporation events close to the SMBH. △ Less

Submitted 7 February, 2024; v1 submitted 23 January, 2024; originally announced January 2024.

Comments: 51 pages, 27 figures, 7 tables. Accepted for publication in The Astrophysical Journal. Abstract revised for final manuscript

High-Resolution, Mid-Infrared Color Temperature Mapping of the Central 10 Arcseconds of the Galaxy

Authors: Cuc K. Dinh, Anna Cuirlo, Mark R. Morris, Rainer Schödel, Andrea Ghez, Tuan Do

Abstract: The neighborhood of the Galactic black hole boasts a plethora of extended interstellar gas and dust features as well as populations of compact (unresolved, or marginally resolved) features such as the G objects. Most are well manifested in the infrared. To disentangle and characterize the infrared structure of extended features and identify compact sources, we used 3.8~$μ$m (L' filter) data from t… ▽ More The neighborhood of the Galactic black hole boasts a plethora of extended interstellar gas and dust features as well as populations of compact (unresolved, or marginally resolved) features such as the G objects. Most are well manifested in the infrared. To disentangle and characterize the infrared structure of extended features and identify compact sources, we used 3.8~$μ$m (L' filter) data from the NIRC2 imager at the Keck Observatory and 8.6~$μ$m (PAH1 filter) data from the VISIR imager at the Very Large Telescope (VLT) to produce the highest-resolution mid-IR color-temperature map of the inner half-parsec of the Galactic Center to date. From this map, we compile a catalog of features that stand out from their background. In particular, we identify 33 compact sources that stand out against the local background temperature, 11 of which are newly identified and are candidates for being members of the G objects population. Additionally, we resolve and newly characterize the morphology of several known extended features. These results prepare the way for ongoing and future JWST studies that have access to a greater range of mid-infrared wavelengths, and thus will allow for refined estimation of the trends of dust temperatures. △ Less

Submitted 27 November, 2023; originally announced November 2023.

Comments: 22 pages, 14 figures

Magnetogenesis from Anisotropic Universe

Authors: Sourav Pal, Debaprasad Maity, Tuan Q. Do

Abstract: The existence of large-scale anisotropy can not be ruled out by the cosmic microwave background (CMB) radiation. Over the years, several models have been proposed in the context of anisotropic inflation to account for CMB's cold spot and hemispheric asymmetry. However, any small-scale anisotropy, if exists during inflation, is not constrained due to its nonlinear evolution in the subsequent phase.… ▽ More The existence of large-scale anisotropy can not be ruled out by the cosmic microwave background (CMB) radiation. Over the years, several models have been proposed in the context of anisotropic inflation to account for CMB's cold spot and hemispheric asymmetry. However, any small-scale anisotropy, if exists during inflation, is not constrained due to its nonlinear evolution in the subsequent phase. This small-scale anisotropy during inflation can play a non-trivial role in giving rise to the cosmic magnetic field, which is the subject of our present study. Assuming a particular phenomenological form of an anisotropic inflationary universe, we have shown that it can generate a large-scale magnetic field at $1$-Mpc scale with a magnitude $\sim 4\times 10^{-20}~G$, within the observed bound. Because of the anisotropy, the conformal flatness property is lost, and the Maxwell field is generated even without explicit coupling. This immediately resolves the strong coupling problem in the standard magnetogenesis scenario. In addition, assuming very low conductivity during the reheating era, we can further observe the evolution of the electromagnetic field with the equation of state (EoS) $ω_{eff}$ and its effects on the present-day magnetic field. △ Less

Submitted 13 November, 2023; originally announced November 2023.

Comments: 15 pages, 5 figures. Comments are welcome

The JWST Galactic Center Survey -- A White Paper

Authors: Rainer Schoedel, Steve Longmore, Jonny Henshaw, Adam Ginsburg, John Bally, Anja Feldmeier, Matt Hosek, Francisco Nogueras Lara, Anna Ciurlo, Mélanie Chevance, J. M. Diederik Kruijssen, Ralf Klessen, Gabriele Ponti, Pau Amaro-Seoane, Konstantina Anastasopoulou, Jay Anderson, Maria Arias, Ashley T. Barnes, Cara Battersby, Giuseppe Bono, Lucía Bravo Ferres, Aaron Bryant, Miguel Cano Gonzáalez, Santi Cassisi, Leonardo Chaves-Velasquez , et al. (85 additional authors not shown)

Abstract: The inner hundred parsecs of the Milky Way hosts the nearest supermassive black hole, largest reservoir of dense gas, greatest stellar density, hundreds of massive main and post main sequence stars, and the highest volume density of supernovae in the Galaxy. As the nearest environment in which it is possible to simultaneously observe many of the extreme processes shaping the Universe, it is one of… ▽ More The inner hundred parsecs of the Milky Way hosts the nearest supermassive black hole, largest reservoir of dense gas, greatest stellar density, hundreds of massive main and post main sequence stars, and the highest volume density of supernovae in the Galaxy. As the nearest environment in which it is possible to simultaneously observe many of the extreme processes shaping the Universe, it is one of the most well-studied regions in astrophysics. Due to its proximity, we can study the center of our Galaxy on scales down to a few hundred AU, a hundred times better than in similar Local Group galaxies and thousands of times better than in the nearest active galaxies. The Galactic Center (GC) is therefore of outstanding astrophysical interest. However, in spite of intense observational work over the past decades, there are still fundamental things unknown about the GC. JWST has the unique capability to provide us with the necessary, game-changing data. In this White Paper, we advocate for a JWST NIRCam survey that aims at solving central questions, that we have identified as a community: i) the 3D structure and kinematics of gas and stars; ii) ancient star formation and its relation with the overall history of the Milky Way, as well as recent star formation and its implications for the overall energetics of our galaxy's nucleus; and iii) the (non-)universality of star formation and the stellar initial mass function. We advocate for a large-area, multi-epoch, multi-wavelength NIRCam survey of the inner 100\,pc of the Galaxy in the form of a Treasury GO JWST Large Program that is open to the community. We describe how this survey will derive the physical and kinematic properties of ~10,000,000 stars, how this will solve the key unknowns and provide a valuable resource for the community with long-lasting legacy value. △ Less

Submitted 14 March, 2024; v1 submitted 18 October, 2023; originally announced October 2023.

Comments: This White Paper will be updated when required (e.g. new authors joining, editing of content). Most recent update: 24 Oct 2023

The development of HISPEC for Keck and MODHIS for TMT: science cases and predicted sensitivities

Authors: Quinn M. Konopacky, Ashley D. Baker, Dimitri Mawet, Michael P. Fitzgerald, Nemanja Jovanovic, Charles Beichman, Garreth Ruane, Rob Bertz, Hiroshi Terada, Richard Dekany, Larry Lingvay, Marc Kassis, David Anderson, Motohide Tamura, Bjorn Benneke, Thomas Beatty, Tuan Do, Shogo Nishiyama, Peter Plavchan, Jason Wang, Ji Wang, Adam Burgasser, Jean-Baptiste Ruffio, Huihao Zhang, Aaron Brown , et al. (50 additional authors not shown)

Abstract: HISPEC is a new, high-resolution near-infrared spectrograph being designed for the W.M. Keck II telescope. By offering single-shot, R=100,000 between 0.98 - 2.5 um, HISPEC will enable spectroscopy of transiting and non-transiting exoplanets in close orbits, direct high-contrast detection and spectroscopy of spatially separated substellar companions, and exoplanet dynamical mass and orbit measureme… ▽ More HISPEC is a new, high-resolution near-infrared spectrograph being designed for the W.M. Keck II telescope. By offering single-shot, R=100,000 between 0.98 - 2.5 um, HISPEC will enable spectroscopy of transiting and non-transiting exoplanets in close orbits, direct high-contrast detection and spectroscopy of spatially separated substellar companions, and exoplanet dynamical mass and orbit measurements using precision radial velocity monitoring calibrated with a suite of state-of-the-art absolute and relative wavelength references. MODHIS is the counterpart to HISPEC for the Thirty Meter Telescope and is being developed in parallel with similar scientific goals. In this proceeding, we provide a brief overview of the current design of both instruments, and the requirements for the two spectrographs as guided by the scientific goals for each. We then outline the current science case for HISPEC and MODHIS, with focuses on the science enabled for exoplanet discovery and characterization. We also provide updated sensitivity curves for both instruments, in terms of both signal-to-noise ratio and predicted radial velocity precision. △ Less

Submitted 19 September, 2023; originally announced September 2023.

Comments: 25 pages, 9 figures. To appear in the Proceedings of SPIE: Techniques and Instrumentation for Detection of Exoplanets XI, vol. 12680 (2023)

Anisotropic power-law inflation for models of non-canonical scalar fields non-minimally coupled to a two-form field

Authors: Tuyen M. Pham, Duy H. Nguyen, Tuan Q. Do, W. F. Kao

Abstract: In this paper, we investigate the validity of the so-called cosmic no-hair conjecture in the framework of anisotropic inflation models of non-canonical scalar fields non-minimally coupled to a two-form field. In particular, we focus on two typical {\it k}-inflation and Dirac-Born-Infeld inflation models, in which we find a set of exact anisotropic power-law inflationary solutions. Interestingly, t… ▽ More In this paper, we investigate the validity of the so-called cosmic no-hair conjecture in the framework of anisotropic inflation models of non-canonical scalar fields non-minimally coupled to a two-form field. In particular, we focus on two typical {\it k}-inflation and Dirac-Born-Infeld inflation models, in which we find a set of exact anisotropic power-law inflationary solutions. Interestingly, these solutions are shown to be stable and attractive during an inflationary phase using the dynamical system analysis. The obtained results indicate that the non-minimal coupling between the scalar and two-form fields acts as a non-trivial source of generating stable spatial anisotropies during the inflationary phase and therefore violates the prediction of the cosmic no-hair conjecture, even when the scalar field is of non-canonical forms. In connection with the Planck 2018 data, tensor-to-scalar ratios of these anisotropic solutions are investigated. As a result, it appears that the tensor-to-scalar ratio of the anisotropic power-law inflationary solution of {\it k}-inflation model turns out to be more highly consistent with the Planck 2018 data than that of Dirac-Born-Infeld model. △ Less

Submitted 4 February, 2024; v1 submitted 5 September, 2023; originally announced September 2023.

Comments: 24 pages, 8 figures. Final version with investigations on the tensor-to-scalar ratio added. Some discussions and relevant references are also added. All calculations remain unchanged. Accepted for publication in EPJC. Comments are welcome

Journal ref: Eur. Phys. J. C 84, 105 (2024)

Near-Infrared Flux Distribution of Sgr A* from 2005-2022: Evidence for an Enhanced Accretion Episode in 2019

Authors: Grant C. Weldon, Tuan Do, Gunther Witzel, Andrea M. Ghez, Abhimat K. Gautam, Eric E. Becklin, Mark R. Morris, Gregory D. Martinez, Shoko Sakai, Jessica R. Lu, Keith Matthews, Matthew W. Hosek Jr., Zoë Haggard

Abstract: Sgr A* is the variable electromagnetic source associated with accretion onto the Galactic center supermassive black hole. While the near-infrared (NIR) variability of Sgr A* was shown to be consistent over two decades, unprecedented activity in 2019 challenges existing statistical models. We investigate the origin of this activity by re-calibrating and re-analyzing all of our Keck Observatory Sgr… ▽ More Sgr A* is the variable electromagnetic source associated with accretion onto the Galactic center supermassive black hole. While the near-infrared (NIR) variability of Sgr A* was shown to be consistent over two decades, unprecedented activity in 2019 challenges existing statistical models. We investigate the origin of this activity by re-calibrating and re-analyzing all of our Keck Observatory Sgr A* imaging observations from 2005-2022. We present light curves from 69 observation epochs using the NIRC2 imager at 2.12 $μ$m with laser guide star adaptive optics. These observations reveal that the mean luminosity of Sgr A* increased by a factor of $\sim$3 in 2019, and the 2019 light curves had higher variance than in all time periods we examined. We find that the 2020-2022 flux distribution is statistically consistent with the historical sample and model predictions, but with fewer bright measurements above 0.6 mJy at the $\sim$2$σ$ level. Since 2019, we have observed a maximum $K_s$ (2.2 $μ$m) flux of 0.9 mJy, compared to the highest pre-2019 flux of 2.0 mJy and highest 2019 flux of 5.6 mJy. Our results suggest that the 2019 activity was caused by a temporary accretion increase onto Sgr A*, possibly due to delayed accretion of tidally-stripped gas from the gaseous object G2 in 2014. We also examine faint Sgr A* fluxes over a long time baseline to search for a quasi-steady quiescent state. We find that Sgr A* displays flux variations over a factor of $\sim$500, with no evidence for a quiescent state in the NIR. △ Less

Submitted 18 August, 2023; originally announced August 2023.

Comments: Accepted to ApJ Letters, 26 pages, 15 figures

Constraining a companion of the galactic center black hole, Sgr A*

Authors: Clifford M. Will, Smadar Naoz, Aurélien Hees, Alexandria Tucker, Eric Zhang, Tuan Do, Andrea Ghez

Abstract: We use 23 years of astrometric and radial velocity data on the orbit of the star S0-2 to constrain a hypothetical intermediate-mass black hole orbiting the massive black hole Sgr A* at the Galactic center. The data place upper limits on variations of the orientation of the stellar orbit (inclination, nodal angle, and pericenter) at levels between 0.02 and 0.07 degrees per year. We use a combinatio… ▽ More We use 23 years of astrometric and radial velocity data on the orbit of the star S0-2 to constrain a hypothetical intermediate-mass black hole orbiting the massive black hole Sgr A* at the Galactic center. The data place upper limits on variations of the orientation of the stellar orbit (inclination, nodal angle, and pericenter) at levels between 0.02 and 0.07 degrees per year. We use a combination of analytic estimates and full numerical integrations of the orbit of S0-2 in the presence of a black-hole binary. For a companion IMBH whose semi-major axis $a_c$ is larger than that of S0-2 (1020 a.u.), we find that in the region between 1000 and 4000 a.u., a companion black hole with mass $m_c$ between $10^3$ and $10^5 M_\odot$ is excluded, with a boundary behaving as $a_c \sim m_c^{1/3}$. For a companion with $a_c < 1020$ a.u., we find that a black hole with mass between $10^3$ and $10^5 \, M_\odot$ is again excluded, with a boundary behaving as $a_c \sim m_c^{-1/2}$. These bounds arise from quadrupolar perturbations of the orbit of S0-2. However, significantly stronger bounds on the mass of an inner companion arise from the fact that the location of S0-2 is measured relative to the bright emission of Sgr A*. As a consequence, that separation is perturbed by the ``wobble'' of Sgr A* about the center of mass between it and the companion, leading to ``apparent'' perturbations of S0-2's orbit that also include a dipole component. The result is a set of bounds as small as $400 \, M_\odot$ at 200 a.u.; the numerical simulations suggest a bound from these effects varying as $a_c \sim m_c^{-1}$. We compare and contrast our results with those from a recent analysis by the GRAVITY collaboration. △ Less

Submitted 11 December, 2023; v1 submitted 31 July, 2023; originally announced July 2023.

Comments: 9 pages, 4 figures; minor revisions to agree with published version

Journal ref: Astrophys. J. 959, 58 (2023)

Improving Photometric Redshift Estimation for Cosmology with LSST using Bayesian Neural Networks

Authors: Evan Jones, Tuan Do, Bernie Boscoe, Jack Singal, Yujie Wan, Zooey Nguyen

Abstract: We present results exploring the role that probabilistic deep learning models can play in cosmology from large-scale astronomical surveys through photometric redshift (photo-z) estimation. Photo-z uncertainty estimates are critical for the science goals of upcoming large-scale surveys such as LSST, however common machine learning methods typically provide only point estimates and lack uncertaintie… ▽ More We present results exploring the role that probabilistic deep learning models can play in cosmology from large-scale astronomical surveys through photometric redshift (photo-z) estimation. Photo-z uncertainty estimates are critical for the science goals of upcoming large-scale surveys such as LSST, however common machine learning methods typically provide only point estimates and lack uncertainties on predictions. We turn to Bayesian neural networks (BNNs) as a promising way to provide accurate predictions of redshift values with uncertainty estimates. We have compiled a galaxy data set from the Hyper Suprime-Cam Survey with grizy photometry, which is designed to be a smaller scale version of large surveys like LSST. We use this data set to investigate the performance of a neural network (NN) and a probabilistic BNN for photo-z estimation and evaluate their performance with respect to LSST photo-z science requirements. We also examine the utility of photo-z uncertainties as a means to reduce catastrophic outlier estimates. The BNN outputs the estimate in the form of a Gaussian probability distribution. We use the mean and standard deviation as the redshift estimate and uncertainty. We find that the BNN can produce accurate uncertainties. Using a coverage test, we find excellent agreement with expectation -- 67.2$\%$ of galaxies between $0 < 2.5$ have 1-$σ$ uncertainties that cover the spectroscopic value. We also include a comparison to alternative machine learning models using the same data. We find the BNN meets two out of three of the LSST photo-z science requirements in the range $0 < z < 2.5$. △ Less

Submitted 18 March, 2024; v1 submitted 22 June, 2023; originally announced June 2023.

Comments: 17 pages, 23 figures, 4 tables

MSC Class: 1

The Galactic Center with Roman

Authors: Sean K. Terry, Matthew W. Hosek Jr., Jessica R. Lu, Casey Lam, Natasha Abrams, Arash Bahramian, Richard Barry, Jean-Phillipe Beaulieu, Aparna Bhattacharya, Devin Chu, Anna Ciurlo, Will Clarkson, Tuan Do, Kareem El-Badry, Ryan Felton, Matthew Freeman, Abhimat Gautam, Andrea Ghez, Daniel Huber, Jason Hunt, Macy Huston, Tharindu Jayasinghe, Naoki Koshimoto, Madeline Lucey, Florian Peißker , et al. (9 additional authors not shown)

Abstract: We advocate for a Galactic center (GC) field to be added to the Galactic Bulge Time Domain Survey (GBTDS). The new field would yield high-cadence photometric and astrometric measurements of an unprecedented ${\sim}$3.3 million stars toward the GC. This would enable a wide range of science cases, such as finding star-compact object binaries that may ultimately merge as LISA-detectable gravitational… ▽ More We advocate for a Galactic center (GC) field to be added to the Galactic Bulge Time Domain Survey (GBTDS). The new field would yield high-cadence photometric and astrometric measurements of an unprecedented ${\sim}$3.3 million stars toward the GC. This would enable a wide range of science cases, such as finding star-compact object binaries that may ultimately merge as LISA-detectable gravitational wave sources, constraining the mass function of stars and compact objects in different environments, detecting populations of microlensing and transiting exoplanets, studying stellar flares and variability in young and old stars, and monitoring accretion onto the central supermassive black hole. In addition, high-precision proper motions and parallaxes would open a new window into the large-scale dynamics of stellar populations at the GC, yielding insights into the formation and evolution of galactic nuclei and their co-evolution with the growth of the supermassive black hole. We discuss the possible trade-offs between the notional GBTDS and the addition of a GC field with either an optimal or minimal cadence. Ultimately, the addition of a GC field to the GBTDS would dramatically increase the science return of Roman and provide a legacy dataset to study the mid-plane and innermost regions of our Galaxy. △ Less

Submitted 21 June, 2023; originally announced June 2023.

Comments: 19 pages, 3 figures. Submitted to the NASA Roman Core Community Surveys White Paper Call

The Orbital Eccentricities of Directly Imaged Companions Using Observable-Based Priors: Implications for Population-level Distributions

Authors: Clarissa R. Do Ó, Kelly K. O'Neil, Quinn M. Konopacky, Tuan Do, Gregory D. Martinez, Jean-Baptiste Ruffio, Andrea M. Ghez

Abstract: The eccentricity of a sub-stellar companion is an important tracer of its formation history. Directly imaged companions often present poorly constrained eccentricities. A recently developed prior framework for orbit fitting called ''observable-based priors'' has the advantage of improving biases in derived orbit parameters for objects with minimal phase coverage, which is the case for the majority… ▽ More The eccentricity of a sub-stellar companion is an important tracer of its formation history. Directly imaged companions often present poorly constrained eccentricities. A recently developed prior framework for orbit fitting called ''observable-based priors'' has the advantage of improving biases in derived orbit parameters for objects with minimal phase coverage, which is the case for the majority of directly imaged companions. We use observable-based priors to fit the orbits of 21 exoplanets and brown dwarfs in an effort to obtain the eccentricity distributions with minimized biases. We present the objects' individual posteriors compared to their previously derived distributions, showing in many cases a shift toward lower eccentricities. We analyze the companions' eccentricity distribution at a population level, and compare this to the distributions obtained with the traditional uniform priors. We fit a Beta distribution to our posteriors using observable-based priors, obtaining shape parameters $α= 1.09^{+0.30}_{-0.22}$ and $β= 1.42^{+0.33}_{-0.25}$. This represents an approximately flat distribution of eccentricities. The derived $α$ and $β$ parameters are consistent with the values obtained using uniform priors, though uniform priors lead to a tail at high eccentricities. We find that separating the population into high and low mass companions yields different distributions depending on the classification of intermediate mass objects. We also determine via simulation that the minimal orbit coverage needed to give meaningful posteriors under the assumptions made for directly imaged planets is $\approx$ 15% of the inferred period of the orbit. △ Less

Submitted 18 July, 2023; v1 submitted 6 June, 2023; originally announced June 2023.

Comments: Published in The Astronomical Journal, Volume 166, Issue 2, id.48, 22 pp. Update includes updated references and ORCID profiles

Anisotropic power-law inflation for a generalized model of two scalar and two vector fields

Authors: Tuan Q. Do, W. F. Kao

Abstract: Cosmological implication of a generalized model of two scalar and two vector fields, in which both scalar fields are non-minimally coupled to each vector field, is studied in this paper. In particular, we will seek a set of new anisotropic power-law inflationary solutions to this model. Additionally, the stability of the obtained solutions will be examined by using the dynamical system approach. A… ▽ More Cosmological implication of a generalized model of two scalar and two vector fields, in which both scalar fields are non-minimally coupled to each vector field, is studied in this paper. In particular, we will seek a set of new anisotropic power-law inflationary solutions to this model. Additionally, the stability of the obtained solutions will be examined by using the dynamical system approach. As a result, we will show that this set of solutions turns out to be stable and attractive during the inflationary phase as expected due to the existence of the unusual couplings between two scalar and two vector fields. Notably, we will point out that the existence of phantom field will lead to an instability of the corresponding anisotropic power-law inflation. △ Less

Submitted 4 February, 2024; v1 submitted 18 April, 2023; originally announced April 2023.

Comments: 26 pages, 6 figures. Final version with some discussions and relevant references added; typos fixed. All calculations remain unchanged. Accepted for publication in Physica Scripta. Comments are welcome

Journal ref: Phys. Scr. 99, 015002 (2024)

The formation history of our Galaxy's nuclear stellar disc constrained from HST observations of the Quintuplet field

Authors: R. Schoedel, F. Nogueras-Lara, M. Hosek Jr, T. Do, J. Lu, A. Martínez Arranz, A. Ghez, R. M. Rich, A. Gardini, E. Gallego-Cano, M. Cano-González, A. T. Gallego-Calvente

Abstract: Until recently it was thought that the nuclear stellar disc at the centre of our Galaxy was formed via quasi-continuous star formation over billions of years. However, an analysis of GALACTICNUCLEUS survey data indicates that >80% of the mass of the stellar disc formed at least 8 Gyr ago and about 5% roughly 1 Gyr ago. Our aim is to derive new constraints on the formation history of the nuclear st… ▽ More Until recently it was thought that the nuclear stellar disc at the centre of our Galaxy was formed via quasi-continuous star formation over billions of years. However, an analysis of GALACTICNUCLEUS survey data indicates that >80% of the mass of the stellar disc formed at least 8 Gyr ago and about 5% roughly 1 Gyr ago. Our aim is to derive new constraints on the formation history of the nuclear stellar disc. We analysed a catalogue of HST/WFC3-IR observations of the Quintuplet cluster field. From this catalogue, we selected about 24000 field stars that probably belong to the nuclear stellar disc. We used red clump giants to deredden the sample and fit the resulting F153M luminosity function with a linear combination of theoretical luminosity functions created from different stellar evolutionary models. We find that >70% of the stellar population in the nuclear disc probably formed more than 10 Gyr ago, while ~15% formed in an event (or series of events) ~1Gyr ago. Up to 10% of the stars appear to have formed in the past tens to hundreds of Myr. These results do not change significantly for reasonable variations in the assumed mean metallicity, sample selection, reddening correction, or stellar evolutionary models. We confirm previous work that changed the formation paradigm for stars in the Galactic Centre. The nuclear stellar disc is indeed a very old structure. There seems to have been little star formation activity between its formation and about 1 Gyr ago. △ Less

Submitted 4 April, 2023; originally announced April 2023.

Comments: Accepted for publication in A&A

Journal ref: A&A 672, L8 (2023)

Stability investigations of isotropic and anisotropic exponential inflation in the Starobinsky-Bel-Robinson gravity

Authors: Tuan Q. Do, Duy H. Nguyen, Tuyen M. Pham

Abstract: In this paper, we would like to examine whether a novel Starobinsky-Bel-Robinson gravity model admits stable exponential inflationary solutions with or without spatial anisotropies. As a result, we are able to derive an exact de Sitter inflationary to this Starobinsky-Bel-Robinson model. Furthermore, we observe that an exact Bianchi type I inflationary solution does not exist in the Starobinsky-Be… ▽ More In this paper, we would like to examine whether a novel Starobinsky-Bel-Robinson gravity model admits stable exponential inflationary solutions with or without spatial anisotropies. As a result, we are able to derive an exact de Sitter inflationary to this Starobinsky-Bel-Robinson model. Furthermore, we observe that an exact Bianchi type I inflationary solution does not exist in the Starobinsky-Bel-Robinson model. However, we find that a modified Starobinsky-Bel-Robinson model, in which the sign of coefficient of $R^2$ term is flipped from positive to negative, can admit the corresponding Bianchi type I inflationary solution. Unfortunately, stability analysis using the dynamical system approach indicates that both of these inflationary solutions turn out to be unstable. Interestingly, we show that a stable de Sitter inflationary solution can be obtained in the modified Starobinsky-Bel-Robinson gravity. △ Less

Submitted 29 November, 2023; v1 submitted 30 March, 2023; originally announced March 2023.

Comments: 32 pages, 3 figures. V3 with many discussions, calculations, and relevant references added. Accepted for publication by International Journal of Modern Physics D. Comments are welcome

Journal ref: International Journal of Modern Physics D 32, 2350087 (2023)

Evidence of a Decreased Binary Fraction for Massive Stars Within 20 Milliparsecs of the Supermassive Black Hole at the Galactic Center

Authors: Devin S. Chu, Tuan Do, Andrea Ghez, Abhimat K. Gautam, Anna Ciurlo, Kelly Kosmo O'Neil, Matthew W. Hosek Jr., Aurélien Hees, Smadar Naoz, Shoko Sakai, Jessica R. Lu, Zhuo Chen, Rory O. Bentley, Eric E. Becklin, Keith Matthews

Abstract: We present the results of the first systematic search for spectroscopic binaries within the central 2 x 3 arcsec$^2$ around the supermassive black hole at the center of the Milky Way galaxy. This survey is based primarily on over a decade of adaptive optics-fed integral-field spectroscopy (R$\sim$4000), obtained as part of the Galactic Center Orbits Initiative at Keck Observatory, and has a limiti… ▽ More We present the results of the first systematic search for spectroscopic binaries within the central 2 x 3 arcsec$^2$ around the supermassive black hole at the center of the Milky Way galaxy. This survey is based primarily on over a decade of adaptive optics-fed integral-field spectroscopy (R$\sim$4000), obtained as part of the Galactic Center Orbits Initiative at Keck Observatory, and has a limiting $K$'-band magnitude of 15.8, which is at least 4 magnitudes deeper than previous spectroscopic searches for binaries at larger radii within the central nuclear star cluster. From this primary dataset, over 600 new radial velocities are extracted and reported, increasing by a factor of 3 the number of such measurements. We find no significant periodic signals in our sample of 28 stars, of which 16 are massive, young (main-sequence B) stars and 12 are low-mass, old (M and K giant) stars. Using Monte Carlo simulations, we derive upper limits on the intrinsic binary star fraction for the young star population at 47% (at 95% confidence) located $\sim$20 mpc from the black hole. The young star binary fraction is significantly lower than that observed in the field (70%). This result is consistent with a scenario in which the central supermassive black hole drives nearby stellar binaries to merge or be disrupted and may have important implications for the production of gravitational waves and hypervelocity stars. △ Less

Submitted 29 March, 2023; originally announced March 2023.

Comments: Accepted by ApJ. 22 pages, 14 figures

Stellar Populations in the Central 0.5 pc of Our Galaxy III: The Dynamical Sub-structures

Authors: Siyao Jia, Ningyuan Xu, Jessica R. Lu, D. S Chu, K. Kosmo O'Neil, W. B. Drechsler, M. W. Hosek Jr., S. Sakai, T. Do, A. Ciurlo, A. K. Gautam, A. M. Ghez, E. Becklin, M. R. Morris, R. O. Bentley

Abstract: We measure the 3D kinematic structures of the young stars within the central 0.5 parsec of our Galactic Center using the 10 m telescopes of the W.~M.~Keck Observatory over a time span of 25 years. Using high-precision measurements of positions on the sky, and proper motions and radial velocities from new observations and the literature, we constrain the orbital parameters for each young star. Our… ▽ More We measure the 3D kinematic structures of the young stars within the central 0.5 parsec of our Galactic Center using the 10 m telescopes of the W.~M.~Keck Observatory over a time span of 25 years. Using high-precision measurements of positions on the sky, and proper motions and radial velocities from new observations and the literature, we constrain the orbital parameters for each young star. Our results show two statistically significant sub-structures: a clockwise stellar disk with 18 candidate stars, as has been proposed before, but with an improved disk membership; a second, almost edge-on plane of 10 candidate stars oriented East-West on the sky that includes at least one IRS 13 star. We estimate the eccentricity distribution of each sub-structure and find that the clockwise disk has <$e$> = 0.39 and the edge-on plane has <$e$> = 0.68. We also perform simulations of each disk/plane with incompleteness and spatially-variable extinction to search for asymmetry. Our results show that the clockwise stellar disk is consistent with a uniform azimuthal distribution within the disk. The edge-on plane has an asymmetry that cannot be explained by variable extinction or incompleteness in the field. The orientation, asymmetric stellar distribution, and high eccentricity of the edge-on plane members suggest that this structure may be a stream associated with the IRS 13 group. The complex dynamical structure of the young nuclear cluster indicates that the star formation process involved complex gas structures and dynamics and is inconsistent with a single massive gaseous disk. △ Less

Submitted 3 February, 2023; originally announced February 2023.

Comments: 41 pages, 26 figures, 13 tables, 2 appendices. Accepted for publication in ApJ

The Swansong of the Galactic Center Source X7: An Extreme Example of Tidal Evolution near the Supermassive Black Hole

Authors: Anna Ciurlo, Randall D. Campbell, Mark R. Morris, Tuan Do, Andrea M. Ghez, Eric E. Becklin, Rory O. Bentley, Devin S. Chu, Abhimat K. Gautam, Yash A. Gursahani, Aurelien Hees, Kelly Kosmo O'Neil, Jessica R. Lu, Gregory D. Martinez, Smadar Naoz, Shoko Sakai, Rainer Schoedel

Abstract: We present two decades of new high-angular-resolution near-infrared data from the W. M. Keck Observatory that reveal extreme evolution in X7, an elongated dust and gas feature, presently located half an arcsecond from the Galactic Center supermassive black hole. With both spectro-imaging observations of Br-γ line-emission and Lp (3.8 μm) imaging data, we provide the first estimate of its orbital p… ▽ More We present two decades of new high-angular-resolution near-infrared data from the W. M. Keck Observatory that reveal extreme evolution in X7, an elongated dust and gas feature, presently located half an arcsecond from the Galactic Center supermassive black hole. With both spectro-imaging observations of Br-γ line-emission and Lp (3.8 μm) imaging data, we provide the first estimate of its orbital parameters and quantitative characterization of the evolution of its morphology and mass. We find that the leading edge of X7 appears to be on a mildly eccentric (e~0.3), relatively short-period (170 years) orbit and is headed towards periapse passage, estimated to occur in ~2036. Furthermore, our kinematic measurements rule out the earlier suggestion that X7 is associated with the stellar source S0-73 or with any other point source that has overlapped with X7 during our monitoring period. Over the course of our observations, X7 has (1) become more elongated, with a current length-to-width ratio of 9, (2) maintained a very consistent long-axis orientation (position angle of 50 deg), (3) inverted its radial velocity differential from tip to tail from -50 to +80 km/sec, and (4) sustained its total brightness (12.8 Lp magnitudes at the leading edge) and color temperature (425 K), which suggest a constant mass of ~50 MEarth. We present a simple model showing that these results are compatible with the expected effect of tidal forces exerted on it by the central black hole and we propose that X7 is the gas and dust recently ejected from a grazing collision in a binary system. △ Less

Submitted 16 January, 2023; originally announced January 2023.

Comments: accepted by ApJ, 20 pages, 16 figures

The Star Formation History of the Milky Way's Nuclear Star Cluster

Authors: Zhuo Chen, Tuan Do, Andrea M. Ghez, Matthew Hosek Jr., Anja Feldmeier-Krause, Devin Chu, Rory Bentley, Jessica R. Lu, Mark R. Morris

Abstract: We report the first star formation history study of the Milky Way's nuclear star cluster (NSC) that includes observational constraints from a large sample of stellar metallicity measurements. These metallicity measurements were obtained from recent surveys from Gemini and VLT of 770 late-type stars within the central 1.5 pc. These metallicity measurements, along with photometry and spectroscopical… ▽ More We report the first star formation history study of the Milky Way's nuclear star cluster (NSC) that includes observational constraints from a large sample of stellar metallicity measurements. These metallicity measurements were obtained from recent surveys from Gemini and VLT of 770 late-type stars within the central 1.5 pc. These metallicity measurements, along with photometry and spectroscopically derived temperatures, are forward modeled with a Bayesian inference approach. Including metallicity measurements improves the overall fit quality, as the low-temperature red giants that were previously difficult to constrain are now accounted for, and the best fit favors a two-component model. The dominant component contains 93%$\pm$3% of the mass, is metal-rich ($\overline{[M/H]}\sim$0.45), and has an age of 5$^{+3}_{-2}$ Gyr, which is $\sim$3 Gyr younger than earlier studies with fixed (solar) metallicity; this younger age challenges co-evolutionary models in which the NSC and supermassive black holes formed simultaneously at early times. The minor population component has low metallicity ($\overline{[M/H]}\sim$ -1.1) and contains $\sim$7% of the stellar mass. The age of the minor component is uncertain (0.1 - 5 Gyr old). Using the estimated parameters, we infer the following NSC stellar remnant population (with $\sim$18% uncertainty): 1.5$\times$10$^5$ neutron stars, 2.5$\times$10$^5$ stellar mass black holes (BHs) and 2.2$\times$10$^4$ BH-BH binaries. These predictions result in 2-4 times fewer neutron stars compared to earlier predictions that assume solar metallicity, introducing a possible new path to understand the so-called "missing pulsar problem". Finally, we present updated predictions for the BH-BH merger rates (0.01-3 Gpc$^{-3}$yr$^{-1}$). △ Less

Submitted 2 December, 2022; originally announced December 2022.

Comments: 37 pages, 23 figures, Accepted to ApJ

Elements of effective machine learning datasets in astronomy

Authors: Bernie Boscoe, Tuan Do, Evan Jones, Yunqi Li, Kevin Alfaro, Christy Ma

Abstract: In this work, we identify elements of effective machine learning datasets in astronomy and present suggestions for their design and creation. Machine learning has become an increasingly important tool for analyzing and understanding the large-scale flood of data in astronomy. To take advantage of these tools, datasets are required for training and testing. However, building machine learning datase… ▽ More In this work, we identify elements of effective machine learning datasets in astronomy and present suggestions for their design and creation. Machine learning has become an increasingly important tool for analyzing and understanding the large-scale flood of data in astronomy. To take advantage of these tools, datasets are required for training and testing. However, building machine learning datasets for astronomy can be challenging. Astronomical data is collected from instruments built to explore science questions in a traditional fashion rather than to conduct machine learning. Thus, it is often the case that raw data, or even downstream processed data is not in a form amenable to machine learning. We explore the construction of machine learning datasets and we ask: what elements define effective machine learning datasets? We define effective machine learning datasets in astronomy to be formed with well-defined data points, structure, and metadata. We discuss why these elements are important for astronomical applications and ways to put them in practice. We posit that these qualities not only make the data suitable for machine learning, they also help to foster usable, reusable, and replicable science practices. △ Less

Submitted 29 November, 2022; v1 submitted 25 November, 2022; originally announced November 2022.

Comments: 5 pages, 1 figure, accepted to the peer-reviewed NeurIPS Machine Learning in the Physical Sciences Workshop, 2022

Anisotropic constant-roll k-inflation model

Authors: Duy H. Nguyen, Tuyen M. Pham, Thien D. Le, Tuan Q. Do

Abstract: In this paper, we would like to figure out whether a k-inflation model admits the Bianchi type I metric as its inflationary solution under a constant-roll condition in the presence of the supergravity motivated coupling between scalar and vector fields, $f^2(φ)F_{μν}F^{μν}$. As a result, some novel anisotropic inflationary solutions are shown to appear along with a power-law one in this scenario.… ▽ More In this paper, we would like to figure out whether a k-inflation model admits the Bianchi type I metric as its inflationary solution under a constant-roll condition in the presence of the supergravity motivated coupling between scalar and vector fields, $f^2(φ)F_{μν}F^{μν}$. As a result, some novel anisotropic inflationary solutions are shown to appear along with a power-law one in this scenario. Furthermore, these solutions are numerically confirmed to be attractive, in contrast to the prediction of the Hawking's cosmic no-hair conjecture. △ Less

Submitted 15 November, 2022; originally announced November 2022.

Comments: 15 pages, 11 figures. To appear in Communications in Physics (2023), a well-established physics journal of Vietnam. Comments are welcome

Journal ref: Communications in Physics 33, 15 (2023)

AIROPA II: Modeling Instrumental Aberrations for Off-Axis Point Spread Functions in Adaptive Optics

Authors: Anna Ciurlo, Paolo Turri, Gunther Witzel, Jessica R. Lu, Tuan Do, Breann N. Sitarski, Michael P. Fitzgerald, Andrea M. Ghez, Carlos Alvarez, Sean K. Terry, Greg Doppmann, James E. Lyke, Sam Ragland, Randall Campbell, Keith Matthews

Abstract: Images obtained with single-conjugate adaptive optics (AO) show spatial variation of the point spread function (PSF) due to both atmospheric anisoplanatism and instrumental aberrations. The poor knowledge of the PSF across the field of view strongly impacts the ability to take full advantage of AO capabilities. The AIROPA project aims to model these PSF variations for the NIRC2 imager at the Keck… ▽ More Images obtained with single-conjugate adaptive optics (AO) show spatial variation of the point spread function (PSF) due to both atmospheric anisoplanatism and instrumental aberrations. The poor knowledge of the PSF across the field of view strongly impacts the ability to take full advantage of AO capabilities. The AIROPA project aims to model these PSF variations for the NIRC2 imager at the Keck Observatory. Here, we present the characterization of the instrumental phase aberrations over the entire NIRC2 field of view and we present a new metric for quantifying the quality of the calibration, the fraction of variance unexplained (FVU). We used phase diversity measurements obtained on an artificial light source to characterize the variation of the aberrations across the field of view and their evolution with time. We find that there is a daily variation of the wavefront error (RMS of the residuals is 94~nm) common to the whole detector, but the differential aberrations across the field of view are very stable (RMS of the residuals between different epochs is 59~nm). This means that instrumental calibrations need to be monitored often only at the center of the detector, and the much more time-consuming variations across the field of view can be characterized less frequently (most likely when hardware upgrades happen). Furthermore, we tested AIROPA's instrumental model through real data of the fiber images on the detector. We find that modeling the PSF variations across the field of view improves the FVU metric by 60\% and reduces the detection of fake sources by 70\%. △ Less

Submitted 19 October, 2022; originally announced October 2022.

Journal ref: Journal of Astronomical Telescopes, Instruments, and Systems, Vol. 8, Issue 3, 038007 (September 2022)

AIROPA IV: Validating Point Spread Function Reconstruction on Various Science Cases

Authors: Sean K. Terry, Jessica R. Lu, Paolo Turri, Anna Ciurlo, Abhimat Gautam, Tuan Do, Michael P. Fitzgerald, Andrea Ghez, Matthew Hosek Jr., Gunther Witzel

Abstract: We present an analysis of six independent on-sky datasets taken with the Keck-II/NIRC2 instrument. Using the off-axis point spread function (PSF) reconstruction software AIROPA, we extract stellar astrometry, photometry, and other fitting metrics in order to characterize the performance of this package. We test the effectiveness of AIROPA to reconstruct the PSF across the field of view in varying… ▽ More We present an analysis of six independent on-sky datasets taken with the Keck-II/NIRC2 instrument. Using the off-axis point spread function (PSF) reconstruction software AIROPA, we extract stellar astrometry, photometry, and other fitting metrics in order to characterize the performance of this package. We test the effectiveness of AIROPA to reconstruct the PSF across the field of view in varying atmospheric conditions, number and location of PSF reference stars, stellar crowding and telescope position angle (PA). We compare the astrometric precision and fitting residuals between a static PSF model and a spatially varying PSF model that incorporates instrumental aberrations and atmospheric turbulence during exposures. Most of the fitting residuals we measure show little to no improvement in the variable-PSF mode over the single-PSF mode. For one of the data sets, we find photometric performance is significantly improved (by ${\sim}10\times$) by measuring the trend seen in photometry as a function of off-axis location. For nearly all other metrics we find comparable astrometric and photometric precision across both PSF modes, with a ${\sim}13$% smaller astrometric uncertainty in variable-PSF mode in the best case. We largely confirm that the spatially variable PSF does not significantly improve the astrometric and other PSF fitting residuals over the static PSF for on-sky observations. We attribute this to unaccounted instrumental aberrations that are not characterized through afternoon adaptive optics (AO) bench calibrations. △ Less

Submitted 12 September, 2022; originally announced September 2022.

Comments: 23 pages, 13 figures, submitted to JATIS

Measuring the Orbits of the Arches and Quintuplet Clusters using HST and Gaia: Exploring Scenarios for Star Formation Near the Galactic Center

Authors: Matthew W. Hosek Jr., Tuan Do, Jessica R. Lu, Mark R. Morris, Andrea M. Ghez, Gregory D. Martinez, Jay Anderson

Abstract: We present new absolute proper motion measurements for the Arches and Quintuplet clusters, two young massive star clusters near the Galactic center. Using multi-epoch HST observations, we construct proper motion catalogs for the Arches ($\sim$35,000 stars) and Quintuplet ($\sim$40,000 stars) fields in ICRF coordinates established using stars in common with the Gaia EDR3 catalog. The bulk proper mo… ▽ More We present new absolute proper motion measurements for the Arches and Quintuplet clusters, two young massive star clusters near the Galactic center. Using multi-epoch HST observations, we construct proper motion catalogs for the Arches ($\sim$35,000 stars) and Quintuplet ($\sim$40,000 stars) fields in ICRF coordinates established using stars in common with the Gaia EDR3 catalog. The bulk proper motions of the clusters are measured to be ($μ_{α*}$, $μ_δ$) = (-0.80 $\pm$ 0.032, -1.89 $\pm$ 0.021) mas/yr for the Arches and ($μ_{α*}$, $μ_δ$) = (-0.96 $\pm$ 0.032, -2.29 $\pm$ 0.023) mas/yr for the Quintuplet, achieving $\sim$5x higher precision than past measurements. We place the first constraints on the properties of the cluster orbits that incorporate the uncertainty in their current line-of-sight distances. The clusters will not approach closer than $\sim$25 pc to SgrA*, making it unlikely that they will inspiral into the Nuclear Star Cluster within their lifetime. Further, the cluster orbits are not consistent with being circular; the average value of r$_{apo}$ / r$_{peri}$ is $\sim$1.9 (equivalent to eccentricity of $\sim$0.31) for both clusters. Lastly, we find that the clusters do not share a common orbit, challenging one proposed formation scenario in which the clusters formed from molecular clouds on the open stream orbit derived by Kruijssen et al. (2015). Meanwhile, our constraints on the birth location and velocity of the clusters offer mild support for a scenario in which the clusters formed via collisions between gas clouds on the x1 and x2 bar orbit families. △ Less

Submitted 17 August, 2022; originally announced August 2022.

Comments: Accepted for publication in ApJ. 38 pages, 25 figures. Proper motion catalogs included in ancillary materials

AIROPA III: Testing Simulated and On-Sky Data

Authors: Paolo Turri, Jessica R. Lu, Gunther Witzel, Anna Ciurlo, Tuan Do, Andrea M. Ghez, Matthew C. Britton, Sam Ragland, Sean K. Terry

Abstract: Adaptive optics images from the W. M. Keck Observatory have delivered numerous influential scientific results, including detection of multi-system asteroids, the supermassive black hole at the center of the Milky Way, and directly imaged exoplanets. Specifically, the precise and accurate astrometry these images yield was used to measure the mass of the supermassive black hole using orbits of the s… ▽ More Adaptive optics images from the W. M. Keck Observatory have delivered numerous influential scientific results, including detection of multi-system asteroids, the supermassive black hole at the center of the Milky Way, and directly imaged exoplanets. Specifically, the precise and accurate astrometry these images yield was used to measure the mass of the supermassive black hole using orbits of the surrounding star cluster. Despite these successes, one of the major obstacles to improved astrometric measurements is the spatial and temporal variability of the point-spread function delivered by the instruments. AIROPA is a software package for the astrometric and photometric analysis of adaptive optics images using point-spread function fitting together with the technique of point-spread function reconstruction. In adaptive optics point-spread function reconstruction, the knowledge of the instrument performance and of the atmospheric turbulence is used to predict the long-exposure point-spread function of an observation. In this paper we present the results of our tests using AIROPA on both simulated and on-sky images of the Galactic Center. We find that our method is very reliable in accounting for the static aberrations internal to the instrument, but it does not improve significantly the accuracy on sky, possibly due to uncalibrated telescope aberrations. △ Less

Submitted 1 July, 2022; originally announced July 2022.

Searching for a Hypervelocity White Dwarf Companion: A Proper Motion Survey of SN 1006

Authors: Joshua V. Shields, Wolfgang Kerzendorf, Matthew W. Hosek Jr., Ken J. Shen, Armin Rest, Tuan Do, Jessica R. Lu, Andrew G. Fullard, Giovanni Strampelli, Alfredo Zenteno

Abstract: Type Ia Supernovae (SNe Ia) are securely understood to come from the thermonuclear explosion of a white dwarf as a result of binary interaction, but the nature of that binary interaction and the secondary object is uncertain. Recently, a double white dwarf model known as the dynamically driven double-degenerate double-detonation (D6) model has become a promising explanation for these events. One r… ▽ More Type Ia Supernovae (SNe Ia) are securely understood to come from the thermonuclear explosion of a white dwarf as a result of binary interaction, but the nature of that binary interaction and the secondary object is uncertain. Recently, a double white dwarf model known as the dynamically driven double-degenerate double-detonation (D6) model has become a promising explanation for these events. One realization of this scenario predicts that the companion may survive the explosion and reside within the remnant as a fast moving ($V_{peculiar} >1000$ km s$^{-1}$), overluminous ($L > 0.1 L_\odot$) white dwarf. Recently, three objects which appear to have these unusual properties have been discovered in the Gaia survey. We obtained photometric observations of the SN Ia remnant SN 1006 with the Dark Energy Camera over four years to attempt to discover a similar star. We present a deep, high precision astrometric proper motion survey of the interior stellar population of the remnant. We rule out the existence of a high proper motion object consistent with our tested realization of the D6 scenario ($V_{transverse} > 600$ km s$^{-1}$ with $m_r < 21$ corresponding to an intrinsic luminosity of $L > 0.0176 L_\odot$). We conclude that such a star does not exist within the remnant, or is hidden from detection by either strong localized dust or the unlikely possibility of ejection from the binary system near parallel to the line of sight. △ Less

Submitted 8 June, 2022; originally announced June 2022.

Comments: 15 pages, 10 figures

Journal ref: https://iopscience.iop.org/article/10.3847/2041-8213/ac7950

Regularized 3D spectroscopy with CubeFit: method and application to the Galactic Center Circumnuclear disk

Authors: Thibaut Paumard, Anna Ciurlo, Mark R. Morris, Tuan Do, Andrea M. Ghez

Abstract: The Galactic Center black hole and the nuclear star cluster are surrounded by a clumpy ring of gas and dust (the circumnuclear disk, CND) that rotates about them at a standoff distance of ~1.5 pc. The mass and density of individual clumps in the CND are disputed. We seek to use H$_2$ to characterize the clump size distribution and to investigate the morphology and dynamics of the interface between… ▽ More The Galactic Center black hole and the nuclear star cluster are surrounded by a clumpy ring of gas and dust (the circumnuclear disk, CND) that rotates about them at a standoff distance of ~1.5 pc. The mass and density of individual clumps in the CND are disputed. We seek to use H$_2$ to characterize the clump size distribution and to investigate the morphology and dynamics of the interface between the ionized interior layer of the CND and the molecular reservoir lying further out (corresponding to the inner rim of the CND, illuminated in ultraviolet light by the central star cluster). We have observed two fields of approximately 20"x20" in the CND at near-infrared wavelengths with the OSIRIS spectro-imager at the Keck Observatory. These two fields, located at the approaching and receding nodes of the CND, best display this interface. Our data cover two H$_2$ lines as well as the Br$γ$ line (tracing H ii). We have developed the tool CubeFit, an original method to extract maps of continuous physical parameters (such as velocity field and velocity dispersion) from integral-field spectroscopy data, using regularization to largely preserve spatial resolution in regions of low signal-to-noise ratio. This original method enables us to isolate compact, bright features in the interstellar medium of the CND. Several clumps in the southwestern field assume the appearance of filaments, many of which are parallel to each other. We conclude that these clumps cannot be self-gravitating. △ Less

Submitted 25 April, 2022; originally announced April 2022.

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

Journal ref: A&A 664, A97 (2022)

Multi-wavelength Variability of Sagittarius A* in July 2019

Authors: H. Boyce, D. Haggard, G. Witzel, S. von Fellenberg, S. P. Willner, E. E. Becklin, T. Do, A. Eckart, G. G. Fazio, M. A. Gurwell, J. L. Hora, S. Markoff, M. R. Morris, J. Neilsen, M. Nowak, H. A. Smith, S. Zhang

Abstract: We report timing analysis of near-infrared (NIR), X-ray, and sub-millimeter (submm) data during a three-day coordinated campaign observing Sagittarius A*. Data were collected at 4.5 micron with the Spitzer Space Telescope, 2-8 keV with the Chandra X-ray Observatory, 3-70 keV with NuSTAR, 340 GHz with ALMA, and at 2.2 micron with the GRAVITY instrument on the Very Large Telescope Interferometer. Tw… ▽ More We report timing analysis of near-infrared (NIR), X-ray, and sub-millimeter (submm) data during a three-day coordinated campaign observing Sagittarius A*. Data were collected at 4.5 micron with the Spitzer Space Telescope, 2-8 keV with the Chandra X-ray Observatory, 3-70 keV with NuSTAR, 340 GHz with ALMA, and at 2.2 micron with the GRAVITY instrument on the Very Large Telescope Interferometer. Two dates show moderate variability with no significant lags between the submm and the infrared at 99% confidence. July 18 captured a moderately bright NIR flare (F_K ~ 15 mJy) simultaneous with an X-ray flare (F ~ 0.1 cts/s) that most likely preceded bright submm flux (F ~ 5.5 Jy) by about +34 (+14 -33) minutes at 99% confidence. The uncertainty in this lag is dominated by the fact that we did not observe the peak of the submm emission. A synchrotron source cooled through adiabatic expansion can describe a rise in the submm once the synchrotron-self-Compton NIR and X-ray peaks have faded. This model predicts high GHz and THz fluxes at the time of the NIR/X-ray peak and electron densities well above those implied from average accretion rates for Sgr A*. However, the higher electron density postulated in this scenario would be in agreement with the idea that 2019 was an extraordinary epoch with a heightened accretion rate. Since the NIR and X-ray peaks can also be fit by a non-thermal synchrotron source with lower electron densities, we cannot rule out an unrelated chance coincidence of this bright submm flare with the NIR/X-ray emission. △ Less

Submitted 24 March, 2022; originally announced March 2022.

Comments: Accepted for publication in The Astrophysical Journal

Photometric Redshifts for Cosmology: Improving Accuracy and Uncertainty Estimates Using Bayesian Neural Networks

Authors: Evan Jones, Tuan Do, Bernie Boscoe, Yujie Wan, Zooey Nguyen, Jack Singal

Abstract: We present results exploring the role that probabilistic deep learning models can play in cosmology from large scale astronomical surveys through estimating the distances to galaxies (redshifts) from photometry. Due to the massive scale of data coming from these new and upcoming sky surveys, machine learning techniques using galaxy photometry are increasingly adopted to predict galactic redshifts… ▽ More We present results exploring the role that probabilistic deep learning models can play in cosmology from large scale astronomical surveys through estimating the distances to galaxies (redshifts) from photometry. Due to the massive scale of data coming from these new and upcoming sky surveys, machine learning techniques using galaxy photometry are increasingly adopted to predict galactic redshifts which are important for inferring cosmological parameters such as the nature of dark energy. Associated uncertainty estimates are also critical measurements, however, common machine learning methods typically provide only point estimates and lack uncertainty information as outputs. We turn to Bayesian neural networks (BNNs) as a promising way to provide accurate predictions of redshift values. We have compiled a new galaxy training dataset from the Hyper Suprime-Cam Survey, designed to mimic large surveys, but over a smaller portion of the sky. We evaluate the performance and accuracy of photometric redshift (photo-z) predictions from photometry using machine learning, astronomical and probabilistic metrics. We find that while the Bayesian neural network did not perform as well as non-Bayesian neural networks if evaluated solely by point estimate photo-z values, BNNs can provide uncertainty estimates that are necessary for cosmology △ Less

Submitted 14 February, 2022; originally announced February 2022.

Comments: 5 pages, 3 figures. Machine Learning and the Physical Sciences, NeurIPS 2021

Kinematics and Feedback in H II regions in the Dwarf Starburst Galaxy IC 10

Authors: Maren Cosens, Shelley A. Wright, Norman Murray, Lee Armus, Karin Sandstrom, Tuan Do, Kirsten Larson, Gregory Martinez, Sanchit Sabhlok, Andrey Vayner, James Wiley

Abstract: We present a survey of the central region of the nearest starburst galaxy, IC 10, using the W. M. Keck Observatory Keck Cosmic Web Imager (KCWI) at high spectral and spatial resolution. We map the central starburst of IC 10 to sample the kinematic and ionization properties of the individual star-forming regions. Using the low spectral resolution mode of KCWI we map the oxygen abundance and with th… ▽ More We present a survey of the central region of the nearest starburst galaxy, IC 10, using the W. M. Keck Observatory Keck Cosmic Web Imager (KCWI) at high spectral and spatial resolution. We map the central starburst of IC 10 to sample the kinematic and ionization properties of the individual star-forming regions. Using the low spectral resolution mode of KCWI we map the oxygen abundance and with the high spectral resolution mode we identify 46 individual H II regions. These H II regions have an average radius of 4.0 pc, star formation rate $\sim1.3\times10^{-4}$ M$_\odot$ yr$^{-1}$, and velocity dispersion $\sim$16 km s$^{-1}$. None of the H II regions appear to be virialized ($\rm α_{vir}>>1$), and, on average, they show evidence of ongoing expansion. IC 10's H II regions are offset from the star forming region size-luminosity scaling relationships, as well as Larson's Law that relates size and velocity dispersion. We investigate the balance of inward and outward pressure, $\rm P_{in}$ and $\rm P_{out}$, finding $\rm P_{out}>P_{in}$ in 89% of H II regions, indicating feedback driven expansion even in these low mass H II regions. We find warm gas pressure ($\rm P_{gas}$) provides the dominant contribution to the outward pressure ($\rm P_{out}$). This counteracts the inward pressure which is dominated by turbulence in the surrounding gas rather than self-gravity. Five H II regions show evidence of outflows which are most likely supported by either stellar winds (2 regions) or champagne flows (3 regions). These observations provide new insights into the state of the star-forming regions in IC 10 and negative feedback from low mass clusters. △ Less

Submitted 8 February, 2022; originally announced February 2022.

Comments: 57 pages, 35 figures, 9 tables, accepted for publication in The Astrophysical Journal

Machine Learning Classification to Identify Catastrophic Outlier Photometric Redshift Estimates

Authors: J. Singal, G. Silverman, E. Jones, T. Do, B. Boscoe, Y. Wan

Abstract: We present results of using a basic binary classification neural network model to identify likely catastrophic outlier photometric redshift estimates of individual galaxies, based only on the galaxies' measured photometric band magnitude values. We find that a simple implementation of this classification can identify a significant fraction of galaxies with catastrophic outlier photometric redshift… ▽ More We present results of using a basic binary classification neural network model to identify likely catastrophic outlier photometric redshift estimates of individual galaxies, based only on the galaxies' measured photometric band magnitude values. We find that a simple implementation of this classification can identify a significant fraction of galaxies with catastrophic outlier photometric redshift estimates while falsely categorizing only a much smaller fraction of non-outliers. These methods have the potential to reduce the errors introduced into science analyses by catastrophic outlier photometric redshift estimates. △ Less

Submitted 23 February, 2022; v1 submitted 14 December, 2021; originally announced December 2021.

Comments: 7 pages, 2 figure blocks. Updated to ApJ accepted version

Measuring the alpha-abundance of subsolar-metallicity stars in the Milky Way's central half-parsec: testing globular cluster and dwarf galaxy infall scenarios

Authors: Rory O. Bentley, Tuan Do, Wolfgang Kerzendorf, Devin S. Chu, Zhuo Chen, Quinn Konopacky, Andrea Ghez

Abstract: While the Milky Way Nuclear star cluster has been studied extensively, how it formed is uncertain. Studies have shown it contains a solar and supersolar metallicity population that may have formed in-situ, along with a subsolar metallicity population that may have formed via mergers of globular clusters and dwarf galaxies. Stellar abundance measurements are critical to differentiate between format… ▽ More While the Milky Way Nuclear star cluster has been studied extensively, how it formed is uncertain. Studies have shown it contains a solar and supersolar metallicity population that may have formed in-situ, along with a subsolar metallicity population that may have formed via mergers of globular clusters and dwarf galaxies. Stellar abundance measurements are critical to differentiate between formation scenarios. We present new measurements of [$M/H$] and $α$-element abundances [$α/Fe$] of two subsolar-metallicity stars in the Galactic Center. These observations were taken with the adaptive-optics assisted high-resolution (R=24,000) spectrograph NIRSPEC in the K-band (1.8 - 2.6 micron). These are the first $α$-element abundance measurements of sub-solar metallicity stars in the Milky Way nuclear star cluster. We measure [$M/H$]=$-0.59\pm 0.11$, [$α/Fe$]=$0.05\pm 0.15$ and [$M/H$]= $-0.81\pm 0.12$, [$α/Fe$]= $0.15\pm 0.16$ for the two stars at the Galactic center; the uncertainties are dominated by systematic uncertainties in the spectral templates. The stars have an [$α/Fe$] in-between the [$α/Fe$] of globular clusters and dwarf galaxies at similar [$M/H$] values. Their abundances are very different than the bulk of the stars in the nuclear star cluster. These results indicate that the sub-solar metallicity population in the Milky Way nuclear star cluster likely originated from infalling dwarf galaxies or globular clusters and are unlikely to have formed in-situ. △ Less

Submitted 5 December, 2021; originally announced December 2021.

Comments: Accepted by ApJ

Anisotropic hyperbolic inflation for a model of two scalar and two vector fields

Authors: Tuan Q. Do, W. F. Kao

Abstract: In this paper, we extend a recent proposed model of two scalar and two vector fields to a hyperbolic inflation scenario, in which the field space of two scalar fields is a hyperbolic space instead of a flat space. In this model, one of the scalar fields is assumed to be a radial field, while the other is set as an angular field. Furthermore, both scalar fields will be coupled to two different vect… ▽ More In this paper, we extend a recent proposed model of two scalar and two vector fields to a hyperbolic inflation scenario, in which the field space of two scalar fields is a hyperbolic space instead of a flat space. In this model, one of the scalar fields is assumed to be a radial field, while the other is set as an angular field. Furthermore, both scalar fields will be coupled to two different vector fields, respectively. As a result, we are able to obtain a set of exact Bianchi type I solutions to this model. Stability analysis is also performed to show that this set of anisotropic solutions is indeed stable and attractive during the inflationary phase. This result indicates that the cosmic no-hair conjecture is extensively violated in this anisotropic hyperbolic inflation model. △ Less

Submitted 9 February, 2022; v1 submitted 26 October, 2021; originally announced October 2021.

Comments: 15 pages, 2 figures. Matches the published version. All calculations are not modified. Comments are welcome

Journal ref: Eur. Phys. J. C 82, 123 (2022)

Anisotropic constant-roll inflation for the Dirac-Born-Infeld model

Authors: Duy H. Nguyen, Tuyen M. Pham, Tuan Q. Do

Abstract: In this paper, we study a non-canonical extension of a supergravity-motivated model acting as a vivid counterexample to the cosmic no-hair conjecture due to its unusual coupling between scalar and electromagnetic fields. In particular, a canonical scalar field is replaced by the string-inspired Dirac-Born-Infeld one in this extension. As a result, exact anisotropic inflationary solutions for this… ▽ More In this paper, we study a non-canonical extension of a supergravity-motivated model acting as a vivid counterexample to the cosmic no-hair conjecture due to its unusual coupling between scalar and electromagnetic fields. In particular, a canonical scalar field is replaced by the string-inspired Dirac-Born-Infeld one in this extension. As a result, exact anisotropic inflationary solutions for this Dirac-Born-Infeld model are figured out under a constant-roll condition. Furthermore, numerical calculations are performed to verify that these anisotropic constant-roll solutions are indeed attractive during their inflationary phase. △ Less

Submitted 22 September, 2021; v1 submitted 29 July, 2021; originally announced July 2021.

Comments: 14 pages, 14 figures. Matches the published version. All results are not changed. Comments are welcome

Journal ref: Eur. Phys. J. C 81, 839 (2021)

k-Gauss-Bonnet inflation

Authors: Tuyen M. Pham, Duy H. Nguyen, Tuan Q. Do

Abstract: We propose a novel {\it k}-Gauss-Bonnet inflationary model, in which a kinetic term of scalar field is allowed to non-minimally couple to the Gauss-Bonnet topological invariant in the absence of a potential of scalar field. As a result, this model is shown to admit an isotropic power-law inflation provided that the scalar field is phantom. Furthermore, stability analysis based on the dynamical sys… ▽ More We propose a novel {\it k}-Gauss-Bonnet inflationary model, in which a kinetic term of scalar field is allowed to non-minimally couple to the Gauss-Bonnet topological invariant in the absence of a potential of scalar field. As a result, this model is shown to admit an isotropic power-law inflation provided that the scalar field is phantom. Furthermore, stability analysis based on the dynamical system method is performed to indicate that this inflation solution is indeed stable and attractive. More interestingly, a gradient instability in tensor perturbations is shown to disappear in this model. △ Less

Submitted 15 February, 2022; v1 submitted 13 July, 2021; originally announced July 2021.

Comments: 9 pages, single column, 2 figures. Updated version, in which many discussions along with relevant references are added. All calculations are not modified. Comments are welcome

Anisotropic power-law inflation for a model of two scalar and two vector fields

Authors: Tuan Q. Do, W. F. Kao

Abstract: Inspired by an interesting counterexample to the cosmic no-hair conjecture found in a supergravity-motivated model recently, we propose a multi-field extension, in which two scalar fields are allowed to non-minimally couple to two vector fields, respectively. This model is shown to admit an exact Bianchi type I power-law solution. Furthermore, stability analysis based on the dynamical system metho… ▽ More Inspired by an interesting counterexample to the cosmic no-hair conjecture found in a supergravity-motivated model recently, we propose a multi-field extension, in which two scalar fields are allowed to non-minimally couple to two vector fields, respectively. This model is shown to admit an exact Bianchi type I power-law solution. Furthermore, stability analysis based on the dynamical system method is performed to show that this anisotropic solution is indeed stable and attractive if both scalar fields are canonical. Nevertheless, if one of the two scalar fields is phantom then the corresponding anisotropic power-law inflation turns unstable as expected. △ Less

Submitted 21 June, 2021; v1 submitted 29 April, 2021; originally announced April 2021.

Comments: 13 pages, 2 figures. Updated version, in which several typos are corrected and one reference is added. All results are not changed. Matches the published version. Comments are welcome

Journal ref: Eur. Phys. J. C 81, 525 (2021)

Upper Limit on Brackett-gamma Emission from the Immediate Accretion Flow onto the Galactic Black Hole

Authors: Anna Ciurlo, Mark R. Morris, Randall D. Campbell, Andrea M. Ghez, Tuan Do, Devin S. Chu

Abstract: We present the first observational constraint on the Br-gamma recombination line emission associated with the supermassive black hole at the center of our Galaxy, known as Sgr A*. By combining 13 years of data with the Adaptive Optics fed integral field spectrograph OSIRIS at the W. M. Keck Observatory obtained as part of the Galactic Center Orbits Initiative, we extract the near-infrared spectrum… ▽ More We present the first observational constraint on the Br-gamma recombination line emission associated with the supermassive black hole at the center of our Galaxy, known as Sgr A*. By combining 13 years of data with the Adaptive Optics fed integral field spectrograph OSIRIS at the W. M. Keck Observatory obtained as part of the Galactic Center Orbits Initiative, we extract the near-infrared spectrum within ~0.2'' of the black hole and we derive an upper limit on the Br-gamma flux. The aperture was set to match the size of the disk-like structure that was recently reported based on millimeter-wave ALMA observations of the hydrogen recombination line, H30-alpha. Our stringent upper limit is at least a factor of 80 (and up to a factor of 245) below what would be expected from the ALMA measurements and strongly constrains possible interpretation of emission from this highly under-luminous supermassive black hole. △ Less

Submitted 25 February, 2021; originally announced February 2021.

Comments: Accepted by ApJ

No-go theorem for inflation in an extended Ricci-inverse gravity model

Authors: Tuan Q. Do

Abstract: In this paper, we propose an extension of the Ricci-inverse gravity, which has been proposed recently as a very novel type of fourth-order gravity, by introducing a second order term of the so-called anticurvature scalar as a correction. The main purpose of this paper is that we would like to see whether the extended Ricci-inverse gravity model admits the homogeneous and isotropic Friedmann-Lemait… ▽ More In this paper, we propose an extension of the Ricci-inverse gravity, which has been proposed recently as a very novel type of fourth-order gravity, by introducing a second order term of the so-called anticurvature scalar as a correction. The main purpose of this paper is that we would like to see whether the extended Ricci-inverse gravity model admits the homogeneous and isotropic Friedmann-Lemaitre-Robertson-Walker metric as its stable inflationary solution. However, a no-go theorem for inflation in this extended Ricci-inverse gravity is shown to appear through a stability analysis based on the dynamical system method. As a result, this no-go theorem implies that it is impossible to have such stable inflation in this extended Ricci-inverse gravity model. △ Less

Submitted 10 January, 2022; v1 submitted 21 January, 2021; originally announced January 2021.

Comments: 14 pages, 3 figures. Updated version with additional discussions on the Minkowskian limit. All calculations are not modified. Matches the published version. Comments are welcome

Journal ref: Eur. Phys. J. C 82, 15 (2022)

Rapid Variability of Sgr A* across the Electromagnetic Spectrum

Authors: G. Witzel, G. Martinez, S. P. Willner, E. E. Becklin, 4 H. Boyce, T. Do, A. Eckart, G. G. Fazio, A. Ghez, M. A. Gurwell, D. Haggard, R. Herrero-Illana, J. L. Hora, Z. Li, J. Liu, N. Marchili, Mark R. Morris, Howard A. Smith, M. Subroweit, J. A. Zensus

Abstract: Sagittarius A* (Sgr A*) is the variable radio, near-infrared (NIR), and X-ray source associated with accretion onto the Galactic center black hole. We have analyzed a comprehensive submillimeter (including new observations simultaneous with NIR monitoring), NIR, and 2-8 keV dataset. Submillimeter variations tend to lag those in the NIR by $\sim$30 minutes. An approximate Bayesian computation (ABC)… ▽ More Sagittarius A* (Sgr A*) is the variable radio, near-infrared (NIR), and X-ray source associated with accretion onto the Galactic center black hole. We have analyzed a comprehensive submillimeter (including new observations simultaneous with NIR monitoring), NIR, and 2-8 keV dataset. Submillimeter variations tend to lag those in the NIR by $\sim$30 minutes. An approximate Bayesian computation (ABC) fit to the X-ray first-order structure function shows significantly less power at short timescales in the X-rays than in the NIR. Less X-ray variability at short timescales combined with the observed NIR-X-ray correlations means the variability can be described as the result of two strictly correlated stochastic processes, the X-ray process being the low-pass-filtered version of the NIR process. The NIR--X-ray linkage suggests a simple radiative model: a compact, self-absorbed synchrotron sphere with high-frequency cutoff close to NIR frequencies plus a synchrotron self-Compton scattering component at higher frequencies. This model, with parameters fit to the submillimeter, NIR, and X-ray structure functions, reproduces the observed flux densities at all wavelengths, the statistical properties of all light curves, and the time lags between bands. The fit also gives reasonable values for physical parameters such as magnetic flux density $B\approx13$ G, source size $L \approx2.2R_{S}$, and high-energy electron density $n_{e}\approx4\times10^{7}$ cm$^{-3}$. An animation illustrates typical light curves, and we make public the parameter chain of our Bayesian analysis, the model implementation, and the visualization code. △ Less

Submitted 7 June, 2021; v1 submitted 18 November, 2020; originally announced November 2020.

Comments: Accepted by ApJS, June 3, 2021. The journal version of figure 21 is animated. The animation can also be found here: https://doi.org/10.17617/1.kctx3s25. This version (version 2) has been revised according to the referee's suggestions and includes optimized figures and some changes and corrections of the text; no scientific conclusions have changed

General relativistic MHD simulations of non-thermal flaring in Sagittarius A*

Authors: Koushik Chatterjee, Sera Markoff, Joseph Neilsen, Ziri Younsi, Gunther Witzel, Alexander Tchekhovskoy, Doosoo Yoon, Adam Ingram, Michiel van der Klis, Hope Boyce, Tuan Do, Daryl Haggard, Michael Nowak

Abstract: Sagittarius A* exhibits regular variability in its multiwavelength emission, including daily X-ray flares and roughly continuous near-infrared (NIR) flickering. The origin of this variability is still ambiguous since both inverse Compton and synchrotron emission are possible radiative mechanisms. The underlying particle distributions are also not well constrained, particularly the non-thermal cont… ▽ More Sagittarius A* exhibits regular variability in its multiwavelength emission, including daily X-ray flares and roughly continuous near-infrared (NIR) flickering. The origin of this variability is still ambiguous since both inverse Compton and synchrotron emission are possible radiative mechanisms. The underlying particle distributions are also not well constrained, particularly the non-thermal contribution. In this work, we employ the GPU-accelerated general relativistic magnetohydrodynamics (GRMHD) code H-AMR perform a study of flare flux distributions, including the effect of particle acceleration for the first time in high-resolution 3D simulations of Sgr A*. For the particle acceleration, we use the general relativistic ray-tracing (GRRT) code BHOSS to perform the radiative transfer, assuming a hybrid thermal+non-thermal electron energy distribution. We extract ~60 h lightcurves in the sub-millimetre, NIR and X-ray wavebands, and compare the power spectra and the cumulative flux distributions of the lightcurves to statistical descriptions for Sgr A* flares. Our results indicate that non-thermal populations of electrons arising from turbulence-driven reconnection in weakly magnetised accretion flows lead to moderate NIR and X-ray flares and reasonably describe the X-ray flux distribution while fulfilling multiwavelength flux constraints. These models exhibit high rms% amplitudes, >~150% both in the NIR and the X-rays, with changes in the accretion rate driving the 230~GHz flux variability, in agreement with Sgr A* observations. △ Less

Submitted 15 September, 2021; v1 submitted 17 November, 2020; originally announced November 2020.

Comments: Accepted by MNRAS

No-go theorem for inflation in Ricci-inverse gravity

Authors: Tuan Q. Do

Abstract: We study a recent proposed Ricci-inverse gravity, which is a very novel type of fourth-order gravity. In particular, we are able to figure out both isotropically and anisotropically inflating universes to this model. More interestingly, these solutions are free from a singularity problem. However, stability analysis based on the dynamical system method shows that both isotropic and anisotropic inf… ▽ More We study a recent proposed Ricci-inverse gravity, which is a very novel type of fourth-order gravity. In particular, we are able to figure out both isotropically and anisotropically inflating universes to this model. More interestingly, these solutions are free from a singularity problem. However, stability analysis based on the dynamical system method shows that both isotropic and anisotropic inflation of this model turn out to be unstable against field perturbations. This result implies a no-go theorem for both isotropic and anisotropic inflation in the Ricci-inverse gravity. △ Less

Submitted 10 May, 2021; v1 submitted 14 September, 2020; originally announced September 2020.

Comments: 15 pages, 2 figures. This is a minor revision, in which the title and abstract are modified in order to express more clearly the main result of this paper. In addition, discussions on some important results investigated by other people are added. Some relevant references are also added. All calculations and conclusions are not changed. Accepted for publication in EPJC. Comments are welcome

Journal ref: Eur. Phys. J. C 81, 431 (2021)

Revealing the Formation of the Milky Way Nuclear Star Cluster via Chemo-Dynamical Modeling

Authors: Tuan Do, Gregory David Martinez, Wolfgang Kerzendorf, Anja Feldmeier-Krause, Manuel Arca Sedda, Nadine Neumayer, Alessia Gualandris

Abstract: The Milky Way nuclear star cluster (MW NSC) has been used as a template to understand the origin and evolution of galactic nuclei and the interaction of nuclear star clusters with supermassive black holes. It is the only nuclear star cluster with a supermassive black hole where we can resolve individual stars to measure their kinematics and metal abundance to reconstruct its formation history. Her… ▽ More The Milky Way nuclear star cluster (MW NSC) has been used as a template to understand the origin and evolution of galactic nuclei and the interaction of nuclear star clusters with supermassive black holes. It is the only nuclear star cluster with a supermassive black hole where we can resolve individual stars to measure their kinematics and metal abundance to reconstruct its formation history. Here, we present results of the first chemo-dynamical model of the inner 1 pc of the MW NSC using metallicity and radial velocity data from the KMOS spectrograph on the Very Large Telescope. We find evidence for two kinematically and chemically distinct components in this region. The majority of the stars belong to a previously known super-solar metallicity component with a rotation axis perpendicular to the Galactic plane. However, we identify a new kinematically distinct sub-solar metallicity component which contains about 7\% of the stars and appears to be rotating faster than the main component with a rotation axis that may be misaligned. This second component may be evidence for an infalling star cluster or remnants of a dwarf galaxy, merging with the MW NSC. These measurements show that the combination of chemical abundances with kinematics is a promising method to directly study the MW NSC's origin and evolution. △ Less

Submitted 4 September, 2020; originally announced September 2020.

Comments: 6 pages, 5 figures, accepted to ApJ Letters

On the origin of a rotating metal-poor stellar population in the Milky Way Nuclear Cluster

Authors: Manuel Arca Sedda, Alessia Gualandris, Tuan Do, Anja Feldmeier-Krause, Nadine Neumayer, Denis Erkal

Abstract: We explore the origin of a population of stars recently detected in the inner parsec of the Milky Way Nuclear Cluster (NC), which exhibit sub-solar metallicity and a higher rotation compared to the dominant population. Using state-of-the-art $N$-body simulations, we model the infall of a massive stellar system into the Galactic center, both of Galactic and extra-galactic origin. We show that the n… ▽ More We explore the origin of a population of stars recently detected in the inner parsec of the Milky Way Nuclear Cluster (NC), which exhibit sub-solar metallicity and a higher rotation compared to the dominant population. Using state-of-the-art $N$-body simulations, we model the infall of a massive stellar system into the Galactic center, both of Galactic and extra-galactic origin. We show that the newly discovered population can either be the remnant of a massive star cluster formed a few kpc away from the Galactic center (Galactic scenario) or be accreted from a dwarf galaxy originally located at 10-100 kpc (extragalactic scenario) and that reached the Galactic center 3-5 Gyr ago. A comparison between our models and characteristic Galactocentric distance and metallicity distributions of Milky Way satellites and globular clusters favours the Galactic scenario. A comparison with clusters associated with the Enceladus-Sausage, Sequoia, Sagittarius and Canis Major structures suggests that the progenitor of the observed metal-poor substructure formed in-situ rather than being accreted. △ Less

Submitted 28 September, 2020; v1 submitted 4 September, 2020; originally announced September 2020.

Comments: 16 pages, 9 figures, 1 table, 2 online data files available upon publication. published in ApJL

On Socially Distant Neighbors: Using Binaries to Constrain the Density of Objects in the Galactic Center

Authors: Sanaea C. Rose, Smadar Naoz, Abhimat K. Gautam, Andrea M. Ghez, Tuan Do, Devin Chu, Eric Becklin

Abstract: Stars often reside in binary configurations. The nuclear star cluster surrounding the supermassive black hole (SMBH) in the Galactic Center (GC) is expected to include a binary population. In this dense environment, a binary frequently encounters and interacts with neighboring stars. These interactions vary from small perturbations to violent collisions. In the former case, weak gravitational inte… ▽ More Stars often reside in binary configurations. The nuclear star cluster surrounding the supermassive black hole (SMBH) in the Galactic Center (GC) is expected to include a binary population. In this dense environment, a binary frequently encounters and interacts with neighboring stars. These interactions vary from small perturbations to violent collisions. In the former case, weak gravitational interactions unbind a soft binary over the evaporation timescale, which depends on the binary properties as well as the density of surrounding objects and velocity dispersion. Similarly, collisions can also unbind a binary, and the collision rate depends on the density. Thus, the detection of a binary with known properties can constrain the density profile in the GC with implications for the number of compact objects, which are otherwise challenging to detect. We estimate the density necessary to unbind a binary within its lifetime for an orbit of arbitrary eccentricity about the SMBH. We find that the eccentricity has a minimal impact on the density constraint. In this proof of concept, we demonstrate that this procedure can probe the density in the GC using hypothetical young and old binaries as examples. Similarly, a known density profile provides constraints on the binary orbital separation. Our results highlight the need to consider multiple dynamical processes in tandem. In certain cases, often closer to the SMBH, the collision timescale rather than the evaporation timescale gives the more stringent density constraint, while other binaries farther from the SMBH provide unreliable density constraints because they migrate inwards due to mass segregation. △ Less

Submitted 2 November, 2020; v1 submitted 14 August, 2020; originally announced August 2020.

Comments: Accepted to ApJ, 20 pages, 11 figures

Stable small spatial hairs in a power-law k-inflation model

Authors: Tuan Q. Do

Abstract: In this paper, we extend our investigation of the validity of the cosmic no-hair conjecture within non-canonical anisotropic inflation. As a result, we are able to figure out an exact Bianchi type I solution to a power-law {\it k}-inflation model in the presence of unusual coupling between scalar and electromagnetic fields as $-f^2(φ)F_{μν}F^{μν}/4$. Furthermore, stability analysis based on the dy… ▽ More In this paper, we extend our investigation of the validity of the cosmic no-hair conjecture within non-canonical anisotropic inflation. As a result, we are able to figure out an exact Bianchi type I solution to a power-law {\it k}-inflation model in the presence of unusual coupling between scalar and electromagnetic fields as $-f^2(φ)F_{μν}F^{μν}/4$. Furthermore, stability analysis based on the dynamical system method indicates that the obtained solution does admit stable and attractive hairs during an inflationary phase and therefore violates the cosmic no-hair conjecture. Finally, we show that the corresponding tensor-to-scalar ratio of this model turns out to be highly consistent with the observational data of the Planck 2018. △ Less

Submitted 24 November, 2020; v1 submitted 9 July, 2020; originally announced July 2020.

Comments: 15 pages, 2 figures. Substantial revision with discussions on seminal papers of Starobinsky related to the cosmic no-hair conjecture and investigations on the corresponding tensor-to-scalar ratio. Main results and conclusions are not changed. Comments are welcome

Journal ref: Eur. Phys. J. C 81, 77 (2021)

The former companion of the hyper-velocity star S5-HVS1

Authors: Wenbin Lu, Jim Fuller, Yael Raveh, Hagai B. Perets, Ting S. Li, Matthew W. Hosek Jr., Tuan Do

Abstract: The hyper-velocity star S5-HVS1, ejected 5 Myr ago from the Galactic Center at 1800 km/s, was most likely produced by tidal break-up of a tight binary by the supermassive black hole SgrA*. Taking a Monte Carlo approach, we show that the former companion of S5-HVS1 was likely a main-sequence star between 1.2 and 6 solar masses and was captured into a highly eccentric orbit with pericenter distance… ▽ More The hyper-velocity star S5-HVS1, ejected 5 Myr ago from the Galactic Center at 1800 km/s, was most likely produced by tidal break-up of a tight binary by the supermassive black hole SgrA*. Taking a Monte Carlo approach, we show that the former companion of S5-HVS1 was likely a main-sequence star between 1.2 and 6 solar masses and was captured into a highly eccentric orbit with pericenter distance in the range 1-10 AU and semimajor axis about $10^3$ AU. We then explore the fate of the captured star. We find that the heat deposited by tidally excited stellar oscillation modes leads to runaway disruption if the pericenter distance is smaller than about 3 AU. Over the past 5 Myr, its angular momentum has been significantly modified by orbital relaxation, which may stochastically drive the pericenter inwards below 3 AU and cause tidal disruption. We find an overall survival probability in the range 5% to 50%, depending on the local relaxation time in the close environment of the captured star, and the initial pericenter at capture. The pericenter distance of the surviving star has migrated to 10-100 AU, making it potentially the most extreme member of the S-star cluster. From the ejection rate of S5-HVS1-like stars, we estimate that there may currently be a few stars in such highly eccentric orbits. They should be detectable (typically Ks < 18.5 mag) by the GRAVITY instrument and by future Extremely Large Telescopes and hence provide an extraordinary probe of the spin of SgrA*. △ Less

Submitted 13 February, 2021; v1 submitted 25 May, 2020; originally announced May 2020.

Comments: MNRAS accepted after minor revision