Source author record

Jessica R. Lu

Jessica R. Lu appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

32works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

32 published item(s)

preprint2022arXiv

AIROPA III: Testing Simulated and On-Sky Data

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.

preprint2022arXiv

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

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.

preprint2022arXiv

An isolated mass gap black hole or neutron star detected with astrometric microlensing

We present the analysis of five black hole candidates identified from gravitational microlensing surveys. Hubble Space Telescope astrometric data and densely sampled lightcurves from ground-based microlensing surveys are fit with a single-source, single-lens microlensing model in order to measure the mass and luminosity of each lens and determine if it is a black hole. One of the five targets (OGLE-2011-BLG-0462/MOA-2011-BLG-191 or OB110462 for short) shows a significant $>1$ mas coherent astrometric shift, little to no lens flux, and has an inferred lens mass of 1.6 - 4.4 $M_\odot$. This makes OB110462 the first definitive discovery of a compact object through astrometric microlensing and it is most likely either a neutron star or a low-mass black hole. This compact object lens is relatively nearby (0.70-1.92 kpc) and has a slow transverse motion of $<$30 km/s. OB110462 shows significant tension between models well-fit to photometry vs. astrometry, making it currently difficult to distinguish between a neutron star and a black hole. Additional observations and modeling with more complex system geometries, such as binary sources are needed to resolve the puzzling nature of this object. For the remaining four candidates, the lens masses are $<2 M_\odot$ and they are unlikely to be black holes; two of the four are likely white dwarfs or neutron stars. We compare the full sample of five candidates to theoretical expectations on the number of black holes in the Milky Way ($\sim 10^8$) and find reasonable agreement given the small sample size.

preprint2022arXiv

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

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

preprint2021arXiv

MOA-2009-BLG-319Lb: A Sub-Saturn Planet Inside the Predicted Mass Desert

We present an adaptive optics (AO) analysis of images from the Keck-II telescope NIRC2 instrument of the planetary microlensing event MOA-2009-BLG-319. The $\sim$10 year baseline between the event and the Keck observations allows the planetary host star to be detected at a separation of $66.5\pm 1.7\,$mas from the source star, consistent with the light curve model prediction. The combination of the host star brightness and light curve parameters yield host star and planet masses of M_host = 0.514 $\pm$ 0.063M_Sun and m_p = 66.0 $\pm$ 8.1M_Earth at a distance of $D_L = 7.0 \pm 0.7\,$kpc. The star-planet projected separation is $2.03 \pm 0.21\,$AU. The planet-star mass ratio of this system, $q = (3.857 \pm 0.029)\times 10^{-4}$, places it in the predicted "planet desert" at $10^{-4} < q < 4\times 10^{-4}$ according to the runaway gas accretion scenario of the core accretion theory. Seven of the 30 planets in the Suzuki et al. (2016) sample fall in this mass ratio range, and this is the third with a measured host mass. All three of these host stars have masses of 0.5 $\leq$ M_host/M_Sun $\leq$ 0.7, which implies that this predicted mass ratio gap is filled with planets that have host stars within a factor of two of 1M_Sun. This suggests that runaway gas accretion does not play a major role in determining giant planet masses for stars somewhat less massive than the Sun. Our analysis has been accomplished with a modified DAOPHOT code that has been designed to measure the brightness and positions of closely blended stars. This will aid in the development of the primary method that the Nancy Grace Roman Space Telescope mission will use to determine the masses of microlens planets and their hosts.

preprint2021arXiv

The 3-D Kinematics of the Orion Nebula Cluster: NIRSPEC-AO Radial Velocities of the Core Population

The kinematics and dynamics of stellar and substellar populations within young, still-forming clusters provides valuable information for constraining theories of formation mechanisms. Using Keck II NIRSPEC+AO data, we have measured radial velocities for 56 low-mass sources within 4' of the core of the Orion Nebula Cluster (ONC). We also re-measure radial velocities for 172 sources observed with SDSS/APOGEE. These data are combined with proper motions measured using $HST$ ACS/WFPC2/WFC3IR and Keck II NIRC2, creating a sample of 135 sources with all three velocity components. The velocities measured are consistent with a normal distribution in all three components. We measure intrinsic velocity dispersions of ($σ_{v_α}$, $σ_{v_δ}$, $σ_{v_r}$) = ($1.64\pm0.12$, $2.03\pm0.13$, $2.56^{+0.16}_{-0.17}$) km s$^{-1}$. Our computed intrinsic velocity dispersion profiles are consistent with the dynamical equilibrium models from Da Rio et al. (2014) in the tangential direction, but not in the line of sight direction, possibly indicating that the core of the ONC is not yet virialized, and may require a non-spherical potential to explain the observed velocity dispersion profiles. We also observe a slight elongation along the north-south direction following the filament, which has been well studied in previous literature, and an elongation in the line of sight to tangential velocity direction. These 3-D kinematics will help in the development of realistic models of the formation and early evolution of massive clusters.

