The development of numerical N -body simulations have allowed to study formation process and evolution of galaxies at different scales. This paper presents the fundamental concepts of N-body systems applied to the cosmological evolution of the $\Lambda$-Cold Dark Matter ($\Lambda$CDM) model. In order to perform structure formation in the Universe, we provide an introduction to the basic equations and their implementation on the GADGET-2 software. We also present a simple guide to modify this code. First, we briefly describe the dark matter in the Universe as well as the theoretical and experimental basis of the $\Lambda$CDM model. Then, we focus on the simulation codes and provide the equations that govern most of the N-body simulations to model the dark matter. We describe the Smoothed Particle Hydrodynamics method used for simulating the gas, star dynamics and structure formation in these simulations. Then, cautiously, we guide the reader to the installation of GADGET-2 on a Linux-based computer, as well as to carry out a couple of examples to operate the code. Finally, by using a computational cluster, we show several results of a large structure simulation, analyse the outputs to display the matter power spectrum, and compare the outcome with theoretical predictions.

In this paper, we present the classical Bondi accretion theory for the case of non-isothermal accretion processes onto a supermassive black hole (SMBH), including the effects of X-ray heating and the radiation force due to electron scattering and spectral lines. The radiation field is calculated by considering an optically thick, geometrically thin, standard accretion disk as the emitter of UV photons and a spherical central object as a source of X-ray emission. In the present analysis, the UV emission from the accretion disk is assumed to have an angular dependence, while the X-ray/central object radiation is assumed to be isotropic. This allows us to build streamlines in any angular direction we need to. The influence of both types of radiation is evaluated for different flux fractions of the X-ray and UV emissions with and without the effects of spectral line driving. We find that the radiation emitted near the SMBH interacts with the infalling matter and modifies the accretion dynamics. In the presence of line driving, a transition resembles from pure type 1 & 2 to type 5 solutions (see Fig2.1 of Frank etal. 2002), which takes place regardless of whether or not the UV emission dominates over the X-ray emission. We compute the radiative factors at which this transition occurs, and discard type 5 solution from all our models. Estimated values of the accretion radius and accretion rate in terms of the classical Bondi values are also given. The results are useful for the construction of proper initial conditions for time-dependent hydrodynamical simulations of accretion flows onto SMBH at the centre of galaxies.

In this paper we derive a novel circular velocity relation for a test particle in a 3D gravitational potential applicable to every system of curvilinear coordinates, suitable to be reduced to orthogonal form. As an illustration of the potentiality of the determined circular velocity expression we perform the rotation curves analysis of UGC 8490 and UGC 9753 and we estimate the total and dark matter mass of these two galaxies under the assumption that their respective dark matter halos have spherical, prolate and oblate spheroidal mass distributions. We employ stellar population synthesis models and the total HI density map to obtain the stellar and HI+He+metals rotation curves of both galaxies. The subtraction of the stellar plus gas rotation curves from the observed rotation curves of UGC 8490 and UGC 9753 generates the dark matter circular velocity curves of both galaxies. We fit the dark matter rotation curves of UGC 8490 and UGC 9753 through the newly established circular velocity formula specialised to the spherical, prolate and oblate spheroidal mass distributions, considering the Navarro, Frenk and White, Burkert, Di Cintio, Einasto and Stadel dark matter halos. Our principal findings are the following: globally, cored dark matter profiles Burkert and Einasto prevail over cuspy Navarro, Frenk and White and Di Cintio. Also, spherical/oblate dark matter models fit better the dark matter rotation curves of both galaxies than prolate dark matter halos.

We present three-dimensional calculations of spherically symmetric Bondi accretion onto a stationary supermassive black hole (SMBH) of mass $10^{8}$ $M_{\odot}$ within a radial range of $0.02-10$ pc, using a modified version of the smoothed particle hydrodynamics (SPH) \gad \sp code, which ensures approximate first-order consistency (i.e., second-order accuracy) for the particle approximation. First-order consistency is restored by allowing the number of neighbours, $n_{\rm neigh}$, and the smoothing length, $h$, to vary with the total number of particles, $N$, such that the asymptotic limits $n_{\rm neigh}\to\infty$ and $h\to 0$ hold as $N\to\infty$. The ability of the method to reproduce the isothermal ($\gamma =1$) and adiabatic ($\gamma =5/3$) Bondi accretion is investigated with increased spatial resolution. In particular, for the isothermal models the numerical radial profiles closely match the Bondi solution, except near the accretor, where the density and radial velocity are slightly underestimated. However, as $n_{\rm neigh}$ is increased and $h$ is decreased, the calculations approach first-order consistency and the deviations from the Bondi solution decrease. The density and radial velocity profiles for the adiabatic models are qualitatively similar to those for the isothermal Bondi accretion. Steady-state Bondi accretion is reproduced by the highly resolved consistent models with a percent relative error of $\lesssim 1$\% for $\gamma =1$ and $\sim 9$\% for $\gamma =5/3$, with the adiabatic accretion taking longer than the isothermal case to reach steady flow. The performance of the method is assessed by comparing the results with those obtained using the standard Gadget and the Gizmo codes.

