Fast and precise propagation of satellite orbits is required for mission design, orbit determination in support of operations and payload data analysis. This demand must also comply with the different accuracy requirements set by a growing variety of scientific and service missions. This contribution proposes a method to improve the computational performance of orbit propagators through an efficient numerical integration that meets the accuracy requirements set by the specific application. This is achieved by appropriately tuning the parameters of the numerical propagator (relative tolerance and maximum time step), establishing a threshold for the perturbing accelerations (Earth's gravitational potential, atmospheric drag, solar radiation pressure, third-body perturbations, relativistic correction to gravity) below which they can be neglected without altering the quality of the results and implementing an efficient and precise algorithm for the harmonic synthesis of the geo-potential and its first-order gradient. In particular, when performing the harmonic synthesis, the number of spherical harmonics to retain (i.e., the expansion degree) is determined by the accuracy requirement. Given that higher-order harmonics decay rapidly with altitude, the expansion degree necessary to meet the target accuracy decreases with height. To improve the computational efficiency, the number of degrees to retain is determined dynamically while the trajectory is being computed. The optimum expansion degree for each altitude is determined by ensuring that the truncation error of the harmonic synthesis is below the threshold acceleration. The work is a generalization to arbitrary orbits of a previous study that focused on communication satellites in geosynchronous inclined orbits. The method is presented and a set of test cases is analysed and discussed.

The Tisserand graph (TG) is a graphical tool commonly employed in the preliminary design of gravity-assisted trajectories. The TG is a two-dimensional map showing essential orbital information regarding the Keplerian orbits resulting from the close passage by one or more massive bodies, given the magnitude of the hyperbolic excess speed ($v_{\infty}$) and the minimum allowed pericenter height for each passage. Contours of constant $v_{\infty}$ populate the TG. Intersections between contours allow to link consecutive flybys and build sequences of encounters en route to a selected destination. When the number of perturbing bodies is large and many $v_{\infty}$ levels are considered, the identification of all the possible sequences of encounters through the visual inspection of the TG becomes a laborious task. Besides, if the sequences are used as input for a numerical code for trajectory design and optimization, an automated examination of the TG is desirable. This contribution describes an automatic technique to explore the TG and find all the encounter paths. The technique is based on a tree search method, and the intersections between contours are found using the regula-falsi scheme. The method is validated through comparisons with solutions available in the open literature. Examples are given of application to interplanetary mission scenarios, including the coupling with a trajectory optimizer.

We present photometric and spectroscopic observations of SN 2013aa and SN 2017cbv, two nearly identical type Ia supernovae (SNe Ia) in the host galaxy NGC 5643. The optical photometry has been obtained using the same telescope and instruments used by the Carnegie Supernova Project. This eliminates most instrumental systematics and provides light curves in a stable and well-understood photometric system. Having the same host galaxy also eliminates systematics due to distance and peculiar velocity, providing an opportunity to directly test the relative precision of SNe Ia as standard candles. The two SNe have nearly identical decline rates, negligible reddening, and remarkably similar spectra and, at a distance of $\sim 20$ Mpc, are ideal as potential calibrators for the absolute distance using primary indicators such as Cepheid variables. We discuss to what extent these two SNe can be considered twins and compare them with other supernova "siblings" in the literature and their likely progenitor scenarios. Using 12 galaxies that hosted 2 or more SNe~Ia, we find that when using SNe~Ia, and after accounting for all sources of observational error, one gets consistency in distance to 3 percent.

Orbit insertion at Saturn requires a large impulsive manoeuver due to the velocity difference between the spacecraft and the planet. This paper presents a strategy to reduce dramatically the hyperbolic excess speed at Saturn by means of deep-space electric propulsion. The interplanetary trajectory includes a gravity assist at Jupiter, combined with low-thrust maneuvers. The thrust arc from Earth to Jupiter lowers the launch energy requirement, while an ad hoc steering law applied after the Jupiter flyby reduces the hyperbolic excess speed upon arrival at Saturn. This lowers the orbit insertion impulse to the point where capture is possible even with a gravity assist with Titan. The control-law algorithm, the benefits to the mass budget and the main technological aspects are presented and discussed. The simple steering law is compared with a trajectory optimizer to evaluate the quality of the results and possibilities for improvement.

The Dark Energy Camera is a new imager with a 2.2-degree diameter field of view mounted at the prime focus of the Victor M. Blanco 4-meter telescope on Cerro Tololo near La Serena, Chile. The camera was designed and constructed by the Dark Energy Survey Collaboration, and meets or exceeds the stringent requirements designed for the wide-field and supernova surveys for which the collaboration uses it. The camera consists of a five element optical corrector, seven filters, a shutter with a 60 cm aperture, and a CCD focal plane of 250 micron thick fully-depleted CCDs cooled inside a vacuum Dewar. The 570 Mpixel focal plane comprises 62 2kx4k CCDs for imaging and 12 2kx2k CCDs for guiding and focus. The CCDs have 15 microns x15 microns pixels with a plate scale of 0.263 arc sec per pixel. A hexapod system provides state-of-the-art focus and alignment capability. The camera is read out in 20 seconds with 6-9 electrons readout noise. This paper provides a technical description of the camera's engineering, construction, installation, and current status.

(Abridged) We study predictions for galaxy cluster observables that can test the statistics of dark matter halo shapes expected in a flat LCDM universe. We present a simple analytical model for the prediction of cluster-scale X-ray observations, approximating clusters as isothermal systems in hydrostatic equilibrium, and dark matter haloes as ellipsoids with uniform axial ratios. We test the model against high-resolution, hydrodynamic cluster simulations to gauge its reliability. We find that this simple prescription does a good job of predicting the distribution of cluster X-ray ellipticities compared to the simulations as long as one focuses on cluster regions that are less sensitive to recent mergers. Based on this simple model, the distribution of cluster-size halo shapes expected in the concordance LCDM cosmology implies an X-ray ellipticity distribution with a mean of 0.32 +- 0.01 and a scatter of 0.14 +- 0.01 for the mass range (1-4)x10^{14} Msun/h. We find it important to include the mass dependence of halo shape to make comparisons to observational samples that contain many, very massive clusters. We analyse the systematics of four observational samples of cluster ellipticities and find that our results are statistically compatible with observations. In particular, we find remarkably good agreement between two recent ROSAT samples and LCDM predictions that DO NOT include gas cooling. We also test how well our analytical model can predict Sunyaev-Zel'dovich decrement maps and find that it is less successful although still useful; the model does not perform as well as a function of flux level in this case because of the changing triaxiality of dark matter haloes as a function of radial distance. Both this effect and the changing alignment of isodensity shells of dark matter haloes leave an imprint on cluster gas...

Using six high resolution dissipationless simulations with a varying box size in a flat LCDM universe, we study the mass and redshift dependence of dark matter halo shapes for M_vir = 9.0e11 - 2.0e14, over the redshift range z=0-3, and for two values of sigma_8=0.75 and 0.9. Remarkably, we find that the redshift, mass, and sigma_8 dependence of the mean smallest-to-largest axis ratio of halos is well described by the simple power-law relation <s> = (0.54 +- 0.02)(M_vir/M_*)^(-0.050 +- 0.003), where s is measured at 0.3 R_vir and the z and sigma_8 dependences are governed by the characteristic nonlinear mass, M_*=M_*(z,sigma_8). We find that the scatter about the mean s is well described by a Gaussian with sigma ~ 0.1, for all masses and redshifts. We compare our results to a variety of previous works on halo shapes and find that reported differences between studies are primarily explained by differences in their methodologies. We address the evolutionary aspects of individual halo shapes by following the shapes of the halos through ~100 snapshots in time. We determine the formation scalefactor a_c as defined by Wechsler et al. (2002) and find that it can be related to the halo shape at z = 0 and its evolution over time.

We present a numerical Hydro+N-body model of A3266 whose X-ray surface brightness, temperature distribution, and galaxy spatial and velocity distribution data are consistent with the A3266 data. The model is an old (~3 Gyr), off-axis merger having a mass ratio of ~2.5:1. The less massive subcluster in the model is moving on a trajectory from southwest to northeast passing on the western side of the dominant cluster while moving into the plane of the sky at ~45 degrees. Off-axis mergers such as this one are an effective mechanism for transferring angular momentum to the intracluster medium (ICM), making possible a large scale rotation of the ICM. We demonstrate here that the ICM rotation predicted by our fully 3-dimensional model of A3266 is observable with current technology. As an example, we present simulated observations assuming the capabilities of the high resolution X-ray spectrometer (XRS) which was to have flown on Astro-E.

The degeneracy between the disk and the dark matter contribution to galaxy rotation curves remains an important uncertainty in our understanding of disk galaxies. Here we discuss a new method for breaking this degeneracy using gravitational lensing by spiral galaxies, and apply this method to the spiral lens B1600+434 as an example. The combined image and lens photometry constraints allow models for B1600+434 with either a nearly singular dark matter halo, or a halo with a sizable core. A maximum disk model is ruled out with high confidence. Further information, such as the circular velocity of this galaxy, will help break the degeneracies. Future studies of spiral galaxy lenses will be able to determine the relative contribution of disk, bulge, and halo to the mass in the inner parts of galaxies.

The frequency with which background galaxies appear as long arcs as a result of gravitational lensing by foreground clusters of galaxies has recently been found to be a very sensitive probe of cosmological models by Bartelmann et al. (1998). They have found that such arcs would be expected far less frequently than observed (by an order of magnitude) in the currently favored model for the universe, with a large cosmological constant $\Omega_\Lambda \sim 0.7$. Here we analyze whether including the effect of cluster galaxies on the likelihood of clusters to generate long-arc images of background galaxies can change the statistics. Taking into account a variety of constraints on the properties of cluster galaxies, we find that there are not enough sufficiently massive galaxies in a cluster for them to significantly enhance the cross section of clusters to generate long arcs. We find that cluster galaxies typically enhance the cross section by only $\lesssim 15%$.

We present results of simple N-body simulations that strengthen the suggestion that A3266 is composed of two subunits of comparable mass that have recently merged. Both the real cluster and the N-body dark-matter cluster show mixed signals of substructure under statistical tests. However, in a decidedly non-statistical approach allowed by the wide-area coverage and large number of redshifts they measured in A3266, Quintana, Ramirez, and Way (1996; QRW) sliced the real cluster in redshift space to uncover a peculiar spatial distribution of galaxies that they suggested was the result of a recent merger. In our simulations, a similar distribution is the result of an ongoing merger between two comparable-mass units that started about 2x10^9 years ago in the N-body simulations. We also find that the distribution of emission line galaxies in A3266 traces the same structure. We discuss further tests of our merger hypothesis, and speculate on the possibility that a similar process might be occurring in other, apparently-relaxed clusters at the present epoch.

Using spectrophotometry of a large sample of galaxies in 19 Abell Clusters, we have selected 42 candidate AGN by the criteria used by Dressler et.al. 1985 (DTS) in their analysis of the statistics of 22 AGN in 14 rich cluster fields, which are based on the equivalent width of [OII]3727, Hbeta, and [OIII]5007 emission. We have then discriminated AGN from HII region-like galaxies (hereafter HII galaxies) in the manner developed by Veilleux & Osterbrok 1987 using the additional information provided by Halpha and [NII]6583 or Halpha and [SII]6716+6731 emission, in order to test the reliability of the selection criteria used by DTS. We find that our sample is very similar to that of DTS before we discriminate AGN from HII galaxies, and would lead to similar conclusions. However, we find that their method inevitably mixes HII galaxies with AGN, even for the most luminous objects in our sample. Other authors have attempted to quantify the relative fraction of cluster to field AGN since the study of DTS (Hill & Oegerle 1993, Biviano et al. 1997) reaching similar conclusions, but using similar criteria to theirs to select AGN (or using the [OIII]5007/Hbeta flux ratio test that also mixes HII galaxies with AGN). Our sample of true AGN is left too small to reach statistically meaningful conclusions, therefore a new study with the more time-consuming method that includes the other lines will be required to quantify the true relative fraction of cluster to field AGN.

Spheroidal components of spiral galaxies have been considered the only dynamically important component in gravitational lensing studies thus far. Here we point out that including the disk component can have a significant effect, depending on the disk inclination, on a variety of lensing properties that are relevant to present studies and future surveys. As an example, we look at the multiple image system B1600+434, recently identified as being lensed by a spiral galaxy. We find that including the disk component one can understand the fairly large image separation as being due to the inclination of a typical spiral, rather than the presence of a very massive halo. The fairly low magnification ratio can also be readily understood if the disk is included. We also discuss how such lensed systems might allow one to constrain parameters of spiral galaxies such as a disk-to-halo mass ratio, and disk mass scale length. Another example we consider is the quasar multiple-lensing cross section, which we find can increase many-fold at high inclination for a typical spiral. Finally, we discuss the changes in the gravitational lensing effects on damped Lyman alpha systems (DLAS) when disk lensing is included.

Many multiply--imaged quasars have been found over the years, but none so far with image separation in excess of $8\arcsec$. The absence of such large splittings has been used as a test of cosmological models: the standard Cold Dark Matter model has been excluded on the basis that it predicts far too many large--separation double images. These studies assume that the lensing structure has the mass profile of a singular isothermal sphere. However, such large splittings would be produced by very massive systems such as clusters of galaxies, for which other gravitational lensing data suggest less singular mass profiles. Here we analyze two cases of mass profiles for lenses: an isothermal sphere with a finite core radius (density $\rho \propto (r^2+r_{core}^2)^{-1})$, and a Hernquist profile ($\rho \propto r^{-1}(r+a)^{-3}$). We find that small core radii $r_{core} \sim 30 h^{-1}$ kpc, as suggested by the cluster data, or large $a \gsim 300 h^{-1}$ kpc, as needed for compatibility with gravitational distortion data, would reduce the number of large--angle splittings by an order of magnitude or more. Thus, it appears that these tests are sensitive both to the cosmological model (number density of lenses) and to the inner lens structure, which is unlikely to depend sensitively on the cosmology, making it difficult to test the cosmological models by large--separation quasar lensing until we reliably know the structure of the lenses themselves.

The distribution of dark matter around galactic or cluster halos has usually been assumed to be approximately isothermal with a non-zero core radius, which is expected to be of the order of the size of the visible matter distribution. Recently, the possibility has been raised that dark matter halos might be singular in the sense that the dark matter density $\rho$ could increase monotonically with radius $r$ down to a very small distance from the center of galaxies or clusters. Such central cusps in the dark matter density could lead to a high flux of gamma rays from WIMP dark matter annihilation. Here we analyze two possibilities that have been discussed in the literature, $\rho \propto r^{-n}$ with $n \approx 1\ {\rm or}\ 2$, and point out that such density profiles are excluded by gravitional lensing analyses on cluster scales and by the rotation curves of gas-rich, halo-dominated dwarf spirals on small scales. We also point out that if spiral galaxies form by gas infall inside dark matter halos, as they are expected to do in any hierarchical clustering model, such profiles almost always lead to falling rotation curves after infall, contrary to observations.