Source author record

Joonas Nättilä

Joonas Nättilä appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

8works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

8 published item(s)

preprint2022arXiv

Multimessenger constraints for ultra-dense matter

Recent rapid progress in neutron-star (NS) observations offers great potential to constrain the properties of strongly interacting matter under the most extreme conditions. In order to fully exploit the current observational inputs and to study the impact of future observations of NS masses, radii, and tidal deformabilities, we analyze a large ensemble of randomly generated viable NS-matter equations of state (EoSs) and the corresponding rotating stellar structures. We discuss the compatibility and impact of various hypotheses and measurements on the EoS, including those involving the merger product of the gravitational-wave (GW) event GW170817, the binary merger components in GW190814, and radius measurements of the pulsar PSR J0740+6620. We obtain an upper limit for the dimensionless spin of a rigidly rotating NS, $|χ| < 0.81$, an upper limit for the compactness of a NS, $GM/(Rc^2) < 0.33$, and find that the conservative hypothesis that the remnant in GW170817 ultimately collapsed to a black hole strongly constrains the EoS and the maximal mass of NSs, implying $M_\text{TOV} < 2.53 M_\odot$ (or $M_\text{TOV} < 2.19 M_\odot$ if we assume that a hypermassive NS was created). Additionally, we derive a novel lower limit for the tidal deformability as a function of NS mass and provide fitting formulae that can be used to set priors for parameter estimation and to discern whether neutron stars or other compact objects are involved in future low-mass GW events. Finally, we find that the recent NICER results for the radius of the massive NS PSR J0740+6620 place strong constraints for the behavior of the EoS, and that the indicated radius values $R(2M_\odot) \gtrsim 11$ km are compatible with moderate speeds of sound in NS matter and thus with the existence of quark matter cores in massive NSs.

preprint2021arXiv

A fully-kinetic model for orphan gamma-ray flares in blazars

Blazars emit a highly-variable non-thermal spectrum. It is usually assumed that the same non-thermal electrons are responsible for the IR-optical-UV emission (via synchrotron) and the gamma-ray emission (via inverse Compton). Hence, the light curves in the two bands should be correlated. Orphan gamma-ray flares (i.e., lacking a luminous low-frequency counterpart) challenge our theoretical understanding of blazars. By means of large-scale two-dimensional radiative particle-in-cell simulations, we show that orphan gamma-ray flares may be a self-consistent by-product of particle energization in turbulent magnetically-dominated pair plasmas. The energized particles produce the gamma-ray flare by inverse Compton scattering an external radiation field, while the synchrotron luminosity is heavily suppressed since the particles are accelerated nearly along the direction of the local magnetic field. The ratio of inverse Compton to synchrotron luminosity is sensitive to the initial strength of turbulent fluctuations (a larger degree of turbulent fluctuations weakens the anisotropy of the energized particles, thus increasing the synchrotron luminosity). Our results show that the anisotropy of the non-thermal particle population is key to modeling the blazar emission.

preprint2020arXiv

Combining Electromagnetic and Gravitational-Wave Constraints on Neutron-Star Masses and Radii

We perform a joint Bayesian inference of neutron-star mass and radius constraints based on GW170817, observations of quiescent low-mass X-ray binaries (QLMXBs), photospheric radius expansion X-ray bursts (PREs), and X-ray timing observations of J0030+0451. With this data set, the form of the prior distribution still has an impact on the posterior mass-radius (MR) curves and equation of state (EOS), but this impact is smaller than recently obtained when considering QLMXBs alone. We analyze the consistency of the electromagnetic data by including an "intrinsic scattering" contribution to the uncertainties, and find only a slight broadening of the posteriors. This suggests that the gravitational-wave and electromagnetic observations of neutron-star structure are providing a consistent picture of the neutron-star mass-radius curve and the EOS.

preprint2020arXiv

Kilohertz quasi-periodic oscillations from neutron star spreading layers

When the accretion disc around a weakly magnetised neutron star (NS) meets the stellar surface, it should brake down to match the rotation of the NS, forming a boundary layer. As the mechanisms potentially responsible for this braking are apparently inefficient, it is reasonable to consider this layer as a spreading layer (SL) with negligible radial extent and structure. We perform hydrodynamical 2D spectral simulations of an SL, considering the disc as a source of matter and angular momentum. Interaction of new, rapidly rotating matter with the pre-existing, relatively slow material co-rotating with the star leads to instabilities capable of transferring angular momentum and creating variability on dynamical timescales. For small accretion rates, we find that the SL is unstable for heating instability that disrupts the initial latitudinal symmetry and produces large deviations between the two hemispheres. This instability also results in breaking of the axial symmetry as coherent flow structures are formed and escape from the SL intermittently. At enhanced accretion rates, the SL is prone to shearing instability and acts as a source of oblique waves that propagate towards the poles, leading to patterns that again break the axial symmetry. We compute artificial light curves of an SL viewed at different inclination angles. Most of the simulated light curves show oscillations at frequencies close to 1kHz. We interpret these oscillations as inertial modes excited by shear instabilities near the boundary of the SL. Their frequencies, dependence on flux, and amplitude variations can explain the high-frequency pair quasi-periodic oscillations observed in many low-mass X-ray binaries.