Recent observations of young stellar systems with the Atacama Large Millimeter/submillimeter Array (ALMA) and the Karl G. Jansky Very Large Array (VLA) are helping to cement the idea that close companion stars form via fragmentation of a gravitationally unstable disk around a protostar early in the star formation process. As the disk grows in mass, it eventually becomes gravitationally unstable and fragments, forming one or more new protostars in orbit with the first at mean separations of 100 astronomical units (AU) or even less. Here we report direct numerical calculations down to scales as small as $\sim 0.1$ AU, using a consistent Smoothed Particle Hydrodynamics (SPH) code, that show the large-scale fragmentation of a cloud core into two protostars accompanied by small-scale fragmentation of their circumstellar disks. Our results demonstrate the two dominant mechanisms of star formation, where the disk forming around a protostar, which in turn results from the large-scale fragmentation of the cloud core, undergoes eccentric ($m=1$) fragmentation to produce a close binary. We generate two-dimensional emission maps and simulated ALMA 1.3 mm continuum images of the structure and fragmentation of the disks that can help explain the dynamical processes occurring within collapsing cloud cores.

We report the detection of the 1.3 mm continuum and the molecular emission of the disks of the young triple system SR24 by analyzing ALMA (The Atacama Large Millimeter/Submillimter Array) subarcsecond archival observations. We estimate the mass of the disks (0.025 solar masses and 4E-5 earth masses for SR24S and SR24N, respectively) and the dynamical mass of the protostars (1.5 and 1.1 solar masses). A kinematic model of the SR24S disk to fit its C18O (2-1) emission allows us to develop an observational method to learn what is the tilt of a rotating and accreting disk. We derive the size, the inclination, the position angle and the sense of rotation of each disk, finding that they are strongly misaligned (108 degrees) and possibly rotate in opposite directions as seen from Earth, in projection. We compare the ALMA observations with 12CO SMA archival observations, which are more sensitive to extended structures. We find three extended structures and estimate their masses: a molecular bridge joining the disks of the system, a molecular gas reservoir associated with SR24N and a gas streamer associated with SR24S. Finally we discuss on the possible origin of the misaligned SR24 system, pointing out that a closer inspection of the northern gas reservoir is needed to better understand it.

In this work we study the mass distribution of two irregular galaxies, UGC 6446 and UGC 7524, by means of HI rotation curves derived from high resolution HI velocity fields obtained through the Westerbork Synthesis Radio Telescope data archive. We constrain the stellar and gas content of both galaxies with stellar population synthesis models and by deriving the HI+He+metals rotation curves from the total HI surface density maps, respectively. The discrepancy between the circular velocity maxima of the stellar plus the HI+He+metals rotation curves and the observed HI rotation curves of both galaxies requires the inclusion of a substantial amount of dark matter. We explore the Navarro Frenk and White, Burkert, Di Cintio, Einasto and Stadel dark matter halo models. We obtain acceptable fits to the observed HI rotation curves of UGC 6446 and UGC 7524 with the cored Burkert, Einasto and Stadel dark matter halos. In particular, Einasto and Stadel models prove to be an appropriate alternative to the Burkert dark matter halo. This result should increase the empirical basis that justify the usage of dark matter exponential models to adjust the observed rotation curves of real galaxies.

