Search

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

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\%.

We present an analysis of the long-term performance of the W. M. Keck observatory laser guide star adaptive optics (LGS-AO) system and explore factors that influence the overall AO performance most strongly. Astronomical surveys can take years or decades to finish, so it is worthwhile to characterize the AO performance on such timescales in order to better understand future results. The Keck telescopes have two of the longest-running LGS-AO systems in use today, and as such they represent an excellent test-bed for processing large amounts of AO data. We use a Keck-II near infrared camera 2 (NIRC2) LGSAO surve of the Galactic Center (GC) from 2005 to 2019 for our analysis, combining image metrics with AO telemetry files, multiaperture scintillation sense/differential imaging motion monitor turbulence profiles, seeing information, weather data, and temperature readings in a compiled dataset to highlight areas of potential performance improvement. We find that image quality trends downward over time, despite multiple improvements made to Keck-II and its AO system, resulting in a 9 mas increase in the average full width at half maximum (FWHM) and a 3% decrease in the average Strehl ratio over the course of the survey. Image quality also trends upward with ambient temperature, possibly indicating the presence of uncorrected turbulence in the beam path. Using nine basic features from our dataset, we train a simple machine learning (ML) algorithm to predict the delivered image quality of NIRC2 given current atmospheric conditions, which could eventually be used for real-time observation planning and exposure time adjustments. A random forest algorithm trained on this data can predict the Strehl ratio of an image to within 18% and the FWHM to within 7%, which is a solid baseline for future applications involving more advanced ML techniques. The assembled dataset and coding tools are released to the public as a resource for testing new predictive control and point spread function-reconstruction algorithms.

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

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

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., Comment: 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

The InfraRed Imaging Spectrograph (IRIS) is a first-light instrument for the Thirty Meter Telescope (TMT) that will be used to sample the corrected adaptive optics field by the Narrow-Field Infrared Adaptive Optics System (NFIRAOS) with a near-infrared (0.8 - 2.4 µm) imaging camera and integral field spectrograph. To better understand IRIS science specifications we use the IRIS data simulator to characterize relative photometric precision and accuracy across the IRIS imaging camera 34”x34” field of view. Because the Point Spread Function (PSF) varies due to the effects of anisoplanatism, we use the Anisoplanatic and Instrumental Reconstruction of Off-axis PSFs for AO (AIROPA) software package to conduct photometric measurements on simulated frames using PSF-fitting as the PSF varies in single-source, binary, and crowded field use cases. We report photometric performance of the imaging camera as a function of the instrumental noise properties including dark current and read noise. Using the same methods, we conduct comparisons of photometric performance with reconstructed PSFs, in order to test the veracity of the current PSF-Reconstruction algorithms for IRIS/TMT., Ground-based and Airborne Instrumentation for Astronomy VIII, December 14-18, 2020, Online Only, United States, Series: Proceedings of SPIE

We present the status and plans for the Keck All sky Precision Adaptive optics (KAPA) program. KAPA includes four key science programs, an upgrade to the Keck I laser guide star (LGS) adaptive optics (AO) facility to improve image quality and sky coverage, AO telemetry based point spread function (PSF) estimates for all science exposures, and an educational component focused on broadening the participation of women and underrepresented groups in instrumentation. For the purpose of this conference we will focus on the AO facility upgrade which includes implementation of a new laser, wavefront sensor and real-time controller to support laser tomography, the laser tomography system itself, and modifications to an existing near-infrared tip-tilt sensor to support multiple natural guide star (NGS) and focus measurements.

