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 the problem of consistency of smoothed particle hydrodynamics (SPH) is solved. A novel error analysis is developed in $n$-dimensional space using the Poisson summation formula, which enables the treatment of the kernel and particle approximation errors in combined fashion. New consistency integral relations are derived for the particle approximation which correspond to the cosine Fourier transform of the classically known consistency conditions for the kernel approximation. The functional dependence of the error bounds on the SPH interpolation parameters, namely the smoothing length $h$ and the number of particles within the kernel support ${\cal{N}}$ is demonstrated explicitly from which consistency conditions are seen to follow naturally. As ${\cal{N}}\to\infty$, the particle approximation converges to the kernel approximation independently of $h$ provided that the particle mass scales with $h$ as $m\propto h^{\beta}$, with $\beta >n$. This implies that as $h\to 0$, the joint limit $m\to 0$, ${\cal{N}}\to\infty$, and $N\to\infty$ is necessary for complete convergence to the continuum, where $N$ is the total number of particles. The analysis also reveals the presence of a dominant error term of the form $(\ln {\cal{N}})^{n}/{\cal{N}}$, which tends asymptotically to $1/{\cal{N}}$ when ${\cal{N}}\gg 1$, as it has long been conjectured based on the similarity between the SPH and the quasi-Monte Carlo estimates.

In this paper we describe the FIT\textit{spec} code, a data mining tool for the automatic fitting of synthetic stellar spectra. The program uses a database of 27\,000 {\sc cmfgen} models of stellar atmospheres arranged in a six-dimensional (6D) space, where each dimension corresponds to one model parameter. From these models a library of 2\,835\,000 synthetic spectra were generated covering the ultraviolet, optical, and infrared region of the electromagnetic spectrum. Using FIT\textit{spec} we adjust the effective temperature and the surface gravity. From the 6D array we also get the luminosity, the metallicity, and three parameters for the stellar wind: the terminal velocity ($v_\infty$), the $\beta$ exponent of the velocity law, and the clumping filling factor ($F_{\rm cl}$). Finally, the projected rotational velocity ($v\cdot\sin i$) can be obtained from the library of stellar spectra. Validation of the algorithm was performed by analyzing the spectra of a sample of eight O-type stars taken from the {\sc iacob} spectroscopic survey of Northern Galactic OB stars. The spectral lines used for the adjustment of the analyzed stars are reproduced with good accuracy. In particular, the effective temperatures calculated with the FIT\textit{spec} are in good agreement with those derived from spectral type and other calibrations for the same stars. The stellar luminosities and projected rotational velocities are also in good agreement with previous quantitative spectroscopic analyses in the literature. An important advantage of FIT\textit{spec} over traditional codes is that the time required for spectral analyses is reduced from months to a few hours.

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.

With the use of simulated supernova catalogs, we show that the statefinder parameters turn out to be poorly and biased estimated by standard cosmography. To this end, we compute their standard deviations and several bias statistics on cosmologies near the concordance model, demonstrating that these are very large, making standard cosmography unsuitable for future and wider compilations of data. To overcome this issue, we propose a new method that consists in introducing the series of the Hubble function into the luminosity distance, instead of considering the usual direct Taylor expansions of the luminosity distance. Moreover, in order to speed up the numerical computations, we estimate the coefficients of our expansions in a hierarchical manner, in which the order of the expansion depends on the redshift of every single piece of data. In addition, we propose two hybrids methods that incorporates standard cosmography at low redshifts. The methods presented here perform better than the standard approach of cosmography both in the errors and bias of the estimated statefinders. We further propose a one-parameter diagnostic to reject non-viable methods in cosmography.

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.

The problem of consistency of smoothed particle hydrodynamics (SPH) has demanded considerable attention in the past few years due to the ever increasing number of applications of the method in many areas of science and engineering. A loss of consistency leads to an inevitable loss of approximation accuracy. In this paper, we revisit the issue of SPH kernel and particle consistency and demonstrate that SPH has a limiting second-order convergence rate. Numerical experiments with suitably chosen test functions validate this conclusion. In particular, we find that when using the root mean square error as a model evaluation statistics, well-known corrective SPH schemes, which were thought to converge to second, or even higher order, are actually first-order accurate, or at best close to second order. We also find that observing the joint limit when $N\to\infty$, $h\to 0$, and $n\to\infty$, as was recently proposed by Zhu et al., where $N$ is the total number of particles, $h$ is the smoothing length, and $n$ is the number of neighbor particles, standard SPH restores full $C^{0}$ particle consistency for both the estimates of the function and its derivatives and becomes insensitive to particle disorder.

