Showing 1–50 of 67 results for author: Witzel, G

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

Ordered magnetic fields around the 3C 84 central black hole

Authors: G. F. Paraschos, J. -Y. Kim, M. Wielgus, J. Röder, T. P. Krichbaum, E. Ros, I. Agudo, I. Myserlis, M. Moscibrodzka, E. Traianou, J. A. Zensus, L. Blackburn, C. -K. Chan, S. Issaoun, M. Janssen, M. D. Johnson, V. L. Fish, K. Akiyama, A. Alberdi, W. Alef, J. C. Algaba, R. Anantua, K. Asada, R. Azulay, U. Bach , et al. (258 additional authors not shown)

Abstract: 3C84 is a nearby radio source with a complex total intensity structure, showing linear polarisation and spectral patterns. A detailed investigation of the central engine region necessitates the use of VLBI above the hitherto available maximum frequency of 86GHz. Using ultrahigh resolution VLBI observations at the highest available frequency of 228GHz, we aim to directly detect compact structures a… ▽ More 3C84 is a nearby radio source with a complex total intensity structure, showing linear polarisation and spectral patterns. A detailed investigation of the central engine region necessitates the use of VLBI above the hitherto available maximum frequency of 86GHz. Using ultrahigh resolution VLBI observations at the highest available frequency of 228GHz, we aim to directly detect compact structures and understand the physical conditions in the compact region of 3C84. We used EHT 228GHz observations and, given the limited (u,v)-coverage, applied geometric model fitting to the data. We also employed quasi-simultaneously observed, multi-frequency VLBI data for the source in order to carry out a comprehensive analysis of the core structure. We report the detection of a highly ordered, strong magnetic field around the central, SMBH of 3C84. The brightness temperature analysis suggests that the system is in equipartition. We determined a turnover frequency of $ν_m=(113\pm4)$GHz, a corresponding synchrotron self-absorbed magnetic field of $B_{SSA}=(2.9\pm1.6)$G, and an equipartition magnetic field of $B_{eq}=(5.2\pm0.6)$G. Three components are resolved with the highest fractional polarisation detected for this object ($m_\textrm{net}=(17.0\pm3.9)$%). The positions of the components are compatible with those seen in low-frequency VLBI observations since 2017-2018. We report a steeply negative slope of the spectrum at 228GHz. We used these findings to test models of jet formation, propagation, and Faraday rotation in 3C84. The findings of our investigation into different flow geometries and black hole spins support an advection-dominated accretion flow in a magnetically arrested state around a rapidly rotating supermassive black hole as a model of the jet-launching system in the core of 3C84. However, systematic uncertainties due to the limited (u,v)-coverage, however, cannot be ignored. △ Less

Submitted 1 February, 2024; originally announced February 2024.

Comments: 15 pages, 6 figures, published in A&A

Journal ref: Issue: A&A Volume 682, February 2024; Article number: L3; Number of pages: 15

A Collection of German Science Interests in the Next Generation Very Large Array

Authors: M. Kadler, D. A. Riechers, J. Agarwal, A. -K. Baczko, H. Beuther, F. Bigiel, T. Birnstiel, B. Boccardi, D. J. Bomans, L. Boogaard, T. T. Braun, S. Britzen, M. Brüggen, A. Brunthaler, P. Caselli, D. Elsässer, S. von Fellenberg, M. Flock, C. M. Fromm, L. Fuhrmann, P. Hartogh, M. Hoeft, R. P. Keenan, Y. Kovalev, K. Kreckel , et al. (66 additional authors not shown)

Abstract: The Next Generation Very Large Array (ngVLA) is a planned radio interferometer providing unprecedented sensitivity at wavelengths between 21 cm and 3 mm. Its 263 antenna element array will be spatially distributed across North America to enable both superb low surface brightness recovery and sub-milliarcsecond angular resolution imaging. The project was developed by the international astronomy com… ▽ More The Next Generation Very Large Array (ngVLA) is a planned radio interferometer providing unprecedented sensitivity at wavelengths between 21 cm and 3 mm. Its 263 antenna element array will be spatially distributed across North America to enable both superb low surface brightness recovery and sub-milliarcsecond angular resolution imaging. The project was developed by the international astronomy community under the lead of the National Radio Astronomy Observatory (NRAO), and is anticipated to be built between 2027 and 2037. Two workshops have been held in 2022 and 2023 with the goal to discuss and consolidate the scientific interests in the ngVLA within the German astronomical community. This community paper constitutes a collection of 48 science ideas which the German community aims to pursue with the ngVLA in the 2030s. This is not a complete list and the ideas are not developed at the level of a "Science Book", such that the present document is mainly meant provide a basis for further discussion within the community. As such, additional contributions are welcome, and will be considered for inclusion in future revisions. △ Less

Submitted 18 June, 2024; v1 submitted 16 November, 2023; originally announced November 2023.

Comments: Version 2.0 (status June 18, 2024): 169 pages, comments and future contributions welcome [v2.0: 7 new science cases added, some minor revisions to other chapters]

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

A search for pulsars around Sgr A* in the first Event Horizon Telescope dataset

Authors: Pablo Torne, Kuo Liu, Ralph P. Eatough, Jompoj Wongphechauxsorn, James M. Cordes, Gregory Desvignes, Mariafelicia De Laurentis, Michael Kramer, Scott M. Ransom, Shami Chatterjee, Robert Wharton, Ramesh Karuppusamy, Lindy Blackburn, Michael Janssen, Chi-kwan Chan, Geoffrey B. Crew, Lynn D. Matthews, Ciriaco Goddi, Helge Rottmann, Jan Wagner, Salvador Sanchez, Ignacio Ruiz, Federico Abbate, Geoffrey C. Bower, Juan J. Salamanca , et al. (261 additional authors not shown)

Abstract: The Event Horizon Telescope (EHT) observed in 2017 the supermassive black hole at the center of the Milky Way, Sagittarius A* (Sgr A*), at a frequency of 228.1 GHz ($λ$=1.3 mm). The fundamental physics tests that even a single pulsar orbiting Sgr A* would enable motivate searching for pulsars in EHT datasets. The high observing frequency means that pulsars - which typically exhibit steep emission… ▽ More The Event Horizon Telescope (EHT) observed in 2017 the supermassive black hole at the center of the Milky Way, Sagittarius A* (Sgr A*), at a frequency of 228.1 GHz ($λ$=1.3 mm). The fundamental physics tests that even a single pulsar orbiting Sgr A* would enable motivate searching for pulsars in EHT datasets. The high observing frequency means that pulsars - which typically exhibit steep emission spectra - are expected to be very faint. However, it also negates pulse scattering, an effect that could hinder pulsar detections in the Galactic Center. Additionally, magnetars or a secondary inverse Compton emission could be stronger at millimeter wavelengths than at lower frequencies. We present a search for pulsars close to Sgr A* using the data from the three most-sensitive stations in the EHT 2017 campaign: the Atacama Large Millimeter/submillimeter Array, the Large Millimeter Telescope and the IRAM 30 m Telescope. We apply three detection methods based on Fourier-domain analysis, the Fast-Folding-Algorithm and single pulse search targeting both pulsars and burst-like transient emission; using the simultaneity of the observations to confirm potential candidates. No new pulsars or significant bursts were found. Being the first pulsar search ever carried out at such high radio frequencies, we detail our analysis methods and give a detailed estimation of the sensitivity of the search. We conclude that the EHT 2017 observations are only sensitive to a small fraction ($\lesssim$2.2%) of the pulsars that may exist close to Sgr A*, motivating further searches for fainter pulsars in the region. △ Less

Submitted 29 August, 2023; originally announced August 2023.

Comments: 33 pages, 7 figures, 6 Tables. Accepted for publication in ApJ

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

Comparison of Polarized Radiative Transfer Codes used by the EHT Collaboration

Authors: Ben S. Prather, Jason Dexter, Monika Moscibrodzka, Hung-Yi Pu, Thomas Bronzwaer, Jordy Davelaar, Ziri Younsi, Charles F. Gammie, Roman Gold, George N. Wong, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Uwe Bach, Anne-Kathrin Baczko, David Ball, Mislav Baloković, John Barrett, Michi Bauböck, Bradford A. Benson, Dan Bintley , et al. (248 additional authors not shown)

Abstract: Interpretation of resolved polarized images of black holes by the Event Horizon Telescope (EHT) requires predictions of the polarized emission observable by an Earth-based instrument for a particular model of the black hole accretion system. Such predictions are generated by general relativistic radiative transfer (GRRT) codes, which integrate the equations of polarized radiative transfer in curve… ▽ More Interpretation of resolved polarized images of black holes by the Event Horizon Telescope (EHT) requires predictions of the polarized emission observable by an Earth-based instrument for a particular model of the black hole accretion system. Such predictions are generated by general relativistic radiative transfer (GRRT) codes, which integrate the equations of polarized radiative transfer in curved spacetime. A selection of ray-tracing GRRT codes used within the EHT collaboration is evaluated for accuracy and consistency in producing a selection of test images, demonstrating that the various methods and implementations of radiative transfer calculations are highly consistent. When imaging an analytic accretion model, we find that all codes produce images similar within a pixel-wise normalized mean squared error (NMSE) of 0.012 in the worst case. When imaging a snapshot from a cell-based magnetohydrodynamic simulation, we find all test images to be similar within NMSEs of 0.02, 0.04, 0.04, and 0.12 in Stokes I, Q, U , and V respectively. We additionally find the values of several image metrics relevant to published EHT results to be in agreement to much better precision than measurement uncertainties. △ Less

Submitted 21 March, 2023; originally announced March 2023.

Comments: Accepted for publication in ApJ

General relativistic effects and the near-infrared and X-ray variability of Sgr A* I

Authors: Sebastiano D. von Fellenberg, Gunther Witzel, Michi Bauböck, Hui-Hsuan Chung, Nicolás Aimar, Matteo Bordoni, Antonia Drescher, Frank Eisenhauer, Reinhard Genzel, Stefan Gillessen, Nicola Marchili, Thibaut Paumard, Guy Perrin, Thomas Ott, Diogo Ribeiro, Eduardo Ros, Frédéric Vincent, Felix Widmann, S. P. Willner, J. Anton Zensus

Abstract: The near-infrared (NIR) and X-ray emission of Sagittarius A* shows occasional bright flares that are assumed to originate from the innermost region of the accretion flow. We identified $25$ $4.5 μm$ and $24$ X-ray flares in archival data obtained with the \textit{Spitzer} and \textit{Chandra} observatories. With the help of general relativistic ray-tracing code, we modeled trajectories of ``hot sp… ▽ More The near-infrared (NIR) and X-ray emission of Sagittarius A* shows occasional bright flares that are assumed to originate from the innermost region of the accretion flow. We identified $25$ $4.5 μm$ and $24$ X-ray flares in archival data obtained with the \textit{Spitzer} and \textit{Chandra} observatories. With the help of general relativistic ray-tracing code, we modeled trajectories of ``hot spots'' and studied the light curves of the flares for signs of the effects of general relativity. Despite their apparent diversity in shape, all flares share a common, exponential impulse response, a characteristic shape that is the building block of the variability. This shape is symmetric, that is, the rise and fall times are the same. Furthermore, the impulse responses in the NIR and X-ray are identical within uncertainties, with an exponential time constant $τ\sim 15$ minute. The observed characteristic flare shape is inconsistent with hot-spot orbits viewed edge-on. Individually modeling the light curves of the flares, we derived constraints on the inclination of the orbital plane of the hot spots with respect to the observer ($i \sim 30^{\circ} , < 75^{\circ} $) and on the characteristic timescale of the intrinsic variability (tens of minutes). △ Less