A small cluster of massive stars residing in the Galactic center, collectively known as IRS13E, is of special interest due to its close proximity to Sgr A* and the possibility that an embedded intermediate-mass black hole (IMBH) binds its member stars. It has been suggested that colliding winds from two member stars, both classified as Wolf-Rayet type, are responsible for the observed X-ray, infrared and radio emission from IRS13E. We have conducted an in-depth study of the X-ray spatial, temporal and spectral properties of IRS13E, based on 5.6 Ms of ultra-deep Chandra observations obtained over 20 years. These X-ray observations show no significant evidence for source variability. We have also explored the kinematics of the cluster members, using Keck near-infrared imaging and spectroscopic data on a 14-yr baseline that considerably improve the accuracy of stars' proper motions. The observations are interpreted using 3-dimensional hydrodynamical simulations of colliding winds tailored to match the physical conditions of IRS13E, leading us to conclude that the observed X-ray spectrum and morphology can be well explained by the colliding wind scenario, in the meantime offering no support for the presence of a putative IMBH. An IMBH more massive than a few $10^3{\rm~M_\odot}$ is also strongly disfavored by the stellar kinematics., 24 pages, 9 figures, Accepted for publication in Astrophysical Journal

The hierarchical nature of galaxy formation suggests that a supermassive black hole binary could exist in our galactic center. We propose a new approach to constraining the possible orbital configuration of such a binary companion to the galactic center black hole Sgr A* through the measurement of stellar orbits. Focusing on the star S0-2, we show that requiring its orbital stability in the presence of a companion to Sgr A* yields stringent constraints on the possible configurations of such a companion. Furthermore, we show that precise measurements of {\it time variations} in the orbital parameters of S0-2 could yield stronger constraints. Using existing data on S0-2 we derive upper limits on the binary black hole separation as a function of the companion mass. For the case of a circular orbit, we can rule out a 10^5 M_sun companion with a semimajor axis greater than 170 astronomical units or 0.8 mpc. This is already more stringent than bounds obtained from studies of the proper motion of Sgr A*. Including other stars orbiting the galactic center should yield stronger constraints that could help uncover the presence of a companion to Sgr A*. We show that a companion can also affect the accretion process, resulting in a variability which may be consistent with the measured infrared flaring timescales and amplitudes. Finally, if such a companion exists, it will emit gravitational wave radiation, potentially detectable with LISA., Comment: 9 pages, 6 figures, accepted to ApJ-Lett

We present $L$'-band imaging of the PDS 70 planetary system with Keck/NIRC2 using the new infrared pyramid wavefront sensor. We detected both PDS 70 b and c in our images, as well as the front rim of the circumstellar disk. After subtracting off a model of the disk, we measured the astrometry and photometry of both planets. Placing priors based on the dynamics of the system, we estimated PDS 70 b to have a semi-major axis of $20^{+3}_{-4}$~au and PDS 70 c to have a semi-major axis of $34^{+12}_{-6}$~au (95\% credible interval). We fit the spectral energy distribution (SED) of both planets. For PDS 70 b, we were able to place better constraints on the red half of its SED than previous studies and inferred the radius of the photosphere to be 2-3~$R_{Jup}$. The SED of PDS 70 c is less well constrained, with a range of total luminosities spanning an order of magnitude. With our inferred radii and luminosities, we used evolutionary models of accreting protoplanets to derive a mass of PDS 70 b between 2 and 4 $M_{\textrm{Jup}}$ and a mean mass accretion rate between $3 \times 10^{-7}$ and $8 \times 10^{-7}~M_{\textrm{Jup}}/\textrm{yr}$. For PDS 70 c, we computed a mass between 1 and 3 $M_{\textrm{Jup}}$ and mean mass accretion rate between $1 \times 10^{-7}$ and $5 \times~10^{-7} M_{\textrm{Jup}}/\textrm{yr}$. The mass accretion rates imply dust accretion timescales short enough to hide strong molecular absorption features in both planets' SEDs., Comment: 20 pages, 5 figures, Accepted to AJ. Updated author list from original version. Fixed equation typo

Searching for space-time variations of the constants of Nature is a promising way to search for new physics beyond General Relativity and the standard model motivated by unification theories and models of dark matter and dark energy. We propose a new way to search for a variation of the fine-structure constant using measurements of late-type evolved giant stars from the S-star cluster orbiting the supermassive black hole in our Galactic Center. A measurement of the difference between distinct absorption lines (with different sensitivity to the fine structure constant) from a star leads to a direct estimate of a variation of the fine structure constant between the star's location and Earth. Using spectroscopic measurements of 5 stars, we obtain a constraint on the relative variation of the fine structure constant below $10^{-5}$. This is the first time a varying constant of Nature is searched for around a black hole and in a high gravitational potential. This analysis shows new ways the monitoring of stars in the Galactic Center can be used to probe fundamental physics., Comment: 7 pages + 10 pages appendix, 3 figures, version accepted for publication

