In this talk we will discuss the recent advances in describing heavy-quark dynamics in the quark-gluon plasma (QGP), which evolves hydrodynamically. Special emphasis is put on the collective flow of the heavy-quarks with the medium constituents, for which we present our latest results obtained within the MC@sHQ+EPOS2 model at $\sqrt{s}=5$~TeV.

The ultracompact low-mass X-ray binary 4U 1820-30 situated in the globular cluster NGC 6624 has an orbital period of only $\approx$11.4 min which likely implies a white dwarf companion. The observed X-ray bursts demonstrate a photospheric radius expansion phase and therefore are believed to reach the Eddington luminosity allowing us to estimate the mass and the radius of the neutron star (NS) in this binary. Here we re-analyse all Rossi X-ray Timing Explorer observations of the system and confirm that almost all the bursts took place during the hard persistent state of the system. This allows us to use the recently developed direct cooling tail method to estimate the NS mass and radius. However, because of the very short, about a second, duration of the cooling tail phases that can be described by the theoretical atmosphere models, the obtained constraints on the NS radius are not very strict. Assuming a pure helium NS atmosphere we found that the NS radius is in the range 10-12 km, if the NS mass is below 1.7 $M_\odot$, and in a wider range of 8-12 km for a higher 1.7-2.0 $M_\odot$ NS mass. The method also constrains the distance to the system to be 6.5$\pm$0.5 kpc, which is consistent with the distance to the cluster. For the solar composition atmosphere, the NS parameters are in strong contradiction with the generally accepted range of possible NS masses and radii.

We present an analysis of the comparison between 2nd and 3rd flow harmonics of $\pi$, protons and D mesons produced in ultrarelativistic heavy ion collisions. We advocate that such an analysis could turn out to be a good strategy in order to quantify the off-equilibrium dynamics of heavy flavor at lower transverse momentum $p_T$.

Determining neutron star (NS) radii and masses can help to understand the properties of matter at supra-nuclear densities. Thermal emission during thermonuclear X-ray bursts from NSs in low-mass X-ray binaries provides a unique opportunity to study NS parameters, because of the high fluxes, large luminosity variations and the related changes in the spectral properties. The standard cooling tail method uses hot NS atmosphere models to convert the observed spectral evolution during cooling stages of X-ray bursts to the Eddington flux F_Edd and the stellar angular size \Omega. These are then translated to the constraints on the NS mass M and radius R. Here we present the improved, direct cooling tail method that generalises the standard approach. First, we adjust the cooling tail method to account for the bolometric correction to the flux. Then, we fit the observed dependence of the blackbody normalization on flux with a theoretical model directly on the M-R plane by interpolating theoretical dependences to a given gravity, hence ensuring only weakly informative priors for M and R instead of F_Edd and \Omega. The direct cooling method is demonstrated using a photospheric radius expansion burst from SAX J1810.8--2609, which has happened when the system was in the hard state. Comparing to the standard cooling tail method, the confidence regions are shifted by 1 sigma towards larger radii, giving R=11.5-13.0 km at M=1.3-1.8 M_sun for this NS.

Electromagnetic probes such as direct photons and dileptons are crucial to study the properties of a Quark-Gluon Plasma (QGP) created in heavy ion collisions. Based on the (3+1)-dimensional event-by-event viscous hydrodynamic model EPOS 3.102, we calculated the anisotropic emission of thermal photons and dileptons in Au+Au collisions at the Relativistic Heavy Ion Collider (RHIC) energy $\sqrt{s_{NN}}=200$~GeV. Thermal emissions from both QGP phase and hadronic gas phase are considered, with AMY rate for photons and Lattice QCD based rates for dileptons in QGP phase. For emission from hadron gas phase, rates based vector meson dominant model are used for both photons and dileptons. Non-thermal contribution to direct photons is calculated with next to leading order QCD. STAR cocktail data are directly used for non-thermal contribution to dileptons. With the same space-time evolution of the collision systems, the two penetrating probes, photons and dileptonsshow some consistency, ie, the emission of both thermal photons and dileptons are underestimated, compared with the transverse momentum spectra of direct photons measured by PHENIX collaboration and the invariant mass spectra of dileptons measured by STAR collaboration, for all centrality classes. With a good constraint of anisotropy of the plasma via the elliptic flow and triangular flow of charged hadrons, the resulted elliptic flow and triangular flow of direct photons agree with PHENIX measurements reasonably well. Thus we made predictions to the flows of thermal dileptons. The elliptic flow of thermal dileptons is predicted larger than the available results from other models, and comparable to the STAR measurement referring to all dileptons (thermal + cocktail) for minimal bias collisions.

