Satellites of Simulated Galaxies: survival, merging, and their relation to the dark and stellar halos

We study the population of satellite galaxies formed in a suite of N-body/gasdynamical simulations of galaxy formation in a LCDM universe. We find little spatial or kinematic bias between the dark matter and the satellite population. The velocity dispersion of the satellites is a good indicator of the virial velocity of the halo: \sigma_{sat}/V_{vir}=0.9 +/- 0.2. Applied to the Milky Way and M31 this gives V_{vir}^{MW}=109 +/- 22$ km/s and V_{vir}^{M31} = 138 +/- 35 km/s, respectively, substantially lower than the rotation speed of their disk components. The detailed kinematics of simulated satellites and dark matter are also in good agreement. By contrast, the stellar halo of the simulated galaxies is kinematically and spatially distinct from the population of surviving satellites. This is because the survival of a satellite depends on mass and on time of accretion; surviving satellites are biased toward low-mass systems that have been recently accreted by the galaxy. Our results support recent proposals for the origin of the systematic differences between stars in the Galactic halo and in Galactic satellites: the elusive ``building blocks'' of the Milky Way stellar halo were on average more massive, and were accreted (and disrupted) earlier than the population of dwarfs that has survived self-bound until the present.Comment: 13 pages, 11 figures, MNRAS in press. Accepted version with minor changes. Version with high resolution figures available at: http://www.astro.uvic.ca/~lsales/SatPapers/SatPapers.htm

Cosmic M\'enage \`a Trois: The Origin of Satellite Galaxies On Extreme Orbits

We examine the orbits of satellite galaxies identified in a suite of N-body/gasdynamical simulations of the formation of $L_*$ galaxies in a LCDM universe. Most satellites follow conventional orbits; after turning around, they accrete into their host halo and settle on orbits whose apocentric radii are steadily eroded by dynamical friction. However, a number of outliers are also present, we find that ~1/3 of satellites identified at $z=0$ are on unorthodox orbits, with apocenters that exceed their turnaround radii. This population of satellites on extreme orbits consists typically of the faint member of a satellite pair that has been ejected onto a highly-energetic orbit during its first approach to the primary. Since the concurrent accretion of multiple satellite systems is a defining feature of hierarchical models of galaxy formation, we speculate that this three-body ejection mechanism may be the origin of (i) some of the newly discovered high-speed satellites around M31 (such as Andromeda XIV); (ii) some of the distant fast-receding Local Group members, such as Leo I; and (iii) the oddly isolated dwarf spheroidals Cetus and Tucana in the outskirts of the Local Group. Our results suggest that care must be exercised when using the orbits of the most weakly bound satellites to place constraints on the total mass of the Local Group.Comment: 10 pages, 6 figures, MNRAS in press. Accepted version with minor changes. Version with high resolution figures available at: http://www.astro.uvic.ca/~lsales/SatPapers/SatPapers.htm

Simulations of galaxy formation in a Λ cold dark matter universe : I : dynamical and photometric properties of a simulated disk galaxy.

We present a detailed analysis of the dynamical and photometric properties of a disk galaxy simulated in the cold dark matter (CDM) cosmogony. The galaxy is assembled through a number of high-redshift mergers followed by a period of quiescent accretion after z1 that lead to the formation of two distinct dynamical components: a spheroid of mostly old stars and a rotationally supported disk of younger stars. The surface brightness profile is very well approximated by the superposition of an R1/4 spheroid and an exponential disk. Each photometric component contributes a similar fraction of the total luminosity of the system, although less than a quarter of the stars form after the last merger episode at z1. In the optical bands the surface brightness profile is remarkably similar to that of Sab galaxy UGC 615, but the simulated galaxy rotates significantly faster and has a declining rotation curve dominated by the spheroid near the center. The decline in circular velocity is at odds with observation and results from the high concentration of the dark matter and baryonic components, as well as from the relatively high mass-to-light ratio of the stars in the simulation. The simulated galaxy lies 1 mag off the I-band Tully-Fisher relation of late-type spirals but seems to be in reasonable agreement with Tully-Fisher data on S0 galaxies. In agreement with previous simulation work, the angular momentum of the luminous component is an order of magnitude lower than that of late-type spirals of similar rotation speed. This again reflects the dominance of the slowly rotating, dense spheroidal component, to which most discrepancies with observation may be traced. On its own, the disk component has properties rather similar to those of late-type spirals: its luminosity, its exponential scale length, and its colors are all comparable to those of galaxy disks of similar rotation speed. This suggests that a different form of feedback than adopted here is required to inhibit the efficient collapse and cooling of gas at high redshift that leads to the formation of the spheroid. Reconciling, without fine-tuning, the properties of disk galaxies with the early collapse and high merging rates characteristic of hierarchical scenarios such as CDM remains a challenging, yet so far elusive, proposition