preprint2020arXiv

Gravitational Microlensing Event Statistics for the Zwicky Transient Facility

Microlensing surveys have discovered thousands of events with almost all events discovered within the Galactic bulge or toward the Magellanic clouds. The Zwicky Transient Facility (ZTF), while not designed to be a microlensing campaign, is an optical time-domain survey that observes the entire northern sky every few nights including the Galactic plane. ZTF observes $\sim10^9$ stars in g-band and r-band and can significantly contribute to the observed microlensing population. We predict that ZTF will observe $\sim$1100 microlensing events in three years of observing within $10^\circ$ degrees latitude of the Galactic plane, with $\sim$500 events in the outer Galaxy ($\ell \geq 10^\circ$). This yield increases to $\sim$1400 ($\sim$800) events by combining every three ZTF exposures, $\sim$1800 ($\sim$900) events if ZTF observes for a total of five years, and $\sim$2400 ($\sim$1300) events for a five year survey with post-processing image stacking. Using the microlensing modeling software PopSyCLE, we compare the microlensing populations in the Galactic bulge and the outer Galaxy. We also present an analysis of the microlensing event ZTF18abhxjmj to demonstrate how to leverage these population statistics in event modeling. ZTF will constrain Galactic structure, stellar populations, and primordial black holes through photometric microlensing.

preprint2020arXiv

Keck Observations Confirm a Super-Jupiter Planet Orbiting M-dwarf OGLE-2005-BLG-071L

We present adaptive optics imaging from the NIRC2 instrument on the Keck-2 telescope that resolves the exoplanet host (and lens) star as it separates from the brighter source star. These observations yield the $K$-band brightness of the lens and planetary host star, as well as the lens-source relative proper motion, $μ_{\rm rel,H}$. in the heliocentric reference frame. The $μ_{\rm rel,H}$ measurement allows determination of the microlensing parallax vector, $π_E$, which had only a single component determined by the microlensing light curve. The combined measurements of $μ_{\rm rel,H}$ and $K_L$ provide the masses of the host stat, $M_{\rm host} = 0.426\pm 0.037 M_\odot$, and planet, $m_p = 3.27 \pm 0.32 M_{\rm Jup}$ with a projected separation of $3.4\pm 0.5\,$AU. This confirms the tentative conclusion of a previous paper (Dong et al. 2009) that this super-Jupiter mass planet, OGLE-2005-BLG-071Lb, orbits an M-dwarf. Such planets are predicted to be rare by the core accretion theory and have been difficult to find with other methods, but there are two such planets with firm mass measurements from microlensing, and an additional 11 planetary microlens events with host mass estimates $< 0.5M_\odot$ and planet mass estimates $> 2$ Jupiter masses that could be confirmed by high angular follow-up observations. We also point out that OGLE-2005-BLG-071L has separated far enough from its host star that it should be possible to measure the host star metallicity withspectra from a high angular resolution telescope such as Keck, the VLT, the Hubble Space Telescope or the James Webb Space Telescope.

preprint2020arXiv

MOA-2007-BLG-400 A Super-Jupiter Mass Planet Orbiting a Galactic BulgeK-dwarf Revealed by Keck Adaptive Optics Imaging

We present Keck/NIRC2 adaptive optics imaging of planetary microlensing event MOA-2007-BLG-400 that resolves the lens star system from the source. We find that the MOA-2007-BLG-400L planetary system consists of a $1.71\pm 0.27 M_{\rm Jup}$ planet orbiting a $0.69\pm 0.04M_{\odot}$ K-dwarf host star at a distance of $6.89\pm 0.77\,$kpc from the Sun. So, this planetary system probably resides in the Galactic bulge. The planet-host star projected separation is only weakly constrained due to the close-wide light curve degeneracy; the 2$σ$ projected separation range is 0.6--$7.2\,$AU. This host mass is at the top end of the range of masses predicted by a standard Bayesian analysis that assumes that all stars have an equal chance of hosting a star of the observed mass ratio. This and the similar result for event MOA-2013-BLG-220 suggests that more massive stars may be more likely to host planets with a mass ratio in the $0.002 < q < 0.004$ range that orbit beyond the snow line. These results also indicate the importance of host star mass measurements for exoplanets found by microlensing. The microlensing survey imaging data from NASA's Nancy Grace Roman Space Telescope (formerly WFIRST) mission will be doing mass measurements like this for a huge number of planetary events. This host lens is the highest contrast lens-source detected in microlensing mass measurement analysis (the lens being 10$\times$ fainter than the source). We present an improved method of calculating photometry and astrometry uncertainties based on the Jackknife method, which produces more accurate errors that are $\sim$$2.5 \times$ larger than previous estimates.

