preparations for #LGAstat

At an undisclosed location in the Catskills, I worked on getting a talk and slides ready for the Local Group Astrostatistics meeting next week in Ann Arbor. Since the whole Stream Team will be there to talk about our stuff on streams, I plan on talking about Bovy and my April-Fools paper on the Solar System (and what that means for dynamical inference in the Milky Way) and The Cannon (because I want to encourage the crew to think about modeling the stars, too). I am not ready for my talk yet!

noisy photometric redshift posteriors

The close reader (there's a reader, I know, but a close reader?) of this blog will have noticed that CampHogg is in a transition. We used to think that projects producing "catalogs" from telescope data ought to produce likelihood information. Now we think this is probably impossible in general, and we will have to live with (at best) posterior probability information, under some priors. We discussed this in group meeting, in particular related to Malz's project and LSST. The photometric redshift system LSST expects to create using cross-correlations (with, say, quasars) will (if the method works) produce posterior probability distribution function estimates (not values, but estimates, which is scary).

Key questions we identified include the following: Understand what the effective priors are for that procedure. Understand what the deep assumptions are behind the cross-correlations; I think they have to be at large scales to work properly (linear bias and all that). And understand whether it has been demonstrated to work, empirically. I got all confused at the end of group meeting about the fact that the method generates noisy estimates of a pdf, which is a strange meta issue. What do you do about a noisily known posterior pdf?

classification, collisions, and calibration

José Garmilla (Princeton) showed up for group meeting, and talked to us about star–galaxy separation, which is of great interest to us at CampHogg. He has been trying many things—including many of the methods we at CampHogg been advocating—all with HyperSuprimeCam data, and he has been hitting many of the things Fadely and I have been hitting over the years: Things break down at faint magnitudes, even the existence of a truth table or training set is hard to establish once you get below about 25th magnitude.

In that same group meeting, Blanton introduced the idea of SDSS fiber collisions and the methods that are used to mitigate them. This is a seemingly trivial question: What to do with the galaxies that don't get redshifts because of proximity to another galaxy that does get a redshift? But when it comes to two-point statistics, there is no trivial answer.

After lunch, Kopytova and I pair-coded some of her spectral-fitting code, implementing the ideas we had about fitting simultaneously for a bunch of calibration or normalization parameters along with the stellar parameters. The code we (very hastily) wrote seems to work well, but of course it has model-complexity parameters (or, really, hyperparameters). That's to think about tomorrow!

nuthin'

All administration, all day; no research to speak of! Except a bit of loose talk at lunch with CampHogg.

long-period exoplanets

My research day started late morning, when I met with Heather Knutson (Caltech), Marta Bryan (Caltech), and Henry Ngo (Caltech) to discuss population inferences involving long-period Jupiter-like planets. For each star, they get a velocity trend ("acceleration", I would call it) and a curvature ("jerk", I would call it), and they use these to constrain the properties of any long-period companion. We discussed what they are doing, and also my ideal vision of what is the best thing that can be done with these data.

One interesting technical point: When you have an acceleration and jerk that are both consistent with zero, do you treat them as "non-detections" or do you treat them identically with detections, and just compute the likelihood of the data given a long-period companion model? In the context of a full hierarchical model, there is no need (I think) to treat the detections and the non-detections differently. The key thing is to permit the model to include zero-mass planets or non-planets or models of differing complexity.

spectral asteroseismology

In my small amount of research time today I worked on a title, abstract, and introduction for a possible paper with Ness on inferring asteroseismological parameters (nu-max and delta-nu) from stellar spectra, using a modification of The Cannon. The project started when a preliminary study suggested that we can get ages. That is, we can train with stars with "known" ages and then infer the ages of other stars. Since the ages ultimately come from asteroseismological parameters, it makes sense to see if we can infer them first. If we can, it might have a big impact on what we think influences stellar spectra.

#ArloFest, day 1

Today was the first day of Landolt Standards & 21st Century Photometry in Baton Rouge, organized by Pagnotta (AMNH) and Clayton (LSU). I came to speak about self-calibration. The day started with a historical overview by Bessel (MSSSO), who gave a lovely talk filled with profiles of the many people who contributed to the development of photometric calibration and magnitude systems. Many of the people he talked about (including himself) have filter systems or magnitude systems named after them! Among the many interesting things he touched on was this paper by Johnson, which I have yet to carefully read, but apparently contains some of the philosophy behind standard-star systems. He also discussed the filter choices for the Skymapper project, which seem very considered.

Suntzeff (TAMU) gave an excellent talk about the limitations of the supernova cosmology projects; his main point is that systematic issues with the photometric calibration system are the dominant term in the uncertainty budget. This is important in thinking about where to apportion new resources. He made a great case for understanding physically every part of the photometric measurement system (and that includes the stars, the atmosphere, the telescope, and the detector pixels, among other things). I couldn't agree more!

Grindlay (CfA) blew us away with the scale and content of the DASCH plate-scanning project at Harvard. It is just awesome, in time span, cadence, and sky coverage. Anyone not searching these data is making a mistake! And, as we were recovering from that, Kafka (AAVSO) blew us away again with the scale and scope of the APASS survey, which was designed, built, operated, reduced, and delivered to the public almost entirely by citizen scientists. It is dramatic; we are not worthy!

There were many other great contributions—too many to mention them all—but the day ended with a crawfish boil and then Josh Peek (STScI) and I at the bar discussing recent explosive conversations in the astronomical community around TMT and development in Hawaii.

One last thing I should say: Arlo Landolt (LSU) has had a huge impact on astronomy; his work has enabled countless projects and scientific measurements and discoveries. The development and stewardship of photometric standards and systems, and all the attention to detail it requires, is unglamorous and time-consuming work, ill-suited to most of the community, and yet absolutely essential to everything we do. I can't thank Landolt—and his collaborators and the whole community of photometrists—enough.

self-calibration for Arlo

I got no research time in today until I headed to the airport, where I worked on my presentation to the Landolt Standards & 21st Century Photometry Symposium in honor of Arlo Landolt (LSU). I will be talking about self-calibration, the set of methods that have replaced Landoldt's (amazingly productive and important) standard-star catalogs for the calibration of huge surveys, which can (in principle) now be self-calibrated. I will talk about both the theory and the practice of the self-calibration arts.

probabilistic cosmology

Boris Leistedt (UCL) showed up for two days of chatting in preparation for his arrival at NYU this coming Fall. We discussed a set of projects in probabilistic cosmology. In one (which I have discussed previously with Fadely), we imagine what it would look like to infer galaxy redshifts from imaging without either a training set or models of galaxy spectral energy distributions. I feel like it might be possible, with the thought experiment being: What if you were given the photometry of the SDSS Luminous Red Galaxies? Couldn't you figure out their redshifts without being told their underlying spectral energy distributions? At least up to some geometric factor? Leistedt predicts that any method that doesn't use either models or training must do worse than any method that has good training data or models. However, it might make a framework for then including both model (theory) and training (data) information in in a principled way.

In another class of subjects, we talked about inference of the density field or the initial conditions. In the first place, if we could infer the density field in patches, we could use that to inform photometric redshifts or reconstruction. In the second, we could perhaps infer the initial conditions (and I mean the phases, not just the power spectrum); this is interesting because the Baryon Acoustic Feature ought to be much stronger and sharper in the initial conditions than it is today. We discussed some conceivably tractable approaches, some based on likelihood-free inference, some based on Gaussian Processes, and some based on machine learning (using simulations as training data!).

exoplanet dynamics, sucking

In group meeting, Dun Wang talked about astrometric calibration of the GALEX Satellite, and Kat Deck (Caltech) talked about the dynamical evolution of exoplanetary systems. She pointed out that we naively expect lots of planets close in period to be locked in resonances, but in fact such resonances are rare, empirically in the Kepler sample. She has explanations for this involving the evolving proto-planetary disk.

After lunch, Deck gave the astro seminar, on planetary system stability and the Kepler planets. She discussed chaos, stability, and heuristic stability criteria. One interesting thing is that there really is no non-heuristic stability criterion: We think of a planetary system as "stable" if there are no catastrophic, order-unity changes to any of the orbital osculating elements. That's not really an equation! And at the talk there was some discussion of the point (counter-intuitive and important) that a system can be stable (by our astronomer definition) for far, far longer than the Lyapunov time. Awesome and important.

At the end of the day, Foreman-Mackey and I made the (astonishing) decision to abort and fail on our NASA proposal: We just ran out of time. I am disappointed; we have an eminently fundable pitch. That said, we just didn't start early enough to make that pitch at the level we wanted. Not sure how to feel about it, but I sure need to catch up on sleep!

finding planets with TTVs

In the tiny bit of research time available today, I spoke with Kat Deck (Caltech) and Foreman-Mackey about finding planets with large transit-timing residuals. These signals aren't precisely periodic, so some searches could miss them entirely. Deck has simple models (based on perturbation theory) for the variations, so we can in principle add only a few parameters and capture a large range in possible variations. This might make us much more sensitive to extremely important and interesting systems hiding in existing data.

The rest of the day was spent writing our NASA proposal, with the exception of lunch with Yann LeCun (facebook) and Léon Bottou (facebook) at the NYC office of facebook Research.

ultra-short-period planets, 2

At group meeting, Roberto Sanchis-Ojeda (Berkeley) and Kat Deck (Caltech) were around to argue exoplanets. We had Sanchis-Ojeda tell us about his ultra-short-period (USP) planet work. He showed us his evidence for a disrupting (or something) USP planet; to me it looks like evidence for a technological civilization making power from a star! But, for some reason, no-one else thought so! We discussed tidal mechanisms and stripping mechanisms, but we concluded that there is no real reason to think of the USPs as being related to hot jupiters or any other kind of planet in any mechanistic way.

ultra-short-period planets

Foreman-Mackey and I spent quality time with Roberto Sanchis-Ojeda (Berkeley) figuring out what we are going to do together, related to ultra-short-period planets. He was pulling for things related to false-positive evaluation, I was arguing for search at lower signal-to-noise ratios. I feel like there are lots of undiscovered systems. The rest of our science time was spent proposal-writing for NASA!

talks all day; Dr Gorbenko

Roberto Sanchis-Ojeda (Berkeley) arrived for a few days of hacking. He is responsible for finding some very short-period transiting exoplanets. We didn't get a chance to set our goals for the day, because the day was packed with talks:

Victor Gorbenko (NYU) gave an absolutely excellent PhD defense. Congratulations Dr Gorbenko. Among other things, he was considering the difference between the theoretical (analytic) action of massless modes on the 1+1 worldsheet of a string (not a 0+1 wordline but s 1+1 worldsheet) and a lattice calculation. In this comparison, he found massive modes. He very nicely connected this work to larger questions in string theory, but also to the origins of string theory in QCD.

At lunch, Guido D'Amico (NYU) spoke about axions, axion–photon interactions, and cosmic transparency. This is related to work I did with Bovy and with More in the optical. He advertised some clever ideas he has about about using the Planck data to constrain this sector using millimeter wavelengths. He has been implementing some of these ideas at #astrohackny.

At the end of the day, Neil Lawrence (Sheffield) spoke about deep machine-learning models built from nested Gaussian Processes. The methods show a lot of promise and connect well to things we are thinking about in CampHogg.

JPL, day 3

I spent the morning with Leonidas Moustakas (JPL), a bit of Foreman-Mackey (in Pasadena for the Sagan Fellowship Symposium), and a bit of Andrew Romero-Wolf (JPL) and Curtis McCully (LCOGT) by phone, discussing projects related to strong-lensing time-delay measurements. We discussed two challenging projects. The first is to determine (from as many as we can construct) the best model for quasar time-domain variability. There are claims in the literature that the damped random walk is the best model, but that (very sensible) model hasn't really been competed against all that much. We know how to do this, lets do this! The reason they want good probabilistic generative models is that they want to determine time delays as precisely as possible, using a probabilistic approach.

The second project is to perform high-quality photometry on the (overlapping from the ground) images of a strongly lensed quasar. In this case—when the point-spread functions overlap—you have to do your photometry by simultaneous fitting, but with the variable (and badly known) PSF of ground-based astronomy, I have never seen such photometry that really looks right: There are always fitting-induced covariances of the overlapping-source light curves. I think this is caused by model mismatch (under-fitting), but I don't really know. Romero-Wolf and McCully pointed out that image differencing methods work well in crowded fields, so I formulated an image-modeling approach to the photometry that is as close to image differencing as possible. I promised to write it up into a document. I am kind-of excited about it; it is still just image modeling, but it makes use of the power of image differencing technologies to get flexibility to fit the real PSF as it is.

Late in the day, Adam Miller (JPL) showed me the JPL Mission Control center and a few high bays, filled with awesome stuff (including some fake Mars!).