Non-Abelian Chern-Simons-Higgs vortices with a quartic potential

We have constructed numerically non-Abelian vortices in an SU(2) Chern-Simons-Higgs theory with a quartic Higgs potential. We have analyzed these solutions in detail by means of improved numerical codes and found some unexpected features we did not find when a sixth-order Higgs potential was used. The generic non-Abelian solutions have been generated by using their corresponding Abelian counterparts as initial guess. Typically, the energy of the non-Abelian solutions is lower than that of the corresponding Abelian one (except in certain regions of the parameter space). Regarding the angular momentum, the Abelian solutions possess the maximal value, although there exist non-Abelian solutions which reach that maximal value too. In order to classify the solutions it is useful to consider the non-Abelian solutions with asymptotically vanishing $A_t$ component of the gauge potential, which may be labelled by an integer number $m$. For vortex number $n=3$ and above, we have found uniqueness violation: two different non-Abelian solutions with all the global charges equal. Finally, we have investigated the limit of infinity Higgs self-coupling parameter and found a piecewise Regge-like relation between the energy and the angular momentum.Comment: 9 pages, 13 figure

Locoregional hyperthermia of deep-seated tumours applied with capacitive and radiative systems. A simulation study

Background: Locoregional hyperthermia is applied to deep-seated tumours in the pelvic region. Two very different heating techniques are often applied: capacitive and radiative heating. In this paper, numerical simulations are applied to compare the performance of both techniques in heating of deep-seated tumours. Methods: Phantom simulations were performed for small (30 × 20 × 50 cm 3 ) and large (45 × 30 × 50 cm 3 ), homogeneous fatless and inhomogeneous fat-muscle, tissue-equivalent phantoms with a central or eccentric target region. Radiative heating was simulated with the 70 MHz AMC-4 system and capacitive heating was simulated at 13.56 MHz. Simulations were performed for small fatless, small (i.e. fat layer typically 3 cm) patients with cervix, prostate, bladder and rectum cancer. Temperature distributions were simulated using constant hyperthermic-level perfusion values with tissue constraints of 44 °C and compared for both heating techniques. Results: For the small homogeneous phantom, similar target heating was predicted with radiative and capacitive heating. For the large homogeneous phantom, most effective target heating was predicted with capacitive heating. For inhomogeneous phantoms, hot spots in the fat layer limit adequate capacitive heating, and simulated target temperatures with radiative heating were 2–4 °C higher. Patient simulations predicted therapeutic target temperatures with capacitive heating for fatless patients, but radiative heating was more robust for all tumour sites and patient sizes, yielding target temperatures 1–3 °C higher than those predicted for capacitive heating. Conclusion: Generally, radiative locoregional heating yields more favourable simulated temperature distributions for deep-seated pelvic tumours, compared with capacitive heating. Therapeutic temperatures are predicted for capacitive heating in patients with (almost) no fat

Isospin and density dependences of nuclear matter symmetry energy coefficients II

Symmetry energy coefficients of explicitly isospin asymmetric nuclear matter at variable densities (from .5$\rho_0$ up to 2 $\rho_0$) are studied as generalized screening functions. An extended stability condition for asymmetric nuclear matter is proposed. We find the possibility of obtaining stable asymmetric nuclear matter even in some cases for which the symmetric nuclear matter limit is unstable. Skyrme-type forces are extensively used in analytical expressions of the symmetry energy coefficients derived as generalized screening functions in the four channels of the particle hole interaction producing alternative behaviors at different $\rho$ and $b$ (respectively the density and the asymmetry coefficient). The spin and spin-isospin coefficients, with corrections to the usual Landau Migdal parameters, indicate the possibility of occurring instabilities with common features depending on the nuclear density and n-p asymmetry. Possible relevance for high energy heavy ions collisions and astrophysical objects is discussed.Comment: 16 pages (latex) plus twelve figures in four eps files, to be published in I.J.M.P.

Asymptotic iteration method for eigenvalue problems

An asymptotic interation method for solving second-order homogeneous linear differential equations of the form y'' = lambda(x) y' + s(x) y is introduced, where lambda(x) \neq 0 and s(x) are C-infinity functions. Applications to Schroedinger type problems, including some with highly singular potentials, are presented.Comment: 14 page

Double percolation effects and fractal behavior in magnetic/superconducting hybrids

Perpendicular magnetic anisotropy ferromagnetic/ superconducting (FM/SC) bilayers with a labyrinth domain structure are used to study nucleation of superconductivity on a fractal network, tunable through magnetic history. As clusters of reversed domains appear in the FM layer, the SC film shows a percolative behavior that depends on two independent processes: the arrangement of initial reversed domains and the fractal geometry of expanding clusters. For a full labyrinth structure, the behavior of the upper critical field is typical of confined superconductivity on a fractal network.Comment: 15 pages, 5 figure