We describe development and application of a Global Astrometric Solution (GAS) to the problem of Pan-STARRS1 (PS1) astrometry. Current PS1 astrometry is based on differential astrometric measurements using 2MASS reference stars, thus PS1 astrometry inherits the errors of the 2MASS catalog. The GAS, based on a single, least squares adjustment to approximately 750k grid stars using over 3000 extragalactic objects as reference objects, avoids this catalog-to-catalog propagation of errors to a great extent. The GAS uses a relatively small number of Quasi-Stellar Objects (QSOs, or distant AGN) with very accurate (<1 mas) radio positions, referenced to the ICRF2. These QSOs provide a hard constraint in the global least squares adjustment. Solving such a system provides absolute astrometry for all the stars simultaneously. The concept is much cleaner than conventional astrometry but is not easy to perform for large catalogs. In this paper we describe our method and its application to Pan-STARRS1 data. We show that large-scale systematic errors are easily corrected but our solution residuals for position (~60 mas) are still larger than expected based on simulations (~10 mas). We provide a likely explanation for the reason the small-scale residual errors are not corrected in our solution as would be expected.

CMOS-hybrid arrays have recently surfaced as competitive optical detectors for use in ground- and space-based astronomy. One source of error in these detectors that does not appear in more traditional CCD arrays is the inter-pixel capacitance component of crosstalk. In this paper we use a single pixel reset method to model inter-pixel capacitance (IPC). We combine this IPC model with a model for charge diffusion to estimate the total crosstalk on H4RG arrays. Finally, we compare our model results to Fe55 data obtained using an astrometric camera built to test the H4RG-B0 generation detectors.

The non-orthogonality of algebraic polynomials of field coordinates traditionally used to model field-dependent corrections to astrometric measurements, gives rise to subtle adverse effects. In particular, certain field dependent perturbations in the observational data propagate into the adjusted coefficients with considerable magnification. We explain how the worst perturbation, resulting in the largest solution error, can be computed for a given non-orthogonal distortion model. An algebraic distortion model of full rank can be converted into a fully orthonormal model based on the Zernike polynomials for a circular field of view, or a basis of functions constructed from the original model by a variant of the Gram-Schmidt orthogonalization process for a rectangular field of view. The relative significance of orthonormal distortion terms is assessed simply by the numerical values of the corresponding coefficients. Orthonormal distortion models are easily extendable when the distribution of residuals indicate the presence of higher order terms.