Neutrino-nucleus reactions cross sections, obtained for neutrino energies in the range $\varepsilon_{\nu}\leq 100-120$ MeV (low- and intermediate-energy range), which refer to promising neutrino detection targets of current terrestrial neutrino experiments, are presented and discussed. At first, we evaluated original cross sections for elastic scattering of neutrinos produced from various astrophysical and laboratory neutrino sources with the most abundant Cd isotopes $^{112}$Cd, $^{114}$Cd and $^{116}$Cd. These isotopes constitute the main material of the COBRA detector aiming to search for neutrinoless double beta decay events and neutrino-nucleus scattering events at the Gran Sasso laboratory (LNGS). The coherent $\nu$-nucleus reaction channel addressed with emphasis here, dominates the neutral current $\nu$-nucleus scattering, events of which have only recently been observed for a first time in the COHERENT experiment at Oak Ridge. Subsequently, simulated $\nu$-signals expected to be recorded at Cd detectors are derived through the application of modern simulation techniques and employment of reliable neutrino distributions of astrophysical $\nu$-sources (as the solar, supernova and Earth neutrinos), as well as laboratory neutrinos (like the reactor neutrinos, the neutrinos produced from pion-muon decay at rest and the $\beta$-beam neutrinos produced from the acceleration of radioactive isotopes at storage rings as e.g. at CERN).

In this work, we simulate $\gamma$-rays created in the hadronic jets of the compact object in binary stellar systems known as microquasars. We utilize as main computational tool the $3$-D relativistic magneto-hydro-dynamical code PLUTO combined with in house derived codes. Our simulated experiments refer to the SS433 X-ray binary, a stellar system in which hadronic jets have been observed. We examine two new model configurations that employ hadron-based emission mechanisms. The simulations aim to explore the dependence of the $\gamma$-ray emissions on the dynamical as well as the radiative properties of the jet (hydrodynamic parameters of the mass-flow density, gas-pressure, temperature of the ejected matter, high energy proton population inside the jet plasma, etc.). The results of the two new scenarios of initial conditions for the micro-quasar stellar system studied, are compared to those of previously considered scenarios.

The hadronic jets in a micro-quasar stellar system are modelled with the relativistic hydrocode PLUTO. We focus on neutrino emission from such jets produced by fast proton (non-thermal) collisions on thermal ones within the hadronic jet. We adopt a semi-analytical approximation for the description of the secondary particles produced from p-p collisions and develop appropriate algorithms using the aforementioned injected protons as input. As a concrete example we consider the SS-433 X-ray binary system for which several observations have been performed the last decades. In contrast to the pre-set distribution of the fast protons along the jet employed in our previous works, in the present paper we simulated it by using a power-law fast proton distribution along the PLUTO hydro-code. This distribution gradually sweeps aside the surrounding winds, during the jet advance through the computational grid. As a first step, in the present work the neutrino energy spectrum is extracted from the model jet, facilitating a range of potential dynamical simulations in currently interesting microquasar jet systems.

In this work, the benefits of the phase fitting technique are embedded in high order discrete Lagrangian integrators. The proposed methodology creates integrators with zero phase lag in a test Lagrangian in a similar way used in phase fitted numerical methods for ordinary differential equations. Moreover, an efficient method for frequency evaluation is proposed based on the eccentricities of the moving objects. The results show that the new method dramatically improves the accuracy and total energy behaviour in Hamiltonian systems. Numerical tests for the 2-body problem with ultra high eccentricity up to 0.99 for 1000000 periods and to the Henon-Heiles Hamiltonian system with chaotic behaviour, show the efficiency of the proposed approach.