preprint2019arXiv

PopSyCLE: A New Population Synthesis Code for Compact Object Microlensing Events

We present a new Milky Way microlensing simulation code, dubbed PopSyCLE (Population Synthesis for Compact object Lensing Events). PopSyCLE is the first resolved microlensing simulation to include a compact object distribution derived from numerical supernovae explosion models and both astrometric and photometric microlensing effects. We demonstrate the capabilities of PopSyCLE by investigating the optimal way to find black holes (BHs) with microlensing. Candidate BHs have typically been selected from wide-field photometric microlensing surveys, such as OGLE, by selecting events with long Einstein crossing times ($t_E>120$ days). These events can be selected at closest approach and monitored astrometrically in order to constrain the mass of each lens; PopSyCLE predicts a BH detection rate of ~40% for such a program. We find that the detection rate can be enhanced to ~85% by selecting events with both $t_E>120$ days and a microlensing parallax of $π_E<0.08$. Unfortunately, such a selection criterion cannot be applied during the event as $π_E$ requires both pre- and post-peak photometry. However, historical microlensing events from photometric surveys can be revisited using this new selection criteria in order to statistically constrain the abundance of BHs in the Milky Way. The future WFIRST microlensing survey provides both precise photometry and astrometry and will yield individual masses of $\mathcal{O}(100-1000) black holes, which is at least an order of magnitude more than is possible with individual candidate follow-up with current facilities. The resulting sample of BH masses from WFIRST will begin to constrain the shape of the black hole present-day mass function, BH multiplicity, and BH kick velocity distributions.

preprint2016arXiv

Radial Trends in IMF-Sensitive Absorption Features in Two Early-Type Galaxies: Evidence for Abundance-Driven Gradients

Samples of early-type galaxies show a correlation between stellar velocity dispersion and the stellar initial mass function (IMF) as inferred from gravity-sensitive absorption lines in the galaxies' central regions. To search for spatial variations in the IMF, we have observed two early-type galaxies with Keck/LRIS and measured radial gradients in the strengths of absorption features from 4000-5500 Å$\,$ and 8000-10,000 Å. We present spatially resolved measurements of the dwarf-sensitive spectral indices NaI (8190 Å) and Wing-Ford FeH (9915 Å), as well as indices for species of H, C$_2$, CN, Mg, Ca, TiO, and Fe. Our measurements show a metallicity gradient in both objects, and Mg/Fe consistent with a shallow gradient in α-enhancement, matching widely observed trends for massive early-type galaxies. The NaI index and the CN$_1$ index at 4160 Å$\,$ exhibit significantly steeper gradients, with a break at $r \sim 0.1 r_{\rm eff}$ ($r \sim 300$ pc). Inside this radius NaI strength increases sharply toward the galaxy center, consistent with a rapid central rise in [Na/Fe]. In contrast, the ratio of FeH to Fe index strength decreases toward the galaxy center. This behavior cannot be reproduced by a steepening IMF inside $0.1 r_{\rm eff}$ if the IMF is a single power law. While gradients in the mass function above $\sim 0.4 M_\odot$ may occur, exceptional care is required to disentangle these IMF variations from the extreme variations in individual element abundances near the galaxies' centers.

preprint2016arXiv

The InfraRed Imaging Spectrograph (IRIS) for TMT: latest science cases and simulations

The Thirty Meter Telescope (TMT) first light instrument IRIS (Infrared Imaging Spectrograph) will complete its preliminary design phase in 2016. The IRIS instrument design includes a near-infrared (0.85 - 2.4 micron) integral field spectrograph (IFS) and imager that are able to conduct simultaneous diffraction-limited observations behind the advanced adaptive optics system NFIRAOS. The IRIS science cases have continued to be developed and new science studies have been investigated to aid in technical performance and design requirements. In this development phase, the IRIS science team has paid particular attention to the selection of filters, gratings, sensitivities of the entire system, and science cases that will benefit from the parallel mode of the IFS and imaging camera. We present new science cases for IRIS using the latest end-to-end data simulator on the following topics: Solar System bodies, the Galactic center, active galactic nuclei (AGN), and distant gravitationally-lensed galaxies. We then briefly discuss the necessity of an advanced data management system and data reduction pipeline.

preprint2016arXiv

The Rapid Transient Surveyor