preprint2020arXiv

Magnetospheric return-current-heated atmospheres of rotation-powered millisecond pulsars

We computed accurate atmosphere models of rotation-powered millisecond pulsars in which the polar caps of a neutron star (NS) are externally heated by magnetospheric return currents. The external ram pressure, energy losses, and stopping depth of the penetrating charged particles were computed self-consistently with the atmosphere model, instead of assuming a simplified deep-heated atmosphere in radiative equilibrium. We used exact Compton scattering formalism to model the properties of the emergent X-ray radiation. The deep-heating approximation was found to be valid only if most of the heat originates from ultra-relativistic bombarding particles with Lorentz factors of $γ\gtrsim 100$. In the opposite regime, the atmosphere attains a distinct two-layer structure with an overheated optically thin skin on top of an optically thick cool plasma. The overheated skin strongly modifies the emergent radiation: it produces a Compton-upscattered high-energy tail in the spectrum and alters the radiation beaming pattern from limb darkening to limb brightening for emitted hard X-rays. This kind of drastic change in the emission properties can have a significant impact on the inferred NS pulse profile parameters as performed, for example, by Neutron star Interior Composition ExploreR. Finally, the connection between the energy distribution of the return current particles and the atmosphere emission properties offers a new tool to probe the exact physics of pulsar magnetospheres.

preprint2015arXiv

Models of neutron star atmospheres enriched with nuclear burning ashes

Low-mass X-ray binaries hosting neutron stars (NS) exhibit thermonuclear (type-I) X-ray bursts, which are powered by unstable nuclear burning of helium and/or hydrogen into heavier elements deep in the NS "ocean". In some cases the burning ashes may rise from the burning depths up to the NS photosphere by convection, leading to the appearance of the metal absorption edges in the spectra, which then force the emergent X-ray burst spectra to shift toward lower energies. These effects may have a substantial impact on the color correction factor $f_c$ and the dilution factor $w$, the parameters of the diluted blackbody model $F_E \approx w B_E(f_c T_{eff})$ that is commonly used to describe the emergent spectra from NSs. The aim of this paper is to quantify how much the metal enrichment can change these factors. We have developed a new NS atmosphere modeling code, which has a few important improvements compared to our previous code required by inclusion of the metals. The opacities and the internal partition functions (used in the ionization fraction calculations) are now taken into account for all atomic species. In addition, the code is now parallelized to counter the increased computational load. We compute a detailed grid of atmosphere models with different exotic chemical compositions that mimic the presence of the burning ashes. From the emerging model spectra we compute the color correction factors $f_c$ and the dilution factors $w$ that can then be compared to the observations. We find that the metals may change $f_c$ by up to about 40%, which is enough to explain the scatter seen in the blackbody radius measurements. The presented models open up the possibility for determining NS mass and radii more accurately, and may also act as a tool to probe the nuclear burning mechanisms of X-ray bursts.

preprint2014arXiv

The effect of accretion on the measurement of neutron star mass and radius in the low-mass X-ray binary 4U 1608-52

Spectral measurements of thermonuclear (type-I) X-ray bursts from low mass X-ray binaries have been used to measure neutron star (NS) masses and radii. A number of systematic issues affect such measurements and have raised concerns as to the robustness of the methods. We present analysis of the X-ray emission from bursts observed from 4U 1608-52 at various persistent fluxes. We find a strong dependence of the burst properties on the flux and spectral hardness of the persistent emission before burst. Bursts occurring during the low-accretion rate (hard) state exhibit evolution of the black body normalisation consistent with the theoretical predictions of NS atmosphere models. However, bursts occurring during the high-accretion rate (soft) state show roughly constant normalisation, which is inconsistent with the NS atmosphere models and therefore these bursts cannot be easily used to determine NS parameters. We analyse the hard-state burst to put the lower limit on the neutron star radius in 4U 1608-52 of 13 km (for masses 1.2-2.4 M$_\odot$). The best agreement with the theoretical NS mass-radius relations is achieved for source distances in the range 3.1-3.7 kpc. We expect that the radius limit will be 10 per cent lower if spectral models including rapid rotation are used instead.

preprint2014arXiv

The influence of accretion geometry on the spectral evolution during thermonuclear (type-I) X-ray bursts

Neutron star (NS) masses and radii can be estimated from observations of photospheric radius-expansion X-ray bursts, provided the chemical composition of the photosphere, the spectral colour-correction factors in the observed luminosity range, and the emission area during the bursts are known. By analysing 246 X-ray bursts observed by the Rossi X-ray Timing Explorer from 11 low-mass X-ray binaries, we find a dependence between the persistent spectral properties and the time evolution of the black body normalisation during the bursts. All NS atmosphere models predict that the colour-correction factor decreases in the early cooling phase when the luminosity first drops below the limiting Eddington value, leading to a characteristic pattern of variability in the measured blackbody normalisation. However, the model predictions agree with the observations for most bursts occurring in hard, low-luminosity, 'island' spectral states, but rarely during soft, high-luminosity, 'banana' states. The observed behaviour may be attributed to the accretion flow, which influences cooling of the NS preferentially during the soft state bursts. This result implies that only the bursts occurring in the hard, low-luminosity spectral states can be reliably used for NS mass and radius determination.