The electromagnetic counterpart to the Galactic center supermassive black hole, Sgr A∗, has been observed in the near-infrared for over 20 yr 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′ band, Sgr A∗ reached flux levels of ∼6 mJy, twice the level of the previously observed peak flux from >13,000 measurements over 130 nights with the Very Large Telescope and Keck Telescopes. We also observe a factor of 75 change in flux over a 2 hr time span with no obvious color changes between 1.6 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 of observing flux levels that are similar to all four 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.© 2019. The American Astronomical Society. All rights reserved.., We thank the anonymous referee for helpful comments. We thank the staff and astronomers at Keck Observatory especially Jim Lyke, Randy Campbell, Sherry Yeh, Greg Doppmann, Cynthia Wilburn, Terry Stickel, and Alan Hatakeyama. Support for this work was provided by NSF AAG grant AST-1412615, the W.M. Keck Foundation, the HeisingSimons Foundation, the Gordon and Betty Moore Foundation, the Levine-Leichtman Family Foundation, the Preston Family Graduate Fellowship (held by A.G.), and the UCLA Galactic Center Star Society. R. S. has received funding from the European Research Council under the European Union's Seventh Framework Programme (FP7/2007-2013)/ERC grant agreement No. 614922. R.S. acknowledges financial support from the State Agency for Research of the Spanish MCIU through the >Center of Excellence Severo Ochoa> award for the Instituto de Astrofisica de Andalucia (SEV-2017-0709). This research was based on data products from the Galactic Center Orbits Initiative, which is hosted at UCLA and is a key science program of the Galactic Center Collaboration. The W.M. Keck Observatory is operated as a scientific partnership among the California Institute of Technology, the University of California, and the National Aeronautics and Space Administration. The Observatory was made possible by the generous financial support of the W.M. Keck Foundation. The authors wish to recognize that the summit of Maunakea has always held a very significant cultural role for the indigenous Hawaiian community. We are most fortunate to have the opportunity to observe from this mountain.

The central 0.1parsecs of the Milky Way host a supermassive black hole identified with the position of the radio and infrared source Sagittarius A* (refs. 1,2), a cluster of young, massive stars (the S stars3) and various gaseous features4,5. Recently, two unusual objects have been found to be closely orbiting Sagittarius A*: the so-called G sources, G1 and G2. These objects are unresolved (having a size of the order of 100astronomical units, except at periapse, where the tidal interaction with the black hole stretches them along the orbit) and they show both thermal dust emission and line emission from ionized gas6-10. G1 and G2 have generated attention because they appear to be tidally interacting with the supermassive Galactic black hole, possibly enhancing its accretion activity. No broad consensus has yet been reached concerning their nature: the G objects show the characteristics of gas and dust clouds but display the dynamical properties of stellar-mass objects. Here we report observations of four additional G objects, all lying within 0.04parsecs of the black hole and forming a class that is probably unique to this environment. The widely varying orbits derived for the six G objects demonstrate that they were commonly but separately formed.

