Starting from the single graphics processing unit (GPU) version of the Smoothed Particle Hydrodynamics (SPH) code DualSPHysics, a multi-GPU SPH program is developed for free-surface flows. The approach is based on a spatial decomposition technique, whereby different portions (sub-domains) of the physical system under study are assigned to different GPUs. Communication between devices is achieved with the use of Message Passing Interface (MPI) application programming interface (API) routines. The use of the sorting algorithm radix sort for inter-GPU particle migration and sub-domain halo building (which enables interaction between SPH particles of different subdomains) is described in detail. With the resulting scheme it is possible, on the one hand, to carry out simulations that could also be performed on a single GPU, but they can now be performed even faster than on one of these devices alone. On the other hand, accelerated simulations can be performed with up to 32 million particles on the current architecture, which is beyond the limitations of a single GPU due to memory constraints. A study of weak and strong scaling behaviour, speedups and efficiency of the resulting programis presented including an investigation to elucidate the computational bottlenecks. Last, possibilities for reduction of the effects of overhead on computational efficiency in future versions of our scheme are discussed.

Experiments showing a seasonal variation of the nuclear decay rates of a number of different nuclei, and decay anomalies apparently related to solar flares and solar rotation, have suggested that the Sun may somehow be influencing nuclear decay processes. Recently, Cooper searched for such an effect in $^{238}$Pu nuclei contained in the radioisotope thermoelectric generators (RTGs) on board the Cassini spacecraft. In this paper we modify and extend Cooper's analysis to obtain constraints on anomalous decays of $^{238}$Pu over a wider range of models, but these limits cannot be applied to other nuclei if the anomaly is composition-dependent. We also show that it may require very high sensitivity for terrestrial experiments to discriminate among some models if such a decay anomaly exists, motivating the consideration of future spacecraft experiments which would require less precision.

Differences in the stellar populations of galaxies can be used to quantify the effect of environment on the star formation history. We target a sample of early-type galaxies from the Sloan Digital Sky Survey in two different environmental regimes: close pairs and a general sample where environment is measured by the mass of their host dark matter halo. We apply a blind source separation technique based on principal component analysis, from which we define two parameters that correlate, respectively, with the average stellar age (eta) and with the presence of recent star formation (zeta) from the spectral energy distribution of the galaxy. We find that environment leaves a second order imprint on the spectra, whereas local properties - such as internal velocity dispersion - obey a much stronger correlation with the stellar age distribution.

(Abridged) We present a detailed study of the stellar populations of a volume-limited sample of early-type galaxies from SDSS, across a range of environments -- defined as the mass of the host dark matter halo. The stellar populations are explored through the SDSS spectra, via projection onto a set of two spectral vectors determined from Principal Component Analysis. We find the velocity dispersion of the galaxy to be the main driver behind the different star formation histories of early-type galaxies. However, environmental effects are seen to play a role (although minor). Galaxies populating the lowest mass halos have stellar populations on average ~1Gyr younger than the rest of the sample. The fraction of galaxies with small amounts of recent star formation is also seen to be truncated when occupying halos more massive than 3E13Msun. The sample is split into satellite and central galaxies for a further analysis of environment. Satellites are younger than central galaxies of the same stellar mass. The younger satellite galaxies in 6E12Msun halos have stellar populations consistent with the central galaxies found in the lowest mass halos of our sample (i.e. 1E12Msun). This result is indicative of galaxies in lower mass halos being accreted into larger halos.

We consider ion heating by turbulent Alfven waves (AWs) and kinetic Alfven waves (KAWs) with perpendicular wavelengths comparable to the ion gyroradius and frequencies smaller than the ion cyclotron frequency. When the turbulence amplitude exceeds a certain threshold, an ion's orbit becomes chaotic. The ion then interacts stochastically with the time-varying electrostatic potential, and the ion's energy undergoes a random walk. Using phenomenological arguments, we derive an analytic expression for the rates at which different ion species are heated, which we test by simulating test particles interacting with a spectrum of randomly phased AWs and KAWs. We find that the stochastic heating rate depends sensitively on the quantity epsilon = dv/vperp, where vperp is the component of the ion velocity perpendicular to the background magnetic field B0, and dv (dB) is the rms amplitude of the velocity (magnetic-field) fluctuations at the gyroradius scale. In the case of thermal protons, when epsilon << eps1, where eps1 is a constant, a proton's magnetic moment is nearly conserved and stochastic heating is extremely weak. However, when epsilon > eps1, the proton heating rate exceeds the cascade power that would be present in strong balanced KAW turbulence with the same value of dv, and magnetic-moment conservation is violated. For the random-phase waves in our test-particle simulations, eps1 is approximately 0.2. For protons in low-beta plasmas, epsilon is approximately dB/B0 divided by the square root of beta, and epsilon can exceed eps1 even when dB/B0 << eps1. At comparable temperatures, alpha particles and minor ions have larger values of epsilon than protons and are heated more efficiently as a result. We discuss the implications of our results for ion heating in coronal holes and the solar wind.

