A density-dependent two-nucleon potential has been derived in the formalism of correlated basis function. The effects of 3-particle interactions has been included by integrating out the degrees of freedom of the third nucleon. The potential can be easily employed in nuclear matter calculations. It yields results in agreement with those obtained from the underlying three-body potential. The use of the density dependent potential allowed us to study the effects of three-nucleon interactions in symmetric nuclear matter within the Auxiliary Field Diffusion Monte Carlo (AFDMC) computational scheme.

The problem of accurately determining the equation of state of nuclear and neutron matter at density near and beyond saturation is still an open challenge. In this paper we will review the most recent progress made by means of Quantum Monte Carlo calculations, which are at present the only ab-inito method capable to treat a sufficiently large number of particles to give meaningful estimates depending only on the choice of the nucleon-nucleon interaction. In particular, we will discuss the introduction of density-dependent interactions, the study of the temperature dependence of the equation of state, and the possibility of accurately studying the effect of the onset of hyperons by developing an accurate hyperon-nucleon and hyperon-nucleon-nucleon interaction.

We present results for neutron star models constructed with a new equation of state for nuclear matter at zero temperature. The ground state is computed using the Auxiliary Field Diffusion Monte Carlo (AFDMC) technique, with nucleons interacting via a semi-phenomenological Hamiltonian including a realistic two-body interaction. The effect of many-body forces is included by means of additional density-dependent terms in the Hamiltonian. In this letter we compare the properties of the resulting neutron-star models with those obtained using other nuclear Hamiltonians, focusing on the relations between mass and radius, and between the gravitational mass and the baryon number.

We calculate the equation of state of neutron matter at zero temperature by means of the auxiliary field diffusion Monte Carlo method (AFDMC) combined with a fixed-phase approximation. The calculation of the energy is carried out by simulating up to 114 neutrons in a periodic box. Special attention was made to reduce finite size effects at the energy evaluation by adding to the interaction the effect due to the truncation of the simulation box, and by performing several simulations using different number of neutrons. The finite size effects due to the kinetic energy were also checked by employing the twist--averaged boundary conditions. We considered a realistic nuclear Hamiltonian containing modern two-- and three--body interactions of the Argonne and Urbana family. The equation of state can be used to compare and to calibrate other many-body calculations and to predict properties of neutron stars.

We present a Quantum Monte Carlo study of the zero temperature equation of state of neutron matter and the computation of the $^1S_0$ pairing gap in the low-density regime with $\rho<0.04$ fm$^{-3}$. The system is described by a non-relativistic nuclear Hamiltonian including both two-- and three--nucleon interactions of the Argonne and Urbana type. This model interaction provides very accurate results in the calculation of the binding energy of light nuclei. A suppression of the gap with respect to the pure BCS theory is found, but sensibly weaker than in other works that attempt to include polarization effects in an approximate way.

The $^1S_0$ pairing in neutron matter has been investigated in presence of realistic two-- and three--nucleon interactions. We have adopted the Argonne $v_{8^\prime}$ NN and the Urbana IX 3N potentials. Quantum Monte Carlo theory, specifically the Auxiliary Field Diffusion Monte Carlo method, and Correlated Basis Function theory are employed in order to get quantitative and reliable estimates of the gap. They both fully take into account the medium modifications due to the interaction induced correlations. The two methods are in good agreement up to the maximum gap density and both point to a slight reduction with respect to the standard BCS value. In fact, the maximum gap is about $2.5 \text{MeV}$ at $k_F \sim 0.8 \text{fm}^{-1}$ in BCS and 2.3--$2.4 \text{MeV}$ at $k_F \sim 0.6 \text{fm}^{-1}$ in correlated matter. At higher densities the Quantum Monte Carlo gap becomes close to BCS. In general, the computed medium polarization effects are much smaller than those previously estimated within \emph{all theories}. Truncations of Argonne $v_{8^\prime}$ to simpler forms give the same gaps in BCS, provided the truncated potentials have been refitted to the same NN data set. Differences among the models appear in the correlated theories, most of the reduction being attributable to the tensor force. The three--nucleon interaction provides an additional increase of the gap of about 0.35 MeV.

Gas falling quasi-spherically onto a black hole forms an inner accretion disc if its specific angular momentum $l$ exceeds $\lmin\sim r_gc$ where $r_g$ is the Schwarzschild radius. The standard disc model assumes $l\gg\lmin$. We argue that, in many black-hole sources, accretion flows have angular momenta just above the threshold for disc formation, $l\simgt\lmin$, and assess the accretion mechanism in this regime. In a range $\lmin<l<\lcr$, a small-scale disc forms in which gas spirals fast into the black hole without any help of horizontal viscous stresses. Such an `inviscid' disc, however, interacts inelastically with the feeding infall. The disc-infall interaction determines the dynamics and luminosity of the accretion flow. The inviscid disc radius can be as large as $14 r_g$, and the energy release peaks at $2r_g$. The disc emits a Comptonized X-ray spectrum with a break at $\sim 100$ keV. This accretion regime is likely to take place in wind-fed X-ray binaries and is also possible in active galactic nuclei.

In wind-fed X-ray binaries the accreting matter is Compton cooled and falls freely onto the compact object. The matter has a modest angular momentum $l$ and accretion is quasi-spherical at large distances from the compact object. Initially small non-radial velocities grow in the converging supersonic flow and become substantial in the vicinity of the accretor. The streamlines with $l>(GMR_*)^{1/2}$ (where $M$ and $R_*$ are the mass and radius of the compact object) intersect outside $R_*$ and form a two-dimensional caustic which emits X-rays. The streamlines with low angular momentum, $l<(GMR_*)^{1/2}$, run into the accretor. If the accretor is a neutron star, a large X-ray luminosity results. We show that the distribution of accretion rate/luminosity over the star surface is sensitive to the angular momentum distribution of the accreting matter. The apparent luminosity depends on the side from which the star is observed and can change periodically with the orbital phase of the binary. The accretor then appears as a `Moon-like' X-ray source.