The Rapid Transient Surveyor (RTS) is a proposed rapid-response, high-cadence adaptive optics (AO) facility for the UH 2.2-m telescope on Maunakea. RTS will uniquely address the need for high-acuity and sensitive near-infrared spectral follow-up observations of tens of thousands of objects in mere months by combining an excellent observing site, unmatched robotic observational efficiency, and an AO system that significantly increases both sensitivity and spatial resolving power. We will initially use RTS to obtain the infrared spectra of ~4,000 Type Ia supernovae identified by the Asteroid Terrestrial-Impact Last Alert System over a two year period that will be crucial to precisely measuring distances and mapping the distribution of dark matter in the z < 0.1 universe. RTS will comprise an upgraded version of the Robo-AO laser AO system and will respond quickly to target-of-opportunity events, minimizing the time between discovery and characterization. RTS will acquire simultaneous-multicolor images with an acuity of 0.07-0.10" across the entire visible spectrum (20% i'-band Strehl in median conditions) and <0.16" in the near infrared, and will detect companions at 0.5" at contrast ratio of ~500. The system will include a high-efficiency prism integral field unit spectrograph: R = 70-140 over a total bandpass of 840-1830 nm with an 8.7" by 6.0" field of view (0.15" spaxels). The AO correction boosts the infrared point-source sensitivity of the spectrograph against the sky background by a factor of seven for faint targets, giving the UH 2.2-m the H-band sensitivity of a 5.7-m telescope without AO.

preprint2015arXiv

Circumstellar discs in Galactic centre clusters: Disc-bearing B-type stars in the Quintuplet and Arches clusters

We investigate the circumstellar disc fraction as determined from L-band excess observations of the young, massive Arches and Quintuplet clusters residing in the central molecular zone of the Milky Way. The Quintuplet cluster was searched for L-band excess sources for the first time. We find a total of 26 excess sources in the Quintuplet cluster and 21 in the Arches cluster, of which 13 are new detections. With the aid of proper motion membership samples, the disc fraction of the Quintuplet cluster was derived for the first time to be 4.0 +/- 0.7%. There is no evidence for a radially varying disc fraction in this cluster. In the case of the Arches cluster, a disc fraction of 9.2 +/- 1.2% approximately out to the cluster's predicted tidal radius, r < 1.5 pc, is observed. This excess fraction is consistent with our previously found disc fraction in the cluster in the radial range 0.3 < r < 0.8 pc. In both clusters, the host star mass range covers late A- to early B-type stars, 2 < M < 15 Msun, as derived from J-band photospheric magnitudes. We discuss the unexpected finding of dusty circumstellar discs in these UV intense environments in the context of primordial disc survival and formation scenarios of secondary discs. We consider the possibility that the L-band excess sources in the Arches and Quintuplet clusters could be the high-mass counterparts to T Tauri pre-transitional discs. As such a scenario requires a long pre-transitional disc lifetime in a UV intense environment, we suggest that mass transfer discs in binary systems are a likely formation mechanism for the B-star discs observed in these starburst clusters.

preprint2015arXiv

HST/WFC3 Observations of an Off-Nuclear Superbubble in Arp 220

We present a high spatial resolution optical and infrared study of the circumnuclear region in Arp 220, a late-stage galaxy merger. Narrowband imaging using HST/WFC3 has resolved the previously observed peak in H$α$+[NII] emission into a bubble-shaped feature. This feature measures 1.6" in diameter, or 600 pc, and is only 1" northwest of the western nucleus. The bubble is aligned with the western nucleus and the large-scale outflow axis seen in X-rays. We explore several possibilities for the bubble origin, including a jet or outflow from a hidden active galactic nucleus (AGN), outflows from high levels of star formation within the few hundred pc nuclear gas disk, or an ultraluminous X-ray source. An obscured AGN or high levels of star formation within the inner $\sim$100 pc of the nuclei are favored based on the alignment of the bubble and energetics arguments.

preprint2015arXiv

Mapping the Outer Edge of the Young Stellar Cluster in the Galactic Center

We present new near-infrared spectroscopic observations of the outer edges of the young stellar cluster around the supermassive black hole at the Galactic center. The observations show a break in the surface-density profile of young stars at approximately 13 arcsec (0.52 pc). These observations spectroscopically confirm previous suggestions of a break based on photometry. Using Gemini North's Near-Infrared Integral Field Spectrometer (NIFS) we are able to detect and separate early- and late-type stars with a 75% completeness at Ks = 15.5. We sample a region with radii between 7" to 23" (0.28 pc to 0.92 pc) from Sgr A*, and present new spectral classifications of 144 stars brighter than Ks = 15.5, where 140 stars are late-type (> 1 Gyr) and only four stars are early-type (young, 4-6 Myr). A broken power-law fit of the early-type surface-density matches well with our data and previously published values. The projected surface-density of late-type stars is also measured and found to be consistent with previous results. We find that the observed early-type surface-density profile is inconsistent with the theory of the young stars originating from a tightly bound infalling cluster, as no significant trail of young stars is found at radii above 13". We also note that either a simple disk instability criterion or a cloud-cloud collision could explain the location of the outer edge, though we lack information to make conclusive remarks on either alternative. If this break in surface-density represents an edge to the young stellar cluster it would set an important scale for the most recent episode of star formation at the Galactic center.

preprint2015arXiv