Heavy-quark observables in ultrarelativistic heavy-ion collisions, like the nuclear modification factor and the elliptic flow, give insight into the mechanisms of high-momentum suppression and low-momentum thermalization of heavy quarks. Here, we present a global study of these two observables within a coupled approach of the heavy-quark propagation in a realistic fluid dynamical medium, MC@sHQ+EPOS2, and compare to experimental data from RHIC and LHC experiments. The heavy quarks scatter elastically and inelastically with the quasiparticles of the quark-gluon plasma (QGP), which are represented consistently with the underlying equation of state. We examine two scenarios: first, we interpret the lattice QCD equation of state as a sum of partonic and hadronic contributions, and second, as a gas of massive partonic quasiparticles. It is observed that independent of their momentum the energy loss of heavy quarks depends strongly on how the lattice QCD equation of state is translated into degrees of freedom of the QGP.

Dileptons, as an electromagnetic probe, are crucial to study the properties of a Quark-Gluon Plasma (QGP) created in heavy ion collisions. We calculated the invariant mass spectra and the anisotropic emission of thermal dielectrons from Au+Au collisions at the Relativistic Heavy Ion Collider (RHIC) energy $\sqrt{s_{NN}}=200$~GeV based on EPOS3. This approach provides a realistic (3+1)-dimensional event-by-event viscous hydrodynamic description of the expanding hot and dense matter with a very particular initial condition, and a large set of hadron data and direct photons (besides $v_{2}$ and $v_{3}$ !) can be successfully reproduced. Thermal dilepton emission from both the QGP phase and the hadronic gas are considered, with the emission rates based on Lattice QCD and a vector meson model, respectively. We find that the computed invariant mass spectra (thermal contribution + STAR cocktail) can reproduce the measured ones from STAR at different centralities. Different compared to other model predictions, the obtained elliptic flow of thermal dileptons is larger than the STAR measurement referring to all dileptons. We observe a clear centrality dependence of thermal dilepton not only for elliptic flow $v_{2}$ but also for higher orders. At a given centrality, $v_{n}$ of thermal dileptons decreases monotonically with $n$ for $2\leq n\leq5$.

The observed large elliptic flow of direct photons is a puzzle in relativistic heavy ion physics. Our previous work provided a possibility to this puzzle with the delayed formation of the quark gluon plasma (QGP) in relativistic heavy ion collisions. In this work, we got the measured transverse momentum spectra, elliptic flow v2 and triangular flow v3 of direct photons from Au+Au collisions at the Relativistic Heavy Ion Collider (RHIC) energy \sqrt s_NN = 200 GeV with all centrality classes explained at the same time without any additional parameter to EPOS3, a hydrodynamic model which can successfully reproduce hadronic data such as rapidity distributions, transverse momenta, elliptic flows and triangular flows from variuos collision systems. The key point is that, EPOS3 has an initial space eccentricity similar to other models, however, a stronger radial flow and a weaker momentum eccentrcity made a good description to the flows of both hadrons and direct photons. The created QGP matter seems discrete than a soup.

Gaia will identify several 1e5 white dwarfs, most of which will be in the solar neighborhood at distances of a few hundred parsecs. Ground-based optical follow-up spectroscopy of this sample of stellar remnants is essential to unlock the enormous scientific potential it holds for our understanding of stellar evolution, and the Galactic formation history of both stars and planets.

Fallback disks are expected to form around young neutron stars. The presence of these disks can be revealed by their blackbody spectrum in the infrared, optical and UV bands. We present a re-reduction of the archival optical and infrared data of the Vela pulsar, together with the existing infrared and UV spectrum of Vela, and model their unpulsed components with the blackbody spectrum of a supernova debris disk. We invoke the quiescent disk solution of Sunyaev and Shakura for the description of the disk in the propeller stage and find the inner radius of the disk to be inside the light cylinder radius. We perform a high resolution X-ray analysis with XMM-Newton and find a narrow absorption feature at 0.57 keV which can be interpreted as K$_{\alpha}$ line of He-like oxygen (OVII). The strength of the line indicates an element over-abundance in our line of sight exceeding the amounts that would be expected from ISM. The spectral feature may originate from the pulsar wind nebula and may be partly caused by the reprocessed X-ray radiation by the fallback disk. We discuss the lower-than-3 braking index of Vela as partially due to the contribution of the propeller torques. Our results suggest that the pulsar mechanism can work simultaneously with the propeller processes and that the debris disks can survive the radiation pressure for at least $\sim 10^4$ years. As Vela is a relatively close object, and a prototypical pulsar, the presence of a disk, if confirmed, may indicate the ubiquity of debris disks around young neutron stars.

After a short presentation of the event generator EPOS, we discuss the production of heavy quarks and prompt photons which has been recently implemented. Whereas we have satisfying results for the charm, work on photons is still in progress.