We present improved relative astrometry for stars within the central half parsec of our Galactic Center (GC) based on data obtained with the 10 m W. M. Keck Observatory from 1995 to 2017. The new methods used to improve the astrometric precision and accuracy include correcting for local astrometric distortions, applying a magnitude-dependent additive error, and more carefully removing instances of stellar confusion. Additionally, we adopt jackknife methods to calculate velocity and acceleration uncertainties. The resulting median proper motion uncertainty is 0.05 mas yr -1 for our complete sample of 1184 stars in the central 10″ (0.4 pc). We have detected 24 accelerating sources, 2.6 times more than the number of previously published accelerating sources, which extend out to 4″ (0.16 pc) from the black hole. Based on S0-2's orbit, our new astrometric analysis has reduced the systematic error of the supermassive black hole (SMBH) by a factor of 2. The linear drift in our astrometric reference frame is also reduced in the north-south direction by a factor of 4. We also find the first potential astrometric binary candidate S0-27 in the GC. These astrometric improvements provide a foundation for future studies of the origin and dynamics of the young stars around the SMBH, the structure and dynamics of the old nuclear star cluster, the SMBH's properties derived from orbits, and tests of general relativity in a strong gravitational field. © 2019. The American Astronomical Society. All rights reserved., We acknowledge support from the W. M. Keck Foundation, the Heising Simons Foundation, and the National Science Foundation (AST-1412615, AST1518273). The research leading to these results has received funding from the European Research Council under the European Union's Seventh Framework Programme (FP7/2007-2013)/ERC grant agreement No. [614922]. The W. M. Keck Observatory is operated as a scientific partnership among the California Institute of Technology, the University of California, and the National Aeronautics and Space Administration. The Observatory was made possible by the generous financial support of the W. M. Keck Foundation.

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., 17 pages, 10 figures

The general theory of relativity predicts that a star passing close to a supermassive black hole should exhibit a relativistic redshift. In this study, we used observations of the Galactic Center star S0-2 to test this prediction.We combined existing spectroscopic and astrometric measurements from 1995-2017, which cover S0-2's 16-year orbit, with measurements from March to September 2018, which cover three events during S0-2's closest approach to the black hole. We detected a combination of special relativistic and gravitational redshift, quantified using the redshift parameter ϒ. Our result, ϒ = 0.88 ± 0.17, is consistent with general relativity (ϒ = 1) and excludes a Newtonian model (U = 0) with a statistical significance of 5ω.Copyright © 2019 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science, This research was supported by the National Key Research and Development Program of China (2016YFD0100404), the National Natural Science Foundation of China (91535108), the Recruitment Program of Global Experts, and the Fundamental Research Funds for the Central Universities.

In various fields of physics and astronomy, access to experimental facilities or to telescopes is becoming more and more competitive and limited. It becomes therefore important to optimize the type of measurements and their scheduling to reach a given scientific objective and to increase the chances of success of a scientific project. In this communication, extending the work of Ford (2008) and of Loredo et al. (2012), we present an efficient adaptive scheduling tool aimed at prioritzing measurements in order to reach a scientific goal. The algorithm, based on the Fisher matrix, can be applied to a wide class of measurements. We present this algorithm in detail and discuss some practicalities such as systematic errors or measurements losses due to contigencies (such as weather, experimental failure, ...). As an illustration, we consider measurements of the short-period star S0-2 in our Galactic Center. We show that the radial velocity measurements at the two turning points of the radial velocity curve are more powerful for detecting the gravitational redshift than measurements at the maximal relativistic signal. We also explicitly present the methodology that was used to plan measurements in order to detect the relativistic redshift considering systematics and possible measurements losses. For the future, we identify the astrometric turning points to be highly sensitive to the relativistic advance of the periastron. Finally, we also identify measurements particularly sensitive to the distance to our Galactic Center: the radial velocities around periastron and the astrometric measurements just before closest approach and at the maximal right ascension astrometric turning point., 12 pages, 6 figures, submitted

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., Comment: Published in AJ. 23 pages, 14 figures. Monte Carlo chains are available in the published article, or are available upon request

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., Comment: 22 pages, 20 figures, accepted to ApJ Letters