Submitted 6 January, 2023; originally announced January 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)

A multiband study and exploration of the radio wave - $γ$-ray connection in 3C 84

Authors: G. F. Paraschos, V. Mpisketzis, J. -Y. Kim, G. Witzel, T. P. Krichbaum, J. A. Zensus, M. A. Gurwell, A. Lähteenmäki, M. Tornikoski, S. Kiehlmann, A. C. S. Readhead

Abstract: Total intensity variability light curves offer a unique insight into the ongoing debate about the launching mechanism of jets. For this work, we utilise the availability of radio and $γ$-ray light curves over a few decades of the radio source 3C 84 (NGC 1275). We calculate the multiband time lags between the flares identified in the light curves via discrete cross-correlation and Gaussian process… ▽ More Total intensity variability light curves offer a unique insight into the ongoing debate about the launching mechanism of jets. For this work, we utilise the availability of radio and $γ$-ray light curves over a few decades of the radio source 3C 84 (NGC 1275). We calculate the multiband time lags between the flares identified in the light curves via discrete cross-correlation and Gaussian process regression. We find that the jet particle and magnetic field energy densities are in equipartition ($k_\textrm{r} = 1.08\pm0.18$). The jet apex is located $z_\textrm{91.5 GHz}= 22 - 645$ $R_\textrm{s}$ ($2 - 20 \times 10^{-3}$ pc) upstream of the 3 mm radio core; at that position, the magnetic field amplitude is $B_\textrm{core}^\textrm{91.5 GHz}= 3 - 10$ G. Our results are in good agreement with earlier studies, which utilised very-long-baseline interferometry. Furthermore, we investigate the temporal relation between the ejection of radio and $γ$-ray flares. Our results are in favour of the $γ$-ray emission being associated with the radio emission. We are able to tentatively connect the ejection of features identified at 43 and 86 GHz to prominent $γ$-ray flares. Finally, we compute the multiplicity parameter $λ$ and the Michel magnetisation $σ_\textrm{M}$ and find that they are consistent with a jet launched by the Blandford & Znajek 1977 mechanism. △ Less

Submitted 18 October, 2022; originally announced October 2022.

Comments: 13 pages, 6 figures, accepted for publication in A&A

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

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

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.

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

Flaremodel: An open-source Python package for one-zone numerical modelling of synchrotron sources

Authors: Y. Dallilar, S. von Fellenberg, M. Bauböck, P. T. de Zeeuw, A. Drescher, F. Eisenhauer, R. Genzel, S. Gillessen, M. Habibi, T. Ott, G. Ponti, J. Stadler, O. Straub, F. Widmann, G. Witzel, A. Young

Abstract: Synchrotron processes, the radiative processes associated with the interaction of energetic charged particles with magnetic field, are of interest in many areas in astronomy, from the interstellar medium to extreme environments near compact objects. Consequently, observations of synchrotron sources carry information on the physical properties of the sources themselves and those of their close vici… ▽ More Synchrotron processes, the radiative processes associated with the interaction of energetic charged particles with magnetic field, are of interest in many areas in astronomy, from the interstellar medium to extreme environments near compact objects. Consequently, observations of synchrotron sources carry information on the physical properties of the sources themselves and those of their close vicinity. In recent years, novel observations of such sources with multi-wavelength collaborations reveal complex features and peculiarities, especially near black holes. Exploring the nature of these sources in more detail necessitates numerical tools complementary to analytical one-zone modelling efforts. In this paper, we introduce an open-source Python package tailored to this purpose, flaremodel. The core of the code consists of low-level utility functions to describe physical processes relevant to synchrotron sources, which are written in C for performance and parallelised with OpenMP for scalability. The Python interface provides access to these functions and built-in source models are provided as a guidance. At the same time, the modular design of the code and the generic nature of these functions enable users to build a variety of source models applicable to many astrophysical synchrotron sources. We describe our methodology and the structure of our code along with selected examples demonstrating capabilities and options for future modelling efforts. △ Less

Submitted 30 November, 2021; originally announced November 2021.

Journal ref: A&A 658, A111 (2022)

Second Scale Submillimeter Variability of Sagittarius A* during flaring activity of 2019: On the Origin of Bright Near Infrared Flares

Authors: Lena Murchikova, Gunther Witzel