The Arches Cluster: Extended Structure and Tidal Radius

At a projected distance of ~26 pc from Sgr A*, the Arches cluster provides insight to star formation in the extreme Galactic Center (GC) environment. Despite its importance, many key properties such as the cluster's internal structure and orbital history are not well known. We present an astrometric and photometric study of the outer region of the Arches cluster (R > 6.25") using HST WFC3IR. Using proper motions we calculate membership probabilities for stars down to F153M = 20 mag (~2.5 M_sun) over a 120" x 120" field of view, an area 144 times larger than previous astrometric studies of the cluster. We construct the radial profile of the Arches to a radius of 75" (~3 pc at 8 kpc), which can be well described by a single power law. From this profile we place a 3-sigma lower limit of 2.8 pc on the observed tidal radius, which is larger than the predicted tidal radius (1 - 2.5 pc). Evidence of mass segregation is observed throughout the cluster and no tidal tail structures are apparent along the orbital path. The absence of breaks in the profile suggests that the Arches has not likely experienced its closest approach to the GC between ~0.2 - 1 Myr ago. If accurate, this constraint indicates that the cluster is on a prograde orbit and is located front of the sky plane that intersects Sgr A*. However, further simulations of clusters in the GC potential are required to interpret the observed profile with more confidence.

preprint2014arXiv

Astrometric performance of the Gemini multi-conjugate adaptive optics system in crowded fields

The Gemini Multi-conjugate adaptive optics System (GeMS) is a facility instrument for the Gemini-South telescope. It delivers uniform, near-diffraction-limited image quality at near-infrared wavelengths over a 2 arcminute field of view. Together with the Gemini South Adaptive Optics Imager (GSAOI), a near-infrared wide field camera, GeMS/GSAOI's combination of high spatial resolution and a large field of view will make it a premier facility for precision astrometry. Potential astrometric science cases cover a broad range of topics including exo-planets, star formation, stellar evolution, star clusters, nearby galaxies, black holes and neutron stars, and the Galactic center. In this paper, we assess the astrometric performance and limitations of GeMS/GSAOI. In particular, we analyze deep, mono-epoch images, multi-epoch data and distortion calibration. We find that for single-epoch, un-dithered data, an astrometric error below 0.2 mas can be achieved for exposure times exceeding one minute, provided enough stars are available to remove high-order distortions. We show however that such performance is not reproducible for multi-epoch observations, and an additional systematic error of ~0.4 mas is evidenced. This systematic multi-epoch error is the dominant error term in the GeMS/GSAOI astrometric error budget, and it is thought to be due to time-variable distortion induced by gravity flexure.

preprint2014arXiv

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

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).

preprint2014arXiv

Properties of the Remnant Clockwise Disk of Young Stars in the Galactic Center

We present new kinematic measurements and modeling of a sample of 116 young stars in the central parsec of the Galaxy in order to investigate the properties of the young stellar disk. The measurements were derived from a combination of speckle and laser guide star adaptive optics imaging and integral field spectroscopy from the Keck telescopes. Compared to earlier disk studies, the most important kinematic measurement improvement is in the precision of the accelerations in the plane of the sky, which have a factor of six smaller uncertainties (~10 uas/yr/yr). We have also added the first radial velocity measurements for 8 young stars, increasing the sample at the largest radii (6"-12") by 25%. We derive the ensemble properties of the observed stars using Monte-Carlo simulations of mock data. There is one highly significant kinematic feature (~20 sigma), corresponding to the well-known clockwise disk, and no significant feature is detected at the location of the previously claimed counterclockwise disk. The true disk fraction is estimated to be ~20%, a factor of ~2.5 lower than previous claims, suggesting that we may be observing the remnant of what used to be a more densely populated stellar disk. The similarity in the kinematic properties of the B stars and the O/WR stars suggests a common star formation event. The intrinsic eccentricity distribution of the disk stars is unimodal, with an average value of <e> = 0.27 +/- 0.07, which we show can be achieved through dynamical relaxation in an initially circular disk with a moderately top-heavy mass function.

preprint2014arXiv

Second generation Robo-AO instruments and systems

The prototype Robo-AO system at the Palomar Observatory 1.5-m telescope is the world's first fully automated laser adaptive optics instrument. Scientific operations commenced in June 2012 and more than 12,000 observations have since been performed at the ~0.12" visible-light diffraction limit. Two new infrared cameras providing high-speed tip-tilt sensing and a 2' field-of-view will be integrated in 2014. In addition to a Robo-AO clone for the 2-m IGO and the natural guide star variant KAPAO at the 1-m Table Mountain telescope, a second generation of facility-class Robo-AO systems are in development for the 2.2-m University of Hawai'i and 3-m IRTF telescopes which will provide higher Strehl ratios, sharper imaging, ~0.07", and correction to λ = 400 nm.