We test here the first stage of a route of modifications to be applied to the public GADGET2 code for dynamically identifying accretion centers during the collision process of two adjacent and identical gas cores. Each colliding core has a uniform density profile and rigid body rotation; its mass and size have been chosen to represent the observed core $L1544$; for the thermal and rotational energy ratios with respect to the potential energy, we assume the values $\alpha=0.3$ and $\beta=0.1$, respectively. These values favor the gravitational collapse of the core. We here study cases of both -head-on and off-center collisions, in which the pre-collision velocity increases the initial sound speed of the barotropic gas by up to several times. In a simulation the accretion centers are formed by the highest density particles, so we here report their location and properties in order to realize the collision effects on the collapsing and colliding cores. In one of the models we observe a roughly spherical distribution of accretion centers located at the front wave of the collision. In a forthcoming publication we will apply the full modified GADGET code to study the collision of turbulent cores.

We investigate the thermodynamics of general fluids that have the constriction that their pressure is constant. For example, this happens in the case of pure dust matter, for which the pressure vanishes and also in the case of standard dark matter phenomenology. Assuming a finite non-zero pressure, the corresponding dynamics is richer than one naively would expect. In particular, it can be considered as a unified description of dark energy and dark matter. We first consider the more general thermodynamic properties of this class of fluids finding the important result that for them adiabatic and isothermal processes should coincide. We therefore study their behaviors in curved space-times where local thermal equilibrium can be appealed. Thus, we show that this dark fluid degenerates with the dark sector of the LCDM model only in the case of adiabatic evolution. We demonstrate that, adding dissipative processes, a phantom behavior can occur and finally we further highlight that an arbitrary decomposition of the dark sector, into ad hoc dark matter and dark energy terms, may give rise to phantom dark energy, whereas the whole dark sector remains non-phantom.

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.

We review the role of fluids in cosmology by first introducing them in General Relativity and then by applying them to a FRW Universe's model. We describe how relativistic and non-relativistic components evolve in the background dynamics. We also introduce scalar fields to show that they are able to yield an inflationary dynamics at very early times (inflation) and late times (quintessence). Then, we proceed to study the thermodynamical properties of the fluids and, lastly, its perturbed kinematics. We make emphasis in the constrictions of parameters by recent cosmological probes.

We examine the problem of the collapse and fragmentation of molecular clouds with a Gaussian density distribution with high resolution, double precision numerical simulations using the GADGET-2 code. To describe the thermodynamic properties of the cloud during the collapse -to mimic the rise of temperature predicted by radiative transfer- we use a barotropic equation of state that introduces a critical density to separate the isothermal and adiabatic regimes. We discuss the effects of this critical density in the formation of multiple systems. We confirm the tendency found for Plummer and Gaussian models that if the collapse changes from isothermal to adiabatic at earlier times that occurs for the models with a lower critical density, the collapse is slowed down, and this enhances the fragments' change to survive. However, this effect happens up to a threshold density below which single systems tend to form. On the other hand, by setting a bigger initial perturbation amplitude, the collapse is faster and in some cases a final single object is formed.

The first generation of stars was formed from primordial gas. Numerical simulations suggest that the first stars were predominantly very massive, with typical masses M > 100 Mo. These stars were responsible for the reionization of the universe, the initial enrichment of the intergalactic medium with heavy elements, and other cosmological consequences. In this work, we study the structure of Zero Age Main Sequence stars for a wide mass and metallicity range and the evolution of 100, 150, 200, 250 and 300 Mo galactic and pregalactic Pop III very massive stars without mass loss, with metallicity Z=10E-6 and 10E-9, respectively. Using a stellar evolution code, a system of 10 equations together with boundary conditions are solved simultaneously. For the change of chemical composition, which determines the evolution of a star, a diffusion treatment for convection and semiconvection is used. A set of 30 nuclear reactions are solved simultaneously with the stellar structure and evolution equations. Several results on the main sequence, and during the hydrogen and helium burning phases, are described. Low metallicity massive stars are hotter and more compact and luminous than their metal enriched counterparts. Due to their high temperatures, pregalactic stars activate sooner the triple alpha reaction self-producing their own heavy elements. Both galactic and pregalactic stars are radiation pressure dominated and evolve below the Eddington luminosity limit with short lifetimes. The physical characteristics of the first stars have an important influence in predictions of the ionizing photon yields from the first luminous objects; also they develop large convective cores with important helium core masses which are important for explosion calculations.

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 simulations within the framework of scalar-tensor theories, in the Newtonian limit, to investigate the influence of massive scalar fields on the dynamics of the collision of two equal spherical clouds. We employ a SPH code modified to include the scalar field to simulate two initially non-rotating protogalaxies that approach each other, and as a result of the tidal interaction, intrinsic angular momentum is generated. We have obtained sufficient large values of J/M to suggest that intrinsic angular momentum can be the result of tidal interactions.

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.