Abstract: In 2019, Sgr A* - the supermassive black hole in the Galactic Center - underwent unprecedented flaring activity in the near infrared (NIR), brightening by up to a factor of 100 compared to quiescent values. Here we report ALMA observations of Sgr A*'s continuum variability at 1.3 mm (230 GHz) - a tracer of the accretion rate - conducted one month after the brightest detected NIR flare and in the m… ▽ More In 2019, Sgr A* - the supermassive black hole in the Galactic Center - underwent unprecedented flaring activity in the near infrared (NIR), brightening by up to a factor of 100 compared to quiescent values. Here we report ALMA observations of Sgr A*'s continuum variability at 1.3 mm (230 GHz) - a tracer of the accretion rate - conducted one month after the brightest detected NIR flare and in the middle of the flaring activity of 2019. We develop an innovative light curve extraction technique which (together with ALMA's excellent sensitivity) allows us to obtain the light curves which are simultaneously of high time resolution (2 seconds) and high signal-to-noise ratio (~ 500). We construct an accurate intrinsic structure function of the Sgr A* submm variability, improving on previous studies by about two orders of magnitude in timescale and one order of magnitude in sensitivity. We compare the June 2019 variability behavior with that of 2001-2017, and suggest that the most likely cause of the bright NIR flares is magnetic reconnection. △ Less

Submitted 7 September, 2021; v1 submitted 23 July, 2021; originally announced July 2021.

Comments: v2: several comments, clarifications, references, and a derivation added. 9 pages, 2 figures, 1 table. Accepted to ApJL. ; v1: 8 pages, 2 figures, 1 table

Constraining particle acceleration in Sgr A* with simultaneous GRAVITY, Spitzer, NuSTAR and Chandra observations

Authors: R. Abuter, A. Amorim, M. Bauböck, F. Baganoff, J. P. Berge, H. Boyce, H. Bonnet, W. Brandner, Y. Clénet, R. Davies, P. T. de Zeeuw, J. Dexter, Y. Dallilar, A. Drescher, A. Eckart, F. Eisenhauer, G. G. Fazio, N. M. Förster Schreiber, K. Foster, C. Gammie, P. Garcia, F. Gao, E. Gendron, R. Genzel, G. Ghisellini , et al. (59 additional authors not shown)

Abstract: We report the time-resolved spectral analysis of a bright near-infrared and moderate X-ray flare of Sgr A*. We obtained light curves in the $M$-, $K$-, and $H$-bands in the mid- and near-infrared and in the $2-8~\mathrm{keV}$ and $2-70~\mathrm{keV}$ bands in the X-ray. The observed spectral slope in the near-infrared band is $νL_ν\propto ν^{0.5\pm0.2}$; the spectral slope observed in the X-ray ban… ▽ More We report the time-resolved spectral analysis of a bright near-infrared and moderate X-ray flare of Sgr A*. We obtained light curves in the $M$-, $K$-, and $H$-bands in the mid- and near-infrared and in the $2-8~\mathrm{keV}$ and $2-70~\mathrm{keV}$ bands in the X-ray. The observed spectral slope in the near-infrared band is $νL_ν\propto ν^{0.5\pm0.2}$; the spectral slope observed in the X-ray band is $νL_ν\propto ν^{-0.7\pm0.5}$. We tested synchrotron and synchrotron self-Compton (SSC) scenarios. The observed near-infrared brightness and X-ray faintness, together with the observed spectral slopes, pose challenges for all models explored. We rule out a scenario in which the near-infrared emission is synchrotron emission and the X-ray emission is SSC. A one-zone model in which both the near-infrared and X-ray luminosity are produced by SSC and a model in which the luminosity stems from a cooled synchrotron spectrum can explain the flare. In order to describe the mean SED, both models require specific values of the maximum Lorentz factor $γ_{max}$, which however differ by roughly two orders of magnitude: the SSC model suggests that electrons are accelerated to $γ_{max}\sim 500$, while cooled synchrotron model requires acceleration up to $γ_{max}\sim5\times 10^{4}$. The SSC scenario requires electron densities of $10^{10}~\mathrm{cm^{-3}}$ much larger than typical ambient densities in the accretion flow, and thus require in an extraordinary accretion event. In contrast, assuming a source size of $1R_s$, the cooled synchrotron scenario can be realized with densities and magnetic fields comparable with the ambient accretion flow. For both models, the temporal evolution is regulated through the maximum acceleration factor $γ_{max}$, implying that sustained particle acceleration is required to explain at least a part of the temporal evolution of the flare. △ Less

Submitted 2 July, 2021; originally announced July 2021.

Comments: accepted for publication in Astronomy & Astrophysics; preview abstract shortened due to arXiv requirements

Journal ref: A&A 654, A22 (2021)

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

Kinematic Structure of the Galactic Center S-cluster

Authors: Basel Ali, Daria Paul, Andreas Eckart, Marzieh Parsa, Michal Zajacek, Florian Peißker, Matthias Subroweit, Monica Valencia-S., Lauritz Thomkins, Gunther Witzel

Abstract: We present a detailed analysis of the kinematics of 112 stars that mostly comprise the high velocity S-cluster and orbit the super massive black hole SgrA* at the center of the Milky Way. For 39 of them orbital elements are known, for the remainder we know proper motions. The distribution of inclinations, and proper motion flight directions deviate significantly from a uniform distribution which o… ▽ More We present a detailed analysis of the kinematics of 112 stars that mostly comprise the high velocity S-cluster and orbit the super massive black hole SgrA* at the center of the Milky Way. For 39 of them orbital elements are known, for the remainder we know proper motions. The distribution of inclinations, and proper motion flight directions deviate significantly from a uniform distribution which one expects if the orientation of the orbits are random. Across the central arcseconds the S-cluster stars are arranged in two almost edge on disks that are located at a position angle approximately +-45 o with respect to the Galactic plane. The angular momentum vectors for stars in each disk point in both directions, i.e. the stars in a given disk rotate in opposite ways. The poles of this structure are located only about 25 o from the line of sight. This structure may be the result of a resonance process that started with the formation of the young B-dwarf stars in the cluster about 6 Myr ago. Alternatively, it indicated the presence of a disturber at a distance from the center comparable to the distance of the compact stellar association IRS13. △ Less

Submitted 19 June, 2020; originally announced June 2020.

Comments: published in Astrophysical Journal 23 pages, 6 tabels, 16 figures, ENHANCED GRAPHICS in ApJ in fullorbits.gif circularized.gif normalized.gif

Consistency of the Infrared Variability of Sgr A* over 22 years

Authors: Zhuo Chen, E. Gallego-Cano, T. Do, G. Witzel, A. M. Ghez, R. Schödel, B. N. Sitarski, E. E. Becklin, J. Lu, M. R. Morris, A. Dehghanfar, A. K. Gautam, A. Hees, M. W. Hosek Jr, S. Jia, A. C. Mangian, K. Matthews

Abstract: We report new infrared measurements of the supermassive black hole at the Galactic Center, Sgr A*, over a decade that was previously inaccessible at these wavelengths. This enables a variability study that addresses variability timescales that are ten times longer than earlier published studies. Sgr A* was initially detected in the near-infrared with adaptive optics observations in 2002. While ear… ▽ More We report new infrared measurements of the supermassive black hole at the Galactic Center, Sgr A*, over a decade that was previously inaccessible at these wavelengths. This enables a variability study that addresses variability timescales that are ten times longer than earlier published studies. Sgr A* was initially detected in the near-infrared with adaptive optics observations in 2002. While earlier data exists in form of speckle imaging (1995 - 2005), Sgr A* was not detected in the initial analysis. Here, we improved our speckle holography analysis techniques. This has improved the sensitivity of the resulting speckle images by up to a factor of three. Sgr A* is now detectable in the majority of epochs covering 7 years. The brightness of Sgr A* in the speckle data has an average observed K magnitude of 16.0, which corresponds to a dereddened flux density of $3.4$ mJy. Furthermore, the flat power spectral density (PSD) of Sgr A* between $\sim$80 days and 7 years shows its uncorrelation in time beyond the proposed single power-law break of $\sim$245 minutes. We report that the brightness and its variability is consistent over 22 years. This analysis is based on simulations using Witzel et al. (2018) model to characterize infrared variability from 2006 to 2016. Finally, we note that the 2001 periapse of the extended, dusty object G1 had no apparent effect on the near-infrared emission from accretion flow onto Sgr A*. The result is consistent with G1 being a self-gravitating object rather than a disrupting gas cloud. △ Less

Submitted 21 August, 2019; originally announced August 2019.

Comments: 22 pages, 20 figures, accepted to ApJ Letters

Unprecedented variability of Sgr A* in NIR

Authors: Tuan Do, Gunther Witzel, Abhimat K. Gautam, Zhuo Chen, Andrea M. Ghez, Mark R. Morris, Eric E. Becklin, Anna Ciurlo, Matthew Hosek Jr., Gregory D. Martinez, Keith Matthews, Shoko Sakai, Rainer Schödel

Abstract: The electromagnetic counterpart to the Galactic center supermassive black hole, Sgr A*, has been observed in the near-infrared for over 20 years and is known to be highly variable. We report new Keck Telescope observations showing that Sgr A* reached much brighter flux levels in 2019 than ever measured at near-infrared wavelengths. In the K$^\prime$ band, Sgr A* reached flux levels of $\sim6$ mJy,… ▽ More The electromagnetic counterpart to the Galactic center supermassive black hole, Sgr A*, has been observed in the near-infrared for over 20 years and is known to be highly variable. We report new Keck Telescope observations showing that Sgr A* reached much brighter flux levels in 2019 than ever measured at near-infrared wavelengths. In the K$^\prime$ band, Sgr A* reached flux levels of $\sim6$ mJy, twice the level of the previously observed peak flux from $>13,000$ measurements over 130 nights with the VLT and Keck Telescopes. We also observe a factor of 75 change in flux over a 2-hour time span with no obvious color changes between 1.6 $μ$m and 2.1 $μ$m. The distribution of flux variations observed this year is also significantly different than the historical distribution. Using the most comprehensive statistical model published, the probability of a single night exhibiting peak flux levels observed this year, given historical Keck observations, is less than $0.3\%$. The probability to observe the flux levels similar to all 4 nights of data in 2019 is less than $0.05\%$. This increase in brightness and variability may indicate a period of heightened activity from Sgr A* or a change in its accretion state. It may also indicate that the current model is not sufficient to model Sgr A* at high flux levels and should be updated. Potential physical origins of Sgr A*'s unprecedented brightness may be from changes in the accretion-flow as a result of the star S0-2's closest passage to the black hole in 2018 or from a delayed reaction to the approach of the dusty object G2 in 2014. Additional multi-wavelength observations will be necessary to both monitor Sgr A* for potential state changes and to constrain the physical processes responsible for its current variability. △ Less

Submitted 5 August, 2019; originally announced August 2019.

Comments: 9 pages, 4 figures, accepted to ApJ Letters

Relativistic redshift of the star S0-2 orbiting the Galactic center supermassive black hole

Authors: Tuan Do, Aurelien Hees, Andrea Ghez, Gregory D. Martinez, Devin S. Chu, Siyao Jia, Shoko Sakai, Jessica R. Lu, Abhimat K. Gautam, Kelly Kosmo O'Neil, Eric E. Becklin, Mark R. Morris, Keith Matthews, Shogo Nishiyama, Randy Campbell, Samantha Chappell, Zhuo Chen, Anna Ciurlo, Arezu Dehghanfar, Eulalia Gallego-Cano, Wolfgang E. Kerzendorf, James E. Lyke, Smadar Naoz, Hiromi Saida, Rainer Schödel , et al. (4 additional authors not shown)

Abstract: General Relativity predicts that a star passing close to a supermassive black hole should exhibit a relativistic redshift. We test this using observations of the Galactic center star S0-2. We combine existing spectroscopic and astrometric measurements from 1995-2017, which cover S0-2's 16-year orbit, with measurements in 2018 March to September which cover three events during its closest approach… ▽ More General Relativity predicts that a star passing close to a supermassive black hole should exhibit a relativistic redshift. We test this using observations of the Galactic center star S0-2. We combine existing spectroscopic and astrometric measurements from 1995-2017, which cover S0-2's 16-year orbit, with measurements in 2018 March to September which cover three events during its closest approach to the black hole. We detect the combination of special relativistic- and gravitational-redshift, quantified using a redshift parameter, $Υ$. Our result, $Υ=0.88 \pm 0.17$, is consistent with General Relativity ($Υ=1$) and excludes a Newtonian model ($Υ=0$ ) with a statistical significance of 5 $σ$. △ Less

Submitted 24 July, 2019; originally announced July 2019.

Comments: 78 pages, 21 figures, accepted to Science, astrometry and radial velocity measurements, and posterior chains available here, http://www.astro.ucla.edu/~ghezgroup/gc/gr/data.zip

Envisioning the next decade of Galactic Center science: a laboratory for the study of the physics and astrophysics of supermassive black holes

Authors: Tuan Do, Andrea Ghez, Jessica R. Lu, Mark Morris, Matthew Hosek Jr., Aurelien Hees, Smadar Naoz, Anna Ciurlo, Philip J. Armitage, Rachael L Beaton, Eric Becklin, Andrea Bellini, Rory O. Bentley, Joss Bland-Hawthorn, Sukanya Chakrabarti, Zhuo Chen, Devin S. Chu, Arezu Dehghanfar, Charles F. Gammie, Abhimat K. Gautam, Reinhard Genzel, Jenny Greene, Daryl Haggard, Joseph Hora, Wolfgang E. Kerzendorf , et al. (16 additional authors not shown)

Abstract: As the closest example of a galactic nucleus, the Galactic center (GC) presents an exquisite laboratory for learning about supermassive black holes (SMBH) and their environment. We describe several exciting new research directions that, over the next 10 years, hold the potential to answer some of the biggest scientific questions raised in recent decades: Is General Relativity (GR) the correct desc… ▽ More As the closest example of a galactic nucleus, the Galactic center (GC) presents an exquisite laboratory for learning about supermassive black holes (SMBH) and their environment. We describe several exciting new research directions that, over the next 10 years, hold the potential to answer some of the biggest scientific questions raised in recent decades: Is General Relativity (GR) the correct description for supermassive black holes? What is the nature of star formation in extreme environments? How do stars and compact objects dynamically interact with the supermassive black hole? What physical processes drive gas accretion in low-luminosity black holes? We describe how the high sensitivity, angular resolution, and astrometric precision offered by the next generation of large ground-based telescopes with adaptive optics will help us answer these questions. First, it will be possible to obtain precision measurements of stellar orbits in the Galaxy's central potential, providing both tests of GR in the unexplored regime near a SMBH and measurements of the extended dark matter distribution that is predicted to exist at the GC. Second, we will probe stellar populations at the GC to significantly lower masses than are possible today, down to brown dwarfs. Their structure and dynamics will provide an unprecedented view of the stellar cusp around the SMBH and will distinguish between models of star formation in this extreme environment. This increase in depth will also allow us to measure the currently unknown population of compact remnants at the GC by observing their effects on luminous sources. Third, uncertainties on the mass of and distance to the SMBH can be improved by a factor of $\sim$10. Finally, we can also study the near-infrared accretion onto the black hole at unprecedented sensitivity and time resolution, which can reveal the underlying physics of black hole accretion. △ Less

Submitted 12 March, 2019; originally announced March 2019.

Comments: 9 pages, 4 figures, submitted for Astro2020 White Paper

The Galactic Center: An Improved Astrometric Reference Frame for Stellar Orbits around the Supermassive Black Hole

Authors: Shoko Sakai, Jessica R. Lu, Andrea Ghez, Siyao Jia, Tuan Do, Gunther Witzel, Abhimat K. Gautam, Aurelien Hees, E. Becklin, K. Matthews, M. W. Hosek Jr

Abstract: Precision measurements of the stars in short-period orbits around the supermassive black hole at the Galactic Center are now being used to constrain general relativistic effects, such as the gravitational redshift and periapse precession. One of the largest systematic uncertainties in the measured orbits has been errors in the astrometric reference frame, which is derived from seven infrared-brigh… ▽ More Precision measurements of the stars in short-period orbits around the supermassive black hole at the Galactic Center are now being used to constrain general relativistic effects, such as the gravitational redshift and periapse precession. One of the largest systematic uncertainties in the measured orbits has been errors in the astrometric reference frame, which is derived from seven infrared-bright stars associated with SiO masers that have extremely accurate radio positions, measured in the Sgr A*-rest frame. We have improved the astrometric reference frame within 14'' of the Galactic Center by a factor of 2.5 in position and a factor of 5 in proper motion. In the new reference frame, Sgr A* is localized to within a position of 0.645 mas and proper motion of 0.03 mas/yr. We have removed a substantial rotation (2.25 degrees per decade), that was present in the previous less-accurate reference frame used to measure stellar orbits in the field. With our improved methods and continued monitoring of the masers, we predict that orbital precession predicted by General Relativity will become detectable in the next ~5 years. △ Less

Submitted 24 January, 2019; originally announced January 2019.

Comments: 17 pages, 10 figures

Simultaneous X-ray and Infrared Observations of Sagittarius A*'s Variability

Authors: H. Boyce, D. Haggard, G. Witzel, S. P. Willner, J. Neilsen, J. L. Hora, S. Markoff, G. Ponti, F. Baganoff, E. Becklin, G. Fazio, P. Lowrance, M. R. Morris, H. A. Smith

Abstract: Emission from Sgr A* is highly variable at both X-ray and infrared (IR) wavelengths. Observations over the last ~20 years have revealed X-ray flares that rise above a quiescent thermal background about once per day, while faint X-ray flares from Sgr A* are undetectable below the constant thermal emission. In contrast, the IR emission of Sgr A* is observed to be continuously variable. Recently, sim… ▽ More Emission from Sgr A* is highly variable at both X-ray and infrared (IR) wavelengths. Observations over the last ~20 years have revealed X-ray flares that rise above a quiescent thermal background about once per day, while faint X-ray flares from Sgr A* are undetectable below the constant thermal emission. In contrast, the IR emission of Sgr A* is observed to be continuously variable. Recently, simultaneous observations have indicated a rise in IR flux density around the same time as every distinct X-ray flare, while the opposite is not always true (peaks in the IR emission may not be coincident with an X-ray flare). Characterizing the behaviour of these simultaneous X-ray/IR events and measuring any time lag between them can constrain models of Sgr A*'s accretion flow and the flare emission mechanism. Using 100+ hours of data from a coordinated campaign between the Spitzer Space Telescope and the Chandra X-ray Observatory, we present results of the longest simultaneous IR and X-ray observations of Sgr A* taken to date. The cross-correlation between the IR and X-ray light curves in this unprecedented dataset, which includes four modest X-ray/IR flares, indicates that flaring in the X-ray may lead the IR by approximately 10-20 minutes with 68% confidence. However, the 99.7% confidence interval on the time-lag also includes zero, i.e., the flaring remains statistically consistent with simultaneity. Long duration and simultaneous multiwavelength observations of additional bright flares will improve our ability to constrain the flare timing characteristics and emission mechanisms, and must be a priority for Galactic Center observing campaigns. △ Less

Submitted 13 December, 2018; originally announced December 2018.

Comments: 12 pages, 6 figures, 3 tables, accepted for publication in ApJ

An Adaptive Optics Survey of Stellar Variability at the Galactic Center

Authors: Abhimat Krishna Gautam., Tuan Do, Andrea M. Ghez, Mark R. Morris, Gregory D. Martinez, Matthew W. Hosek Jr., Jessica R. Lu, Shoko Sakai, Gunther Witzel, Siyao Jia, Eric E. Becklin, Keith Matthews

Abstract: We present a $\approx 11.5$ year adaptive optics (AO) study of stellar variability and search for eclipsing binaries in the central $\sim 0.4$ pc ($\sim 10''$) of the Milky Way nuclear star cluster. We measure the photometry of 563 stars using the Keck II NIRC2 imager ($K'$-band, $λ_0 = 2.124 \text{ } μ\text{m}$). We achieve a photometric uncertainty floor of $Δm_{K'} \sim 0.03$ ($\approx 3\%$), c… ▽ More We present a $\approx 11.5$ year adaptive optics (AO) study of stellar variability and search for eclipsing binaries in the central $\sim 0.4$ pc ($\sim 10''$) of the Milky Way nuclear star cluster. We measure the photometry of 563 stars using the Keck II NIRC2 imager ($K'$-band, $λ_0 = 2.124 \text{ } μ\text{m}$). We achieve a photometric uncertainty floor of $Δm_{K'} \sim 0.03$ ($\approx 3\%$), comparable to the highest precision achieved in other AO studies. Approximately half of our sample ($50 \pm 2 \%$) shows variability. $52 \pm 5\%$ of known early-type young stars and $43 \pm 4 \%$ of known late-type giants are variable. These variability fractions are higher than those of other young, massive star populations or late-type giants in globular clusters, and can be largely explained by two factors. First, our experiment time baseline is sensitive to long-term intrinsic stellar variability. Second, the proper motion of stars behind spatial inhomogeneities in the foreground extinction screen can lead to variability. We recover the two known Galactic center eclipsing binary systems: IRS 16SW and S4-258 (E60). We constrain the Galactic center eclipsing binary fraction of known early-type stars to be at least $2.4 \pm 1.7\%$. We find no evidence of an eclipsing binary among the young S-stars nor among the young stellar disk members. These results are consistent with the local OB eclipsing binary fraction. We identify a new periodic variable, S2-36, with a 39.43 day period. Further observations are necessary to determine the nature of this source. △ Less

Submitted 12 November, 2018; originally announced November 2018.

Comments: 69 pages, 28 figures, 12 tables. Accepted for publication in The Astrophysical Journal

Journal ref: The Astrophysical Journal, 871:103 (32pp), 2019 January 20

Improving Orbit Estimates for Incomplete Orbits with a New Approach to Priors -- with Applications from Black Holes to Planets

Authors: K. Kosmo O'Neil, G. D. Martinez, A. Hees, A. M. Ghez, T. Do, G. Witzel, Q. Konopacky, E. E. Becklin, D. S. Chu, J. Lu, K. Matthews, S. Sakai

Abstract: We propose a new approach to Bayesian prior probability distributions (priors) that can improve orbital solutions for low-phase-coverage orbits, where data cover less than approximately 40% of an orbit. In instances of low phase coverage such as with stellar orbits in the Galactic center or with directly-imaged exoplanets, data have low constraining power and thus priors can bias parameter estimat… ▽ More We propose a new approach to Bayesian prior probability distributions (priors) that can improve orbital solutions for low-phase-coverage orbits, where data cover less than approximately 40% of an orbit. In instances of low phase coverage such as with stellar orbits in the Galactic center or with directly-imaged exoplanets, data have low constraining power and thus priors can bias parameter estimates and produce under-estimated confidence intervals. Uniform priors, which are commonly assumed in orbit fitting, are notorious for this. We propose a new observable-based prior paradigm that is based on uniformity in observables. We compare performance of this observable-based prior and of commonly assumed uniform priors using Galactic center and directly-imaged exoplanet (HR 8799) data. The observable-based prior can reduce biases in model parameters by a factor of two and helps avoid under-estimation of confidence intervals for simulations with less than about 40% phase coverage. Above this threshold, orbital solutions for objects with sufficient phase coverage such as S0-2, a short-period star at the Galactic center with full phase coverage, are consistent with previously published results. Below this threshold, the observable-based prior limits prior influence in regions of prior dominance and increases data influence. Using the observable-based prior, HR 8799 orbital analyses favor lower eccentricity orbits and provide stronger evidence that the four planets have a consistent inclination around 30 degrees to within 1-sigma. This analysis also allows for the possibility of coplanarity. We present metrics to quantify improvements in orbital estimates with different priors so that observable-based prior frameworks can be tested and implemented for other low-phase-coverage orbits. △ Less

Submitted 10 June, 2019; v1 submitted 14 September, 2018; originally announced September 2018.

Comments: Published in AJ. 23 pages, 14 figures. Monte Carlo chains are available in the published article, or are available upon request

Journal ref: AJ, 158, 4 (2019)

Multiwavelength Light Curves of Two Remarkable Sagittarius A* Flares

Authors: G. G. Fazio, J. L. Hora, G. Witzel, S. P. Willner, M. L. N. Ashby, F. Baganoff, E. Becklin, S. Carey, D. Haggard, C. Gammie, A. Ghez, M. A. Gurwell, J. Ingalls, D. Marrone, M. R. Morris, H. A. Smith

Abstract: Sgr A*, the supermassive black hole (SMBH) at the center of our Milky Way Galaxy, is known to be a variable source of X-ray, near-infrared (NIR), and submillimeter (submm) radiation and therefore a prime candidate to study the electromagnetic radiation generated by mass accretion flow onto a black hole and/or a related jet. Disentangling the power source and emission mechanisms of this variability… ▽ More Sgr A*, the supermassive black hole (SMBH) at the center of our Milky Way Galaxy, is known to be a variable source of X-ray, near-infrared (NIR), and submillimeter (submm) radiation and therefore a prime candidate to study the electromagnetic radiation generated by mass accretion flow onto a black hole and/or a related jet. Disentangling the power source and emission mechanisms of this variability is a central challenge to our understanding of accretion flows around SMBHs. Simultaneous multiwavelength observations of the flux variations and their time correlations can play an important role in obtaining a better understanding of possible emission mechanisms and their origin. This paper presents observations of two flares that both apparently violate the previously established patterns in the relative timing of submm/NIR/X-ray flares from Sgr A*. One of these events provides the first evidence of coeval structure between NIR and submm flux increases, while the second event is the first example of the sequence of submm/X-ray/NIR flux increases all occurring within ~1 hr. Each of these two events appears to upend assumptions that have been the basis of some analytic models of flaring in Sgr A*. However, it cannot be ruled out that these events, even though unusual, were just coincidental. These observations demonstrate that we do not fully understand the origin of the multiwavelength variability of Sgr A*, and show that there is a continued and important need for long-term, coordinated, and precise multiwavelength observations of Sgr A* to characterize the full range of variability behavior. △ Less

Submitted 19 July, 2018; originally announced July 2018.

Comments: 9 pages, 3 figures, to appear in the Astrophysical Journal

Variability Timescale and Spectral Index of Sgr A* in the Near Infrared: Approximate Bayesian Computation Analysis of the Variability of the Closest Supermassive Black Hole

Authors: G. Witzel, G. Martinez, J. Hora, S. P. Willner, M. R. Morris, C. Gammie, E. E. Becklin, M. L. N. Ashby, F. Baganoff, S. Carey, T. Do, G. G. Fazio, A. Ghez, W. J. Glaccum, D. Haggard, R. Herrero-Illana, J. Ingalls, R. Narayan, H. A. Smith

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 present an analysis of the most comprehensive NIR variability dataset of Sgr A* to date: eight 24-hour epochs of continuous monitoring of Sgr A* at 4.5 $μ$m with the IRAC instrument on the Spitzer Space Telescope, 93 epochs of 2.18 $μ$m data from Na… ▽ 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 present an analysis of the most comprehensive NIR variability dataset of Sgr A* to date: eight 24-hour epochs of continuous monitoring of Sgr A* at 4.5 $μ$m with the IRAC instrument on the Spitzer Space Telescope, 93 epochs of 2.18 $μ$m data from Naos Conica at the Very Large Telescope, and 30 epochs of 2.12 $μ$m data from the NIRC2 camera at the Keck Observatory, in total 94,929 measurements. A new approximate Bayesian computation method for fitting the first-order structure function extracts information beyond current Fast Fourier Transformation (FFT) methods of power spectral density (PSD) estimation. With a combined fit of the data of all three observatories, the characteristic coherence timescale of Sgr A* is $τ_{b} = 243^{+82}_{-57}$ minutes ($90\%$ credible interval). The PSD has no detectable features on timescales down to 8.5 minutes ($95\%$ credible level), which is the ISCO orbital frequency for a dimensionless spin parameter $a = 0.92$. One light curve measured simultaneously at 2.12 and 4.5 $μ$m during a low flux-density phase gave a spectral index $α_s = 1.6 \pm 0.1$ ($F_ν\propto ν^{-α_s}$). This value implies that the Sgr A* NIR color becomes bluer during higher flux-density phases. The probability densities of flux densities of the combined datasets are best fit by log-normal distributions. Based on these distributions, the Sgr A* spectral energy distribution is consistent with synchrotron radiation from a non-thermal electron population from below 20 GHz through the NIR. △ Less

Submitted 5 September, 2018; v1 submitted 1 June, 2018; originally announced June 2018.

Comments: Accepted for publication in ApJ on May 30, 2018. A machine readable version of the light curve data is included in the journal's online publication. Version 2 includes proof corrections

Off-axis point spread function characterisation in laser-guide star adaptive optics systems

Authors: O. Beltramo-Martin, C. M. Correia, E. Mieda, B. Neichel, T. Fusco, G. Witzel, J. Lu, J. -P. Véran

Abstract: Adaptive optics (AO) restore the angular resolution of ground-based telescopes, but at the cost of delivering a time- and space-varying point spread function (PSF) with a complex shape. PSF knowledge is crucial for breaking existing limits on the measured accuracy of photometry and astrometry in science observations. In this paper, we concentrate our analyses on anisoplanatism signature only onto… ▽ More Adaptive optics (AO) restore the angular resolution of ground-based telescopes, but at the cost of delivering a time- and space-varying point spread function (PSF) with a complex shape. PSF knowledge is crucial for breaking existing limits on the measured accuracy of photometry and astrometry in science observations. In this paper, we concentrate our analyses on anisoplanatism signature only onto PSF: for large-field observations (20") with single- conjugated AO, PSFs are strongly elongated due to anisoplanatism that manifests itself as three different terms for Laser-guide star (LGS) systems: angular, focal and tilt. We propose a generalized model that relies on a point-wise decomposition of the phase and encompasses the non-stationarity of LGS systems. We demonstrate it is more accurate and less computationally demanding than existing models: it agrees with end-to-end physical-optics simulations to within 0.1% of PSF measurables, such as Strehl-ratio, FWHM and fraction of variance unexplained.Secondly, we study off-axis PSF modelling is with respect to $C_n^2(h)$ profile (heights and fractional weights). For 10m class telescope, PSF morphology is estimated at 1% level as long as we model the atmosphere with at least 7 layers whose heights and weights are known respectively with 200m and 10% precision. As a verification test we used the Canada's NRC-Herzberg HeNOS testbed data, featuring four lasers. We highlight capability of retrieving off-axis PSF characteristics within 10% of fraction of variance unexplained, which complies with the expected range from the sensitivity analysis. Our new off-axis PSF modelling method lays the ground-work for testing on-sky in the near future. △ Less

Submitted 14 April, 2018; originally announced April 2018.

Confusing binaries: the role of stellar binaries in biasing disk properties in the Galactic Center

Authors: Smadar Naoz, Andrea M. Ghez, Aurelien Hees, Tuan Do, Gunther Witzel, Jessica R. Lu

Abstract: The population of young stars near the Supermassive black hole (SMBH) in the Galactic Center (GC) has presented an unexpected challenge to theories of star formation. Kinematics measurements of these stars have revealed a stellar disk structure (with an apparent 20% disk membership) that has provided important clues to the origin of these mysterious young stars. However many of the apparent disk p… ▽ More The population of young stars near the Supermassive black hole (SMBH) in the Galactic Center (GC) has presented an unexpected challenge to theories of star formation. Kinematics measurements of these stars have revealed a stellar disk structure (with an apparent 20% disk membership) that has provided important clues to the origin of these mysterious young stars. However many of the apparent disk properties are difficult to explain, including the low disk membership fraction and the high eccentricities, given the youth of this population. Thus far, all efforts to derive the properties of this disk have made the simplifying assumption that stars at the GC are single stars. Nevertheless, stellar binaries are prevalent in our Galaxy, and recent investigations suggested that they may also be abundant in the Galactic Center. Here we show that binaries in the disk can largely alter the apparent orbital properties of the disk. The motion of binary members around each other adds a velocity component, which can be comparable to the magnitude of the velocity around the SMBH in the GC. Thus neglecting the contribution of binaries can significantly vary the inferred stars' orbital properties. While the disk orientation is unaffected the apparent disk's 2D width is increased to about 11.2deg, similar to the observed width. For a population of stars orbiting the SMBH with zero eccentricity, unaccounted for binaries will create a wide apparent eccentricity distribution with an average of 0.23.This is consistent with the observed average eccentricity of the stars' in the disk. We suggest that this high eccentricity value, which poses a theoretical challenge, may be an artifact of binary stars. Finally our results suggest that the actual disk membership might be significantly higher than the one inferred by observations that ignore the contribution of binaries, alleviating another theoretical challenge. △ Less

Submitted 11 January, 2018; originally announced January 2018.

Comments: Accepted to ApJ Lett, 5 figures, 8 pages

Investigating the Binarity of S0-2: Implications for its Origins and Robustness as a Probe of the Laws of Gravity around a Supermassive Black Hole

Authors: Devin S. Chu, Tuan Do, Aurelien Hees, Andrea Ghez, Smadar Naoz, Gunther Witzel, Shoko Sakai, Samantha Chappell, Abhimat K. Gautam, Jessica R. Lu, Keith Matthews

Abstract: The star S0-2, which orbits the supermassive black hole (SMBH) in our Galaxy with a period of 16 years, provides the strongest constraint on both the mass of the SMBH and the distance to the Galactic center. S0-2 will soon provide the first measurement of relativistic effects near a SMBH. We report the first limits on the binarity of S0-2 from radial velocity monitoring, which has implications for… ▽ More The star S0-2, which orbits the supermassive black hole (SMBH) in our Galaxy with a period of 16 years, provides the strongest constraint on both the mass of the SMBH and the distance to the Galactic center. S0-2 will soon provide the first measurement of relativistic effects near a SMBH. We report the first limits on the binarity of S0-2 from radial velocity monitoring, which has implications for both understanding its origin and robustness as a probe of the central gravitational field. With 87 radial velocity measurements, which include 12 new observations presented, we have the data set to look for radial velocity variations from S0-2's orbital model. Using a Lomb-Scargle analysis and orbit fitting for potential binaries, we detect no radial velocity variation beyond S0-2's orbital motion and do not find any significant periodic signal. The lack of a binary companion does not currently distinguish between different formation scenarios for S0-2. The upper limit on the mass of a companion star ($M_{\text{comp}}$) still allowed by our results has a median upper limit of $M_{\text{comp}}$ $\sin i \leq$ 1.6 M$_{\odot}$ for periods between 1 and 150 days, the longest period to avoid tidal break up of the binary. We also investigate the impact of the remaining allowed binary system on the measurement of the relativistic redshift at S0-2's closest approach in 2018. While binary star systems are important to consider for this experiment, we find plausible binaries for S0-2 will not alter a 5$σ$ detection of the relativistic redshift. △ Less

Submitted 18 December, 2017; v1 submitted 14 September, 2017; originally announced September 2017.

Comments: 12 pages, 7 figures, accepted to ApJ

The Post-Periapse Evolution of Galactic Center Source G1: The second case of a resolved tidal interaction with a supermassive black hole

Authors: Gunther Witzel, Breann. N. Sitarski, Andrea M. Ghez, Mark R. Morris, Aurelien Hees, Tuan Do, Jessica R. Lu, Smadar Naoz, Anna Boehle, Greg Martinez, Samantha Chappell, Rainer Schödel, Leo Meyer, Sylvana Yelda, Ereic E. Becklin, Keith Matthews

Abstract: We present new Adaptive Optics (AO) imaging and spectroscopic measurements of Galactic Center source G1 from W. M. Keck Observatory. Our goal is to understand its nature and relationship to G2, which is the first example of a spatially-resolved object interacting with the supermassive black hole (SMBH). Both objects have been monitored with AO for the past decade (2003 - 2014) and are comparativel… ▽ More We present new Adaptive Optics (AO) imaging and spectroscopic measurements of Galactic Center source G1 from W. M. Keck Observatory. Our goal is to understand its nature and relationship to G2, which is the first example of a spatially-resolved object interacting with the supermassive black hole (SMBH). Both objects have been monitored with AO for the past decade (2003 - 2014) and are comparatively close to the black hole ($a_{\rm{min}} \sim$200-300 AU) on very eccentric orbits ($e_{\rm{G1}}\sim$0.99; $e_{\rm{G2}}\sim$0.96). While G2 has been tracked before and during periapse passage ($T_{0} \sim$ 2014.2), G1 has been followed since soon after emerging from periapse ($T_{0} \sim$ 2001.3). Our observations of G1 double the previously reported observational time baseline, which improves its orbital parameter determinations. G1's orbital trajectory appears to be in the same plane as that of G2, but with a significantly different argument of periapse ($Δω$ = 21$\pm$4 degrees). This suggests that G1 is an independent object and not part of a gas stream containing G2 as has been proposed. Furthermore, we show for the first time that: (1) G1 is extended in the epochs closest to periapse along the direction of orbital motion and (2) G1 becomes significantly smaller over time, (450 AU in 2004 to less than 170 AU in 2009). Based on these observations, G1 appears to be the second example of an object tidally interacting with a SMBH. G1's existence 14 years after periapse, along with its compactness in epochs further from the time of periapse, suggest that this source is stellar in nature. △ Less

Submitted 7 July, 2017; originally announced July 2017.

Comments: submitted to ApJ

Testing the gravitational theory with short-period stars around our Galactic Center

Authors: A. Hees, A. M. Ghez, T. Do, J. R. Lu, M. R. Morris, E. E. Becklin, G. Witzel, A. Boehle, S. Chappell, Z. Chen, D. Chu, A. Ciurlo, A. Dehghanfar, E. Gallego-Cano, A. Gautam, S. Jia, K. Kosmo, G. D. Martinez, K. Matthews, S. Naoz, S. Sakai, R. Schödel

Abstract: Motion of short-period stars orbiting the supermassive black hole in our Galactic Center has been monitored for more than 20 years. These observations are currently offering a new way to test the gravitational theory in an unexplored regime: in a strong gravitational field, around a supermassive black hole. In this proceeding, we present three results: (i) a constraint on a hypothetical fifth forc… ▽ More Motion of short-period stars orbiting the supermassive black hole in our Galactic Center has been monitored for more than 20 years. These observations are currently offering a new way to test the gravitational theory in an unexplored regime: in a strong gravitational field, around a supermassive black hole. In this proceeding, we present three results: (i) a constraint on a hypothetical fifth force obtained by using 19 years of observations of the two best measured short-period stars S0-2 and S0-38 ; (ii) an upper limit on the secular advance of the argument of the periastron for the star S0-2 ; (iii) a sensitivity analysis showing that the relativistic redshift of S0-2 will be measured after its closest approach to the black hole in 2018. △ Less

Submitted 30 May, 2017; originally announced May 2017.

Comments: 4 pages, 2 figures, proceedings of the 52nd Rencontres de Moriond, Gravitation Session

Testing General Relativity with stellar orbits around the supermassive black hole in our Galactic center

Authors: A. Hees, T. Do, A. M. Ghez, G. D. Martinez, S. Naoz, E. E. Becklin, A. Boehle, S. Chappell, D. Chu, A. Dehghanfar, K. Kosmo, J. R. Lu, K. Matthews, M. R. Morris, S. Sakai, R. Schödel, G. Witzel

Abstract: In this Letter, we demonstrate that short-period stars orbiting around the supermassive black hole in our Galactic Center can successfully be used to probe the gravitational theory in a strong regime. We use 19 years of observations of the two best measured short-period stars orbiting our Galactic Center to constrain a hypothetical fifth force that arises in various scenarios motivated by the deve… ▽ More In this Letter, we demonstrate that short-period stars orbiting around the supermassive black hole in our Galactic Center can successfully be used to probe the gravitational theory in a strong regime. We use 19 years of observations of the two best measured short-period stars orbiting our Galactic Center to constrain a hypothetical fifth force that arises in various scenarios motivated by the development of a unification theory or in some models of dark matter and dark energy. No deviation from General Relativity is reported and the fifth force strength is restricted to an upper 95% confidence limit of $\left|α\right| < 0.016$ at a length scale of $λ=$ 150 astronomical units. We also derive a 95% confidence upper limit on a linear drift of the argument of periastron of the short-period star S0-2 of $\left|\dot ω_\textrm{S0-2} \right|< 1.6 \times 10^{-3}$ rad/yr, which can be used to constrain various gravitational and astrophysical theories. This analysis provides the first fully self-consistent test of the gravitational theory using orbital dynamic in a strong gravitational regime, that of a supermassive black hole. A sensitivity analysis for future measurements is also presented. △ Less

Submitted 22 May, 2017; originally announced May 2017.

Comments: 8 pages, 3 figures, version accepted for publication in Phys. Rev. Letters. Original title "Seeking a fifth force with..." changed upon request of the PRL editors

Journal ref: Phys. Rev. Letters 118, 211101, 2017

Constraining the Variability and Binary Fraction of Galactic Center Young Stars

Authors: Abhimat K. Gautam, Tuan Do, Andrea M. Ghez, Jessica R. Lu, Mark R. Morris, Shoko Sakai, Gunther Witzel, Breann N. Sitarski, Samantha Chappell

Abstract: We present constraints on the variability and binarity of young stars in the central 10 arcseconds (~0.4 pc) of the Milky Way Galactic Center (GC) using Keck Adaptive Optics data over a 12 year baseline. Given our experiment's photometric uncertainties, at least 36% of our sample's known early-type stars are variable. We identified eclipsing binary systems by searching for periodic variability. In… ▽ More We present constraints on the variability and binarity of young stars in the central 10 arcseconds (~0.4 pc) of the Milky Way Galactic Center (GC) using Keck Adaptive Optics data over a 12 year baseline. Given our experiment's photometric uncertainties, at least 36% of our sample's known early-type stars are variable. We identified eclipsing binary systems by searching for periodic variability. In our sample of spectroscopically confirmed and likely early-type stars, we detected the two previously discovered GC eclipsing binary systems. We derived the likely binary fraction of main sequence, early-type stars at the GC via Monte Carlo simulations of eclipsing binary systems, and find that it is at least 32% with 90% confidence. △ Less

Submitted 7 October, 2016; originally announced October 2016.

Comments: Accepted for publication in Proceedings of IAU Symposium 322: The Multi-Messenger Astrophysics of the Galactic Centre, 2 pages, 1 figure

An Improved Distance and Mass Estimate for Sgr A* from a Multistar Orbit Analysis

Authors: A. Boehle, A. M. Ghez, R. Schödel, L. Meyer, S. Yelda, S. Albers, G. D. Martinez, E. E. Becklin, T. Do, J. R. Lu, K. Matthews, M. R. Morris, B. Sitarski, G. Witzel

Abstract: We present new, more precise measurements of the mass and distance of our Galaxy's central supermassive black hole, Sgr A*. These results stem from a new analysis that more than doubles the time baseline for astrometry of faint stars orbiting Sgr A*, combining two decades of speckle imaging and adaptive optics data. Specifically, we improve our analysis of the speckle images by using information a… ▽ More We present new, more precise measurements of the mass and distance of our Galaxy's central supermassive black hole, Sgr A*. These results stem from a new analysis that more than doubles the time baseline for astrometry of faint stars orbiting Sgr A*, combining two decades of speckle imaging and adaptive optics data. Specifically, we improve our analysis of the speckle images by using information about a star's orbit from the deep adaptive optics data (2005 - 2013) to inform the search for the star in the speckle years (1995 - 2005). When this new analysis technique is combined with the first complete re-reduction of Keck Galactic Center speckle images using speckle holography, we are able to track the short-period star S0-38 (K-band magnitude = 17, orbital period = 19 years) through the speckle years. We use the kinematic measurements from speckle holography and adaptive optics to estimate the orbits of S0-38 and S0-2 and thereby improve our constraints of the mass ($M_{bh}$) and distance ($R_o$) of Sgr A*: $M_{bh} = 4.02\pm0.16\pm0.04\times10^6~M_{\odot}$ and $7.86\pm0.14\pm0.04$ kpc. The uncertainties in $M_{bh}$ and $R_o$ as determined by the combined orbital fit of S0-2 and S0-38 are improved by a factor of 2 and 2.5, respectively, compared to an orbital fit of S0-2 alone and a factor of $\sim$2.5 compared to previous results from stellar orbits. This analysis also limits the extended dark mass within 0.01 pc to less than $0.13\times10^{6}~M_{\odot}$ at 99.7% confidence, a factor of 3 lower compared to prior work. △ Less

Submitted 19 July, 2016; originally announced July 2016.

Comments: 56 pages, 14 figures, accepted to ApJ

Using infrared/X-ray flare statistics to probe the emission regions near the event horizon of Sgr A*

Authors: Salome Dibi, Sera Markoff, Renaud Belmont, Julien Malzac, Joey Neilsen, Gunther Witzel

Abstract: The supermassive black hole at the centre of the Galaxy flares at least daily in the infrared (IR) and X-ray bands, yet the process driving these flares is still unknown. So far detailed analysis has only been performed on a few bright flares. In particular, the broadband spectral modelling suffers from a strong lack of simultaneous data. However, new monitoring campaigns now provide data on thous… ▽ More The supermassive black hole at the centre of the Galaxy flares at least daily in the infrared (IR) and X-ray bands, yet the process driving these flares is still unknown. So far detailed analysis has only been performed on a few bright flares. In particular, the broadband spectral modelling suffers from a strong lack of simultaneous data. However, new monitoring campaigns now provide data on thousands of flaring events, allowing a statistical analysis of the flare properties. In this paper, we investigate the X-ray and IR flux distributions of the flare events. Using a self-consistent calculation of the particle distribution, we model the statistical properties of the flares. Based on a previous work on single flares, we consider two families of models: pure synchrotron models and synchrotron self-Compton (SSC) models. We investigate the effect of fluctuations in some relevant parameters (e.g. acceleration properties, density, magnetic field) on the flux distributions. The distribution of these parameters is readily derived from the flux distributions observed at different wavelengths. In both scenarios, we find that fluctuations of the power injected in accelerated particles plays a major role. This must be distributed as a power-law (with different indices in each model). In the synchrotron dominated scenario, we derive the most extreme values of the acceleration power required to reproduce the brightest flares. In that model, the distribution of the acceleration slope fluctuations is constrained and in the SSC scenario we constrain the distributions of the correlated magnetic field and flow density variations. △ Less

Submitted 21 June, 2016; originally announced June 2016.

Comments: 9 pages, 3 tables, 6 figures, MNRAS, June 2016

Merging Binaries in the Galactic Center: The eccentric Kozai-Lidov mechanism with stellar evolution

Authors: Alexander P. Stephan, Smadar Naoz, Andrea M. Ghez, Gunther Witzel, Breann N. Sitarski, Tuan Do, Bence Kocsis

Abstract: Most, if not all, stars in the field are born in binary configurations or higher multiplicity systems. In dense stellar environment such as the Galactic Center (GC), many stars are expected to be in binary configurations as well. These binaries form hierarchical triple body systems, with the massive black hole (MBH) as the third, distant object. The stellar binaries are expected to undergo large a… ▽ More Most, if not all, stars in the field are born in binary configurations or higher multiplicity systems. In dense stellar environment such as the Galactic Center (GC), many stars are expected to be in binary configurations as well. These binaries form hierarchical triple body systems, with the massive black hole (MBH) as the third, distant object. The stellar binaries are expected to undergo large amplitude eccentricity and inclination oscillations via the so-called "eccentric Kozai-Lidov" (EKL) mechanism. These eccentricity excitations, combined with post main sequence stellar evolution, can drive the inner stellar binaries to merge. We study the mergers of stellar binaries in the inner 0.1 pc of the GC caused by gravitational perturbations due to the MBH. We run a large set of Monte Carlo simulations that include the secular evolution of the orbits, general relativistic precession, tides, and post-main-sequence stellar evolution. We find that about 13 % of the initial binary population will have merged after a few million years and about 29 % after a few billion years. These expected merged systems represent a new class of objects at the GC and we speculate that they are connected to G2-like objects and the young stellar population. △ Less

Submitted 25 May, 2016; v1 submitted 8 March, 2016; originally announced March 2016.

Comments: 12 pages, 8 figures, accepted by MNRAS and available through 'Advance Access'

The X-ray Flux Distribution of Sagittarius A* as Seen by Chandra

Authors: J. Neilsen, S. Markoff, M. A. Nowak, J. Dexter, G. Witzel, N. Barrière, Y. Li, F. K. Baganoff, N. Degenaar, P. C. Fragile, C. Gammie, A. Goldwurm, N. Grosso, D. Haggard

Abstract: We present a statistical analysis of the X-ray flux distribution of Sgr A* from the Chandra X-ray Observatory's 3 Ms Sgr A* X-ray Visionary Project (XVP) in 2012. Our analysis indicates that the observed X-ray flux distribution can be decomposed into a steady quiescent component, represented by a Poisson process with rate $Q=(5.24\pm0.08)\times10^{-3}$ cts s$^{-1},$ and a variable component, repre… ▽ More We present a statistical analysis of the X-ray flux distribution of Sgr A* from the Chandra X-ray Observatory's 3 Ms Sgr A* X-ray Visionary Project (XVP) in 2012. Our analysis indicates that the observed X-ray flux distribution can be decomposed into a steady quiescent component, represented by a Poisson process with rate $Q=(5.24\pm0.08)\times10^{-3}$ cts s$^{-1},$ and a variable component, represented by a power law process ($dN/dF\propto F^{-ξ},$ $ξ=1.92_{-0.02}^{+0.03}$). This slope matches our recently-reported distribution of flare luminosities. The variability may also be described by a log-normal process with a median unabsorbed 2-8 keV flux of $1.8^{+0.9}_{-0.6}\times10^{-14}$ erg s$^{-1}$ cm$^{-2}$ and a shape parameter $σ=2.4\pm0.2,$ but the power law provides a superior description of the data. In this decomposition of the flux distribution, all of the intrinsic X-ray variability of Sgr A* (spanning at least three orders of magnitude in flux) can be attributed to flaring activity, likely in the inner accretion flow. We confirm that at the faint end, the variable component contributes ~10% of the apparent quiescent flux, as previously indicated by our statistical analysis of X-ray flares in these Chandra observations. Our flux distribution provides a new and important observational constraint on theoretical models of Sgr A*, and we use simple radiation models to explore the extent to which a statistical comparison of the X-ray and infrared can provide insights into the physics of the X-ray emission mechanism. △ Less

Submitted 9 December, 2014; originally announced December 2014.

Comments: 13 pages, 7 figures, accepted for publication in ApJ. Comments welcome

Polarized light from Sgr A* in the near-infrared $K_{s}$-band

Authors: B. Shahzamanian, A. Eckart, M. Valencia-S., G. Witzel, M. Zamaninasab, N. Sabha, M. García-Marín, V. Karas, G. D. Karssen, A. Borkar, M. Dovuciak, D. Kunneriath, M. Bursa, R. Buchholz, J. Moultaka, C. Straubmeier

Abstract: We present a statistical analysis of polarized near-infrared (NIR) light from Sgr A*, the radio source associated with the supermassive black hole at the center of the Milky Way. The observations have been carried out using the adaptive optics instrument NACO at the VLT UT4 in the infrared $K_\mathrm{s}$-band from 2004 to 2012. Several polarized flux excursions were observed during these years. Li… ▽ More We present a statistical analysis of polarized near-infrared (NIR) light from Sgr A*, the radio source associated with the supermassive black hole at the center of the Milky Way. The observations have been carried out using the adaptive optics instrument NACO at the VLT UT4 in the infrared $K_\mathrm{s}$-band from 2004 to 2012. Several polarized flux excursions were observed during these years. Linear polarization at 2.2 $μm$, its statistics and time variation, can be used constrain the physical conditions of the accretion process onto this supermassive black hole. With an exponent of about 4 for the number density histogram of fluxes above 5~mJy, the distribution of polarized flux density is closely linked to the single state power-law distribution of the total $K_\mathrm{s}$-band flux densities reported earlier. We find typical polarization degrees of the order of 20\%$\pm$10\% and a preferred polarization angle of 13$^o$$\pm$15$^o$. Simulations show the uncertainties under a total flux density of $\sim 2\,{\rm mJy}$ are probably dominated by observational effects. At higher flux densities there are intrinsic variations of polarization degree and angle within rather well constrained ranges. Since the emission is most likely due to optically thin synchrotron radiation, this preferred polarization angle we find is very likely coupled to the intrinsic orientation of the Sgr A* system i.e. a disk or jet/wind scenario associated with the super massive black hole. If they are indeed linked to structural features of the source the data imply a rather stable geometry and accretion process for the Sgr A* system. △ Less

Submitted 31 October, 2014; originally announced November 2014.

Comments: 19 pages, 18 figures, 4 tables, submitted to Astronomy & Astrophysics

Journal ref: A&A 576, A20 (2015)

Detection of Galactic Center source G2 at 3.8 $μ$m during periapse passage

Authors: Gunther Witzel, Andrea M. Ghez, Mark R. Morris, Breann N. Sitarski, Anna Boehle, Smadar Naoz, Randall Campbell, Eric E. Becklin, Gabriela Canalizo, Samantha Chappell, Tuan Do, Jessica R. Lu, Keith Matthews, Leo Meyer, Alan Stockton, Peter Wizinowich, Sylvana Yelda

Abstract: We report new observations of the Galactic Center source G2 from the W. M. Keck Observatory. G2 is a dusty red object associated with gas that shows tidal interactions as it nears closest approach with the Galaxy's central black hole. Our observations, conducted as G2 passed through periapse, were designed to test the proposal that G2 is a 3 earth mass gas cloud. Such a cloud should be tidally dis… ▽ More We report new observations of the Galactic Center source G2 from the W. M. Keck Observatory. G2 is a dusty red object associated with gas that shows tidal interactions as it nears closest approach with the Galaxy's central black hole. Our observations, conducted as G2 passed through periapse, were designed to test the proposal that G2 is a 3 earth mass gas cloud. Such a cloud should be tidally disrupted during periapse passage. The data were obtained using the Keck II laser guide star adaptive optics system (LGSAO) and the facility near-infrared camera (NIRC2) through the K' [2.1 $μ$m] and L' [3.8 $μ$m] broadband filters. Several results emerge from these observations: 1) G2 has survived its closest approach to the black hole as a compact, unresolved source at L'; 2) G2's L' brightness measurements are consistent with those over the last decade; 3) G2's motion continues to be consistent with a Keplerian model. These results rule out G2 as a pure gas cloud and imply that G2 has a central star. This star has a luminosity of $\sim$30 $L_{\odot} $ and is surrounded by a large ($\sim$2.6 AU) optically thick dust shell. The differences between the L' and Br-$γ$ observations can be understood with a model in which L' and Br-$γ$ emission arises primarily from internal and external heating, respectively. We suggest that G2 is a binary star merger product and will ultimately appear similar to the B-stars that are tightly clustered around the black hole (the so-called S-star cluster). △ Less

Submitted 15 October, 2014; v1 submitted 7 October, 2014; originally announced October 2014.

Comments: Accepted by ApJ Letters, 2014 October 14

Spitzer/IRAC Observations of the Variability of Sgr A* and the Object G2 at 4.5 microns

Authors: J. L. Hora, G. Witzel, M. L. N. Ashby, E. E. Becklin, S. Carey, G. G. Fazio, A. Ghez, J. Ingalls, L. Meyer, M. R. Morris, H. A. Smith, S. P. Willner

Abstract: We present the first detection from the Spitzer Space Telescope of 4.5 micron variability from Sgr A*, the emitting source associated with the Milky Way's central black hole. The >23 hour continuous light curve was obtained with the IRAC instrument in 2013 December. The result characterizes the variability of Sgr A* prior to the closest approach of the G2 object, a putative infalling gas cloud tha… ▽ More We present the first detection from the Spitzer Space Telescope of 4.5 micron variability from Sgr A*, the emitting source associated with the Milky Way's central black hole. The >23 hour continuous light curve was obtained with the IRAC instrument in 2013 December. The result characterizes the variability of Sgr A* prior to the closest approach of the G2 object, a putative infalling gas cloud that orbits close to Sgr A*. The high stellar density at the location of Sgr A* produces a background of ~250 mJy at 4.5 microns in each pixel with a large pixel-to-pixel gradient, but the light curve for the highly variable Sgr A* source was successfully measured by modeling and removing the variations due to pointing wobble. The observed flux densities range from the noise level of ~0.7 mJy rms in a 6.4-s measurement to ~10 mJy. Emission was seen above the noise level ~34% of the time. The light curve characteristics, including the flux density distribution and structure function, are consistent with those previously derived at shorter infrared wavelengths. We see no evidence in the light curve for activity attributable to the G2 interaction at the observing epoch, ~100 days before the expected G2 periapsis passage. The IRAC light curve is more than a factor of two longer than any previous infrared observation, improving constraints on the timescale of the break in the power spectral distribution of Sgr A* flux densities. The data favor the longer of the two previously published values for the timescale. △ Less

Submitted 8 August, 2014; originally announced August 2014.

Comments: 13 pages, 10 figures, 2 tables, accepted for publication in the ApJ

The Intrinsic Two-Dimensional Size of Sagittarius A*

Authors: Geoffrey C. Bower, Sera Markoff, Andreas Brunthaler, Casey Law, Heino Falcke, Dipankar Maitra, M. Clavel, A. Goldwurm, M. R. Morris, Gunther Witzel, Leo Meyer, A. M. Ghez

Abstract: We report the detection of the two-dimensional structure of the radio source associated with the Galactic Center black hole, Sagittarius A*, obtained from Very Long Baseline Array (VLBA) observations at a wavelength of 7mm. The intrinsic source is modeled as an elliptical Gaussian with major axis size 35.4 x 12.6 R_S in position angle 95 deg East of North. This morphology can be interpreted in the… ▽ More We report the detection of the two-dimensional structure of the radio source associated with the Galactic Center black hole, Sagittarius A*, obtained from Very Long Baseline Array (VLBA) observations at a wavelength of 7mm. The intrinsic source is modeled as an elliptical Gaussian with major axis size 35.4 x 12.6 R_S in position angle 95 deg East of North. This morphology can be interpreted in the context of both jet and accretion disk models for the radio emission. There is supporting evidence in large angular-scale multi-wavelength observations for both source models for a preferred axis near 95 deg. We also place a maximum peak-to-peak change of 15% in the intrinsic major axis size over five different epochs. Three observations were triggered by detection of near infrared (NIR) flares and one was simultaneous with a large X-ray flare detected by NuSTAR. The absence of simultaneous and quasi-simultaneous flares indicates that not all high energy events produce variability at radio wavelengths. This supports the conclusion that NIR and X-ray flares are primarily due to electron excitation and not to an enhanced accretion rate onto the black hole. △ Less

Submitted 6 May, 2014; originally announced May 2014.

Comments: accepted for publication in ApJ

A formal method for identifying distinct states of variability in time-varying sources: SgrA* as an example

Authors: L. Meyer, G. Witzel, F. A. Longstaff, A. M. Ghez

Abstract: Continuously time variable sources are often characterized by their power spectral density and flux distribution. These quantities can undergo dramatic changes over time if the underlying physical processes change. However, some changes can be subtle and not distinguishable using standard statistical approaches. Here, we report a methodology that aims to identify distinct but similar states of tim… ▽ More Continuously time variable sources are often characterized by their power spectral density and flux distribution. These quantities can undergo dramatic changes over time if the underlying physical processes change. However, some changes can be subtle and not distinguishable using standard statistical approaches. Here, we report a methodology that aims to identify distinct but similar states of time variability. We apply this method to the Galactic supermassive black hole, where 2.2 um flux is observed from a source associated with SgrA*, and where two distinct states have recently been suggested. Our approach is taken from mathematical finance and works with conditional flux density distributions that depend on the previous flux value. The discrete, unobserved (hidden) state variable is modeled as a stochastic process and the transition probabilities are inferred from the flux density time series. Using the most comprehensive data set to date, in which all Keck and a majority of the publicly available VLT data have been merged, we show that SgrA* is sufficiently described by a single intrinsic state. However the observed flux densities exhibit two states: a noise-dominated and a source-dominated one. Our methodology reported here will prove extremely useful to assess the effects of the putative gas cloud G2 that is on its way toward the black hole and might create a new state of variability. △ Less

Submitted 20 March, 2014; originally announced March 2014.

Comments: Submitted to ApJ; 33 pages, 4 figures; comments welcome

The Keplerian orbit of G2

Authors: L. Meyer, A. M. Ghez, G. Witzel, T. Do, K. Phifer, B. N. Sitarski, M. R. Morris, A. Boehle, S. Yelda, J. R. Lu, E. Becklin

Abstract: We give an update of the observations and analysis of G2 - the gaseous red emission-line object that is on a very eccentric orbit around the Galaxy's central black hole and predicted to come within 2400 Rs in early 2014. During 2013, the laser guide star adaptive optics systems on the W. M. Keck I and II telescopes were used to obtain three epochs of spectroscopy and imaging at the highest spatial… ▽ More We give an update of the observations and analysis of G2 - the gaseous red emission-line object that is on a very eccentric orbit around the Galaxy's central black hole and predicted to come within 2400 Rs in early 2014. During 2013, the laser guide star adaptive optics systems on the W. M. Keck I and II telescopes were used to obtain three epochs of spectroscopy and imaging at the highest spatial resolution currently possible in the near-IR. The updated orbital solution derived from radial velocities in addition to Br-Gamma line astrometry is consistent with our earlier estimates. Strikingly, even ~6 months before pericenter passage there is no perceptible deviation from a Keplerian orbit. We furthermore show that a proposed "tail" of G2 is likely not associated with it but is rather an independent gas structure. We also show that G2 does not seem to be unique, since several red emission-line objects can be found in the central arcsecond. Taken together, it seems more likely that G2 is ultimately stellar in nature, although there is clearly gas associated with it. △ Less

Submitted 5 December, 2013; originally announced December 2013.

Comments: Proceedings of IAU Symposium #303, "The Galactic Center: Feeding and Feedback in a Normal Galactic Nucleus"; 2013 September 30 - October 4, Santa Fe New Mexico (USA)

Ks- and Lp-band polarimetry on stellar and bow-shock sources in the Galactic center

Authors: Rainer M. Buchholz, Gunther Witzel, Rainer Schoedel, Andreas Eckart

Abstract: Infrared observations of the Galactic center (GC) provide a unique opportunity to study stellar and bow-shock polarization effects in a dusty environment. The goals of this work are to present new Ks- and Lp-band polarimetry on an unprecedented number of sources in the central parsec of the GC, thereby expanding our previous results in the H- and Ks-bands. We use AO-assisted Ks- and Lp-band observ… ▽ More Infrared observations of the Galactic center (GC) provide a unique opportunity to study stellar and bow-shock polarization effects in a dusty environment. The goals of this work are to present new Ks- and Lp-band polarimetry on an unprecedented number of sources in the central parsec of the GC, thereby expanding our previous results in the H- and Ks-bands. We use AO-assisted Ks- and Lp-band observations, obtained at the ESO VLT. High precision photometry and the new polarimetric calibration method for NACO allow us to map the polarization in a region of 8" x 25" (Ks) resp. 26" x 28" (Lp). These are the first polarimetric observations of the GC in the Lp-band in 30 years, with vastly improved spatial resolution compared to previous results. This allows resolved polarimetry on bright bow-shock sources in this area for the first time at this wavelength. We find foreground polarization to be largely parallel to the Galactic plane (Ks-band: 6.1% at 20 degrees, Lp-band: 4.5% at 20 degrees, in good agreement with our previous findings and with older results. The previously described Lp-band excess in the foregound polarization towards the GC could be confirmed here for a much larger number of sources. The bow-shock sources contained in the FOV seem to show a different relation between the polarization in the observed wavelength bands than what was determined for the foreground. This points to the different relevant polarization mechanisms. The resolved polarization patterns of IRS 5 and 10W match the findings we presented earlier for IRS~1W. Additionally, intrinsic Lp-band polarization was measured for IRS 1W and 21, as well as for other, less prominent MIR-excess sources (IRS 2S, 2L, 5NE). The new data offer support for the presumed bow-shock nature of several of these sources (1W, 5, 5NE, 10W, 21) and for the model of bow-shock polarization presented in our last work. △ Less

Submitted 5 August, 2013; originally announced August 2013.

Comments: 19 pages, 18 figures

Keck Observations of the Galactic Center Source G2: Gas Cloud or Star?

Authors: K. Phifer, T. Do, L. Meyer, A. M. Ghez, G. Witzel, S. Yelda, A. Boehle, J. R. Lu, M. R. Morris, E. E. Becklin, K. Matthews

Abstract: We present new observations and analysis of G2 - the intriguing red emission-line object which is quickly approaching the Galaxy's central black hole. The observations were obtained with the laser guide star adaptive optics systems on the W. M. Keck I and II telescopes and include spectroscopy (R ~ 3600) centered on the Hydrogen Br-gamma line as well as K' (2.1 micrometer) and L' (3.8 micrometer)… ▽ More We present new observations and analysis of G2 - the intriguing red emission-line object which is quickly approaching the Galaxy's central black hole. The observations were obtained with the laser guide star adaptive optics systems on the W. M. Keck I and II telescopes and include spectroscopy (R ~ 3600) centered on the Hydrogen Br-gamma line as well as K' (2.1 micrometer) and L' (3.8 micrometer) imaging. Analysis of these observations shows the Br-gamma line emission has a positional offset from the L' continuum. This offset is likely due to background source confusion at L'. We therefore present the first orbital solution derived from Br-gamma line astrometry, which when coupled with radial velocity measurements, results in a later time of closest approach (2014.21 +/- 0.14), closer periastron (130 AU, 1900Rs), and higher eccentricity (0.9814 +/- 0.0060) compared to a solution using L' astrometry. The new orbit casts doubt on previous associations of G2 and a low surface brightness "tail". It is shown that G2 has no K' counterpart down to K' ~ 20 mag. G2's L' continuum and the Br-gamma line-emission is unresolved in almost all epochs; however it is marginally extended in our highest quality Br-gamma data set from 2006 and exhibits a clear velocity gradient at that time. While the observations altogether suggest that G2 has a gaseous component which is tidally interacting with the central black hole, there is likely a central star providing the self-gravity necessary to sustain the compact nature of this object. △ Less

Submitted 23 April, 2013; v1 submitted 18 April, 2013; originally announced April 2013.

Comments: 14 pages, 3 figures, submitted to ApJ. For higher resolution figures see http://www.astro.ucla.edu/~ghezgroup/gc/G2_2013.pdf

Source-intrinsic near-infrared properties of Sgr A*: Total intensity measurements

Authors: G. Witzel, A. Eckart, M. Bremer, M. Zamaninasab, B. Shahzamanian, M. Valencia-S., R. Schödel, V. Karas, R. Lenzen, N. Marchili, N. Sabha, M. Garcia-Marin, R. M. Buchholz, D. Kunneriath, C. Straubmeier

Abstract: We present a comprehensive data description for Ks-band measurements of Sgr A*. We characterize the statistical properties of the variability of Sgr A* in the near-infrared, which we find to be consistent with a single-state process forming a power-law distribution of the flux density. We discover a linear rms-flux relation for the flux-density range up to 12 mJy on a timescale of 24 minutes. This… ▽ More We present a comprehensive data description for Ks-band measurements of Sgr A*. We characterize the statistical properties of the variability of Sgr A* in the near-infrared, which we find to be consistent with a single-state process forming a power-law distribution of the flux density. We discover a linear rms-flux relation for the flux-density range up to 12 mJy on a timescale of 24 minutes. This and the power-law flux density distribution implies a phenomenological, formally non-linear statistical variability model with which we can simulate the observed variability and extrapolate its behavior to higher flux levels and longer timescales. We present reasons why data with our cadence cannot be used to decide on the question whether the power spectral density of the underlying random process shows more structure at timescales between 25 min and 100 min compared to what is expected from a red noise random process. △ Less

Submitted 29 August, 2012; originally announced August 2012.

Comments: Accepted to ApJS, August 27, 2012

Near-infrared proper motions and spectroscopy of infrared excess sources at the Galactic Center

Authors: A. Eckart, K. Muzic, S. Yazici, B. Shahzamanian, G. Witzel, N. Sabha, L. Moser, M. Garcia-Marin, M. Valencia-S., B. Jalali, M. Bremer, C. Straubmeier, D. Kunneriath, J. Moultaka

Abstract: There are a number of faint compact infrared excess sources in the central stellar cluster of the Milky Way. Their nature and origin is unclear. In addition to several isolated objects of this kind we find a small but dense cluster of co-moving sources (IRS13N) about 3" west of SgrA* just 0.5" north of the bright IRS13E cluster of WR and O-type stars. Based on their color and brightness, there are… ▽ More There are a number of faint compact infrared excess sources in the central stellar cluster of the Milky Way. Their nature and origin is unclear. In addition to several isolated objects of this kind we find a small but dense cluster of co-moving sources (IRS13N) about 3" west of SgrA* just 0.5" north of the bright IRS13E cluster of WR and O-type stars. Based on their color and brightness, there are two main possibilities: (1) they may be dust embedded stars older than few Myr, or (2) extremely young, dusty stars with ages less than 1Myr. We present fist H- and Ks-band identifications or proper motions of the IRS13N members, the high velocity dusty S-cluster object (DSO), and other infrared excess sources in the central field. We also present results of NIR H- and Ks-band ESO-SINFONI integral field spectroscopy of ISR13N. We show that within the uncertainties, the proper motions of the IRS13N sources in Ks- and L'-band are identical. This indicates that the bright L'-band IRS13N sources are indeed dust enshrouded stars rather than core-less dust clouds. The proper motions show that the IRS13N sources are not strongly gravitationally bound to each other implying that they have been formed recently. We also present a first H- and Ks-band identification as well as proper motions and HKsL'-colors of a fast moving DSO which was recently found in the cluster of high speed S-stars that surround the super-massive black hole Sagittarius A* (SgrA*). Most of the compact L'-band excess emission sources have a compact H- or Ks-band counterpart and therefore are likely stars with dust shells or disks. Our new results and orbital analysis from our previous work favor the hypothesis that the infrared excess IRS13N members and other dusty sources close to SgrA* are very young dusty stars and that star formation at the GC is a continuously ongoing process. △ Less

Submitted 22 December, 2012; v1 submitted 9 August, 2012; originally announced August 2012.

Comments: 20 pages, 18 figures, 4 tables plus appendix with 16 figures and 3 tables accepted by A&A