preprint2014arXiv

The InfraRed Imaging Spectrograph (IRIS) for TMT: Overview of innovative science programs

IRIS (InfraRed Imaging Spectrograph) is a first light near-infrared diffraction limited imager and integral field spectrograph being designed for the future Thirty Meter Telescope (TMT). IRIS is optimized to perform astronomical studies across a significant fraction of cosmic time, from our Solar System to distant newly formed galaxies (Barton et al. [1]). We present a selection of the innovative science cases that are unique to IRIS in the era of upcoming space and ground-based telescopes. We focus on integral field spectroscopy of directly imaged exoplanet atmospheres, probing fundamental physics in the Galactic Center, measuring 10^4 to 10^10 Msun supermassive black hole masses, resolved spectroscopy of young star-forming galaxies (1 < z < 5) and first light galaxies (6 < z < 12), and resolved spectroscopy of strong gravitational lensed sources to measure dark matter substructure. For each of these science cases we use the IRIS simulator (Wright et al. [2], Do et al. [3]) to explore IRIS capabilities. To highlight the unique IRIS capabilities, we also update the point and resolved source sensitivities for the integral field spectrograph (IFS) in all five broadband filters (Z, Y, J, H, K) for the finest spatial scale of 0.004" per spaxel. We briefly discuss future development plans for the data reduction pipeline and quicklook software for the IRIS instrument suite.

preprint2014arXiv

The orbital motion of the Quintuplet cluster - a common origin for the Arches and Quintuplet clusters?

We investigate the orbital motion of the Quintuplet cluster near the Galactic center with the aim of constraining formation scenarios of young, massive star clusters in nuclear environments. Three epochs of adaptive optics high-angular resolution imaging with Keck/NIRC2 and VLT/NACO were obtained over a time baseline of 5.8 years, delivering an astrometric accuracy of 0.5-1 mas/yr. Proper motions were derived in the cluster reference frame and were used to distinguish cluster members from the majority of field stars. Fitting the cluster and field proper motion distributions with 2D gaussian models, we derive the orbital motion of the cluster for the first time. The Quintuplet is moving with a 2D velocity of 132 +/- 15 km/s with respect to the field along the Galactic plane, which yields a 3D orbital velocity of 167 +/- 15 km/s when combined with the previously known radial velocity. From a sample of 119 stars measured in three epochs, we derive an upper limit to the velocity dispersion in the core of the Quintuplet cluster of sigma_1D < 10 km/s. Knowledge of the three velocity components of the Quintuplet allows us to model the cluster orbit in the potential of the inner Galaxy. Comparing the Quintuplet's orbit with the Arches orbit, we discuss the possibility that both clusters originated in the same area of the central molecular zone. [abridged]

preprint2013arXiv

Keck Adaptive Optics Observations of the Protostellar Disk around Radio Source I in the Orion Kleinmann-Low Nebula

We have made the first detection of a near-infrared counterpart associated with the disk around Radio Source "I," a massive protostar in the Kleinmann-Low Nebula in Orion using imaging with laser guide star adaptive optics on the Keck II telescope. The infrared emission is evident in images acquired using L' (3.8 microns) and Ms (4.7 microns) filters and is not detectable at K' (2.1 microns). The observed morphology strongly suggests that we are seeing some combination of scattered and thermal light emanating from the disk. The disk is also manifest in the L'/Ms flux ratio image. We interpret the near-infrared emission as the illuminated surface of a nearly edge-on disk, oriented so that only the northern face is visible; the opposite surface remains hidden by the disk. We do not see infrared radiation associated directly with the star proposed to be associated with Source "I." The data also suggest that there is a cavity above and below the disk that is oriented perpendicular to the disk, and is sculpted by the known, strong outflow from the inner disk of Source I. We compare our data to models of a protostar with a surrounding disk, envelope, and wind-blown cavity in order to elucidate the nature of the disk around Radio Source I.

preprint2013arXiv

Stellar Populations in the Central 0.5 pc of the Galaxy I: A New Method for Constructing Luminosity Functions and Surface-Density Profiles

We present new high angular resolution near-infrared spectroscopic observations of the nuclear star cluster surrounding the Milky Way's central supermassive black hole. Using the integral-field spectrograph OSIRIS on Keck II behind the laser-guide-star adaptive optics system, this spectroscopic survey enables us to separate early-type (young, 4-6 Myr) and late-type (old, >1 Gyr) stars with a completeness of 50% down to K' = 15.5 mag, which corresponds to ~10 \msun for the early-type stars. This work increases the radial extent of reported OSIRIS/Keck measurements by more than a factor of 3 from 4" to 14" (0.16 pc to 0.56 pc), along the projected disk of young stars. For our analysis, we implement a new method of completeness correction using a combination of star-planting simulations and Bayesian inference. We assign probabilities for the spectral type of every source detected in deep imaging down to K' = 15.5 mag using information from spectra, simulations, number counts, and the distribution of stars. The inferred radial surface-density profiles, $Σ(R) \propto R^{-Γ}$, for the young stars and late-type giants are consistent with earlier results ($Γ_{early} = 0.93 \pm 0.09$, $Γ_{late} = 0.16 \pm 0.07$). The late-type surface-density profile is approximately flat out to the edge of the survey. While the late-type stellar luminosity function is consistent with the Galactic bulge, the completeness-corrected luminosity function of the early-type stars has significantly more young stars at faint magnitudes compared to previous surveys with similar depth. This luminosity function indicates that the corresponding mass function of the young stars is likely less top-heavy than that inferred from previous surveys.