We study the consistency and convergence of smoothed particle hydrodynamics (SPH), as a function of the interpolation parameters, namely the number of particles $N$, the number of neighbors $n$, and the smoothing length $h$, using simulations of the collapse and fragmentation of protostellar rotating cores. The calculations are made using a modified version of the GADGET-2 code that employs an improved scheme for the artificial viscosity and power-law dependences of $n$ and $h$ on $N$, as was recently proposed by Zhu et al., which comply with the combined limit $N\to\infty$, $h\to 0$, and $n\to\infty$ with $n/N\to 0$ for full SPH consistency, as the domain resolution is increased. We apply this realization to the "standard isothermal test case" in the variant calculated by Burkert & Bodenheimer and the Gaussian cloud model of Boss to investigate the response of the method to adaptive smoothing lengths in the presence of large density and pressure gradients. The degree of consistency is measured by tracking how well the estimates of the consistency integral relations reproduce their continuous counterparts. In particular, $C^{0}$ and $C^{1}$ particle consistency is demonstrated, meaning that the calculations are close to second-order accuracy. As long as $n$ is increased with $N$, mass resolution also improves as the minimum resolvable mass $M_{\rm min}\sim n^{-1}$. This aspect allows proper calculation of small-scale structures in the flow associated with the formation and instability of protostellar disks around the growing fragments, which are seen to develop a spiral structure and fragment into close binary/multiple systems as supported by recent observations.

We present digital tables for the radiative terms that appear in the energy and momentum equations used to simulate the accretion onto supermassive black holes (SMBHs) in the center of galaxies. Cooling and heating rates and radiative accelerations are calculated with two different Spectral Energy Distributions (SEDs). One SED is composed of an accretion disk + [X-ray]-powerlaw, while the other is made of an accretion disk + [Corona]-bremsstrahlung with $T_X=1.16 \times 10^8$ K, where precomputed conditions of adiabatic expansion are included. Quantification of different physical mechanisms at operation are presented, showing discrepancies and similarities between both SEDs in different ranges of fundamental physical parameters (i.e., ionization parameter, density, and temperature). With the recent discovery of outflows originating at sub-parsec scales, these tables may provide a useful tool to model gas accretion processes onto a SMBH.

In this article we investigate the outer and inner mass distributions of the irregular galaxies UGC 4284 and UGC 11861, taking advantage of published HI and H{\alpha} high resolution rotation curves and constraining the stellar disk of both galaxies throughout stellar population synthesis studies. In addition we take into account the gas content of both galaxies deriving the HI+He rotation curve. The deduced baryonic rotation curves (star+gas) are inadequate to account for the total mass of UGC 4284 and UGC 11861, for that reason we examine the possibility of dark matter to explain the incongruity between the observed HI and H{\alpha} rotation curves of UGC 4284 and UGC 11861 and the derived baryonic rotation curves. We consider NFW, Burkert, DiCintio, Einasto, and the Stadel dark matter halos, to analyse the dark matter content of UGC 4284 and UGC 11861. The principal results of this work are that cored dark matter models better reproduce the dark matter H{\alpha} and HI rotation curves of UGC 11861 and the dark matter HI rotation curve of UGC 4284, while, the H{\alpha} rotation curve of UGC 4284 is better reproduced by a cuspy DiCintio DM model. In general, cored exponential two-parameters models Einasto and Stadel, give better fits than Burkert. This trend, as well as to confirm past results, presents for the first time a comparison between two different exponential dark matter models, Einasto and Stadel, in an attempt to better constrain the range of possible exponential dark matter models applied to real galaxies.

Using the data obtained previously from Fabry-Perot interferometry, we study the orbital characteristics of the interacting pair of galaxies KPG 302 with the aim to estimate a possible interaction history, the conditions necessary for the spiral arms formation and initial satellite mass. We found by performing N-body/SPH simulations of the interaction that a single passage can produce a grand design spiral pattern in less than 1 Gyr. Althought we reproduce most of the features with the single passage, the required satellite to host mass ratio should be 1:5, which is not confirmed with the dynamical mass estimate made from the measured rotation curve. We conclude that a more realistic interaction scenario would require several passages in order to explain the mass ratio discrepancy.

The method of Smoothed Particle Hydrodynamics (SPH) has been widely studied and implemented for a large variety of problems, ranging from astrophysics to fluid dynamics and elasticity problems in solids. However, the method is known to have several deficiencies and discrepancies in comparison with traditional mesh-based codes. In particular, there has been a discussion about its ability to reproduce the Kelvin-Helmholtz Instability in shearing flows. Several authors reported that they were able to reproduce correctly the instability by introducing some improvements to the algorithm. In this contribution, we compare the results of Read et al. (2010) implementation of the SPH algorithm with the original Gadget-2 N-body/SPH code.