(Abridged) We study the stellar populations of 14 elliptical galaxies in the Virgo cluster. We propose an alternative approach to the standard side-band method to measure equivalent widths (EWs). Our Boosted Median Continuum maps the EWs more robustly than the side-band method, minimising the effect from neighbouring absorption lines and reducing the age-metallicity degeneracy. We concentrate on Balmer lines (Hbeta,Hgamma,Hdelta), the G band and the 4000A break as age-sensitive indicators, and on the combination [MgFe] as the main metallicity indicator. We go beyond the standard comparison of the observations with simple stellar populations (SSP) and consider various models to describe the star formation histories, either with a continuous star formation rate or with a mixture of two different SSPs. Composite models are found to give more consistent fits among individual line strengths and agree with an independent estimate using the spectral energy distribution. Our age and metallicity estimates correlate well with stellar mass or velocity dispersion, with a significant threshold around 5E10 Msun above which galaxies are uniformly old and metal rich. In a more speculative way, our models suggest that it is formation **epoch** and not formation timescale what drives the Mass-Age relationship of elliptical galaxies.

We extract from the Sloan Digital Sky Survey a sample of 347 systems involving early type galaxies separated by less than 30 kpc, in projection, and 500 km/s in radial velocity. These close pairs are likely progenitors of dry mergers. The (optical) spectra is used to determine how the interaction affects the star formation history and nuclear activity of the galaxies. The emission lines (or lack thereof) are used to classify the sample into AGN, star forming or quiescent. Increased AGN activity and reduced star formation in early-type pairs that already appear to be interacting indicate that the merging process changes the nature of nebular activity, a finding that is also supported by an increase in AGN luminosity with decreasing pair separation. Recent star formation is studied on the absorption line spectra, both through principal component analysis as well as via a comparison of the spectra with composite stellar population models. We find that the level of recent star formation in close pairs is raised relative to a control sample of early-type galaxies. This excess of residual star formation is found throughout the sample of close pairs and does not correlate with pair separation or with visual signs of interaction. Our findings are consistent with a scenario whereby the first stage of the encounter (involving the outer parts of the halos) trigger residual star formation, followed by a more efficient inflow towards the centre -- switching to an AGN phase -- after which the systems are quiescent.

(Abridged) We apply Principal Component Analysis (PCA) to a sample of early-type galaxies from the Sloan Digital Sky Survey (SDSS) in order to infer differences in their star formation histories. Out of the first few principal components (PC), we study four which give information about stellar populations and velocity dispersion. We construct two parameters (eta and zeta) as linear combinations of PC1 and PC2. We find zeta to be most sensitive to recent episodes of star formation, and eta to be strongly dependent on the average age of the stellar populations. The distribution of the eta component of the composites appear to be indistinguishable between high and low density regions, whereas the distribution of zeta parameters have a significant skew towards lower values for galaxies in low density regions. This result suggests that galaxies in lower density environments are less likely to present weak episodes of recent star formation. In contrast, a significant number of galaxies from our high density subsample -- which includes clusters (both outer regions and centres) and groups -- underwent small but detectable recent star formation at high metallicity, in agreement with recent estimates targeting elliptical galaxies in Hickson Compact Groups and in the field (Ferreras et al.).

Principal Component Analysis (PCA) is a well-known multivariate technique used to decorrelate a set of vectors. PCA has been extensively applied in the past to the classification of stellar and galaxy spectra. Here we apply PCA to the optical spectra of early-type galaxies, with the aim of extracting information about their star formation history. We consider two different data sets: 1) a reduced sample of 30 elliptical galaxies in Hickson compact groups and in the field, and 2) a large volume-limited (z<0.1) sample of ~7,000 galaxies from the Sloan Digital Sky Survey. Even though these data sets are very different, the homogeneity of the populations results in a very similar set of principal components. Furthermore, most of the information (in the sense of variance) is stored into the first few components in both samples. The first component (PC1) can be interpreted as an old population and carries over 99% of the variance. The second component (PC2) is related to young stars and we find a correlation with NUV flux from GALEX. Model fits consistently give younger ages for those galaxies with higher values of PC2.

We present particle-in-cell simulations of collisionless magnetic reconnection in a system (like the magnetopause) with a large density asymmetry across the current layer. In the presence of an ambient component of the magnetic field perpendicular to the reconnection plane the gradient creates a diamagnetic drift that advects the X-line with the electron diamagnetic velocity. When the relative drift between the ions and electrons is of the order the Alfven speed the large scale outflows from the X-line necessary for fast reconnection cannot develop and the reconnection is suppressed. We discuss how these effects vary with both the plasma beta and the shear angle of the reconnecting field and discuss observational evidence for diamagnetic stabilization at the magnetopause.