OSIRIS is a near-infrared (1.0--2.4 $\mu$m) integral field spectrograph operating behind the adaptive optics system at Keck Observatory, and is one of the first lenslet-based integral field spectrographs. Since its commissioning in 2005, it has been a productive instrument, producing nearly half the laser guide star adaptive optics (LGS AO) papers on Keck. The complexity of its raw data format necessitated a custom data reduction pipeline (DRP) delivered with the instrument in order to iteratively assign flux in overlapping spectra to the proper spatial and spectral locations in a data cube. Other than bug fixes and updates required for hardware upgrades, the bulk of the DRP has not been updated since initial instrument commissioning. We report on the first major comprehensive characterization of the DRP using on-sky and calibration data. We also detail improvements to the DRP including characterization of the flux assignment algorithm; exploration of spatial rippling in the reduced data cubes; and improvements to several calibration files, including the rectification matrix, the bad pixel mask, and the wavelength solution. We present lessons learned from over a decade of OSIRIS data reduction that are relevant to the next generation of integral field spectrograph hardware and data reduction software design., Comment: 18 pages, 16 figures; accepted for publication in AJ

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

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

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, $\lambda_0 = 2.124 \text{ } \mu \text{m}$). We achieve a photometric uncertainty floor of $\Delta 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., Comment: 69 pages, 28 figures, 12 tables. Accepted for publication in The Astrophysical Journal

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., 9 pages, 3 figures, to appear in the Astrophysical Journal

Knowledge of the point spread function (PSF) is critical to many astronomical science cases. However, the PSF can be very difficult to estimate for cases where there are many crowded point sources or for observations of extended objects. Additionally, for adaptive optics observations, the PSF can be very complex with both spatial and temporal variability in the PSF. Integral-field spectroscopy behind adaptive optics is especially challenging because the fields of view are typically too small to sample the halo for even a single PSF. Here, we present a method for semi-empirical PSF reconstruction for integral field spectrographs using a combination of point source observations on a parallel imager, instrumental aberration measurements, and atmospheric turbulence profiles. This work builds upon the PSF reconstruction project AIROPA designed for imaging and extending it to IFU work (AIROPA-IFU). By using empirical calibrators from the parallel imager, which has a much larger field of view, and accounting for anisoplantic effects and instrumental aberrations, we can predict the PSF on the spectrograph. An important aspect is being able to predict the PSF at many different wavelengths based on observations from broad-band imaging. Here, we discuss how science cases such as observations of stars at the Galactic center can benefit from this method. We also establish metrics to quantitatively assess the performance of PSF reconstruction. We show that for bright stars, AIROPA-IFU can produce spectra with signal to noise ratio 50% higher than with simple aperture extraction of a data cube.

The integral field spectrograph OSIRIS at Keck I has been used to measure the motion of the stars around the supermassive black hole at the Center of the Galaxy. The small field of view provided and the crowding of the region prevent any good PSF estimate. A parallel imager can be used simultaneously to the IFU. However, its distance of 19 arcseconds prevents the observed PSF to be directly applied to the IFU because of anisoplanatism and instrumental aberrations. The Galactic Center Group at UCLA has developed an algorithms to predict PSF variability for Keck AO images (Off-axis PSF reconstruction, AIROPA software package). AIROPA allows us to use the parallel imager to correctly predict the IFU’s PSF. We modified this package to adapt it to the case of OSIRIS imager and IFU (AIROPA-IFU) and characterized the instrumental aberrations of both detectors. Here, we present preliminary results of the application of this post-processing tool to OSIRIS datasets of the Galactic Center.

A precise extinction law is a critical input when interpreting observations of highly reddened sources such as young star clusters and the Galactic Center (GC). We use Hubble Space Telescope observations of a region of moderate extinction and a region of high extinction to measure the optical and near-infrared extinction law (0.8 $\mu$m -- 2.2 $\mu$m). The moderate extinction region is the young massive cluster Westerlund 1 (Wd1; A$_{Ks} \sim$ 0.6 mag), where 453 proper motion-selected main-sequence stars are used to measure the shape of the extinction law. To quantify the shape we define the parameter $\mathcal{S}_{1/\lambda}$, which behaves similarly to a color excess ratio but is continuous as a function of wavelength. The high extinction region is the GC (A$_{Ks} \sim$ 2.5 mag), where 819 red clump stars are used to determine the normalization of the law. The best-fit extinction law is able to reproduce the Wd1 main sequence colors, which previous laws misestimate by 10%-30%. The law is inconsistent with a single power law, even when only the near-infrared filters are considered, and has A$_{F125W}$/A$_{Ks}$ and A$_{F814W}$/A$_{Ks}$ values that are 18% and 24% larger than the commonly used \citet{Nishiyama:2009fc} law, respectively. Using the law we recalculate the Wd1 distance to be 3896 $\pm$ 328 pc from published observations of eclipsing binary W13. This new extinction law should be used for highly reddened populations in the Milky Way, such as the Quintuplet cluster and Young Nuclear Cluster. A python code is provided to generate the law for future use., Comment: 32 pages, 20 figures; submitted to ApJ

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., Comment: Accepted to ApJ Lett, 5 figures, 8 pages

