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.

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.

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.