We present a kinematical study of the nearly edge-on galaxy ESO 379-G006 that shows the existence of extraplanar ionized gas. With Fabry-Perot spectroscopy at H-alpha, we study the kinematics of ESO 379-G006 using velocity maps and position-velocity diagrams parallel to the major and to the minor axis of the galaxy. We build the rotation curve of the disk and discuss the role of projection effects due to the fact of viewing this galaxy nearly edge-on. The twisting of the isovelocities in the radial velocity field of the disk of ESO 379-G006 as well as the kinematic asymmetries found in some position-velocity diagrams parallel to the minor axis of the galaxy suggest the existence of deviations to circular motions in the disk that can be modeled and explained with the inclusion of a radial inflow probably generated by a bar or by spiral arms. We succeeded in detecting extraplanar Diffuse Ionized Gas in this galaxy. At the same time, from the analysis of position-velocity diagrams, we found some evidence that the extraplanar gas could lag in rotation velocity with respect to the midplane rotation.

In this work we extend the study on the mass distribution of the spiral galaxy NGC 5278, performing 1D and 2D bulge-disk decomposition to determine which components constitute the baryonic mass in this galaxy. Our analysis does not detect any bulge, instead we find a bright source, probably related with the central AGN, and an exponential disk. We fix the stellar disk contribution to the rotation curve (RC) with broad band photometric observations and population synthesis models, to obtain 2D mass distribution of the stellar disk. In the particular case of NGC 5278, we find that the typical assumption of considering the mass-to-luminosity ratio (M/L) of the disk as constant along the galactocentric radius is not valid. We also extract a baryonic RC from the mass profile, to determine the inability of this baryonic RC and also of the baryonic RC with more and less 30% disk mass (in order to consider the disk mass errors) to fit the entire RC. We perform the RC decomposition of NGC 5278 considering the baryonic RC and four types of dark matter halos: Hernquist; Burkert; Einasto and Navarro, Frenk & White. Our results show that Hernquist halo models better our observed RC in the case of determined disk mass ($M_d=5.6\times 10^{10}$ solar masses) and also with 30% less disk mass. In the case of 30% more disk mass the cored Einasto (n < 4) halo is the best fitting model.