Dong, Hui et. al., Because of strong and spatially highly variable interstellar extinction and extreme source crowding, the faint (K ≥ 15) stellar population in the Milky Way's nuclear star cluster is still poorly studied. RR Lyrae stars provide us with a tool to estimate the mass of the oldest, relative dim stellar population. Recently, we analysed HST/WFC3/IR observations of the central 2.3 × 2.3 arcmin2 of the Milky Way and found 21 variable stars with periods between 0.2 and 1 d. Here, we present a further comprehensive analysis of these stars. The period- luminosity relationship of RR Lyrae is used to derive their extinctions and distances. Using multiple approaches, we classify our sample as 4 RRc stars, 4 RRab stars, 3 RRab candidates and 10 binaries. Especially, the four RRab stars show sawtooth light curves and fall exactly on to the Oosterhoff I division in the Bailey diagram. Compared to the RRab stars reported by Minniti et al., our new RRab stars have higher extinction (AK > 1.8) and should be closer to the Galactic Centre. The extinction and distance of one RRab stars match those for the Milky Way's nuclear star cluster given in previous works. We perform simulations and find that after correcting for incompleteness, there could be not more than 40 RRab stars within the Milky Way's nuclear star cluster and in our field of view. Through comparing with the known globular clusters of the Milky Way, we estimate that if there exists an old, metal-poor (-1.5 < [Fe/H] < -1) stellar population in the Milky Way nuclear star cluster on a scale of 5 × 5 pc, then it contributes at most 4.7 × 105 M⊙, i.e. ~18 per cent of the stellar mass. © 2017 The Authors., The research leading to these results has received funding from the European Research Council under the European Union's Seventh Framework Programme (FP7/2007-2013)/ERC grant agreement no [614922]. The work is also supported partly by NASA via the grant GO-14589, provided by the Space Telescope Science Institute. FN-L acknowledges financial support from an MECD predoctoral contract, code FPU14/01700. This work uses observations made with the NASA/ESA HST and the data archive at the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc. under NASA contract NAS 5-26555.

We used 4-yr baseline Hubble Space Telescope/Wide Field Camera 3 IR observations of the Galactic Centre in the F153M band (1.53 μm) to identify variable stars in the central ~2.3 arcmin × 2.3 arcmin field.We classified 3845 long-term (periods from months to years) and 76 short-term (periods of a few days or less) variables among a total sample of 33 070 stars. For 36 of the latter ones, we also derived their periods (< 3 d). Our catalogue not only confirms bright long period variables and massive eclipsing binaries identified in previous works but also contains many newly recognized dim variable stars. For example, we found δ Scuti and RR Lyrae stars towards the Galactic Centre for the first time, as well as one BL Her star (period < 1.3 d). We cross-correlated our catalogue with previous spectroscopic studies and found that 319 variables have well-defined stellar types, such as Wolf-Rayet, OB main sequence, supergiants and asymptotic giant branch stars. We used colours and magnitudes to infer the probable variable types for those stars without accurately measured periods or spectroscopic information. We conclude that the majority of unclassified variables could potentially be eclipsing/ellipsoidal binaries and Type II Cepheids. Our source catalogue will be valuable for future studies aimed at constraining the distance, star formation history and massive binary fraction of the Milky Way nuclear star cluster. © 2017 The Authors., The research leading to these results has received funding from the European Research Council under the European Union's Seventh Framework Programme (FP7/2007-2013) / ERC grant agreement no [614922]. The work is also supported by NASA via the grant AR-14589, provided by the Space Telescope Science Institute. FN-L acknowledges financial support from an MECD predoctoral contract, code FPU14/01700.

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|\alpha\right| < 0.016$ at a length scale of $\lambda=$ 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 \omega_\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., Comment: 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

