• The gyrokinetic toroidal code (GTC) has been upgraded for global simulations by coupling the core and scrape-off layer (SOL) regions across the separatrix with field-aligned particle-grid interpolations. A fully kinetic particle pusher for high frequency waves (ion cyclotron frequency and beyond) and a guiding center pusher for low frequency waves have been implemented using cylindrical coordinates in a global toroidal geometry. The two integrators correctly capture the particle orbits and agree well with each other, conserving energy and canonical angular momentum. As a verification and application of this new capability, ion guiding center simulations have been carried out to study ion orbit losses at the edge of the DIII-D tokamak for single null magnetic separatrix discharges. The ion loss conditions are examined as a function of the pitch angle for cases without and with a radial electric field. The simulations show good agreement with past theoretical results and with experimentally observed feature in which high energy ions flow out along the ion drift orbits and then hit the divertor plates. A measure of the ion direct orbit loss fraction shows that the loss fraction increases with the ion energy for DIII-D in the initial velocity space. Finally, as a further verification of the capability of the new code, self-consistent simulations of zonal flows in the core region of the DIII-D tokamak were carried out. All DIII-D simulations were perfomed in the absence of turbulence.
  • One-dimensional particle-in-cell simulation is used to simulate the capacitively coupled argon plasma for a range of driving frequency from 13.56 MHz to 100 MHz. The argon chemistry set can, selectively, include two metastable levels enabling multi-step ionization and metastable pooling. The results show that the plasma density decreases when metastable atoms are included with higher discrepancy at higher excitation frequency. The contribution of multistep ionization to overall density increases with excitation frequency. The electron temperature increases with the inclusion of metastable atoms and decreases with excitation frequency. At lower excitation frequency, the density of Ar** (3p5 4p, 13.1 eV) is higher than Ar* (3p5 4s, 11.6 eV), whereas, at higher excitation frequencies the Ar* (3p5 4s, 11.6 eV) is the dominant metastable atom. The metastable and electron temperature profile evolve from a parabolic profile at lower excitation frequency to a saddle type profile at higher excitation frequency. With metastable, the electron energy distribution function (EEDF) changes its shape from Druyvesteyn type, at low excitation frequency, to bi-Maxwellian, at high frequency plasma excitation, however a three-temperature EEDF is observed without metastable atoms.
  • We present criteria for the photometric selection of M-dwarfs using all-sky photometry, with a view to identifying M-dwarf candidates for inclusion in the input catalogues of upcoming all-sky surveys, including TESS and FunnelWeb. The criteria are based on Gaia, WISE and 2MASS all-sky photometry, and deliberately do not rely on astrometric information. In the lead-up to the availability of truly distance-limited samples following the release of Gaia DR2, this approach has the significant benefit of delivering a sample unbiased with regard to space velocity. Our criteria were developed by using Galaxia synthetic galaxy model predictions to evaluate both M-dwarf completeness and false-positive detections (i.e. non-M-dwarf contamination rates). In addition to the previously known sensitivity of J-H colour for giant-dwarf discrimination at cool temperatures, we find the WISE W1-W2 colour is also effective at discriminating M-dwarfs from cool giants. We have derived two sets of Gaia G > 14.5 criteria - a "high-completeness" set that contains 78,340 stars, of which 30.7-44.4% are expected to be M-dwarfs and contains 99.3% of the total number of expected M-dwarfs; and a "low-contamination" set that prioritises the stars most likely to be M-dwarfs at a cost of a reduction in completeness. This subset contains 40,505 stars and is expected to be comprised of 58.7-64.1% M-dwarfs, with a completeness of 98%. Comparison of the high-completeness set with the TESS Input Catalogue has identified 234 stars not currently in that catalogue, which preliminary analysis suggests could be useful M-dwarf targets for TESS. We also compared the criteria to selection via absolute magnitude and a combination of both methods. We found that colour selection in combination with an absolute magnitude limit provides the most effective way of selecting M-dwarfs en masse.
  • This article reports on the feasibility of testing of the symmetry under reversal in time in a purely leptonic system constituted by positronium atoms using the J-PET detector. The present state of T symmetry tests is discussed with an emphasis on the scarcely explored sector of leptonic systems. Two possible strategies of searching for manifestations of T violation in non-vanishing angular correlations of final state observables in the decays of metastable triplet states of positronium available with J-PET are proposed and discussed. Results of a pilot measurement with J-PET and assessment of its performance in reconstruction of three-photon decays are shown along with an analysis of its impact on the sensitivity of the detector for the determination of T -violation sensitive observables.
  • The overlap between the spectroscopic Galactic Archaeology with HERMES (GALAH) survey and $Gaia$ provides a high-dimensional chemodynamical space of unprecedented size. We present a first analysis of a subset of this overlap, of 7066 dwarf, turn-off, and sub-giant stars. These stars have spectra from the GALAH survey and high parallax precision from the $Gaia$ DR1 Tycho-$Gaia$ Astrometric Solution. We investigate correlations between chemical compositions, ages, and kinematics for this sample. Stellar parameters and elemental abundances are derived from the GALAH spectra with the spectral synthesis code Spectroscopy Made Easy. We determine kinematics and dynamics, including action angles, from the $Gaia$ astrometry and GALAH radial velocities. Stellar masses and ages are determined with Bayesian isochrone matching, using our derived stellar parameters and absolute magnitudes. Our stars span mainly metallicities $-0.7 < \mathrm{[Fe/H]} < +0.5$. We report measurements of Li, C, O, Na, Mg, Al, Si, K, Ca, Sc, Ti, V, Cr, Mn, Co, Ni, Cu, Zn, Y, and Ba and note that for some elements (Li, O, Al, and Fe), we employ non-LTE calculations. These elements represent several nucleosynthesis channels including $\alpha$, iron-peak, and s-process. We find strong correlations between stellar age and abundance, including for $\alpha$-elements and the s-process elements Y and Ba. Furthermore, in the high-$\alpha$ metal-rich regime, age can be used in combination with chemistry to separate stars from the extended high-$\alpha$ sequence and the low-$\alpha$ one. We find that the low-$\alpha$ metal-poor stars in the solar neighbourhood have significantly higher angular momentum $L_z$ than the Sun, which implies that they are on eccentric orbits and originate from the outer disk. Stars of the high-$\alpha$-sequence on the other hand have lower $L_z$ than the Sun, implying they are visitors from the inner disk.
  • We experimentally identify the magnetostatic modes active for Brillouin light scattering in the optical whispering gallery modes of a yttrium iron garnet sphere. Each mode is identified by magnetic field dispersion of ferromagnetic-resonance spectroscopy and coupling strength to the known field distribution of the microwave drive antenna. Our optical measurements confirm recent predictions that higher-order magnetostatic modes can also generate optical scattering, according to the selection rules derived from the axial symmetry. From this we summarize the selection rules for Brillouin light scattering. We give experimental evidence that the optomagnonic coupling to non-uniform magnons can be higher than that of the uniform Kittel mode.
  • Microscopic spin-dependent charge transfer, backtransfer and scattering processes at interfaces are crucial for a broad range of dynamic spin phenomena in ferromagnetic heterostructures. We disentangle the fundamental quantum processes in optically induced spin transfer across a model epitaxial Co/Cu(001) interface by means of a direct comparison between femtosecond time-resolved interface-sensitive magnetization-induced second harmonic generation and parameter-free ab-initio time-dependent density functional theory. Due to agreement between theory and experiment, we identify transfer of majority electrons from the Co to the Cu, and simultaneous minority electron backtransfer from the Cu to the Co interface layers, as driving the spin dynamics in the first 30 fs after laser excitation. We separate this interfacial spin transfer from subsequently dominating spin flip processes induced by spin-orbit coupling, which reduce the total spin polarization in the heterostructure. Our combined experimental-theoretical approach shows the potential of a predictive theoretical description to elucidate spin injection and accumulation at ferromagnet/metal interfaces in conjunction with interface-sensitive magneto-optical experiments.
  • The possibility to separate signals caused by 511 keV photons created in annihilation of electron-positron pairs and the so-called prompt photons from nuclei de- excitation is investigated. It could potentially be used to improve the quality of reconstructed images in the J-PET scanner in 2+1 photon tomography. Firstly, a research is conducted for several radioisotopes that decay via \b{eta}+ decay followed by de-excitation of an excited nucleus. Efficiency, purity and false positive rate are calculated for each isotope as a function of energy deposited threshold, with a hypothesis that signals caused by 511 keV photons deposit smaller values of energy than 1 z 13the selected threshold, while prompt photons deposit larger energy than the threshold. Analysis of the results accompanied with physical properties of radioisotopes suggests using 44 Sc, which is the most promising candidate for medical applications. With the use of GATE and J-POS simulation software, in-phantom scattering was introduced and the best energy deposited threshold value was estimated to be approximately 375 keV. It corresponds to almost 100% efficiency for 511 keV signals, 75% purity for 511 keV photons, and approximately 70% efficiency and purity for prompt photons.
  • A recently proposed exchange-correlation functional (Dewhurst et al. 10.1021/acs.jctc.7b01049) with in density functional theory, which ensures that the exchange-correlation magnetic field is source-free, is shown to give non-zero internal spin-torque. This spin-torque is identically zero for all conventional local and semi-local functionals. Extension of this source-free functional to the time domain is used to study the effect of the internal spin-torque on the laser induced spin-dynamics in bulk Co, Ni and interfaces of these metals with Pt. It is shown that the internal spin-torque contribute significantly to spin-dynamics only when the magneto crystalline anisotropy energy is small, as in the case of cubic bulk materials. For surfaces or interfaces, where the anisotropy energy is large, these torques are too small to cause any significant precession of spins in early times ($<$ 100fs). Further more it is shown that the spin-dynamics caused by the internal spin-torque is slow compared to the inter-site spin transfer and spin-orbit mediated spin-flips.
  • The normalized $^{238}$U(n,f)/$^{235}$U(n,f) cross section ratio has been measured using the NIFFTE fission Time Projection Chamber from the reaction threshold to $30$~MeV. The fissionTPC is a two-volume MICROMEGAS time projection chamber that allows for full three-dimensional reconstruction of fission-fragment ionization profiles from neutron-induced fission. The measurement was performed at the Los Alamos Neutron Science Center, where the neutron energy is determined from neutron time-of-flight. The $^{238}$U(n,f)/$^{235}$U(n,f) ratio reported here is the first cross section measurement made with the fissionTPC, and will provide new experimental data for evaluation of the $^{238}$U(n,f) cross section, an important standard used in neutron-flux measurements. Use of a development target in this work prevented the determination of an absolute normalization, to be addressed in future measurements. Instead, the measured cross section ratio has been normalized to ENDF/B-VIII.$\beta$5 at 14.5 MeV.
  • In recent years the Tsallis statistics is gaining popularity in describing charged particle produc- tion and their properties, in particular pT spectra and the multiplicities in high energy particle collisions. Motivated by its success, an analysis of the LHC data of proton-proton collisions at ener- gies ranging from 0.9 TeV to 7 TeV in different rapidity windows for charged particle multiplicities has been done. A comparative analysis is performed in terms of the Tsallis distribution, the Gamma distribution and the shifted-Gamma distribution. An interesting observation on the inapplicability of these distributions at sqrt{s}=7 TeV in the lower rapidity windows is intriguing. The non-extensive nature of the Tsallis statistics is studied by determining the entropic index and its energy depen- dence. The analysis is extrapolated to predict the multiplicity distribution at sqrt{s}=14 TeV for one rapidity window, |y| < 1.5 with the Tsallis function.
  • We use data from the Radial Velocity Experiment (RAVE) and the Tycho-Gaia astrometric solution catalogue (TGAS) to compute the velocity fields yielded by the radial (VR), azimuthal (Vphi) and vertical (Vz) components of associated Galactocentric velocity. We search in particular for variation in all three velocity components with distance above and below the disc midplane, as well as how each component of Vz (line-of-sight and tangential velocity projections) modifies the obtained vertical structure. To study the dependence of velocity on proper motion and distance we use two main samples: a RAVE sample including proper motions from the Tycho-2, PPMXL and UCAC4 catalogues, and a RAVE-TGAS sample with inferred distances and proper motions from the TGAS and UCAC5 catalogues. In both samples, we identify asymmetries in VR and Vz. Below the plane we find the largest radial gradient to be dVR / dR = -7.01+- 0.61 km\s kpc, in agreement with recent studies. Above the plane we find a similar gradient with dVR / dR= -9.42+- 1.77 km\s kpc. By comparing our results with previous studies, we find that the structure in Vz is strongly dependent on the adopted proper motions. Using the Galaxia Milky Way model, we demonstrate that distance uncertainties can create artificial wave-like patterns. In contrast to previous suggestions of a breathing mode seen in RAVE data, our results support a combination of bending and breathing modes, likely generated by a combination of external or internal and external mechanisms.
  • TIFR Near Infrared Imaging Camera-II is a closed-cycle Helium cryo-cooled imaging camera equipped with a Raytheon 512 x 512 pixels InSb Aladdin III Quadrant focal plane array having sensitivity to photons in the 1-5 microns wavelength band. In this paper, we present the performance of the camera on the newly installed 3.6-m Devasthal Optical Telescope (DOT) based on the calibration observations carried out during 2017 May 11-14 and 2017 October 7-31. After the preliminary characterization, the camera has been released to the Indian and Belgian astronomical community for science observations since 2017 May. The camera offers a field-of-view of ~86.5 arcsec x 86.5 arcsec on the DOT with a pixel scale of 0.169 arcsec. The seeing at the telescope site in the near-infrared bands is typically sub-arcsecond with the best seeing of ~0.45 arcsec realized in the near-infrared K-band on 2017 October 16. The camera is found to be capable of deep observations in the J, H and K bands comparable to other 4-m class telescopes available world-wide. Another highlight of this camera is the observational capability for sources up to Wide-field Infrared Survey Explorer (WISE) W1-band (3.4 microns) magnitudes of 9.2 in the narrow L-band (nbL; lambda_{cen} ~3.59 microns). Hence, the camera could be a good complementary instrument to observe the bright nbL-band sources that are saturated in the Spitzer-Infrared Array Camera ([3.6] <= 7.92 mag) and the WISE W1-band ([3.4] <= 8.1 mag). Sources with strong polycyclic aromatic hydrocarbon (PAH) emission at 3.3 microns are also detected. Details of the observations and estimated parameters are presented in this paper.
  • Momentum-resolved scattering experiments with laser-cooled atomic targets have been performed since almost two decades with MOTRIMS (Magneto-Optical Trap Recoil Ion Momentum Spectroscopy) setups. Compared to experiments with gas-jet targets, MOTRIMS features significantly lower target temperatures allowing for an excellent recoil ion momentum resolution. However, the coincident and momentum-resolved detection of electrons was long rendered impossible due to incompatible magnetic field requirements. Here we report on a novel experimental approach which is based on an all-optical $^6$Li atom trap that -- in contrast to magneto-optical traps -- does not require magnetic field gradients in the trapping region. Atom temperatures of about 2 mK and number densities up to 10$^9$ cm$^{-3}$ make this trap ideally suited for momentum-resolved electron-ion coincidence experiments. The overall configuration of the new trap is very similar to conventional magneto-optical traps. It mainly requires small modifications of laser beam geometries and polarization which makes it easily implementable in other existing MOTRIMS experiments.
  • The Jagiellonian Positron Emission Tomograph (J-PET) project carried out in the Institute of Physics of the Jagiellonian University is focused on construction and tests of the first prototype of PET scanner for medical diagnostic which allows for the simultaneous 3D imaging of the whole human body using organic scintillators. The J-PET prototype consists of 192 scintillator strips forming three cylindrical layers which are optimized for the detection of photons from the electron-positron annihilation with high time- and high angular-resolutions. In this article we present time calibration and synchronization of the whole J-PET detection system by irradiating each single detection module with a 22Na source and a small detector providing common reference time for synchronization of all the modules.
  • Positron Annihilation Lifetime Spectroscopy (PALS) has shown to be a powerful tool to study the nanostructures of porous materials. Positron Emissions Tomography (PET) are devices allowing imaging of metabolic processes e.g. in human bodies. A newly developed device, the J-PET (Jagiellonian PET), will allow PALS in addition to imaging, thus combining both analyses providing new methods for physics and medicine. In this contribution we present a computer program that is compatible with the J-PET software. We compare its performance with the standard program LT 9.0 by using PALS data from hexane measurements at different temperatures. Our program is based on an iterative procedure, and our fits prove that it performs as good as LT 9.0.
  • The influence of magnetic fields (B-fields) in the formation and evolution of bipolar bubbles, due to the expanding ionization fronts (I-fronts) driven by the Hii regions that are formed and embedded in filamentary molecular clouds, has not been well-studied yet. In addition to the anisotropic expansion of I-fronts into a filament, B-fields are expected to introduce an additional anisotropic pressure which might favor expansion and propagation of I-fronts to form a bipolar bubble. We present results based on near-infrared polarimetric observations towards the central $\sim$8'$\times$8' area of the star-forming region RCW57A which hosts an Hii region, a filament, and a bipolar bubble. Polarization measurements of 178 reddened background stars, out of the 919 detected sources in the JHKs-bands, reveal B-fields that thread perpendicular to the filament long axis. The B-fields exhibit an hour-glass morphology that closely follows the structure of the bipolar bubble. The mean B-field strength, estimated using the Chandrasekhar-Fermi method, is 91$\pm$8 {\mu}G. B-field pressure dominates over turbulent and thermal pressures. Thermal pressure might act in the same orientation as those of B-fields to accelerate the expansion of those I-fronts. The observed morphological correspondence among the B-fields, filament, and bipolar bubble demonstrate that the B-fields are important to the cloud contraction that formed the filament, gravitational collapse and star formation in it, and in feedback processes. The latter include the formation and evolution of mid-infrared bubbles by means of B-field supported propagation and expansion of I-fronts. These may shed light on preexisting conditions favoring the formation of the massive stellar cluster in RCW57A.
  • Multiplicity distributions of charged particles produced in the $e^{+}e^{-}$ collisions at LEP2 energies ranging from 91 to 206 GeV in full phase space, are compared with predictions from Tsallis $q$-statistics and the recently proposed Weibull distribution functions.~The analysis uses data from two LEP experiments, L3 and OPAL.~It is shown that Tsallis $q$-statistics explains the data in a statistically acceptable manner in full phase space at all energies, while the Weibull distribution fails to explain the underlying properties of the data.~Modifications to the distributions proposed earlier, are applied to uncover manifold improvements in explaining the data characteristics.
  • We present results for the ratios of mean ($M_B$), variance ($\sigma_B^2$), skewness ($S_B)$ and kurtosis ($\kappa_B$) of net baryon-number fluctuations obtained in lattice QCD calculations with a physical light to strange quark mass ratio. Using next-to-leading order Taylor expansions in baryon chemical potential we find that qualitative features of these ratios closely resemble the corresponding experimentally measured cumulants ratios of net proton-number fluctuations for beam energies down to $\sqrt{s_{_{NN}}} \ge 19.6$ GeV. We show that the difference in cumulant ratios for the mean net baryon-number, $M_B/\sigma_B^2=\chi_1^B(T,\mu_B)/\chi_2^B(T,\mu_B)$ and the normalized skewness, $S_B\sigma_B=\chi_3^B(T,\mu_B)/\chi_2^B(T,\mu_B)$, naturally arises in QCD thermodynamics. Moreover, we establish a close relation between skewness and kurtosis ratios, $S_B\sigma_B^3/M_B=\chi_3^B(T,\mu_B)/\chi_1^B(T,\mu_B)$ and $\kappa_B\sigma_B^2=\chi_4^B(T,\mu_B)/\chi_2^B(T,\mu_B)$, valid at small values of the baryon chemical potential.
  • Deciphering the assembly history of the Milky Way is a formidable task, which becomes possible only if one can produce high-resolution chrono-chemo-kinematical maps of the Galaxy. Data from large-scale astrometric and spectroscopic surveys will soon provide us with a well-defined view of the current chemo-kinematical structure of the Milky Way, but will only enable a blurred view on the temporal sequence that led to the present-day Galaxy. As demonstrated by the (ongoing) exploitation of data from the pioneering photometric missions CoRoT, Kepler, and K2, asteroseismology provides the way forward: solar-like oscillating giants are excellent evolutionary clocks thanks to the availability of seismic constraints on their mass and to the tight age-initial-mass relation they adhere to. In this paper we identify five key outstanding questions relating to the formation and evolution of the Milky Way that will need precise and accurate ages for large samples of stars to be addressed, and we identify the requirements in terms of number of targets and the precision on the stellar properties that are needed to tackle such questions. By quantifying the asteroseismic yields expected from PLATO for red-giant stars, we demonstrate that these requirements are within the capabilities of the current instrument design, provided that observations are sufficiently long to identify the evolutionary state and allow robust and precise determination of acoustic-mode frequencies. This will allow us to harvest data of sufficient quality to reach a 10% precision in age. This is a fundamental pre-requisite to then reach the more ambitious goal of a similar level of accuracy, which will only be possible if we have to hand a careful appraisal of systematic uncertainties on age deriving from our limited understanding of stellar physics, a goal which conveniently falls within the main aims of PLATO's core science.
  • Spatial and time scales of spin density fluctuations (SDF) as well as the fluctuation-dissipation theorem (FDT) were analyzed in 3d ferromagnets using ab initio linear response calculations of complete wavevector and energy dependence of the dynamic spin susceptibility tensor. We demonstrate that SDF are spread continuously over the entire Brillouin zone and while majority of them reside within the 3d bandwidth, a significant amount comes from much higher energies. A justification of adiabatic approximation in spin dynamics is provided. FDT is shown to have two main constituents: a minor low-energy spin wave contribution and a much larger high-energy component from local SDF. The obtained universal dependence of the on-site spin correlator (SC) and the related effective fluctuating moment on the 3d band population is further discussed.
  • Spin-dependent exchange-correlation energy functionals in use today depend on the charge density and the magnetization density: $E_{\rm xc}[\rho,{\bf m}]$. However, it is also correct to define the functional in terms of the curl of ${\bf m}$ for physical external fields: $E_{\rm xc}[\rho,\nabla\times{\bf m}]$. The exchange-correlation magnetic field, ${\bf B}_{\rm xc}$, then becomes source-free. We study this variation of the theory by uniquely removing the source term from local and generalized gradient approximations to the functional. By doing so, the total Kohn-Sham moments are improved for a wide range of materials for both functionals. Significantly, the moments for the pnictides are now in good agreement with experiment. We also predict dramatic differences in the spatial geometry of ${\bf B}_{\rm xc}$ for the pnictides. Our source-free method is simple to implement in all existing density functional theory codes.
  • The first transiting planetesimal orbiting a white dwarf was recently detected in K2 data of WD1145+017 and has been followed up intensively. The multiple, long, and variable transits suggest the transiting objects are dust clouds, probably produced by a disintegrating asteroid. In addition, the system contains circumstellar gas, evident by broad absorption lines, mostly in the u'-band, and a dust disc, indicated by an infrared excess. Here we present the first detection of a change in colour of WD1145+017 during transits, using simultaneous multi-band fast-photometry ULTRACAM measurements over the u'g'r'i'-bands. The observations reveal what appears to be 'bluing' during transits; transits are deeper in the redder bands, with a u'-r' colour difference of up to ~-0.05 mag. We explore various possible explanations for the bluing. 'Spectral' photometry obtained by integrating over bandpasses in the spectroscopic data in- and out-of-transit, compared to the photometric data, shows that the observed colour difference is most likely the result of reduced circumstellar absorption in the spectrum during transits. This indicates that the transiting objects and the gas share the same line-of-sight, and that the gas covers the white dwarf only partially, as would be expected if the gas, the transiting debris, and the dust emitting the infrared excess, are part of the same general disc structure (although possibly at different radii). In addition, we present the results of a week-long monitoring campaign of the system.
  • A joint theoretical and experimental investigation is performed to understand the underlying physics of laser-induced demagnetization in Ni and Co films. Experimentally dynamics of spins is studied by determining the time-dependent amplitude of the Voigt vector and theoretically ab-initio calculations are performed using time-dependent density functional theory. We demonstrate that overall spin-orbit induced spin-flips are the most significant contributors with super-diffusive spin transport playing very limited to no role. Our study highlights the material dependent nature of the demagnetization during the process of thermalization of non-equilibrium spins.
  • We calculated the QCD equation of state using Taylor expansions that include contributions from up to sixth order in the baryon, strangeness and electric charge chemical potentials. Calculations have been performed with the Highly Improved Staggered Quark action in the temperature range $T\in [135~{\rm MeV}, 330~{\rm MeV}]$ using up to four different sets of lattice cut-offs corresponding to lattices of size $N_\sigma^3\times N_\tau$ with aspect ratio $N_\sigma/N_\tau=4$ and $N_\tau =6-16$. The strange quark mass is tuned to its physical value and we use two strange to light quark mass ratios $m_s/m_l=20$ and $27$, which in the continuum limit correspond to a pion mass of about $160$ MeV and $140$ MeV espectively. Sixth-order results for Taylor expansion coefficients are used to estimate truncation errors of the fourth-order expansion. We show that truncation errors are small for baryon chemical potentials less then twice the temperature ($\mu_B\le 2T$). The fourth-order equation of state thus is suitable for the modeling of dense matter created in heavy ion collisions with center-of-mass energies down to $\sqrt{s_{NN}}\sim 12$ GeV. We provide a parametrization of basic thermodynamic quantities that can be readily used in hydrodynamic simulation codes. The results on up to sixth order expansion coefficients of bulk thermodynamics are used for the calculation of lines of constant pressure, energy and entropy densities in the $T$-$\mu_B$ plane and are compared with the crossover line for the QCD chiral transition as well as with experimental results on freeze-out parameters in heavy ion collisions. These coefficients also provide estimates for the location of a possible critical point. We argue that results on sixth order expansion coefficients disfavor the existence of a critical point in the QCD phase diagram for $\mu_B/T\le 2$ and $T/T_c(\mu_B=0) > 0.9$.