We present the photometric and kinematic characterization of two groups, USGC U268 and USGC U376 located in different regions of the Leo cloud. U268, composed of 10 catalogued members and 11 new added members, has a small fraction (~24%) of early-type galaxies (ETGs). U376 has 16 plus 8 new added members, with ~38% of ETGs. We find the presence of significant substructures in both groups suggesting that they are likely accreting galaxies. U268 is located in a more loose environment than U376. For each member galaxy, broad band integrated and surface photometry have been obtained in far-UV and near-UV with GALEX, and in u,g, r, i, z (SDSS) bands. H_alpha imaging and 2D high resolution kinematical data have been obtained using PUMA Scanning Fabry-Perot interferometer at the 2.12 m telescope in San Pedro M\'artir, (Baja California, M\'exico). We improved the galaxy classification and we detected morphological and kinematical distortions that may be connected to either on-going and/or past interaction/accretion events or environmental induced secular evolution. U268 appears more active than U376, with a large fraction of galaxies showing interaction signatures (60% vs. 13%). The presence of bars among late-type galaxies is ~10% in U268 and ~$29% in U376. The cumulative distribution of (FUV - NUV) colours of galaxies in U268 is significantly different than that in U376 with galaxies in U268 bluer than those in U376. In the (FUV-r vs. M_r) and (NUV-r vs. M_r) planes no members of U268 are found in the `red sequence', even early-type galaxies lie in the `blue sequence' or in the `green valley'. Most (80%) of the early-type members in U376 inhabits the `red sequence, a large fraction of galaxies, of different morphological types, are located in the `green valley', while the `blue sequence' is under-populated with respect to U268.

In this work we present scanning Fabry-Perot H$\alpha$ observations of the isolated interacting galaxy pair NGC 5278/79 obtained with the PUMA Fabry-Perot interferometer. We derived velocity fields and rotation curves for both galaxies. For NGC 5278 we also obtained the residual velocity map to investigate the non-circular motions, and estimated its mass by fitting the rotation curve with a disk+halo components. We test three different types of halo (pseudo-isothermal, Hernquist and Navarro Frenk White) and obtain satisfactory fits to the rotation curve for all profiles. The amount of dark matter required by pseudo-isothermal profile is about ten times smaller than, that for the other two halo distributions. Finally, our kinematical results together with the analysis of dust lanes distribution and of surface brightness profiles along the minor axis allowed us to determine univocally that both components of the interacting pair are trailing spirals.

In the context of the AMIGA project, we used Fabry-Perot observations in order to study the dynamics of the ionised gas in the isolated galaxy CIG 0314. From the Halpha observations, we could obtain the velocity field and rotation curve of the galaxy. A detail analysis of the velocity field is done in order to understand the kinematics of the gas to gather clues on the mechanisms which favour or inhibit star formation, in particular along the bar. The visible and dark matter content can be reached, as well as an estimation of the mass of the galaxy.

An important dynamic parameter of barred galaxies is the bar pattern speed. Among several methods that are used for the determination of the pattern speed the Tremaine-Weinberg method has the advantage of model independency and accuracy. In this work we apply the method to a simulated bar including gas dynamics and study the effect of 2D spectroscopy data quality on robustness of the method. We added a white noise and a Gaussian random field to the data and measured the corresponding errors in the pattern speed. We found that a signal to noise ratio in surface density ~5 introduces errors of ~20% for the Gaussian noise, while for the white noise the corresponding errors reach ~50%. At the same time the velocity field is less sensitive to contamination. On the basis of the performed study we applied the method to the NGC 3367 spiral galaxy using H{\alpha} Fabry-Perot interferometry data. We found for the pattern speed 43 \pm 6 km/s/kpc for this galaxy.

We study bar formation in galactic disks as a consequence of the collision of two spiral galaxies under the influence of a potential which is obtained from the Newtonian limit of a scalar--tensor theory of gravity. We found that dynamical effects depend on parameters ($\alpha$, $\lambda$) of the theory. In particular, we observe that the bar is shorter for weaker tidal perturbations, which in turn corresponds to smaller values of $\lambda$ used in our numerical experiments.

We present $\Lambda$CDM $N$-body cosmological simulations in the framework of a static general scalar-tensor theory of gravity. Due to the influence of the non-minimally coupled scalar field, the gravitational potential is modified by a Yukawa type term, yielding a new structure formation dynamics. We present some preliminary results and, in particular, we compute the density and velocity profiles of the most massive group.

We use the Newtonian limit of a general scalar-tensor theory around a background field to study astrophysical effects. The gravitational theory modifies the standard Newtonian potential by adding a Yukawa term to it, which is quantified by two theoretical parameters: $\lambda$, the lenghtscale of the gravitational interaction and its strength, $\alpha$. Within this formalism we firstly present a numerical study on the formation of bars in isolated galaxies. We have found for positive $\alpha$ that the modified gravity destabilizes the galactic discs and leads to rapid bar formation in isolated galaxies. Values of $\lambda$ in the range $\approx 8$ -- 14 kpc produce strongest bars in isolated models. Then, we extent this work by considering tidal effects due to interacting galaxies. We send two spirals to collide and study the bar properties of the remnant. We characterize the bar kinematical properties in terms of our parameters ($\lambda, \alpha$).

The joint influence of numerical parameters such as the number of particles N, the gravitational softening length $\epsilon$ and the time-step $\Delta t$ is investigated in the context of galaxy simulations. For isolated galaxy models we have performed a convergence study and estimated the numerical parameters ranges for which the relaxed models do not deviate significantly from its initial configuration. By fixing N, we calculate the range of the mean interparticle separation $\lambda(r)$ along the disc radius. We have found that in the simulations with N=1310720 particles $\lambda$ varies by a factor of 6, and the corresponding final Toomre's parameters Q change by only about 5 per cent. By decreasing N, the $\lambda$ and Q ranges broaden. Large $\epsilon$ and small N cause an earlier bar formation. For a given set of parameters the disc heating is smaller with the Plummer softening than with the spline softening. For galaxy collision models numerical simulations indicate that the properties of the formed bars strongly depend upon the selection of N and $\epsilon$. Large values of the gravitational softening parameter and a small number of particles results in the rapid formation of a well defined, slowly rotating bar. On the other hand, small values of $\epsilon$ produce a small, rapidly rotating disc with tightly wound spiral arms, and subsequently a weak bar emerges. We have found that by increasing N, the bar properties converge and the effect of the softening parameter diminishes. Finally, in some cases short spiral arms are observed at the ends of the bar that change periodically from trailing to leading and vice-versa - the wiggle.