Over two decades of astrometric and radial velocity data of short period stars in the Galactic center have the potential to provide unprecedented tests of General Relativity and insight into the astrophysics of supermassive black holes. Fundamental to this is understanding the underlying statistical issues of fitting stellar orbits. Unintended prior effects can obscure actual physical effects from General Relativity and the underlying extended mass distribution. At the heart of this is dealing with large parameter spaces inherent to multi star fitting and ensuring acceptable coverage properties of the resulting confidence intervals within the Bayesian framework. This proceeding will detail some of the UCLA Galactic Center Group's analysis and work in addressing these statistical issues., Comment: 2 pages, 1 figure. Proceedings IAUS 322

Optimal error estimation is key to achieve accurate photometry and astrometry. Stellar fluxes and positions in high angular resolution images are typically measured with PSF fitting routines, such as StarFinder. However, the formal uncertainties computed by these software packages tend to seriously underestimate the relevant uncertainties. We present a new approach to deal with this problem using a resampling method to obtain robust and reliable uncertainties without loss of sensitivity.

The Galactic center offers us a unique opportunity to test General Relativity (GR) with the orbits of stars around a supermassive black hole. Observations of these stars have been one of the great successes of adaptive optics on 8-10 m telescopes, driving the need for the highest angular resolution and astrometric precision. New tests of gravitational physics in the strong gravity regime with stellar orbits will be made possible through the leap in angular resolution and sensitivity from the next generation of extremely large ground-based telescopes. We present new simulations of specific science cases such as the detection of the GR precession of stars, the measurement of extended dark mass, and the distance to the Galactic center. We use realistic models of the adaptive optics system for TMT and the IRIS instrument to simulate these science cases. In additions, the simulations include observational issues such as the impact of source confusion on astrometry and radial velocities in the dense environment of the Galactic center. We qualitatively show how improvements in sensitivity, astrometric and spectroscopic precision, and increasing the number of stars affect the science with orbits at the Galactic center. We developed a tool to determine the constraints on physical models using a joint fit of over 100 stars that are expected to be observable with TMT. These science cases require very high astrometric precision and stability, thus they provide some of the most stringent constraints on the planned instruments and adaptive optics systems., Comment: 12 pages, 6 figures

The Galactic Center (GC) is dominated by the gravity of a super-massive black hole (SMBH), Sagittarius A$^*$, and is suspected to contain a sizable population of binary stars. Such binaries form hierarchical triples with the SMBH, undergoing Eccentric Kozai-Lidov (EKL) evolution, which can lead to high eccentricity excitations for the binary companions' mutual orbit. This effect can lead to stellar collisions or Roche-lobe crossings, as well as orbital shrinking due to tidal dissipation. In this work we investigate the dynamical and stellar evolution of such binary systems, especially with regards to the binaries' post-main-sequence evolution. We find that the majority of binaries (~75%) is eventually separated into single stars, while the remaining binaries (~25%) undergo phases of common-envelope evolution and/or stellar mergers. These objects can produce a number of different exotic outcomes, including rejuvenated stars, G2-like infrared-excess objects, stripped giant stars, Type Ia supernovae (SNe), cataclysmic variables (CVs), symbiotic binaries (SBs), or compact object binaries. We estimate that, within a sphere of 250 Mpc radius, about 7.5 to 15 Type Ia SNe per year should occur in galactic nuclei due to this mechanism, potentially detectable by ZTF and ASAS-SN. Likewise we estimate that, within a sphere of 1 Gpc$^3$ volume, about 10 to 20 compact object binaries form per year that could become gravitational wave sources. Based on results of EKL-driven compact object binary mergers in galactic nuclei by Hoang at al. (2018), this compact object binary formation rate translates to about 15 to 30 events per year detectable by Advanced LIGO., 8 pages, 3 figures, accepted by ApJ