PSR B1259-63 is the only known binary system with a radio pulsar from which the non-pulsed radio and X-ray emission was detected. The companion star in this system is a Be star SS 2883. A rapidly rotating radio pulsar is expected to produce a wind of relativistic particles. Be stars are known to produce highly asymmetric mass loss. Due to the interaction of the pulsar wind and the Be star wind the system of two shocks between the pulsar and the Be star forms. In this paper we show that the observed non-pulsed radio emission from the system is a result of the synchrotron emission of the relativistic particles in the outflow beyond the shock wave and that the non-pulsed X-ray emission is due to the inverse Compton scattering of the Be star photons on this particles.

Observed spectra of Active Galactic Nuclei (AGN) and luminous X-ray binaries in our Galaxy suggest that both hot (~10^9 K) and cold (~10^6 K) plasma components exist close to the central accreting black hole. Hard X-ray component of the spectra is usually explained by Compton upscattering of optical/UV photons from optically thick cold plasma by hot electrons. Observations also indicate that some of these objects are quite efficient in converting gravitational energy of accretion matter into radiation. Existing theoretical models have difficulties in explaining the two plasma components and high intensity of hard X-rays. Most of the models assume that the hot component emerges from the cold one due to some kind of instability, but no one offers a satisfactory physical explanation for this. Here we propose a solution to these difficulties that reverses what was imagined previously: in our model the hot component forms first and afterward it cools down to form the cold component. In our model, accretion flow has initially a small angular momentum, and thus it has a quasi-spherical geometry at large radii. Close to the black hole, the accreting matter is heated up in shocks that form due to the action of the centrifugal force. The hot post-shock matter is very efficiently cooled down by Comptonization of low energy photons and condensates into a thin and cold accretion disk. The thin disk emits the low energy photons which cool the hot component.

Consider binary system with a millisecond pulsar ejecting relativistic particles and an optical star emitting soft photons with the energy $\omega \simeq 1-10$ eV. These low-energy photons are scattered by the relativistic electrons and positrons of the pulsar's wind. The scattered photons forms a wide spectrum from hard X-ray band up to gamma band $\epsilon \simeq 1-1000$GeV >.When the pulsar wind is isotropic the luminosity of gamma radiation $L_\gamma =L_\gamma (\psi)$ depends heavily on the angle $\psi $ between the directions to the optical star and to the observer from the pulsar. During the orbital motion this angle varies periodically giving rise to the periodical change of the observed intensity of the gamma radiation and its spectrum. We calculated the spectral shape of the scattered hard photons. Under the assumption that the energy losses of the relativistic particles are small we receive analytical formulas. We apply our results to the binary system PSR B1259-63 and show that if the wind from the Be star is accounted for then it is possible to reproduce the observed spectrum.

A new statistical approach is presented to study the thermal instability process of optically thin unmagnetized plasma. In this approach the time evolution of mass distribution function over temperature is calculated. This function characterizes the statistical properties of the multiphase medium of arbitrary spaced three-dimensional structure of arbitrary temperature perturbations. We construct our theory under the isobarical condition (P=const over space), which is satisfied in the short wavelength limit. The developed theory is illustrated in the case of thermal instability of a slowly expanding interstellar cloud. Numerical solutions of equations of the statistical theory are constucted and compared with hydrodynamical solutions. The results of both approaches are identical in the short wavelength range when the isobarity condition is satisfied. Also the limits of applicability of the statistical theory are estimated. The possible evolution of initial spectrum of perturbations is discussed. The proposed theory and numerical models can be relevant to the formation of the two-phases medium in the ~1pc region around quasars. Then small warm (T~10000K) clouds are formed as the result of thermal instability in an expanded gas fragment, which is a product of either a star-star or star-accretion disk collision.

As first shown by Shvartsman (1970), a neutron star accreting close to the Eddington limit must acquire a positive charge in order for electrons and protons to move at the same speed. The resulting electrostatic field may contribute to accelerating positrons produced near the star surface in conjunction with the radiative force. We reconsider the balance between energy gains and losses, including inverse Compton (IC), bremsstrahlung and non--radiative scatterings. It is found that, even accounting for IC losses only, the maximum positron energy never exceeds $\approx 400$ keV. The electrostatic field alone may produce energies $\approx 50$ keV at most. We also show that Coulomb collisions and annihilation with accreting electrons severely limit the number of positrons that escape to infinity.

It is generally supposed that when the "compactness" $l \equiv L\sigma_T/(r m_e c^3)$ in photons above the pair-production threshold is large, few $\gamma$-rays can escape. We demonstrate that even when $l \gg 1$, if the high energy and low energy photons are produced in geometrically-separated regions, many of the $\gamma$-rays can, in fact, escape. Pair-production along a thin surface separating the two sources creates enough Compton optical depth to deflect most of the low energy photons away from the high energy ones. Those few low-energy photons which penetrate the shielding surface are reduced in opacity by advection to large distance and small density, by relativistic beaming along the inner edge of the surface, and by Compton upscattering to higher energies inside the surface. The pairs in this surface flow outward relativistically, forming a structure resembling a pair-dominated mildly relativistic jet.