We investigate the elliptic and the triangular flow of heavy mesons in ultrarelativistic heavy-ion collisions at RHIC and the LHC. The dynamics of heavy quarks is coupled to the locally thermalized and fluid dynamically evolving quark-gluon plasma. The elliptic flow of $D$ mesons and the centrality dependence measured at the LHC is well reproduced for purely collisional and bremsstrahlung interactions. Due to the event-by-event fluctuating initial conditions from the EPOS2 model, the $D$ meson triangular flow is predicted to be nonzero at $\sqrt{s}=200$ GeV and $\sqrt{s}=2.76$ TeV. We study the centrality dependence and quantify the contributions stemming from flow of the light bulk event and the hadronization process. The flow coefficients as response to the initial eccentricities behave differently for heavy mesons than for light hadrons due to their inertia. Higher-order flow coefficients of heavy flavor become important in order to quantify the degree of thermalization.

We investigate the charm-quark propagation in the QGP media produced in ultrarelativistic heavy-ion collisions at RHIC and the LHC. Purely collisional and radiative processes lead to a significant suppression of final $D$-meson spectra at high transverse momentum and a finite flow of heavy quarks inside the fluid dynamical evolution of the light partons. The $D$-meson nuclear modification factor and the elliptic flow are studied at two collision energies. We further propose to measure the triangular flow of $D$ mesons, which we find to be nonzero in non-central collisions.

We extend the Gunion Bertsch calculation of gluon radiation in single scattering to the case of finite mass quarks. This case applies to the radiative energy-loss of heavy quarks of intermediate energies propagating in a quark gluon plasma. We discuss more specifically the dead cone effect as well as the mass hierarchy of the collisional and radiative energy loss and provide some predictions for observables sensitive to the mass hierarchy of energy loss in ultrarelativistic heavy ion collisions.

In this paper we study the azimuthal correlations of heavy quarks in Pb+Pb collisions with $\sqrt{s}=2.76$ TeV at LHC. Due to the interaction with the medium heavy quarks and antiquarks are deflected from their original direction and the initial correlation of the pair is broadened. We investigate this effect for different transverse momentum classes. Low-momentum heavy-quark pairs lose their leading order back-to-back initial correlation, while a significant residual correlation survives at large momenta. Due to the larger acquired average deflection from their original directions the azimuthal correlations of heavy-quark pairs are broadened more efficiently in a purely collisional energy loss mechanism compared to including radiative corrections. This discriminatory feature survives when next-to-leading order production processes are included.

Observables of heavy-quark azimuthal correlations in heavy-ion collisions are a new and promising tool for the investigation of the in-medium energy loss. We explore the potential of these observables to discriminate the collisional and radiative contributions within a hybrid EPOS+MC@sHQ transport approach.

We report the discovery of pulsations in the spectroscopic PG 1159 type pre-white dwarf SDSS J075415.12+085232.18. Analysis of the spectrum by Werner, Rauch and Kepler (2014) indicated Teff=120 000+/-10 000 K, log g=7.0+/-0.3, mass M=0.52+/-0.02 Msun, C/He=0.33 by number. We obtained time-series images with the SOAR 4.1 m telescope and 2.1 m Otto Struve telescope at McDonald Observatory and show the star is also a variable PG 1159 type star, with dominant period of 525 s.

The suppression of D mesons in non-central heavy-ion collisions is investigated. The anisotropy in collisions at finite impact parameter leads to an ordering of all-angle, in- and out-of-plane nuclear modification factors due to the different in-medium path lengths. Within our MC@sHQ+EPOS model of heavy-quark propagation in the QGP we demonstrate that fluctuating initial conditions lead to an effective reduction of the energy loss of heavy quarks, which is seen in a larger nuclear modification factor at intermediate and high transverse momenta. The elliptic flow at small transverse momenta is reduced.

Only a few percent of cool, old white dwarfs (WDs) have infrared excesses interpreted as originating in small hot disks due to the infall and destruction of single asteroids that come within the star's Roche limit. Infrared excesses at 24 micron were also found to derive from the immediate vicinity of younger, hot WDs, most of which are still central stars of planetary nebulae (CSPN). The incidence of CSPN with this excess is 18%. The Helix CSPN, with a 24 micron excess, has been suggested to have a disk formed from collisions of Kuiper belt-like objects (KBOs). In this paper, we have analyzed an additional sample of CSPN to look for similar infrared excesses. These CSPN are all members of the PG 1159 class and were chosen because their immediate progenitors are known to often have dusty environments consistent with large dusty disks. We find that, overall, PG 1159 stars do not present such disks more often than other CSPN, although the statistics (5 objects) are poor. We then consider the entire sample of CSPN with infrared excesses, and compare it to the infrared properties of old WDs, as well as cooler post-AGB stars. We conclude with the suggestion that the infrared properties of CSPN more plausibly derive from AGB-formed disks rather than disks formed via the collision of KBOs, although the latter scenario cannot be ruled out. We finally remark that there seems to be an association between CSPN with a 24 micron excess and confirmed or possible binarity of the central star.