preprint2013arXiv

Stellar Populations in the Central 0.5 pc of the Galaxy II: The Initial Mass Function

The supermassive black hole at the center of the Milky Way plays host to a massive, young cluster that may have formed in one of the most inhospitable environments in the Galaxy. We present new measurements of the global properties of this cluster, including the initial mass function (IMF), age, and cluster mass. These results are based on Keck laser-guide-star adaptive optics observations used to identify the young stars and measure their Kp-band luminosity function as presented in Do et al. 2013. A Bayesian inference methodology is developed to simultaneously fit the global properties of the cluster utilizing the observations and extensive simulations of synthetic star clusters. We find that the slope of the mass function for this cluster is α= 1.7 +/- 0.2, which is steeper than previously reported, but still flatter than the traditional Salpeter slope of 2.35. The age of the cluster is between 2.5-5.8 Myr with 95% confidence, which is a younger age than typically adopted but consistent within the uncertainties of past measurements. The exact age of the cluster is difficult to determine since our results show two distinct age solutions (3.9 Myr and 2.8 Myr) due to model degeneracies in the relative number of Wolf-Rayet and OB stars. The total cluster mass is between 14,000 - 37,000 \msun above 1 \msun and it is necessary to include multiple star systems in order to fit the observed luminosity function and the number of observed Wolf-Rayet stars. The new IMF slope measurement is now consistent with X-ray observations indicating a factor of 10 fewer X-ray emitting pre-main-sequence stars than expected when compared with a Salpeter IMF. The young cluster at the Galactic center is one of the few definitive examples of an IMF that deviates significantly from the near-universal IMFs found in the solar neighborhood.

preprint2010arXiv

Improving Galactic Center Astrometry by Reducing the Effects of Geometric Distortion

We present significantly improved proper motion measurements of the Milky Way's central stellar cluster. These improvements are made possible by refining our astrometric reference frame with a new geometric optical distortion model for the W. M. Keck II 10 m telescope's Adaptive Optics camera (NIRC2) in its narrow field mode. For the first time, this distortion model is constructed from on-sky measurements, and is made available to the public. When applied to widely dithered images, it produces residuals in the separations of stars that are a factor of ~3 smaller compared to the outcome using previous models. By applying this new model, along with corrections for differential atmospheric refraction, to widely dithered images of SiO masers at the Galactic center, we improve our ability to tie into the precisely measured radio Sgr A*-rest frame. The resulting infrared reference frame is ~2-3 times more accurate and stable than earlier published efforts. In this reference frame, Sgr A* is localized to within a position of 0.6 mas and a velocity of 0.09 mas/yr, or ~3.4 km/s at 8 kpc (1 sigma). While earlier proper motion studies defined a reference frame by assuming no net motion of the stellar cluster, this approach is fundamentally limited by the cluster's intrinsic dispersion and therefore will not improve with time. We define a reference frame with SiO masers and this reference frame's stability should improve steadily with future measurements of the SiO masers in this region. This is essential for achieving the necessary reference frame stability required to detect the effects of general relativity and extended mass on short-period stars at the Galactic center.

preprint2010arXiv

Increasing the Scientific Return of Stellar Orbits at the Galactic Center

We report a factor of $\sim$3 improvement in Keck Laser Guide Star Adaptive Optics (LGSAO) astrometric measurements of stars near the Galaxy's supermassive black hole (SMBH). By carrying out a systematic study of M92, we have produced a more accurate model for the camera's optical distortion. Updating our measurements with this model, and accounting for differential atmospheric refraction, we obtain estimates of the SMBH properties that are a factor of $\sim$2 more precise, and most notably, increase the likelihood that the black hole is at rest with respect to the nuclear star cluster. These improvements have also allowed us to extend the radius to which stellar orbital parameter estimates are possible by a factor of 2.

preprint2010arXiv

Recent Results and Perspectives for Precision Astrometry and Photometry with Adaptive Optics