The Quintuplet star cluster is one of only three known young ($10^4$ M$_\odot$) clusters within $\sim100$ pc of the Galactic Center. In order to explore star cluster formation and evolution in this extreme environment, we analyze the Quintuplet's dynamical structure. Using the HST WFC3-IR instrument, we take astrometric and photometric observations of the Quintuplet covering a $120''\times120''$ field-of-view, which is $19$ times larger than those of previous proper motion studies of the Quintuplet. We generate a catalog of the Quintuplet region with multi-band, near-infrared photometry, proper motions, and cluster membership probabilities for $10,543$ stars. We present the radial density profile of $715$ candidate Quintuplet cluster members with $M\gtrsim4.7$ M$_\odot$ out to $3.2$ pc from the cluster center. A $3��$ lower limit of $3$ pc is placed on the tidal radius, indicating the lack of a tidal truncation within this radius range. Only weak evidence for mass segregation is found, in contrast to the strong mass segregation found in the Arches cluster, a second and slightly younger massive cluster near the Galactic Center. It is possible that tidal stripping hampers a mass segregation signature, though we find no evidence of spatial asymmetry. Assuming that the Arches and Quintuplet formed with comparable extent, our measurement of the Quintuplet's comparatively large core radius of $0.62^{+0.10}_{-0.10}$ pc provides strong empirical evidence that young massive clusters in the Galactic Center dissolve on a several Myr timescale., 25 pages (21-page main text, 4-page appendix), 18 figures, submitted to ApJ

As a young massive cluster in the Central Molecular Zone, the Arches cluster is a valuable probe of the stellar Initial Mass Function (IMF) in the extreme Galactic Center environment. We use multi-epoch Hubble Space Telescope observations to obtain high-precision proper motion and photometric measurements of the cluster, calculating cluster membership probabilities for stars down to 1.8 M$_{\odot}$ between cluster radii of 0.25 pc -- 3.0 pc. We achieve a cluster sample with just ~8% field contamination, a significant improvement over photometrically-selected samples which are severely compromised by the differential extinction across the field. Combining this sample with K-band spectroscopy of 5 cluster members, we forward model the Arches cluster to simultaneously constrain its IMF and other properties (such as age and total mass) while accounting for observational uncertainties, completeness, mass segregation, and stellar multiplicity. We find that the Arches IMF is best described by a 1-segment power law that is significantly top-heavy: $\alpha$ = 1.80 $\pm$ 0.05 (stat) $\pm$ 0.06 (sys), where dN/dm $\propto$ m$^{-\alpha}$, though we cannot discount a 2-segment power law model with a high-mass slope only slightly shallower than local star forming regions ($\alpha$ = 2.04$^{+0.14}_{-0.19}$ $\pm$ 0.04) but with a break at 5.8$^{+3.2}_{-1.2}$ $\pm$ 0.02 M$_{\odot}$. In either case, the Arches IMF is significantly different than the standard IMF. Comparing the Arches to other young massive clusters in the Milky Way, we find tentative evidence for a systematically top-heavy IMF at the Galactic Center., Comment: 31 pages, 25 figures. Accepted to ApJ; Revised extinction law, (very) minor changes to results, no impact on conclusions

The central parsec of our Galaxy hosts not only a supermassive black hole, but also a large population of young stars (age

We discuss observations of Sagittarius A* with NACO@VLT in K-band and recent synchronous observations with NIRC2@Keck II and OSIRIS@Keck I in L′-band and H-band, respectively. The variability of Sagittarius A* in the near infrared is a continuous one-state process that can be described by a pure red-noise process having a timescale of a few hours. We describe this process and its properties in detail. Our newest observations with the Keck telescopes represent the first truly synchronous high cadence data set to test for time variability of the spectral index within the near infrared. We discovered a time-variable spectral index that might be interpreted as a time lag of the L′-band with respect to the H-band.