We investigate how the possible existence of hadronic bound states above the deconfinement transition temperature $T_c$ affects heavy-quark observables like the nuclear modification factor, the elliptic flow and azimuthal correlations. Lattice QCD calculations suggest that above $T_c$ the effective degrees of freedom might not be exclusively partonic but that a certain fraction of hadronic degrees of freedom might already form at higher temperatures. This is an interesting questions by itself but also has a strong influence on other probes of the strongly interacting matter produced in ultrarelativistic heavy-ion collisions. A substantial fraction of hadronic bound states above $T_c$ reduces on average the interaction of the heavy quarks with colored constituents of the medium. We find that all studied observables are highly sensitive to the active degrees of freedom in the quark-gluon plasma.

Disk accretion onto stars leads to the formation of a Boundary Layer (BL) near the stellar surface where the disk makes contact with the star. Albeit a large fraction of the total luminosity of the system originates from this tiny layer connecting the accretion disk and the accreting object, its structure has not been fully understood yet. It is the aim of this work, to obtain more insight into the Boundary Layer around the white dwarf in compact binary systems. There are still many uncertainties concerning the extent and temperature of the BL and the rotation rate of the white dwarf. We perform numerical hydrodynamical simulations, where the problem is treated in a one-dimensional, radial approximation (slim disk). The turbulence is described by the alpha parameter viscosity. We include both cooling from the disk surfaces and radial radiation transport. The radiation energy is treated in the one-temperature approximation. For a given M_dot our results show a strong dependence on the stellar mass and rotation rate. Both, the midplane and the effective temperature rise considerably with increasing stellar mass or decreasing stellar rotation rate. Our simulations further show, that the radiation energy and pressure are indeed important in the BL. However, some models show a low optical depth in the BL, making it necessary to find a better representation for optically thin regions. The combination of a high mass and a small radius, such as in white dwarfs, can lead to an enormous energy release in the BL, provided the WD rotates slowly. Since the radial extent of BLs is typically very small (about 0.02 to 0.05 R_star), this leads to surface temperatures of a few hundred thousand Kelvin. All of our models showed subsonic infall velocities with Mach numbers of < 0.4 at most.

Radio detection of cosmic-ray-induced air showers has come to a flight the last decade. Along with the experimental efforts, several theoretical models were developed. The main radio-emission mechanisms are established to be the geomagnetic emission due to deflection of electrons and positrons in Earth's magnetic field and the charge-excess emission due to a net electron excess in the air shower front. It was only recently shown that Cherenkov effects play an important role in the radio emission from air showers. In this article we show the importance of these effects to extract quantitatively the position of the shower maximum from the radio signal, which is a sensitive measure for the mass of the initial cosmic ray. We also show that the relative magnitude of the charge-excess and geomagnetic emission changes considerably at small observer distances where Cherenkov effects apply.

We present a macroscopic calculation of coherent electro-magnetic radiation from air showers initiated by ultra-high energy cosmic rays, based on currents obtained from three-dimensional Monte Carlo simulations of air showers in a realistic geo-magnetic field. We discuss the importance of a correct treatment of the index of refraction in air, given by the law of Gladstone and Dale, which affects the pulses enormously for certain configurations, compared to a simplified treatment using a constant index. We predict in particular a geomagnetic Cherenkov radiation, which provides strong signals at high frequencies (GHz), for certain geometries together with "normal radiation" from the shower maximum, leading to a double peak structure in the frequency spectrum. We also provide some information about the numerical procedures referred to as EVA 1.0.

In this contribution, we present some predictions for the production of D and B mesons in ultrarelativistic heavy ion collisions at RHIC and LHC energies and confront them with experimental results obtained so far by the STAR, PHENIX, ALICE and CMS collaborations. We next discuss some preliminary results obtained with an improved description of the medium based on EPOS initial conditions, and its possible implications on the nuclear modification factor and on the elliptic flow of heavy quarks.

The physics of radio emission from cosmic-ray induced air showers is shortly summarized. It will be shown that the radio signal at different distances from the shower axis provides complementary information on the longitudinal shower evolution, in particular the early part, and on the distribution of the electrons in the shower core. This complements the information obtained from surface, fluorescence, and muon detectors and is very useful in getting a comprehensive picture of an air shower.