disrupted and perturbed

I spent much of the day working on an old-school perturbation-theory calculation. I consider a cold tidal stream in a host galaxy potential, perturbed by the close passage of a point mass. In the limit of small perturbations of very cold streams, this calculation has only two free parameters: the angle between the direction of the stream and the velocity of the perturber (in the comoving frame of the stream or equivalent), and the time since the impulse (in some scaled time units related to the mass and velocity of the perturber or amplitude of the perturbation). This is all very idealized, but actually, to zeroth order, I think I may have exhaustively described all possible perturbations to cold streams. Now the idea is to use this to constrain the substructure in our Galaxy with observations of cold streams.

finished writing

I finished the faint-source proper-motion paper. It still needs to be vetted by collaborators, but I am stoked. Here is the abstract:

The near future of astrophysics involves many large solid-angle, multi-epoch, multi-band imaging surveys. These surveys will, at their faint limits, have data on large numbers of sources that are too faint to detect at any individual epoch. Here we show that it is possible to measure in multi-epoch data not only the fluxes and positions, but also the parallaxes and proper motions of sources that are too faint to detect at any individual epoch. The method involves fitting a model of a moving point source simultaneously to all imaging, taking account of the noise and point-spread function in each image. By this method it is possible—in well-understood data—to measure the proper motion of a point source with an uncertainty (found after marginalizing over flux, mean position, and parallax) roughly equal to the minimum possible uncertainty given the information in the data, which is limited by the point-spread function, the distribution of observation times, and the total signal-to-noise in the combined data. We demonstrate our technique on artificial data and on multi-epoch Sloan Digital Sky Survey imaging of the SDSS Southern Stripe. With the SDSSSS data we show that with this technique it is possible to distinguish very red brown dwarfs from very high-redshift quasars and from resolved galaxies more than 1.6 mag fainter than by the traditional technique. Proper motions distinguish faint brown dwarfs from faint quasars with better fidelity than multi-band imaging alone; we present 16 new candidate brown dwarfs in the SDSSSS, identified on the basis of high proper motion. They are likely to be halo stars because none has a significantly measured parallax.

post-starbursts in GALEX, writing

Continued writing on the faint-motion paper. It is funny how much there is left to do when a project is done!

Wu began working on the GALEX properties of post-starburst galaxies: They aren't detected at the MIS depth. That is good, because star-forming galaxies are detected. So a GALEX selection is likely to work, at some level. The question is, how good will it be? We would like to have a GALEX-based selection of the post-starburst galaxies so we can perform emission-line studies without worrying about the fact that the post-starbursts are selected on the basis of emission lines. Fortuitously, Christy Tremonti (Arizona) showed up at the MPIA today for a month, so she may be able to help out.

writing not finished

I tried, but I didn't finish the faint-proper-motion paper today.

submission and resubmission

Bovy, with some help from Moustakas and me, got ready the galaxy-cluster transparency paper for resubmission in response to referee. The referee really made a big difference to the paper, because he or she recommended averaging the samples in a better way, which improved the results. I, with help from Barron and Roweis, got ready the Blind Date paper (on estimating image dates using proper motions) for resubmission. And I, with help from Lang, have promised my co-authors I will get the paper on faint-source proper motions ready for submission by the end of the week. That end approaches fast.

bimodality, transparency

I switched my search for low kurtosis directions in spectrum space into a search for bimodal directions. That is, I wrote down a scalar (which has to do with k-means with k=2) that decreases as a distribution becomes more bimodal. Then I searched Tsalmantza's high-variance PCA components for directions in the space that are most bimodal. I find three, perpendicular bimodal directions! Of course each one is a different version of the red–blue galaxy bimodality, of which I have been an unheard critic. More on this as I understand it better.

Surhud More (MPIA) and I began working this week on the monopole term in the opacity of the Universe, using the consistency of baryon-acoustic and supernovae measures of the expansion history to check the phase-space conservation of photons. This test is (nearly) independent of world model, as it depends almost entirely on purely special-relativistic considerations. We hope to have an LPU (least publishable unit) on the subject soon.

null hypothesis

Wu and I are back onto looking at the processes that lead to post-starburst galaxies, this time with Frank van den Bosch and Anna Pasquali here at MPIA. The first step was to create comparison samples to make null hypotheses; because the catalog we are using has redshift and flux dependencies, we built comparison samples to be exactly matched in redshift and brightness (stellar mass). Wu finished that today.

kurtosis minimization

With Tsalmantza's help, I got the kurtosis minimization working on high-variance directions in the SDSS spectral space. In the minimal kurtosis directions (within the high-variance subspace), the star-forming and non-star-forming galaxies separate very clearly, and there are other tantalizing structures. I think this technique may have legs.

blind date

I worked on the response to referee for our blind date paper. We should be able to resubmit this week.

finishing paper

Lang and I discussed finishing the faint-source proper-motion paper, among other things.

kurtosis

A few days ago I bashed PCA on various grounds, in particular that it ranks components by their contribution to the data variance, and it is rarely the data variance about which one cares. Today in discussions with Tsalmantza I realized that one could rank components by the kurtosis of their amplitudes (rather than the variance), and lowest first. This has a number of advantages, but one is that (uninteresting) data artifacts and outliers tend to create high-kurtosis directions in data space, and another is that if there are directions that are multi-modal, they tend also to be low in kurtosis (think the color distribution of galaxies, which is bimodal and low in kurtosis). It is still a very frequentist approach, but a search for minimal kurtosis directions in data space might be productive. Tsalmantza and I hope to give it a shot next week.

cross-correlation issues

I made plots of my quasar–photon and white-dwarf–photon cross-correlations, and the random samples which are supposed to be equivalent. There must be some kind of bug, because the random has a negative feature at the center! So I will spend the rest of this week de-bugging.

mean white dwarf in GALEX

Schiminovich and I hope to detect intergalactic scattering with our quasar–photon cross-correlations. In order to make this detection, which will require precision, we need to create a differential experiment. The first difference is between mean quasars and mean white dwarfs. The white dwarfs are so close, they should have essentially no scattering (or scattering local to the observatory that is shared with the quasars). I made the white-dwarf–photon cross-correlations this weekend.

The mean white-dwarf image is probably interesting in its own right, if I broke down the white dwarfs by type and temperature, because they would provide extremely high signal-to-noise GALEX information. Does anyone want that?

cross-correlate and visualize

Wrote code to combine the jackknife subsamples of my quasar–photon cross-correlation functions and visualize the output. The redshift-dependence of the ultraviolet flux from quasars was not as strong as I expected. Either I have a bug in my code or a bug in my thinking.

ultraviolet vs redshift

I got back up to speed on the quasar-photon cross-correlation, splitting the SDSS photometric quasar sample by redshift. This is the first test: Does the ultraviolet flux depend on redshift as we expect? It ought to drop out at the redshifts at which the two GALEX bandpasses cross Lyman-alpha and the Lyman limit. Hope to have results for tomorrow.