Large ground-based telescopes equipped with adaptive optics (AO) systems have ushered in a new era of high-resolution infrared photometry and astrometry. Relative astrometric accuracies of <0.2 mas have already been demonstrated from infrared images with spatial resolutions of 55-95 mas resolution over 10-20'' fields of view. Relative photometric accuracies of 3% and absolute photometric accuracies of 5%-20% are also possible. I will review improvements and current limitations in astrometry and photometry with adaptive optics of crowded stellar fields. These capabilities enable experiments such as measuring orbits for brown dwarfs and exoplanets, studying our Galaxy's supermassive black hole and its environment, and identifying individual stars in young star clusters, which can be used test the universality of the initial mass function.

preprint2010arXiv

Testing stellar cusp formation theories with observations of the Milky Way nuclear star cluster

(Abridged) We report on the structure of the nuclear star cluster in the innermost 0.16 pc of the Galaxy as measured by the number density profile of late-type giants. Using laser guide star adaptive optics in conjunction with the integral field spectrograph, OSIRIS, at the Keck II telescope, we are able to differentiate between the older, late-type (~ 1 Gyr) stars, which are presumed to be dynamically relaxed, and the unrelaxed young (~ 6 Myr) population. This distinction is crucial for testing models of stellar cusp formation in the vicinity of a black hole, as the models assume that the cusp stars are in dynamical equilibrium in the black hole potential. We find that contamination from young stars is significant, with more than twice as many young stars as old stars in our sensitivity range (K < 15.5) within the central arcsecond. Based on the late-type stars alone, the surface stellar number density profile, is flat, with a projected power law slope of -0.26+-0.24. These results are consistent with the nuclear star cluster having no cusp, with a core profile that is significantly flatter than predicted by most cusp formation theories. Here, we also review the methods for further constraining the true three-dimensional radial profile using kinematic measurements. Precise acceleration measurements in the plane of the sky as well as along the line of sight has the potential to directly measure the density profile to establish whether there is a "hole" in the distribution of late-type stars in the inner 0.1 pc.

preprint2009arXiv

High angular resolution integral-field spectroscopy of the Galaxy's nuclear cluster: a missing stellar cusp?

We report on the structure of the nuclear star cluster in the innermost 0.16 pc of the Galaxy as measured by the number density profile of late-type giants. Using laser guide star adaptive optics in conjunction with the integral field spectrograph, OSIRIS, at the Keck II telescope, we are able to differentiate between the older, late-type ($\sim$ 1 Gyr) stars, which are presumed to be dynamically relaxed, and the unrelaxed young ($\sim$ 6 Myr) population. This distinction is crucial for testing models of stellar cusp formation in the vicinity of a black hole, as the models assume that the cusp stars are in dynamical equilibrium in the black hole potential. Based on the late-type stars alone, the surface stellar number density profile, $Σ(R) \propto R^{-Γ}$, is flat, with $Γ= -0.27\pm0.19$. Monte Carlo simulations of the possible de-projected volume density profile, n(r) $\propto r^{-γ}$, show that $γ$ is less than 1.0 at the 99.73 % confidence level. These results are consistent with the nuclear star cluster having no cusp, with a core profile that is significantly flatter than predicted by most cusp formation theories, and even allows for the presence of a central hole in the stellar distribution. Of the possible dynamical interactions that can lead to the depletion of the red giants observable in this survey -- stellar collisions, mass segregation from stellar remnants, or a recent merger event -- mass segregation is the only one that can be ruled out as the dominant depletion mechanism. The lack of a stellar cusp around a supermassive black hole would have important implications for black hole growth models and inferences on the presence of a black hole based upon stellar distributions.

preprint2008arXiv

A near-IR variability study of the Galactic black hole: a red noise source with no detected periodicity

We present the results of near-infrared (2 and 3 microns) monitoring of Sgr A*-IR with 1 min time sampling using the natural and laser guide star adaptive optics (LGS AO) system at the Keck II telescope. Sgr A*-IR was observed continuously for up to three hours on each of seven nights, between 2005 July and 2007 August. Sgr A*-IR is detected at all times and is continuously variable, with a median observed 2 micron flux density of 0.192 mJy, corresponding to 16.3 magnitude at K'. These observations allow us to investigate Nyquist sampled periods ranging from about 2 minutes to an hour. Using Monte Carlo simulations, we find that the variability of Sgr A* in this data set is consistent with models based on correlated noise with power spectra having frequency dependent power law slopes between 2.0 to 3.0, consistent with those reported for AGN light curves. Of particular interest are periods of ~20 min, corresponding to a quasi-periodic signal claimed based upon previous near-infrared observations and interpreted as the orbit of a 'hot spot' at or near the last stable orbit of a spinning black hole. We find no significant periodicity at any time scale probed in these new observations for periodic signals. This study is sensitive to periodic signals with amplitudes greater than 20% of the maximum amplitude of the underlying red noise component for light curves with duration greater than ~2 hours at a 98% confidence limit.